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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Probing High-Pressure Structural Evolution in Polyurea with In Situ Energy-Dispersive X-ray Diffraction and Molecular Dynamics Simulations

Polyurea, an elastomer with a phase-segregated microstructure, has been proven as an effective coating in defense applications. To gain a more complete understanding of the high-pressure atomic-level morphology of these phases and to validate molecular dynamics (MD) simulations, multi-angle energy-dispersive X-ray diffraction experiments were performed in situ up to pressures of ~6 GPa at room temperature. Structure factors were obtained and compared to MD simulations with an average error of less than 5% between major peak positions. The first sharp diffraction peak shifted from 4.56 Å to lower d-spacing with pressure, indicating compression between hard segments. This was further supported by the behavior of a peak at ~3.86 Å from the pair distribution function (PDF), suspected to represent π-stacking and separation between soft segments. Compression within the hard segments themselves is minimal as low-r peaks in the PDF are not greatly affected by pressure.

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

Influence of diluent concentration in localized high concentration electrolytes: elucidation of hidden diluent-Li + interactions and Li + transport mechanism

Localized high concentration electrolytes (LHCE) offer a viable dilution strategy for high concentration electrolytes (HCE) as the dilution process barely impacts the enhanced reductive/oxidative behavior of the HCE formulation but significantly lowers the overall viscosity and, in most cases, increases the ionic conductivity. On the other hand, experimental studies indicate that fluorinated ether electrolytes such as 1,1,2,2-tetrafluoroethylene 2,2,3,3-tetrafluoropropyl ether (TTE) help grow enhanced passivation layers on Ni-rich NMC cathodes. In this work, we study LHCE formulations based on lithium bis(fluorosulfonyl)imide (LiFSI), dimethyl carbonate (DMC), and TTE as the diluent. Here we use molecular dynamics methodologies, and Raman spectra measurements, to address to what extent the diluent content impacts the coordination behavior of the aggregated structures in the electrolyte and to evaluate the Li + transport properties under the influence of an external electric field. In contrast to other fluorinated ethers, we find that TTE interacts with Li + via fluorine atoms, partially limiting the DMC–Li + interactions hence altering the Li + solvation coordination. This competitive interaction with Li+ between the organic solvent and the TTE diluent influences the electrolyte's reductive/oxidative behavior. Nevertheless, the bonding strength of the Li + –F TTE is much weaker than those of the Li + –O DMC and Li + –OFSI – . Therefore, the existence of Li + –F TTE is in a transient state rather than in a steady state. These results provide plausible guiding rules for future dilution strategies of HCE electrolytes. We also demonstrate that Li + ions drift under the electric field's influence via repeated ion dissociation/association processes following a hopping conduction mechanism. Li + ions jump between aggregated networks where Li–O interactions dominate via diluent-enriched phases, a process in which the solvation shells temporarily mutate to a Li–F dominated coordination structure. We expect our results to contribute an improved atomic-level understanding of the solvation structure and dynamics of LHCE electrolytes.

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

Influence of isoelectronic mass modulation on phonon decay and transport in layered pnictogen-carbon systems

This study systematically investigates the structural, electronic, and thermal transport properties of layered pnictogen–carbon (Pn 2 C 2 ; Pn = P, As, Sb, Bi) monolayers, revealing a profound influence of isoelectronic mass modulation on phonon dynamics and thermal conductivity. Through first-principles calculations and lattice dynamics analysis, we demonstrate that increasing the atomic mass of pnictogens leads to elongated Pn-C bonds, weakened interatomic interactions, and enhanced phonon-phonon scattering, resulting in a dramatic reduction in lattice thermal conductivity (κ). Specifically, κ decreases exponentially from 65.6 W/m-K in P 2 C 2 to an ultralow 0.37 W/m-K in Bi 2 C 2 , driven by suppressed acoustic-optical phonon gaps Δ$^{Γ}_{A-0}$, increased anharmonicity, and reduced phonon lifetimes and group velocities. The transition from semiconducting to metallic behavior in heavier pnictogen-based systems further enhances phonon-electron scattering, contributing to thermal conductivity suppression. Structural analysis highlights the strengthening of out-of-plane C—C bonds and the contraction of C—Pn—C bond angles, which disrupt in-plane phonon transport. Further, these findings establish a design framework for engineering ultralow-κ materials through chemical substitution and structural tuning, offering significant potential for thermoelectric and thermal management applications. This work provides fundamental insights into phonon decay mechanisms and thermal transport in 2D materials, paving the way for advanced phononic and energy-efficient devices.

