Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “peak finding”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Interactions Enhance Ramp Reversal Memory in Locally Phase Separated Materials

The ramp-reversal memory (RRM) effect in metal–insulator transition metal oxides (TMOs), a non-volatile resistance change induced by repeated temperature cycling, has attracted considerable interest in neuromorphic computing and non-volatile memory devices. Our previous defect motion model successfully explained RRM in vanadium dioxide (VO 2 ), capturing observed critical temperature shifts and memory accumulation throughout the sample. However, this approach lacked interactions between metallic and insulating domains. Here, we extend our model by combining a correlated Random Field Ising Model with defect diffusion-segregation, enabling accurate hysteresis modeling while predicting the relationship between RRM and domain interactions. Our simulations demonstrate that the maximum RRM occurs when the turnaround temperature approaches the inflection point. This peak in RRM vs. turnaround temperature is consistent with prior transport measurements, as well as our own optical measurements reported here. Significantly, we find that increasing nearest-neighbor interactions enhances the maximum memory effect, thus providing a clear mechanism for optimizing RRM performance. Since our model employs minimal assumptions, we predict that RRM should be a widespread phenomenon in materials exhibiting patterned phase coexistence of electronic domains. This work not only advances fundamental understanding of memory behavior in TMOs but also establishes a much-needed theoretical framework for optimizing device applications.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Homogenization of Dendritic Structures and the High-Temperature Strength of the Refractory High-Entropy Alloy MoNbTaVW

MoNbTaVW, a pioneering refractory high-entropy alloy (RHEA), is renowned for its exceptional strength at elevated temperatures. Like most RHEAs, it solidifies into a dendritic microstructure with steep concentration gradients, necessitating heat treatment for equilibration. Steep concentration gradients typically complicate mechanical property analysis and modeling, key challenges for high-throughput alloy discovery. This study examines the effects of high-temperature homogenization (1800 °C for 8 hours) on the microstructure and mechanical properties of MoNbTaVW across temperatures up to 1200 °C. Postmortem microstructural analysis, coupled with chemical, crystallographic, and mechanical assessments, revealed that yield strength was surprisingly insensitive to homogenization, with less than a 3 pct difference between as-cast and heat-treated materials. Nanoindentation confirmed minimal nanohardness changes across the dendritic structure at room temperature. However, homogenization significantly enhanced high-temperature work hardening, producing higher peak compressive strength. Both as-cast and homogenized MoNbTaVW exhibited room-temperature strengths of 1400 MPa, exceeding previously reported values. These findings demonstrate that high-temperature treatment enables microstructural homogenization without compromising strength, a unique behavior among RHEAs. In conclusion, while not universally applicable, this insight highlights the importance of understanding microstructural development in experimental alloys and offers a pathway to design alloys that achieve homogeneous-like properties without heat treatment.

