Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Random Effects”

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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Mineral Deposition on the Rough Walls of a Fracture

Modeling carbonate growth in fractures and pores is important for understanding carbon sequestration in the environment or when supersaturated solutions are injected into rocks. Here, we study the simple but nontrivial problem of calcite growth on fractures with rough walls of the same mineral using kinetic Monte Carlo simulations of attachment and detachment of molecules and scaling approaches. First, we consider wedge-shaped fracture walls whose upper terraces are in the same low-energy planes and show that the valleys are slowly filled by the propagation of parallel monolayer steps in the wedge sides. The growth ceases when the walls reach these low-energy configurations so that a gap between the walls may not be filled. Second, we consider fracture walls with equally separated monolayer steps (vicinal surfaces with roughness below 1 nm) and show that growth by step propagation will eventually clog the fracture gap. In both cases, scaling approaches predict the times to attain the final configurations as a function of the initial geometry and the step-propagation velocity, which is set by the saturation index. The same reasoning applied to a random wall geometry shows that step propagation leads to lateral filling of surface valleys until the wall reaches the low-energy crystalline plane that has the smallest initial density of molecules. Thus, the final configurations of the fracture walls are much more sensitive to the crystallography than to the roughness or the local curvature. The framework developed here may be used to determine those configurations, the times to reach them, and the mass of deposited mineral. Effects of transport limitations are discussed when the fracture gap is significantly narrowed.

calcite↗
DOE OSTI2021

Investigation of theoretical scaling laws using large eddy simulations for airborne spreading of viral contagion from sneezing and coughing

Using a set of large eddy point-particle simulations, we explore the fluid dynamics of an ejected puff resulting from a cough/sneeze. The ejection contains over 61,000 potentially virus-laden droplets at an injection Reynolds number of about 46,000, comparable to an actual cough/ sneeze. We observe that global puff properties, such as centroid, puff volume, momentum, and buoyancy vary little across realizations. Other properties, such as maximum extent, shape, and edge velocity of the puff, may exhibit substantial variation. In many realizations, a portion of the puff splits off and advances along a random direction, while keeping airborne droplet nuclei afloat. This peeled-off portion provides a mechanism for virus-laden droplets to travel over large distances in a short amount of time. We also observe that the vast majority of droplets remain suspended within the puff after all liquid has evaporated. The main objectives of the study are to (i) evaluate assumptions of Balachandar’s et al. theory [Int. J. Multiphase Flow 132, 103439 (2020)], which include buoyancy effects, shape of the puff, and droplet evaporation rate, (ii) obtain values of closure parameters, which include location and time of the virtual origin, and puff entrainment and drag coefficients, and (iii) evaluate the accuracy of the theory in predicting the shape, size, and location of the puff, as well as droplet number density long after ejection. The theory adequately predicts global puff properties including size, velocity, and distance traveled, the largest size of droplets that exit the puff due to settling, and the droplet size distribution within the puff long after ejection.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Characterization and thermometry of dissipatively stabilized steady states

In this work we study the properties of dissipatively stabilized steady states of noisy quantum algorithms, exploring the extent to which they can be well approximated as thermal distributions, and proposing methods to extract the effective temperature T. We study an algorithm called the relaxational quantum eigensolver (RQE), which is one of a family of algorithms that attempt to find ground states and balance error in noisy quantum devices. In RQE, we weakly couple a second register of auxiliary ‘shadow’ qubits to the primary system in Trotterized evolution, thus engineering an approximate zero-temperature bath by periodically resetting the auxiliary qubits during the algorithm’s runtime. Balancing the infinite temperature bath of random gate error, RQE returns states with an average energy equal to a constant fraction of the ground state. We probe the steady states of this algorithm for a range of base error rates, using several methods for estimating both T and deviations from thermal behavior. In particular, we both confirm that the steady states of these systems are often well-approximated by thermal distributions, and show that the same resources used for cooling can be adopted for thermometry, yielding a fairly reliable measure of the temperature. These methods could be readily implemented in near-term quantum hardware, and for stabilizing and probing Hamiltonians where simulating approximate thermal states is hard for classical computers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Data-Driven Preemptive Voltage Monitoring and Control Using Probabilistic Voltage Sensitivities

