Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “atomic 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 OSTI2026

Theoretical Insights into Reaction-Induced Transformation and Tuning of Catalytic Behavior in Heterogenous Catalysis

Reaction-induced transformations in heterogenous catalysis represent diverse phenomena that challenge traditional views of static catalyst surfaces. From surface adsorbate dynamics, atomic rearrangements, to composition and phase transitions, these processes reveal the profound differences between idealized model systems under ultrahigh vacuum and the complex, evolving interfaces that govern real catalytic behaviors under reaction conditions. Here, this perspective reviews recent theoretical efforts to provide atomic-level mechanistic insights into significant reaction-induced transformations and their impact on catalytic activity and selectivity. It underscores the need for an integrated framework that combines predictive simulations with operando characterization to uncover active sites and mechanisms under realistic operating conditions. Achieving this requires accelerating existing simulations to fully capture diverse reaction-induced surface dynamics, enabling scalable and accurate modeling of catalysts as condition-dependent, dynamically evolving systems. Such approaches are critical to bridge the gap between theory and practice, offering a pathway to more impactful and predictive catalyst design.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Coordination Sphere of Lanthanide Aqua Ions Resolved with Ab Initio Molecular Dynamics and X-ray Absorption Spectroscopy

To resolve the fleeting structures of lanthanide Ln 3+ aqua ions in solution, we (i) performed the first ab initio molecular dynamics (AIMD) simulations of the entire series of Ln 3+ aqua ions in explicit water solvent using pseudopotentials and basis sets recently optimized for lanthanides and (ii) measured the symmetry of the hydrating waters about Ln 3+ ions (Nd 3+ , Dy 3+ , Er 3+ , Lu 3+ ) for the first time with extended X-ray absorption fine structure (EXAFS). EXAFS spectra were measured experimentally and generated from AIMD trajectories to directly compare simulation, which concurrently considers the electronic structure and the atomic dynamics in solution, with experiment. In this study, we performed a comprehensive evaluation of EXAFS multiple-scattering analysis (up to 6.5 Å) to measure Ln–O distances and angular correlations (i.e., symmetry) and elucidate the molecular geometry of the first hydration shell. This evaluation, in combination with symmetry-dependent L 3 - and L 1 -edge spectral analysis, shows that the AIMD simulations remarkably reproduces the experimental EXAFS data. The error in the predicted Ln–O distances is less than 0.07 Å for the later lanthanides, while we observed excellent agreement with predicted distances within experimental uncertainty for the early lanthanides. Our analysis revealed a dynamic, symmetrically disordered first coordination shell, which does not conform to a single molecular geometry for most lanthanides. This work sheds critical light on the highly elusive coordination geometry of the Ln 3+ aqua ions.

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

Real-space local self-motion of protonated and deuterated water

Here, we report on the self-part of the Van Hove correlation function, the correlation function describing the dynamics of a single molecule, of water and deuterated water. The correlation function is determined by transforming inelastic scattering spectra of neutrons or x rays over a wide range of momentum transfer Q and energy transfer E to space R and time t. The short-range diffusivity is estimated from the Van Hove correlation function in the framework of the Gaussian approximation. The diffusivity has been found to be different from the long-range macroscopic diffusivity, providing information about local atomic dynamics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Heat transport in liquid water from first-principles and deep neural network simulations

In this work, we compute the thermal conductivity of water within linear response theory from equilibrium molecular dynamics simulations, by adopting two different approaches. In one, the potential energy surface (PES) is derived on the fly from the electronic ground state of density functional theory (DFT) and the corresponding analytical expression is used for the energy flux. In the other, the PES is represented by a deep neural network (DNN) trained on DFT data, whereby the PES has an explicit local decomposition and the energy flux takes a particularly simple expression. By virtue of a gauge invariance principle, established by Marcolongo, Umari, and Baroni, the two approaches should be equivalent if the PES were reproduced accurately by the DNN model. We test this hypothesis by calculating the thermal conductivity, at the GGA (PBE) level of theory, using the direct formulation and its DNN proxy, finding that both approaches yield the same conductivity, in excess of the experimental value by approximately 60%. Besides being numerically much more efficient than its direct DFT counterpart, the DNN scheme has the advantage of being easily applicable to more sophisticated DFT approximations, such as meta-GGA and hybrid functionals, for which it would be hard to derive analytically the expression of the energy flux. We find in this way that a DNN model, trained on meta-GGA (SCAN) data, reduces the deviation from experiment of the predicted thermal conductivity by about 50%, leaving the question open as to whether the residual error is due to deficiencies of the functional, to a neglect of nuclear quantum effects in the atomic dynamics, or, likely, to a combination of the two.

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

