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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Molecular Dynamics Simulations of Supercritical Carbon Dioxide and Water using TraPPE and SWM4-NDP Force Fields

The increased levels of carbon dioxide (CO 2 ) emissions due to the combustion of fossil fuels and the consequential impact on global climate change have made CO 2 capture, storage, and utilization a significant area of focus for current research. In most electrochemical CO 2 applications, water is used as a proton donor due to its high availability and mobility and use as a polar solvent. Additionally, supercritical CO 2 is a promising avenue for electrochemical applications due to its unique chemical and physical properties. Consequently, understanding the interactions between water and supercritical CO 2 is of great importance for future electrochemical applications. Molecular dynamics (MD) simulation is a powerful tool that enables atomistic-resolution dynamics of molecular systems, which can complement and guide future experimental investigations. This study employed atomistic MD to study the cosolubilities, codiffusivities, and structure of supercritical CO 2 and water systems, with a polarizable water model (SWM4-NDP) and a nonpolarizable CO 2 model (TraPPE). Additionally, ab initio MD simulations were used to better understand how atomistic polarizable/nonpolarizable models compare to explicit modeling of electron densities. The polarizable water model exhibited substantial improvement in water-associated properties. In conclusion, we anticipate the development of a compatible polarizable CO 2 model to yield similar improvement, providing a pathway for realizing novel high-pressure electrochemical systems.

25 ENERGY STORAGE↗
DOE OSTI2020

Minimal Experimental Bias on the Hydrogen Bond Greatly Improves Ab Initio Molecular Dynamics Simulations of Water

Experiment Directed Simulations (EDS) is a method within a class of techniques seeking to improve molecular simulations by minimally biasing the system Hamiltonian to reproduce certain experimental observables. In a previous application of EDS to ab initio molecular dynamics (AIMD) simulation based on electronic density functional theory (DFT), the AIMD simulations of water were biased to reproduce its experimentally derived solvation structure. In particular, by solely biasing the O-O pair correlation functions, other structural and dynamical properties that were not biased were improved. In this work, the hypothesis is tested that directly biasing the OH pair correlation (and hence the H-O∙∙∙H hydrogen bonding), will provide an even better improvement of DFT-based water properties in AIMD simulations. The logic behind this hypothesis is that for most electronic DFT descriptions of water the hydrogen bonding is known to be deficient due to anomalous charge transfer and over polarization in the DFT. Using recent advances to the EDS learning algorithm, we thus train a minimal bias on AIMD water that reproduces the O-H radial distribution function derived from the highly accurate MB-pol model of water. Finally, it is then confirmed that biasing the O-H pair correlation alone can lead to improved AIMD water properties, with structural and dynamical properties in even closer to experiment than the previous EDS-AIMD model.

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

Modeling protein structures from predicted contacts with modern molecular dynamics potentials: accuracy, sensitivity, and refinement

Protein structure prediction has become increasingly popular and successful in recent years. An essential step for fragment-free, template-free methods is the generation of a final three-dimensional protein model from a set of predicted amino acid contacts that are often described by interresidue pairwise atomic distances. Here we explore the use of modern, open-source molecular dynamics (MD) engines, which have been continually developed over the last three decades with all-atom Hamiltonians to model biomolecular structure and dynamics, to generate accurate protein structures starting from a set of inferred pairwise distances. Additionally, the ability of MD empirical physical potentials to correct inaccuracies in the predicted geometries is tested. We rigorously characterize the effect of modeling parameters on results, the effect of different amounts of error in the predicted distances on the final structures, and test the ability of post-processing analysis to sort the best models out of a set of statistical replicas. We find that with exact distances and with noisy distances, the method can produce excellent structural models, and that the molecular dynamics force field seems to help correct errors in distance predictions, resisting the effects of applied noise.

Davidson, Russ↗
NASA NTRS1987

Molecular dynamics of the water liquid-vapor interface

The results of molecular dynamics calculations on the equilibrium interface between liquid water and its vapor at 325 K are presented. For the TIP4P model of water intermolecular pair potentials, the average surface dipole density points from the vapor to the liquid. The most common orientations of water molecules have the C2 nu molecular axis roughly parallel to the interface. The distributions are quite broad and therefore compatible with the intermolecular correlations characteristic of bulk liquid water. All near-neighbor pairs in the outermost interfacial layers are hydrogen bonded according to the common definition adopted here. The orientational preferences of water molecules near a free surface differ from those near rigidly planar walls which can be interpreted in terms of patterns found in hexagonal ice 1. The mean electric field in the interfacial region is parallel to the mean polarization which indicates that attention cannot be limited to dipolar charge distributions in macroscopic descriptions of the electrical properties of this interface. The value of the surface tension obtained is 132 +/- 46 dyn/cm, significantly different from the value for experimental water of 68 dyn/cm at 325 K.