Increased penetration levels of distributed variable renewable generation can cause random voltage fluctuations and violations at multiple nodes. Traditional methods of voltage control typically involve reactionary responses of capacitor banks, tap changers, and recently even smart inverters. But because of the lack of foresight in voltage violations, these controls are ineffective to completely mitigate the issue. Therefore, new methods of predicting voltage violations subject to random power injection changes in the distribution network are needed, which can be used to guide optimal and dynamic methods of voltage control. This work lays the foundation for such preemptive voltage monitoring and control by proposing an analytical and sensor data-driven voltage sensitivity analysis method. Driven by stochastic data and forecasts, the method can be used to develop probabilistic voltage sensitivities and consequently to predict system nodes with high likelihood of voltage limit violations. The effectiveness of this method is tested on IEEE 69-node distribution system integrated with distributed solar. The results demonstrate the proposed method's ability to successfully predict nodes with high probability of voltage violations for a specific time-series simulation. The results also demonstrate the ability to guide timely power injection control actions to mitigate future voltage violations.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

Mu2e resonant extraction regulation system simulation in delivery ring

Mu2e is an upcoming experiment at Fermilab that relies on the slowly extracted 8 GeV proton beam from the Delivery Ring. The experiment imposes strong requirements on the spill uniformity. To address these requirements, the fast spill regulations system is being developed and commissioned. To inform this development and optimize the system performance we are carrying out the detailed simulations of the regulation process. The simulation includes the effect of six harmonic sextupoles that excite the third-integer resonance and three fast ramping quadrupoles that drive the horizontal tune to 29/3. The components of spill regulation system are designed to mitigate long-term drifts in the beam, ensuring stable operation over extended timescales, as well as addresses rapid variations within single spill. In this study, we review the regulation system design, simulation of the slow regulation, and the fast regulation PID regulation to curtail random variations in the extraction rate that could occur within a single spill.

Narayanan, Aakaash [Fermilab]↗
DOE OSTI2024

Emulation and detection of physical faults and cyber-attacks on building energy systems through real-time hardware-in-the-loop experiments

The increasing use of remote or mobile access, integrated wearable technologies, data exchange, and cloud-based data analytics in modern smart buildings is steering the building industry towards open communication technologies. The increased connectivity and accessibility could lead to more cyber-attacks in smart buildings. On the other hand, physical faults (e.g., HVAC -heating, ventilation, and air-conditioning faults) may have similar adverse impacts as those from the cyber-attacks on building energy systems, such as occupant discomfort, energy wastage, and equipment downtime. However, current physical behavior-based anomaly detection methods fail to differentiate between cyber-attacks and physical faults in building energy systems. Moreover, the challenge in collecting real-world threat data with ground truth has led researchers to rely on numerical models with user-defined assumptions, which may not accurately reflect real-world conditions due to the lack of in-situ experimental datasets. To address these challenges and gaps, this paper presents a flexible hardware-in-the-loop (HIL) testbed for generating cyber-attack and physical fault datasets and demonstrating threat detection algorithms in a real building automation system (BAS) environment. This testbed combines hardware (i.e., real BAS with local HVAC controllers and a physical network) with software (i.e., high-fidelity models to represent behaviors of building envelope and HVAC energy systems), enabling emulations of realistic threats. Five HIL experiments, including one baseline without any threats, two with physical faults, and two with cyber-attacks, were conducted to generate datasets containing detailed network traffic and system states. A joint classification framework, incorporating a network analyzer and a physical HVAC fault detector, was proposed to automatically detect cyber-physical abnormalities on BAS at both the network and the physical HVAC levels. The network analyzer comprises a conditional random fields (CRF) based command validator and a statistics-based detection strategy. The fault detector employs a weather and schedule-based pattern matching and feature-based principal component analysis (WPM-FPCA) method. Evaluation of the classification using four metrics from the multi-class confusion matrix revealed an average accuracy of 90.2%, recall of 89.7%, precision of 88.5% and F1-score of 89.2%. Finally, these results demonstrate that the proposed joint classification framework can effectively differentiate between specific types of cyber-attacks (e.g., device reinitialization attack, network Denial-of-Service attack) and physical faults (e.g., air handling unit operational fault, cooling coil valve stuck) in real time for improved building energy management.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Impact of chemical ordering on thermodynamic properties of point defects and Xe substitutional in U-10Mo

