Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nonequilibrium molecular 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 OSTI2020

Transport phenomena in electrolyte solutions: Nonequilibrium thermodynamics and statistical mechanics

The theory of transport phenomena in multicomponent electrolyte solutions is presented here through the integration of continuum mechanics, electromagnetism, and nonequilibrium thermodynamics. The governing equations of irreversible thermodynamics, including balance laws, Maxwell's equations, internal entropy production, and linear laws relating the thermodynamic forces and fluxes, are derived. Green–Kubo relations for the transport coefficients connecting electrochemical potential gradients and diffusive fluxes are obtained in terms of the flux–flux time correlations. Furthermore, the relationship between the derived transport coefficients and those of the Stefan–Maxwell and infinitely dilute frameworks are presented, and the connection between the transport matrix and experimentally measurable quantities is described. To exemplify the application of the derived Green–Kubo relations in molecular simulations, the matrix of transport coefficients for lithium and chloride ions in dimethyl sulfoxide is computed using classical molecular dynamics and compared with experimental measurements.

42 ENGINEERING↗
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↗
NASA NTRS1991

Algorithm Computes Hypersonic Flow Of Air

Algorithm simulates three-dimensional hypersonic flow of air or another mixture of gases around blunt body. Accounts for nonequilibrium thermochemical effects via two-temperature mathematical model of molecular excitations and equations for finite-rate chemical reactions. These equations fully coupled to fluid dynamical equations.

Palmer, Grant↗
DOE OSTI2020

Nonequilibrium Solvent Effects during Photodissociation in Liquids: Dynamical Energy Surfaces, Caging and Chemical Identity

In the gas phase, potential energy surfaces can be used to provide insight into the details of photochemical reaction dynamics. In solution, however, it is unclear what potential energy surfaces, if any, can be used to describe even simple chemical reactions such as the photodissociation of a diatomic solute. In this paper, we use mixed quantum/classical (MQC) molecular dynamics (MD) to study the photodissociation of $Na$ $^+_2$ in both liquid Ar and liquid tetrahydrofuran (THF). We examine both the gas-phase potential surfaces and potentials of mean force (PMF), which assume that the solvent remains at equilibrium with the solute throughout the photodissociation process and show that neither resemble a nonequilibrium dynamical energy surface that is generated by taking the time integral of work. For the photodissociation of $Na$ $^+_2$ in liquid Ar, the dynamical energy surface shows clear signatures of solvent caging, and the degree of caging is directly related to the mass of the solvent atoms. For $Na$ $^+_2$ in liquid THF, local specific interactions between the solute and solvent lead to changes in chemical identity that create a kinetic trap that effectively prevents the molecule from dissociating. Here, the results show that nonequilibrium effects play an important role even in simple solution-phase reactions, requiring the use of dynamical energy surface to understand such chemical events.

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

Observation of Fundamental Mechanisms in Compression-Induced Phase Transformations Using Ultrafast X-ray Diffraction

As theoretically hypothesized for several decades in group IV transition metals, we have discovered a dynamically stabilized body-centered cubic (bcc) intermediate state in Zr under uniaxial loading at sub-nanosecond timescales. Under ultrafast shock wave compression, rather than the transformation from α-Zr to the more disordered hex-3 equilibrium ω-Zr phase, in its place we find the formation of a previously unobserved nonequilibrium bcc metastable intermediate. We probe the compression-induced phase transition pathway in zirconium using time-resolved sub-picosecond x-ray diffraction analysis at the Linac Coherent Light Source. We also present molecular dynamics simulations using a potential derived from first-principles methods which independently predict this intermediate phase under ultrafast shock conditions. In contrast with experiments on longer timescale (> 10 ns) where the phase diagram alone is an adequate predictor of the crystalline structure of a material, our recent study highlights the importance of metastability and time dependence in the kinetics of phase transformations.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Single-shot ultrafast dynamics of nanosecond pulsed plasmas: Transition from ps–ns nonequilibrium to near-full ionization

The ultrafast multi-stage evolution of state-defining properties in atmospheric-pressure nanosecond pulsed plasmas is quantified using single-discharge, jitter-free, continuous streak-sweep spectroscopy of N 2 (C → B) molecular and N + /O + ionic species with a single-shot time resolution as short as ∼160 ps. An early period of extreme nonequilibrium is observed with vibrational temperatures [T V (C)] dropping from ∼8000 K at <200 ps after plasma breakdown to <4000 K within ∼1 ns, with near-ambient rotational temperatures [T R (C)] ∼ 300 K due to limited collisional energy transfer. This early-time T V (C) trend is representative of a direct and unquenched look into the high-energy tail of the electron energy distribution function; thus, it tracks real-time changes in the mean electron energy via the N2(C) emission signatures. This is followed by the rapid onset of N + /O + ionic emission after a distinct time delay of ∼14.9 ns due to multistep chemical kinetics. The ionic emission enables determination of electron densities (n e ) >2 × 10 19 cm −3 and electron temperatures (T e ) >36 000 K, indicating the transition to a nearly fully ionized regime. The ps–ns temporal dynamics are also compared between air and N 2 plasmas to assess the influence of collisional partners, as well as across the anode, cathode, and central gap regions to identify spatial variations in plasma behavior. Finally, this work demonstrates versatile continuously probing approach capable of spectro-temporal and spatially resolved single-shot measurements of key state variables in the ps–ns evolution of atmospheric-pressure nanosecond pulsed plasmas.