NASA Center ARC↗
NASA NTRS2016

Molecular Dynamics Modeling of Carbon Nanotube Composite Fracture Using ReaxFF

Carbon nanotube (CNT) fiber reinforced composites with specific tensile strengths and moduli approaching those of aerospace grade carbon fiber composites have recently been reported. This achievement was enabled by the emerging availability of high N/tex yarns in kilometer-scale quantities. While the production of this yarn is an impressive advance, its strength is still much lower than that of the individual CNTs comprising the yarn. Closing this gap requires understanding load transfer between CNTs at the nanometer dimensional scale. This work uses reactive molecular dynamics simulations to gain an understanding at the nanometer scale of the key factors that determine CNT nanocomposite mechanical performance, and to place more realistic upper bounds on the target properties. While molecular dynamics simulations using conventional force fields can predict elastic properties, the ReaxFF reactive forcefield can also model fracture behavior because of its ability to accurately describe bond breaking and formation during a simulation. The upper and lower bounds of CNT composite properties are investigated by comparing systems composed of CNTs continuously connected across the periodic boundary with systems composed of finite length CNTs. These lengths, effectively infinite for the continuous tubes and an aspect ratio of 13 for the finite length case, result from practical limitations on the number of atoms that can be included in a simulation. Experimentally measured aspect ratios are typically on the order of 100,000, so the calculated results should represent upper and lower limits on experimental mechanical properties. Finally, the effect of various degrees of covalent crosslinking between the CNTs and amorphous carbon matrix is considered to identify the amount of CNT-matrix covalent bonding that maximizes overall composite properties.

Jensen, Benjamin D.↗
DOE OSTI2021

Nonadiabatic Molecular Dynamics Study of the Relaxation Pathways of Photoexcited Cyclooctatetraene

In the current study, we present nonadiabatic (NAMD) and adiabatic molecular dynamics simulations of the transition-state dynamics of photoexcited cyclooctatetraene (COT). The equilibrium-state structure and absorption spectra are analyzed using the semiempirical Austin Model 1 potential. The NAMD simulations are obtained by a surface-hopping algorithm. We analyzed in detail an active excited to ground state relaxation pathway accompanied by an S 2 /S 3 (D 2d ) → S 1 (D 8h ) → S o (D 4h ) → S o (D 2d ) double-bond shifting mechanism. The simulated excitation lifetime is in good agreement with experiment. The first excited singlet state S1 plays a crucial role in the photochemistry. Here, the obtained critical molecular conformations, energy barrier, and transition-state lifetime results will provide a basis for further investigations of the bond-order inversion and photoswitching process of COT.

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

Stochastic symplectic reduced-order modeling for model-form uncertainty quantification in molecular dynamics simulations in various statistical ensembles

Here, this work focuses on the representation of model-form uncertainties in molecular dynamics simulations in various statistical ensembles. In prior contributions, the modeling of such uncertainties was formalized and applied to quantify the impact of, and the error generated by, pair-potential selection in the microcanonical ensemble (NVE). In this work, we extend this formulation and present a linear-subspace reduced-order model for the canonical (NVT) and isobaric (NPT) ensembles. The symplectic reduced-order basis is randomized on the tangent space of the Stiefel manifold to provide topological relationships and capture model-form uncertainty. Using the Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS), we assess the relevance of these stochastic reduced-order atomistic models on canonical problems involving a Lennard-Jones fluid and an argon crystal melt.

42 ENGINEERING↗
DOE OSTI2023

Near-surface damage and mixing in Si-Cl 2 -Ar atomic layer etching processes: Insights from molecular dynamics simulations

Silicon-chlorine-argon (Si-Cl 2 -Ar) atomic layer etching (ALE) is simulated using classical molecular dynamics (MD). The simulations provide a detailed view into the near-surface region during ALE processing. Bombardment of Ar + ions creates a mixed amorphous region that significantly differs from the picture of ideal ALE. There is also a significant change in the Si etch yield and the etch product distribution as a function of Ar + ion fluence. The Si etch yield is the highest at the beginning of the bombardment step but eventually decays to the physical sputtering yield. Atomic Cl and silicon chlorides are major etch products at the start of an ion bombardment step, but quickly decay. Atomic Si yields remain relatively constant as a function of Ar + ion fluence. As a result, a new schematic of Si-Cl 2 -Ar ALE is presented in order to emphasize the complex behavior observed in MD simulations.

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