The accurate knowledge of defect energetics is critical to understanding the aging and irradiation behavior of U-10Mo nuclear fuel, which is selected as the fuel type for conversion of the United States High-Performance Research Reactors (HPRRs). Furthermore, using hybrid molecular dynamics and Monte Carlo (MDMC) simulation, we studied the impact of chemical ordering on the formation energies of vacancies, interstitials, and the solution energy of the Xe substitutionals. Instead of forming a random solid solution (RSS), substantial short-range-order (SRO) develops in U-10Mo, particularly at low temperatures. Mo atoms are found to repel each other and prefer U-rich local atomic environments within the 1st nearest neighbor (1NN) cutoff. Compared to the case of a RSS, the state with equilibrated Mo ordering shifts the distributions of vacancy and interstitial formation energies due to the dependence of defect energies on the local atomic environment, without a clear effect on Xe solution energy. In the operation temperature range (100–250 °C) of U-Mo fuels, neglecting SRO can lead to an inaccurate estimate of thermal equilibrium point defect concentrations by over an order of magnitude and incorrectly predict the preference among different types of dumbbells, highlighting the critical importance of accounting for the impact of chemical ordering for accurate atomistic calculations of defect properties.

Molecular dynamics and Monte Carlo (MDMC)↗
DOE OSTI2022

Redox Defect Thermochemistry of FeAl 2 O 4 Hercynite in Water Splitting from First-Principles Methods

Solar thermochemical hydrogen (STCH) production is a promising route to produce fuels from sunlight via high-temperature water splitting. However, efficient and technologically viable implementations of this process only allow a narrow window of thermodynamic boundary conditions that can be used to cycle the system, thus limiting the design space for suitable metal oxide redox mediators. An oxygen defect redox mechanism can contribute a favorable reduction entropy to expand this window, and computational evaluation of materials with high oxygen defect entropies could play a pivotal role in guiding the discovery and design of suitable oxides. This study employs first-principles calculations to investigate the redox mediating defect mechanism of the STCH candidate material, hercynite (FeAl 2 O 4 ). We compare the results of total energy calculations from density functional theory (DFT) with beyond-DFT approaches, including hybrid functionals and the random phase approximation, which are among the most advanced methods currently feasible for supercell defect calculations. Using the predicted formation energies, we perform thermodynamic modeling of FeAl2O4 reduction and oxidation via free energy minimization that incorporates ideal gas, configurational, and vibrational entropy contributions evaluated within the quasi-harmonic approximation. Special attention is devoted to understanding interactions among co-existing defects, such as the association of pairs and complexes of O vacancies and cation antisite defects, and the effect of mutually compensating defect charges. Our results corroborate the notion that the details of defect interactions can be decisive for the viability of hydrogen production within the desirable STCH process window.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Molecular bond stabilization in the strong-field dissociation of O 2 +

We theoretically examine the rotational and vibrational dynamics of O 2 + molecular ions exposed to intense, short laser pulses for conditions realized in contemporary pump-probe experiments. We solve the time-dependent Schrödinger equation within the Born-Oppenheimer approximation for an initial distribution of randomly aligned molecular ions. For fixed peak intensities, our numerical results show that total, angle-integrated O 2 + → O( 3 P) + O + ( 4 S 0 ) dissociation yields do not monotonically increase with increasing infrared-probe pulse duration. We find this pulse-duration-dependent stabilization to be consistent with the transient trapping of nuclear probability density in a light-induced (bond-hardening) potential-energy surface and robust against rotational excitation. We analyze this stabilization effect and its underlying bond-hardening mechanism (i) in the time domain, by following the evolution of partial nuclear probability densities associated with the dipole-coupled O 2 + (a 4 Π u ) and O 2 + (f 4 Π g ) cationic states, and (ii) in the frequency domain, by examining rovibrational quantum-beat spectra for the evolution of the partial nuclear probability densities associated with these states. In this work, our analysis reveals the characteristic timescale for the bond-hardening mechanism in O 2 + and explains the onset of bond stabilization for sufficiently long pulse durations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Performance of plane wavefront Fizeau interferometers in power spectral density measurements with tilted plane optics

Binary pseudo-random array (BPRA) test samples are useful devices for calibrating the instrument transfer function (ITF) of Fizeau interferometers, interferometric microscopes, and other optical and non-optical surface and wavefront metrology instruments. The intrinsic white noise character of the power spectral density (PSD) function of the BPRA pattern simplifies the extraction of the ITF from the measured PSD. The ITF determined in a dedicated calibration experiment can be used to reconstruct the surface height profile from the measured data, effectively enhancing the instrument’s spatial resolution. For a high confidence reconstruction procedure, a reliable analytical model of the IFT is desirable. Usually, the model accounts for the contributions to the ITF related the imperfections of the instrument’s optical and detector systems. Here, we experimentally demonstrate that in the case of surface height metrology with Fizeau interferometers, the PSD measurements and, therefore, the efficacy of the ITF calibration of the tool, are strongly affected by the instrument data acquisition and processing procedures, as well as by the shape of the optic under test and its alignment with respect to the interferometer.