Electric discharges↗
DOE OSTI2026

Single-Molecule Electron Transport in Peptoids

Peptoids are structural analogs of peptides in which side chains are appended to the backbone nitrogen rather than the α-carbon. The sequence-defined modularity of peptoids enables precise control over structure−function relationships, enabling applications in energy storage and biomedical materials. Despite recent progress, the role of sequence and conformation on electron transport in peptoid molecules is not fully understood. Here, we synthesize a library of peptoid oligomers and characterize their molecular electronic properties using the scanning tunneling microscope-break junction (STM-BJ) technique. Our results show well-defined electron transport behavior for peptoid sequences containing aromatic side groups lacking hydrogen bonds (H-bonds) and without chemical substitutions at the N−C α position. This behavior fundamentally differs from electron transport in peptides, where H-bond interactions give rise to higher conductance states. All-atom molecular dynamics (MD) simulations are used to understand the conformational heterogeneity of peptoids, and molecular conformations obtained from MD simulations are used in quantum mechanical calculations based on the nonequilibrium Green’s function−density functional theory (NEGF-DFT) formalism. In all cases, computational results are in reasonable qualitative agreement with experiments. Our work demonstrates that the conductance behavior of peptoids depends on monomer identity, including side-chain aromaticity and substitution at the N−C α position. Overall, this work provides new insights into the structure−function relationships governing electron transport in peptoid-based materials and establishes design rules for peptoid-based molecular junctions.

Charge transport↗
NASA NTRS1993

Molecular gas dynamics applied to low-thrust propulsion

The Direct Simulation Monte Carlo method is currently being applied to study flowfields of small thrusters, including both the internal nozzle and the external plume flow. The DSMC method is employed because of its inherent ability to capture nonequilibrium effects and proper boundary physics in low-density flow that are not readily obtained by continuum methods. Accurate prediction of both the internal and external nozzle flow is important in determining plume expansion which, in turn, bears directly on impingement and contamination effects.

Zelesnik, Donna↗
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 OSTI2023

Cations Control Lipid Bilayer Memcapacitance Associated with Long-Term Potentiation

Phospholipid bilayers can be described as capacitors whose capacitance per unit area (specific capacitance, C m ) is determined by their thickness and dielectric constant–independent of applied voltage. It is also widely assumed that the C m of membranes can be treated as a “biological constant”. Recently, using droplet interface bilayers (DIBs), it was shown that zwitterionic phosphatidylcholine (PC) lipid bilayers can act as voltage-dependent, nonlinear memory capacitors, or memcapacitors. When exposed to an electrical “training” stimulation protocol, capacitive energy storage in lipid membranes was enhanced in the form of long-term potentiation (LTP), which enables biological learning and long-term memory. LTP was the result of membrane restructuring and the progressive asymmetric distribution of ions across the lipid bilayer during training, which is analogous, for example, to exponential capacitive energy harvesting from self-powered nanogenerators. Here, we describe how LTP could be produced from a membrane that is continuously pumped into a nonequilibrium steady state, altering its dielectric properties. During this time, the membrane undergoes static and dynamic changes that are fed back to the system’s potential energy, ultimately resulting in a membrane whose modified molecular structure supports long-term memory storage and LTP. Here, we also show that LTP is very sensitive to different salts (KCl, NaCl, LiCl, and TmCl 3 ), with LiCl and TmCl 3 having the most profound effect in depressing LTP, relative to KCl. This effect is related to how the different cations interact with the bilayer zwitterionic PC lipid headgroups primarily through electric-field-induced changes to the statistically averaged orientations of water dipoles at the bilayer headgroup interface.

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

Visualizing interfacial collective reaction behaviour of Li–S batteries