2D materials↗
DOE OSTI2023

Structure–Thermodynamic Relationship of a Polysaccharide Gel (Alginate) as a Function of Water Content and Counterion Type (Na vs Ca)

Biofilms are the predominant mode of microbial life on Earth, and so a deep understanding of microbial communities-and their impacts on environmental processes-requires a firm understanding of biofilm properties. Because of the importance of biofilms to their microbial inhabitants, microbes have evolved different ways of engineering and reconfiguring the matrix of extracellular polymeric substances (EPS) that constitute the main non-living component of biofilms. This ability makes it difficult to distinguish between the biotic and abiotic origins of biofilm properties. An important route toward establishing this distinction has been the study of simplified models of the EPS matrix. This study builds on such efforts by using atomistic simulations to predict the nanoscale (≤10 nm scale) structure of a model EPS matrix and the sensitivity of this structure to interpolymer interactions and water content. To accomplish this, we use replica exchange molecular dynamics (REMD) simulations to generate all-atom configurations of ten 3.4 kDa alginate polymers at a range of water contents and Ca–Na ratios. Simulated systems are solvated with explicitly modeled water molecules, which allows us to capture the discrete structure of the hydrating water and to examine the thermodynamic stability of water in the gels as they are progressively dehydrated. Our primary findings are that (i) the structure of the hydrogels is highly sensitive to the identity of the charge-compensating cations, (ii) the thermodynamics of water within the gels (specific enthalpy and free energy) are, surprisingly, only weakly sensitive to cation identity, and (iii) predictions of the differential enthalpy and free energy of hydration include a short-ranged enthalpic term that promotes hydration and a longer-ranged (presumably entropic) term that promotes dehydration, where short and long ranges refer to distances shorter or longer than ~0.6 nm between alginate strands.

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

Hydrogen Spillover Is Regulating Minority Rh 1 Active Sites on TiO 2 in Room-Temperature Ethylene Hydrogenation

The complicated dynamics of active sites on single-atom catalysts under reducing conditions limits their applications in hydrogenation reactions and mechanistic understanding. Herein, we report that on Rh 1 /TiO 2 , *H spillover during room-temperature ethylene hydrogenation hydroxylates and reduces TiO 2 , enhancing the intrinsic activity of Rh 1 by 9-fold. Spectroscopic and kinetic evidence suggests that the spillover of *H is suppressed by their facile reaction with C 2 H 4 , most of the spilled *H are nonreactive spectators, and >99% turnovers occur on a small subset (<20%) of exposed “active Rh 1 ”. Steady-state kinetics indicates competitive adsorption between H and C 2 H 4 , H 2 dissociation is the rate-determining step, and the apparent activation barrier (E a,app ) of the reaction is ~48 kJ/mol. The evolution of Rh 1 under H 2 was further tracked by spectroscopic and microscopic techniques at elevated temperatures. At 200 °C, more Rh 1 are exposed, but these Rh 1 are at least 5-fold less active than that of the “active Rh 1 ”. At 300 °C, Rh clusters derived from Rh 1 become the main active sites, shifting E a,app to 62 kJ/mol, characteristic of Rh nanoparticles. At ≥400 °C, larger and more active Rh particles in the strong metal–support interaction state are created. In conclusion, this work revealed the unexpected regulation effects of *H spillover on M 1 active sites under ambient conditions, differentiated the minority active M 1 sites, and demonstrated how the stability of M 1 under reducing atmospheres affects hydrogenation catalysis.

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