Fragile-to-strong transition in liquid As 2 S 3 under pressure: The effect of melt metallization

The well-known classification of glass-forming melts into fragile and strong liquids has several notable exceptions, including water, silica, and certain phase-change materials (PCMs). These exceptional fluid systems exhibit a fragile-to-strong transition (FST) upon cooling: a transformation from a high-temperature liquid with fast atomic dynamics, low viscosity, and low flow activation energy, to a viscous supercooled melt with high energy barriers near the glass transition temperature T g . This behavior is critically important for non-volatile memories, photonic tensor cores, reconfigurable metamaterials, and other devices, that use PCMs, enabling nanosecond-scale crystallization in the fragile regime and long data retention in the strong regime near or below T g . A significant structural transformation is expected between these two viscosity regimes, along with a semiconductor-metal (SC-M) transition upon heating, driven by high internal pressure and associated density increase. By applying high external pressure to the canonical low-conducting chalcogenide melt As 2 S 3 , we observed both the FST and the SC-M transition, occurring simultaneously within the same domain of the P, T−phase space. These findings suggest that the FST is not limited to a few exceptional liquids but is a common phenomenon, at least in systems that exhibit melt metallization within specific regions of their P, T−phase diagrams.

first-principles molecular dynamics↗
DOE OSTI2022

Large negative thermal expansion in GdFe(CN) 6 driven by unusual low-frequency modes

Understanding the negative thermal expansion (NTE) mechanism is of great importance. In this work, we consider the new NTE compound GdFe(CN) 6 (αv = -34.2×10 -6 K -1 ) as a case study to investigate the NTE mechanism from the perspective of the lattice vibrational dynamics. The atomic mean-square displacements suggest that the NTE of GdFe(CN) 6 comes from the strong tension effect induced by the transverse vibrations of the atomic –Fe–Ctriple bondN–Gd– linkages, with the largest contribution given by N atoms. Lattice dynamics calculations show that three low-frequency optical modes at about 50 cm -1 show the largest negative Grüneisen parameters thus providing the largest contribution to the NTE. In conclusion, the existence of these unusual low-frequency vibrational modes can be ascribed to the presence of GdN 6 trigonal prisms in the framework structure of GdFe(CN) 6 .

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

Transcritical diffuse-interface hydrodynamics of propellants in high-pressure combustors of chemical propulsion systems

Rocket engines and high-power new generations of gas-turbine jet engines and diesel engines oftentimes involve the injection of one or more reactants at subcritical temperatures into combustor environments at high pressures, and more particularly at pressures higher than those corresponding to the critical points of the individual components of the mixture, which typically range from 13 to 50 bars for most propellants. This class of trajectories in the thermodynamic space has been traditionally referred to as transcritical. However, the fundamental understanding of fuel atomization, vaporization, mixing, and combustion processes at such high pressures remains elusive. In particular, whereas fuel sprays are relatively well characterized at normal pressures, analyses of dispersion of fuel in high-pressure combustors are hindered by the limited experimental diagnostics and theoretical formulations available. The description of the thermodynamics of hydrocarbon-fueled mixtures employed in chemical propulsion systems is complex and involves mixing-induced phenomena, including an elevation of the critical point whereby the coexistence region of the mixture extends up to pressures much larger than the critical pressures of the individual components. As a result, interfaces subject to surface-tension forces may persist in multicomponent systems despite the high pressures, and may give rise to unexpected spray-like atomization dynamics that are otherwise absent in monocomponent systems above their critical point. In this article, the current understanding of this phenomenon is reviewed within the context of propulsion systems fueled by heavy hydrocarbons. Emphasis is made on analytical descriptions at mesoscopic scales of interest for computational fluid dynamics. In particular, a set of modifications of the constitutive laws in the Navier–Stokes equations for multicomponent flows, supplemented with a high-pressure equation of state and appropriate redefinitions of the thermodynamic potentials, are introduced in this work based on an extended version of the diffuse-interface theory of van der Waals. The resulting formulation involves revisited forms of the stress tensor and transport fluxes of heat and species, and enables a description of the mesoscopic volumetric effects induced by transcritical interfaces consistently with the thermodynamic phase diagram of the mixture at high pressures. Additionally, applications of the theory are illustrated in canonical problems, including dodecane/nitrogen transcritical interfaces in non-isothermal systems. The results indicate that a transcritical interface is formed between the propellant streams that persists downstream of the injection orifice over distances of the same order as the characteristic thermal-entrance length of the fuel stream. The transcritical interface vanishes at an edge that gives rise to a fully supercritical mixing layer.

