Search NASA⌕ Search

SEARCH · Search NASA

Results for “ergodic”

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.

34 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Coherence-Induced Deep Thermalization Transition in Random Permutation Quantum Dynamics

We report a phase transition in the projected ensemble—the collection of postmeasurement wave functions of a local subsystem obtained by measuring its complement. The transition emerges in systems undergoing random permutation dynamics, a type of quantum time evolution wherein computational basis states are shuffled without creating superpositions. It separates a phase exhibiting deep thermalization, where the projected ensemble is distributed over Hilbert space in a maximally entropic fashion (Haar random), from a phase where it is minimally entropic (“classical bit-string ensemble”). Crucially, this deep thermalization transition is invisible to the subsystem’s density matrix, which always exhibits thermalization to infinite temperature across the phase diagram. Through a combination of analytical arguments and numerical simulations, we show that the transition is tuned by the total amount of injected by the input state and the measurement basis, and is exhibited robustly across different microscopic models. Our findings represent a novel form of ergodicity-breaking universality in quantum many-body dynamics, characterized not by a failure of regular thermalization, but rather by a failure of deep thermalization.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Generic Hilbert space fragmentation in Kogut-Susskind lattice gauge theories

At the heart of quantum many-body physics lies the understanding of mechanisms that avoid quantum thermalization in an isolated system quenched far from equilibrium. A prominent example is Hilbert space fragmentation, which has recently emerged as an ergodicity-breaking mechanism in constrained spin models. Here, we show that Kogut-Susskind formulations of lattice gauge theories in d+1D (d spatial and one temporal dimensions) give rise to Hilbert space fragmentation, and discuss possible implications for understanding continuum physics. Lastly, our findings not only prove that lattice gauge theories are a natural platform for Hilbert space fragmentation, they also serve as a guide to the conditions under which these models can be faithfully used to infer the thermalization properties of quantum chromodynamics.

Eigenstate thermalization↗
DOE OSTI2024

Quantum fragmentation in the extended quantum breakdown model

We introduce a one-dimensional (1D) extended quantum breakdown model comprising a fermionic and a spin degree of freedom per site, and featuring a spatially asymmetric breakdown-type interaction between the fermions and spins. Furthermore, our model resembles the breakdown process of particles incident into a cloud chamber with nonzero quantum amplitudes of both exciting and not exciting the local vapor atoms. We analytically show that, in the absence of any magnetic field for the spins, the model exhibits Hilbert space fragmentation within each symmetry sector into exponentially many Krylov subspaces and hence displays nonthermal dynamics. Here, we demonstrate that the fragmentation naturally occurs in an entangled basis and thus provides an example of “quantum fragmentation.” Besides establishing the nature of fragmentation analytically, we also study the long-time behavior of the entanglement entropy and its deviation from the expected Page value as a probe of ergodicity in the system. Upon introducing a magnetic field for the spins, most of the Krylov subspaces merge and the model becomes chaotic. Finally, we study the effects of strong randomness on the system and observe behavior similar to that of many-body localized systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Operator dynamics in Floquet many-body systems

We study operator dynamics in many-body quantum systems, focusing on generic features of systems that are ergodic, spatially extended, and lack conserved densities. Quantum circuits of various types provide simple models for such systems. We focus on Floquet quantum circuits, comparing their behavior with what has been found previously for circuits that are random in time. Floquet circuits, which have discrete time-translation symmetry, represent an intermediate case between circuits that are random in time and lack any symmetry, and systems with a time-independent Hamiltonian and continuous time-translation invariance. By making this comparison, one of our aims is to identify signatures of time-translation symmetry in Floquet operator dynamics. To characterize behavior we examine a variety of quantities in solvable models and numerically: operator autocorrelation functions; the partial spectral form factor; the out-of-time-order correlator (OTOC); and the paths in operator space that make the dominant contributions to the ensemble-averaged autocorrelation functions. Our most striking result is that ensemble-averaged autocorrelation functions show behavior that is distinctively different in Floquet systems compared to systems in which successive time-steps are independent. Specifically, while average autocorrelation functions decay on a microscopic timescale for circuits that are random in time, in Floquet systems they have a late-time tail with a duration that grows parametrically with the size of the operator support. In the simplest models this tail is separated from the initial decay by a minimum, so that the average autocorrelation function has an intermediate-time peak. The existence of these tails provides a way to understand deviations of the spectral form factor from random matrix behavior at times shorter than the Thouless time. In contrast to this feature in autocorrelation functions, we find no new aspects to the behavior of OTOCs for Floquet models compared to random-in-time circuits. We show that this difference between averaged autocorrelation functions and OTOCs can be understood in terms of the paths in operator space that contribute to the two quantities: paths for the former retain a limited support at late times, while paths for the latter are dominated by operator spreading. Published by the American Physical Society 2025