Constant-potential molecular dynamics simulations of molten salt double layers for FLiBe and FLiNaK

We report the results of constant-potential molecular dynamics simulations of the double layer interface between molten 2LiF–BeF 2 (FLiBe) and 23LiF–6NaF–21KF (FLiNaK) fluoride mixtures and idealized solid electrodes. Employing methods similar to those used in studies of chloride double layers, we compute the structure and differential capacitance of molten fluoride electric double layers as a function of applied voltage. The role of molten salt structure is probed through comparisons between FLiBe and FLiNaK, which serve as models for strong and weak associate-forming salts, respectively. In FLiBe, screening involves changes in Be–F–Be angles and alignment of the oligomers parallel to the electrode, while in FLiNaK, the electric field is screened mainly by rearrangement of individual ions, predominantly the polarizable potassium cation.

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

Evaluating the active site-substrate interplay between x-ray crystal structure and molecular dynamics in chorismate mutase

Designing realistic quantum mechanical (QM) models of enzymes is dependent on reliably discerning and modeling residues, solvents, and cofactors important in crafting the active site microenvironment. Interatomic van der Waals contacts have previously demonstrated usefulness toward designing QM-models, but their measured values (and subsequent residue importance rankings) are expected to be influenceable by subtle changes in protein structure. Using chorismate mutase as a case study, this work examines the differences in ligand-residue interatomic contacts between an x-ray crystal structure and structures from a molecular dynamics simulation. Select structures are further analyzed using symmetry adapted perturbation theory to compute ab initio ligand-residue interaction energies. The findings of this study show that ligand-residue interatomic contacts measured for an x-ray crystal structure are not predictive of active site contacts from a sampling of molecular dynamics frames. In addition, the variability in interatomic contacts among structures is not correlated with variability in interaction energies. However, the results spotlight using interaction energies to characterize and rank residue importance in future computational enzymology workflows.

Chemistry↗
DOE OSTI2023

Probing the Structure and Dynamics of Fluid Mixtures in Porous Materials Through Ultrafast Vibrational Spectro-Microscopy and Many-Body Molecular Dynamics

This research program focused on the characterization of both molecular structure and dynamics of aqueous solutions in various metal-organic frameworks (MOFs) and organic polymers that have recently been proposed for applications in water treatment technologies. This was accomplished by integrating ultrafast vibrational spectroscopy and microscopy, which are sensitive to local environments, molecular orientation, and dynamical couplings, with many-body molecular dynamics (MB-MD) simulations, which enable realistic modeling of aqueous systems and porous materials under different thermodynamic conditions. Our studies allowed for identifying the physical mechanisms and characterizing the underlying molecular interactions between water and host materials that determine the adsorption and transport processes of water in various MOFs and β-cyclodextrin polymers – prototypical examples of porous materials and polymers, that have recently been proposed for technological applications in water purification. By combining bulk- and surface-sensitive, spatially resolved ultrafast vibrational spectroscopy with MB-MD simulations, it was possible to gain broad insight into the structure and mobility of water under heterogeneous confinement at various length scales. In particular, measurements of surface- and bulk-sensitive vibrational spectra at single crystal level, integrated with MB-MD simulations were used to characterize domain-specific, structure–binding affinity relationships. Besides providing specific information about the behavior of aqueous solutions in prototypical MOFs and organic polymers for potential applications in water treatment technologies, our studies also provided fundamental insights into the properties of aqueous solutions in confined heterogeneous environments, which have broad implications in many areas, ranging from heterogeneous catalysis to oil recovery, ion transport processes, and reaction-diffusion processes in crowded systems.

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

Ab initio molecular dynamics (AIMD) simulations of NaCl, UCl 3 molten salts

Ab initio molecular dynamics (AIMD) simulations are used to calculate select thermophysical and thermodynamic properties of NaCl, UCl and NaCl-UCl 3 molten salts. Following established approaches, the AIMD simulations include a model for Van der Waals interactions. The Langreth & Lundqvist (vDW-DF2), DFT-D3, and density-dependent energy correction (DFT-dDsC) dispersion models are first tested for molten NaCl in order to assess their accuracy for density and heat capacity predictions across a range of temperatures. Based on the NaCl results, the vdW-DF2 and DFT-dDsC methods are extended to the UCl 3 system and compared to available experimental data. Next, mixtures of NaCl-UCl 3 are investigated and analyzed with respect to the deviation from ideal solution behavior. For the DFT-dDsC simulations, density deviates by up to 2% from an ideal mixture, with the maximum occurring close to the eutectic composition. The mixing energy also deviates from an ideal solution and exhibits a minimum of about -0.074 to eV per formula unit, again close to the eutectic composition. Finally, the compressibility and species diffusivity of the pure and mixed salt systems are calculated. The diffusivities are slightly reduced in the mixed compared to the pure systems and the compressibilities loosely follow a linear correlation as a function of the UCl 3 composition. The trends observed for mixing properties are rationalized by correlating them to the coordination chemistry, which emphasizes the importance of maintaining the extended network formed by U and Cl ions as the NaCl concentration increases. The concentration at which breakdown of the extended network occurs, roughly coincides with the minimum of both the mixing energy and the deviation from ideal solution behavior for density.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Molecular dynamics simulations of displacement cascades in LiAlO2 and LiAl5O8 ceramics

