Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nonequilibrium systems”

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

Ultrafast dynamics of a fermion chain in a terahertz field-driven optical cavity

In this article, we study the effect of a terahertz field-driven single cavity mode for ultrafast control of a fermion chain with dissipation-induced nonlinearity and quadratic coupling to an infrared-active phonon mode. Unlike the first-order phase transition in the nonequilibrium steady state of the system without a cavity or with a high-quality cavity, we find that a realistic dissipation in a low-quality driven cavity process prevents such a transition. Without realistic photon loss from the cavity, the transition only survives for the lower phonon-polariton branch with strong drive; a weak laser field fails to induce the phase transition and renders the polaritons symmetrical. The ability to control the phase transition is crucial for realizing strongly modulated steady states; we propose experimentally feasible regimes where this occurs.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

$${\mathcal{P}}{\mathcal{T}}$$-Symmetric Effective Model for Nonequilibrium Phase Transitions in a Dissipative Fermionic Mott Insulator Chain

Nonequilibrium phase transitions in open dissipative systems can be described as instabilities in the spectra and wavefunctions of effective non-Hermitian Hamiltonians invariant under simultaneous parity (P) and time-reversal (T) transformations. The degree of non-Hermiticity reflects the strength of the external drive and dissipation, and the transition is described as a loss of the PT symmetry of the solutions corresponding to stationary low-drive dynamics. This approach has been successfully applied to spin, superconducting, and Mott insulator systems. However, the microscopic foundations for the employed phenomenological models are currently lacking. Here we propose a microscopic mechanism leading to the PT-symmetric effective model in the context of the nonequilibrium Mott transition in a dissipative Hubbard chain. Our model comprises a half-filled fermionic Hubbard chain subject to a constant electric field. The dissipation is introduced via the electron-phonon coupling. We obtain the explicit expressions for the non-Hermitian parameter in terms of the electron-phonon coupling strength and driving field. Analyzing the implications of microscopic model, we find a re-entrant Mott insulator with the increasing electric field for phonon density of states that increases slower than the square of the energy (such as in one or two dimensions), or varies non-monotonously with energy.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Pump-induced motion of an interface between competing orders

In this work, we study the motion of an interface separating two regions with different electronic orders following a short-duration pump that drives the system out of equilibrium. Using a generalized Ginzburg-Landau approach and assuming that the main effect of the nonequilibrium drive is to transiently heat the system we address the question of the direction of interface motion—in other words, which ordered region expands and which contracts after the pump. Our analysis includes the effects of differences in free-energy landscape and in order parameter dynamics and identifies circumstances in which the drive may act to increase the volume associated with the subdominant order, for example when the subdominant order has a second-order free-energy landscape while the dominant order has a first-order one.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Unified, Geometric Framework for Nonequilibrium Protocol Optimization

Controlling thermodynamic cycles to minimize the dissipated heat is a long-standing goal in thermodynamics, and more recently, a central challenge in stochastic thermodynamics for nanoscale systems. Here, we introduce a theoretical and computational framework for optimizing nonequilibrium control protocols that can transform a system between two distributions in a minimally dissipative fashion. These protocols optimally transport a system along paths through the space of probability distributions that minimize the dissipative cost of a transformation. Furthermore, we show that the thermodynamic metric—determined via a linear response approach—can be directly derived from the same objective function that is optimized in the optimal transport problem, thus providing a unified perspective on thermodynamic geometries. As a result, we investigate this unified geometric framework in two model systems and observe that our procedure for optimizing control protocols is robust beyond linear response.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Transport, correlations, and chaos in a classical disordered anharmonic chain

