Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Classical 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 OSTI2022

Interfacial water asymmetry at ideal electrochemical interfaces

Controlling electrochemical reactivity requires a detailed understanding of the charging behavior and thermodynamics of the electrochemical interface. Experiments can independently probe the overall charge response of the electrochemical double layer by capacitance measurements and the thermodynamics of the inner layer with potential of maximum entropy measurements. Relating these properties by computational modeling of the electrochemical interface has so far been challenging due to the low accuracy of classical molecular dynamics (MD) for capacitance and the limited time and length scales of ab initio MD. Here, we combine large ensembles of long-time-scale classical MD simulations with charge response from electronic density functional theory to predict the potential-dependent capacitance of a family of ideal aqueous electrochemical interfaces with different peak capacitances. We show that while the potential of maximum capacitance varies, this entire family exhibits an electrode charge of maximum capacitance (CMC) between -2.9 and -2.2 μC/cm 2 , regardless of the details in the electronic response. Simulated heating of the same interfaces reveals that the entropy peaks at a charge of maximum entropy (CME) of -5.1 ± 0.6 μC/cm 2 , in agreement with experimental findings for metallic electrodes. The CME and CMC both indicate asymmetric response of interfacial water that is stronger for negatively charged electrodes, while the difference between CME and CMC illustrates the richness in behavior of even the ideal electrochemical interface.

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

The effect of hydroxyl spacing in diols on the solvation structure, dynamics, and transport properties of choline chloride-based deep eutectic solvents

Deep eutectic solvents (DESs) are a class of liquids that offer great potential in alleviating some of the challenges present in today's long-term energy storage methods because they have physical properties that are favorable for storable electrolyte solutions. In this work, a series of glycols (ethylene glycol, 1,3-propanediol, 1,4-butanediol, and 1,5-pentanediol) were studied as potential hydrogen bond donors (HBD) with a common choline chloride (ChCl) as the hydrogen bond acceptor (HBA). The solvation dynamics of the prepared systems were studied by measuring the solvent reorganization response using femtosecond transient absorption spectroscopy (fs-TA). Conductivity, viscosity, density, ET(30) polarity, and dynamics of the prepared DESs were analyzed, with a particular interest in determining the effect of HBD chain length on these parameters. Here, classical molecular dynamics simulations were employed to investigate how the local liquid structure, solvent dynamics, and bulk solvent properties vary with changes in glycol chain length.

Conductivity↗
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 OSTI2025

Phonon second harmonic generation in NaBr studied by inelastic neutron scattering and computer simulation

The phenomenon of second harmonic generation (SHG) was found for phonons in anharmonic NaBr by inelastic neutron scattering. The temperature dependence of this phonon SHG was measured from 300 K to 650 K. At 300 K the second harmonic (SH) is seen as a high-energy branch around 33 meV, nearly independent of $\overrightarrow{Q}$. The temperature effective potential (TDEP) method and classical molecular dynamics (MD) simulation with machine learning interatomic potential were able to reproduce the SH, and showed that SHG occurs with the flat transverse optical (TO) phonon branch. A classical model of a nonlinear medium explains the intensity and lifetime of the SH, compared to those of the TO modes. Also successful was a quantum model based on the Heisenberg-Langevin equation for interacting phonons coupled to a thermal bath, which also predicts a spectral distribution of the SH. In conclusion, the measured temperature dependence of the intensity of the second harmonic showed that it follows the Planck distribution of a one-phonon quasiparticle, and not two TO phonons.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A molecular dynamics study of laser-excited gold