02 PETROLEUM↗
DOE OSTI2021

Investigation of post-injection strategies for diesel engine Catalyst Heating Operation using a vapor-liquid-equilibrium-based spray model

Most multidimensional engine simulations spend much time solving for non-equilibrium spray dynamics (atomization, collision, vaporization). However, their accuracy is limited by significant grid dependency, and the need for extensive calibration. This is critical for modeling cold-start diesel fuel post injections, which occur at low temperatures and pressures, far from typical model validation ranges. At the same time, resolving micron-scale spray phenomena would render full Eulerian multiphase calculations prohibitive. In this study, an improved phase equilibrium based approach was implemented and assessed for simulating diesel catalyst heating operation strategies. A phase equilibrium solver based on the model by Yue and Reitz [1] was implemented: a fully multiphase CFD solver is employed with an engineering-size engine grid, and fuel injection is modeled using the standard Lagrangian parcels approach. Mass and energy from the liquid parcels are released to the Eulerian multiphase mixture according to an equilibrium-based liquid jet model. An improved phase equilibrium solver was developed to handle large real-gas mixtures such as those from accurate chemical kinetics mechanisms. The liquid-jet model was improved such that momentum transfer to the Eulerian solver better reproduces the physical spray jet structure. Validation of liquid/vapor penetration predictions showed that the model yields accurate results with very limited tuning and low sensitivity to the few calibration constants. In-cylinder simulations of diesel catalyst heating operation strategies showed that capturing spray structure is paramount when short, transient injection pulses and low temperatures are present. Furthermore, the EP model provides improved predictions of post-injection spray structure and ignitability, while conventional spray modeling does not capture the increase of liquid penetration during the expansion stroke. Finally, the only important EP model calibration constant, Cliq, does not affect momentum transfer, but it changes the local charge cooling distribution through the local energy transfer, which makes it candidate to additional research. The results confirm that non-equilibrium spray processes do not need to be resolved in engineering simulations of high-pressure diesel sprays.

42 ENGINEERING↗
DOE OSTI2020

Negatively Charged Lipids Exhibit Negligible Effects on the Water Repellency of Montmorillonite Films

Amphiphilic molecules can alter the wettability of soil minerals. To determine how the headgroup chemistry of amphiphiles determines these effects, we investigate a system of the clay montmorillonite with long-chain phospholipids. We use phosphatidylglycerol (PG) phospholipids to contrast with our previous work using phosphatidylethanolamine lipids (PE). Zwitterionic PE lipids can sorb to the negatively-charged montmorillonite surface, whereas negatively-charged PG lipids cannot. Employing a suite of techniques from molecular dynamics, atomic force microscopy, fluorescence microscopy, and contact angle measurements, we define sample characteristics from molecular scale structure to the macroscopic wettability. We find that PG lipids do not significantly alter montmorillonite wetting characteristics, such as the contact angle, flow viscosity, and the characteristic timescale for droplet imbibition. In comparing PE and PG lipid/clay films, we find that, among the phospholipids compared, they must have three characteristics in order to change clay/lipid film wettability: they must bind to the mineral surface, be solid at room temperature, and have a relatively continuous distribution throughout the film.

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

The critical role of hot carrier cooling in optically excited structural transitions