Here, we explore transport properties in a disordered nonlinear chain of classical harmonic oscillators, and thereby identify a regime exhibiting behavior analogous to that seen in quantum many-body-localized systems. Through extensive numerical simulations of this system connected at its ends to heat baths at different temperatures, we computed the heat current and the temperature profile in the nonequilibrium steady state as a function of system size $\textit{N}$, disorder strength Δ, and temperature $\textit{T}$. The conductivity $κ_N$, obtained for finite length ($\textit{N}$), saturates to a value $κ_∞$ > 0 in the large $\textit{N}$ limit, for all values of disorder strength Δ and temperature $\textit{T}$ > 0. We show evidence that for any Δ > 0 the conductivity goes to zero faster than any power of $\textit{T}$ in the ($\textit{T}$/Δ) → 0 limit, and find that the form $κ_∞$ ~ $e^{{–B|\text{ln}(CΔ/T)|}^{3}}$ fits our data. This form has earlier been suggested by a theory based on the dynamics of multioscillator chaotic islands. The finite-size effect can be $κ_N < κ_∞$ due to boundary resistance when the bulk conductivity is high (the weak disorder case), or $κ_N > κ_∞$ due to direct bath-to-bath coupling through bulk localized modes when the bulk is weakly conducting (the strong disorder case). We also present results on equilibrium dynamical correlation functions and on the role of chaos on transport properties. Finally, we explore the differences in the growth and propagation of chaos in the weak and strong chaos regimes by studying the classical version of the out-of-time-ordered commutator.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Active Colloids as Models, Materials, and Machines

Active colloids use energy input at the particle level to propel persistent motion and direct dynamic assemblies. We consider three types of colloids animated by chemical reactions, time-varying magnetic fields, and electric currents. For each type, we review the basic propulsion mechanisms at the particle level and discuss their consequences for collective behaviors in particle ensembles. These microscopic systems provide useful experimental models of nonequilibrium many-body physics in which dissipative currents break time-reversal symmetry. Freed from the constraints of thermodynamic equilibrium, active colloids assemble to form materials that move, reconfigure, heal, and adapt. Colloidal machines based on engineered particles and their assemblies provide a basis for mobile robots with increasing levels of autonomy. This review provides a conceptual framework for understanding and applying active colloids to create material systems that mimic the functions of living matter. We highlight opportunities for chemical engineers to contribute to this growing field.

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

Superconducting-like response in driven systems near the Mott transition

Here we point out that fractionalized bosonic charge excitations can explain the recently discovered photoinduced superconducting-like response in $\kappa$-(ET) 2 Cu [N(CN) 2 ] Br, an organic metal close to the Mott transition. The pump laser exerts a periodic drive on the fractionalized field, creating a nonequilibrium condensate, which gives a Drude peak much narrower than the equilibrium scattering rate, hence a superconducting-like response. Our proposal illuminates alternative possibilities of detecting fractionalization and can be readily tested in spin liquid candidates and in cold atom systems.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Influence of Markovianity and self-consistency on time-resolved spectral functions of driven quantum systems

We present a systematic comparison of the real-time Dyson expansion (RTDE) with established nonequilibrium Green's function (GF) approaches for simulating driven, interacting quantum systems. Focusing on density matrix dynamics, time-off-diagonal GFs, and time-resolved photoemission spectra, we benchmark RTDE against fully self-consistent Kadanoff-Baym equation (KBE) calculations, the generalized Kadanoff-Baym ansatz, and exact diagonalization for small systems using second-order many-body perturbation theory. Using a driven two-band Hubbard model, we show that mean-field single-particle density matrix trajectories provide a reliable baseline for RTDE across a broad range of interaction strengths and excited-carrier populations. Further, RTDE accurately captures correlation effects in the GFs, including long-lived oscillations and revivals that are strongly suppressed by the overdamping inherent to self-consistent KBE schemes. As a consequence, RTDE resolves rich nonequilibrium spectral structure in time-resolved photoemission, such as interaction- and population-dependent quasiparticle splittings and band gap renormalization, which are largely washed out in self-consistent approaches yet are present in exact solutions. Furthermore, our results demonstrate that RTDE bridges the gap between mean-field propagation and full two-time KBE simulations, retaining favorable linear scaling while capturing essential dynamical correlations relevant for ultrafast spectroscopy.

Electronic structure↗
DOE OSTI2024

Robust 3D multi-material hydrodynamics using discontinuous Galerkin methods

A high-order discontinuous Galerkin (DG) method is presented for nonequilibrium multi-material (m ≥ 2) flow with sharp interfaces. Material interfaces are reconstructed using the algebraic THINC approach, resulting in a sharp interface resolution. The system assumes stiff velocity relaxation and pressure nonequilibrium. The presented DG method uses Dubiner's orthogonal basis functions on tetrahedral elements. This results in a unique combination of sharp multimaterial interfaces and high-order accurate solutions in smooth single-material regions. A novel shock indicator based on the interface conservation condition is introduced to mark regions with discontinuities. Slope limiting techniques are applied only in these regions so that nonphysical oscillations are eliminated while maintaining high-order accuracy in smooth regions. A local projection is applied on the limited solution to ensure discrete closure law preservation. The effectiveness of this novel limiting strategy is demonstrated for complex three-dimensional multi-material problems, where robustness of the method is critical. The presented numerical problems demonstrate that more accurate and efficient multi-material solutions can be obtained by the DG method, as compared to second-order finite volume methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Energy-conserving contact dynamics of nonspherical rigid-body particles

