Search NASA⌕ Search

SEARCH · Search NASA

Results for “ergodicity”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Regional-Scale 3D Ground-Motion Simulations of M w 7 Earthquakes on the Hayward Fault, Northern California Resolving Frequencies 0–10 Hz and Including Site-Response Corrections

Large earthquake ground-motion simulations in 3D Earth models provide constraints on site-specific shaking intensities but have suffered from limited frequency resolution and ignored site response in soft soils. In this work we report new regional-scale 3D simulations for moment magnitude 7.0 scenario earthquakes on the Hayward Fault, northern California with SW4. Simulations resolved significantly broader band frequencies (0–10 Hz) than previous studies and represent the highest resolution simulations for any such earthquake to date. Seismic waves were excited by a kinematic rupture following Graves and Pitarka (2016) and obeyed wave propagation in a 3D Earth model with topography from the U.S. Geological Survey (USGS) assuming a minimum shear wavespeed, VSmin, of 500 m/s. We corrected motions for linear and nonlinear site response for the shear wavespeed, VS, from the USGS 3D model, using a recently developed ground-motion model (GMM) for Fourier amplitude spectra . At soft soil locations subjected to strong shaking, the site-corrected intensities reflect the competing effects of linear amplification by low V S material, reduction of stiffness during nonlinear deformation, and damping of high frequencies. Sites with near-surface V S of 500 m/s or greater require no linear site correction but can experience amplitude reduction due to nonlinear response. Averaged over all sites, we obtained reasonable agreement with empirical ergodic median GMMs currently used for seismic hazard and design ground motions (epsilon less than 1), with marked improvement at soft sedimentary sites. At specific locations, the simulated shaking intensities show systematic differences from the GMMs that reveal path and site effects not captured in these ergodic models. Results suggest how next generation regional-scale earthquake simulations can provide higher spatial and frequency resolution while including effects of soft soils that are commonly ignored in scenario earthquake ground-motion simulations.

58 GEOSCIENCES↗
DOE OSTI2025

Spectral anomalies and broken symmetries in maximally chaotic quantum maps

Spectral statistics such as the level spacing statistics and spectral form factor (SFF) are widely expected to accurately identify “ergodicity,” including the presence of underlying macroscopic symmetries, in generic quantum systems ranging from quantized chaotic maps to interacting many-body systems. By studying various quantizations of maximally chaotic maps that break a discrete classical symmetry upon quantization, we demonstrate that this approach can be misleading and fail to detect macroscopic symmetries. Notably, the same classical map can exhibit signatures of different random matrix symmetry classes in short-range spectral statistics depending on the quantization. While the long-range spectral statistics encoded in the early time ramp of the SFF are more robust and correctly identify macroscopic symmetries in several common quantizations, we also demonstrate analytically and numerically that the presence of Berry-like phases in the quantization leads to spectral anomalies, which break this correspondence. Finally, we provide numerical evidence that long-range spectral rigidity remains directly correlated with ergodicity in the quantum dynamical sense of visiting a complete orthonormal basis.

