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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

KCl-UCl 3 molten salts investigated by Ab Initio Molecular Dynamics (AIMD) simulations: A comparative study with three dispersion models

Ab Initio Molecular Dynamics (AIMD) simulations are performed on molten KCl-UCl 3 salt mixtures to determine energies, heat capacities, and densities. The density-dependent energy correction (DFT-dDsC), Grimme et al.’s DFT-D3, and Langreth & Lundqvist (vdW-cx) models are used for dispersion forces and combined with the Perdew-Burke-Ernzerhof (PBE) exchange-correlation potential with a Hubbard U param eter for the 5$f$ electrons of uranium. After validating predictions for the end-member systems to literature data, KCl-UCl 3 mixtures are studied at select temperatures. Densities and energies both deviate from ideal solution behavior, with the maximum deviation occurring around 36% UCl 3 for mixing energies and slightly lower (29% UCl 3 ) for densities. Compared to the NaCl-UCl 3 system, which was previously investigated using the same simulation methodologies, the KCl-UCl 3 density and mixing energy deviations from ideal solution behavior are larger by almost a factor of two. No deviation from ideal solution behavior for heat capacity was observed. The AIMD predictions for mixing energies and densities agree qualitatively with experimental data, though the spread in data obtained from the various dispersion force models utilized, measurements, and empirical estimates makes strong conclusions difficult. The dependence of thermodynamic and thermophysical properties on composition is correlated with the local chemistry of the solution phase, in particular, the tendency of UCl 3 to form network structures.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1981

Cluster adsorption on amorphous and crystalline surfaces - A molecular dynamics study of model Pt on Cu and model Pd on Pt

Molecular dynamics was used to study the structure, dispersion and short-time behavior of ten-atom clusters adsorbed onto amorphous and crystalline substrates, in which the cluster atoms differed from the substrate atoms. Two adatom-substrate model systems were chosen; one, in which the interaction energy between adatom pairs was greater than that between substrate pairs, and the other, in which the reverse was true. At relatively low temperature ranges, increased dispersion of cluster atoms occurred: (a) on the amorphous substrate as compared to the FCC(100) surface, (b) with increasing reduced temperature, and (c) with adatom-substrate interaction energy stronger than adatom-adatom interaction. Two-dimensional clusters (rafts) on the FCC(100) surface displayed migration of edge atoms only, indicating a mechanism for the cluster rotation and shape changes found in experimental studies.

Garofalini, S. H.↗
NASA NTRS1994

Molecular Dynamics of a Water-Lipid Bilayer Interface

We present results of molecular dynamics simulations of a glycerol 1-monooleate bilayer in water. The total length of analyzed trajectories is 5ns. The calculated width of the bilayer agrees well with the experimentally measured value. The interior of the membrane is in a highly disordered fluid state. Atomic density profile, orientational and conformational distribution functions, and order parameters indicate that disorder increases toward the center of the bilayer. Analysis of out-of-plane thermal fluctuations of the bilayer surfaces occurring at the time scale of the present calculations reveals that the distribution of modes agrees with predictions of the capillary wave model. Fluctuations of both bilayer surfaces are uncorrelated, yielding Gaussian distribution of instantaneous widths of the membrane. Fluctuations of the width produce transient thinning defects in the bilayer which occasionally span almost half of the membrane. The leading mechanism of these fluctuations is the orientational and conformational motion of head groups rather than vertical motion of the whole molecules. Water considerably penetrates the head group region of the bilayer but not its hydrocarbon core. The total net excess dipole moment of the interfacial water points toward the aqueous phase, but the water polarization profile is non-monotonic. Both water and head groups significantly contribute to the surface potential across the interface. The calculated sign of the surface potential is in agreement with that from experimental measurements, but the value is markedly overestimated. The structural and electrical properties of the water-bilayer system are discussed in relation to membrane functions, in particular transport of ions and nonelectrolytes across membranes.

Wilson, Michael A.↗
DOE OSTI2024

Neutron scattering and neural-network quantum molecular dynamics investigation of the vibrations of ammonia along the solid-to-liquid transition

