Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Shadow molecular dynamics for flexible multipole models

Shadow molecular dynamics provide an efficient and stable atomistic simulation framework for flexible charge models with long-range electrostatic interactions. Shadow molecular dynamics simulations are driven by approximate “shadow” Born–Oppenheimer potentials for which the exact charges and forces are directly accessible without relying on costly (and approximate) iterative solvers. While previous implementations have been limited to atomic monopole charge distributions, we extend this approach to flexible multipole models. We derive detailed expressions for the shadow energy functions, potentials, and force terms, explicitly incorporating monopole–monopole, dipole–monopole, and dipole–dipole interactions. In our formulation, both atomic monopoles and atomic dipoles are treated as extended dynamical variables alongside the propagation of the nuclear degrees of freedom. We demonstrate that introducing the additional dipole degrees of freedom preserves the stability and accuracy previously seen in monopole-only shadow molecular dynamics simulations. In addition, we present a shadow molecular dynamics scheme where the monopole charges are held fixed while the dipoles remain flexible. Our extended shadow dynamics provide a framework for stable, computationally efficient, and versatile molecular dynamics simulations involving long-range interactions between flexible multipoles. This is of particular current interest in combination with machine-learned interatomic potentials, including long-range electrostatic interactions.

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

Towards sharp error analysis of extended Lagrangian molecular dynamics

The extended Lagrangian molecular dynamics (XLMD) method provides a useful framework for reducing the computational cost of a class of molecular dynamics simulations with constrained latent variables. The XLMD method relaxes the constraints by introducing a fictitious mass ε for the latent variables and solving a set of singularly perturbed ordinary differential equations. While favorable numerical performance of XLMD has been demonstrated in several different contexts in the past decade, mathematical analysis of the method remains scarce. Here, we propose the first error analysis of the XLMD method in the context of a classical polarizable force field model. While the dynamics with respect to the atomic degrees of freedom are general and nonlinear, the key mathematical simplification of the polarizable force field model is that the constraints on the latent variables are given by a linear system of equations. We prove that when the initial value of the latent variables is compatible in a sense that we define, XLMD converges as the fictitious mass ε is made small with $\mathscr{O}$(ε) error for the atomic degrees of freedom and with $\mathscr{O}$($\sqrt{ε}$) error for the latent variables, when the dimension of the latent variable d' is 1. Furthermore, when the initial value of the latent variables is improved to be optimally compatible in a certain sense, we prove that the convergence rate can be improved to $\mathscr{O}$(ε) for the latent variables as well. Numerical results verify that both estimates are sharp not only for d'=1, but also for arbitrary d'. In the setting of general d', we do obtain convergence, but with the non-sharp rate of $\mathscr{O}$($\sqrt{ε}$) for both the atomic and latent variables.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Dynamics of plasma atomic layer etching: Molecular dynamics simulations and optical emission spectroscopy

Atomic layer etching is intrinsically dynamic as it involves sequential and repeated exposures of a surface to be etched with different species at different energies. The composition and structure of the near surface region change in both time and depth. Full understanding of this process requires resolving both temporal and spatial variations. In this work, we consider silicon (Si) atomic layer etching (ALE) by alternating exposure to chlorine gas (Cl 2 ) and argon ions (Ar + ). Molecular dynamics (MD) simulations are compared to experimental measurements with the aim of better understanding the dynamics of ALE and to test the simulation procedure. The simulations help to more fully interpret the experimental measurements. Further, optical emission measured just above the surface being etched can be related to etch products and can, therefore, be directly compared to simulation predictions. The simulations capture the measured initial product distribution leaving the surface and match the measured etch per cycle reasonably well. While simulations demonstrate the importance of ion-induced surface damage and mixing into a layer below the surface, the depth of which depends mainly on ion energy, the experiments suggest there is more Cl mixed into the layer than the MD procedure predicts.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Li-Metal Anode in Dilute Electrolyte LiFSI/TMP: Electrochemical Stability Using Ab Initio Molecular Dynamics