Benefiting from high energy density (2,600 Wh kg -1 )) and low cost, lithium-sulfur (Li-S) batteries are considered promising candidates for advanced energy-storage systems. Despite tremendous efforts in suppressing the long-standing shuttle effect of lithium polysulfides, understanding of the interfacial reactions of lithium polysulfides at the nanoscale remains elusive. This is mainly because of the limitations of in situ characterization tools in tracing the liquid-solid conversion of unstable lithium polysulfides at high temporal-spatial resolution. There is an urgent need to understand the coupled phenomena inside Li-S batteries, specifically, the dynamic distribution, aggregation, deposition and dissolution of lithium polysulfides. Here, by using in situ liquid-cell electrochemical transmission electron microscopy, we directly visualized the transformation of lithium polysulfides over electrode surfaces at the atomic scale. Notably, an unexpected gathering-induced collective charge transfer of lithium polysulfides was captured on the nanocluster active-centre-immobilized surface. It further induced an instantaneous deposition of nonequilibrium Li2S nanocrystals from the dense liquid phase of lithium polysulfides. Without mediation of active centres, the reactions followed a classical single-molecule pathway, lithium polysulfides transforming into Li 2 S 2 and Li 2 S step by step. Molecular dynamics simulations indicated that the long-range electrostatic interaction between active centres and lithium polysulfides promoted the formation of a dense phase consisting of Li + and S n 2- (2 < n ≤ 6), and the collective charge transfer in the dense phase was further verified by ab initio molecular dynamics simulations. Finally, the collective interfacial reaction pathway unveils a new transformation mechanism and deepens the fundamental understanding of Li-S batteries.

25 ENERGY STORAGE↗
NASA NTRS1997

Aerodynamics of Stardust Sample Return Capsule

Successful return of interstellar dust and cometary material by the Stardust Sample Return Capsule requires an accurate description of the Earth entry vehicle's aerodynamics. This description must span the hypersonic-rarefied, hypersonic-continuum, supersonic, transonic, and subsonic flow regimes. Data from numerous sources are compiled to accomplish this objective. These include Direct Simulation Monte Carlo analyses, thermochemical nonequilibrium computational fluid dynamics, transonic computational fluid dynamics, existing wind tunnel data, and new wind tunnel data. Four observations are highlighted: 1) a static instability is revealed in the free-molecular and early transitional-flow regime due to aft location of the vehicle s center-of-gravity, 2) the aerodynamics across the hypersonic regime are compared with the Newtonian flow approximation and a correlation between the accuracy of the Newtonian flow assumption and the sonic line position is noted, 3) the primary effect of shape change due to ablation is shown to be a reduction in drag, and 4) a subsonic dynamic instability is revealed which will necessitate either a change in the vehicle s center-of-gravity location or the use of a stabilizing drogue parachute.

Mitcheltree, R. A.↗
DOE OSTI2025

Nonequilibrium effects in high-gain inertial confinement fusion

Recent experimental demonstrations of ignition and target gain in inertial confinement fusion (ICF) have stimulated interest in exploring the fundamental physics of violent deuterium-tritium (DT) burn in high-gain ICF targets. A significant DT-burn fraction is a necessary condition for high energy gain and large neutron yields (>100MJ). Using classical molecular-dynamics (MD) simulations and a hybrid fluid-kinetic model, we examine how a large fraction of low-energy 𝛼 particles can kick D and T ions out of equilibrium in high-gain ICF targets. The MD results suggest that (1) temperatures of 𝑇 𝐷 and 𝑇 𝑇 can differ by as much as ∼20% of their mean temperature and (2) the deviation of the DT energy distribution from the Maxwell-Boltzmann function can exceed ∼30%. Some of these MD observations, such as the preferential heating of D ions by low-energy 𝛼 particles and the temperature separation, can be explained by a proposed hybrid fluid-kinetic model. Furthermore, the implication of such nonequilibrium effects on the DT reactivity is also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Steric effects of central dogma processes on the compaction and segregation of bacterial nucleoids

The bacterial cytoplasm is characterized by a distinctive membrane-less organelle, the nucleoid, which harbors the chromosomal DNA. Here, we investigate the steric effects of dynamic processes associated with transcription and translation on the structure of this organelle using coarse-grained molecular dynamics simulations that incorporate out-of-equilibrium reactions. Our model captures the scale of the entire cell and incorporates a reaction-diffusion system for ribosomes and polyribosomes, coupling their nonequilibrium kinetics to DNA excluded volume interactions. Our findings demonstrate that out-of-equilibrium reactions increase the size of the nucleoid and the number of ribosomes in this subcellular region. In addition, we show that nucleoid size is proportional to transcriptional activity. Our model reproduces the time-dependent change in nucleoid size observed in rifampicin treatment experiments, where the pool of polyribosomes is depleted. Furthermore, we find that these active processes are essential for complete sister chromosome separation and correct nucleoid positioning within the cell. Overall, our study reveals the effects of the central dogma processes on the internal organization and localization of bacterial nucleoids.

