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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Atomic Mixing of Initially Separated Mixtures

We consider the atomic mixing of two materials, initially separated in space, at fixed temperature– pressure conditions [T, p], and zero mass-mean velocity. Each material is composed of a specified set of atomic species whose density corresponds to the specieswise real Equation of State (EOS) $\tilde{&#961}$ i [T, p]. Because there are no spatial gradients in the mean [T, p] state, the only mechanism for mixing is a specieswise nonequilibrium thermal velocity distribution. The recipe for nonequilibrium species mass fluxes is taken directly from Hirschfelder, Curtiss, and Bird, who express the diffusion velocities in terms of gradients in the molar concentrations and pairwise coefficients of diffusion. EOS data from the SESAME library is used for all species. Data for the specieswise self diffusion coefficients is adopted from the Orbital Free Molecular Dynamics (OFMD) simulations of Ticknor et al.; the so-called Darken approximation is used to express the mutual diffusion coefficients of a mixture in the manner used by White et al. A one-dimensional spatial domain is used; with the left half occupied by one multispecies material and the right half occupied by the other, and with zero–gradient boundary conditions at the ends of the domain. Hence the total mass in the domain is a constant. At time t = 0 the materials are separated. Accordingly, spatial gradients at the separation point are infinite and the species fluxes are also infinite. As t → $\textbf{∞}$ the spatial gradients approach zero and a uniform set of species concentrations is obtained; at this point the system is said to be at “thermodynamic equilibrium.” Analytic solution of the problem is not possible in general due to the nonlinear concentration dependence of the species fluxes. Hence we utilize the numerical solution apparatus in Mathematica for obtaining accurate time-dependent solutions for the mixing of various initial species masses for conditions of interest in laser-driven Inertial Confinement Fusion. Concentration profiles and the speed of propagation of the mixing fronts are discussed in light of the nonlinearity of the problem.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Three-state harmonic models for photoinduced charge transfer

A widely used strategy for simulating the charge transfer between donor and acceptor electronic states in an all-atom anharmonic condensedphase system is based on invoking linear response theory to describe the system in terms of an effective spin-boson model Hamiltonian. Extending this strategy to photoinduced charge transfer processes requires also taking into consideration the ground electronic state in addition to the excited donor and acceptor electronic states. In this paper, we revisit the problem of describing such nonequilibrium processes in terms of an effective three-state harmonic model. We do so within the framework of nonequilibrium Fermi’s golden rule (NE-FGR) in the context of photoinduced charge transfer in the carotenoid–porphyrin–C 60 (CPC 60 ) molecular triad dissolved in explicit tetrahydrofuran (THF). To this end, we consider different ways for obtaining a three-state harmonic model from the equilibrium autocorrelation functions of the donor–acceptor, donor–ground, and acceptor–ground energy gaps, as obtained from all-atom molecular dynamics simulations of the CPC 60 /THF system. The quantum-mechanically exact time-dependent NE-FGR rate coefficients for two different charge transfer processes in two different triad conformations are then calculated using the effective three-state model Hamiltonians as well as a hierarchy of more approximate expressions that lead to the instantaneous Marcus theory limit. Our results show that the photoinduced charge transfer in CPC 60 /THF can be described accurately by the effective harmonic three-state models and that nuclear quantum effects are small in this system.

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

Kibble–Zurek scaling of nonequilibrium phase transition in barium titanate

Far-from-equilibrium phase transition dynamics is one of the grand challenges in modern materials science. A theoretical landmark is the Kibble–Zurek (KZ) scaling to describe the relationship between the temperature quenching rate and the resulting defect density in the vicinity of symmetry-breaking phase transformations. Despite the confirmation of the KZ scaling in ferroic perovskite materials and macroscopic simulations, its atomistic mechanisms remain elusive. Here, we demonstrate the KZ scaling using all-atom molecular dynamics simulations for a prototypical ferroelectric perovskite, barium titanate, with the scaling exponent corresponding to the theoretical prediction for rapid quenching. Simulated diffuse neutron scattering data are presented to guide future experiments.

Physics↗
DOE OSTI2023

Athermal annealing of pre-existing defects in crystalline silicon

Systematic investigations of electronic energy loss (S e ) effects on pre-existing defects in crystalline silicon (Si) are crucial to provide reliance on the use of ionizing irradiation to anneal pre-existing defects, leading to successful implementation of this technology in the fabrication of Si-based devices. In this regard, the S e effects on nonequilibrium defect evolution in pre-damaged Si single crystals at 300 K has been investigated using intermediate-energy ions (12 MeV O and Si ions) that interact with the pre-damaged surface layers of Si mainly by ionization, except at the end of their range where the nuclear energy loss (S n ) is no longer negligible. Furthermore, under these irradiation conditions, experimental results and molecular dynamics simulations have revealed that pre-existing disorder in Si can be almost fully annealed by subsequent irradiation with intermediate-energy incident ions with S e values as low as 1.5–3.0 keV/nm. Selective annealing of pre-existing defect levels in Si at room temperature can be considered as an effective strategy to mediate the transient enhanced diffusion of dopants in Si. This approach is more desirable than the regular thermal annealing, which is not compatible with the processing requirements that fall below the typical thermal budget.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Amphiphilic Baskets for Supramolecular Nanoarchitectures at Interfaces: Inverted Monolayer Formation on Water

Interfacial chemistry of molecular baskets remains poorly understood despite their promise for supramolecular applications of detection and sequestration of toxic molecules including those of illicit drugs, organophosphorus compounds, and anticancer agents. We present a fundamental investigation of the interfacial behavior of three amphiphilic supramolecular baskets (ASB 4, 8, and 12), having increasingly longer yet linear alkyl chains at the top of their bowl-shaped cavity. The studies were completed at the air−water interface to elucidate surface activity, interfacial stability, self-assembly, and monolayer organization that drive inverted monolayer formation, in which the molecular arms orient toward the aqueous phase in a configuration opposite to that typically observed for lipids. Herein, surface pressure−area isotherms of ASB 4, 8, 12, deposited on a water surface, were performed in tandem with nonequilibrium relaxation experiments to quantify surface activity, thermodynamic stability, and monolayer compressibility of the baskets’ monolayer assembly. Brewster angle microscopy enabled direct visualization of morphological evolution, aggregation, and packing at the interface. We show that systematic extension of the hydrocarbon arms, from four to 12 methylene groups, progressively modifies intermolecular packing, drives distinct two-dimensional aggregation pathways, and increases number densities at the air−water interface. Atomistic molecular dynamics simulations corroborate many of these experimentally observed trends and provide mechanistic detail on the cooperative roles of basket topology and interfacial concentration in regulating the hydration structure and dynamics within the cavities generated by surface-adsorbing baskets, consistent with observed variations in surface activity and packing. Our results establish how the topology of these unique supramolecules and their concentration govern interfacial organization and offer a rational framework for designing amphiphiles with predictable behavior at soft interfaces.

Hydrocarbons↗
DOE OSTI2021

Metastable–solid phase diagrams derived from polymorphic solidification kinetics

Nonequilibrium processes during solidification can lead to kinetic stabilization of metastable crystal phases. A general framework for predicting the solidification conditions that lead to metastable-phase growth is developed and applied to a model face-centered cubic (fcc) metal that undergoes phase transitions to the body-centered cubic (bcc) as well as the hexagonal close-packed phases at high temperatures and pressures. Large-scale molecular dynamics simulations of ultrarapid freezing show that bcc nucleates and grows well outside of the region of its thermodynamic stability. An extensive study of crystal–liquid equilibria confirms that at any given pressure, there is a multitude of metastable solid phases that can coexist with the liquid phase. We define for every crystal phase, a solid cluster in liquid (SCL) basin, which contains all solid clusters of that phase coexisting with the liquid. Additionally, a rigorous methodology is developed that allows for practical calculations of nucleation rates into arbitrary SCL basins from the undercooled melt. It is demonstrated that at large undercoolings, phase selections made during the nucleation stage can be undone by kinetic instabilities amid the growth stage. On these bases, a solidification–kinetic phase diagram is drawn for the model fcc system that delimits the conditions for macroscopic grains of metastable bcc phase to grow from the melt. We conclude with a study of unconventional interfacial kinetics at special interfaces, which can bring about heterogeneous multiphase crystal growth. A first-order interfacial phase transformation accompanied by a growth-mode transition is examined.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Electron-Induced Proton Transfer

The pathway of activationless proton transfer induced by an electron-transfer reaction is studied theoretically. Long-range electron transfer produces highly nonequilibrium medium polarization that can drive proton transfer through an activationless transition during the process of thermalization, dynamically altering the screening of the electron–proton Coulomb interaction by the medium. The cross electron–proton reorganization energy is the main energy parameter of the theory, which exceeds in magnitude the proton-transfer reorganization energy roughly by the ratio of the electron-transfer to proton-transfer distance. This parameter, which can be either positive or negative, is related to the difference in pK a values in two electron-transfer states. The relaxation time of the medium is on the (sub)picosecond time scale, which establishes the characteristic time for activationless proton transfer. Microscopic calculations predict substantial retardation of the collective relaxation dynamics compared to the continuum estimates due to the phenomenology analogous to de Gennes narrowing. As a result, nonequilibrium medium configuration promoting proton transfer can be induced by either thermal or photoinduced charge transfer.

14 SOLAR ENERGY↗
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 OSTI2024

Quantum Simulations of Radiation Damage in a Molecular Polyethylene Analog

Abstract An atomic‐level understanding of radiation‐induced damage in simple polymers like polyethylene is essential for determining how these chemical changes can alter the physical and mechanical properties of important technological materials such as plastics. Ensembles of quantum simulations of radiation damage in a polyethylene analog are performed using the Density Functional Tight Binding method to help bind its radiolysis and subsequent degradation as a function of radiation dose. Chemical degradation products are categorized with a graph theory approach, and occurrence rates of unsaturated carbon bond formation, crosslinking, cycle formation, chain scission reactions, and out‐gassing products are computed. Statistical correlations between product pairs show significant correlations between chain scission reactions, unsaturated carbon bond formation, and out‐gassing products, though these correlations decrease with increasing atom recoil energy. The results present relatively simple chemical descriptors as possible indications of network rearrangements in the middle range of excitation energies. Ultimately, the work provides a computational framework for determining the coupling between nonequilibrium chemistry in polymers and potential changes to macro‐scale properties that can aid in the interpretation of future radiation damage experiments on plastic materials.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2023

An ab initio free energy study of the reaction mechanism and rate-limiting steps of CO 2 capture by aqueous glycine

Aqueous amino acids are promising absorbents for direct air capture (DAC) of CO 2 . Herein, we investigate the possibility of kinetic control of CO 2 absorption with aqueous anionic glycine (GLY - ) by employing extensive ab initio molecular dynamics simulations, free energy analysis, and reaction rate theory. We find that first GLY - binds to CO 2 by overcoming a barrier (7.4 kcal/mol) to form a zwitterion intermediate, which then releases a proton by overcoming a similar barrier. Despite the similarity in the barrier, zwitterion formation appears to be the rate-limiting step because it is two orders of magnitude slower (microseconds) than the proton release step. This is predominantly due to stronger nonequilibrium solvent effects for the former that cause many barrier-recrossing events and effectively slow down the reaction rate. Such a detailed fundamental understanding of the amino acid-based CO 2 -absorption mechanism and rates is key to improving the kinetic efficiency of DAC technology.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1997

Data Base for CFD Validation

The flow behind the shock wave formed around objects which fly at hypervelocity behaves differently from that of a perfect gas. Molecules become vibrationally excited, dissociated, and ionized. The hot gas may emit or absorb radiation. When the atoms produced by dissociation reach the wall surface, chemical reactions, including recombination, may occur. The thermochemical phenomena of vibration, dissociation, ionization, surface chemical reaction, and radiation are referred to commonly as high-temperature real-gas phenomena. The phenomena cause changes in the dynamic behavior of the flow and the surface pressure and heat transfer distribution around the object. The character of a real gas is described by the internal degrees of freedom and state of constituent molecules; nitrogen and oxygen for air. The internal energy states, rotation, vibration and electronic, of the molecules are excited and, in the limit, the molecular bonds are exceeded and the gas dissociated into atomic and, possibly, ionic constituents. The process of energy transfer causing excitation, dissociation and recombination is a rate process controlled by particle collisions. Binary, two-body, collisions are sufficient to cause internal excitation, dissociation and ionization while three-body collisions are required to recombine the particles into molecular constituents. If the rates of energy transfer are fast with respect to the local fluid dynamic time scale the gas is in, or nearly in, equilibrium. If the energy transfer rates are very slow the gas can be described as frozen. In all other instances, wherein any of the energy exchange rates are comparable to the local fluid time scale, the gas will be thermally or chemically reacting and out of equilibrium. Real gas thermochemical nonequilibrium processes are important in the determination of aerodynamic heating; both convective (including wall catalytic effects) and radiative heating. To illustrate this we consider the hypervelocity flow over a bluff body typical of an atmospheric entry vehicle or an aerospace transfer vehicle (ASTV). The qualitative aspects of a hypersonic flow field over a bluff body are discussed in two parts, forebody and afterbody, with attention to which particular physical effects must be included in an analysis. This will indicate what type of numerical modeling will be adequate in each region of the flow. A bluff forebody flow field is dominated by the presence of the strong bow shock wave and the consequent heating, and chemical reaction of the gas. At high altitude hypersonic flight conditions the thermal excitation and chemical reaction of the gas occur slowly enough that a significant portion of the flow field is in a state of thermochemical nonequilibrium. A second important effect is the presence of the thick boundary layer along the forebody surface. In this region there are large thermal and chemical species gradients due to the interaction of the gas with the wall. Also at high altitudes the shock wave and the boundary layer may become so thick that they merge; in this case the entire shock layer is dominated by viscous effects.

Deiwert, George S.↗
DOE OSTI2022

Fabricating Graphene Oxide/h-BN Metal Insulator Semiconductor Diodes by Nanosecond Laser Irradiation

To employ graphene’s rapid conduction in 2D devices, a heterostructure with a broad bandgap dielectric that is free of traps is required. Within this paradigm, h-BN is a good candidate because of its graphene-like structure and ultrawide bandgap. We show how to make such a heterostructure by irradiating alternating layers of a-C and a-BN film with a nanosecond excimer laser, melting and zone-refining constituent layers in the process. With Raman spectroscopy and ToF-SIMS analyses, we demonstrate this localized zone-refining into phase-pure h-BN and rGO films with distinct Raman vibrational modes and SIMS profile flattening after laser irradiation. Furthermore, in comparing laser-irradiated rGO-Si MS and rGO/h-BN/Si MIS diodes, the MIS diodes exhibit an increased turn-on voltage (4.4 V) and low leakage current. The MIS diode I-V characteristics reveal direct tunneling conduction under low bias and Fowler-Nordheim tunneling in the high-voltage regime, turning the MIS diode ON with improved rectification and current flow. This study sheds light on the nonequilibrium approaches to engineering h-BN and graphene heterostructures for ultrathin field effect transistor device development.

2D materials↗
DOE OSTI2020

Addition of Short Polymer Chains Mechanically Reinforces Glassy Poly(2-vinylpyridine)–Silica Nanoparticle Nanocomposites

The addition of hard fillers to a polymer matrix is a well-known process for achieving mechanical reinforcement. With a decrease in the size of the fillers, the contribution from polymer–particle nanometer-sized interfaces become significant, and these interfaces affect the mechanical performance of polymer nanocomposites (PNCs) beyond the limits established for conventional composites. However, the molecular mechanisms underlying the improvement in the mechanical performance of glassy PNCs remain unresolved, necessitating a deeper understanding of the structure–property relationships in these intrinsically heterogeneous systems. In this effort, by using Brillouin light scattering (BLS) and dynamic mechanical analysis (DMA), we demonstrated that adding shorter chains to a PNC prepared with high molecular weight polymers significantly improved the mechanical properties of the PNC in the glassy state. The strongest enhancement of mechanical properties occurred at an optimum concentration of short chains. This is in contrast to behavior of the glass transition temperature of PNCs which shows a monotonic decrease with an increase in the concentration of shorter chains. Using experimental data and coarse-grained molecular dynamics (MD) simulations, we have identified the molecular mechanism leading to the observed non-monotonic changes in mechanical reinforcement. Here, this mechanism includes changes in the nanoscale organization at the interface, combined with chain stretching amplified by the addition of the short chains. Overall, our approach paves a simple, cost-effective pathway to fabricating glassy PNCs with significantly improved mechanical properties that will fill various practical needs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1992

Radiative energy transfer in molecular gases

Basic formulations, analyses, and numerical procedures are presented to study radiative interactions in gray as well as nongray gases under different physical and flow conditions. After preliminary fluid-dynamical considerations, essential governing equations for radiative transport are presented that are applicable under local and nonlocal thermodynamic equilibrium conditions. Auxiliary relations for relaxation times and spectral absorption models are also provided. For specific applications, several simple gaseous systems are analyzed. The first system considered consists of a gas bounded by two parallel plates having the same temperature. Within the gas there is a uniform heat source per unit volume. For this system, both vibrational nonequilibrium effects and radiation conduction interactions are studied. The second system consists of fully developed laminar flow and heat transfer in a parallel plate duct under the boundary condition of a uniform surface heat flux. For this system, effects of gray surface emittance are studied. With the single exception of a circular geometry, the third system is considered identical to the second system. Here, the influence of nongray walls is also studied.

Tiwari, Surendra N.↗
DOE OSTI2021

Laser-Induced Cooperative Transition in Molecular Electronic Crystal

The competing and non-equilibrium phase transitions, involving dynamic tunability of cooperative electronic and magnetic states in strongly correlated materials, show great promise in quantum sensing and information technology. To date, the stabilization of transient states is still in the preliminary stage, particularly with respect to molecular electronic solids. Here, a dynamic and cooperative phase in potassium-7,7,8,8-tetracyanoquinodimethane (K-TCNQ) with the control of pulsed electromagnetic excitation is demonstrated. Simultaneous dynamic and coherent lattice perturbation with 8 ns pulsed laser (532 nm, 15 MW cm -2 , 10 Hz) in such a molecular electronic crystal initiates a stable long-lived (over 400 days) conducting paramagnetic state (≈42 Ωcm), showing the charge–spin bistability over a broad temperature range from 2 to 360 K. Comprehensive noise spectroscopy, in situ high-pressure measurements, electron spin resonance (ESR), theoretical model, and scanning tunneling microscopy/spectroscopy (STM/STS) studies provide further evidence that such a transition is cooperative, requiring a dedicated charge–spin–lattice decoupling to activate and subsequently stabilize nonequilibrium phase. Finally, the cooperativity triggered by ultrahigh-strain-rate (above 10 6 s -1 ) pulsed excitation offers a collective control toward the generation and stabilization of strongly correlated electronic and magnetic orders in molecular electronic solids and offers unique electro-magnetic phases with technological promises.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Quantum Computation of Dynamical Quantum Phase Transitions and Entanglement Tomography in a Lattice Gauge Theory

Strongly coupled gauge theories far from equilibrium may exhibit unique features that could illuminate the physics of the early universe and of hadron and ion colliders. Studying real-time phenomena has proven challenging with classical-simulation methods but is a natural application of quantum simulation. To demonstrate this prospect, we quantum compute nonequal-time correlation functions and perform entanglement tomography of nonequilibrium states of a simple lattice gauge theory, the Schwinger model, using a trapped-ion quantum computer by IonQ Inc. As an ideal target for near-term devices, a recently predicted dynamical quantum phase transition in this model is studied by preparing, quenching, and tracking the subsequent nonequilibrium dynamics in three ways: (i) overlap echos signaling dynamical transitions, (ii) nonequal-time correlation functions with an underlying topological nature, and (iii) the entanglement structure of nonequilibrium states, including entanglement Hamiltonians. These results constitute the first observation of a dynamical quantum phase transition in a lattice gauge theory on a quantum computer and are a first step toward investigating topological phenomena in nuclear and high-energy physics using quantum technologies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Ab initio electron dynamics in high electric fields: Accurate prediction of velocity-field curves

Electron dynamics in external electric fields governs the behavior of solid-state electronic devices. First-principles calculations enable precise predictions of charge transport in low electric fields. However, studies of high-field electron dynamics remain elusive due to a lack of accurate and broadly applicable methods. In this work, we develop an efficient approach to solve the real-time Boltzmann transport equation with both the electric field term and ab initio electron-phonon collisions. These simulations provide field-dependent electronic distributions in the time domain, allowing us to investigate both transient and steady-state transport in electric fields ranging from low to high (> 10 kV/cm). The broad capabilities of our approach are shown by computing nonequilibrium electron occupations and velocity-field curves in Si, GaAs, and graphene, obtaining results in quantitative agreement with experiment. Our approach sheds light on microscopic details of transport in high electric fields, including the dominant scattering mechanisms and valley occupation dynamics. Our results demonstrate quantitatively accurate calculations of electron dynamics in low to high electric fields, with broad application to power and micro-electronics, optoelectronics, and sensing.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