Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonequilibrium dynamics”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Modeling Equilibrium Solid–Liquid Interfaces under Effective Constant Chemical Potential Using Machine Learning Interatomic Potentials

The chemical potential (μ) of species in solution is essential for understanding various chemical processes at interfaces. Molecular dynamics (MD) simulations, constrained by fixed compositions, cannot maintain constant chemical potential with reference to a targeted concentration or chemical potential under nonequilibrium or dynamic conditions, as solute species can migrate to the interface and deplete (or enrich) the bulk due to solute-interface interactions. In this study, we introduce a simple and computationally efficient approach named iterative quasi-constant chemical potential molecular dynamics (iqCμMD) simulation, which helps simulate targeted molar concentrations of species in solution. iqCμMD overcomes the limitations of conventional MD by adjusting the number of species in the solution to reach a target bulk concentration (chemical potential), which allows simulation of the interface under the bulk conditions comparable to experiment. We demonstrate our approach using machine learning interatomic potential (MLIP)-based MD simulations of the Na 2 SO 4,aq –graphene interface, and to show the transferability of our approach, we also perform classical force field-based MD simulations of NaCl aq –air and NaCl aq –graphite interfaces, which produce comparable results to previous CμMD simulations. Our results also show that the iqCμMD approach efficiently achieves the desired bulk ion concentration within two iterations, and by utilizing MLIPs, we can achieve converged results using relatively small-scale simulations compared to previous CμMD simulations. By combining iqCμMD with MLIP-driven simulations, solid–liquid interfaces can be modeled under an effective constant chemical potential with DFT-level accuracy. Here, we show that iqCμMD offers a robust and simple computational framework for constant chemical potential simulations, as its only requirement is to be able to converge interfacial simulations with a measurable bulk region.

Chemical structure↗
DOE OSTI2023

Real-time evolution of Anderson impurity models via tensor network influence functionals

In this work, we present and analyze two tensor network-based influence functional approaches for simulating the real-time dynamics of quantum impurity models such as the Anderson model. Via comparison with recent numerically exact simulations, we show that such methods accurately capture the long-time nonequilibrium quench dynamics. The two parameters that must be controlled in these tensor network influence functional approaches are a time discretization (Trotter) error and a bond dimension (tensor network truncation) error. We show that the actual numerical uncertainties are controlled by an intricate interplay of these two approximations, which we demonstrate in different regimes. Our work opens the door to using these tensor network influence functional methods as general impurity solvers.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Random Quantum Circuits

Quantum circuits—built from local unitary gates and local measurements—are a new playground for quantum many-body physics and a tractable setting to explore universal collective phenomena far from equilibrium. These models have shed light on longstanding questions about thermalization and chaos, and on the underlying universal dynamics of quantum information and entanglement. In addition, such models generate new sets of questions and give rise to phenomena with no traditional analog, such as dynamical phase transitions in quantum systems that are monitored by an external observer. Quantum circuit dynamics is also topical in view of experimental progress in building digital quantum simulators that allow control of precisely these ingredients. Randomness in the circuit elements allows a high level of theoretical control, with a key theme being mappings between real-time quantum dynamics and effective classical lattice models or dynamical processes. Many of the universal phenomena that can be identified in this tractable setting apply to much wider classes of more structured many-body dynamics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Nonthermal Bonding Origin of a Novel Photoexcited Lattice Instability in SnSe

Lattice dynamics measurements are often crucial tools for understanding how materials transform between different structures. Here, we report time-resolved x-ray scattering-based measurements of the nonequilibrium lattice dynamics in SnSe, a monochalcogenide reported to host a novel photoinduced lattice instability. By fitting interatomic force models to the fluence dependent excited-state dispersion, we determine the nonthermal origin of the lattice instability to be dominated by changes of interatomic interactions along a bilayer-connecting bond, rather than of an intralayer bonding network that is of primary importance to the lattice instability in thermal equilibrium.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Ultrafast visualization of phase transitions in nonequilibrium warm dense matter

Intense ultrafast laser excitation brings materials into highly nonequilibrium states with complex solid–liquid phase transitions. These ultrafast processes yield warm dense matter (WDM) conditions characterized by comparable thermal and Fermi energies, and strong ion–ion coupling. Here, we review recent studies employing mega-electron-volt ultrafast electron diffraction to resolve the structural dynamics of nonequilibrium WDM created by femtosecond laser heating of solids. Furthermore, for the study of warm dense gold, the results showed homogeneous melting that occurs within tens of picoseconds at absorbed energy densities of 0.4–1.4 MJ/kg. These results constrained the electron–ion coupling rate and revealed the melting sensitivity to nucleation seeds. For the study of radiation-damaged tungsten, the results showed a melting transition below the melting temperature. Molecular dynamics simulations suggested that this melting behavior is driven by vacancy defect clusters from radiation damage. These studies provided atomic-level insights into the melting behavior of materials under extreme conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Tuning between continuous time crystals and many-body scars in long-range XYZ spin chains.