Uncovering fast solid-acid proton conductors based on dynamics of polyanion groups and proton bonding strength

Achieving high proton conductivity in inorganic solids is key for advancing many electrochemical technologies, including low-energy nano-electronics and energy-efficient fuel cells and electrolyzers. A quantitative understanding of the physical traits of a material that regulate proton diffusion is necessary for accelerating the discovery of fast proton conductors. In this work, we have mapped the structural, chemical and dynamic properties of solid acids to the elementary steps of the Grotthuss mechanism of proton diffusion. Our approach combines ab initio molecular dynamics simulations, analysis of phonon spectra and atomic structure calculations. We have identified the donor–hydrogen bond lengths and the acidity of polyanion groups as key descriptors of local proton transfer and the vibrational frequencies of the cation framework as the key descriptor of lattice flexibility. The latter facilitates rotations of polyanion groups and long-range proton migration in solid acid proton conductors. The calculated lattice flexibility also correlates with the experimentally reported superprotonic transition temperatures. Using these descriptors, we have screened the Materials Project database and identified potential solid acid proton conductors with monovalent, divalent and trivalent cations, including Ag + , Sr 2+ , Ba 2+ and Er 3+ cations, which go beyond the traditionally considered monovalent alkali cations (Cs + , Rb + , K + , and NH 4 + ) in solid acids.

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

Reaction processes at step edges on S-decorated Cu(111) and Ag(111) surfaces: MD analysis utilizing machine learning derived potentials

We report a variety of complexation, reconstruction, and sulfide formation processes can occur at step edges on the {111} surfaces of coinage metals (M) in the presence of adsorbed S under ultra-high vacuum conditions. Given the cooperative many-atom nature of these reaction processes, Molecular Dynamics (MD) simulation of the associated dynamics is instructive. However, only quite restricted Density Functional Theory (DFT)-level ab initio MD is viable. Thus, for M = Ag and Cu, we instead utilize the DeePMD framework to develop machine-learning derived potentials, retaining near-DFT accuracy for the M–S systems, which should have broad applicability. These potentials are validated by comparison with DFT predictions for various key quantities related to the energetics of S on M(111) surfaces. The potentials are then utilized to perform extensive MD simulations elucidating the above diverse restructuring and reaction processes at step edges. Key observations from MD simulations include the formation of small metal–sulfur complexes, especially MS 2 ; development of a local reconstruction at A-steps featuring an S-decorated {100} motif; and 3D sulfide formation. Additional analysis yields further information on the kinetics for metal–sulfur complex formation, where these complexes can strongly enhance surface mass transport, and on the propensity for sulfide formation.

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

Advances in Operando Electrocatalysis with High-Resolution Hard X-Ray Spectroscopy

Electrocatalysis unfolds at the interface—where solids, liquids, and gases meet to exchange charge, transform molecules, and set the foundation for technologies like CO 2 conversion, and water splitting. These interfacial regions are inherently dynamic, chemically diverse, and structurally heterogeneous. How atoms rearrange, oxidize, coordinate ligands, or adsorb reaction intermediates under electro-chemical potential defines a catalyst’s performance. Furthermore, capturing these changes with element specificity and electronic sensitivity under operando conditions remains one of the most critical challenges in the field.

