Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Classical molecular dynamic simulations”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Damage initiation and evolution in Al-Si layered microstructures under shock loading conditions at atomic scales

Designing materials with microstructural features like multiphase interfaces have shown significant promise for usage in the next generation defense and nuclear applications. The dynamic response of these interfaces in extreme environments is observed to vary with the deformability of individual phases, which could subsequently alter the favored damage nucleation sites related to spall failure. Here, we investigate the role of spacing of interfaces in a nanocrystalline layered Aluminum-Silicon system on the shock wave propagation behavior, and microstructural and defect evolution using classical molecular dynamics simulations. The simulations are carried out with different hypothetical Al/Si microstructures including variations in the distribution of Si grains as layers in a nanocrystalline Al matrix. The molecular dynamics simulations suggest that spall failure is preferentially initiated at Al/Si interfaces. In addition, for the same system volume and same concentration of Si, a higher number of layers is marked by an increase in shock wave velocity and reduced resistance to spall failure.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Optimized utilization of COMB3 reactive potentials in LAMMPS

An investigation to optimize the application of the third-generation charge optimized many-body (COMB3) interatomic potential and associated input parameters was carried out through the study of solid–liquid interactions in classical molecular dynamics simulations. The rates of these molecular interactions are understood through the wetting rates of water nano-droplets on a bare copper (111) surface. Implementing the Langevin thermostat, the influence of simulation time step, the number of atoms in the system, the frequency at which charge equilibration is performed, and the temperature relaxation rate are all examined. The results indicate that time steps of 0.4 fs are possible when using longer relaxation times for the system temperature, which is almost double the typical time step used for reactive potentials. Additionally, the use of the charge equilibration allows for a fewer atomic layers to be used in the Cu slab. In addition, charge equilibrium schemes do not need to be performed every time step to ensure accurate charge transfer. Interestingly, the rate of wetting for the nanodroplets is dominantly dependent on the temperature relaxation time, which is predicted to significantly change the viscosity of the water droplets. This work provides a pathway for optimizing simulations using the COMB3 reactive interatomic potential.

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

Vapor–liquid equilibrium of water with the MB-pol many-body potential

Among the many existing molecular models of water, the MB-pol many-body potential has emerged as a remarkably accurate model, capable of reproducing thermodynamic, structural, and dynamic properties across water's solid, liquid, and vapor phases. In this work, we assessed the performance of MB-pol with respect to an important set of properties related to vapor-liquid coexistence and interfacial behavior. Through direct coexistence classical molecular dynamics simulations at temperatures of 400 K < T < 600 K, we calculated properties such as equilibrium coexistence densities, vapor-liquid interfacial tension, vapor pressure, and enthalpy of vaporization and compared the MB-pol results to experimental data. We also compared rigid vs fully flexible variants of the MB-pol model and evaluated system size effects for the properties studied. We found that the MB-pol model predictions are in good agreement with experimental data, even for temperatures approaching the vapor-liquid critical point; this agreement was largely insensitive to system sizes or the rigid vs flexible treatment of the intramolecular degrees of freedom. Furthermore, these results attest to the chemical accuracy of MB-pol and its high degree of transferability, thus enabling MB-pol's application across a large swath of water's phase diagram.

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

Energy, angle, and temperature dependencies of the sticking of D atoms on Li surfaces

Detailed experimental and computational information on the response of lithium surfaces to irradiation by slow hydrogenic particles (ions, atoms, molecules) is sparse and mainly speculative. In this work, we present a computational study of the reflection and retention of deuterium (D) atoms at crystalline and amorphous lithium surfaces at 300 and 500 K, where the D atoms have an impact energy in the range of 0.025–5 eV and incident angles of 0° (perpendicular incidence) or 85° (near-grazing incidence). Classical molecular dynamics simulations are performed with the reactive bond-order force field (ReaxFF) potentials. This study provides quantitative information on the deuterium sticking probability and recycling coefficient for lithium surfaces. Further, our results support the ongoing work at the Lithium Tokamak eXperiment-β fusion experiment as well as relevant experiments in the laboratory setting.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

DEHP− extractant binding to trivalent lanthanide Er3+: Fast binding accompanied by concerted angular motions of hydration water