Ab initio molecular dynamics simulations were performed for Li + conducting electrolytes based on trimethyl phosphates (TMP) and lithium bis(fluorosulfonyl)imide (Li + FSI – ) salt in contact with a Li-metal electrode. We focused on the transient-state behavior at the electrolyte, interfacial electrolyte–Li-metal electrode, and lithium reference electrode–electrolyte–Li-metal electrode to study dynamics and activation energy barriers of the Li + ion, electrochemical and thermal stability of the interface electrode-electrolyte, and potential behavior of the Li-metal electrode, respectively. Our results show that in the most stable state, Li + ions are tetrahedrally coordinated to three TMPs and one FSI – . The inner solvation shell of a Li-ion is composed of three TMP and one FSI – in one contact ion-pair and four TMPs in a solvent-separated ion-pair. On the other hand, Li-ions transport through electrolyte cages takes place when they are coordinated with three or less molecules that could be a combination of TMPs and FSI – . The decomposition pathway of the LiFSI salt when in direct contact with the Li-metal anode starts with defluorination of FSI – , rapidly losing F– to the lithium surface, forming LiF species. The remaining FSO 2 NSO 2 –2 with the addition of 2e – from the Li-metal decomposes into SO 2 –2 and NFSO 2 –2 . SO 2 –2 deposits on the Li-surface and decomposes into Li 2 O and Li 2 S. The remaining NFSO 2 –2 defluorinates, losing F – ion to the lithium surface, resulting in LiF and the remaining NSO 2 –1 deposits on the lithium surface and decomposes in the following picoseconds, forming several binary compounds such as Li 3 N, Li 2 S and Li 2 O. In contrast, when the salt is solvated by the TMP molecules, avoiding a direct contact with the Li metal electrode, only one defluorination occurs, decomposing the FSI – into FSO 2 NSO 2 –2 and a F – . The two anions remain stable as they are solvated by the TMP molecules. Furthermore, we also analyzed the open circuit potential energy (OCPE) of the Li-metal electrode during the SEI formation. OCPE is calculated from the average local potential profile difference within the Li-metal electrode and a pristine Li-crystal reference electrode (LRE). When no SEI is formed, the Li-metal electrode has an average OCPE of +0.36 eV vs LRE. Due to the formation of a SEI, the Li-metal electrode has an average OCPE between -0.07 and -0.21 eV vs LRE. The OCPE of the Li-metal electrode decreases by ~0.42 eV when a SEI is formed.

25 ENERGY STORAGE↗
DOE OSTI2020

Temperature and composition dependent screw dislocation mobility in austenitic stainless steels from large-scale molecular dynamics

Abstract Extensive molecular dynamics simulations are performed to determine screw dislocation mobility in austenitic Fe 0.7 Ni x Cr 0.3-x stainless steels as a function of temperature ranging from 100 to 1300 K, resolved shear stress from 30 to 140 MPa, and Ni composition from 0.0 to 30.0 at%. These mobility data are fitted to a linear mobility law with a nonzero stress offset, referred to as the threshold stress. We find that both the linear drag coefficient and the threshold stress increase with Ni composition. The drag coefficient increases with temperature, whereas the threshold stress decreases with temperature. Based on these calculations, we determine fitting functions for the linear solute drag coefficient as a function of temperature and composition. The mobility laws determined in this study may serve to inform dislocation dynamics simulations pertinent to dislocation network evolution at elevated temperatures for a wide composition range of austenitic stainless steels.

Chu, Kevin (ORCID:000000031929710X)↗
DOE OSTI2023

Equation of state for He bubbles in W and model of He bubble growth and bursting near W{100} surfaces derived from molecular dynamics simulations