The hot carrier cooling occurs in most photoexcitation-induced phase transitions (PIPTs), but its role has often been neglected in many theoretical simulations as well as in proposed mechanisms. Here, by including the previously ignored hot carrier cooling in real-time time-dependent density functional theory (rt-TDDFT) simulations, we investigated the role of hot carrier cooling in PIPTs. Taking IrTe2 as an example, we reveal that the cooling of hot electrons from the higher energy levels of spatially extended states to the lower energy levels of the localized Ir–Ir dimer antibonding states strengthens remarkably the atomic driving forces and enhances atomic kinetic energy. These two factors combine to dissolute the Ir–Ir dimers on a timescale near the limit of atomic motions, thus initiating a deterministic kinetic phase transition. We further demonstrate that the subsequent cooling induces nonradiative recombination of photoexcited electrons and holes, leading to the ultrafast recovery of the Ir–Ir dimers observed experimentally. These findings provide a complete picture of the atomic dynamics in optically excited structural phase transitions.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Temperature Control of Crystal Optics for Ultrahigh-Resolution Applications

The temperature control of crystal optics is critical for ultrahigh resolution applications such as those used in meV-resolved Inelastic Scattering. Due to the low count rate and long acquisition time of these experiments, for 1-meV energy resolution, the absolute temperature stability of the crystal optics must be maintained below 4 mK to ensure the required stability of lattice constant, thereby ensuring the energy stability of the optics. Furthermore, the temperature control with sub-mK precision enables setting the absolute temperature of individual crystal, making it possible to align the reflection energy of each crystal’s rocking curve in sub-meV resolution thereby maximizing the combined efficiency of the crystal optics. In this contribution, we report the details of an EPICS control system using PT1000 sensors, Keithley 3706A 7.5 digits sensor scanner, and Wiener MPOD LV power supply for the analyzer crystals of the Inelastic X-ray Scattering (IXS) beamline 10-ID at NSLS-II**. We were able to achieve absolute temperature stability below 1 mK and sub-meV energy alignment for several asymmetrically cut analyzer crystals. The EPICS ePID record was used for the control of the power supplies based on the PT1000 sensor input that was read with 7.5 digits accuracy from the Keithley 3706A scanner. The system enhances the performance of the meV-resolved IXS spectrometer with currently a 1.4 meV total energy resolution and unprecedented spectral sharpness for studies of atomic dynamics in a broad range of materials.

Gofron, Kaz↗
DOE OSTI2025

Time-resolved atomic-resolution Brownian tomography of single nanocrystals reveals size-dependent dynamics

Atomic-resolution structure identification of nanocrystals by graphene liquid cell electron microscopy (GLC-EM) has revealed that small, solubilized platinum nanocrystals consist of an ordered crystalline core surrounded by mobile surface atoms, which dissociate during oxidative etching, resulting in distinct temporal structural states. Requirements imposed by the 3D reconstruction algorithm limit the number of structural states that can be resolved. We introduce a regularized 3D reconstruction algorithm that exploits the redundancy inherent in the experimental data, allowing us to improve the time resolution. Our developments provide a comprehensive molecular picture at unprecedented spatial and temporal resolution of the nonlinear, linear, and fluctuating dynamic phenomena that single nanocrystals undergo during the GLC-EM experiment. We determined atomic structures of 66 temporal structural states, extracted from 15 time trajectories of individual nanocrystals. Large (478 to 698 atoms) and small (<300 atoms) nanocrystals show etching that preserves a stable core, whereas mid-sized (351 to 571 atoms) nanocrystals present dynamics that change the coordination of the core.

Meana-Pañeda, Rubén↗
DOE OSTI2024

Compiling Quantum Circuits for Dynamically Field-Programmable Neutral Atoms Array Processors

Dynamically field-programmable qubit arrays (DPQA) have recently emerged as a promising platform for quantum information processing. In DPQA, atomic qubits are selectively loaded into arrays of optical traps that can be reconfigured during the computation itself. Leveraging qubit transport and parallel, entangling quantum operations, different pairs of qubits, even those initially far away, can be entangled at different stages of the quantum program execution. Such reconfigurability and non-local connectivity present new challenges for compilation, especially in the layout synthesis step which places and routes the qubits and schedules the gates. In this paper, we consider a DPQA architecture that contains multiple arrays and supports 2D array movements, representing cutting-edge experimental platforms. Within this architecture, we discretize the state space and formulate layout synthesis as a satisfiability modulo theories problem, which can be solved by existing solvers optimally in terms of circuit depth. For a set of benchmark circuits generated by random graphs with complex connectivities, our compiler OLSQ-DPQA reduces the number of two-qubit entangling gates on small problem instances by 1.7x compared to optimal compilation results on a fixed planar architecture. To further improve scalability and practicality of the method, we introduce a greedy heuristic inspired by the iterative peeling approach in classical integrated circuit routing. Using a hybrid approach that combined the greedy and optimal methods, we demonstrate that our DPQA-based compiled circuits feature reduced scaling overhead compared to a grid fixed architecture, resulting in 5.1X less two-qubit gates for 90 qubit quantum circuits. These methods enable programmable, complex quantum circuits with neutral atom quantum computers, as well as informing both future compilers and future hardware choices.

