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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
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 OSTI2020

A Machine Learning Assisted Development of a Model for the Populations of Convective and Stratiform Clouds

Abstract Traditional parameterizations of the interaction between convection and the environment have relied on an assumption that the slowly varying large‐scale environment is in statistical equilibrium with a large number of small and short‐lived convective clouds. They fail to capture nonequilibrium transitions such as the diurnal cycle and the formation of mesoscale convective systems as well as observed precipitation statistics and extremes. Informed by analysis of radar observations, cloud‐permitting model simulation, theory, and machine learning, this work presents a new stochastic cloud population dynamics model for characterizing the interactions between convective and stratiform clouds, with the goal of informing the representation of these interactions in global climate models. Fifteen wet seasons of precipitating cloud observations by a C‐band radar at Darwin, Australia are fed into a machine learning algorithm to obtain transition functions that close a set of coupled equations relating large‐scale forcing, mass flux, the convective cell size distribution, and the stratiform area. Under realistic large‐scale forcing, the derived transition functions show that, on the one hand, interactions with stratiform clouds act to dampen the variability in the size and number of convective cells and therefore in the convective mass flux. On the other, for a given convective area fraction, a larger number of smaller cells is more favorable for the growth of stratiform area than a smaller number of larger cells. The combination of these two factors gives rise to solutions with a few convective cells embedded in a large stratiform area, reminiscent of mesoscale convective systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Collapse of the waiting time effect in a spin glass

Waiting time measurements are reported here for a bulk polycrystalline sample of the spin glass Cu 0.95 Mn 0.05 . We observe a suppression of aging as the glass temperature is approached from below. These results are interpreted in terms of a spin-glass correlation length approaching and reaching crystallite sizes. They provide a direct experimental connection between the correlation length scale and relaxation times in the spin-glass state. In particular, these measurements track the value of the dynamic critical exponent as the system transitions from the critical temperature regime to its low-temperature, quasistatic nonequilibrium phase.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Overcoming immiscibility toward bimetallic catalyst library

Bimetallics are emerging as important materials that often exhibit distinct chemical properties from monometallics. However, there is limited access to homogeneously alloyed bimetallics because of the thermodynamic immiscibility of the constituent elements. Overcoming the inherent immiscibility in bimetallic systems would create a bimetallic library with unique properties. Here, we present a nonequilibrium synthesis strategy to address the immiscibility challenge in bimetallics. As a proof of concept, we synthesize a broad range of homogeneously alloyed Cu-based bimetallic nanoparticles regardless of the thermodynamic immiscibility. The nonequilibrated bimetallic nanoparticles are further investigated as electrocatalysts for carbon monoxide reduction at commercially relevant current densities (>100 mA cm -2 ), in which Cu 0.9 Ni 0.1 shows the highest multicarbon product Faradaic efficiency of ~76% with a current density of ~93 mA cm -2 . The ability to overcome thermodynamic immiscibility in multimetallic synthesis offers freedom to design and synthesize new functional nanomaterials with desired chemical compositions and catalytic properties.

25 ENERGY STORAGE↗
DOE OSTI2020

ACCESS: Confirmation of No Potassium in the Atmosphere of WASP-31b

We present a new optical (400–950 nm) transmission spectrum of the hot Jupiter WASP-31b (M = 0.48 M {sub J}; R = 1.54 R {sub J}; P = 3.41 days), obtained by combining four transit observations. These transits were observed with IMACS on the Magellan Baade Telescope at Las Campanas Observatory as part of the ACCESS project. We investigate the presence of clouds/hazes in the upper atmosphere of this planet, as well as the contribution of stellar activity on the observed features. In addition, we search for absorption features of the alkali elements Na i and K i, with particular focus on K i, for which there have been two previously published disagreeing results. Observations with Hubble Space Telescope (HST)/STIS detected K i, whereas ground-based low- and high-resolution observations did not. We use equilibrium and nonequilibrium chemistry retrievals to explore the planetary and stellar parameter space of the system with our optical data combined with existing near-IR observations. Our best-fit model is that with a scattering slope consistent with a Rayleigh slope (α=5.3{sub −3.1}{sup +2.9}), high-altitude clouds at a log cloud top pressure of −3.6{sub −2.1}{sup +2.7} bars, and possible muted H{sub 2}O features. We find that our observations support other ground-based claims of no K i. Clouds are likely why signals like H{sub 2}O are extremely muted and Na or K cannot be detected. We then juxtapose our Magellan/IMACS transmission spectrum with existing VLT/FORS2, HST/WFC3, HST/STIS, and Spitzer observations to further constrain the optical-to-infrared atmospheric features of the planet. We find that a steeper scattering slope (α = 8.3 ± 1.5) is anchored by STIS wavelengths blueward of 400 nm and only the original STIS observations show significant potassium signal.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Protein–Ligand Binding Free-Energy Calculations with ARROW–A Purely First-Principles Parameterized Polarizable Force Field

