Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nonequilibrium statistical mechanics”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS1971

The closure approximation in the hierarchy equations.

The expectation of the solution process in a stochastic operator equation can be obtained from averaged equations only under very special circumstances. Conditions for validity are given and the significance and validity of the approximation in widely used hierarchy methods and the ?self-consistent field' approximation in nonequilibrium statistical mechanics are clarified. The error at any level of the hierarchy can be given and can be avoided by the use of the iterative method.

Adomian, G.↗
DOE OSTI2023

Quantum bath augmented stochastic nonequilibrium atomistic simulations for molecular heat conduction

Classical molecular dynamics (MD) has been shown to be effective in simulating heat conduction in certain molecular junctions since it inherently takes into account some essential methodological components which are lacking in the quantum Landauer-type transport model, such as many-body full force-field interactions, anharmonicity effects and nonlinear responses for large temperature biases. However, the classical MD reaches its limit in the environments where the quantum effects are significant (e.g. with low-temperatures substrates, presence of extremely high frequency molecular modes). Here, we present an atomistic simulation methodology for molecular heat conduction that incorporates the quantum Bose–Einstein statistics into an “effective temperature” in the form of a modified Langevin equation. We show that the results from such a quasi-classical effective temperature MD method deviates drastically when the baths temperature approaches zero from classical MD simulations and the results converge to the classical ones when the bath approaches the high-temperature limit, which makes the method suitable for full temperature range. In addition, we show that our quasi-classical thermal transport method can be used to model the conducting substrate layout and molecular composition (e.g. anharmonicities, high-frequency modes). Anharmonic models are explicitly simulated via the Morse potential and compared to pure harmonic interactions to show the effects of anharmonicities under quantum colored bath setups. Finally, the chain length dependence of heat conduction is examined for one-dimensional polymer chains placed in between quantum augmented baths.

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

Universal energy-speed-accuracy trade-offs in driven nonequilibrium systems

The connection between measure theoretic optimal transport and dissipative nonequilibrium dynamics provides a language for quantifying nonequilibrium control costs, leading to a collection of thermodynamic speed limits, which rely on the assumption that the target probability distribution is perfectly realized. This is almost never the case in experiments or numerical simulations, so here we address the situation in which the external controller is imperfect. We obtain a lower bound for the dissipated work in generic nonequilibrium control problems that (1) is asymptotically tight and (2) matches the thermodynamic speed limit in the case of optimal driving. Along with analytically solvable examples, we refine this imperfect driving notion to systems in which the controlled degrees of freedom are slow relative to the nonequilibrium relaxation rate, and identify independent energy contributions from fast and slow degrees of freedom. Furthermore, we develop a strategy for optimizing minimally dissipative protocols based on optimal transport flow matching, a generative machine learning technique. Furthermore, this latter approach ensures the scalability of both the theoretical and computational framework we put forth. Crucially, we demonstrate that we can compute the terms in our bound numerically using efficient algorithms from the computational optimal transport literature and that the protocols we learn saturate the bound.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Simulations of classical three-body thermalization in one dimension

One-dimensional systems, such as nanowires or electrons moving along strong magnetic field lines, have peculiar thermalization physics. The binary collision of pointlike particles, typically the dominant process for reaching thermal equilibrium in higher-dimensional systems, cannot thermalize a 1D system. We study how dilute classical 1D gases thermalize through three-body collisions. We consider a system of identical classical point particles with pairwise repulsive inverse power-law potential V ij ∝ 1/|x i –x j | n or the pairwise Lennard-Jones potential. Using Monte Carlo methods, we compute a collision kernel and use it in the Boltzmann equation to evolve a perturbed thermal state with temperature T toward equilibrium. We explain the shape of the kernel and its dependence on the system parameters. Additionally, we implement molecular dynamics simulations of a many-body gas and show agreement with the Boltzmann evolution in the low-density limit. For the inverse power-law potential, the rate of thermalization is proportional to ρ 2 ⁢T$\frac{1}{2}$ – $\frac{1}{n}$, where ρ is the number density. Furthermore, the corresponding proportionality constant decreases with increasing n.

1-dimensional systems↗
DOE OSTI2023

Full Counting Statistics of Charge in Chaotic Many-Body Quantum Systems

We investigate the full counting statistics of charge transport in U(1)-symmetric random unitary circuits. We consider an initial mixed state prepared with a chemical potential imbalance between the left and right halves of the system and study the fluctuations of the charge transferred across the central bond in typical circuits. Using an effective replica statistical mechanics model and a mapping onto an emergent classical stochastic process valid at large on-site Hilbert space dimension, we show that charge transfer fluctuations approach those of the symmetric exclusion process at long times, with subleading t –1/2 quantum corrections. Here, we discuss our results in the context of fluctuating hydrodynamics and macroscopic fluctuation theory of classical nonequilibrium systems and check our predictions against direct matrix-product state calculations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Monitored Fluctuating Hydrodynamics