Shou, Laura [Univ. of Maryland, College Park, MD (↗
NASA NTRS1976

Isomeric structures and structural transformations in small clusters

The temperature dependence of the energies of the isomers of a seven-particle system is studied with a view toward understanding ergodicity problems in Monte Carlo simulations. It is found that the phase space of particles in a cluster is not ergodic at lower temperatures.

Kaelberer, J. B.↗
NASA NTRS1978

A method for determining the weak statistical stationarity of a random process

A method for determining the weak statistical stationarity of a random process is presented. The core of this testing procedure consists of generating an equivalent ensemble which approximates a true ensemble. Formation of an equivalent ensemble is accomplished through segmenting a sufficiently long time history of a random process into equal, finite, and statistically independent sample records. The weak statistical stationarity is ascertained based on the time invariance of the equivalent-ensemble averages. Comparison of these averages with their corresponding time averages over a single sample record leads to a heuristic estimate of the ergodicity of a random process. Specific variance tests are introduced for evaluating the statistical independence of the sample records, the time invariance of the equivalent-ensemble autocorrelations, and the ergodicity. Examination and substantiation of these procedures were conducted utilizing turbulent velocity signals.

Sadeh, W. Z.↗
NASA NTRS1983

Ion trajectories in a space-charge wave on a relativistic electron beam

The motion of an ion in a space-charge wave on a strongly magnetized electron beam is investigated. The motions of the ion perpendicular and parallel to the beam direction are coupled by a nonlinear term in the ion Hamiltonian that is proportional to the wave amplitude, and this coupling causes the motion of the ion to deviate significantly from that of a linear harmonic oscillator in certain resonant regions of phase space. A sequence of canonical transformations is employed to investigate the motion of the ion in these regions. It is determined that wave amplitudes that are too small to trap beam electrons are too small to cause these resonances to overlap. When this overlap does not occur, the motion is found not to be discernibly ergodic in any three-dimensional subspace of the energy hypersurface because there exists a third constant of the motion in addition to the total energy and angular momentum. The numerically integrated ion trajectories are studied using surface-of-section techniques in order to verify these findings. It is found that the third constant of the motion constrains an ion initially trapped in a potential well of the wave to remain trapped in this well. It is concluded that ergodic behavior poses no threat at attempts at collective ion acceleration in space-charge waves on an electron beam.

Russell, D. A.↗
NASA NTRS2003

On the Confidence Limit of Hilbert Spectrum

Confidence limit is a routine requirement for Fourier spectral analysis. But this confidence limit is established based on ergodic theory: For stationary process, temporal average equals the ensemble average. Therefore, one can divide the data into n-sections and treat each section as independent realization. Most natural processes in general, and climate data in particular, are not stationary; therefore, there is a need for the Hilbert Spectral analysis for such processes. Here ergodic theory is no longer applicable. We propose to use various adjustable parameters in the shifting processes of the Empirical Mode Decomposition (EMD) method to obtain an ensemble of Intrinsic Mode Function 0 sets. Based on such an ensemble, we introduce a statistical measure in. a form of confidence limits for the Intrinsic Mode Functions, and consequently, the Hilbert spectra. The criterion of selecting the various adjustable parameters is based on the orthogonality test of the resulting M F sets. Length-of-day data from 1962 to 2001 will be used to illustrate this new approach. Its implication in climate data analysis will also be discussed.

Huang, Norden↗
NASA NTRS2011

CUDA Simulation of Homogeneous, Incompressible Turbulence

We discuss very fast Compute Unified Device Architecture (CUDA) simulations of ideal homogeneous incompressible turbulence based on Fourier models. These models have associated statistical theories that predict that Fourier coefficients of fluid velocity and magnetic fields (if present) are zero-mean random variables. Prior numerical simulations have shown that certain coefficients have a non-zero mean value that can be very large compared to the associated standard deviation. We review the theoretical basis of this "broken ergodicity" as applied to 2-D and 3-D fluid and magnetohydrodynamic simulations of homogeneous turbulence. Our new simulations examine the phenomenon of broken ergodicity through very long time and large grid size runs performed on a state-of-the-art CUDA platform. Results comparing various CUDA hardware configurations and grid sizes are discussed. NS and MHD results are compared.

Morin, Lee↗
NASA NTRS2017

Time and Measurement Days

Questions in data analysis involving the concepts of time and measurement are often pushed into the background or reserved for a philosophical discussion. Some examples are: a) Is causality a consequence of the laws of physics, or can the arrow of time be reversed? b) Can we determine the arrow of time of an event? c) Do we need the continuum hypothesis for the underlying function in any measurement process? d) Can we say anything about the analyticity of the underlying process of an event? e) Would it be valid to model a non-analytical process as function of time? f) What are the implications of all these questions for classical Fourier techniques? However, in the age of big data gathered either from space missions supplying ultra-precise long time series, or e.g. LIGO data from the ground, the moment to bring these questions to the foreground seems arrived. The limitations of our understanding of some fundamental processes is emphasized by the lack of solution for problems open for more than 2 decades, such as the non-detection of solar g-modes, or the modal identification of main sequence stellar pulsators like delta Scuti stars. Flicker noise or 1/f noise, for example, attributed in solar-like stars to granulation, is analyzed mostly only to apply noise reduction techniques, neither considering the classical problem of 1/f noise that was introduced a 100 years ago, nor taking into account ergodic or non-ergodic solutions that make inapplicable spectral analysis techniques in practice. This topic was discussed by Nicholas W. Watkins during the ITISE meeting held in Granada in 2016. There he presented preliminary results of his research on Mandelbrot's related work. We reproduce here his quotation of Mandelbrot (1999) "There is a sharp contrast between a highly anomalous ("non-white") noise that proceeds in ordinary clock time and a noise whose principal anomaly is that it is restricted to fractal time", suggesting a connection with the above proposed topics that could be phrased as the following additional questions:a) Is self-organized criticality (SOC) frequent in astrophysical phenomena? b) Could all fractals in nature be considered stochastic? c) Could we establish mathematical/physical relationships between chaotic and fractal behaviors in time series? d) Could the differences between fractals and chaos in terms of analyticity be used to understand the residuals of the fitting of stellar light curves? In this meeting we would like to approximate these problems from a holistic and multidisciplinary perspective, taking into account not only technical issues but also the deeper implications. In particular the concept of connectivity (introduced in Pascual-Granado et al. A&A, 2015) could be used to implement, within the framework of ARMA processes, an "arrow of time" (see attached document), and so studying the possible implications in the concept of time as envisaged by Watkins.

