Search NASA⌕ Search

SEARCH · Search NASA

Results for “Fluctuation theorems”

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

Normal modes, rotational inertia, and thermal fluctuations of trapped ion crystals

The normal modes of a trapped ion crystal are derived using an approach based on the Hermitian properties of the system’s dynamical matrix. This method is equivalent to the standard Bogoliubov method, but for classical systems it is arguably simpler and more general in that canonical coordinates are not necessary. The theory is developed for stable, unstable, and neutrally-stable systems. The method is then applied to ion crystals in a Penning trap. Reduced eigenvalue problems for the case of large applied magnetic field are developed, for which the spectrum breaks into E x B drift modes, axial modes, and cyclotron modes. Thermal fluctuation levels in these modes are analyzed and shown to be consistent with the Bohr-van-Leeuwen theorem, provided that neutrally-stable modes associated with crystal rotations are included in the analysis. An expression for the rotational inertia of the crystal is derived, and a magnetic contribution to this inertia, which dominates in large magnetic fields, is described. An unusual limit is discovered for the special case of spherically-symmetric confinement, in which the rotational inertia does not exist and changes in angular momentum leave the rotation frequency unaffected.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Dispersion and the speed-limited particle-in-cell algorithm

This paper discusses temporally continuous and discrete forms of the speed-limited particle-in-cell (SLPIC) method first treated by Werner et al. [Phys. Plasmas 25, 123512 (2018)]. The dispersion relation for a 1D1V electrostatic plasma whose fast particles are speed-limited is derived and analyzed. By examining the normal modes of this dispersion relation, we show that the imposed speed-limiting substantially reduces the frequency of fast electron plasma oscillations while preserving the correct physics of lower-frequency plasma dynamics (e.g. ion acoustic wave dispersion and damping). We then demonstrate how the timestep constraints of conventional electrostatic particle-in-cell methods are relaxed by the speed-limiting approach, thus enabling larger timesteps and faster simulations. Here, these results indicate that the SLPIC method is a fast, accurate, and powerful technique for modeling plasmas wherein electron kinetic behavior is nontrivial (such that a fluid/Boltzmann representation for electrons is inadequate) but evolution is on ion timescales.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

On the symmetries of cosmological perturbations

The space of inflationary models is vast, containing wide varieties of mechanisms, symmetries, and spectra of particles. Consequently, the space of observational signatures is similarly complex. Hence, it is natural to look for boundaries of the space of models and their signatures. In this paper, we explore the possible symmetries associated with the primordial cosmological perturbations and their correlators in the asymptotic future. Assuming the observed homogeneity, isotropy and (approximate) scale invariance, we prove three main results. First, correlation functions of scalar metric fluctuations are uniquely characterized by soft theorems and are free from ambiguity under field redefinitions. Second, whatever the particle content and interactions, when the standard soft theorems apply, invariance under de Sitter boosts (linearly realized conformal invariance) is only possible if all connected correlators vanish identically, i.e. if the theory is free. Third, conformal invariance is the largest set of linearly realized (bosonic) symmetries of the correlators of any single scalar, irrespectively of any soft theorems or particle content.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Influence of inhomogeneous stochasticity on the falsifiability of mean-field theories and examples from accretion disc modelling

ABSTRACT Despite spatial and temporal fluctuations in turbulent astrophysical systems, mean-field theories can be used to describe their secular evolution. However, observations taken over time scales much shorter than dynamical time scales capture a system in a single state of its turbulence ensemble. Comparing with mean-field theory can falsify the latter only if the theory is additionally supplied with a quantified precision. The central limit theorem provides appropriate estimates to the precision only when fluctuations contribute linearly to an observable and with constant coherent scales. Here, we introduce an error propagation formula that relaxes both limitations, allowing for non-linear functional forms of observables and inhomogeneous coherent scales and amplitudes of fluctuations. The method is exemplified in the context of accretion disc theories, where inhomogeneous fluctuations in the surface temperature are propagated to the disc emission spectrum – the latter being a non-linear and non-local function of the former. The derived precision depends non-monotonically on emission frequency. Using the same method, we investigate how binned spectral fluctuations in telescope data change with the spectral resolving power. We discuss the broader implications for falsifiability of a mean-field theory.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Finite-time correlations boost large voltage angle fluctuations in electric power grids