Abstract Vibrational spectroscopy allows us to understand complex physical and chemical interactions of molecular crystals and liquids such as ammonia, which has recently emerged as a strong hydrogen fuel candidate to support a sustainable society. We report inelastic neutron scattering measurement of vibrational properties of ammonia along the solid-to-liquid phase transition with high enough resolution for direct comparisons to ab-initio simulations. Theoretical analysis reveals the essential role of nuclear quantum effects (NQEs) for correctly describing the intermolecular spectrum as well as high energy intramolecular N-H stretching modes. This is achieved by training neural network models using ab-initio path-integral molecular dynamics (PIMD) simulations, thereby encompassing large spatiotemporal trajectories required to resolve low energy dynamics while retaining NQEs. Our results not only establish the role of NQEs in ammonia but also provide general computational frameworks to study complex molecular systems with NQEs.

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

Deformation behavior of annealed Cu64Zr36 metallic glass via molecular dynamics simulations

In the present work, we investigate the isothermal annealing process of Cu 64 Zr 36 metallic glass (MG) by means of molecular dynamics (MD) simulations, and characterize the compression properties under different temperatures and compressive rates. The pair distribution functions, the bond-orientational order, and the cluster-type index method are modeled and analyzed to characterize changes in the structure. Results of the modeling and analysis reveal that the MgZn 2 -type Laves, MgCu 2 -type Laves and five-fold local symmetry structures can be formed during the isothermal holding at 950 K. Study on the compression properties shows that at a fixed strain rate, the yield strength and yield drop all increase with decreasing temperature. At a fixed temperature, however, the increase of the strain rate leads to a noticeable increase in the yield strength at 600 K, but has little effect on yield strength at 10 K. Moreover, the modeling and analysis of the structure at a temperature of 10 K and strain-rate of 5 × 10 7 s -1 demonstrate that the order-disorder transformation initiates the shear band.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Autonomous alloy composition optimization using molecular dynamics guided by a large language model

Here, we present an autonomous materials discovery framework that couples a large language model (LLM) with molecular dynamics (MD) simulations to optimize Fe–Cr–Mn alloy compositions for tensile strength. Starting from six distinct compositions, the LLM operated as an intelligent agent, iteratively proposing changes based on prior simulation results and constraints. Over 50 iterations per case, the LLM adaptively explored the composition space, identifying high-strength regions, not easily accessible by conventional methods. The highest strength, 18.7 GPa, was achieved with Fe 71 Cr 25 Mn 4 composition, identified from a Fe 75 Cr 20 Mn 5 starting point. The LLM autonomously adjusted its strategy in real time, demonstrating closed-loop decision-making using commodity hardware. This approach showcases the potential of LLMs as scientific co-pilots, capable of accelerating materials discovery and generalizable to other domains like biology and drug design.

Autonomy↗
DOE OSTI2022

Investigation of γ - ( U , Z r ) structural properties and its interfacial properties with liquid sodium using ab initio molecular dynamics

In this study, the elastic properties, structural parameters, sound velocity, and Debye temperature of γ–(U,Zr) were computed using ab initio molecular dynamics (AIMD) at temperatures between 1000 K and 1400 K and for Zr content between 0 at. % and 100 at. %. UZr is used as a metallic fuel for Sodium Fast Reactors (SFRs). The study of the mechanical and thermal behavior of these alloys leads to a better data-informed fuel design. The bulk modulus, shear modulus, Young's modulus, and Poisson's ratio were calculated from the elastic constants and their dependence on Zr content and temperature was investigated, comparing the results with previous computational work and the available experimental data in the literature. Interfacial properties between UZr (up to 32 at. % which typically exists in nuclear fuel) and liquid sodium are also of interest due to the presence of a sodium bond between the fuel and the cladding in metallic nuclear fuel. The interfacial energy between γ–(U,Zr) and liquid sodium, the surface tension of liquid sodium, and the work of adhesion were computed at different temperatures and Zr concentrations. It was demonstrated that γ–(U,Zr) is completely wetted by liquid sodium at all the investigated temperatures and Zr concentrations. Finally, this work provides the basis for the determination of interfacial resistances in SFRs and their implementation into heat transfer fuel performance simulations, which will be the subject of future work.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Broadband NMR Relaxometry as a Powerful Technique to Study Molecular Dynamics of Ionic Liquids

