Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Molecular dynamics simulation”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Machine learning predictions of diffusion in bulk and confined ionic liquids using simple descriptors

Ionic liquids have many intriguing properties and widespread applications such as separations and energy storage. However, ionic liquids are complex fluids and predicting their behavior is difficult, particularly in confined environments. We introduce fast and computationally efficient machine learning (ML) models that can predict diffusion coefficients and ionic conductivity of bulk and nanoconfined ionic liquids over a wide temperature range (350–500 K). The ML models are trained on molecular dynamics simulation data for 29 unique ionic liquids as bulk fluids and confined in graphite slit pores. This model is based on simple physical descriptors of the cations and anions such as molecular weight and surface area. Here, we also demonstrate that accurate results can be obtained using only descriptors derived from SMILES (simplified molecular-input line-entry system) codes for the ions with minimal computational effort. This offers a fast and efficient method for estimating diffusion and conductivity of nanoconfined ionic liquids at various temperatures without the need for expensive molecular dynamics simulations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Direct observation of coherent femtosecond solvent reorganization coupled to intramolecular electron transfer

It is well known that the solvent plays a critical role in ultrafast electron-transfer reactions. However, solvent reorganization occurs on multiple length scales, and selectively measuring short-range solute–solvent interactions at the atomic level with femtosecond time resolution remains a challenge. Here we report femtosecond X-ray scattering and emission measurements following photoinduced charge-transfer excitation in a mixed-valence bimetallic (Fe ii Ru iii ) complex in water, and their interpretation using non-equilibrium molecular dynamics simulations. Combined experimental and computational analysis reveals that the charge-transfer excited state has a lifetime of 62 fs and that coherent translational motions of the first solvation shell are coupled to the back electron transfer. Overall, our molecular dynamics simulations identify that the observed coherent translational motions arise from hydrogen bonding changes between the solute and nearby water molecules upon photoexcitation, and have an amplitude of tenths of ångströms, 120–200 cm -1 frequency and ~100 fs relaxation time. This study provides an atomistic view of coherent solvent reorganization mediating ultrafast intramolecular electron transfer.

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

Iodine recombination in xenon solvent: Clusters in the gas to liquid-like state transition

Supercritical fluids (SCFs) have attracted significant attention as solvents for chemical reactions due to their unique properties, such as high diffusivity, low viscosity, and tunable solvation properties. These properties profoundly influence reaction kinetics and are often attributed to the formation of molecular clusters within SCFs. To study the effect of supercritical solvent on chemical reactivity and dynamics of reactions, one needs to understand the dynamics of clusters in supercritical fluid. Extensive experiments on the photodissociation and recombination of iodine in supercritical fluids served as a model system for understanding these effects. Experimental studies have been complemented by theoretical and computational investigations, which mostly employ Monte Carlo or empirical molecular dynamics simulations. However, computational studies using non-reactive force fields and ab initio approaches present challenges in capturing reactive processes at larger scales within supercritical fluids. Here, in this work, we developed the ReaxFF parameters by training against quantum mechanics data. ReaxFF reactive force field based molecular dynamics simulations were performed, studying the dynamics of a xenon solvent and cage effect at different thermodynamic conditions for the iodine recombination reaction. We show that the conditions near the critical point are the optimal conditions to study the cage effect. We show that the average lifetime of xenon clusters ranging between 5 and 11 ps is comparable to iodine geminate recombination. Our simulation results of iodine recombination in xenon solvent demonstrate the higher probability of iodine molecule formation in the presence of xenon clusters. Finally, we show that the supercritical condition exhibits the highest recombination rate for iodine atoms.

Cage effect↗
DOE OSTI2024

Probing the active sites of oxide encapsulated electrocatalysts with controllable oxygen evolution selectivity

Electrocatalysts encapsulated by nanoscopic overlayers can control the rate of redox reactions at the outer surface of the overlayer or at the buried interface between the overlayer and the active catalyst, leading to complex behavior in the presence of two competing electrochemical reactions. This study investigated oxide encapsulated electrocatalysts (OECs) comprised of iridium (Ir) thin films coated with an ultrathin (2–10 nm thick) silicon oxide (SiO x ) or titanium oxide (TiO x ) overlayer. The performance of SiO x |Ir and TiO x |Ir thin film electrodes towards the oxygen evolution reaction (OER) and Fe(II)/Fe(III) redox reactions were evaluated. An improvement in selectivity towards the OER was observed for all OECs. Overlayer properties, namely ionic and electronic conductivity, were assessed using a combination of electroanalytical methods and molecular dynamics simulations. SiO x and TiO x overlayers were found to be permeable to H 2 O and O 2 such that the OER can occur at the MO x |Ir (M = Ti, Si) buried interface, which was further supported with molecular dynamics simulations of model SiO 2 coatings. In contrast, Fe(II)/Fe(III) redox reactions occur to the same degree with TiO x overlayers having thicknesses less than 4 nm as bare electrocatalyst, while SiO x overlayers inhibit redox reactions at all thicknesses. This observation is attributed to differences in electronic transport between the buried interface and outer overlayer surface, as measured with through-plane conductivity measurements of wetted overlayer materials. These findings reveal the influence of oxide overlayer properties on the activity and selectivity of OECs and suggest opportunities to tune these properties for a wide range of electrochemical reactions.