Solvent extraction of trivalent rare earth metal ions by organophosphorus extractants proceeds via binding of phosphoric acid headgroups to the metal ion. Water molecules in the tightly bound first hydration shell of the metal ions must be displaced by oxygen atoms from phosphoric acid headgroups. Here, we use classical molecular dynamics simulations to explore the event in which a fully hydrated Er3+ binds to its first phosphoric acid headgroup. Approach of the headgroup into the region between the first and second hydration shells leads to a fast ejection of a water molecule that is accompanied by reordering of the hydration water molecules, including discretization of their angular positions and collective rotation about the metal ion. The water molecule ejected from the first shell is located diametrically opposite from the binding oxygen. Headgroup binding places a headgroup oxygen closer to Er3+ than its first hydration shell and creates a loosely bound water that subsequently exchanges between the first shell and its environment. This second exchange of water also occurs at discrete angular positions. This geometrical aspect of binding may be of relevance to understanding the binding and transport of ion–extractant complexes that are expected to occur at the organic–aqueous liquid–liquid interface used in solvent extraction processes.

Chemistry↗
DOE OSTI2023

Molecular dynamics study on the interfacial properties of mixtures of monomers of polyvinylpyrrolidone (PVP)-based battery binders on graphene and graphite surfaces

This study investigates the behavior of two different mixtures of monomers of polyvinylpyrrolidone (PVP)-based battery binders, polyvinylpyrrolidone:polyvinylidene difluoride (PVP:PVDF) and polyvinylpyrrolidone:polyacrylic acid (PVP:PAA), at graphene and graphite interfaces using classical molecular dynamics simulations. The aim is to identify the best performing monomer binder blend and carbon-based material for the design of battery-optimized energy devices. The PVP:PAA monomer binder blend and graphite are found to have the best interaction energies, densification upon adsorption, and more ordered structure. The adsorption of both monomer binder blends is strongly guided by the higher affinity of PVP and PAA monomeric molecules for the surfaces compared to PVDF. The structure of adsorbed layers of PVP:PVDF monomer binder blend on graphene and graphite develops more quickly than PVP:PAA, indicating faster kinetics. This study complements a previous density functional theory study recently reported by our group and contributes to a better understanding of the nanoscopic features of relevant interfacial regions involving mixtures of monomers of PVP-based battery binders and different carbon-based materials. In conclusion, the effect of a blend of commonly used monomer binders on carbon-based materials is essential for obtaining tightly bound anode and cathode active materials in lithium-ion batteries, which is crucial for designing battery-optimized energy devices.

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

What do far-infrared spectra of solitary water in “water-in-solvent” systems reveal about water’s solvation and dynamics?

Classical molecular dynamics simulations of water in ionic and dipolar solvents were used to interpret the far-infrared (FIR) rotation/libration spectra of “solitary water” in terms of water’s rotational dynamics and interactions with solvents. Seven solvents represented by nonpolarizable all-atom force fields and a series of idealized variable-charge solvents were used to span the range of solvent polarities (hydrogen bonding) studied experimentally. Simulated spectra capture the solvent dependence observed, as well as the relationship between the frequencies of water libration (νL) and OH stretching bands (νOH). In more strongly interacting solvents, simulated νL are ∼20% higher than those of experiment. In all solvents, the simulated spectra are composites of rotational motions about the two axes perpendicular to water’s dipole moment, and the different frequencies of these two motions are responsible for the breadth of the libration band and the bimodal shape observed in halide ionic liquids. Simulations overestimate the separation of these two components in most solvents. The character of water rotational motions changes markedly with solvent polarity, from quasi-free rotation in nonpolar and weakly polar solvents to highly constrained libration in strongly hydrogen bonding environments. The changeover to librational motions dominating the spectrum occurs between solvents such as benzene (νL ∼ 250 cm−1) and acetonitrile (νL ∼ 400 cm−1). For solvents in the latter category, the mean frequency of the experimental FIR band provides a direct measure of mean-squared torques and, therefore, force constants associated with interactions constraining water’s librational motion.

Chemistry↗
DOE OSTI2024

Melting point of iron at high pressure: An assessment of uncertainties and effect of electronic temperature