Understanding the contact dynamics of nonspherical particles is crucial for accurately modeling colloidal and granular systems where shape anisotropy dictates structural organization and transport properties. We here introduce an energy-conserving contact dynamics framework for arbitrary convex rigid-body particles by implementing vertex–boundary interactions in 2D and vertex–surface and edge–edge detection in 3D. The established formulation enables continuous force evaluation and prevents particle overlap while conserving total energy during translational and rotational motion. We demonstrate the framework’s stability and its utility to capture packing behavior, anisotropic diffusion, and equations of state of polygonal and polyhedral particles as examples. The framework establishes a robust and extensible foundation for investigating nonequilibrium dynamics of complex nonspherical particulate systems, enabling enhanced understanding of phenomena across spatiotemporal scales in self- and directed-assembly, granular flows, and hydrodynamics, potentially coupled with interactions that represent underlying physical mechanisms.

Discrete element method↗
DOE OSTI2022

Anomalous phase separation in a correlated electron system: Machine-learning–enabled large-scale kinetic Monte Carlo simulations

Significance Phase separation is crucial to the functionalities of many correlated electron materials with notable examples including colossal magnetoresistance in manganites and high- T c superconductivity in cuprates. However, the nonequilibrium phase-separation dynamics in such systems are poorly understood theoretically, partly because the required multiscale modeling is computationally very demanding. With the aid of machine-learning methods, we have achieved large-scale dynamical simulations in a representative correlated electron system. We observe an unusual relaxation process that is beyond the framework of classical phase-ordering theories. We also uncover a correlation-induced freezing behavior, which could be a generic feature of phase separation in correlated electron systems.

Zhang, Sheng↗
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

Ergodicity breaking in rapidly rotating C 60 fullerenes

Ergodicity, the central tenet of statistical mechanics, requires an isolated system to explore all available phase space constrained by energy and symmetry. Mechanisms for violating ergodicity are of interest for probing nonequilibrium matter and protecting quantum coherence in complex systems. Polyatomic molecules have long served as a platform for probing ergodicity breaking in vibrational energy transport. Here, we report the observation of rotational ergodicity breaking in an unprecedentedly large molecule, 12 C 60 , determined from its icosahedral rovibrational fine structure. The ergodicity breaking occurs well below the vibrational ergodicity threshold and exhibits multiple transitions between ergodic and nonergodic regimes with increasing angular momentum. These peculiar dynamics result from the molecule’s distinctive combination of symmetry, size, and rigidity, highlighting its relevance to emergent phenomena in mesoscopic quantum systems.

Science & Technology - Other Topics↗
DOE OSTI2024

Probabilistic measures for biological adaptation and resilience

This paper introduces an approach to quantifying ecological resilience in biological systems, particularly focusing on noisy systems responding to episodic disturbances with sudden adaptations. Incorporating concepts from nonequilibrium statistical mechanics, we propose a measure termed “ecological resilience through adaptation,” specifically tailored to noisy, forced systems that undergo physiological adaptation in the face of stressful environmental changes. Randomness plays a key role, accounting for model uncertainty and the inherent variability in the dynamical response among components of biological systems. Our measure of resilience is rooted in the probabilistic description of states within these systems and is defined in terms of the dynamics of the ensemble average of a model-specific observable quantifying success or well-being. Our approach utilizes stochastic linear response theory to compute how the expected success of a system, originally in statistical equilibrium, dynamically changes in response to a environmental perturbation and a subsequent adaptation. Importantly, the resulting mathematical derivations allow for the estimation of resilience in terms of ensemble averages of simulated or experimental data. Finally, through a simple but clear conceptual example, we illustrate how our resilience measure can be interpreted and compared to other existing frameworks in the literature. The methodology is general but inspired by applications in plant systems, with the potential for broader application to complex biological processes.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Optimizing Ammonia Separation via Reactive Absorption for Sustainable Ammonia Synthesis