Yoshimura, Takato (ORCID:0000000309159846)↗
DOE OSTI2024

Coordinate versus momentum cuts and effects of collective flow on critical fluctuations

We analyze particle number fluctuations in the crossover region near the critical endpoint of a first-order phase transition by utilizing molecular dynamics simulations of the classical Lennard-Jones fluid. We extend our previous study [V. A. Kuznietsov , ] by incorporating longitudinal collective flow. The scaled variance of particle number distribution inside different coordinate and momentum space acceptances is computed through ensemble averaging and found to agree with earlier results obtained using time averaging, validating the ergodic hypothesis for fluctuation observables. Presence of a sizable collective flow is found to be essential for observing large fluctuations from the critical point in momentum space acceptances. We discuss our findings in the context of heavy-ion collisions. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Emergent hydrodynamic mode on SU(2) plaquette chains and quantum simulation

We search for emergent hydrodynamic modes in real-time Hamiltonian dynamics of 2+1-dimensional SU(2) lattice gauge theory on a quasi-one-dimensional plaquette chain, by numerically computing symmetric correlation functions of energy densities on lattice sizes of about 20 with the local Hilbert space truncated at 𝑗 max = $\frac{1}{2}$. Because of the Umklapp processes, we only find a mode for energy diffusion. The symmetric correlator exhibits transport peak near zero frequency with a width approximately proportional to momentum squared at small momentum, when the system is fully quantum ergodic, as indicated by the eigenenergy level statistics. This transport peak leads to a power-law 𝑡 −$\frac{1}{2}$ decay of the symmetric correlator at late time, also known as the long-time tail, as well as diffusionlike spreading in position space. We also introduce a quantum algorithm for computing the symmetric correlator on a quantum computer and find it gives results consistent with exact diagonalization when tested on the IBM emulator. Finally we discuss the future prospect of searching for the sound modes.

Hamiltonian systems↗
DOE OSTI2026

Scale setting of SU⁡(𝑁) Yang–Mills theory, topology and large-𝑁 volume independence

We set the scale of SU⁡(𝑁) Yang-Mills theories for 𝑁 =3, 5, 8 and in the large-𝑁 limit via gradient flow, as a first step towards the computation of the large-𝑁 Λ-parameter using step scaling. We adopt twisted boundary conditions to achieve large-𝑁 volume reduction and the Parallel Tempering on Boundary Conditions algorithm to tame topological freezing. This setup allows accurate determinations of the gradient-flow scales down to lattice spacings as fine as ∼0.025 fm for all the explored values of 𝑁, a regime that has never been reached with ergodic algorithms. Moreover, we are able to precisely estimate the finite-size systematics related to topological freezing, and to show the suppression of finite-volume effects expected by virtue of large-𝑁 twisted volume reduction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Boundary-induced classical generalized Gibbs ensemble with angular momentum

We investigate how confinement geometry leads to the emergence of a Generalized Gibbs Ensemble (GGE) in classical systems. Unlike the standard Gibbs ensemble, the GGE includes additional conserved quantities, such as angular momentum, that arise from boundary-induced symmetries. Using analytical arguments based on the maximum entropy principle, we show that circular boundaries preserve angular momentum and drive the system toward a chiral, non-ergodic GGE that violates time-reversal symmetry. This ensemble differs fundamentally from the Gibbs case, producing near-boundary condensation and revealing how geometry alone can alter thermal equilibration. To quantify these effects, we introduce an order parameter measuring deviations from Gibbs behavior and demonstrate that conventional Monte Carlo methods must incorporate angular momentum conservation under such conditions. Our study highlights how geometric constraints shape non-equilibrium statistical ensembles and lead to subtle departures from the Bohr-van Leeuwen theorem. These predictions are validated through detailed simulations of confined classical hard-disk gases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Simulation of Physics-Based 0-10Hz Strong Motion Using High Performance Computing Supporting Refinements to Regional Ground Motion Models for the Central Eastern US