Abstract Decarbonization in the energy sector has been accompanied by an increased penetration of new renewable energy sources in electric power systems. Such sources differ from traditional productions in that, first, they induce larger, undispatchable fluctuations in power generation and second, they lack inertia. Recent measurements have indeed reported long, non-Gaussian tails in the distribution of local voltage frequency data. Large frequency deviations may induce grid instabilities, leading in worst-case scenarios to cascading failures and large-scale blackouts. In this article, we investigate how correlated noise disturbances, characterized by the cumulants of their distribution, propagate through meshed, high-voltage power grids. For a single source of fluctuations, we show that long noise correlation times boost non-Gaussian voltage angle fluctuations so that they propagate similarly to Gaussian fluctuations over the entire network. However, they vanish faster, over short distances if the noise fluctuates rapidly. We furthermore demonstrate that a Berry–Esseen theorem leads to the vanishing of non-Gaussianities as the number of uncorrelated noise sources increases. Our predictions are corroborated by numerical simulations on realistic models of power grids.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

On numerical energy conservation for an implicit particle-in-cell method coupled with a binary Monte-Carlo algorithm for Coulomb collisions

We report conventional particle-in-cell (PIC) methods suffer from enhanced numerical heating (explicit PIC) or cooling (semi-implicit PIC) when coupled with a binary Monte-Carlo algorithm for Coulomb collisions. In this work, a fully-implicit θ-PIC scheme (with adjustable time-biasing parameter 1/2 ≤ θ ≤) is considered. The discrete change in energy of a closed system after a time step for this scheme scales with (1/2 - θ)C θ , where C θ is a positive definite quantity that depends on the frequency spectrum of the energy in the fields. Collisions lead to additional energy in the field fluctuations associated with high-frequency light waves produced by a numerical Bremsstrahlung process, which can result in a large increase in the numerical cooling rate for θ > 1/2. However, for θ = 1/2, energy is exactly conserved. The energy in the field fluctuations on long time scales agrees with that calculated using the equipartition theorem for a classical system in thermodynamic equilibrium.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Kinetic theory of particle-in-cell simulation plasma and the ensemble averaging technique

Abstract We derive the kinetic theory of fluctuations in physically and numerically stable particle-in-cell (PIC) simulations of electrostatic plasmas. The starting point is the single-time correlation at the start of the simulation between the statistical fluctuations of the weighted densities of macroparticle centers in the plasma particle phase-space. The fluctuations are associated with different initial conditions, typically due to the random initial conditions (in velocity space) of the macroparticles/simulation plasma, assigned according to their initial distribution of probability. The single-time correlations at all time steps and in each spatial grid cell are then determined from the Laplace–Fourier transforms of the discretized Klimontovich-like equation for the macroparticles and Maxwell’s equations for the fields, as computed by modern PIC codes. We recover the expressions for the electrostatic field and the plasma particle density fluctuation autocorrelation spectra as well as the kinetic equations describing the average evolution of PIC-simulated plasma particles, first derived by Langdon (1970b Proc. 4th Conf. Numerical Simulation of Plasmas ) using a test macroparticle approach perturbing a discretized Vlasovian plasma and then averaging the obtained physical quantity over the initial macroparticle velocity distribution. We generalize and extend these results to the modern algorithms in PIC codes using arbitrary macroparticle weights. Analytical estimates of statistical fluctuation amplitudes are derived as a function of the plasma simulation parameters, using the central limit theorem in the limit of a large number of macroparticles per cell. The theory is then used to analyze the ensemble averaging technique of PIC simulations where statistical averages are performed over ensembles of PIC simulations, modeling the same plasma physics problem but using different statistical realizations of the initial distribution functions of the macroparticles. This method is illustrated by linear Landau damping uncovering (from noise, which is usually considered numerical) the physical fluctuations driven by a single small amplitude electrostatic wave perturbing a PIC simulation plasma in equilibrium.