The structural evolution of laser-excited systems of gold has previously been measured through ultrafast MeV electron diffraction. However, there has been a long-standing inability of atomistic simulations to provide a consistent picture of the melting process, leading to large discrepancies between the predicted threshold energy density for complete melting, as well as the transition between heterogeneous and homogeneous melting. We make use of two-temperature classical molecular dynamics simulations utilizing three highly successful interatomic potentials and reproduce electron diffraction data presented by Mo et al. [Science 360, 1451–1455 (2018)]. We recreate the experimental electron diffraction data, employing both a constant and temperature-dependent electron–ion equilibration rate. In all cases, we are able to match time-resolved electron diffraction data, and find consistency between atomistic simulations and experiments, only by allowing laser energy to be transported away from the interaction region. This additional energy-loss pathway, which scales strongly with laser fluence, we attribute to hot electrons leaving the target on a timescale commensurate with melting.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Structure of ice confined in silica nanopores

Observed anomalous thermodynamic properties of confined water such as deviations in the melting point and freezing point motivate the determination of the structure of confined water as a function of pore size and temperature. Here, we investigate the dynamic evolution of the structure of confined ice in SBA-15 porous materials with pore diameters of 4 nm, 6 nm, and 8 nm at temperatures ranging from 183 K to 300 K using in operando Wide-Angle X-Ray Scattering (WAXS) measurements, X-Ray Partial Distribution Function (PDF) measurements, and classical Molecular Dynamics (MD) simulations. Formation of hexagonal ice structures is noted in all the three pore sizes. In silica nanopores with diameters of 4 nm, cubic ice formation is noted in addition to hexagonal ice. Longer lasting hydrogen bonds and longer residence times of the water molecules in the first coordination shell contribute to observed crystalline organization of ice in confinement. Self-diffusion coefficients of confined liquid water, predicted from classical MD simulations, are four orders of magnitude higher compared to ice formed in confinement. These experimental and simulation results provide comprehensive insights underlying the organization of confined water and ice in silica nanopores and the underlying physico-chemical interactions that contribute to the observed structures.

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

Bespoke Liquid/Liquid Interfaces (Final Technical Report)

The goal of DE-SC0001815 was to advance the basic science of liquid:liquid interface formation, to develop a deeper understanding of the mechanisms of phase separation and the essential relationships between solution composition, organization and dynamics that underlie the kinetic regime of solvent extraction. This included learning how interfacial organization and dynamics alters the properties of the primary coordination sphere of ions and the free energy of transport of ions complexes across a phase boundary. We relied primarily upon classical molecular dynamics studies to determine the equilibrium ensembles of these complex systems, but also utilized ab-initio MD and cluster-based density functional theory (DFT) calculations when more detailed investigation of the electronic structure was needed. We continued development of graph-theory based analyses to elucidate hierarchical correlations and expanded into geometric topology methods to quantify the collectively organized structures that can organize at a liquid/liquid interface during solute transport. One of the main conclusions was from the observation of two distinct mechanisms for solute transport - those that derive from amplifications of interfacial heterogeneity and surface roughness, and those wherein surface roughness has been dampened and instead collectively organized macrostructures work to bring solutes into the organic phase. It was our aim to create a concrete chemical model of the underlying driving forces behind interfacial primary and secondary structure formation and to map out the energetic features of solute transport so that tailored liquid/liquid can be developed that have characteristic kinetic features associated with mass transport.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Atomic scale etching of diamond: insights from molecular dynamics simulations

Diamond is a promising material for multiple applications in quantum information processing and sensing as well as applications in microelectronics. However, diamond devices can be limited by surface defects that compromise charge stability and spin coherence, among others. Improved strategies in plasma etching of diamond could play an important role in minimizing or eliminating these defects. In this work, we explore plasma-assisted atomic scale etching of diamond using argon ions (Ar + ), hydrogen ions (H + ) and hydrogen atoms (H). We employ classical molecular dynamics (MD) simulations and test several interatomic potentials based on the Reactive Empirical Bond Order (REBO) form with comparisons to a variety of published experimental results. We performed MD simulations of low-energy hydrogen ($\leqslant$50 eV) and argon ( $\leqslant$200 eV) ion bombardment of diamond surfaces. Ar + bombardment can be used to locally smooth initially rough diamond surfaces via the formation of an amorphous C layer, the thickness of which increases with argon ion energy. Subsequent exposure with hydrogen ions (or fast neutrals) will selectively etch this amorphous C layer, leaving the underlying diamond layer mostly intact if the H energy is maintained below about 10 eV. The simulations suggest that combining Ar + smoothing with selective, near threshold energy H removal of amorphous C can be an effective strategy for diamond surface engineering, leading to more reliable and sensitive diamond color center devices.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS2020