Rietema, C. J. [Lawrence Livermore National Labora↗
DOE OSTI2025

Abiotic influences on continuous conifer forest structure across a subalpine watershed

Understanding the abiotic drivers of high-elevation forest physiognomy is essential for forecasting how mountain ecosystems will respond to emerging environmental pressures. Most prior studies of these relationships have relied on small samples of the full landscape, resulting in limited power to detect dominant covariates and their interactions. Here we report the first evaluation of abiotic influences on a complement of accurate, wall-to-wall estimates of conifer forest structure and composition at the watershed scale. In a subalpine conifer domain in the Colorado Rocky Mountains (USA), we developed a novel method for deriving stand structure metrics from waveform LiDAR data, which showed high fidelity with field inventory. We quantified the relationships between structural and compositional metrics and climate, topographic, edaphic, and geologic factors. Our results showed that peak snow water equivalent (SWE), snow disappearance rate, and elevation explained most of the variation in forest structure. The highest stand density, basal area, maximum canopy height, and quadratic mean diameter occurred in sites with SWE around one standard deviation below mean, but with long snow residence times. Stand density decreased linearly with elevation, while other metrics peaked between 3000 m and 3200 m. Substrate properties had weaker influence. Continuous mapping of through-canopy forest structure enabled our novel findings of the dominant role of snowpack in explaining structural and compositional variation, and of elevation thresholds. Our reproducible approach facilitates assessment of forest-topoclimate relationships in other conifer-dominated landscapes and improves understanding of the baseline patterns controlling forest structure, which is needed for predicting long-term ecological change.

forest ecology↗
DOE OSTI2026

Pedestal origin and extrapolation of high-density small edge-localised-modes peak parallel energy fluence in ITER and SPARC

Experimental analysis and simulations with the BOUT++ code show that small edge-localised modes (ELMs) in reactor-relevant high-density regimes originate in a region close to the separatrix and only marginally perturb the pedestal structure. The measured divertor peak parallel energy fluence (ε ∥,peak ) for a database of small ELM scenarios in DIII-D and ASDEX Upgrade can be reproduced, within 40 % accuracy on average, if an ad hoc modification of the Eich peak parallel ELM energy fluence model is applied to account for the small ELM pedestal birth location. This allows for first-order extrapolation of small-ELM divertor ε ∥,peak to ITER and SPARC, resulting in values that satisfy the nominal melting threshold of tungsten monoblocks of 12 MJ m −2 . The findings reported in this study, both via modelling and direct measurements, constitute a step forward in assessing small ELMs in high edge-collisionality scenarios as a viable plasma regime for the operation of next-generation fusion machines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Dust dominates high-altitude snow darkening and melt over high-mountain Asia

Westerly-driven, long-range transportation of dust particles in elevated aerosol layers (EALs) is a persistent phenomenon during spring and summer over the Indian subcontinent. During the snow accumulation season, EALs transport substantial dust to the snow-covered slopes of high-mountain Asia (HMA). Here we use multiple satellite-based estimates to demonstrate a robust physical association between the EALs and dust-induced snow darkening over HMA. Results from a fully-coupled atmosphere-chemistry-snow model support these observations, revealing a signature of increasing dust-induced snow darkening with surface elevation over HMA that peaks near 4500 m. Moreover, the influence of dust on snow darkening is greater than that of black carbon above 4000 m. Our findings suggest a discernable role of dust in the observed spatial heterogeneity of snowmelt and snowline trends over HMA and highlight an increasing contribution of dust to snowmelt as the snowline rises with warming.

Elevated aerosol layers (EAL), Light absorbing par↗
DOE OSTI2025

Morphological Characterization of Uranyl Fluoride Particles via Atomic Force Microscopy

Uranium hexafluoride (UF 6 ) undergoes a rapid hydrolysis reaction when exposed to atmospheric water. In addition to producing hazardous HF gas, the hydrolysis reaction produces uranyl fluoride (UO 2 F 2 ), a radioactive solid phase particulate material. Because of the technological utility of UF 6 in the nuclear fuel cycle, understanding the transport properties of UO 2 F 2 aerosol produced via UF 6 hydrolysis is important for accident scenarios. Moreover, the fundamental chemical and physical properties of the UF 6 hydrolysis reaction are not completely understood. Recently, several experiments on the aerosol phase properties of UO 2 F 2 produced in this way have shown that under most relevant conditions, the particle size distribution (PSD) of UO 2 F 2 can be extremely small, approximately 3 to 5 nm, which is well below the threshold that can be routinely observed via scanning electron microscopy (SEM). Although readily observable in the aerosol phase, observation of nanometer-sized particles in the condensed phase (i.e. deposited on surfaces) remains a challenge. Here, in this study, we have used atomic force microscopy (AFM) to study the PSD and morphological characteristics of UO 2 F 2 deposited at low and high concentrations under different humidity conditions, a primary variable in the hydrolysis reaction. Here, we find strong agreement between PSD measured in the aerosol phase via scanning mobility particle sizing and PSD measured via AFM, with particle sizes peaked below 4 nm for low-humidity conditions. At higher humidity, the distribution is centered around 5 to 10 nm but extends up to 20 nm. These results are in stark contrast to previous measurements using SEM that show PSD on the order of 300- to 1000-nm particle sizes; moreover, these are the first direct measurements of individual particles of UO 2 F 2 having been produced via UF 6 hydrolysis deposited on surfaces. These measurements, therefore, open a new avenue for collecting and detecting UO 2 F 2 in the condensed phase and further refine the PSD, which is critical for environmental transport determinations.

Uranium hexafluoride↗
DOE OSTI2024

Warming effects on grassland soil microbial communities are amplified in cool months

Abstract Global warming modulates soil respiration (RS) via microbial decomposition, which is seasonally dependent. Yet, the magnitude and direction of this modulation remain unclear, partly owing to the lack of knowledge on how microorganisms respond to seasonal changes. Here, we investigated the temporal dynamics of soil microbial communities over 12 consecutive months under experimental warming in a tallgrass prairie ecosystem. The interplay between warming and time altered (P < 0.05) the taxonomic and functional compositions of microbial communities. During the cool months (January to February and October to December), warming induced a soil microbiome with a higher genomic potential for carbon decomposition, community-level ribosomal RNA operon (rrn) copy numbers, and microbial metabolic quotients, suggesting that warming stimulated fast-growing microorganisms that enhanced carbon decomposition. Modeling analyses further showed that warming reduced the temperature sensitivity of microbial carbon use efficiency (CUE) by 28.7% when monthly average temperature was low, resulting in lower microbial CUE and higher heterotrophic respiration (Rh) potentials. Structural equation modeling showed that warming modulated both Rh and RS directly by altering soil temperature and indirectly by influencing microbial community traits, soil moisture, nitrate content, soil pH, and gross primary productivity. The modulation of Rh by warming was more pronounced in cooler months compared to warmer ones. Together, our findings reveal distinct warming-induced effects on microbial functional traits in cool months, challenging the norm of soil sampling only in the peak growing season, and advancing our mechanistic understanding of the seasonal pattern of RS and Rh sensitivity to warming.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Material Effects on Electron-Capture Decay in Cryogenic Sensors

Several current searches for physics beyond the standard model are based on measuring the electron-capture (EC) decay of radionuclides implanted into cryogenic high-resolution sensors. The sensitivity of these experiments has already reached the level where systematic effects related to atomic state energy changes from the host material are a limiting factor. One example is a neutrino mass study based on the nuclear EC decay of 7 Be to 7 Li inside cryogenic Ta-based sensors. Here, to understand the material effects at the required level, we use density-functional theory to model the electronic structure of lithium atoms in different atomic environments of the polycrystalline Ta absorber film. The calculations reveal that the Li 1s binding energies can vary by more than 2 eV due to insertion at different lattice sites, at grain boundaries, in disordered Ta, and in the vicinity of various impurities. However, the total range of Li 1s shifts does not exceed 4 eV, even for extreme amorphous disorder. Furthermore, when investigating the effects on the Li 2s levels, we find broadening of more than 5 eV due to hybridization with the Ta band structure. Material effects are shown to contribute significantly to peak broadening in Ta-based sensors that are used to search for physics beyond the standard model in the EC decay of 7 Be, but they do not explain the full extent of observed broadening. Understanding these in-medium effects will be required for current- and future-generation experiments that observe low-energy radiation from the EC decay of implanted isotopes to evaluate potential limitations on the measurement sensitivity.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Carbon complexes in highly C-doped GaN

We investigate the properties of heavily C-doped GaN grown by hydride vapor phase epitaxy using both optical experiments and hybrid density functional theory calculations. Previous work has established that carbon acceptors (C N ) give rise to a yellow luminescence band near 2.2 eV along with a blue luminescence band near 2.9 eV. Photoluminescence measurements show the yellow band shifting as a function of carbon concentration, suggesting a change in the behavior of carbon species as carbon content increases. With hybrid density functional theory we calculate the electrical and optical behavior of carbon centers containing multiple carbon impurities, which may arise in heavily doped material. We compare the behavior of these complexes to the isolated centers, and find that the dicarbon donor-acceptor (C Ga –C N ) complex is a candidate to explain the shift in the yellow luminescence peak. Tricarbon complexes have high formation energies and modest binding energies, and also give rise to optical transitions that are inconsistent with the observed spectra. We also identify the split dicarbon interstitial on the gallium site as a low-energy species with a large binding energy that may act to compensate carbon acceptors. Here, local vibrational modes are calculated for carbon impurity centers, and we compare these results to recent experiments. Dicarbon and tricarbon complexes involving C Ga and C N exhibit modes that are only slightly higher than the isolated species, while carbon interstitials and related complexes give rise to vibrational modes much higher than C Ga and C N .

36 MATERIALS SCIENCE↗
DOE OSTI2025

Cosmological quasiparticles and the cosmological collider

The interplay between cosmology and strongly coupled dynamics can yield transient spectral features that vanish at late times, but which may leave behind phenomenological signatures in the spectrum of primordial fluctuations. Of particular interest are strongly coupled extensions of the standard model featuring approximate conformal invariance. In flat space, the spectral density for a scalar operator in a conformal field theory is characterized by a continuum with scaling law governed by the dimension of the operator, and is otherwise featureless. Anti–de Sitter/conformal field theory arguments suggest that for large 𝑁, in an inflationary background with Hubble rate 𝐻, this continuum is gapped. We demonstrate that there can be additional peak structures that become sharp and particlelike at phenomenologically interesting regions in parameter space, and we estimate their contribution to cosmological observables. We find phenomena that are potentially observable in future experiments that are unique to these models, including displaced oscillatory features in the squeezed limit of the bispectrum. These particles can be either fundamental, and localized to a UV brane, or composite at the Hubble scale, 𝐻, and bound to a horizon in the bulk of the 5D geometry. We comment on how stabilization of conformal symmetry breaking vacua can be correlated with these spectral features and their phenomenology.

conformal field theory↗
DOE OSTI2025

Deep Neural Network Assisted Distributed Strain and Temperature Fiber Sensor System for Natural Gas Pipeline Monitoring

Natural gas pipeline integrity monitoring is crucial to detect potential leaks, find structural issues, and prevent environmental damage. This article presents a system of natural gas pipeline monitoring that uses a specialized double Brillouin peak sensing fiber along with the Brillouin optical time domain analysis (BOTDAs) technique. The calibrated sensing fiber coefficients for strain and temperature are 41.8 kHz/ με and 0.9 MHz/°C for peak 1; and 47.2 kHz/ με , and 1.11 MHz/°C for peak 2, respectively. Initially, lab tests were performed by installing a short section of double Brillouin peak fiber (DBPF) on a 1-in steel pipe under pressure up to 1000 per square inch (psi) at elevated temperatures. Simultaneous distributed measurements of temperature and pressure-induced hoop strain were successfully measured. Considering the long processing speed to extract Brillouin frequency shift (BFS), we employ a novel probabilistic deep neural network (PDNN) framework for rapid BFS prediction. Additionally, using the Finite Element Method, the effects of the pipeline pressure on hoop strain were modeled and compared to the experimental hoop strain under the same set of pipeline conditions. Finally, an actual 4-in outer diameter steel natural gas pipeline was used for pilot-scale tests, where hoop strain was measured at various pressure levels. Leaks were simulated to demonstrate accurate pipeline integrity monitoring. At an internal pipe pressure of 1000 psi, hoop strain of approximately 300 με was observed, and the sensitivity was calculated as 0.28 με /psi. The results of this pilot-scale study demonstrated that the system is capable of performing distributed monitoring sufficient to detect pipeline pressure and the presence of leaks to ensure the safe operation of gas pipelines in the field.

03 NATURAL GAS↗
DOE OSTI2021

Simulation and Modeling of Time-Resolved X-Ray Detector for the Saturn Accelerator

In this work, we present the technology-aided computer design (TCAD) device simulation and modeling of a silicon p-i-n diode for detecting time-dependent X-ray radiation. We show that the simulated forward and reverse breakdown current–voltage characteristics agree well with the measured data under nonradiation environment by only calibrating carrier lifetimes for the forward bias case and avalanche model critical fields for the reverse bias condition. Using the calibrated parameters and other nominal material properties, we simulated the radiation responses of the p-i-n diode and compared with experimental data when the diode was exposed to X-ray radiation at Sandia’s Saturn facility and the Idaho State University (ISU) TriMeV facility. For Saturn’s Gaussian dose-rate pulses, we show three findings from TCAD simulations. First, the simulated photocurrents are in excellent agreement with the measured data for two dose-rate pulses with peak values of 1.16×10 10 and 1.88×10 10 rad(Si)/s. Second, the simulation results of high dose-rate pulses predict increased delayed photocurrents with longer time tails in the diode electrical responses due to excess carrier generation. Third, simulated peak values of diode radiation responses versus peak dose rates at different bias conditions provide useful guidance to determine the dose-rate range that the p-i-n diode can reliably detect in experiment. For TriMeV’s non-Gaussian dose-rate pulse, our simulated diode response is in decent agreement with the measured data without further calibration. We also studied the effects of device geometry, recombination process, and dose-rate enhancement via TCAD simulations to understand the higher measured response in the time after the peak dose-rate radiation for the p-i-n diode exposed to TriMeV irradiation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Combining Astrometry and Elemental Abundances: The Case of the Candidate Pre-Gaia Halo Moving Groups G03-37, G18-39, and G21-22

While most moving groups are young and nearby, a small number have been identified in the Galactic halo. Understanding the origin and evolution of these groups is an important piece of reconstructing the formation history of the halo. Here we report on our analysis of three putative halo moving groups: G03-37, G18-39, and G21-22. Based on Gaia EDR3 data, the stars associated with each group show some scatter in velocity (e.g., Toomre diagram) and integrals of motion (energy, angular momentum) spaces, counter to expectations of moving-group stars. We choose the best candidate of the three groups, G21-22, for follow-up chemical analysis based on high-resolution spectroscopy of six presumptive members. Using a new Python code that uses a Bayesian method to self-consistently propagate uncertainties from stellar atmosphere solutions in calculating individual abundances and spectral synthesis, we derive the abundances of α- (Mg, Si, Ca, Ti), Fe-peak (Cr, Sc, Mn, Fe, Ni), odd-Z (Na, Al, V), and neutron-capture (Ba, Eu) elements for each star. We find that the G21-22 stars are not chemically homogeneous. Based on the kinematic analysis for all three groups and the chemical analysis for G21-22, we conclude the three are not genuine moving groups. The case for G21-22 demonstrates the benefit of combining kinematic and chemical information in identifying conatal populations when either alone may be insufficient. Comparing the integrals of motion and velocities of the six G21-22 stars with those of known structures in the halo, we tentatively associate them with the Gaia-Enceladus accretion event.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Sandblasting the r -process: Spallation of Ejecta from Neutron Star Mergers

Neutron star mergers (NSMs) are rapid neutron-capture (r-process) nucleosynthesis sites that expel matter at high velocities, from 0.1c to as high as 0.6c. Nuclei ejected at these speeds are sufficiently energetic to initiate spallation nuclear reactions with interstellar medium (ISM) particles. We adopt a thick-target model for the propagation of high-speed heavy nuclei in the ISM, similar to the transport of cosmic rays. We find that spallation may create observable perturbations to NSM isotopic abundances, particularly around the low-mass edges of the r-process peaks where neighboring nuclei have very different abundances. The extent to which spallation modifies the final NSM isotopic yields depends on: (1) the ejected abundances, which are determined by the NSM astrophysical conditions and the properties of nuclei far from stability, (2) the ejecta velocity distribution and propagation in interstellar matter, and (3) the spallation cross sections. Observed solar and stellar r-process yields could thus constrain the velocity distribution of ejected neutron star matter, assuming NSMs are the dominant r-process source. We suggest avenues for future work, including measurement of relevant cross sections.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

EDEN: Exploring Disks Embedded in N-body Simulations of Milky Way–mass Halos from Symphony

We investigate the impact of galactic disks on the tidal stripping of cold dark matter subhalos within Milky Way (MW)-mass halos (M vir ∼ 10 12 M ⊙ ) using a new simulation suite, EDEN. By re-simulating 45 MW-mass zoomed-in halos from the N-body Symphony compilation with embedded disk potentials, which evolve according to star formation histories predicted by the UniverseMachine model, we self-consistently tie disk growth to halo accretion rate and significantly expand the range of disk masses and formation histories studied. We use the particle-tracking-based subhalo finder Symfind to enhance the robustness of subhalo tracking. We find that disks near the median disk-to-halo mass ratio of our sample (M *,Disk /M vir,host = 2%) reduce subhalo peak mass functions within 100 kpc by about 10% for peak masses above 10 8 M ⊙ . Heavier, MW/M31-like disks (M *,Disk /M vir,host ≳ 5%) lead to a reduction of more than 40%. Subhalo abundance suppression is more pronounced near halo centers, with particularly enhanced stripping for subhalos accreted over 8 Gyr ago on orbits with pericenters <100 kpc. Suppression is further amplified when disk mass is increased within fixed halo and disk assembly histories. In all cases, the suppression we measure should be interpreted as stripping below the mass resolution limit rather than complete subhalo disruption. This study reshapes our understanding of the MW’s impact on its satellites, suggesting it strips subhalos more efficiently than typical MW-mass galaxies due to its larger disk-to-halo mass ratio and earlier disk formation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Raman spectroscopy of neutron irradiated silicon carbide

The effects of neutron irradiation on microstructural evolution and the resultant changes in physical and mechanical properties are of critical importance for the development of silicon carbide (SiC) materials for nuclear applications. This study neutron-irradiated βSiC under a wide range of conditions at temperatures between 235 and 750°C and neutron doses of 0.01–11.8 displacements per atom, and then evaluated the effects on the SiC structure using Raman spectroscopy. The SiC optical phonon lines were shifted to lower wavenumbers by irradiation. Correlations were found among the wavenumber of the longitudinal optical phonon line, irradiation-induced swelling, and irradiation temperature. The peak shift also correlated indirectly with decreasing thermal conductivity of irradiated SiC. The irradiation-induced peak shift is explained by combinations of lattice strain, reduction of the elastic modulus, and other factors including decreasing coherent domain size. These findings bridge irradiation-induced microstructural changes and property changes and illustrate how Raman spectroscopy is a useful tool for nondestructively assessing irradiated SiC materials for nuclear applications.

Koyanagi, T.↗
DOE OSTI2026

Use and Siting of Electric Vehicle Charging Stations in Juneau, Alaska

This report details a study of electric vehicle (EV) Level 2 charging stations in Juneau, Alaska. Utilization analyses of six public over five years and 250 residential chargers over two years are included, and a composite score is introduced to identify optimal locations for future charging stations that target residents of manufactured and multifamily housing (MMFH) in Juneau. We find that public charging station usage is very location-dependent, with three chargers in use more than 60% of days during the peak hour of the day (which ranges from 10 a.m. to 7 p.m.), including a charger near residential housing, illuminating potential needs for additional public chargers in those areas. Residential charging utilization typically occurs overnight - opposite to most public charging stations analyzed - and spikes after 10 p.m. This suggests that Alaska Electric Light & Power Company's time-of-use charging program, which lowers electricity rates at 10 p.m. to incentivize overnight charging, is very effective. Residential charging data also show that households tend to charge 15 hours per week, or 9% of the time, meaning that multiple households could likely share one charger if one were provided near MMFH locations. This is supported by residential charging session analysis, which shows that the median household has around two night charging sessions per week. The EV siting analysis identifies areas of high housing density, low access to public chargers, and unconstrained feeders. A cluster of MMFH parcels in Douglas demonstrated the highest composite scores considering all factors, being the only area to have a perfect score of 2.25.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Electronic density of states of a U ( 1 ) quantum spin liquid with spinon Fermi surface. I. Orbital magnetic field effects

A quantum spin liquid with a spinon Fermi surface is an exotic insulator that hosts neutral Fermi surfaces inside the insulating gap. In an external magnetic field, it has been pointed out that the neutral Fermi surfaces are Landau quantized to form Landau levels due to the coupling to the induced emergent gauge magnetic field. In this paper, we calculate the electronic density of states (as observed in tunneling experiments) of the quantum spin liquid in an orbital magnetic field. Here, we find that the Landau levels from the neutral Fermi surfaces give rise to a set of steps emerging at the upper and lower Hubbard band edges. Each of the Hubbard band edge steps further develop into a band edge resonance peak when a weak gauge binding arising from the gauge field fluctuations is taken into account. Importantly, each Hubbard band edge step and its resulting resonance peak in the weak gauge binding are found to have a corresponding Landau level from the neutral Fermi surfaces, so the Hubbard band edge steps and the band edge resonance peaks provide signatures to the unique feature that characterizes the Landau quantization of the in-gap neutral Fermi surfaces in the spin liquid. We further consider the strong gauge-binding regime where the band edge resonance peaks move into the Mott gap and develop into true in-gap bound states. In the strong gauge binding regime, we solve the Landau-level spectrum of the in-gap bound states in an orbital magnetic field. For the in-gap bound state with a Mexican-hat-like band dispersion, we find that the envelop energy to have a state excited from the bound state Landau levels decreases quadratically with the magnetic field. The quadratic decrease behavior of the envelop energy is consistent with the intuition that applying magnetic field localizes the states and energetically promotes the in-gap bound states formation. Finally, we discuss the connection of our results to the electronic density of states spectra measured in the layered 1T-TaS 2 . We point out that a quantum spin liquid with a quasibound state in the upper Hubbard band can give the density of states spectra similar to the one measured in the experiment.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