fluctuations correlations↗
DOE OSTI2025

Topologically Protected Flatness in Chiral Moiré Heterostructures

The observation of delicate correlated phases in twisted heterostructures of graphene and transition metal dichalcogenides suggests that moiré flat bands are intrinsically resilient against certain types of disorder. Here, we investigate the robustness of moiré flat bands in the chiral limit of the Bistritzer-MacDonald model—applicable to both platforms in certain limits—and demonstrate drastic differences between the first magic angle and higher magic angles in response to chiral symmetric disorder that arise, for instance, from lattice relaxation. We understand these differences using a hidden constant of motion that permits the decomposition of the non-Abelian gauge field induced by interlayer tunnelings into two decoupled Abelian ones. At all magic angles, the resulting effective magnetic field splits into an anomalous contribution and a fluctuating part. The anomalous field maps the moiré flat bands onto a zeroth Dirac Landau level, whose flatness withstands any chiral symmetric perturbation such as nonuniform magnetic fields due to a topological index theorem—thereby underscoring a topological mechanism for band flatness. Only the first magic angle can fully harness this topological protection due to its weak fluctuating magnetic field. In higher magic angles, the amplitude of fluctuations largely exceeds the anomalous contribution, which we find results in a physically meaningless chiral operator and an extremely large sensitivity to microscopic details and an exponential collapse of the single-particle gap. Through numerical simulations, we further study various types of disorder and identify the scattering processes that are enhanced or suppressed in the chiral limit. Interestingly, we find that the topological suppression of disorder broadening persists away from the chiral limit and is further accentuated by isolating a single sublattice polarized flat band in energy. Our analysis suggests the Berry curvature hot spot at the top of the K and K ′ valence band in the transition metal dichalcogenide monolayers is essential for the stability of its moiré flat bands and their correlated states. Published by the American Physical Society 2025

Crépel, Valentin (ORCID:0000000302403412)↗
DOE OSTI2022

Quasilinear theory for inhomogeneous plasma

Here, this paper presents quasilinear theory (QLT) for a classical plasma interacting with inhomogeneous turbulence. The particle Hamiltonian is kept general; for example, relativistic, electromagnetic and gravitational effects are subsumed. A Fokker–Planck equation for the dressed ‘oscillation-centre’ distribution is derived from the Klimontovich equation and captures quasilinear diffusion, interaction with the background fields and ponderomotive effects simultaneously. The local diffusion coefficient is manifestly positive-semidefinite. Waves are allowed to be off-shell (i.e. not constrained by a dispersion relation), and a collision integral of the Balescu–Lenard type emerges in a form that is not restricted to any particular Hamiltonian. This operator conserves particles, momentum and energy, and it also satisfies the H -theorem, as usual. As a spin-off, a general expression for the spectrum of microscopic fluctuations is derived. For on-shell waves, which satisfy a quasilinear wave-kinetic equation, the theory conserves the momentum and energy of the wave–plasma system. The action of non-resonant waves is also conserved, unlike in the standard version of QLT. Dewar's oscillation-centre QLT of electrostatic turbulence is proven formally as a particular case and given a concise formulation. Also discussed as examples are relativistic electromagnetic and gravitational interactions, and QLT for gravitational waves is proposed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Quantifying unitary flow efficiency and entanglement for many-body localization

Holography tells us that there should be a mapping between a boundary Hamiltonian and bulk degrees of freedom. We probe the bulk geometry of the Wegner Wilson Flow (WWF) in the context of many-body localization by addressing efficiency and bulk entanglement growth measures through approximating upper bounds on the boundary entanglement entropy. Here, we connect these upper bounds to the Fubini-Study metric and clarify how a central quantity, the information fluctuation complexity, distinguishes bulk unitary rotation from entanglement production. We also give a short new proof of the small incremental entangling theorem in the absence of ancillas, achieving a dimension independent, universal factor of c = 2.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Causal horizons, geodesic completeness and stability in slow contraction cosmology