Fast field cycling nuclear magnetic resonance (FFC NMR) relaxometry technique has been demonstrated to be a useful analytical tool to investigate molecular dynamics in very diverse systems during the last decades. Of particular importance has been its application in studying ionic liquids, upon which this review article is based. Furthermore, some of the research carried out on ionic liquids during the last ten years using this technique is highlighted in this article with the aim of promoting the favorable features of FFC NMR applied toward understanding dynamics of complex systems.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Molecular dynamics simulation of seeded crystal growth in glass

The mechanism of non-congruent growth of a crystal from glass has been sought using molecular dynamics simulations. Specifically, as a model of this process, the growth of a lithium niobate (LiNbO 3 ) crystal seed sandwiched between two lithium niobosilicate (LNS) glass slabs has been simulated as a function of time and temperature. The growth of pre-existing crystal is strongly affected by the orientation of crystal seed, temperature, and the SiO 2 concentration in the surrounding LNS glass matrix. The orientation of LiNbO 3 seed surface that has inherently larger interplanar distance results in a relatively slower crystal growth. The addition of SiO 2 to LNS system significantly decreases the crystal growth, which primarily occurs in the region devoid of Si. The suppressive effect of SiO 2 on growth rate can be traced to the existence of defect complex comprising of Si substituted at the Nb site and a nearby Nb vacancy.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Molecular dynamics studies of helium bubble effects on grain boundary fracture vulnerabilities in an Fe 70 Ni 11 Cr 19 –1%H austenitic stainless steel

Comprehensive molecular dynamics tensile test simulations have been performed to study the delamination processes of seven different grain boundaries / cleavage planes (Σ1{111}, Σ3{111}, Σ5{100}, Σ7{111}, Σ9{411}, Σ11{311}, and R{100}/{411}) containing a helium bubble. Combinations of a variety of conditions are explored including different strain rates, system dimensions, bubble density, bubble radius, bubble pressure, and temperature. We found that in general, grain boundaries absorb less energies with decreasing strain rate but increasing bubble areal density, bubble pressure, bubble radius, and temperature. The propensity of grain boundary delamination is sensitive to grain boundary type: The random grain boundary R{100}/{411} is one of the most brittle boundaries whereas the Σ1{111} cleavage plane and the Σ3{111} twin boundary are two of the toughest boundaries. The sorted list of grain boundary fracture vulnerability obtained from our dynamic tensile test simulations differs from the one obtained from our decohesion energy calculations, confirming the important role of plastic deformation during fracture. Finally, detailed mechanistic analyses are performed to interpret the simulated results.

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

Modeling Chemical Reactions in Classical Molecular Dynamics Simulations

An algorithm capable of incorporating multi-step reaction mechanisms into atomistic molecular dynamics (MD) simulations using traditional fixed valence force fields is proposed and implemented within the framework of LAMMPS (Large-scale Atomic Molecular Massively Parallel Simulator). This extension, referred to as fix bond/react, enables bonding topology modifications during a running MD simulation using pre- and post-reaction bonding templates to carry out a pre-specified reaction. Candidate reactants are first identified by interatomic separation, followed by the application of a generalized topology matching algorithm to confirm they match the pre-reaction template. This is followed by a topology conversion to match the post-reaction template and a dynamic relaxation to minimize high energy configurations. Two case studies, the condensation polymerization of nylon 6,6 and the formation of a highly-crosslinked epoxy, are simulated to demonstrate the robustness, stability, and speed of the algorithm. Improvements which could increase its utility are discussed.

Jacob R Gissinger↗
DOE OSTI2023

Evaluation of 6-OxP-CD, an Oxime-based cyclodextrin as a viable medical countermeasure against nerve agent poisoning: Experimental and molecular dynamic simulation studies on its inclusion complexes with cyclosarin, soman and VX