An accurate calculation of the melting point of iron at various pressures in the Earth's core is important for understanding the core structure, geodynamo, and the Earth's history. Previous studies have assessed the melt line of iron at these extreme conditions using various experimental measurement techniques as well as both ab initio and classic molecular dynamics simulations. However, experimental measurements have uncertainties up to several hundred Kelvin, and inconsistencies remain among simulation results. Here in this work, we propose an iterative framework that couples density functional theory (DFT) calculations and molecular dynamics simulations performed using an ensemble of interatomic potentials to assess the effect of electronic temperature on the melting point. We systematically validate the potentials by comparing lattice constants and phonon dispersion curves at 0 K and enthalpy differences between liquid and HCP, FCC, BCC phases of iron close to the melt line at 300 GPa with DFT. Our results show that HCP iron melts at 6144 K (at 300 GPa), BCC phase is thermodynamically unstable, and FCC is metastable at this temperature. The melting points of FCC and BCC phases at 300 GPa are 5858 and 5647 K, respectively.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Mean force emission theory for classical bremsstrahlung in strongly coupled plasmas

This work presents mean force emission theory, which extends the classical theory of bremsstrahlung emission to strongly coupled plasmas. In the high-frequency limit, the theory reduces to solving for the electron trajectory during a binary collision, but where the electron–ion interactions occur through the potential of mean force. In the low-frequency limit, it uses an autocorrelation formalism that captures effects of multiple collisions and strongly correlated motion. The predictions are benchmarked by comparison with first-principles classical molecular dynamics simulations of a fully ionized hydrogen plasma in which all interactions are repulsive. The comparison shows good agreement up to Coulomb coupling strengths of Γ∼30. The theory improves upon traditional models by including strong coupling effects and systematically including the effect of multiple collisions. Furthermore, mean force emission theory provides evidence that the Drude correction factor commonly used in quantum calculations of optical quantities may not be adequate at strong coupling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Solving Newton’s equations of motion with large timesteps using recurrent neural networks based operators

Classical molecular dynamics simulations are based on solving Newton’s equations of motion. Using a small timestep, numerical integrators such as Verlet generate trajectories of particles as solutions to Newton’s equations. We introduce operators derived using recurrent neural networks that accurately solve Newton’s equations utilizing sequences of past trajectory data, and produce energy-conserving dynamics of particles using timesteps up to 4000 times larger compared to the Verlet timestep. We demonstrate significant speedup in many example problems including 3D systems of up to 16 particles.

Newton’s equations↗
DOE OSTI2026

Stoichiometry dependent properties of cerium hydride: An active learning developed interatomic potential study

Cerium hydride has a variety of interesting properties, including a known lattice contraction and densification with increasing hydrogen content. However, precise stoichiometric control is not experimentally straightforward and ab initio approaches are not computationally feasible for many properties such as melting and low temperature diffusion. Therefore, we develop a machine-learned interatomic potential for cerium hydride that is valid for H to Ce ratios from 2.0 to 3.0. A query-by-committee active learning approach is used to develop the training set. Leveraging classical molecular dynamics simulations, we assess a range of properties and provide fundamental mechanisms for the trends with stoichiometry. Finally, a majority of the properties follow the trend of lattice contraction, being governed by the stronger lattice binding induced by adding octahedral atoms.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Diffuse scattering from dynamically compressed single-crystal zirconium following the pressure-induced α → ω phase transition

The prototypical α → ω phase transition in zirconium is an ideal test bed for our understanding of polymorphism under extreme loading conditions. After half a century of study, a consensus had emerged that the transition is realized via one of two distinct displacive mechanisms, depending on the nature of the compression path. However, recent dynamic-compression experiments equipped with diffraction diagnostics performed in the past few years have revealed new transition mechanisms, demonstrating that our understanding of the underlying atomistic dynamics and transition kinetics is in fact far from complete. We present classical molecular dynamics simulations of the α → ω phase transition in single-crystal zirconium shock compressed along the [0001] axis using a machine-learning-class potential. The transition is predicted to proceed primarily via a modified version of the two-stage Usikov-Zilberstein mechanism, whereby the high-pressure ω phase heterogeneously nucleates at boundaries between grains of an intermediate β phase. We further observe the fomentation of atomistic disorder at the junctions between β grains, leading to the formation of highly defective interstitial material between the ω grains. We directly compare synthetic x-ray diffraction patterns generated from our simulations with those obtained using femtosecond diffraction in recent dynamic-compression experiments, and show that the simulations produce the same unique, anisotropic diffuse scattering signal unlike any previously seen from an elemental metal. Our simulations suggest that the diffuse signal arises from a combination of thermal diffuse scattering, nanoparticlelike scattering from residual kinetically stabilized α and β grains, and scattering from interstitial defective structures. Published by the American Physical Society 2024