Chang, Mu-Hung [Univ. of Tennessee, Knoxville, TN ↗
DOE OSTI2021

Collective Vibrational Strong Coupling Effects on Molecular Vibrational Relaxation and Energy Transfer: Numerical Insights via Cavity Molecular Dynamics Simulations**

Abstract For a small fraction of hot CO 2 molecules immersed in a liquid‐phase CO 2 thermal bath, classical cavity molecular dynamics simulations show that forming collective vibrational strong coupling (VSC) between the C=O asymmetric stretch of CO 2 molecules and a cavity mode accelerates hot‐molecule relaxation. This acceleration stems from the fact that polaritons can be transiently excited during the nonequilibrium process, which facilitates intermolecular vibrational energy transfer. The VSC effects on these rates 1) resonantly depend on the cavity mode detuning, 2) cooperatively depend on Rabi splitting, and 3) collectively scale with the number of hot molecules. For larger cavity volumes, the average VSC effect per molecule can remain meaningful for up to N ≈10 4 molecules forming VSC. Moreover, the transiently excited lower polariton prefers to relax by transferring its energy to the tail of the molecular energy distribution rather than distributing it equally to all thermal molecules. As far as the parameter dependence is concerned, the vibrational relaxation data presented here appear analogous to VSC catalysis in Fabry–Pérot microcavities.

Li, Tao E.↗
DOE OSTI2021

Collective Vibrational Strong Coupling Effects on Molecular Vibrational Relaxation and Energy Transfer: Numerical Insights via Cavity Molecular Dynamics Simulations**

Abstract For a small fraction of hot CO 2 molecules immersed in a liquid‐phase CO 2 thermal bath, classical cavity molecular dynamics simulations show that forming collective vibrational strong coupling (VSC) between the C=O asymmetric stretch of CO 2 molecules and a cavity mode accelerates hot‐molecule relaxation. This acceleration stems from the fact that polaritons can be transiently excited during the nonequilibrium process, which facilitates intermolecular vibrational energy transfer. The VSC effects on these rates 1) resonantly depend on the cavity mode detuning, 2) cooperatively depend on Rabi splitting, and 3) collectively scale with the number of hot molecules. For larger cavity volumes, the average VSC effect per molecule can remain meaningful for up to N ≈10 4 molecules forming VSC. Moreover, the transiently excited lower polariton prefers to relax by transferring its energy to the tail of the molecular energy distribution rather than distributing it equally to all thermal molecules. As far as the parameter dependence is concerned, the vibrational relaxation data presented here appear analogous to VSC catalysis in Fabry–Pérot microcavities.

Li, Tao E.↗
DOE OSTI2025

Temperature Evolution of the Activation Barriers Leads to Meyer−Neldel Rules for Structural Relaxation and Transport in Polymers

Understanding activation barriers controlling structural relaxation in glass-forming liquids, molecular transport, and ionic conductivity in amorphous polymers is a grand challenge of fundamental scientific and materials engineering interest across disciplines. Over decades, intriguing but puzzling empirical correlations between the elementary time scale of activated barrier crossing and the apparent Arrhenius activation energy, the so-called Meyer−Neldel (MN) rules, have been discovered in diverse liquids and glasses. Here, in this study, we formulate and successfully apply a new experimental analysis and an explicitly dynamical theoretical framework which provides an understanding of the origin, validity, and failure of such correlations, that bridge and unify the three fields of structural relaxation, molecular transport, and ionic conductivity in liquids and quenched glasses. Distinct quasi-universal laws are predicted in equilibrated liquids and nonequilibrium glasses, consistent with experiments. Our analysis reveals that even if the relaxation appears Arrhenius over a limited temperature range, the physical activation barrier is generally temperature-dependent in polymeric systems even below glass transition temperature. In addition, we show that the approximate validity of classical MN rules hinges on a linear temperature dependence of this barrier and the temperature range probed in experiments. Our findings are relevant for controlling the activation barrier in functional soft polymeric materials relevant to molecular separations, barrier coatings, and charge transport, and also provide new constraints on the theoretical understanding of the mechanism underlying slow activated dynamics in glass-forming condensed matter.

Meyer-Neldel rules↗
DOE OSTI2025

Designing a quantum-accurate machine-learning potential to enable large-scale simulations of deuterium under shock

Large-scale molecular dynamics of deuterium under shock can elucidate kinetic processes vital to the target design in inertial confinement fusion and high-energy-density experiments. However, modeling the complex evolution of this material from an insulating molecular state at ambient pressure to an ionized, atomic fluid under strong shock is beyond the capability of simple pair and even bond order potentials. We thus train a quantum-accurate and broadly transferable machine-learning interatomic potential for deuterium using the Chebyshev Interaction Model for Efficient Simulations framework. We show that due to an improved description of the molecular-to-atomic transition, our model is able to better reproduce the ab initio equation of state, radial distribution functions, and principal Hugoniot than bond order potentials. This represents an important step toward large-scale quantum-accurate and nonequilibrium simulations of complicated systems under dynamic changes including phase transitions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