Abstract Molecular dynamics (MD) simulations are performed to derive an equation of state (EOS) for helium (He) bubbles in tungsten (W) and to study the growth of He bubbles under a W(100) surface until they burst. We study the growth as a function of the initial nucleation depth of the bubbles. During growth, successive loop-punching events are observed, accompanied by shifts in the depth of the bubble towards the surface. Subsequently, the MD data are used to derive models that describe the conditions that cause the loop punching and bursting events. Simulations have been performed at 500, 933, 1500, 2000, and 2500 K to fit the parameters in the models. To compute the pressure in the bubble at the loop punching and bursting events from the models, we derive an EOS for He bubbles in tungsten with an accompanying volume model to compute the bubble volume for a given number of vacancies ( $$N_\text {V}$$ N V ), He atoms ( $$N_\text {He}$$ N He ), and temperature ( T ). To derive the bubble EOS, we firstly derive the EOS for a free He gas. The derived free-gas EOS can accurately predict all MD data included in the analysis (which span up to 54 GPa at 2500 K). Subsequently, the bubble EOS is derived based on the free-gas EOS by correcting the gas density to account for the interaction between He and W atoms. The EOS for the bubbles is fitted to data from MD simulations of He bubbles in bulk W that span a wide range of gas density and sizes up to about 3 nm in diameter. The pressure of subsurface bubbles at the loop punching events as calculated using the bubble-EOS and the volume model agrees well with the pressure obtained directly from the MD simulations. In the loop punching model, for bubbles consisting of $$N_\text {V}$$ N V vacancies and $$N_\text {He}$$ N He helium atoms, the $$N_\text {He}/N_\text {V}$$ N He / N V ratio that causes the event, the resulting increase in $$N_\text {V}$$ N V , and the associated shift of the bubble depth are formulated as a function of $$N_\text {V}$$ N V and T . In the bursting model, a bubble must simultaneously reach a certain depth and $$N_\text {He}/N_\text {V}$$ N He / N V ratio in order to burst. The burst depth and $$N_\text {He}/N_\text {V}$$ N He / N V are also modeled as a function of $$N_\text {V}$$ N V and T . The majority of the loop punching events occur at bubble pressures between 20 and 60 GPa, depending on the bubble size and temperature. The larger the bubble and the higher the temperature, the lower the bubble pressure. Furthermore, our results indicate that at a higher temperature, a bubble can burst from a deeper region.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Full structural ensembles of intrinsically disordered proteins from unbiased molecular dynamics simulations

Abstract Molecular dynamics (MD) simulation is widely used to complement ensemble-averaged experiments of intrinsically disordered proteins (IDPs). However, MD often suffers from limitations of inaccuracy. Here, we show that enhancing the sampling using Hamiltonian replica-exchange MD (HREMD) led to unbiased and accurate ensembles, reproducing small-angle scattering and NMR chemical shift experiments, for three IDPs of varying sequence properties using two recently optimized force fields, indicating the general applicability of HREMD for IDPs. We further demonstrate that, unlike HREMD, standard MD can reproduce experimental NMR chemical shifts, but not small-angle scattering data, suggesting chemical shifts are insufficient for testing the validity of IDP ensembles. Surprisingly, we reveal that despite differences in their sequence, the inter-chain statistics of all three IDPs are similar for short contour lengths (< 10 residues). The results suggest that the major hurdle of generating an accurate unbiased ensemble for IDPs has now been largely overcome.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Predicting 1 H NMR relaxation in Gd 3+ -aqua using molecular dynamics simulations

Atomistic molecular dynamics simulations are used to predict 1 H NMR T 1 relaxation of water from paramagnetic Gd 3+ ions in solution at 25 °C. Simulations of the T 1 relaxivity dispersion function r 1 computed from the Gd 3+ – 1 H dipole–dipole autocorrelation function agree within ≃8% of measurements in the range f 0 ≃ 5 ↔ 500 MHz, without any adjustable parameters in the interpretation of the simulations, and without any relaxation models. Here the simulation results are discussed in the context of the Solomon-Bloembergen-Morgan inner-sphere relaxation model, and the Hwang-Freed outer-sphere relaxation model. Below f 0 ≲ 5 MHz, the simulation overestimates r 1 compared to measurements, which is used to estimate the zero-field electron-spin relaxation time. The simulations show potential for predicting r 1 at high frequencies in chelated Gd 3+ contrast-agents used for clinical MRI.

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

Structure and energetics of model amphiphilic molecules at the water liquid-vapor interface - A molecular dynamics study

A molecular dynamics study of adsorption of p-n-pentylphenol at infinite dilution at the water liquid-vapor interface is reported. The calculated free energy of adsorption is -8.8 +/- 0.7 kcal/mol, in good agreement with the experimental value of -7.3 kcal/mol. The transition between the interfacial region and the bulk solution is sharp and well-defined by energetic, conformational, and orientational criteria. At the water surface, the phenol head group is mostly immersed in aqueous solvent. The most frequent orientation of the hydrocarbon tail is parallel to the interface, due to dispersion interactions with the water surface. This arrangement of the phenol ring and the alkyl chain requires that the chain exhibits a kink. As the polar head group is being moved into the solvent, the chain length increases and the tail becomes increasingly aligned toward the surface normal, such that the nonpolar part of the molecule exposed to water is minimized. The same effect was achieved when phenol was replaced by a more polar head group, phenolate.