Abstract Molecular dynamics was employed to investigate the radiation damage due to collision cascades in LiAlO 2 and LiAl 5 O 8 , the latter being a secondary phase formed in the former during irradiation. Atomic displacement cascades were simulated by initiating primary knock-on atoms (PKA) with energy values = 5, 10 and 15 keV and the damage was quantified by the number of Frenkel pairs formed for each species: Li, Al and O. The primary challenges of modeling an ionic system with and without a core–shell model for oxygen atoms were addressed and new findings on the radiation resistance of these ceramics are presented. The working of a variable timestep function and the kinetics in the background of the simulations have been elaborated to highlight the novelty of the simulation approach. More importantly, the key results indicated that LiAlO 2 experiences much more radiation damage than LiAl 5 O 8 , where the number of Li Frenkel pairs in LiAlO 2 was 3–5 times higher than in LiAl 5 O 8 while the number of Frenkel pairs for Al and O in LiAlO 2 are ~ 2 times higher than in LiAl 5 O 8 . The primary reason is high displacement threshold energies (E d ) in LiAl 5 O 8 for Li cations. The greater E d for Li imparts higher resistance to damage during the collision cascade and thus inhibits amorphization in LiAl 5 O 8 . The presented results suggest that LiAl 5 O 8 is likely to maintain structural integrity better than LiAlO 2 in the irradiation conditions studied in this work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

An examination of the performance of molecular dynamics force fields: Silicon and silicon dioxide reactive ion etching

Classical molecular dynamics (MD) simulations of plasma-surface interactions were performed of physical sputtering and reactive ion etching (RIE), with predictions based on several force fields. In this paper, we focus mainly on SiO 2 but include some results for Si substrates as well. Here, we compare predictions from these MD simulations to experimental studies of SiO 2 physical sputtering (by Ar + ions), RIE of Si, and RIE of SiO 2 (both using F atoms and Ar + ions). MD results using different published force fields are compared to reported yields from published vacuum beam experiments. The near-surface depth profiles predicted using different force fields are compared. One motivation for the present study is to document the nature and magnitude of differences in the predictions for selected systems and conditions of practical interest.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1989

A model for including thermal conduction in molecular dynamics simulations

A technique is introduced for including thermal conduction in molecular dynamics simulations for solids. A model is developed to allow energy flow between the computational cell and the bulk of the solid when periodic boundary conditions cannot be used. Thermal conduction is achieved by scaling the velocities of atoms in a transitional boundary layer. The scaling factor is obtained from the thermal diffusivity, and the results show good agreement with the solution for a continuous medium at long times. The effects of different temperature and size of the system, and of variations in strength parameter, atomic mass, and thermal diffusivity were investigated. In all cases, no significant change in simulation results has been found.

Wu, Yue↗
NASA NTRS1991

Molecular Dynamics study of Pb overlayer on Cu(100)

Isothermal-isobaric Molecular Dynamics (MD) simulation of a submonolayer Pb film in c(2x2) ordered structure adsorbed on a Cu(100) substrate showed retention of order to high T. The Embedded Atom Method (EAM) calculated the energy of atoms of overlayer and substrate. The time-averaged squared modulus of the two dimensional structure factor for the Pb overlayer measured the order of the overlayer. The results are for increasing T only, and require verification by simulated cooling.

Karimi, M.↗
DOE OSTI2025

Data from "Deep Potential Molecular Dynamics Simulations of Low-Temperature Plasma-Surface Interactions"

Data and input files related to the paper "Deep Potential Molecular Dynamics Simulations of Low-Temperature Plasma-Surface Interactions" (https://doi.org/10.1116/6.0004027). This includes the final DP model used in all simulations, training data set, example input files to run DeepMD (with LAMMPS), and data tables summarizing the results obtained from the simulations.

machine learning models↗