Yashchuk, V V↗
DOE OSTI2024

The effect of monitoring complexity on stakeholder acceptance of CO2 geological storage projects in the US gulf coast region

Environmental monitoring at geologic CO 2 storage sites is required by regulations for the purposes of environmental protection and emissions accounting in the case of leakage to surface. However, another very important goal of environmental monitoring is to assure stakeholders that the project is monitored for safety and effectiveness. With current efforts to optimize monitoring for cost-effectiveness, the question remains: will optimization of monitoring approaches degrade stakeholder assurance, or do heavily-instrumented sites communicate higher risk to a stakeholder? We report the results of a stakeholder survey in Gulf Coast states of the US where carbon capture and storage (CCS) is developing quickly. We rely on a 2 by 2 factorial experiment in which we manipulate message complexity (complex v. simple) and social norm (support from scientists v. support from community members). Subjects were randomly assigned to one of four conditions: 1) complex message with scientist support; 2) complex message with community member support; 3) simple message with scientist support; or 4) simple message with community member support. In addition to the experimental stimuli, subjects were also asked about their need for cognition, attitudes toward science and scientists, attitudes about climate change and support for carbon capture and storage (CCS). Our sample is drawn from residents in states bordering the western Gulf of Mexico (Texas, Louisiana, Florida) where CO 2 geologic storage is being planned both onshore and offshore. The results offer important implications for public outreach efforts to key stakeholders.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Text Mining for Process–Structure–Properties Relationships in Metals

With the advent of large language models (LLMs), the vast unstructured text within millions of academic papers is increasingly accessible for materials discovery—although significant challenges remain. While LLMs offer promising few- and zero-shot learning capabilities, particularly valuable in the materials domain where expert annotations are scarce, general-purpose LLMs often fail to address key materials-specific queries without further adaptation. To bridge this gap, fine-tuning LLMs on human-labeled data is essential for effective structured knowledge extraction (Liu in The Importance of Human-Labeled Data in the Era of LLMs, 2023). Here, in this study, we introduce a novel annotation schema designed to extract generic process–structure–properties relationships from scientific literature. We demonstrate the utility of this approach using a dataset of 128 abstracts, with annotations drawn from two distinct domains: high-temperature materials (Domain I) and uncertainty quantification in simulating materials microstructure (Domain II). Initially, we developed a conditional random field (CRF) model based on MatBERT—a domain-specific BERT variant—and evaluated its performance on Domain I. Subsequently, we compared this model with a fine-tuned LLM (GPT-4o from OpenAI) under identical conditions. Our results indicate that fine-tuning LLMs can significantly improve entity extraction performance over the BERT-CRF baseline on Domain I. However, when additional examples from Domain II were incorporated, the performance of the BERT-CRF model became comparable to that of the GPT-4o model. These findings underscore the potential of our schema for structured knowledge extraction and highlight the complementary strengths of both modeling approaches.

Materials science↗
DOE OSTI2024

Analyzing the impact of design factors on solar module thermomechanical durability using interpretable machine learning techniques

Solar modules in utility-scale systems are expected to maintain decades of lifetime to rival conventional energy sources. However, cyclic thermomechanical loading often degrades their long-term performance, highlighting the importance of effective design to mitigate thermal expansion mismatches between module materials. Given the complex composition of solar modules, isolating the impact of individual components on overall durability remains a challenging task. In this work, we analyze a comprehensive data set that comprises bill-of-materials (BOM) and thermal cycling power loss from 251 distinct module designs to identify the predominant design factors and their impacts on the thermomechanical durability of modules. The methodology of our analysis combines machine learning modeling (random forest) and Shapley additive explanation (SHAP) to correlate design factors with power loss and interpret the model’s decision-making. The interpretation reveals that silicon type (monocrystalline or polycrystalline), encapsulant thickness, busbar numbers, and wafer thickness predominantly influence the degradation. With lower power loss of around 0.6% on average in the SHAP analysis, monocrystalline cells present better durability than polycrystalline cells. This finding is further substantiated by statistical testing on our raw data set. The SHAP analysis also demonstrates that while thicker encapsulants lead to reduced power loss, further increasing their thickness over around 0.6 to 0.7 mm does not yield additional benefits, particularly for the front side one. In addition, other important BOM features such as the number of busbars are analyzed. This study provides a blueprint for utilizing explainable machine learning techniques in a complex material system and can potentially guide future research on optimizing the design of solar modules.

