SEARCH · Engineering Papers
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.
In this work, we explore time-resolved dynamics of non-equilibrium magnons in Rb2MnF4, a model two-dimensional square lattice Heisenberg antiferromagnet using a novel laser-neutron pump-probe technique. This technique uses a nanosecond pulsed laser to periodically excite magnons to out-of-equilibrium states and then probe the resulting dynamics using INS.
Nonequilibrium dynamics play an important role in all contexts of physics, both classical and quantum as well as living and nonliving, so it is crucial to develop a foundational understanding of nonequilibrium phase transitions. In this work we investigate an important class of nonequilibrium dynamics in the form of nonreciprocal interactions. In particular we study how nonreciprocal coupling between two 𝑂(𝑛𝑖) order parameters (with 𝑖 = 1,2) affects the universality at a multicritical point, extending the analysis of J. T. Young et al. [Phys. Rev. X 10, 011039 (2020)], which considered the case 𝑛 1 = 𝑛 2 = 1, i.e., a ℤ 2 × ℤ 2 model. We show that nonequilibrium fixed points (NEFPs) emerge for a broad range of 𝑛 1 ,𝑛 2 and exhibit intrinsically nonequilibrium critical phenomena, namely a violation of fluctuation-dissipation relations at all scales and underdamped oscillations near criticality in contrast to the overdamped relaxational dynamics of the corresponding equilibrium models. Furthermore, the NEFPs exhibit an emergent discrete scale invariance in certain physically relevant regimes of 𝑛 1 ,𝑛 2 , but not others, depending on whether the critical exponent 𝜈 is real or complex. The boundary between these two regions is described by an exceptional point in the renormalization group (RG) flow, leading to distinctive features in correlation functions and the phase diagram. Another contrast with the previous work is the number and stability of the NEFPs as well as the underlying topology of the RG flow. Lastly, we investigate an extreme form of nonreciprocity where one order parameter is independent of the other order parameter but not vice versa. Unlike the ℤ 2 × ℤ 2 model, which becomes nonperturbative in this case, we identify a distinct nonequilibrium universality class whose dependent field similarly violates fluctuation-dissipation relations but does not exhibit discrete scale invariance or underdamped oscillations near criticality.
Nonequilibrium states in soft condensed matter require a systematic approach to characterize and model materials, enhancing predictability and applications. Among the tools, X-ray photon correlation spectroscopy (XPCS) provides exceptional temporal and spatial resolution to extract dynamic insight into the properties of the material. However, existing models might overlook intricate details. We introduce an approach for extracting the transport coefficient, denoted as $J(t)$, from the XPCS studies. This coefficient is a fundamental parameter in nonequilibrium statistical mechanics and is crucial for characterizing transport processes within a system. Our method unifies the Green–Kubo formulas associated with various transport coefficients, including gradient flows, particle–particle interactions, friction matrices, and continuous noise. We achieve this by integrating the collective influence of random and systematic forces acting on the particles within the framework of a Markov chain. We initially validated this method using molecular dynamics simulations of a system subjected to changes in temperatures over time. Subsequently, we conducted further verification using experimental systems reported in the literature and known for their complex nonequilibrium characteristics. The results, including the derived $J(t)$ and other relevant physical parameters, align with the previous observations and reveal detailed dynamical information in nonequilibrium states. This approach represents an advancement in XPCS analysis, addressing the growing demand to extract intricate nonequilibrium dynamics. Further, the methods presented are agnostic to the nature of the material system and can be potentially expanded to hard condensed matter systems.
Coherent anti-Stokes Raman scattering (CARS) and nitric oxide molecular tagging velocimetry (NO-MTV) were used to characterize the freestream in Sandia’s Hypersonic Shock Tunnel (HST) using a burst-mode laser operated at 100 kHz. Experiments were performed at nominal freestream velocities of 3 and 4 km/s using both air and N 2 test gas. The CARS diagnostic provides nonequilibrium characterization of the flow by measuring vibrational and rotational temperatures of N 2 and O 2 , which are compared to NO temperatures from separate laser absorption experiments. Simultaneous, collinear freestream velocities were measured using NO MTV along with pitot pressure measurements. Furthermore, this extensive freestream dataset is compared to nonequilibrium CFD that is capable of modeling species-specific, vibrational temperatures throughout the nozzle expansion. Significant nonequilibrium between vibrational and rotational temperatures was measured at each flow condition. N 2 exhibits the most nonequilibrium, followed by O 2 and NO. The CFD model captures this trend, although it consistently overpredicts N 2 and O 2 vibrational temperatures. At 3 km/s, the modeled NO nonequilibrium is underpredicted, whereas it is overpredicted at 4 km/s. Good agreement is seen between CFD and the velocity and rotational temperature measurements. Experiments with water added to the test gas yielded no discernable difference in vibrational relaxation.
Understanding the dynamics of nonequilibrium quantum many-body systems is an important research topic in a wide range of fields across condensed matter physics, quantum optics, and high-energy physics. However, numerical studies of large-scale nonequilibrium phenomena in realistic materials face serious challenges due to intrinsic high-dimensionality of quantum many-body problems and the absence of time-invariance. The nonequilibrium properties of many-body systems can be described by the dynamics of the correlator, or the Green's function of the system, whose time evolution is given by a high-dimensional system of integro-differential equations, known as the Kadanoff–Baym equations (KBEs). The time-convolution term in KBEs, which needs to be recalculated at each time step, makes it difficult to perform long-time numerical simulation. In this paper, we develop an operator-learning framework based on recurrent neural networks (RNNs) to address this challenge. We utilize RNNs to learn the nonlinear mapping between Green's functions and convolution integrals in KBEs. By using the learned operators as a surrogate model in the KBE solver, we obtain a general machine-learning scheme for predicting the dynamics of nonequilibrium Green's functions. Besides significant savings per each time step, the new methodology reduces the temporal computational complexity from $O(N_t^3)$ to $O(N_t)$ where N t is the number of steps taken in a simulation, thereby making it possible to study large many-body problems which are currently infeasible with conventional KBE solvers. Through various numerical examples, we demonstrate the effectiveness of the operator-learning based approach in providing accurate predictions of physical observables such as the reduced density matrix and time-resolved photoemission spectra. Moreover, our framework exhibits clear numerical convergence and can be easily parallelized, thereby facilitating many possible further developments and applications.
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.
Coal is an abundant natural resource and there is motivation to find new uses for it that do not intrinsically involve combustion. One approach is to explore new ways of processing coal, and in this work, we focus on the transformation of coal in a nonequilibrium plasma generated from an equimolar mixture of nitrogen and hydrogen. The outcome of the nonequilibrium plasma reaction is fundamentally different than a thermal control reaction carried out using the same gas composition, pressure, and temperature range. The nonequilibrium plasma produces a gas mixture that is enriched in acetylene and its derivatives. Furthermore, when compared to the thermal control experiment, the solid char byproduct of the nonequilibrium plasma has a very reactive surface and is spontaneously combustible at ambient temperature. Experiments performed to characterize the reaction kinetics of coal in the plasma suggest that the mechanism proceeds through a sequential process by which the coal particle temperature rises to a point where devolatilization can occur, the devolatilization reaction happens, followed by parallel reactions of released organic vapors in the plasma phase and surface activation. In conclusion, the reaction rate appears to be limited by the time it takes for the coal particle temperature to rise, consistent with previous results reported for reactions of coal in thermal plasma.
Ammonia (NH 3 ) is one of the most important industrial chemicals. The conventional NH 3 synthesis method-the Haber–Bosch process-converts atmospheric nitrogen (N 2 ) into NH 3 using H 2 with an iron (Fe) catalyst. However, this process requires high pressures (100–200 atm) and temperatures (700–800 K) near thermal equilibrium. Recently, Fe-based nanocatalysts have been reported to produce promising NH 3 yields under atmospheric pressures and temperature-modulated nonequilibrium conditions. Understanding the mechanism of nonequilibrium catalysis with programmed temperature variation could help to optimize this fully electrified and less energy-intensive process. Although reactive molecular dynamics (RMD) simulations can be a useful tool to model nonequilibrium catalytic processes, they require the development of accurate force fields (i.e., interatomic potentials). Here, we present a machine-learned (ML) force field within the Deep Potential MD (DPMD) framework, trained using periodic density functional theory (DFT) calculations, to model NH 3 synthesis on Fe catalysts with various surface adsorbates such as *N, *H, *N 2 , *H 2 , *NH, *NH 2 , and *NH 3 . Here, we generated the DFT data from static models of elementary reactions on the most stable (110) surface of body-centered cubic Fe, which then were augmented by data from constant number of particles–volume–temperature (NVT) DFT-MD trajectories at various temperatures. Finally, we utilized the fully optimized ML force field to investigate reaction dynamics at an Fe(110) surface at linearly increasing temperatures using NVT-DPMD simulations. Our simulations indicate that pulsed temperature ramping could prove favorable for NH3 synthesis. For example, we conducted ramping under multiple sets of conditions: (i) from 900 to 1200 K over periods of 0.1–0.3 ns for Fe surfaces precovered with N or NH along with H; and (ii) from 300 to 600 K over 0.1–0.3 ns for Fe surfaces precovered with NH 3 . While our simulations so far are limited to short time scales (very rapid heating), these observations shed light on the mechanism of the high NH 3 synthesis rate achieved in a novel temperature-modulated nonequilibrium catalytic reactor using pulsed heating and cooling.
The coupling between quasiparticles and bosonic excitations rules the energy transfer pathways in condensed matter systems. The possibility of inferring the strength of specific coupling channels from their characteristic timescales measured in nonequilibrium experiments is still an open question. Here, in this study, we investigate MgB 2 , in which conventional superconductivity at temperatures as high as 39 K is mediated by the strong coupling between the conduction electrons and the 𝐸 2𝑔 phonon mode. By means of broadband time-resolved optical spectroscopy, we show that this selective electron-phonon coupling dictates the nonequilibrium optical response of MgB 2 at early times (<100 fs) after photoexcitation. Furthermore, based on an effective temperature model analysis, we estimate its contribution to the total electron-boson coupling function extracted from complementary equilibrium spectroscopy approaches, namely, optical reflectivity and angle-resolved photoemission spectroscopy. The coupling strength with the 𝐸 2𝑔 phonon modes is thus estimated to be 𝜆 ≃ 0.56, which is approximately half of the total coupling constant, in agreement with 𝑎𝑏 initio calculations from the literature. As a benchmark, broadband time-resolved optical spectroscopy is performed also on the isostructural and nonsuperconducting compound AlB 2 , showing that the nonequilibrium optical response relaxes on a slower timescale due to the lack of strongly coupled phonon modes. Our findings demonstrate the possibility to resolve and quantify selective electron-phonon coupling from nonequilibrium optical spectroscopy.
Using matrix product state techniques we study the nonequilibrium dynamical response of the half-filled Hubbard ladder when subject to an optical pump. Optical pumping offers a way of producing and manipulating new strongly correlated phenomena by suppressing existing magnetic correlations. The ladder allows the effects of pump directionality to be investigated, and compared to a single chain it has strong spin-charge coupling and a fully gapped excitation spectrum, promising different nonequilibrium physics. We compute time-dependent correlations, including the nonequilibrium dynamical structure factors for spin and charge. By deriving a combined spin-charge sum rule that applies both in and out-of-equilibrium, we show that spectral weight is pumped directly from the antiferromagnetic spin response into a low energy 𝜔 ∼ 0 charge response below the Mott gap. The transfer of weight is pump direction dependent: pumping directed along the legs disrupts magnetic correlations more than pumping in the rung-direction, even if the post-pump energy density is similar. The charge correlation length is dramatically enhanced by the pump, whilst the spin correlations are most strongly suppressed at nearest- and next-nearest-neighbor spacings. After the pump the system is in a nonthermal correlated metallic state, with gapless charge excitations and approximately equal spin and charge correlation lengths, emphasizing the importance of treating these degrees of freedom on an equal footing in nonequilibrium systems.
During his time at Royal Institute of Technology (Kungliga Tekniska högskolan) in Sweden, the present author learned nonequilibrium thermodynamics from Mats Hillert. The key concepts are the separation of internal and external variables of a system and the definitions of potentials and molar quantities. In equilibrium thermodynamics derived by Gibbs, the internal variables are not independent and can be fully evaluated from given external variables. While irreversible thermodynamics led by Onsager focuses on internal variables though often mixed with external variables. Hillert integrated them together by first emphasizing their differences and then examining their connections. His philosophy was reflected by the title of his book “Phase Equilibria, Phase Diagrams and Phase Transformations” that puts equilibrium, nonequilibrium, and internal processes on equal footing. Here, in the present paper honoring Hillert, the present author reflects his experiences with Hillert and his work in last 40 years and expresses his gratitude for all the wisdom and support from him in terms of “Hillert nonequilibrium thermodynamics” and discusses some recent topics that the present author has been working on.
A key objective in nuclear and high-energy physics is to describe nonequilibrium dynamics of matter, e.g., in the early Universe and in particle colliders, starting from the standard model of particle physics. Classical computing methods, via the framework of lattice gauge theory, have experienced limited success in this mission. Quantum simulation of lattice gauge theories holds promise for overcoming computational limitations. Because of local constraints (Gauss’s laws), lattice gauge theories have an intricate Hilbert-space structure. This structure complicates the definition of thermodynamic properties of systems coupled to reservoirs during equilibrium and nonequilibrium processes. We show how to define thermodynamic quantities such as work and heat using strong-coupling thermodynamics, a framework that has recently burgeoned within the field of quantum thermodynamics. Our definitions suit instantaneous quenches, simple nonequilibrium processes undertaken in quantum simulators. To illustrate our framework, we compute the work and heat exchanged during a quench in a Z 2 lattice gauge theory coupled to matter in 1+1 dimensions. Here, the thermodynamic quantities, as functions of the quench parameter, evidence a phase transition. For general thermal states, we derive a simple relation between a quantum many-body system’s entanglement Hamiltonian, measurable with quantum-information-processing tools, and the Hamiltonian of mean force, used to define strong-coupling thermodynamic quantities.
Two-level systems (TLSs) are a dominant decoherence mechanism in superconducting qubits, microwave kinetic inductance detectors, and quantum dots. TLSs are tunneling states commonly found in amorphous materials and interfaces, which interact electromechanically, leading to energy loss. Fluctuations in the frequency of the TLSs are a significant problem for qubit operation and force recalibration every few hours. In this study, we probe the effect of alternating bias at cryogenic temperatures on TLSs in a purpose-built LC oscillator. When an in situ alternating bias is applied, the steady-state spectra of the TLSs disappear. Spectroscopy reveals transient behavior, in which the TLS frequency fluctuates on the order of minutes. Thermal cycling above 10 K reverses these effects, restoring the TLS spectrum to its original state. Importantly, the intrinsic loss tangent of the LC oscillator remains unchanged before and after the application of the alternating bias. We propose that the disappearance of the steady-state spectrum is caused by nonequilibrium energy buildup from strain in the oxide film introduced by the pulsed voltage-bias sequence. Understanding this nonequilibrium energy could shed light on TLS fluctuations and reduce the calibration requirements for qubits.
The evaporation model for water isotopes proposed by Craig and Gordon (1965, https://books.google.co.in/books?id=6wIKAQAAIAAJ) is used in most isotope-enabled atmospheric models for the parameterization of nonequilibrium fractionation during evaporation from the ocean. In this model, one of the most uncertain parameters is the nonequilibrium fractionation factor 𝛼 𝑘 . Many isotope models use the formulation of Merlivat and Jouzel (1979, https://doi.org/10.1029/jc084ic08p05029), which parameterizes 𝛼 𝑘 as a function of wind speed and distinguishes between a smooth and a rough regime to account for the effect of ocean waves. The resulting discontinuity in 𝛼 𝑘 between smooth and rough regimes has been disputed by several empirical studies. Here, we present a new approach to parameterizing 𝛼𝑘 by explicitly accounting for the influence of wave drag on the momentum flux near the surface. Following Cifuentes-Lorenzen et al. (2018, https://doi.org/10.1007/s10546-018-0376-0), we add a third wave-induced component to the total momentum flux, in addition to the viscous and turbulent components, and extend the definition of the eddy viscosity to account for the momentum flux due to waves and turbulent dissipation near the surface. The new scheme predicts a slight decrease of 𝛼 𝑘 with wind speed, similar to the smooth-regime parameterization of Merlivat and Jouzel (1979, https://doi.org/10.1029/jc084ic08p05029). This new parameterization is incorporated into the isotope-enabled Community Atmosphere Model, where it improves the correlation of simulated and measured vapor deuterium excess relative to the default version and a version with constant 𝛼 𝑘 , suggesting that it may be used as a valid representation of fractionation during evaporation from the ocean in future isotope models.
Electronic structure calculations in the time domain provide a deeper understanding of nonequilibrium dynamics in materials. The real-time Boltzmann equation (rt-BTE), used in conjunction with accurate interactions computed from first principles, has enabled reliable predictions of coupled electron and lattice dynamics. However, the timescales and system sizes accessible with this approach are still limited, with two main challenges being the different timescales of electron and phonon interactions and the cost of computing collision integrals. As a result, only a few examples of these calculations exist, mainly for two-dimensional (2D) materials. Here we leverage adaptive and multirate time integration methods to achieve a major step forward in solving the coupled rt-BTEs for electrons and phonons. Relative to conventional (non-adaptive) time-stepping, our approach achieves a 10x speedup for a target accuracy, or greater accuracy by 3–6 orders of magnitude for the same computational cost, enabling efficient calculations in both 2D and bulk materials. This efficiency is showcased by computing the coupled electron and lattice dynamics in graphene up to ~100 ps, as well as modeling ultrafast lattice dynamics and thermal diffuse scattering maps in bulk materials (silicon and gallium arsenide). In addition to improved efficiency, our adaptive method can resolve the characteristic rates of different physical processes, thus naturally bridging different timescales. This enables simulations of longer timescales and provides a framework for modeling multiscale dynamics of coupled degrees of freedom in matter. Our work opens new opportunities for quantitative studies of nonequilibrium physics in materials, including driven lattice dynamics with phonons coupled to electrons, spin, and other degrees of freedom.
Aims.The goal is to develop a database of rate coefficients for rotational state-to-state transitions in H 2 O + H 2 O collisions that is suitable for the modeling of energy transfer in nonequilibrium conditions, in which the distribution of rotational states of H 2 O deviates from local thermodynamic equilibrium. Methods.A two-temperature model was employed that assumed that although there is no equilibrium between all possible degrees of freedom in the system, the translational and rotational degrees of freedom can be expected to achieve their own equilibria independently, and that they can be approximately characterized by Boltzmann distributions at two different temperatures,T kin andT rot . Results.Upon introducing our new parameterization of the collisional rates, taking into account their dependence on bothT kin andT rot , we find a change of up to 20% in the H 2 O rotational level populations for both ortho and para-H 2 O for the part of the cometary coma where the nonequilibrium regime occurs.
Landauer's principle states that erasing a bit of information at fixed temperature T costs at least $k$ B $T$ ln 2 units of work. Here we investigate erasure at varying temperature, to which Landauer's result does not apply. Here we formulate bit erasure as a stochastic nonequilibrium process involving a compression of configuration space, with physical and logical states associated in a symmetric way. Erasure starts and ends at temperature T, but temperature can otherwise vary with time in an arbitrary way. Defined in this way, erasure is governed by a set of nonequilibrium fluctuation relations that show that varying-temperature erasure can done with less work than $k$ B $T$ ln 2. As a result, erasure and the complementary process of bit randomization can be combined to form a work-producing engine cycle.
Using equilibrium fluctuations to understand the response of a physical system to an externally imposed perturbation is the basis for linear response theory, which is widely used to interpret experiments and shed light on microscopic dynamics. For nonequilibrium systems, perturbations cannot be interpreted simply by monitoring fluctuations in a conjugate observable and general response results rely on path ensemble averaging. Furthermore, these techniques do not apply to perturbations that affect the diffusion tensor in a stochastic system. Here, we introduce an “effective” physical process that represents the diffusion perturbed dynamics and enables accurate calculations of responses to a change in the diffusion. Interestingly, the effective dynamics contain an additional drift involving the instantaneous “score” of the system, and we leverage score matching algorithms to carry out nonequilibrium response calculations on systems for which the exact stationary distribution is unknown.