We show that cosmological models with a semi-infinite phase of slow contraction (ekpyrosis) possess a combination of properties that can address several fundamental problems in cosmology, otherwise faced in contracting de Sitter phases or standard big bang expansion. In particular, flat or open slow contraction admits a stable past attractor that asymptotes to Minkowski space and is past geodesically complete, as well as a stable, flat, homogeneous, and isotropic future attractor with negligible Weyl curvature (and, therefore, negligible gravitational entropy). In bouncing cosmologies, this contracting attractor is terminated by a smooth, non-singular bounce that transforms the attractor properties at the end of contraction into the initial conditions for the subsequent expanding phase. Cosmologies incorporating a slow contraction phase have no particle horizon and therefore avoid the causal horizon problem. The past Minkowski attractor also generates an initial spectrum of vacuum-like quantum fluctuations on all wavelengths. Moreover, because the averaged expansion rate along past-directed geodesics is non-positive, models incorporating a semi-infinite phase of slow contraction also evade the Borde–Guth–Vilenkin theorem. By contrast, contracting de Sitter space possesses a finite particle horizon and becomes unstable in the presence of scalar fields, matter, or radiation.

Cosmology and Nongalactic Astrophysics (astro-ph.C↗
DOE OSTI2020

The vecro hypothesis

We consider three fundamental issues in quantum gravity: (a) the black hole information paradox (b) the unboundedness of entropy that can be stored inside a black hole horizon (c) the relation between the black hole horizon and the cosmological horizon. With help from the small corrections theorem, we convert each of these issues into a sharp conflict. We then argue that all three conflicts can be resolved by the following hypothesis: the vacuum wave functional of quantum gravity contains a “vecro” component made of virtual fluctuations of configurations of the same type that arise in the fuzzball structure of black hole microstates. Further, if we assume that causality holds to leading order in gently curved spacetime, then we must have such a vecro component in order to resolve the above conflicts. The term vecro stands for “Virtual Extended Compression-Resistant Object”, and characterizes the nature of the vacuum fluctuations that resolve the puzzles. It is interesting that puzzle (c) may relate the role of quantum gravity in black holes to observations in the sky.

Astronomy & Astrophysics↗
DOE OSTI2022

Validating and optimizing mismatch tolerance of Doppler backscattering measurements with the beam model (invited)

In this work, we use the beam model of Doppler backscattering (DBS), which was previously derived from beam tracing and the reciprocity theorem, to shed light on mismatch attenuation. This attenuation of the backscattered signal occurs when the wavevector of the probe beam’s electric field is not in the plane perpendicular to the magnetic field. Correcting for this effect is important for determining the amplitude of the actual density fluctuations. Previous preliminary comparisons between the model and Mega-Ampere Spherical Tokamak (MAST) plasmas were promising. In this work, we quantitatively account for this effect on DIII-D, a conventional tokamak. We compare the predicted and measured mismatch attenuation in various DIII-D, MAST, and MAST-U plasmas, showing that the beam model is applicable in a wide variety of situations. Finally, we performed a preliminary parameter sweep and found that the mismatch tolerance can be improved by optimizing the probe beam’s width and curvature at launch. This is potentially a design consideration for new DBS systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Quantum-gravitational null Raychaudhuri equation

Abstract We consider a congruence of null geodesics in the presence of a quantized spacetime metric. The coupling to a quantum metric induces fluctuations in the congruence; we calculate the change in the area of a pencil of geodesics induced by such fluctuations. For the gravitational field in its vacuum state, we find that quantum gravity contributes a correction to the null Raychaudhuri equation which is of the same sign as the classical terms. We thus derive a quantum-gravitational focusing theorem valid for linearized quantum gravity.

Physics↗
DOE OSTI2021

Stochastic density functional theory: Real- and energy-space fragmentation for noise reduction

Stochastic density functional theory (sDFT) is becoming a valuable tool for studying ground-state properties of extended materials. The computational complexity of describing the Kohn-Sham orbitals is replaced by introducing a set of random (stochastic) orbitals leading to linear and often sub-linear scaling of certain ground-state observables at the account of introducing a statistical error. Schemes to reduce the noise are essential, for example, for determining the structure using the forces obtained from sDFT. Recently, we have introduced two embedding schemes to mitigate the statistical fluctuations in the electron density and resultant forces on the nuclei. Both techniques were based on fragmenting the system either in real space or slicing the occupied space into energy windows, allowing for a significant reduction in the statistical fluctuations. For chemical accuracy, further reduction of the noise is required, which could be achieved by increasing the number of stochastic orbitals. However, the convergence is relatively slow as the statistical error scales as 1 / N χ according to the central limit theorem, where N χ is the number of random orbitals. In this paper, we combined the embedding schemes mentioned above and introduced a new approach that builds on overlapped fragments and energy windows. The new approach significantly lowers the noise for ground-state properties, such as the electron density, total energy, and forces on the nuclei, as demonstrated for a G-center in bulk silicon.

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

Evaluation of two-particle properties within finite-temperature self-consistent one-particle Green’s function methods: Theory and application to GW and GF2

One-particle Green’s function methods can model molecular and solid spectra at zero or non-zero temperatures. One-particle Green’s functions directly provide electronic energies and one-particle properties, such as dipole moment. However, the evaluation of two-particle properties, such as $\langle$S 2 $\rangle$ and $\langle$N 2 $\rangle$, can be challenging because they require a solution of the computationally expensive Bethe–Salpeter equation to find two-particle Green’s functions. We demonstrate that the solution of the Bethe–Salpeter equation can be completely avoided. Applying the thermodynamic Hellmann–Feynman theorem to self-consistent one-particle Green’s function methods, we derive expressions for two-particle density matrices in a general case and provide explicit expressions for GF2 and GW methods. Such density matrices can be decomposed into an antisymmetrized product of correlated one-electron density matrices and the two-particle electronic cumulant of the density matrix. Cumulant expressions reveal a deviation from ensemble representability for GW, explaining its known deficiencies. We analyze the temperature dependence of $\langle$S 2 $\rangle$ and $\langle$N 2 $\rangle$ for a set of small closed-shell systems. Interestingly, both GF2 and GW show a non-zero spin contamination and a non-zero fluctuation of the number of particles for closed-shell systems at the zero-temperature limit.

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

Electromagnetic energy transport by tearing fluctuations in a self-organized reversed-field pinch plasma

Fluctuation measurements reveal the outward electromagnetic energy flux needed to drive the dynamo electromotive force supporting magnetic self-organization in a reversed-field pinch plasma. The radial Poynting flux due to tearing mode fluctuations is measured with an insertable probe during magnetic relaxation. This flux corresponds to transient power levels much larger than the input power and comparable to the global equilibrium magnetic energy transient loss rate. The probe measurements of this flux are roughly as predicted by a simple Poynting's theorem model upon substitution of equilibrium measurement data.

Physics↗
DOE OSTI2025

Lattice QCD Study of Pion Electroproduction and Weak Production from a Nucleon

Quantum fluctuations in QCD influence nucleon structure and interactions, with pion production serving as a key probe of chiral dynamics. In this Letter, we present a lattice QCD calculation of multipole amplitudes at threshold, related to both pion electroproduction and weak production from a nucleon, using two gauge ensembles near the physical pion mass. We develop a technique for spin projection and construct multiple operators for analyzing the generalized eigenvalue problem in both the nucleon-pion system in the center-of-mass frame and the nucleon system with nonzero momentum. The numerical lattice results are then compared with those extracted from experimental data and predicted by low-energy theorems incorporating one-loop corrections. Published by the American Physical Society 2025

Gao, Yu-Sheng (ORCID:0009000406829247)↗