14 SOLAR ENERGY↗
DOE OSTI2023

MALDI-MS screening of microbial colonies with isomer resolution to select fatty acid desaturase variants

Creating controlled lipid unsaturation locations in oleochemicals can be a key to many bioengineered products. However, evaluating the effects of modifications to the acyl-ACP desaturase on lipid unsaturation is not currently amenable to high-throughput assays, limiting the scale of redesign efforts to <200 variants. Here, we report a rapid MS assay for profiling the positions of double bonds on membrane lipids produced by Escherichia coli colonies after treatment with ozone gas. By MS measurement of the ozonolysis products of Δ6 and Δ8 isomers of membrane lipids from colonies expressing recombinant Thunbergia alata desaturase, we screened a randomly mutagenized library of the desaturase gene at 5 s per sample. Two variants with altered regiospecificity were isolated, indicated by an increase in 16:1 Δ8 proportion. We also demonstrated the ability of these desaturase variants to influence the membrane composition and fatty acid distribution of E. coli strains deficient in the native acyl-ACP desaturase gene, fabA. Lastly, we used the fabA deficient chassis to concomitantly express a non-native acyl- ACP desaturase and a medium-chain thioesterase from Umbellularia californica, demonstrating production of only saturated free fatty acids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Using soil library hyperspectral reflectance and machine learning to predict soil organic carbon: Assessing potential of airborne and spaceborne optical soil sensing

Soil organic carbon (SOC) is a key variable to determine soil functioning, ecosystem services, and global carbon cycles. Spectroscopy, particularly optical hyperspectral reflectance coupled with machine learning, can provide rapid, efficient, and cost-effective quantification of SOC. However, how to exploit soil hyperspectral reflectance to predict SOC concentration, and the potential performance of airborne and satellite data for predicting surface SOC at large scales remain relatively underknown. Here, this study utilized a continental-scale soil laboratory spectral library (37,540 full-pedon 350–2500 nm reflectance spectra with SOC concentration of 0–780 g·kg –1 across the US) to thoroughly evaluate seven machine learning algorithms including Partial-Least Squares Regression (PLSR), Random Forest (RF), K-Nearest Neighbors (KNN), Ridge, Artificial Neural Networks (ANN), Convolutional Neural Networks (CNN), and Long Short-Term Memory (LSTM) along with four preprocessed spectra, i.e. original, vector normalization, continuum removal, and first-order derivative, to quantify SOC concentration. Furthermore, by using the coupled soil-vegetation-atmosphere radiative transfer model, we simulated twelve airborne and spaceborne hyper/multi-spectral remote sensing data from surface bare soil laboratory spectra to evaluate their potential for estimating SOC concentration of surface bare soils. Results show that LSTM achieved best predictive performance of quantifying SOC concentration for the whole data sets (R 2 = 0.96, RMSE = 30.81 g·kg –1 ), mineral soils (SOC ≤ 120 g·kg –1 , R 2 = 0.71, RMSE = 10.60 g·kg –1 ), and organic soils (SOC > 120 g·kg –1 , R 2 = 0.78, RMSE = 62.31 g·kg –1 ). Spectral data preprocessing, particularly the first-order derivative, improved the performance of PLSR, RF, Ridge, KNN, and ANN, but not LSTM or CNN. We found that the SOC models of mineral and organic soils should be distinguished given their distinct spectral signatures. Finally, we identified that the shortwave infrared is vital for airborne and spaceborne hyperspectral sensors to monitor surface SOC. This study highlights the high accuracy of LSTM with hyperspectral/multispectral data to mitigate a certain level of noise (soil moisture <0.4 m 3 ·m –3 , green leaf area < 0.3 m 2 ·m –2 , plant residue <0.4 m 2 ·m –2 ) for quantifying surface SOC concentration. Forthcoming satellite hyperspectral missions like Surface Biology and Geology (SBG) have a high potential for future global soil carbon monitoring, while high-resolution satellite multispectral fusion data can be an alternative.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Self-biased magnetoelectric switching at room temperature in three-phase ferroelectric–antiferromagnetic–ferrimagnetic nanocomposites