Physics↗
DOE OSTI2022

Quantum dynamics of cold atomic gas with SU(1,1) symmetry

Motivated by recent advances in quantum dynamics, we investigate the dynamics of the system with SU(1,1) symmetry. Instead of performing the time-ordered integral for the evolution operator of the time-dependent Hamiltonian, we show that the time evolution operator can be expressed as an SU(1,1) group element. Since the SU(1,1) group describes the “rotation” on a hyperbolic surface, the dynamics can be visualized on a Poincaré disk, a stereographic projection of the upper hyperboloid. As an example, we present the trajectory of the revival of Bose-Einstein condensation and that of the scale-invariant Fermi gas on the Poincaré disk. Further considering quantum gases in an oscillating lattice, we also study the dynamics of the system with time-dependent single-particle dispersion.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS1987

Dynamics of atom-surface interactions

Scattering techniques currently being used to examine the dynamics and energetics of gas-surface energy exchange processes and gas-surface reaction mechanisms are reviewed. In particular, recent inelastic scattering measurements are highlighted which are revealing the microscopic basis for collision-induced gas-surface energy exchange, e.g., which surface vibrational modes actively participate in translational energy accommodation. Reactive scattering and laser desorptive experiments which examine energy disposal in volatile products are also discussed. Finally, an efficient atomic oxygen beam source is described which is suitable for terrestrial studies of gas-surface interactions.

Sibener, Steven J.↗
DOE OSTI2020

Highly selective electrocatalytic CO 2 reduction to ethanol by metallic clusters dynamically formed from atomically dispersed copper

Direct electrochemical conversion of CO 2 to ethanol offers a promising strategy of lowering CO 2 emission while storing energy from renewable electricity. However, current electrocatalysts offer only limited selectivity toward ethanol. Here we report a copper catalyst synthesized by a unique Cu-Li amalgm method over a commercial carbon support that achieved Faradaic efficiency (FE) higher than 91% at -0.7 V (RHE) and the active potential as low as -0.4 V (RHE) during direct electrocatalytic CO 2 -to-ethanol conversion. The catalyst also demonstrated stability over an extended period of operation. A strong correlation between the catalytic selectivity and the initial Cu atoms dispersion was found and Operando X-ray absorption spectroscopy identified a dynamic and reversible transformation from atomically dispersed copper atoms to Cu n (n = 3 and 4) under the electrochemical reaction. Frist-principles calculations further elucidate the possible catalytic mechanism of CO 2 reduction over Cu n .

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

False vacuum decay and nucleation dynamics in neutral atom systems

Metastable states of quantum many-body systems with confinement offer a means to simulate false vacuum phenomenology, including nonequilibrium dynamical processes like decay by nucleation, in truncated limits. Recent work has examined the decay process in one-dimensional (1D) ferromagnetic Ising spins and superfluids. Here, in this paper, we study nucleation dynamics in 1D antiferromagnetic neutral atom chains with Rydberg interactions, using both numerical simulations and analytic modeling. We apply a staggered local detuning field to generate the metastable and ground states. Our efforts focus on two dynamical regimes: decay and annealing. In the first, we corroborate the phenomenological decay rate scaling and determine the associated parameter range for the decay process; in the second, we uncover and elucidate a procedure to anneal the metastable state from the initial to the final system, with intermediate nucleation events. We further propose experimental protocols to prepare the required states and perform quenches on near-term neutral atom quantum simulators, examining the experimental feasibility of our proposed setup and parameter regime.

97 MATHEMATICS AND COMPUTING↗