Persistent oscillatory dynamics in nonequilibrium many-body systems is a tantalizing manifestation of ergodicity breakdown that continues to attract much attention. Recent works have focused on two classes of such systems: discrete time crystals and quantum many-body scars (QMBS). While both systems host oscillatory dynamics, its origin is expected to be fundamentally different: discrete time crystal is a phase of matter which spontaneously breaks the Z 2 symmetry of the external periodic drive, while QMBS span a subspace of nonthermalizing eigenstates forming an su(2) algebra representation. Here, we ask a basic question: is there a physical system that allows us to tune between these two dynamical phenomena? In contrast to much previous work, we investigate the possibility of a continuous time crystal (CTC) in undriven, energy-conserving systems exhibiting prethermalization. We introduce a long-range XYZ spin model and show that it encompasses both a CTC phase as well as QMBS. We map out the dynamical phase diagram using numerical simulations based on exact diagonalization and time-dependent variational principle in the thermodynamic limit. We identify a regime where QMBS and CTC order coexist, and we discuss experimental protocols that reveal their similarities as well as key differences.

Bull, Kieran↗
DOE OSTI2024

Hidden phonon highways promote photoinduced interlayer energy transfer in twisted transition metal dichalcogenide heterostructures

Vertically stacked van der Waals (vdW) heterostructures exhibit unique electronic, optical, and thermal properties that can be manipulated by twist-angle engineering. However, the weak phononic coupling at a bilayer interface imposes a fundamental thermal bottleneck for future two-dimensional devices. Using ultrafast electron diffraction, we directly investigated photoinduced nonequilibrium phonon dynamics in MoS 2 /WS 2 at 4° twist angle and WSe 2 /MoSe 2 heterobilayers with twist angles of 7°, 16°, and 25°. We identified an interlayer heat transfer channel with a characteristic timescale of ~20 picoseconds, about one order of magnitude faster than molecular dynamics simulations assuming initial intralayer thermalization. Atomistic calculations involving phonon-phonon scattering suggest that this process originates from the nonthermal phonon population following the initial interlayer charge transfer and scattering. Our findings present an avenue for thermal management in vdW heterostructures by tailoring nonequilibrium phonon populations.

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

Machine-learning modeling of magnetization dynamics in quasi-equilibrium and driven metallic spin systems

Here, we present a perspective on recent progress in machine-learning (ML) force-field approaches for large-scale Landau–Lifshitz–Gilbert (LLG) simulations of metallic spin systems. Building on a generalization of the Behler–Parrinello (BP) architecture originally developed for quantum molecular dynamics, we develop scalable and transferable ML models that faithfully capture the complex, environment-dependent electron-mediated exchange fields characteristic of itinerant magnets. A central ingredient of this framework is the implementation of symmetry-aware magnetic descriptors based on group-theoretical bispectrum formalisms. Leveraging these ML force fields, LLG simulations faithfully reproduce hallmark non-collinear magnetic orders—such as the 120° and tetrahedral states—on the triangular lattice, and successfully capture the complex spin textures emerging in the mixed-phase states of a square-lattice double-exchange model under thermal quench. We further discuss a generalized potential theory that extends the BP formalism to incorporate both conservative and nonconservative electronic torques, thereby enabling ML models to learn nonequilibrium exchange fields from computationally demanding microscopic approaches such as nonequilibrium Green’s-function techniques. This extension yields quantitatively accurate predictions of voltage-driven domain-wall motion and establishes a foundation for quantum-accurate, multiscale modeling of nonequilibrium spin dynamics and spintronic functionalities.

Descriptors↗
DOE OSTI2026

Slow Quasiparticle Dynamics and Anyonic Statistics in a Fractional Quantum Hall Fabry-Pérot Interferometer