Protein–ligand binding free-energy calculations using molecular dynamics (MD) simulations have emerged as a powerful tool for in silico drug design. Here, we present results obtained with the ARROW force field (FF)–a multipolar polarizable and physics-based model with all parameters fitted entirely to high-level ab initio quantum mechanical (QM) calculations. ARROW has already proven its ability to determine solvation free energy of arbitrary neutral compounds with unprecedented accuracy. The ARROW FF parameterization is now extended to include coverage of all amino acids including charged groups, allowing molecular simulations of a series of protein–ligand systems and prediction of their relative binding free energies. We ensure adequate sampling by applying a novel technique that is based on coupling the Hamiltonian Replica exchange (HREX) with a conformation reservoir generated via potential softening and nonequilibrium MD. ARROW provides predictions with near chemical accuracy (mean absolute error of ~0.5 kcal/mol) for two of the three protein systems studied here (MCL1 and Thrombin). The third protein system (CDK2) reveals the difficulty in accurately describing dimer interaction energies involving polar and charged species. Overall, for all of the three protein systems studied here, ARROW FF predicts relative binding free energies of ligands with a similar accuracy level as leading nonpolarizable force fields.

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

Real-Time Time-Dependent Density Functional Theory for Simulating Nonequilibrium Electron Dynamics

The explicit real-time propagation approach for time-dependent density functional theory (RT-TDDFT) has increasingly become a popular first-principles computational method for modeling various time-dependent electronic properties of complex chemical systems. Here, in this Perspective, we provide a nontechnical discussion of how this first-principles simulation approach has been used to gain novel physical insights into nonequilibrium electron dynamics phenomena in recent years. Following a concise overview of the RT-TDDFT methodology from a practical standpoint, we discuss our recent studies on the electronic stopping of DNA in water and the Floquet topological phase as examples. Our discussion focuses on how RT-TDDFT simulations played a unique role in deriving new scientific understandings. We then discuss existing challenges and some new advances at the frontier of RT-TDDFT method development for studying increasingly complex dynamic phenomena and systems.

Hamiltonians↗
DOE OSTI2021

Inducing thermodynamically blocked atomic ordering via strongly driven nonequilibrium kinetics

Ultrafast light-matter interactions enable inducing exotic material phases by promoting access to kinetic processes blocked in equilibrium. Despite potential opportunities, actively using nonequilibrium kinetics for material discovery is limited by the poor understanding on intermediate states of driven systems. Here, using single-pulse time-resolved imaging with x-ray free-electron lasers, we found intermediate states of photoexcited bismuth nanoparticles that showed kinetically reversed surface ordering during ultrafast melting. This entropy-lowering reaction was further investigated by molecular dynamics simulations to reveal that observed kinetics were thermodynamically buried in equilibrium, which emphasized the critical role of electron-mediated ultrafast free-energy modification in inducing exotic material phases. This study demonstrated that ultrafast photoexcitations of electrons provide an efficient strategy to induce hidden material phases by overcoming thermodynamic barriers via nonequilibrium reaction pathways.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Hyperuniform Active Chiral Fluids with Tunable Internal Structure

Large density fluctuations observed in active systems and hyperuniformity are two seemingly incompatible phenomena. However, the formation of hyperuniform states has been recently predicted in nonequilibrium fluids formed by chiral particles performing circular motion with the same handedness. In this work, we report evidence of hyperuniformity realized in a chiral active fluid comprised of pear-shaped Quincke rollers of arbitrary handedness. We show that hyperunifonnity and large density fluctuations, triggered by dynamic clustering, coexist in this system at different length scales. The system loses its hyperuniformity as the curvature of particles' motion increases, transforming them into localized spinners. Our results experimentally demonstrate a novel hyperuniform active fluid and provide new insights into an interplay between chirality, activity, and hyperuniformity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Strong Hilbert space fragmentation and fractons from subsystem and higher-form symmetries