Pohorille, Andrew↗
DOE OSTI2021

An Ab Initio Multiple Cloning Method for Non-Adiabatic Excited-State Molecular Dynamics in NWChem

The recently developed ab-initio multiple cloning (AIMC) approach based on the multicongurational Ehrenfest (MCE) method provides a powerful and accurate way of describing the excited-state dynamics of molecular systems. The AIMC method is a controlled approximation to non-adiabatic dynamics with a particular strength in the proper description of decoherence effects due to branching of vibrational wavepackets at a level crossing. In this work, we report a new implementation of the AIMC algorithm in the open source NWChem computational chemistry program. The framework combines linear-response time-dependent density functional theory with Ehrenfest mean-eld theory to determine the equations of motion for classical trajectories. The multi-dimensional wave function is decomposed into a superposition of Gaussian coherent states guided by Ehrenfest trajectories (i.e., MCE approach), which can clone with fully quantum mechanical amplitudes and phases. By using an efficient time-derivative based non-adiabatic coupling approach within the AIMC method, all observable are calculated on-the-y in the non-adiabatic molecular dynamics process. As a representative example, we apply our implementation to study the ultrafast photoinduced electronic and vibrational energy transfer in a pyridine molecule. The effects of the cloning procedure on electronic and vibrational coherence, relaxation and unidirectional energy transfer are discussed. This new AIMC implementation provides a high-level non-adiabatic molecular dynamics framework for simulating photoexcited dynamics in complex molecular systems and experimentally relevant ultrafast spectroscopic probes, such as nonlinear coherent optical and X-ray signals.

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

Method of information entropy for convergence assessment of molecular dynamics simulations

The lack of a reliable method to evaluate the convergence of molecular dynamics simulations has contributed to discrepancies in different areas of molecular dynamics. Here, the method of information entropy is introduced to molecular dynamics for stationarity assessment. The Shannon information entropy formalism is used to monitor the convergence of the atom motion to a steady state in a continuous spatial domain and is also used to assess the stationarity of calculated multidimensional fields such as the temperature field in a discrete spatial domain. It is demonstrated in this work that monitoring the information entropy of the atom position matrix provides a clear indicator of reaching steady state in radiation damage simulations, non-equilibrium molecular dynamics thermal conductivity computations, and simulations of Poiseuille and Couette flow in nanochannels. A main advantage of the present technique is that it is non-local and relies on fundamental quantities available in all molecular dynamics simulations. Unlike monitoring average temperature, the technique is applicable to simulations that conserve total energy such as reverse non-equilibrium molecular dynamics thermal conductivity computations and to simulations where energy dissipates through a boundary as in radiation damage simulations. The method is applied to simulations of iron using the Tersoff/ZBL splined potential, silicon using the Stillinger–Weber potential, and to Lennard–Jones fluid. Its applicability to both solids and fluids shows that the technique has potential for generalization to other areas in molecular dynamics.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Neural Network Molecular Dynamics at Scale

Neural network molecular dynamics (NNMD) simulations could revolutionize atomistic modeling of materials with quantum-mechanical accuracy at a fraction of computational cost. However, popular NNMD frameworks are generally implemented for a single computing node, and conventional energy-based NN models still suffer from large time-to-solution (T2S), prohibiting the application of NNMD to challenging materials simulations encompassing large spatiotemporal scales. Consequently, no leadership-scale NNMD simulation has thus far been reported. Here, we present a scalable parallel NNMD software (RXMD-NN) based on our scalable reactive molecular dynamics simulation engine named RXMD. RXMD-NN has achieved high scalability up to 786,432 IBM BlueGene/Q cores involving 1.7 billion atoms. Furthermore, we have achieved 4.6-fold reduction of T2S by using a novel network that directly predicts atomic forces from feature vectors. Reduced T2S has for the first time allowed the study of large-scale off-stoichiometry effect in a widely used phase change material, Ge 2 Se 2 Te 5 , thereby resolving its "firstsharp diffraction peak mystery".