36 MATERIALS SCIENCE↗
DOE OSTI2022

Discovery of multimechanisms of screw dislocation interaction in bcc iron from open-ended saddle point searches

Dislocation motion and interactions determine mechanical properties in body-centered cubic (bcc) metallic materials. However, studying mechanisms for the screw dislocation interaction is fundamentally challenging since many underlying processes involve mesotimescales and atomistic resolution, currently inaccessible by either experimental techniques or continuum theoretical methods. Here, in this paper, we develop a computational capability based on self-evolving atomistic kinetic Monte Carlo (SEAKMC) to sample the critical events and saddle point energies related to screw dislocations and their junctions. The method is first validated by calculating the stress dependence of Peierls barriers and formation energies of kink pairs and cross-slip kink pairs on a single screw dislocation in bcc iron. Then the method is applied to a binary junction of a pair of intersecting screw dislocations, the structure of which is crucial for low-temperature plastic deformation. We identify three important mechanisms: coplanar cross-slipping, jog-pinning, and a previously unknown unzipping mechanism during the evolution of the binary junction. The mechanisms are then further validated using classical molecular dynamics simulations. The computational capability developed in this paper provides an effective tool to evaluate screw dislocation related thermally activated events in complex stress conditions. The mechanisms discovered in this paper provide critical insights into temperature dependence of the anomalous slip, a breakdown of the Schmidt law, during the plastic deformation in bcc iron and can be generalized to other bcc metals.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Fracture mechanisms of sodium silicate glasses

Reactive classical molecular dynamics simulations of sodium silicate glasses, xNa 2 O–(100 – x)SiO 2 (x = 10–30), under quasi-static loading, were performed for the analysis of molecular scale fracture mechanisms. Mechanical properties of the sodium silicate glasses were consistent with experimentally reported values, and the amount of crack propagation varied with reported fracture toughness values. The most crack propagation occurred in NS20 systems (20-mol% Na 2 O) compared with the other simulated compositions. Dissipation via two mechanisms, the first through sodium migration as a lower activation energy process and the second through structural rearrangement as a higher activation energy process, was calculated and accounted for the energy that was not stored elastically or associated with the formation of new fracture surfaces. Here, a correlation between crack propagation and energy dissipation was identified, with systems with higher crack propagation exhibiting less energy dissipation. Sodium silicate glass compositions with lower energy dissipation also exhibited the most sodium movement and structural rearrangement within 10 Å of the crack tip during loading. Therefore, high sodium mobility near the crack tip may enable energy dissipation without requiring formation of structural defects. Therefore, the varying mobilities of the network modifiers near crack tips influence the brittleness and the crack growth rate of modified amorphous oxide systems.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Inelastic relaxation processes in amorphous sodium silicates

During fracture amorphous oxides exhibit irreversible processes, including inelastic and nonrecoverable relaxation effects in the process zone surrounding the crack tip. Here, classical molecular dynamics simulations were used with a reactive forcefield to evaluate inelastic relaxation processes in five amorphous sodium silicate compositions. Overall, the 20% Na 2 O-SiO 2 (NS20) composition exhibited the most inelastic relaxation, followed by the 15% Na 2 O-SiO 2 (NS15) composition, the 25% Na 2 O-SiO 2 (NS25) composition, and finally the 10% (NS10) and 30% (NS30) Na 2 O-SiO 2 compositions. Coordination analysis of the Na + ions identified that during inelastic relaxation the Na+ ions were increasingly coordinated by nonbridging oxygens (NBOs) for the NS10 and NS15 compositions, which was supported by radial analysis of the O-Na-O bond angles surrounding the crack tip. Across the sodium silicate compositional range, two different inelastic relaxation mechanism were identified based on the amount of bridging oxygens (BOs) and NBOs in the Na+ ion coordination shell. At lower (NS10) and higher (NS30) sodium compositions, the entire structured relaxed toward the crack tip. In contrast at intermediate sodium concentrations (NS20) the Na + ion migrates toward the crack tip separately from the network structure. By developing a fundamental understanding of how modified silica systems respond to static stress fields, we will be able to predict how varying amorphous silicate systems exhibit slow crack growth.