08 HYDROGEN↗
DOE OSTI2025

Development and transferability of neural-network models for plasma-surface interactions

Plasma-surface interactions are increasingly critical to modern technologies; yet, accurate molecular dynamics simulations remain limited by the capabilities of interatomic potentials. Deep Potentials (DPs) promise to revolutionize the field by providing a systematic method for producing accurate interatomic potentials. The primary challenge of DP development is selecting a dataset, which efficiently spans the set of atomic environments one expects to encounter in the subsequent molecular dynamics simulations. The computational cost of density functional theory calculations, which are the typical basis for DP development, makes it impossible to directly verify the quality of a given DP. To address this challenge, we explore the development of a deep-learned interatomic potential, “DeepREBO,” trained to reproduce the behavior of the REBO2 empirical potential, enabling direct validation of training methodology and transferability. Using an active learning framework, we begin with a minimal dataset and iteratively expand it to train a Deep Potential-Smooth Edition model that faithfully reproduces REBO2 results for 25 eV hydrogen bombardment of diamond (001), a particularly challenging case. We show that small, carefully curated datasets can outperform large, unguided ones, with effective models requiring fewer than 15 000 snapshots. Subsequent transferability tests demonstrate that while DeepREBO generalizes well to diamond (111) surfaces, performance degrades for amorphous carbon or higher-energy impacts, highlighting the need for use-case-specific training data. We also evaluate methods to improve short-range repulsion. This study outlines best practices for training robust deep potentials and underscores the importance of dataset design for predictive plasma simulations.

Ab-initio molecular dynamics↗
DOE OSTI2021

A kinetic model of friction in strongly coupled strongly magnetized plasmas

Plasmas that are strongly magnetized in the sense that the gyrofrequency exceeds the plasma frequency exhibit novel transport properties that are not well understood. As a representative example, we compute the friction force acting on a massive test charge moving through a strongly coupled and strongly magnetized one-component plasma using a generalized Boltzmann kinetic theory. Recent works studying the weakly coupled regime have shown that strong magnetization leads to a transverse component of the friction force that is perpendicular to both the Lorentz force and velocity of the test charge, in addition to the stopping power component aligned antiparallel to the velocity. Recent molecular dynamics simulations have also shown that strong Coulomb coupling in addition to strong magnetization gives rise to a third component of the friction force in the direction of the Lorentz force. Here, we show that the generalized Boltzmann kinetic theory captures these effects and generally agrees well with the molecular dynamics simulations over a broad range of Coulomb coupling and magnetization strength regimes. The theory is also used to show that the “gyro” component of the friction in the direction of the Lorentz force arises due to asymmetries associated with gyromotion during short-range collisions. Computing the average motion of the test charge through the background plasma, the transverse force is found to strongly influence the trajectory by changing the gyroradius and the gyrofriction force is found to slightly change the gyrofrequency of the test charge resulting in a phase shift.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Visualizing interfacial collective reaction behaviour of Li–S batteries

Benefiting from high energy density (2,600 Wh kg -1 )) and low cost, lithium-sulfur (Li-S) batteries are considered promising candidates for advanced energy-storage systems. Despite tremendous efforts in suppressing the long-standing shuttle effect of lithium polysulfides, understanding of the interfacial reactions of lithium polysulfides at the nanoscale remains elusive. This is mainly because of the limitations of in situ characterization tools in tracing the liquid-solid conversion of unstable lithium polysulfides at high temporal-spatial resolution. There is an urgent need to understand the coupled phenomena inside Li-S batteries, specifically, the dynamic distribution, aggregation, deposition and dissolution of lithium polysulfides. Here, by using in situ liquid-cell electrochemical transmission electron microscopy, we directly visualized the transformation of lithium polysulfides over electrode surfaces at the atomic scale. Notably, an unexpected gathering-induced collective charge transfer of lithium polysulfides was captured on the nanocluster active-centre-immobilized surface. It further induced an instantaneous deposition of nonequilibrium Li2S nanocrystals from the dense liquid phase of lithium polysulfides. Without mediation of active centres, the reactions followed a classical single-molecule pathway, lithium polysulfides transforming into Li 2 S 2 and Li 2 S step by step. Molecular dynamics simulations indicated that the long-range electrostatic interaction between active centres and lithium polysulfides promoted the formation of a dense phase consisting of Li + and S n 2- (2 < n ≤ 6), and the collective charge transfer in the dense phase was further verified by ab initio molecular dynamics simulations. Finally, the collective interfacial reaction pathway unveils a new transformation mechanism and deepens the fundamental understanding of Li-S batteries.