Here, we introduce a route to Hilbert space fragmentation in high dimensions leveraging the group-word formalism. We show that taking strongly fragmented models in one dimension and “lifting” to higher dimensions using subsystem symmetries can yield strongly fragmented dynamics in higher dimensions, with subdimensional (e.g., lineonic) excitations. This provides a route to higher-dimensional strong fragmentation, and also a route to fractonic behavior. Meanwhile, lifting one-dimensional fragmented models to higher dimensions using higher-form symmetries yields models with topologically robust fragmentation. In three or more spatial dimensions, one can also “mix and match” subsystem and higher-form symmetries, leading to canonical fracton models such as X cube. We speculate that this approach could also yield a route to non-Abelian fractons. These constructions unify a number of phenomena that have been discussed in the literature, as well as furnishing models with unique properties.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Minimal two-body quantum absorption refrigerator

Here, we study the phenomenon of absorption refrigeration, where refrigeration is achieved by heating instead of work, in two different setups: a minimal set up based on coupled qubits and two nonlinearly coupled resonators. Considering ZZ interaction between the two qubits, we outline the basic ingredients required to achieve cooling. Using local as well as global master equations, we observe that inclusion of an XX type term in the qubit-qubit coupling is detrimental to cooling. We compare the cooling effect obtained in the qubit case with that of nonlinearly coupled resonators (multilevel system) where the ZZ interaction translates to a Kerr-type nonlinearity. For small to intermediate strengths of nonlinearity, we observe that multilevel quantum systems, for example qutrits, give better cooling effect compared to the qubits. Using Keldysh nonequilibrium Green's function formalism, we go beyond first order sequential tunneling processes and study the effect of higher order processes on refrigeration. We find a reduced cooling effect compared to the master equation calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

On Gibbs Equilibrium and Hillert Nonequilibrium Thermodynamics

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.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Scrambling Transition in a Radiative Random Unitary Circuit

Here we study quantum information scrambling in a random unitary circuit that exchanges qubits with an environment at a rate p. As a result, initially localized quantum information not only spreads within the system, but also spills into the environment. Using the out-of-time-order correlator (OTOC) to characterize scrambling, we find a nonequilibrium phase transition in the directed percolation universality class at a critical swap rate p c : for p < p c the ensemble-averaged OTOC exhibits ballistic growth with a tunable light cone velocity, while for p > p c the OTOC fails to percolate within the system and vanishes uniformly within a finite timescale, indicating that all local operators are rapidly swapped into the environment. To elucidate its information-theoretic consequences, we demonstrate that the transition in operator spreading coincides with a transition in an observer’s ability to decode the system’s initial quantum information from the swapped-out, or “radiated,” qubits. We present a simple decoding scheme which recovers the system’s initial information with perfect fidelity in the nonpercolating phase, and with continuously decreasing fidelity with decreasing swap rate in the percolating phase. Depending on the initial state of the swapped-in qubits, we further observe a corresponding entanglement transition in the coherent information from the system into the radiated qubits.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A Review: Microstructural and Phase Evolution in Alloys during Extended Plastic Deformation

Materials in service and during processing are often subjected to plastic deformation. For multi-phase metallic alloys, simple geometric models and atomistic simulations suggest that two distinctive regimes in these materials’ evolution during deformation exist. At low strains, evolutions are often dominated by the kinetic roughening of interfaces, which results from the superdiffusive transport of matter in sheared crystals. Furthermore, at high strains and temperatures where thermal diffusion is sluggish, on the other hand, shearing- induced forced atomic mixing dominates these evolutions, resulting in significantly enhanced solubility. Distinguishing these two regimes is shown to provide a convenient framework for rationalizing and analyzing recent experiments and simulations on wear of layered structures in the low strain-regime and nonequilibrium phase co-existence and self-organization in highly immiscible or reactive alloy systems in the high-strain regime.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Computing equilibrium free energies through a nonequilibrium quench