Ultrafast Melting at Metal-Semiconductor Interfaces: A Molecular Dynamics Study

Quantifying the thermodynamics of interfacial melting is important for several applications such as laser machining and designing radiation hardened electronic devices. Due to the short time and length scales that are involved in the process, molecular simulations provide an ideal path for understating the interfacial melting phenomena. In this work we investigate tungsten (W) – silicon carbide (SiC) interfaces which are relevant in modern high-power electronic devices. Classical molecular dynamics simulations along with the two-temperature model (TTM) were used to study the radiation induced melting of the interface while a neural network-based machine learning framework was used to track the transient re-crystallization process. The damage endured by the semiconductor at the interface following the radiation event is discussed in detail while the re-crystallization of the metal is described in terms of the classical nucleation theory.

Ashwin Ravichandran↗
DOE OSTI2024

Temperature and density dependent pair potential for deuterium under shock

Large-scale classical molecular dynamics (CMD) simulations naturally include the microscopic physics necessary for atomistic modeling of shock release at the ablator-fuel interface in an inertial confinement fusion (ICF) capsule. Here, the multi-megabar shocks utilized in ICF experiments can drive the deuterium fuel from ambient to electron volt temperatures (T) and multi-fold compression. Modeling interatomic interactions over such an extreme range of conditions is challenging for empirical bond order potentials. We generate a pair potential for deuterium with explicit temperature and mass density dependence from ab initio density functional theory molecular dynamics using the iterative Boltzmann inversion method. This potential accurately reproduces the radial distribution functions and pressures from DFT in CMD equilibrium simulations across a wide range of thermodynamic conditions, yet fails to return the expected Hugoniot relations when used in direct CMD shock simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Structural and Spectroscopic Properties of Butanediol-Modified Boehmite Materials

Glycoboehmite (GB) materials are synthesized by a solvothermal reaction to form layered aluminum oxyhydroxide (boehmite) modified by intercalated butanediol molecules. These hybrid materials offer a platform to design materials with potentially novel sorption, wetting, and catalytic properties. Several synthetic methods have been used, resulting in different structural and spectroscopic properties, but atomistic detail is needed to determine the interlayer structure to explore the synthetic control of GB materials. Here, in this study, we use classical molecular dynamics (MD) simulations to compare the structural properties of GB interlayers containing chemisorbed butanediol molecules as a function of diol loading. Accompanying quantum (density functional theory, DFT) static calculations and MD simulations are used to validate the classical model and compute the infrared spectra of various models. Classical MD results reveal the existence of two unique interlayer environments at higher butanediol loading, corresponding to smaller (cross-linked) and expanded interlayers. DFT-computed infrared spectra reveal the sensitivity of the aluminol O–H stretch frequencies to the interlayer environment, consistent with the spectrum of the synthesized material. Insight from these simulations will aid in the characterization of the newly synthesized GB materials.

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

Elucidating molecular scale interactions underlying the freezing behavior of salt solutions in silica nanopores