data analysis↗
DOE OSTI2024

Local Poincaré algebra from quantum chaos

The local two-dimensional Poincaré algebra near the horizon of an eternal AdS black hole, or in proximity to any bifurcate Killing horizon, is generated by the Killing flow and outward null translations on the horizon. In holography, this local Poincaré algebra is reflected as a pair of unitary flows in the boundary Hilbert space whose generators under modular flow grow and decay exponentially with a maximal Lyapunov exponent. This is a universal feature of many geometric vacua of quantum gravity. To explain this universality, we show that a two-dimensional Poincaré algebra emerges in any quantum system that has von Neumann subalgebras associated with half-infinite modular time intervals (modular future and past subalgebras) in a limit analogous to the near-horizon limit. In ergodic theory, quantum dynamical systems with future or past algebras are called quantum K-systems. The surprising statement is that modular K-systems are always maximally chaotic. Interacting quantum systems in the thermodynamic limit and large N theories above the Hawking-Page phase transition are examples of physical theories with future/past subalgebras. We prove that the existence of (modular) future/past von Neumann subalgebras also implies a second law of (modular) thermodynamics and the exponential decay of (modular) correlators. We generalize our results from the modular flow to any dynamical flow with a positive generator and interpret the positivity condition as quantum detailed balance.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Constraining chaos: Enforcing dynamical invariants in the training of reservoir computers

Drawing on ergodic theory, we introduce a novel training method for machine learning based forecasting methods for chaotic dynamical systems. The training enforces dynamical invariants—such as the Lyapunov exponent spectrum and the fractal dimension—in the systems of interest, enabling longer and more stable forecasts when operating with limited data. The technique is demonstrated in detail using reservoir computing, a specific kind of recurrent neural network. Finally, results are given for the Lorenz 1996 chaotic dynamical system and a spectral quasi-geostrophic model of the atmosphere, both typical test cases for numerical weather prediction.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Perturbative instability of nonergodic phases in non-Abelian quantum chains

An important challenge in the field of many-body quantum dynamics is to identify nonergodic states of matter beyond many-body localization (MBL). Strongly disordered spin chains with non-Abelian symmetry and chains of non-Abelian anyons are natural candidates, as they are incompatible with standard MBL. In such chains, real space renormalization group methods predict a partially localized, nonergodic regime known as a quantum critical glass (a critical variant of MBL). Here, this regime features a treelike hierarchy of integrals of motion and symmetric eigenstates with entanglement entropy that scales as a logarithmically enhanced area law. We argue that such tentative nonergodic states are perturbatively unstable using an analytic computation of the scaling of off-diagonal matrix elements and accessible level spacing of local perturbations. Our results indicate that strongly disordered chains with non-Abelian symmetry display either spontaneous symmetry breaking or ergodic thermal behavior at long times. We identify the relevant length and timescales for thermalization: Even if such chains eventually thermalize, they can exhibit nonergodic dynamics up to parametrically long timescales with a nonanalytic dependence on disorder strength.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Equilibration and “Thermalization” in the Adapted Caldeira–Leggett Model

I explore the processes of equilibration exhibited by the Adapted Caldeira–Leggett (ACL) model, a small unitary “toy model” developed for numerical studies of quantum decoherence between an SHO and an environment. I demonstrate how dephasing allows equilibration to occur in a wide variety of situations. While the finite model size and other “unphysical” aspects prevent the notions of temperature and thermalization from being generally applicable, certain primitive aspects of thermalization can be realized for particular parameter values. I link the observed behaviors to intrinsic properties of the global energy eigenstates, and argue that the phenomena I observe contain elements which might be key ingredients that lead to ergodic behavior in larger more realistic systems. The motivations for this work range from curiosity about phenomena observed in earlier calculations with the ACL model to much larger questions related to the nature of equilibrium, thermalization, and the emergence of physical laws.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Self‐consistent transport simulation of boron dust particle injection in the peripheral plasma in Large Helical Device

Abstract The trajectories and the ablation positions of boron dust particles dropped from an impurity powder dropper in the peripheral plasma in the Large Helical Device (LHD) were calculated using a three‐dimensional edge plasma simulation code (EMC3‐EIRENE) and a dust transport simulation code (DUSTT). The simulation shows that the trajectory of the boron dust particles is deflected at the upper divertor leg due to the effect of the hydrogen plasma flow, and the ablation positions of the dust particles in an ergodic layer change toward the outboard side of the torus for higher plasma densities. The effect of the boron ion flow in the divertor leg on the deflection is investigated by coupling the two codes self‐consistently. The simulation predicts that the boron ions in the divertor leg, which are produced by sputtering on the divertor plates, which do not affect the change in the ablation positions. It also shows that the ablation positions move toward the inboard side and approach the Last Closed Flux Surface (LCFS) in case of increased boron dust drop rates, which is caused by the lowered plasma flow in the upper divertor leg due to the lowered electron temperature by radiation cooling by the dropped dust particles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Effect of Magnetic Islands on Neoclassical Heat Diffusivity in a Global Tokamak Simulation