Anyons are two-dimensional particles with fractional exchange statistics that emerge as elementary excitations of fractional quantum Hall phases. Experimentally, their exchange statistics can be measured in the edge-state Fabry-Pérot interferometer, wherein the presence of 𝑁 𝑞⁢𝑝 localized anyons contributes a phase 𝑁 𝑞⁢𝑝⁢ 𝜃 𝑎 to the interference pattern where 𝜃 𝑎 is twice the exchange phase. Here we report the observation of large, hysteretic phase jumps in a monolayer graphene Fabry-Pérot interferometer at 𝜈 = 1/3. When the filling factor is increased from 𝜈 < 1/3 toward the center of the plateau, we observe phase slips with magnitude Δ⁢𝜃 ≈ 2⁢𝜋/3, consistent with the addition of individual quasiparticles to the interferometer bulk. These phase slips occur as instantaneous jumps in the interference signal, with intervals between the jumps indicating quasiparticle equilibration times exceeding 20 min. We use this long timescale to investigate the effect of changes in interferometer area 𝐴 𝐼 and 𝑁 𝑞⁢𝑝 independently at fixed magnetic field, revealing a striking memory effect in the phase slip magnitude. In particular, as the 𝜈 =1/3 plateau is approached from higher filling, we observed phase slips with Δ⁢𝜃 significantly larger than 2⁢𝜋/3 over the same range of gate voltage where quantized jumps are seen for increasing 𝜈. We discuss this asymmetry in terms of bulk-edge coupling of quasiparticles localized near the edge or in the bulk, and argue that this effect can be qualitatively reconciled with theoretical expectations for strongly interacting quasiparticles in the presence of weak disorder and strongly nonequilibrium charge dynamics. Besides providing a replication of interferometric measurements sensitive to 𝜃 𝑎 , our results highlight the key role played by charge dynamics on signatures of the anyon phase, and demonstrate that fractional quasiparticles can be indefinitely localized in nonequilibrium configurations.

anyons↗
DOE OSTI2025

Hidden domain boundary dynamics toward crystalline perfection

A central paradigm of nonequilibrium physics concerns the dynamics of heterogeneity and disorder, impacting processes ranging from the behavior of glasses to the emergent functionality of active matter. Understanding these complex mesoscopic systems requires probing the microscopic trajectories associated with irreversible processes, the role of fluctuations and entropy growth, and the timescales on which nonequilibrium responses are ultimately maintained. Approaches that illuminate these processes in model systems may enable a more general understanding of other heterogeneous nonequilibrium phenomena, and potentially define ultimate speed and energy cost limits for information processing technologies. Here, we apply ultrafast single-shot X-ray photon correlation spectroscopy to resolve the nonequilibrium, heterogeneous, and irreversible mesoscale dynamics during a light-induced phase transition in a (PbTiO 3 ) 16 /(SrTiO 3 ) 16 superlattice. Such ferroelectric superlattice systems are a useful platform to study phase transitions and topological dynamics due to their high degree of tunability. This provides an approach for capturing the nucleation of the light-induced phase, the formation of transient mesoscale defects at the boundaries of the nuclei, and the eventual annihilation of these defects, even in systems with complex polarization topologies. We identify a nonequilibrium correlation response spanning >10 orders of magnitude in timescales, with multistep behavior similar to the plateaus observed in supercooled liquids and glasses. We further show how the observed time-dependent long-time correlations can be understood in terms of stochastic and non-Markovian dynamics of domain walls, encoded in waiting-time distributions with power-law tails. This work defines possibilities for probing the nonequilibrium and correlated dynamics of disordered and heterogeneous media.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Direct observation of strong momentum-dependent electron-phonon coupling in a metal

Phonon scattering in metals is one of the most fundamental processes in materials science. However, understanding such processes has remained challenging and requires detailed information on interactions between phonons and electrons. We use an ultrafast electron diffuse scattering technique to resolve the nonequilibrium phonon dynamics in femtosecond–laser-excited tungsten in both time and momentum. We determine transient populations of phonon modes which show strong momentum dependence initiated by electron-phonon coupling. For phonons near Brillouin zone border, we observe a transient rise in their population on a timescale of approximately 1 picosecond driven by the strong electron-phonon coupling, followed by a slow decay on a timescale of approximately 8 picosecond governed by the weaker phonon-phonon relaxation process. We find that the exceptional harmonicity of tungsten is needed for isolating the two processes, resulting in long-lived nonequilibrium phonons in a pure metal. Our finding highlights that electron-phonon scattering can be the determinant factor in the phonon thermal transport of metals.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Reaction-induced departures from continuum Navier–Stokes turbulence