Garcia-Esparza, Angel T. [SLAC National Accelerato↗
DOE OSTI2023

Probing reaction channels via reinforcement learning

Abstract Chemical reactions are dynamical processes involving the correlated reorganization of atomic configurations, driving the conversion of an initial reactant into a result product. By virtue of the metastability of both the reactants and products, chemical reactions are rare events, proceeding fleetingly. Reaction pathways can be modelled probabilistically by using the notion of reactive density in the phase space of the molecular system. Such density is related to a function known as the committor function, which describes the likelihood of a configuration evolving to one of the nearby metastable regions. In theory, the committor function can be obtained by solving the backward Kolmogorov equation (BKE), which is a partial differential equation (PDE) defined in the full dimensional phase space. However, using traditional methods to solve this problem is not practical for high dimensional systems. In this work, we propose a reinforcement learning based method to identify important configurations that connect reactant and product states along chemical reaction paths. By shooting multiple trajectories from these configurations, we can generate an ensemble of states that concentrate on the transition path ensemble. This configuration ensemble can be effectively employed in a neural network-based PDE solver to obtain an approximation solution of a restricted BKE, even when the dimension of the problem is very high. The resulting solution provides an approximation for the committor function that encodes mechanistic information for the reaction, paving a new way for understanding of complex chemical reactions and evaluation of reaction rates.

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

Optical nanofiber testbeds for benchmarking membrane-waveguide photonic integrated circuit platforms toward on-chip quantum inertial sensing

Recent advances in cold atom interferometry with optical and magnetic atom guides have set the stage for quantum inertial sensors capable of operating in dynamic environments. In this work, we present three key innovations—evanescent-field (EF) atom guides, optical nanofiber testbeds, and membrane-waveguide photonic integrated circuit (PIC) platforms—to advance EF-guided atom interferometry. First, we demonstrate EF atom guides on optical nanofiber testbeds, which serve as performance benchmarks for our membrane-waveguide PIC platforms. Second, we achieve low-power (⁠ ~ 5 mW) guiding of freely moving, laser-cooled 133 Cs atoms in two-color, traveling-wave EF optical dipole traps at the novel, heat-efficient magic wavelengths of 793 and 937 nm (i.e., “793/937-nm EF atom guides”). Concurrently, we design and fabricate membrane-waveguide PIC platforms for these EF atom guides; in our prior work, we showed that these structures safely accommodate 4–6 times the required optical trap power under vacuum and enable dense cold atom generation via magneto-optical trapping in the vicinity of the optical wavguide for efficient loading. Third, we verify preserved atomic coherence via microwave fields and EF-coupled Doppler-free Raman beams; to our knowledge, this is the first report of coherence fringes driven by co-propagating EF-coupled Raman beams with only 150 nW of total optical power. By providing a direct comparison between optical nanofiber testbeds and membrane-waveguide PIC platforms, our results lay critical groundwork for the on-chip realization of EF-guided atom interferometry and the development of fully integrated, compact, lightweight, and low-power quantum accelerometers and gyroscopes.

Orozco, Adrian [Sandia National Laboratories (SNL-↗
DOE OSTI2021

A.I. for nuclear physics

Nuclear science is concerned with the understanding of the nature of matter, its basic constituents and their interaction to form the elements and the properties we observe. This includes the forms of matter we see around us and also exotic forms such as those that existed in the first moments after the Big Bang and that exist today inside neutron stars. The techniques, tools, and expertise needed for nuclear physics (NP) research are therefore diverse in nature. State-of-the art accelerators are being developed to illuminate the dynamical basis of the core of the atom in terms of the fundamental constituents called quarks and gluons and to increase the number of isotopes with known properties. This scientific infrastructure is reaching scales and complexities that require computational methods for tasks such as anomaly detection in operational data. New methodologies are needed to detect anomalies and to optimize operating parameters, predict failures as well as to discover new optimization algorithms.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Molecular Structure and Electron Affinity of Metal-Solvent Complexes: Insights from Density Functional Theory Simulations

A molecular level understanding of the structure and energetics of the monovalent and divalent metal ion complexes is of great importance for development of next-generation batteries. Here, Density Functional Theory (DFT) simulations at the ωb97xD/6-31 + G(d,p) level of theory are performed to investigate the interaction of metal ions (Li + , Na + , K + , Mg 2+ , Ca 2+ , Zn 2+ ) with 26 organic solvent molecules. The reduction energetics (electron affinity and reduction potential) and structural responses of the solvent molecules and the molecular complexes are discussed. The DFT calculations are carried out to investigate the structure, energetics and electron affinities of chelated complexes of water (H 2 O), tetrahydrofuran (THF) and di-methoxy ethane (DME) solvent molecules. Additionally, ab initio dynamic simulations (AIMD) at 298 K using atom centered density matrix propagation (ADMP) formalism are performed to understand the spontaneous structure formation upon electron attachment of the metal ion-solvent complexes. The ADMP simulations indicate the decomposition of Mg + -(DME) 3 complex via cleavage of C–O bond of one of the three DME molecules indicating irreversible decomposition of DME in the presence of the Mg + radical. We believe that the data collected as part of this investigation serves as a library of fundamental knowledge towards a deeper understanding of the electrode-electrolyte interfacial reactions.

25 ENERGY STORAGE↗
DOE OSTI2022

Large-Scale Atomistic Simulations: Investigating Free Expansion

The experiment investigates free expansion of a supercritical fluid into a two-phase liquid-vapor coexistence region. A huge molecular dynamics simulation (6 billion Lennard-Jones atoms) was run on 5760 GPUs (33% of LLNL Sierra) using LAMMPS/Kokkos software. This improved visualization workflow and started preliminary simulations of aluminum using SNAP machine learning potential.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Neutron Scattering in High Magnetic Fields: Community Recommendations to Seize the Moment for Scientific Breakthroughs

Oak Ridge National Laboratory (ORNL) operates the world’s two most powerful neutron sources: the High Flux Isotope Reactor (HFIR) and the Spallation Neutron Source (SNS). These form a national user facility that provides neutron scattering instruments to probe and image static and dynamic properties of materials down to the atomic scale. Because the neutron carries a magnetic moment and the energies and wavelengths of neutrons match energy and length scales associated with magnetic materials, these facilities provide unique sensitivity to magnetism.

36 MATERIALS SCIENCE↗
NASA NTRS1980

Chemical and isotopic fractionation by grain size separates

Fractionation of refractory elements according to grain size is argued to occur during their growth. Two major modes should exist: (1) during thermal condensation sequences whenever the condensing phase (e.g. Mg2SiO4) does not alloy with the precondensed phase (e.g. MgAl2O4); (2) during accretion of gaseous atoms in the nonequilibrated interstellar medium. Processes dynamically sorting grains according to size (e.g. sedimentation) therefore are potentially capable of achieving fractionations normally attributed to separations of dust and gas. This paper considers the first mode during supernova condensation; however, it also can occur in an equilibrium solar condensation sequence owing to an overlooked freedom in that simplified description.

Clayton, D. D.↗
NASA NTRS1986

Solar concentrator materials development

Materials with potential applications in reflective and refractive solar dynamic concentrators are tested for resistance to atomic oxygen degradation. It is found that inorganic coatings such as MgF2, SiO(x), and ITO provide excellent protection for reflective surfaces while organic materials are much more susceptible to erosion and mass loss. Of the organic polymers tested, the silicones have the highest intrinsic resistance to atomic oxygen degradation.

Morel, D. E.↗
NASA NTRS1992

Stiffnesses by (TtN) ensemble molecular dynamics

Calculation of elastic constants requires initializing a lattice of atoms, an h matrix defining the computational cell boundaries, and the parameters defining the interatomic potentials, then numerically integrating the equations of motion for the atoms and the computational cell boundaries, calculating strains and stresses at each timestep, and updating at each timestep the average quantities required for calculation of fluctuations. Enforcements of the DYNAMO FORTRAN code included addition of documentation, improved memory management and data flow, and enabling use of several interatomic potentials. Documentation included headers for each subroutine, which document the parameters input to the routine, algorithms employed, and output.

Tibbits, Patrick↗