The ability of the cyclodextrin-oxime construct 6-OxP-CD to bind and degrade the nerve agents Cyclosarin (GF), Soman (GD) and S -[2-[Di(propan-2-yl)amino]ethyl] O -ethyl methylphosphonothioate (VX) has been studied using 31 P-nuclear magnetic resonance (NMR) under physiological conditions. While 6-OxP-CD was found to degrade GF instantaneously under these conditions, it was found to form an inclusion complex with GD and significantly improve its degradation (t 1/2 ~ 2 hrs) relative over background (t 1/2 ~ 22 hrs). Consequently, effective formation of the 6-OxP-CD:GD inclusion complex results in the immediate neutralization of GD and thus preventing it from inhibiting its biological target. In contrast, NMR experiments did not find evidence for an inclusion complex between 6-OxP-CD and VX, and the agent’s degradation profile was identical to that of background degradation (t 1/2 ~ 24 hrs). As a complement to this experimental work, molecular dynamics (MD) simulations coupled with Molecular Mechanics-Generalized Born Surface Area (MM-GBSA) calculations have been applied to the study of inclusion complexes between 6-OxP-CD and the three nerve agents. These studies provide data that informs the understanding of the different degradative interactions exhibited by 6-OxP-CD with each nerve agent as it is introduced in the CD cavity in two different orientations (up and down). For its complex with GF, it was found that the oxime in 6-OxP-CD lies in very close proximity (P GF …O Oxime ~ 4–5 Å) to the phosphorus center of GF in the ‘down GF ’ orientation for most of the simulation accurately describing the ability of 6-OxP-CD to degrade this nerve agent rapidly and efficiently. Further computational studies involving the center of masses (COMs) for both components (GF and 6-OxP-CD) also provided some insight on the nature of this inclusion complex. Distances between the COMs (ΔCOM) lie closer in space in the ‘down GF ’ orientation than in the ‘up GF ’ orientation; a correlation that seems to hold true not only for GF but also for its congener, GD. In the case of GD, calculations for the ‘down GD ’ orientation showed that the oxime functional group in 6-OxP-CD although lying in close proximity (P GD …O Oxime ~ 4–5 Å) to the phosphorus center of the nerve agent for most of the simulation, adopts another stable conformation that increase this distance to ~ 12–14 Å, thus explaining the ability of 6-OxP-CD to bind and degrade GD but with less efficiency as observed experimentally (t 1/2 ~ 4 hr. vs. immediate). Lastly, studies on the VX:6-OxP-CD system demonstrated that VX does not form a stable inclusion complex with the oxime-bearing cyclodextrin and as such does not interact in a way that is conducive to an accelerated degradation scenario. Collectively, these studies serve as a basic platform from which the development of new cyclodextrin scaffolds based on 6-OxP-CD can be designed in the development of medical countermeasures against these highly toxic chemical warfare agents.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Stability of Calcium Ion Battery Electrolytes: Predictions from Ab Initio Molecular Dynamics Simulations

Multivalent batteries, such as magnesium-ion, calcium-ion, and zinc-ion batteries, have attracted significant attention as next-generation electrochemical energy storage devices to complement conventional lithium-ion batteries (LIBs). Among them, calcium-ion batteries (CIBs) are the least explored due to the difficult reversible Ca deposition-dissolution. In this work, we examined the stability of four different Ca salts with weakly coordinating anions and three different solvents commonly employed in existing battery technologies to identify suitable candidates for CIBs. By employing Born-Oppenheimer molecular dynamics (BOMD) simulations on salt-Ca and solvent-Ca interfaces, we find that the tetraglyme solvent and carborane salt are promising candidates for CIBs. Due to the strong reducing nature of the calcium surface, the other salts and solvents readily decompose. Further, we explain the microscopic mechanisms of salt/solvent decomposition on the Ca surface using time-dependent projected density of states, time-dependent charge-transfer plots, and climbing-image nudged elastic band calculations. Collectively, this work presents the first mechanistic assessment of the dynamical stability of candidate salts and solvents on a Ca surface using BOMD simulations, and provides a predictive path toward designing stable electrolytes for CIBs.

Born-Oppenheimer molecular dynamics↗
DOE OSTI2021

Molecular dynamics simulations of a dicationic ionic liquid for CO 2 capture

A dicationic ionic liquid ([DBU-PEG][Tf 2 N] 2 ) was studied using classical molecular dynamics simulations to examine its structural and gas separation properties. The dication was designed in an attempt to improve CO 2 solubility by means of tuning the cation-anion interactions of the ionic liquid (IL). The computational model was compared to experimentally obtained density, viscosity, and powder X-ray diffraction spectra. The structure of the IL was further investigated with radial distribution functions and free volume analysis through cavity distributions. It was found that the shape and charge distribution of the dication enhances CO 2 interaction: the CO 2 molecule is hugged by the dication along the PEG linker and close to one of the cationic ends. Here, the geminal design of the dication allows for strong interaction with CO 2 , showing promise as a means of carbon capture.

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