Reactive hydrodynamic turbulence is an inherently multiscale phenomenon, characterized by the separation between energy-containing, viscous, and molecular length and time scales. The separation between the viscous scale (the Kolmogorov scale) and the molecular mean free path ostensibly justifies a macroscopic description of reactive turbulence via the Navier–Stokes (NS) equations. However, here we use molecular-level simulations to demonstrate that exothermic bimolecular reactions can cause the NS description of turbulence to break down in the near-continuum regime. Sufficiently energetic heat-releasing reactive collisions strongly distort the Maxwell–Boltzmann velocity distribution function, modifying not only the macroscopic chemical rate law but the kinetic-energy-transfer processes as well. This translational nonequilibrium ultimately introduces significant departures from the NS description in the kinetic energy spectra at scales orders of magnitude larger than both the molecular mean free path and the Kolmogorov length scale. These departures prove substantial enough to meaningfully alter integrated quantities, including the overall turbulence kinetic energy itself.

molecular gas dynamics↗
DOE OSTI2023

Integrated Quantum-Classical Protocol for the Realistic Description of Solvated Multinuclear Mixed-Valence Transition-Metal Complexes and Their Solvatochromic Properties

Linear cyanide-bridged polymetallic complexes, which undergo photoinduced metal-to-metal charge transfer, represent prototypical systems for studying long-range electron-transfer reactions and understanding the role played by specific solute–solvent interactions in modulating the excited-state dynamics. Here, to tackle this problem, while achieving a statistically meaningful description of the solvent and of its relaxation, one needs a computational approach capable of handling large polynuclear transition-metal complexes, both in their ground and excited states, as well as the ability to follow their dynamics in several environments up to nanosecond time scales. Here, we present a mixed quantum classical approach, which combines large-scale molecular dynamics (MD) simulations based on an accurate quantum mechanically derived force field (QMD-FF) and self-consistent QMD polarized point charges, with IR and UV–vis spectral calculations to model the solvation dynamics and optical properties of a cyano-bridged trinuclear mixed-valence compound (trans-[(NC) 5 Fe III (μ-CN)Ru II (pyridine) 4 (μ-NC)Fe III (CN) 5 ] 4– ). We demonstrate the reliability of the QMD-FF/MD approach in sampling the solute conformational space and capturing the local solute–solvent interactions by comparing the results with higher-level quantum mechanics/molecular mechanics (QM/MM) MD reference data. The IR spectra calculated along the classical MD trajectories in different solvents correctly predict the red shift of the CN stretching band in the aprotic medium (acetonitrile) and the subtle differences measured in water and methanol, respectively. By explicitly including the solvent molecules around the cyanide ligands and calculating the thermal averaged absorption spectra using time-dependent density functional theory calculations within the Tamm–Dancoff approximation, the experimental solvatochromic shift is quantitatively reproduced going from water to methanol, while it is overestimated for acetonitrile. This discrepancy can likely be traced back to the lack of important dispersion interactions between the solvent cyano groups and the pyridine substituents in our micro solvation model. The proposed protocol is applied to the ground state in water, methanol, and acetonitrile and can be flexibly generalized to study excited-state nonequilibrium solvation dynamics.

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

Solving quasi-free and quadratic Lindblad master equations for open fermionic and bosonic systems

The dynamics of Markovian open quantum systems are described by Lindblad master equations. For fermionic and bosonic systems that are quasi-free, i.e. with Hamiltonians that are quadratic in the ladder operators and Lindblad operators that are linear in the ladder operators, we derive the equation of motion for the covariance matrix. This determines the evolution of Gaussian initial states and the steady states, which are also Gaussian. Using ladder super-operators (a.k.a. third quantization), here we show how the Liouvillian can be transformed to a many-body Jordan normal form which also reveals the full many-body spectrum. Extending previous work by Prosen and Seligman, we treat fermionic and bosonic systems on equal footing with Majorana operators, shorten and complete some derivations, also address the odd-parity sector for fermions, give a criterion for the existence of bosonic steady states, cover non-diagonalizable Liouvillians also for bosons, and include quadratic systems. In extension of the quasi-free open systems, quadratic open systems comprise additional Hermitian Lindblad operators that are quadratic in the ladder operators. While Gaussian states may then evolve into non-Gaussian states, the Liouvillian can still be transformed to a useful block-triangular form, and the equations of motion for k-point Green’s functions form a closed hierarchy. Based on this formalism, results on criticality and dissipative phase transitions in such models are discussed in a companion paper.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Reversible Ferroelectric Polarization Modulation of Chiral Molecular Ferroelectrics by Circularly Polarized Light