25 ENERGY STORAGE↗
NASA NTRS2002

Tensile Strength of Carbon Nanotubes Under Realistic Temperature and Strain Rate

Strain rate and temperature dependence of the tensile strength of single-wall carbon nanotubes has been investigated with molecular dynamics simulations. The tensile failure or yield strain is found to be strongly dependent on the temperature and strain rate. A transition state theory based predictive model is developed for the tensile failure of nanotubes. Based on the parameters fitted from high-strain rate and temperature dependent molecular dynamics simulations, the model predicts that a defect free micrometer long single-wall nanotube at 300 K, stretched with a strain rate of 1%/hour, fails at about 9 plus or minus 1% tensile strain. This is in good agreement with recent experimental findings.

Wei, Chen-Yu↗
DOE OSTI2021

Small uranium and oxygen interstitial clusters in UO 2 : An empirical potential study

The aggregation of defect clusters in UO 2 under irradiations creates dislocation loops and cavities, which in turn contribute to swelling of the fuels, leading to detrimental pellet-clad mechanical interactions. The kinetics of the defect clusters, however, have received less attention and are not well understood. As such, using two interatomic potentials, the Morelon potential and a many-body model by Cooper, Rushton, and Grimes (CRG), the energetics and kinetics of small uranium interstitial clusters, with and without bound oxygen interstitials, are investigated with molecular dynamics and accelerated molecular dynamics simulations. For the single bound anti-Schottky interstitial cluster (U i 2O i ), both potentials find the interstitialcy migration mechanism to be a low barrier migration mechanism. However, for the CRG potential, an oxygen Frenkel pair assisted interstitial cluster migration mechanism with a lower migration barrier is identified during accelerated molecular dynamics simulations. For the Morelon potential, the interstitial clusters have geometrically compact structures, and small interstitial clusters (two and three anti-Schottky clusters) are more mobile than the single bound anti-Schottky interstitial cluster; for the CRG potential, the interstitial clusters have geometrically non-compact structures, yielding a more complex energy landscape. The rapid migration of interstitials and interstitial clusters are important for the irradiation response of UO 2 and correlates well with the very low barriers reported in a few historic experimental studies.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2020

Lithium oxidation and electrolyte decomposition at Li-metal/liquid electrolyte interfaces

In this work, we examine the evolution of events occurring when a Li metal surface is in contact with a 2 M solution of a Li salt in a solvent or mixture of solvents, via classical molecular dynamics simulations with a reactive force field allowing bond breaking and bond forming. The main events include Li oxidation and electrolyte reduction along with expansion of the Li surface layers forming a porous phase that is the basis for the formation of the solid-electrolyte interphase (SEI) components. Nucleation of the main SEI components (LiF, Li oxides, and some organics) is characterized. The analysis clearly reveals the details of these physical–chemical events as a function of time, during 20 nanoseconds. The effects of the chemistry of the electrolyte on Li oxidation and dissolution in the liquid electrolyte, and SEI nucleation and structure are identified by testing two salts: LiPF 6 and LiCF 3 SO 3 , and various solvents including ethers and carbonates and mixtures of them. The kinetics and thermodynamics of Li 6 F, the core nuclei in the LiF crystal, are studied by analysis of the MD trajectories, and via density functional theory calculations respectively. The SEI formed in this computational experiment is the “native” film that would form upon contact of the Li foil with the liquid electrolyte. As such, this work is the first in a series of computational experiments that will help elucidate the intricate interphase layer formed during battery cycling using metal anodes.

25 ENERGY STORAGE↗
DOE OSTI2022

Structural evolution in Au- and Pd-based metallic glass forming liquids and the case for improved molecular dynamics force fields