Magnetoelectric systems could be used to develop magnetoelectric random access memory and microsensor devices. One promising system is the two-phase 3-1-type multiferroic nanocomposite in which a one-dimensional magnetic column is embedded in a three-dimensional ferroelectric matrix. However, it suffers from a number of limitations including unwanted leakage currents and the need for biasing with a magnetic field. Here we show that the addition of an antiferromagnet to a 3-1-type multiferroic nanocomposite can lead to a large, self-biased magnetoelectric effect at room temperature. Our three-phase system is composed of a ferroelectric Na 0.5 Bi 0.5 TiO 3 matrix in which ferrimagnetic NiFe 2 O 4 nanocolumns coated with antiferromagnetic p-type NiO are embedded. This system, which is self-assembled, exhibits a magnetoelectric coefficient of up to 1.38 × 10 –9 s m –1 , which is large enough to switch the magnetic anisotropy from the easy axis (K eff = 0.91 × 10 4 J m –3 ) to the easy plane (K eff = –1.65 × 10 4 J m –3 ).

42 ENGINEERING↗
DOE OSTI2025

Anatomy of information scrambling and decoherence in the integrable Sachdev-Ye-Kitaev model

The growth of information scrambling, captured by out-of-time-order correlation functions (OTOCs), is a central indicator of the nature of many-body quantum dynamics. Here, we compute analytically the complete time dependence of the OTOC for an integrable Sachdev-Ye-Kitaev (SYK) model, 𝑁 Majoranas with random two-body interactions of infinite range, coupled to a Markovian bath at finite temperature. In the limit of no coupling to the bath, the time evolution of scrambling experiences different stages. For 𝑡≲√𝑁, after an initial polynomial growth, the OTOC approaches saturation in a power-law fashion with oscillations superimposed. At 𝑡∼√𝑁, the OTOC reverses trend and starts to decrease linearly in time. The reason for this linear decrease is that, despite being a subleading 1/𝑁 effect, the OTOC in this region is governed by the spectral form factor of the antisymmetric couplings of the SYK model. The linear decrease stops at 𝑡∼2⁢𝑁, the Heisenberg time, where saturation occurs. The effect of the environment is an overall exponential decay of the OTOC for times longer than the inverse of the coupling strength to the bath. Here, the oscillations at 𝑡≲√𝑁 indicate lack of thermalization—a desired feature for better performance of quantum information devices.

Open quantum systems & decoherence↗
DOE OSTI2022

Avalanches and many-body resonances in many-body localized systems

Here we numerically study both the avalanche instability and many-body resonances in strongly disordered spin chains exhibiting many-body localization (MBL). Finite-size systems behave like MBL within the MBL regimes, which we divide into the asymptotic MBL phase and the finite-size MBL regime; the latter regime is, however, thermal in the limit of large systems and long times. In both Floquet and Hamiltonian models, we identify some landmarks within the MBL regimes. Our first landmark is an estimate of where the MBL phase becomes unstable to avalanches, obtained by measuring the slowest relaxation rate of a finite chain coupled to an infinite bath at one end. Our estimates indicate that the actual MBL-to-thermal phase transition occurs much deeper in the MBL regimes than has been suggested by most previous studies. Our other landmarks involve systemwide many-body resonances: We find that the effective matrix elements producing eigenstates with systemwide many-body resonances are enormously broadly distributed. This broad distribution means that the onset of such resonances in typical samples occurs quite deep in the MBL regimes, and the first such resonances typically involve rare pairs of eigenstates that are farther apart in energy than the minimum gap. Thus we find that the resonance properties define two landmarks that divide the MBL regimes of finite-size systems into three subregimes: (i) at strongest randomness, typical samples do not have any eigenstates that are involved in systemwide many-body resonances; (ii) there is a substantial intermediate subregime where typical samples do have such resonances but the pair of eigenstates with the minimum spectral gap does not, so the size of the minimum gap agrees with expectations from Poisson statistics; and (iii) in the weaker randomness subregime, the minimum gap is larger than predicted by Poisson level statistics because it is involved in a many-body resonance and thus subject to level repulsion. Nevertheless, even in this third subregime, all but a vanishing fraction of eigenstates remain nonresonant and the system thus still appears MBL in most respects. Based on our estimates of the location of the avalanche instability, it might be that the MBL phase is only part of subregime (i) and the other subregimes are entirely in the thermal phase, even though they look localized in most respects, so are in the finite-size MBL regime.

36 MATERIALS SCIENCE↗