Resonant magnetic perturbations (RMPs) alter magnetic field topology via island formation and can modify core transport in toroidal fusion devices. Using the global gyrokinetic particle‐in‐cell code XGC‐S, originally developed for stellarator geometries, we quantify how island topology affects neoclassical heat transport in a circular tokamak with ($m,n$) = (5,2) RMP‐induced islands. Electron radial heat diffusivity exhibits a two‐peak structure with a dominant enhancement near the O/X‐point region and a secondary peak at the outer island boundary, whereas ion transport remains close to the neoclassical prediction. The observed electron diffusivity is strongly enhanced over the neoclassical baseline but remains well below the Rechester–Rosenbluth estimate, consistent with the non‐ergodic character of the island topology. The density dependence of diffusivity is weaker inside the islands than outside, indicating the presence of topology‐driven transport channels that do not depend strongly on collisionality. Two‐dimensional heat flux maps reveal an up–down antisymmetric structure at the island separatrices consistent with magnetic gradient and curvature drift effects, in contrast to the symmetric temperature flattening inside the islands. These results identify topology‐dependent avenues for electron heat transport in RMP‐perturbed cores and suggest regimes in which island geometry, rather than collisions alone, governs radial transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Spectral form factor of a quantum spin glass

It is widely expected that systems which fully thermalize are chaotic in the sense of exhibiting random-matrix statistics of their energy level spacings, whereas integrable systems exhibit Poissonian statistics. In this paper, we investigate a third class: spin glasses. These systems are partially chaotic but do not achieve full thermalization due to large free energy barriers. We examine the level spacing statistics of a canonical infinite-range quantum spin glass, the quantum p-spherical model, using an analytic path integral approach. We find statistics consistent with a direct sum of independent random matrices, and show that the number of such matrices is equal to the number of distinct metastable configurations — the exponential of the spin glass “complexity” as obtained from the quantum Thouless-Anderson-Palmer equations. We also consider the statistical properties of the complexity itself and identify a set of contributions to the path integral which suggest a Poissonian distribution for the number of metastable configurations. Our results show that level spacing statistics can probe the ergodicity-breaking in quantum spin glasses and provide a way to generalize the notion of spin glass complexity beyond models with a semi-classical limit.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Sub-Cloud Turbulence Explains Cloud-Base Updrafts for Shallow Cumulus Ensembles: First Observational Evidence

Sub-cloud turbulent kinetic energy has been used to parameterize the cloud-base updraft velocity (w b ) in cumulus parameterizations. The validity of this idea has never been proved in observations. Instead, it was challenged by recent Doppler Lidar observations showing a poor correlation between the two. We argue that the low correlation is likely caused by the difficulty of a fixed-point Lidar to measure ensemble properties of cumulus fields. Taking advantage of the stationarity and ergodicity of early afternoon convection, we developed a Lidar sampling methodology to measure w b of a shallow cumulus (ShCu) ensemble (not a single ShCu). By analyzing 128 ShCu ensembles over the Southern Great Plains, we show that the ensemble properties of sub-cloud turbulence explain nearly half of the variability in ensemble-mean wb, demonstrating the ability of sub-cloud turbulence to dictate w b . Lastly, the derived empirical formulas will be useful for developing cumulus parameterizations and satellite inference of w b .

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Realization of two-dimensional discrete time crystals with anisotropic Heisenberg coupling

A discrete time crystal (DTC) is an out-of-equilibrium phase of matter that spontaneously breaks discrete time-translation symmetry. Previous studies have been limited to a set of models with Ising-like couplings - and mostly only in one dimension - thus precluding our understanding of the existence (or not) of DTCs in models with more realistic interactions. In this work, by combining the latest generation of IBM quantum processors with state-of-the-art tensor network methods, we demonstrate the existence of a DTC in a two-dimensional system governed by anisotropic Heisenberg interactions. We uncover a rich phase diagram encompassing spin-glass, ergodic, and time-crystalline phases, and identify the interplay of initialization, interaction anisotropy, and driving protocols in stabilizing the DTC phase. By extending the study of Floquet matter beyond simplified models, we lay the groundwork for exploring how driven systems bridge the gap between quantum coherence and emergent non-equilibrium thermodynamics.

Phase transitions and critical phenomena↗