The results of a combined experimental and computational investigation of the structural evolution of Au 81 Si 19 , Pd 82 Si 18 , and Pd 77 Cu 6 Si 17 metallic glass forming liquids are presented. Electrostatically levitated metallic liquids are prepared, and synchrotron x-ray scattering studies are combined with embedded atom method molecular dynamics simulations to probe the distribution of relevant structural units. Metal–metalloid based metallic glass forming systems are an extremely important class of materials with varied glass forming ability and mechanical processibility. High quality experimental x-ray scattering data are in poor agreement with the data from the molecular dynamics simulations, demonstrating the need for improved interatomic potentials. The first peak in the x-ray static structure factor in Pd 77 Cu 6 Si 17 displays evidence for a Curie–Weiss type behavior but also a peak in the effective Curie temperature. In conclusion, a proposed order parameter distinguishing glass forming ability, 1/(S(T,q 1 )−1, shows a peak in the effective Curie temperature near a crossover temperature established by the behavior of the viscosity, T A .

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

Thermal coupling-decoupling mechanism of heat transfer across van der Waals interfaces in n-eicosane

Progress toward improving thermal transport properties of n-eicosane, a promising phase change material for thermal energy storage near human body temperature, is hindered by lack of a general theoretical framework to model thermal interfacial conductance (TIC) of the van der Waals (vdW) force bonded molecular interfaces, which have little compositional/structural mismatch. Combining molecular dynamics simulations at temperatures up to the melting point and a time-domain analysis of interfacial heat currents in individual molecules, we unveil that (1) the heat flux across a molecular vdW interface has an alternating coupling-decoupling pattern at the single molecule level due to thermal fluctuations, and (2) the net heat transfer during the thermally coupled periods is orders of magnitude larger than that during the thermally decoupled periods. Furthermore, we demonstrate that a vdW interface's TIC correlates almost linearly with two atomistic-scale coupling-decoupling parameters: αc and σc, that represent the duration and strength of the thermal coupling, respectively. This thermal coupling-decoupling mechanism proposes a new paradigm that explicitly accounts for the effects of structural flexibility and dynamical fluctuation on heat transfer across molecular interfaces, and will likely provide new insights into modeling thermal transport properties in a wide-ranging soft materials, including self-assembled monolayers, paraffins or lipids.

42 ENGINEERING↗
DOE OSTI2025

Long timescale solvation dynamics and confinement: The case of non-ionic deep eutectic solvents of lauric acid and N-methylacetamide

Microscopic segregation and molecular heterogeneities in complex liquids are the result of the interplay between different intermolecular forces, all of which contribute to the energy landscape of the system. A consequence of the intricate energy landscape is the nontrivial effect on the solvation dynamics. Here, the impact of molecular heterogeneities on the solvation dynamics is studied using infrared spectroscopies and molecular dynamics simulations. In particular, this study focuses on the dynamical effect of nanoscopic heterogeneities present in deep eutectic solvents (DESs) composed of lauric acid (LA) and N-methylacetamide (NMA). To this end, a molecular probe containing a carbon triple bond is used as an infrared reporter. The results show that the vibrational probe is likely to be located in the NMA polar domains. Furthermore, the probe solvation dynamics derived from the 2DIR spectra presents a slowdown of its timescale with increasing LA concentration in the DES. Kubo modeling of the probe solvation dynamics shows a correlation between the amplitude of its long time component and the presence of molecular heterogeneities in the sample. Semiclassical modeling of the vibrational line shape of the triple bond stretch demonstrates that the heterogeneities affect the whole solvation dynamics of the system through the amplitudes of the frequency fluctuations. Molecular dynamics simulations confirm the experimental results and their interpretation by showing a slowdown of the solvation dynamics when the LA heterogeneities are present. Overall, the study presents a molecular framework to explain the effect of confinement created by nanoscopic LA heterogeneities on the solvation dynamics of the system.

Chemical compounds↗
DOE OSTI2020

Comparison of fixed charge and polarizable models for predicting the structural, thermodynamic, and transport properties of molten alkali chlorides

Results from extensive molecular dynamics simulations of molten LiCl, NaCl, KCl, and RbCl over a wide range of temperatures are reported. Comparison is made between the “Polarizable Ion Model” (PIM) and the non-polarizable “Rigid Ion Model” (RIM). Densities, self-diffusivities, shear viscosities, ionic conductivities, and thermal conductivities are computed and compared with experimental data. In addition, radial distribution functions are computed from ab initio molecular dynamics simulations and compared with the two sets of classical simulations as well as experimental data. Here, the two classical models perform reasonably well at capturing structural and dynamic properties of the four molten alkali chlorides, both qualitatively and often quantitatively. With the singular exception of liquid density, for which the PIM is more accurate than the RIM, there are few clear trends to suggest that one model is more accurate than the other for the four alkali halide systems studied here.

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

Reinforcement learning based hybrid bond-order coarse-grained interatomic potentials for exploring mesoscale aggregation in liquid–liquid mixtures

Exploring mesoscopic physical phenomena has always been a challenge for brute-force all-atom molecular dynamics simulations. Although recent advances in computing hardware have improved the accessible length scales, reaching mesoscopic timescales is still a significant bottleneck. Coarse-graining of all-atom models allows robust investigation of mesoscale physics with a reduced spatial and temporal resolution but preserves desired structural features of molecules, unlike continuum-based methods. Here, we present a hybrid bond-order coarse-grained forcefield (HyCG) for modeling mesoscale aggregation phenomena in liquid–liquid mixtures. The intuitive hybrid functional form of the potential offers interpretability to our model, unlike many machine learning based interatomic potentials. We parameterize the potential with the continuous action Monte Carlo Tree Search (cMCTS) algorithm, a reinforcement learning (RL) based global optimizing scheme, using training data from all-atom simulations. The resulting RL-HyCG correctly describes mesoscale critical fluctuations in binary liquid–liquid extraction systems. cMCTS, the RL algorithm, accurately captures the mean behavior of various geometrical properties of the molecule of interest, which were excluded from the training set. The developed potential model along with the RL-based training workflow could be applied to explore a variety of other mesoscale physical phenomena that are typically inaccessible to all-atom molecular dynamics simulations.

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

The hTERT core promoter forms three parallel G-quadruplexes

The structure of the 68 nt sequence with G-quadruplex forming potential within the hTERT promoter is disputed. One model features a structure with three stacked parallel G-quadruplex units, while another features an unusual duplex hairpin structure adjoined to two stacked parallel and antiparallel quadruplexes. We report here the results of an integrated structural biology study designed to distinguish between these possibilities. As part of our study, we designed a sequence with an optimized hairpin structure and show that its biophysical and biochemical properties are inconsistent with the structure formed by the hTERT wild-type sequence. By using circular dichroism, thermal denaturation, nuclear magnetic resonance spectroscopy, analytical ultracentrifugation, small-angle X-ray scattering, molecular dynamics simulations and a DNase I cleavage assay we found that the wild type hTERT core promoter folds into a stacked, three-parallel G-quadruplex structure. The hairpin structure is inconsistent with all of our experimental data obtained with the wild-type sequence. All-atom models for both structures were constructed using molecular dynamics simulations. These models accurately predicted the experimental hydrodynamic properties measured for each structure. We found with certainty that the wild-type hTERT promoter sequence does not form a hairpin structure in solution, but rather folds into a compact stacked three-G-quadruplex conformation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Molecular dynamics and continuum analyses of the electrokinetic zeta potential in nanostructured slit channels

This work presents a theoretical and numerical study of electrokinetic flow and the zeta potential for the case of slit channels with nanoscale surface features of dimensions comparable to the Debye length, by employing molecular dynamics simulations and continuum-level analyses. Additionally, a simple analytical model for considering the average effect of such surface features producing nanoscale roughness is proposed by employing matched asymptotic solutions for the charge density and fluid flow field, and matching conditions that satisfy electroneutrality and the Onsager reciprocal relation between the electroosmotic flow rate and streaming current. The proposed analytical model quantitatively accounts for results from molecular dynamics simulations that consider the presence of ion solvation shells and surface hydration layers. Our analysis indicates that a simultaneous knowledge of the electroosmotic and pressure-driven flow rate or streaming current can be instrumental to unambiguously determine the zeta potential in the presence of nanoscale surface features.

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

Plasticity and Kinky Chemistry of Carbon Nanotubes

Since their discovery in 1991, carbon nanotubes have been the subject of intense research interest based on early predictions of their unique mechanical, electronic, and chemical properties. Materials with the predicted unique properties of carbon nanotubes are of great interest for use in future generations of aerospace vehicles. For their structural properties, carbon nanotubes could be used as reinforcing fibers in ultralight multifunctional composites. For their electronic properties, carbon nanotubes offer the potential of very high-speed, low-power computing elements, high-density data storage, and unique sensors. In a continuing effort to model and predict the properties of carbon nanotubes, Ames accomplished three significant results during FY99. First, accurate values of the nanomechanics and plasticity of carbon nanotubes based on quantum molecular dynamics simulations were computed. Second, the concept of mechanical deformation catalyzed-kinky-chemistry as a means to control local chemistry of nanotubes was discovered. Third, the ease of nano-indentation of silicon surfaces with carbon nanotubes was established. The elastic response and plastic failure mechanisms of single-wall nanotubes were investigated by means of quantum molecular dynamics simulations.

Srivastava, Deepak↗