Metal halide salts such as magnesium chloride have been shown to be promising candidates for ammonia storage materials for energy storage and agriculture applications due to their ability to incorporate several moles of ammonia per mole of salt. Ammonia exiting a synthesis reactor can be separated from nitrogen and hydrogen by absorption into magnesium chloride. Such an absorption can be more complete and hotter than separation via ammonia condensation, the current standard in the Haber–Bosch process. Here, we discuss the optimal conditions for the cyclic uptake and release of ammonia from the supported magnesium chloride absorbents. An automated system was designed for measuring the nonequilibrium working capacity of the absorbent, as well as the impact of important operating conditions such as absorption and desorption temperature, pressure, and desorption time. Measurements of absorption and desorption kinetics provide insight into the mechanisms involved. The temperatures and pressures during absorption and desorption were designed to use minimal energy input to maximize the uptake and release of ammonia within a reasonable amount of time. In a laboratory-scale bed, absorption has a small unused bed length, so it is largely independent of mass transfer; it is dominated by how fast ammonia is fed to the bed. Yet, desorption is restricted both by the speed of heating the bed and by diffusion out of the absorbent. These measurements provide guidelines for ammonia separations and cycling sorbent materials on a larger scale.

10 SYNTHETIC FUELS↗
DOE OSTI2022

Reversible to irreversible transitions for cyclically driven particles on periodic obstacle arrays

Here, we examine the collective dynamics of disks moving through a square array of obstacles under cyclic square wave driving. Below a critical density, we find that the system organizes into a reversible state in which the disks return to the same positions at the end of every drive cycle. Above this density, the dynamics are irreversible and the disks do not return to the same positions after each cycle. The critical density depends strongly on the angle θ between the driving direction and a symmetry axis of the obstacle array, with the highest critical densities appearing at commensurate angles such as θ = 0° and θ = 45° and the lowest critical densities falling at θ = arctan (0.618), the inverse of the golden ratio, where the flow is the most degenerate. As the density increases, the number of cycles required to reach a reversible state grows as a power law with an exponent near ν = 1.36, similar to what is found in periodically driven colloidal and superconducting vortex systems.

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

Steady-state properties of multi-orbital systems using quantum Monte Carlo

A precise dynamical characterization of quantum impurity models with multiple interacting orbitals is challenging. In quantum Monte Carlo methods, this is embodied by sign problems. A dynamical sign problem makes it exponentially difficult to simulate long times. A multi-orbital sign problem generally results in a prohibitive computational cost for systems with multiple impurity degrees of freedom even in static equilibrium calculations. Here, we present a numerically exact inchworm method that simultaneously alleviates both sign problems, enabling simulation of multi-orbital systems directly in the equilibrium or nonequilibrium steady-state. The method combines ideas from the recently developed steady-state inchworm Monte Carlo framework [Erpenbeck et al., Phys. Rev. Lett. 130, 186301 (2023)] with other ideas from the equilibrium multi-orbital inchworm algorithm [Eidelstein et al., Phys. Rev. Lett. 124, 206405 (2020)]. We verify our method by comparison with analytical limits and numerical results from previous methods.

Chemistry↗
DOE OSTI2026

Field-driven ion pairing dynamics in concentrated electrolytes

Here, we investigate ion pairing dynamics in electrolytes driven far from equilibrium using molecular simulations and nonequilibrium rate theory. Focusing on 0.5M LiPF 6 in water and acetonitrile under uniform electric fields, we compute transition path theory observables, including reactive fluxes and mean first-passage times of ion pairing. Moreover, we introduce a dynamical proxy of free-ion population, where its field-induced change is strongly correlated with the nonlinear enhancement of conductivity, yielding an increase of 40% at 50 mV/Å in acetonitrile, compared to that of less than 10% in aqueous electrolytes. Further kinetic analysis elucidates that Onsager’s classical theory substantially overestimates field-induced enhancement of ion pair dissociation in molecular electrolytes. This discrepancy arises from solvent-mediated dynamical pathways and field-induced dielectric decrement that suppress ion pair dissociation within explicit solvents, highlighting that a faithful description of molecular details is essential. Our results provide a molecular interpretation of nonlinear electrolyte transport beyond continuum theories and establish a general framework for quantifying nonequilibrium reaction kinetics in condensed phase systems.

Moon, Seokjin [University of California, Berkeley,↗