Many methods to accelerate sampling of molecular configurations are based on the idea that temperature can be used to accelerate rare transitions. These methods typically compute equilibrium properties at a target temperature using reweighting or through Monte Carlo exchanges between replicas at higher temperatures. A recent paper [G. M. Rotskoff and E. Vanden-Eijnden, Phys. Rev. Lett. 122, 150602 (2019)] demonstrated that accurate equilibrium densities of states can also be computed through a nonequilibrium “quench” process, where sampling is performed at a higher temperature to encourage rapid mixing and then quenched to lower energy states with dissipative dynamics. Here, we provide an implementation of the quench dynamics in LAMMPS and evaluate a new formulation of nonequilibrium estimators for the computation of partition functions or free energy surfaces (FESs) of molecular systems. We show that the method is exact for a minimal model of N-independent harmonic springs and use these analytical results to develop heuristics for the amount of quenching required to obtain accurate sampling. We then test the quench approach on alanine dipeptide, where we show that it gives an FES that is accurate near the most stable configurations using the quench approach but disagrees with a reference umbrella sampling calculation in high FE regions. We then show that combining quenching with umbrella sampling allows the efficient calculation of the free energy in all regions. Moreover, by using this combined scheme, we obtain the FES across a range of temperatures at no additional cost, making it much more efficient than standard umbrella sampling if this information is required. Finally, we discuss how this approach can be extended to solute tempering and demonstrate that it is highly accurate for the case of solvated alanine dipeptide without any additional modifications.

Chemistry↗
DOE OSTI2024

Local strain inhomogeneities during electrical triggering of a metal–insulator transition revealed by X-ray microscopy

Electrical triggering of a metal–insulator transition (MIT) often results in the formation of characteristic spatial patterns such as a metallic filament percolating through an insulating matrix or an insulating barrier splitting a conducting matrix. When MIT triggering is driven by electrothermal effects, the temperature of the filament or barrier can be substantially higher than the rest of the material. Using X-ray microdiffraction and dark-field X-ray microscopy, we show that electrothermal MIT triggering leads to the development of an inhomogeneous strain profile across the switching device, even when the material does not undergo a pronounced, discontinuous structural transition coinciding with the MIT. Diffraction measurements further reveal evidence of unique features associated with MIT triggering including lattice distortions, tilting, and twinning, which indicate structural nonuniformity of both low- and high-resistance regions inside the switching device. Such lattice deformations do not occur under equilibrium, zero-voltage conditions, highlighting the qualitative difference between states achieved through increasing temperature and applying voltage in nonlinear electrothermal materials. Electrically induced strain, lattice distortions, and twinning could have important contributions in the MIT triggering process and drive the material into nonequilibrium states, providing an unconventional pathway to explore the phase space in strongly correlated electronic systems.

42 ENGINEERING↗
DOE OSTI2021

Ultrafast Neuromorphic Dynamics Using Hidden Phases in the Prototype of Relaxor Ferroelectrics

Materials possessing multiple states are promising to emulate synaptic and neuronic behaviors. Their operation frequency, typically in or below the GHz range, however, limits the speed of neuromorphic computing. Ultrafast THz electric field excitation has been employed to induce nonequilibrium states of matter, called hidden phases in oxides. One may wonder if there are systems for which THz pulses can generate neuronic and synaptic behavior, via the creation of hidden phases. Using atomistic simulations, we discover that relaxor ferroelectrics can emulate all the key neuronic and memristive synaptic features. Their occurrence originates from the activation of many hidden phases of polarization order, resulting from the response of nanoregions to THz pulses. We find such phases further possess different dielectric constants, which is also promising for memcapacitor devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Quenching through the QCD chiral phase transition

We present a detailed numerical and analytical study of the out-of-equilibrium dynamics of Model G, the dynamical universality class relevant to the chiral phase transition. We perform numerical 3D stochastic (Langevin) simulations of the 𝑂⁡(4) critical point for large lattices in the chiral limit. We quench the system from the high-temperature unbroken phase to the broken phase and study the nonequilibrium dynamics of pion fields. Strikingly, the nonequilibrium evolution of the two-point functions exhibits a regime of growth, a parametrically large enhancement, and a subsequent slow relaxation to equilibrium. We analyze our numerical results using dynamic critical scaling and mean-field theory. The growth of the two-point functions is determined by the nonlinear dynamics of an ideal non-Abelian superfluid, which is a limit of Model G that reflects the broken chiral symmetry. We also relate the nonequilibrium two-point functions to a long-lived parametric enhancement of soft pion yields relative to thermal equilibrium following a quench.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