In collaboration with the U.S. Nuclear Regulatory Commission (NRC) the LLNL has developed a computationally efficient simulation platform designed to perform physics-based ground motion simulations for crustal earthquakes in the Stable Continental Regions of Central and Eastern US (CEUS), using high-performance computing. The main objective of the earthquake simulations was to use synthetic ground motion to provide constrains to refinements of existing ergodic Ground Motion Models (GMMs), for large magnitude earthquakes and near-fault distances, for which these models are less reliable. Physics-based broadband (0-10Hz) ground motion simulations were used to estimate the near-fault ground motion amplitudes and within event and between-event variabilities associated with fault rupture characteristics. In our simulations we used a 3D regional velocity model that was based on Saikia’s 1D velocity model (1994). In simulations performed during the first stage of this project the Saikia’s velocity model demonstrated better performance in modelling high frequency regional wave propagation for the CEUS region recorded during the Mw5.0 November 7, 2016, Cushing Oklahoma (Taylor et al., 2017), and Mw5.8 September 3, 2016, Pawnee Oklahoma earthquakes. The proposed regional 3D model includes random perturbations to the 1D background model using the stochastic scheme of Pitarka and Mellors (2021). In addition, validation analysis of the rupture generator and regional wave propagation models, using comparisons with different GMMs for Mw6.5 and Mw7.0 scenario earthquakes in the CEUS region resulted in a very good match between the simulated and empirical ground motion models. For the purposes of seismic hazard assessment at the existing and planned nuclear power plants, NRC is interested in studies aimed at improving the current ground motion models (GMM) for both Stable Continental Regions (SCR) in the Central and Eastern US and Active Crustal Regions (ACR) in the Western US. Due to lack of recorded data, these improvements require synthetic data for short fault distances and large magnitude earthquakes for which the existing recorded data is not enough to uniquely constrain the GMMs. The need for simulations and strong motion data is especially critical for the CEUS region where we do not have recorded data from potentially large damaging earthquakes with moment magnitudes 6.0 and higher. In this the project, we focused on 10Hz simulations of Mw7.0 scenario earthquakes with strike slip and thrust faulting mechanisms. We used more than 50 Mw7.0 earthquake rupture scenarios to investigate the ground motion uncertainty due to unknown earthquake rupture parameters, in particular, the slip distribution, rupture velocity, and faulting mechanism, and their implication on ground motion amplification due to forward rupture directivity effects.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Zentropy Theory for Transformative Functionalities of Magnetic and Superconducting Materials

The proposed research developed the zentropy theory through applications to complex magnetic materials and superconductors under the hypothesis that the emergent properties of complex magnetic materials and superconductors can be predicted by statistical mechanics of ergodic microstates with their partition functions computed from DFT-predicted free energies. The key objective is to develop approaches to systematically determine the types and number of microstates and the supercell size in DFT-based calculations through convergency of macroscopic functionalities, with the incorporation of our mixed-space approach accounting for the interactions between periodic supercells. In addition to use scientific intuitions to guide the design of important microstates, the key innovation of the proposed research is to integrate the domain knowledge and the material-property-descriptor database (MPDD) with 4 million microstates, which is supported by our deep neural network machine learning models (SIPFENN: structure-informed prediction of formation energy using neural networks) and integrated with our high throughput DFT Tool Kit (DFTTK). For complex magnetic materials, one of the objectives is to develop approaches to calculate short-range ordering from the statistical distribution of each microstate. For superconductors, the divergency of quasiparticle effective mass at a quantum critical point will be investigated, and the superconducting and non-superconducting microstates will be delineated through analysis of electronic band structure, density of states, charge density, and Fermi surface.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Derivation of A Representative Elementary Volume (REV) for Upscaled Two-Phase Flow in Porous Media

Relative permeability plays an important role in the upscaling of multiphase flow in porous media from the pore scale to the Darcy scale. The entire concept of relative permeability is contingent on the existence of a representative elementary volume (REV). As we move to smaller samples to measure relative permeability, such as with digital core analysis, the concept of a classical REV has become increasingly unlikely when using the conventional approach to defining a representative volume. The “‘conventional”’ understanding of an REV is that a large enough volume must be considered such that spatial variability averages out. In digital rock methods, such as pore-scale simulations based on micro-computed tomography (CT) images, the domain size is typically 2 to 4 mm. This is approximately the length scale of a single-phase flow REV using the classic REV approach. However, the single-phase perspective does not consider the complex dynamics and fluctuations often observed in multiphase flow systems, even at centimeter-scale experiments and/or simulations. A fundamental question is, therefore, whether the domain size commonly used in digital rock simulations can provide a consistent energy budget such that the concept of relative permeability exists. Based on first principles, relative permeability accounts for the rate of energy dissipated in a stationary process. If the dynamics are fluctuating, the energy dissipated can vary but will average out over a long enough timescale. The key to determining the validity of the relative permeability is the timescale of the measurement, not the spatial scale. The conventional REV theory assumes that spatial, temporal, and ensemble averages are equivalent in an ergodic system, but it does not provide a way to test this assumption. We provide a formal way to identify the timescale where the relative permeability accurately captures energy dissipation as a way to validate relative permeability measurements and quantitatively assess their accuracy. This result will be tested for a practical SCAL test, determining how long a flow experiment needs to be run to accurately characterize the rate of energy dissipation by the flow. The outcome will be a best practice guide for the determination of relative permeability from core-scale experiments and/or digital core simulations that ensure the energy budget is fully accounted for in the relative permeability coefficient.