ReaxFF molecular dynamics simulations of electrolyte–water systems at supercritical temperature

We have performed ReaxFF molecular dynamics simulations of alkali metal–chlorine pairs in different water densities at supercritical temperature (700 K) to elucidate the structural and dynamical properties of the system. The radial distribution function and the angular distribution function explain the inter-ionic structural and orientational arrangements of atoms during the simulation. The coordination number of water molecules in the solvation shell of ions increases with an increase in the radius of ions. We find that the self-diffusion coefficient of metal ions increases with a decrease in density under supercritical conditions due to the formation of voids within the system. The hydrogen bond dynamics has been interpreted by the residence time distribution of various ions, which shows Li+ having the highest water retaining capability. The void distribution within the system has been analyzed by using the Voronoi polyhedra algorithm providing an estimation of void formation within the system at high temperatures. We observe the formation of salt clusters of Na+ and K+ at low densities due to the loss of dielectric constants of ions. The diffusion of ions gets altered dramatically due to the formation of voids and nucleation of ions in the system.

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

Molecular dynamics simulations of inelastic x-ray scattering from shocked copper

By taking the spatial and temporal Fourier transforms of the coordinates of the atoms in molecular dynamics simulations conducted using an embedded-atom-method potential, we calculate the inelastic scattering of x rays from copper single crystals shocked along [001] to pressures of up to 70 GPa. Above the Hugoniot elastic limit, we find that the copious stacking faults generated at the shock front introduce strong quasi-elastic scattering (QES) that competes with the inelastic scattering signal, which remains discernible within the first Brillouin zone; for specific directions in reciprocal space outside the first zone, the QES dominates the inelastic signal overwhelmingly. The synthetic scattering spectra we generate from our Fourier transforms suggest that energy resolutions of order 10 meV would be required to distinguish inelastic from quasi-elastic scattering within the first Brillouin zone of shock-loaded copper. We further note that high-resolution inelastic scattering also affords the possibility of directly measuring particle velocities via the Doppler shift. These simulations are of relevance to future planned inelastic scattering experiments at x-ray Free Electron Laser facilities.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Understanding Hydrogenation Chemistry at MgB 2 Reactive Edges from Ab Initio Molecular Dynamics

Solid-state hydrogen storage materials often operate via transient, multistep chemical reactions at complex interfaces that are difficult to capture. Here, we use direct ab initio molecular dynamics simulations at accelerated temperatures and hydrogen pressures to probe the hydrogenation chemistry of the candidate material MgB 2 without a priori assumption of reaction pathways. Focusing on highly reactive ($10\overline{1}0$) edge planes where initial hydrogen attack is likely to occur, we track mechanistic steps toward the formation of hydrogen-saturated BH 4 – units and key chemical intermediates, involving H 2 dissociation, generation of functionalities and molecular complexes containing BH 2 and BH 3 motifs, and B–B bond breaking. Further, the genesis of higher-order boron clustering is also observed. Different charge states and chemical environments at the B-rich and Mg-rich edge planes are found to produce different chemical pathways and preferred speciation, with implications for overall hydrogenation kinetics. The reaction processes rely on B–H bond polarization and fluctuations between ionic and covalent character, which are critically enabled by the presence of Mg 2+ cations in the nearby interphase region. Our results provide guidance for devising kinetic improvement strategies for MgB 2 -based hydrogen storage materials, while also providing a template for exploring chemical pathways in other solid-state energy storage reactions.

08 HYDROGEN↗
DOE OSTI2020

Cavity molecular dynamics simulations of liquid water under vibrational ultrastrong coupling

We simulate vibrational strong coupling (VSC) and vibrational ultrastrong coupling (V-USC) for liquid water with classical molecular dynamics simulations. When the cavity modes are resonantly coupled to the O–H stretch mode of liquid water, the infrared spectrum shows asymmetric Rabi splitting. The lower polariton (LP) may be suppressed or enhanced relative to the upper polariton (UP) depending on the frequency of the cavity mode. Furthermore, although the static properties and the translational diffusion of water are not changed under VSC or V-USC, we do find the modification of the orientational autocorrelation function of H 2 O molecules especially under V-USC, which could play a role in ground-state chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