Elucidating the influence of nanoscale confinement on the freezing behavior of salt solutions is of fundamental interest to environmental security and materials science and engineering. Specific structural information such as the coordination environment of ions in confined aqueous medium, effects on the extended hexagonal network of water molecules and their mutual non-bonding interactions in confinement are sparse in the literature along with the change in the dynamical characteristics of water in the presence of ions. To address these knowledge gaps, the current study is focused on investigating the influence of reduced dimensionality arising from nanoscale confinement on the structural evolution of salt solutions on freezing and the associated fluid–surface interactions. In this regard, operando wide-angle X-ray scattering (WAXS) measurements and classical molecular dynamics (MD) simulations are conducted with water and 0.5 M CaCl 2 , MgCl 2 and KCl solutions confined in 4 nm sized SBA-15 silica pores upon cooling from 300 K to 200 K. The freezing point of the salt solutions is depressed to 235 K which is about 10 K lower than that of confined water. The translational dynamics of confined salt solutions indicates shifts in the fragile-to-strong dynamical crossover to a lower temperature compared to pure water. The strong electrostatic attraction between the cations and the surrounding water molecules contributes to the freezing point depression in confined salt solutions. These insights unlock the molecular-scale basis and mechanisms underlying the freezing behavior of confined salt solutions.

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

Tracking structural solvent reorganization and recombination dynamics following e – photoabstraction from aqueous I – with femtosecond x-ray spectroscopy and scattering

Here we present a sub-picosecond resolved investigation of the structural solvent reorganization and geminate recombination dynamics following 400 nm two-photon excitation and photodetachment of a valence p electron from the aqueous atomic solute, I–(aq). The measurements utilized time-resolved X-ray Absorption Near Edge Structure (TR-XANES) spectroscopy and X-ray Solution Scattering (TR-XSS) at the Linac Coherent Light Source x-ray free electron laser in a laser pump/x-ray probe experiment. The XANES measurements around the L 1 -edge of the generated nascent iodine atoms (I 0 ) yield an average electron ejection distance from the iodine parent of 7.4 ± 1.5 Å with an excitation yield of about 1/3 of the 0.1M NaI aqueous solution. The kinetic traces of the XANES measurement are in agreement with a purely diffusion-driven geminate iodine–electron recombination model without the need for a long-lived (I 0 :e–) contact pair. Nonequilibrium classical molecular dynamics simulations indicate a delayed response of the caging H 2 O solvent shell and this is supported by the structural analysis of the XSS data: We identify a two-step process exhibiting a 0.1 ps delayed solvent shell reorganization time within the tight H-bond network and a 0.3 ps time constant for the mean iodine–oxygen distance changes. The results indicate that most of the reorganization can be explained classically by a transition from a hydrophilic cavity with a well-ordered first solvation shell (hydrogens pointing toward I–) to an expanded cavity around I 0 with a more random orientation of the H 2 O molecules in a broadened first solvation shell.

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

Better, Cheaper, Faster Molecular Dynamics