Mcclure, James [Virginia Tech, Blacksburg]↗
DOE OSTI2026

The Maximal Entanglement Limit in Statistical and High-energy Physics

These lectures advocate the idea that quantum entanglement provides a unifying foundation for both statistical physics and high-energy interactions. I argue that, at sufficiently long times or high energies, most quantum systems approach a Maximal Entanglement Limit (MEL) in which phases of quantum states become unobservable, reduced density matrices acquire a thermal form, and probabilistic descriptions emerge without invoking ergodicity or classical randomness. Within this framework, the emergence of probabilistic parton model, thermalization in the break-up of confining strings and in high-energy collisions, and the universal small-x behavior of structure functions arise as direct consequences of entanglement and geometry of high-dimensional Hilbert space.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Nonlocal effects on thermal transport in hydrodynamic simulations of unmagnetized MagLIF-relevant gaspipes on NIF

We present simulations of heat flow relevant to gaspipe experiments on the National Ignition Facility to investigate kinetic effects on transport phenomena. D 2 and neopentane (C 5 H 12 ) filled targets are used to study the laser preheat stage of a MagLIF scheme where an axial magnetic field is sometimes applied to the target. Simulations were done with the radiation-MHD code HYDRA with a collision-dominated fluid model and the SNB nonlocal electron thermal conduction model. Using the SNB model to evolve the electron temperature increased the heat front propagation of neopentane gas targets compared to a local model by limiting radial heat flow. This increases electron temperature near the axis, which decreases laser absorption. We find that the effect of heat flow models on temperature profiles and laser propagation is modest. Beyond the SNB model, we utilize HYDRA to initialize plasma conditions for the Vlasov–Fokker–Planck K2 code. We run K2 until a quasi-steady state is reached and examine the impact of kinetic effects on heat transport. Although axial heat flow is well predicted by fluid models, the fluid model consistently overpredicts radial heat flow up to 150% in regions with the largest temperature gradient of D 2 filled gaspipes. On the other hand, the SNB nonlocal electron conduction model is found to be adequate for capturing kinetic heat flow in gaspipes.

Lau, Ryan Y. [Univ. of Colorado, Boulder, CO (Unit↗
DOE OSTI2026

Spatiotemporal control of laser intensity using differentiable programming

Optical techniques for spatiotemporal control can produce laser pulses with custom amplitude, phase, or polarization structure. In nonlinear optics and plasma physics, the use of structured pulses typically follows a forward design approach, in which the efficacy of a known structure is analyzed for a particular application. Inverse approaches, in contrast, enable the discovery of new structures with the potential for superior performance. Here, an implementation of the unidirectional pulse propagation equation that supports automatic differentiation is combined with gradient-based optimization to design structured pulses with features that are advantageous for a range of nonlinear optical and plasma-based applications: (1) a longitudinally uniform intensity over an extended region, (2) a superluminal intensity peak that travels many Rayleigh ranges with constant duration, spot size, and amplitude, and (3) a laser pulse that ionizes a gas to form a uniform column of plasma. In the final case, optimizing the full spatiotemporal structure improves the performance by a factor of 15 compared to optimizing only spatial or only temporal structure, highlighting the advantage of spatiotemporal control.

automatic differentiation↗
DOE OSTI2025

Single-Shot Reconstruction of Electron Beam Longitudinal Phase Space in a Laser Wakefield Accelerator

We report on a single-shot longitudinal phase-space reconstruction diagnostic for electron beams in a laser wakefield accelerator via the experimental observation of distinct periodic modulations in the angularly resolved spectra. Such modulated angular spectra arise as a result of the direct interaction between the ultrarelativistic electron beam and the laser driver in the presence of the wakefield. A constrained theoretical model for the coupled oscillator, assisted by a genetic algorithm, can recreate the experimental electron spectra and, thus, fully reconstructs the longitudinal phase-space distribution of the electron beam with a temporal resolution of approximately 1.3 fs. In particular, it reveals the slice energy spread of the electron beam, which is important to measure for applications such as x-ray free electron lasers. In our experiment, the root-mean-square energy spread retrieved is bounded at 9.9 MeV, corresponding to a 0.9%–3.0% relative spread, despite the overall GeV energy beam having approximately 100% relative energy spread.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