36 MATERIALS SCIENCE↗
DOE OSTI2025

On the role of α-alumina in the origin of life: Surface-driven assembly of amino acids

We investigate the hypothesis that mineral/water interfaces played a crucial catalytic role in peptide formation by promoting the self-assembly of amino acids. Using classical molecular dynamics simulations, we demonstrate that the α-alumina(0001) surface exhibits an affinity of 4k B Tfor individual glycine or GG dipeptide molecules due to hydrogen bonds. In simulations with multiple glycine molecules, surface-bound glycine enhances further adsorption, leading to the formation of long chains connected by hydrogen bonds between the carboxyl and amine groups of glycine molecules. We find that the likelihood of observing chains longer than 10 glycine units increases by at least five orders of magnitude at the surface compared to the bulk. This surface-driven assembly is primarily due to local high density and alignment with the alumina surface pattern. Together, these results propose a model for how mineral surfaces can induce configuration-specific assembly of amino acids, thereby promoting condensation reactions.

Science & Technology - Other Topics↗
DOE OSTI2021

A Theoretical Model for Computing Freezing Point Depression of Lithium-Ion Battery Electrolytes

Reliable prediction of freezing point depression in liquid electrolytes will accelerate the development of improved Li-ion batteries which can operate in low temperature environments. In this work we establish a computational methodology to calculate activity coefficients and liquidus lines for battery-relevant liquid electrolytes. Electronic structure methods are used in conjuction with classical molecular dynamics simulations and theoretical expressions for Born solvation energy, ion-atmosphere effects from Debye-Hückel theory and solvent entropic effects. The framework uses no a priori knowledge beyond neat solvent properties and the concentration of salt. LiPF 6 in propylene carbonate (PC), LiPF 6 in dimethyl carbonate (DMC) and LiClO 4 in DMC are investigated up to 1 molal with accuracy better than 3 °C when compared to experimental freezing point measurements. Finally, we find that the difference in freezing point depression between the propylene carbonate-based electrolyte and the dimethyl carbonate electrolytes originates from the difference in the solvent dielectric constant.

25 ENERGY STORAGE↗
DOE OSTI2024

DRIFT: Diffusivity Regulation of Isotopes by nanoconFinementT

Spontaneous isotope fractionation has been reported under nanoconfinement conditions in naturally occurring systems, but the origin of this phenomena is currently unknown. Two existing hypotheses have been proposed, one based on changes in the solvation environment of the isotopes that reduces the non-mass dependent hydrodynamics contribution to diffusion. The other is that isotopes have mass-dependent surface adsorption, varying their total diffusion through nanoconfined channels. To investigate these hypotheses, benchtop experiments, nuclear magnetic resonance (NMR) spectroscopy, and molecule scale modeling were applied. Classical molecular dynamics simulations identified that the Na + and Cl - hydration shells across the three different salt solutions ( 22 Na 35 Cl, 23 Na 35 Cl, 24 Na 35 Cl) did not vary as a function of the Na + isotope, but that there was a significant pore size effect, with larger hydration shells at larger pore sizes. Additionally, while total adsorption times did not vary as a function of the Na + isotope or pore size, the free ion concentration, or those adsorbed on the surface for <5% of the simulation time did exhibit isotope dependence. Experimentally, challenges occurred developing a repeatable experiment, but NMR characterization of water diffusion rates through ordered alumina membranes was able to identify the existence of two distinct water environments associated with water inside and outside the pore. Further NMR studies could be used to confirm variation in hydration shells and diffusion rates of dissolved ions in water. Ultimately, mass-dependence adsorption is a primary driver of variations in isotope diffusion rates, rather than variation in hydration shells that occur under nanoconfinement.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