fist sharp diffraction peak↗
DOE OSTI2026

From bulk to surface: Structure and dynamics of amorphous alumina from deep potential molecular dynamics

Understanding the atomic-scale structure and dynamics of amorphous oxide surfaces is essential for interpreting their chemical reactivity, mechanical stability, and interfacial behavior, yet direct experimental characterization remains challenging. We employ Deep Potential (DP) molecular dynamics to generate large-scale, ab initio -quality models of amorphous Al 2 O 3 bulk glasses and melt-quenched free surfaces, enabling a quantitative analysis of both structure and relaxation dynamics with statistical confidence inaccessible to direct ab initio simulation. The trained DP model reproduces experimental liquid and glass structure, captures the cooling-rate dependence of the bulk glass transition, and corrects systematic biases in the polyhedral populations predicted by widely used classical force fields. At the free surface, mass density recovers to bulk values over ~10 Å, while local coordination requires a slightly wider subsurface region to fully converge. The outermost layer is oxygen-enriched, exhibits altered polyhedral connectivity with contracted Al–O bonds, and hosts a broad population of under-coordinated motifs (notably AlO 3 and OAl 2 ) whose abundances are governed by glass stability. These under-coordinated surface motifs exhibit distinct vibrational signatures and occur as locally paired Lewis acid and Brønsted base sites consistent with bond-valence compensation, yet remain spatially dispersed rather than aggregating into extended clusters. Despite this pronounced structural heterogeneity, surface relaxation and the glass-transition temperature remain comparable to their bulk counterparts, suggesting that the disordered surface is kinetically stable once formed. Together, these results establish a molecular-level picture of amorphous alumina surfaces and demonstrate the capability of machine-learned potentials to resolve structure–property relationships in disordered oxide interfaces.

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

Dynamics of Nanoscale Grain-Boundary Decohesion in Aluminum by Molecular-Dynamics Simulation

The dynamics and energetics of intergranular crack growth along a flat grain boundary in aluminum is studied by a molecular-dynamics simulation model for crack propagation under steady-state conditions. Using the ability of the molecular-dynamics simulation to identify atoms involved in different atomistic mechanisms, it was possible to identify the energy contribution of different processes taking place during crack growth. The energy contributions were divided as: elastic energy, defined as the potential energy of the atoms in fcc crystallographic state; and plastically stored energy, the energy of stacking faults and twin boundaries; grain-boundary and surface energy. In addition, monitoring the amount of heat exchange with the molecular-dynamics thermostat gives the energy dissipated as heat in the system. The energetic analysis indicates that the majority of energy in a fast growing crack is dissipated as heat. This dissipation increases linearly at low speed, and faster than linear at speeds approaching 1/3 the Rayleigh wave speed when the crack tip becomes dynamically unstable producing periodic dislocation bursts until the crack is blunted.

Yamakov, V.↗
DOE OSTI2021

The role of sub-surface hydrogen on CO 2 reduction and dynamics on Ni(110): An ab initio molecular dynamics study

The catalytic reduction of carbon dioxide is a crucial step in many chemical industrial reactions, such as methanol synthesis, the reverse water-gas shift reaction (rWGS) and formic acid synthesis. Here, we investigate the role of bulk hydrogen, where hydrogen atoms are found deep inside a metal surface as opposed to subsurface ones, upon CO 2 reduction over a Ni(110) surface using density functional theory (DFT) and ab initio molecular dynamics (AIMD) simulations. While it has previously been shown that subsurface hydrogen stabilizes CO 2 and can aid in overcoming reaction barriers, the role of bulk hydrogen is less studied and thus unknown with regard to CO 2 reduction. Furthermore, we find that the presence of bulk hydrogen can significantly alter the electronic structure of the Ni(110) surface, particularly the work function and d-band center, such that CO 2 adsorbs more strongly to the surface and is more easily reduced. Our results show an enhanced CO 2 dissociation in the presence of bulk hydrogen, shedding light on a hitherto underappreciated mechanistic pathway for CO 2 reduction on metal surfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