Abstract The optical modulation of ferroelectric polarization constitutes a transformative, non‐contact strategy for the precise manipulation of ferroelectric properties, heralding advancements in optically stimulated ferroelectric devices. Despite its potential, progress in this domain is constrained by material limitations and the intricate nature of the underlying mechanisms. Recent studies have achieved efficient regulation of ferroelectric polarization and thermal conductivity in chiral ferroelectric thin films through the application of left‐ and right‐handed circularly polarized light (LCP and RCP). Differential absorption of circularly polarized light (CPL) induces nonequilibrium carrier dynamics, generating distinctive interfacial electrostatic fields that enable precise control of ultrathin ferroelectric films. For (R)‐BINOL−DIPASi and (S)‐BINOL−DIPASi (C 26 H 26 O 2 Si), polarization changes surpass 23%, exhibiting opposite response under LCP and RCP excitation. In R chiral films, remnant polarization decreases from 1.05 µC cm − 2 under LCP to 0.85 µC cm − 2 under RCP, whereas in S chiral films, polarization increases from 0.85 µC cm − 2 under LCP to 0.98 µC cm − 2 under RCP. This reversible modulation facilitates reliable switching between ON and OFF states, presenting the potential of chiral ferroelectric materials for flexible, high‐speed integrated photonic sensor technologies.

Chemistry↗
DOE OSTI2025

Field-Driven Simulations to Probe the Impact of Ionic Correlations on Solution Transport Coefficients in Binary, Ternary, and Reciprocal Quaternary Aqueous Electrolytes

Ionic correlations play a critical role in governing transport properties of mixed-salt aqueous electrolyte solutions, yet their quantitative characterization remains challenging, particularly for multicomponent aqueous systems. Here, we develop a general nonequilibrium molecular dynamics framework to efficiently compute Onsager transport coefficients and ionic correlations in mixed-salt aqueous solutions. Using field-driven simulations, we obtain accurate Onsager matrices for LiCl/ KCl, KCl/KBr, and LiBr/KCl electrolyte solutions with significantly reduced computational cost relative to equilibrium Green−Kubo methods. The framework enables direct assessment of how attractive cation−anion and repulsive like-ion interactions contribute to conductivity and salt diffusivity across compositions. While static ion pairing exhibits strong composition dependence, dynamic ion correlations remain nearly invariant, leading to constant deviations from Nernst−Einstein predictions. These results highlight the disconnect between static ion association and dynamic transport correlations, and they establish a transferable approach for analyzing ion transport in complex electrolyte environments relevant to separation processes and electrochemical systems.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Differentiating Thermal Conductances at Semiconductor Nanocrystal/Ligand and Ligand/Solvent Interfaces in Colloidal Suspensions

Infrared-pump, electronic-probe (IPEP) spectroscopy is used to measure heat flow into and out of CdSe nanocrystals suspended in an organic solvent, where the surface ligands are initially excited with an infrared pump pulse. Subsequently, the heat is transferred from the excited ligands to the nanocrystals and in parallel to the solvent. Parallel heat transfer in opposite directions uniquely enables us to differentiate the thermal conductances at the nanocrystal/ligand and ligand/solvent interfaces. Using a novel solution to the heat diffusion equation, we fit the IPEP data to find that the nanocrystal/ligand conductances range from 88 to 135 MW m –2 K –1 and are approximately 1 order of magnitude higher than the ligand/solvent conductances, which range from 7 to 26 MW m –2 K –1 . Transient nonequilibrium molecular dynamics (MD) simulations of nanocrystal suspensions agree with IPEP data and show that ligands bound to the nanocrystal by bidentate bonds have more than twice the per-ligand conductance as those bound by monodentate bonds.

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

A kinetic model of friction in strongly coupled strongly magnetized plasmas

Plasmas that are strongly magnetized in the sense that the gyrofrequency exceeds the plasma frequency exhibit novel transport properties that are not well understood. As a representative example, we compute the friction force acting on a massive test charge moving through a strongly coupled and strongly magnetized one-component plasma using a generalized Boltzmann kinetic theory. Recent works studying the weakly coupled regime have shown that strong magnetization leads to a transverse component of the friction force that is perpendicular to both the Lorentz force and velocity of the test charge, in addition to the stopping power component aligned antiparallel to the velocity. Recent molecular dynamics simulations have also shown that strong Coulomb coupling in addition to strong magnetization gives rise to a third component of the friction force in the direction of the Lorentz force. Here, we show that the generalized Boltzmann kinetic theory captures these effects and generally agrees well with the molecular dynamics simulations over a broad range of Coulomb coupling and magnetization strength regimes. The theory is also used to show that the “gyro” component of the friction in the direction of the Lorentz force arises due to asymmetries associated with gyromotion during short-range collisions. Computing the average motion of the test charge through the background plasma, the transverse force is found to strongly influence the trajectory by changing the gyroradius and the gyrofriction force is found to slightly change the gyrofrequency of the test charge resulting in a phase shift.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