Recent, revolutionary progress in genomics and structural, molecular and cellular biology has created new opportunities for molecular-level computer simulations of biological systems by providing vast amounts of data that require interpretation. These opportunities are further enhanced by the increasing availability of massively parallel computers. For many problems, the method of choice is classical molecular dynamics (iterative solving of Newton's equations of motion). It focuses on two main objectives. One is to calculate the relative stability of different states of the system. A typical problem that has' such an objective is computer-aided drug design. Another common objective is to describe evolution of the system towards a low energy (possibly the global minimum energy), "native" state. Perhaps the best example of such a problem is protein folding. Both types of problems share the same difficulty. Often, different states of the system are separated by high energy barriers, which implies that transitions between these states are rare events. This, in turn, can greatly impede exploration of phase space. In some instances this can lead to "quasi non-ergodicity", whereby a part of phase space is inaccessible on time scales of the simulation. To overcome this difficulty and to extend molecular dynamics to "biological" time scales (millisecond or longer) new physical formulations and new algorithmic developments are required. To be efficient they should account for natural limitations of multi-processor computer architecture. I will present work along these lines done in my group. In particular, I will focus on a new approach to calculating the free energies (stability) of different states and to overcoming "the curse of rare events". I will also discuss algorithmic improvements to multiple time step methods and to the treatment of slowly decaying, log-ranged, electrostatic effects.

Pohorille, Andrew↗
DOE OSTI2023

Quantum bath augmented stochastic nonequilibrium atomistic simulations for molecular heat conduction

Classical molecular dynamics (MD) has been shown to be effective in simulating heat conduction in certain molecular junctions since it inherently takes into account some essential methodological components which are lacking in the quantum Landauer-type transport model, such as many-body full force-field interactions, anharmonicity effects and nonlinear responses for large temperature biases. However, the classical MD reaches its limit in the environments where the quantum effects are significant (e.g. with low-temperatures substrates, presence of extremely high frequency molecular modes). Here, we present an atomistic simulation methodology for molecular heat conduction that incorporates the quantum Bose–Einstein statistics into an “effective temperature” in the form of a modified Langevin equation. We show that the results from such a quasi-classical effective temperature MD method deviates drastically when the baths temperature approaches zero from classical MD simulations and the results converge to the classical ones when the bath approaches the high-temperature limit, which makes the method suitable for full temperature range. In addition, we show that our quasi-classical thermal transport method can be used to model the conducting substrate layout and molecular composition (e.g. anharmonicities, high-frequency modes). Anharmonic models are explicitly simulated via the Morse potential and compared to pure harmonic interactions to show the effects of anharmonicities under quantum colored bath setups. Finally, the chain length dependence of heat conduction is examined for one-dimensional polymer chains placed in between quantum augmented baths.

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

Charged Particle Dynamics in the Lunar Environment

As the Artemis Program streams forward, organizations of scientists and engineers across the country have been coming together to solve the complex network of problems to once again achieve the milestone of successfully touching down on the moon. This project is no exception, and it has been an honor to work with the Electrostatics and Surface Physics Laboratory (ESPL), a lab within the Exploration Research and Technology Programs’ Spaceport Technologies Office (UB-G) located at the National Aeronautics and Space Administration at Kennedy Space Center (NASA KSC). The authors and their mentor James R. Phillips III3, in tandem with researchers from the Astrodynamics and Space Robotics Laboratory (ASRL) at the University of Central Florida (UCF), have been working on creating a state-of-the-art (SOA) granular gas dynamics model for particulate contamination prevention and mitigation purposes. In essence, the underlying objective for this project is to more accurately model the resulting electrodynamic interactions between lunar regolith grains with applications to dust mitigation and rocket engine plume surface interactions. To achieve this goal, the team has been expanding upon existing open source classical molecular dynamics code developed by Sandia National Laboratories (SNL). The following reports on the details of the problem at hand as well as the contributions that the authors have made towards resolution, including but not limited to the encoding of physical attributes and interactions for non-spherical polydisperse particle distributions within various bed geometries and preemptive data analysis implementations. Significant data analysis processes were utilized, and several original algorithms were created to perform critical evaluations, resulting in only a 0.006% error in discrepancy.

Annelisa B Esparza↗
DOE OSTI2023

Molecular Dynamics Simulation of Hygroscopic Aging Effects in Epoxy Polymer

The automobile industry is incorporating more lightweight content in car designs to boost fuel-economy. New structural adhesives are needed to mitigate the corrosion and thermal expansion issues associated with joining dissimilar lightweight materials, but adhesive developers lack a fundamental understanding of the chemistry that occurs in the adhesive as the joint ages. In this study, we developed structural adhesive molecular models and applied classical molecular dynamics simulations and density functional theory calculations to gain molecular insights into the influence of water molecules on the properties of epoxy-based adhesives (DGEBA + Jeffamine (JD230)). The simulations were complemented by experimental synthesis and characterization. Our work underscores the impact of water molecules on the local structure of the epoxy network as well as resulting mechanical properties. Water molecules were mainly coordinated with hydroxyls, primary amines and secondary amines, but also weakly coordinated with ether linkages, which were found most probable to be labile. Simulated stress–strain data indicates that increasing the water content deteriorates the mechanical properties. The Young’s modulus decreased by ~ 30% when the water content increased to 3 wt%. We conclude, this integration of molecular-level chemical insights with mechanical property simulations of the hydrated epoxy system and experimental validation holds the promise to advance lightweight joint technologies.

36 MATERIALS SCIENCE↗
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↗