We introduce a hydrodynamic framework for describing monitored classical stochastic processes. We study the conditional ensembles for these monitored processes—i.e., we compute spacetime correlation functions conditioned on a fixed, typical measurement record. In the presence of global symmetries we show that these conditional ensembles can undergo measurement-induced “sharpening” phase transitions as a function of the monitoring rate; moreover, even weak monitoring can give rise to novel critical phases, derived entirely from a classical perspective. We give a simple hydrodynamic derivation of the known “charge-fuzzy phase” for weakly monitored diffusive many-body quantum systems. We show that although the unmonitored symmetric and asymmetric exclusion processes are in different universality classes of transport, the fluctuations in their conditional ensembles flow to the same fixed point with emergent relativistic invariance under monitoring. On the other hand, weakly monitored systems with non-Abelian symmetries enter a novel strongly coupled fixed point with nontrivial dynamical exponent, which we characterize. Our formalism naturally accounts for monitoring general observables, such as currents or density gradients, and allows for a direct calculation of information-theoretic diagnostics of sharpening transitions, including the Shannon entropy of the measurement record.

classical statistical mechanics↗
DOE OSTI2022

Integrable nonunitary quantum circuits

Here we show that the integrable Lindblad superoperators found recently can be used to build integrable nonunitary quantum circuits with two-site gates by demonstrating that the R(SIC) matrices are completely positive and trace preserving. Using the bond-site transformation, we obtain the corresponding integrable nonunitary quantum circuits with three-site gates. When restricted to the diagonals of the density matrix, these quantum models reduce to integrable classical cellular automata. This approach may pave the way for a systematic construction of circuits.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Quantum Ornstein-Zernike theory for two-temperature two-component plasmas

Laboratory plasma production almost always preferentially heats either the ions or electrons, leading to a two-temperature state. In this state, density functional theory molecular dynamic simulation is the state of the art for modeling bulk material properties. We construct a statistical mechanics model for the two temperature limit that is theoretically consistent with the molecular dynamics method. We proceed to derive the electron-ion multi-temperature quantum Ornstein-Zernike equations for the first time. This allows the construction of a two-temperature two-component plasma model using the average atom from which we can compute bulk material properties at a fraction of the computation time of the two-temperature density functional theory simulation. The accuracy of the model is benchmarked against ion pair correlation and self-diffusion results from ab initio simulation. Here, we proceed to compute the viscosity and ion thermal conductivity as a function of both ion and electron temperature.

Ab initio molecular dynamics↗
DOE OSTI2022

Model-Free Measurement of Local Entropy Production and Extractable Work in Active Matter

Time-reversal symmetry breaking and entropy production are universal features of nonequilibrium phenomena. Despite its importance in the physics of active and living systems, the entropy production of systems with many degrees of freedom has remained of little practical significance because the high dimensionality of their state space makes it difficult to measure. Here, in this work, we introduce a local measure of entropy production and a numerical protocol to estimate it. We establish a connection between the entropy production and extractability of work in a given region of the system and show how this quantity depends crucially on the degrees of freedom being tracked. We validate our approach in theory, simulation, and experiments by considering systems of active Brownian particles undergoing motility-induced phase separation, as well as active Brownian particles and E.coli in a rectifying device in which the time-reversal asymmetry of the particle dynamics couples to spatial asymmetry to reveal its effects on a macroscopic scale.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Comparison between explicit and implicit discretization strategies for a dissipative thermal environment

We investigate strategies for simulating open quantum systems coupled to dissipative baths by comparing explicit wave function-based discretization [via multi-layer multi-configuration time-dependent Hartree (ML-MCTDH)] and the implicit density matrix-based master equation method [via tree tensor network hierarchical equations of motion (TTN-HEOM)]. For dissipative baths characterized by exponentially decaying bath correlation functions, the implicit discretization approach of HEOM—rooted in bath correlation function decompositions—proves significantly more efficient than explicit discretization of the bath into discrete harmonic modes. Explicit methods, like ML-MCTDH, require extensive mode discretization to approximate continuum baths, leading to computational bottlenecks. Case studies for two-level systems and a Fenna–Matthews–Olson complex model highlight TTN-HEOM’s superiority in capturing dissipative dynamics with relaxations with a minimal number of auxiliary modes, while the explicit methods are as exact as the HEOM in pure dephasing regimes. This comparison is enabled by the TENSO package, which has both ML-MCTDH and TTN-HEOM implemented using the same computational structure and propagation strategy.

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

Thermalization of Interacting Quasi-One-Dimensional Systems

Many experimentally relevant systems are quasi-one-dimensional, consisting of nearly decoupled chains. In these systems, there is a natural separation of scales between the strong intrachain interactions and the weak interchain coupling. When the intrachain interactions are integrable, weak interchain couplings play a crucial part in thermalizing the system. Here, we develop a Boltzmann-equation formalism involving a collision integral that is asymptotically exact for any interacting integrable system, and apply it to develop a quantitative theory of relaxation in coupled Bose gases in the experimentally relevant Newton's cradle setup. Here we find that relaxation involves a broad spectrum of timescales. We provide evidence that the Markov process governing relaxation at late times is gapless; thus, the approach to equilibrium is generally nonexponential, even for spatially uniform perturbations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Universal Kardar-Parisi-Zhang Dynamics in Integrable Quantum Systems

Although the Bethe ansatz solution of the spin-1/2 Heisenberg model dates back nearly a century, the anomalous nature of its high-temperature transport dynamics has only recently been uncovered. Indeed, numerical and experimental observations have demonstrated that spin transport in this paradigmatic model falls into the Kardar-Parisi-Zhang (KPZ) universality class. This has inspired the significantly stronger conjecture that KPZ dynamics, in fact, occur in all integrable spin chains with non-Abelian symmetry. Here, in this letter, we provide extensive numerical evidence affirming this conjecture. Moreover, we observe that KPZ transport is even more generic, arising in both supersymmetric and periodically driven models. Motivated by recent advances in the realization of SU(N)-symmetric spin models in alkaline-earth-based optical lattice experiments, we propose and analyze a protocol to directly investigate the KPZ scaling function in such systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

A practical approach to calculating magnetic Johnson noise for precision measurements

Magnetic Johnson noise is an important consideration for many applications involving precision magnetometry, and its significance will only increase in the future with improvements in measurement sensitivity. The fluctuation–dissipation theorem can be utilized to derive analytic expressions for magnetic Johnson noise in certain situations, but when used in conjunction with finite element analysis tools, the combined approach is particularly powerful as it provides a practical means to calculate the magnetic Johnson noise arising from conductors of arbitrary geometry and permeability. In this paper, we demonstrate this method to be one of the most comprehensive approaches presently available to calculate thermal magnetic noise. In particular, its applicability is shown to not be limited to cases where the noise is evaluated at a point in space but also can be expanded to include cases where the magnetic field detector has a more general shape, such as a finite-size loop, a gradiometer, or a detector that consists of a polarized atomic species trapped in a volume. Furthermore, some physics insights gained through studies made using this method are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Subexponential Decay of Local Correlations from Diffusion-Limited Dephasing

Chaotic quantum systems at finite entropy density are expected to act as their own heat baths, rapidly dephasing local quantum superpositions. Here, we argue that in fact this dephasing is generically subexponential in one-dimensional systems with conservation laws: all local correlation functions decay as exp⁡[−𝒪⁡(𝑡 𝛼 )] with 0 ≤ 𝛼 ≤ 2/3, even when the operators are orthogonal to all hydrodynamic modes. The mechanism is diffusion-limited dephasing, in which rare low-entropy regions (“voids”) protect quantum coherences. This intrinsically quantum effect lies beyond standard hydrodynamics and disappears under extrinsic dephasing. In random charge-conserving circuits we find 𝛼 = 1/2, while in generic translation-invariant Floquet systems we bound 𝛼 ≤ 2/3. Our arguments are general, subject principally to the assumption that thermal fluctuations can create regions of zero entropy density. In systems with energy conservation, this assumption is automatically satisfied because of the third law of thermodynamics.

information scrambling↗
DOE OSTI2023

Simulating dielectric spectra: A demonstration of the direct electric field method and a new model for the nonlinear dielectric response

Here we demonstrate a method to compute the dielectric spectra of fluids in molecular dynamics (MD) by directly applying electric fields to the simulation. We obtain spectra from MD simulations with low magnitude electric fields (≈0.01 V/Å) in agreement with spectra from the fluctuation–dissipation method for water and acetonitrile. We examine this method’s trade-off between noise at low field magnitudes and the nonlinearity of the response at higher field magnitudes. We then apply the Booth equation to describe the nonlinear response of both fluids at low frequency (0.1 GHz) and high field magnitude (up to 0.5 V/Å). We develop a model of the frequency-dependent nonlinear response by combining the Booth description of the static nonlinear dielectric response of fluids with the frequency-dependent linear dielectric response of the Debye model. We find good agreement between our model and the MD simulations of the nonlinear dielectric response for both acetonitrile and water.

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

Theory of the ion–electron temperature relaxation rate in strongly magnetized plasmas

Recent works have shown that strongly magnetized plasmas characterized by having a gyrofrequency greater than the plasma frequency exhibit novel transport properties. One example is that the friction force on a test charge shifts, obtaining components perpendicular to its velocity in addition to the typical stopping power component antiparallel to its velocity. Here, we apply a recent generalization of the Boltzmann equation for strongly magnetized plasmas to calculate the ion–electron temperature relaxation rate. Strong magnetization is generally found to increase the temperature relaxation rate perpendicular to the magnetic field and to cause the temperatures parallel and perpendicular to the magnetic field to not relax at equal rates. This, in turn, causes a temperature anisotropy to develop during the equilibration. Strong magnetization also breaks the symmetry of independence of the sign of the charges of the interacting particles on the collision rate, commonly known as the “Barkas effect.” Furthermore, it is found that the combination of oppositely charged interaction and strong magnetization causes the ion–electron parallel temperature relaxation rate to be significantly suppressed, scaling inversely proportional to the magnetic field strength.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