Search NASA⌕ Search

SEARCH · Search NASA

Results for “Molecular 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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
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 OSTI2023

Code modernization strategies for short-range non-bonded molecular dynamics simulations

Modern HPC systems are increasingly relying on greater core counts and wider vector registers. Thus, applications need to be adapted to fully utilize these hardware capabilities. One class of applications that can benefit from this increase in parallelism are molecular dynamics simulations. In this paper, we describe our efforts at modernizing the ESPResSo++ simulation package for molecular dynamics by restructuring its particle data layout for efficient memory accesses and applying vectorization techniques to benefit the calculation of short-range non-bonded forces, which results in an overall three times speedup and serves as a baseline for further optimizations. We also implement fine-grained parallelism for multi-core CPUs through HPX, a C++ runtime system which uses lightweight threads and an asynchronous many-task approach to maximize concurrency. Our goal is to evaluate the performance of an HPX-based approach compared to the bulk-synchronous MPI-based implementation. This requires the introduction of an additional layer to the domain decomposition scheme that defines the task granularity. On spatially inhomogeneous systems, which impose a corresponding load-imbalance in traditional MPI-based approaches, we demonstrate that by choosing an optimal task size, the efficient work-stealing mechanisms of HPX can overcome the overhead of communication resulting in an overall 1.4 times speedup compared to the baseline MPI version.

97 MATHEMATICS AND COMPUTING↗
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 OSTI2025

Initial position optimization in molecular dynamics simulations for a Coulomb system

A new algorithm for molecular dynamics (MD) simulations is developed to optimize plasma particle distributions at given initial temperatures. By combining velocity scaling and reassignment, the method effectively eliminates the initial rise and oscillation in temperatures observed with randomly distributed positions. These rises and oscillations are undesired numerical artifacts observed in conventional plasma MD simulations, arising from unoptimized particle positions. The algorithm demonstrates temperature relaxation without initial rises or oscillations, as well as precise flow velocity relaxation, enabling accurate measurement of relaxation times. The code is accelerated using graphics processing units for parallel processing, enhancing the study of plasma dynamics. The proposed method for distributing physically valid particles in MD simulations enables accurate studies of intrinsic collision processes in plasmas, including the dynamics of strongly coupled plasmas, plasma–wave interactions, and transport phenomena in magnetized plasmas. The paper concludes with a discussion of potential applications and future enhancements to the algorithm.

Jo, Jawon (ORCID:0009000924193285)↗
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↗
DOE OSTI2020

Hydration and reaction mechanisms on sodium silicate glass surfaces from molecular dynamics simulations with reactive force fields

Recent development of reactive force fields have enabled molecular dynamics simulations of interactions between silicate glasses and water at the atomistic scale. While multicomponent silicate glasses encompass a wide variety of compositions and properties, one common structural feature in these glasses is the combination of the network structure that is made up of silica tetrahedra linked through corner sharing interspersed with network modifiers like alkali and alkaline-earth ions that break up the Si–O–Si linkages by forming nonbridging oxygen. In reactions with water, ion exchange between alkali ions in the glass and proton or hydronium in the solution, as well as hydrolysis reaction of the Si–O–Si linkages and subsequent silanol formation, is observed and well documented. Furthermore, we have used a set of recently developed reactive force field to investigate the reactions between water and the surfaces of silica and sodium silicate glasses of different compositions for reactions up to 8 nanoseconds. Our results indicate sodium leaching into water and diffusion of water molecules up to 25 Å into the glass surface. We examined the structural and compositional changes inside the glass and around the diffused ions and use these to explain the rates of silanol formation at the surface. We also observed proton transport in the glass which has an indirect influence on the silanol formation rates. While the surface of the glass was rough to start with, it undergoes further modification into a hydrated gel-like structure in the glass for up to 5 Å in the higher alkali containing glasses. It was found that the leached sodium ions remain close to the interface and that fragments of silicate network from the surface is capable of dislodging from the bulk glass and enter the aqueous solution. These simulations thus provide insights into the formation and structure of an alteration layers commonly observed in multicomponent silicate glasses corroded in aqueous solutions.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Steering Molecular Dynamics Simulations of Membrane-Associated Proteins with Neutron Reflection Results

Herein, we present a novel method to incorporate structural results from surface-sensitive scattering, such as X-ray or neutron reflectometry, into molecular dynamics simulations. While reflectometry techniques generally provide a means to determine the molecular-scale structures of organized interfacial films, they were recently shown to offer the capability to characterize the structures of protein-membrane complexes supported by a solid substrate. One-dimensional information inherent in the experimental results is used in the form of component volume occupancy (CVO) profiles, which describe the distribution of molecular components within an interfacial architecture, to construct real-space constraints in the form of a biasing potential for the simulation that vanishes when the simulated and experimental profiles agree. This approach improves the correspondence between simulation and experiment, as shown in the re-evaluation of an neutron-reflection-derived structure which was approximated by an independent molecular dynamics simulation in earlier work, and it also leads to faster equilibration of ensemble structures. We further show that time averaging the CVO profile that develops in the simulation while biasing with this approach permits fluctuations about the average that are necessary for conformational exploration of the system. This method is particularly valuable for studies of proteins at interfaces that contain disordered regions since the conformation of such regions is difficult to judge from the analysis of one-dimensional experimental profiles and may take prohibitively long to equilibrate in simulations.

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

Large-scale molecular dynamics simulations of bubble collapse in water: Effects of system size, water model, and nitrogen

Molecular dynamics simulations in the microcanonical ensemble are performed to study the collapse of a bubble in liquid water using the single-site mW and the four-site TIP4P/2005 water models. To study system size effects, simulations for pure water systems are performed using periodically replicated simulation boxes with linear dimensions, L, ranging from 32 to 512 nm with the largest systems containing 8.7 × 10 6 and 4.5 × 10 9 molecules for the TIP4P/2005 and mW water models, respectively. The computationally more efficient mW water model allows us to reach converging behavior when the bubble dynamics results are plotted in reduced units, and the limiting behavior can be obtained through linear extrapolation in L –1 . Qualitative differences are observed between simulations with the mW and TIP4P/2005 water models, but they can be explained by the models’ differences in predicted viscosity and surface tension. Although bubble collapse occurs on time scales of only hundreds of picoseconds, the system sizes used here are sufficiently large to obtain bubble dynamics consistent with the Rayleigh–Plesset equation when using the models’ thermophysical properties as input. For the conditions explored here, extreme heating of the interfacial water molecules near the time of collapse is observed for the larger mW water systems (but the model underpredicts the viscosity), whereas heating is less pronounced for the TIP4P/2005 water systems because its larger viscosity contribution slows the collapse dynamics. The presence of nitrogen within the bubble only starts to affect bubble dynamics near the very end of the initial collapse, leading to an incomplete collapse and strong rebound for the mW water model. Although nitrogen is non-condensable at 300 K, it becomes highly compressed and reaches a liquid-like density near the collapse point. We find that the dissolution of nitrogen is much slower than the movement of the collapsing water front, and the re-expansion of the dense nitrogen droplet gives rise to bubble rebound. The incompatibility of the collapse and dissolution time scales should be considered for continuum-scale modeling of bubble dynamics. Finally, we also confirm that the diffusion coefficient for dissolved nitrogen is insensitive to pressure as the liquid transitions from a compressed to a stretched state.

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

The Stress-strain Behavior of Polymer-Nanotube Composites from Molecular Dynamics Simulations

Stress-strain curves of polymer-carbon nanotube composites are derived from molecular dynamics simulations of a single-walled carbon nanotube embedded in polyethylene. A comparison is made between the response to mechanical loading of a composite with a long, continuous nanotube (replicated via periodic boundary conditions) and the response of a composite with a short, discontinuous nanotube. Both composites are mechanically loaded in the direction of and transverse to the NT axis. The long-nanotube composite shows an increase in the stiffness relative to the polymer and behaves anisotropically under the different loading conditions. The short-nanotube composite shows no enhancement relative to the polymer, most probably because of its low aspect ratio. The stress-strain curves are compared with rule-of-mixtures predictions.

Frankland, S. J. V.↗
DOE OSTI2023

Coarse-grained explicit-solvent molecular dynamics simulations of semidilute unentangled polyelectrolyte solutions

In this study, we present results from explicit-solvent coarse-grained molecular dynamics (MD) simulations of fully charged, salt-free, and unentangled polyelectrolytes in semidilute solutions. The inclusion of a polar solvent in the model allows for a more physical representation of these solutions at concentrations, where the assumptions of a continuum dielectric medium and screened hydrodynamics break down. The collective dynamic structure factor of polyelectrolytes, S(q, t), showed that at q > q*, where q* = 2π/ξ is the polyelectrolyte peak in the structure factor S(q) and ξ is the correlation length, the relaxation time obtained from fits to stretched exponential was $\tau$ KWW ~ q -3 , which describes unscreened Zimm-like dynamics. This is in contrast to implicit-solvent simulations using a Langevin thermostat where $\tau$ KWW ~ q -2 . At q < q*, a crossover region was observed that eventually transitions to another inflection point $\tau$ KWW ~ q -2 at length scales larger than ξ for both implicit- and explicit-solvent simulations. The simulation results were also compared to scaling predictions for correlation length, ξ ~ c$-½\atop{p}$, specific viscosity, η sp ~ c$½\atop{p}$, and diffusion coefficient, D ~ c$0\atop{p}$, where c p is the polyelectrolyte concentration. The scaling prediction for ξ holds; however, deviations from the predictions for η sp and D were observed for systems at higher c p , which are in qualitative agreements with recent experimental results. This study highlights the importance of explicit-solvent effects in molecular dynamics simulations, particularly in semidilute solutions, for a better understanding of polyelectrolyte solution behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Molecular dynamics simulation of metallic Al–Ce liquids using a neural network machine learning interatomic potential

Al-rich Al-Ce alloys have the possibility of replacing heavier steel and cast-irons for use in high-temperature applications. Knowledge about the structures and properties of Al-Ce alloys at liquid state is vital for optimizing the manufacture process to produce desired allows. However, reliable molecular dynamics simulation of Al-Ce alloy systems remains a great challenge due to the lack of accurate Al-Ce interatomic potential. In this work, an artificial neural network (ANN) deep machine learning (ML) method is used to develop a reliable interatomic potential for Al-Ce alloy. Ab initio molecular dynamics (AIMD) simulation data on Al-Ce liquid with small unit cell (~200 atoms) and on the known Al-Ce crystalline compounds are collected to train the interatomic potential using ANN-ML. The obtained ANN-ML model reproduces well the energies, forces, and atomic structure of Al 90 Ce 10 liquid and crystalline phases of Al-Ce compounds in comparison with ab initio results. The developed ANN-ML potential is applied in molecular dynamics simulations to study the structures and properties of metallic Al 90 Ce 10 liquid, which would provide useful insight for guiding experimental process to produce desired Al-Ce allows.

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

Molecular Dynamics Simulation of Carbon Nanotube Based Gears

We used molecular dynamics to investigate the properties and design space of molecular gears fashioned from carbon nanotubes with teeth added via a benzyne reaction known to occur with C60. A modified, parallelized version of Brenner's potential was used to model interatomic forces within each molecule. A Leonard-Jones 6-12 potential was used for forces between molecules. One gear was powered by forcing the atoms near the end of the buckytube to rotate, and a second gear was allowed.to rotate by keeping the atoms near the end of its buckytube on a cylinder. The meshing aromatic gear teeth transfer angular momentum from the powered gear to the driven gear. A number of gear and gear/shaft configurations were simulated. Cases in vacuum and with an inert atmosphere were examined. In an extension to molecular dynamics technology, some simulations used a thermostat on the atmosphere while the hydrocarbon gear's temperature was allowed to fluctuate. This models cooling the gears with an atmosphere. Results suggest that these gears can operate at up to 50-100 gigahertz in a vacuum or inert atmosphere at room temperature. The failure mode involves tooth slip, not bond breaking, so failed gears can be returned to operation by lowering temperature and/or rotation rate. Videos and atomic trajectory files in xyz format are presented.

Han, Jie↗
DOE OSTI2020

Fast predictions of liquid-phase acid-catalyzed reaction rates using molecular dynamics simulations and convolutional neural networks

The rates of liquid-phase, acid-catalyzed reactions relevant to the upgrading of biomass into high-value chemicals are highly sensitive to solvent composition and identifying suitable solvent mixtures is theoretically and experimentally challenging. We show that the complex atomistic configurations of reactant–solvent environments generated by classical molecular dynamics simulations can be exploited by 3D convolutional neural networks to enable accurate predictions of Brønsted acid-catalyzed reaction rates for model biomass compounds. We develop a 3D convolutional neural network, which we call SolventNet, and train it to predict acid-catalyzed reaction rates using experimental reaction data and corresponding molecular dynamics simulation data for seven biomass-derived oxygenates in water–cosolvent mixtures. We show that SolventNet can predict reaction rates for additional reactants and solvent systems an order of magnitude faster than prior simulation methods. This combination of machine learning with molecular dynamics enables the rapid, high-throughput screening of solvent systems and identification of improved biomass conversion conditions.

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

A fast recursive algorithm for molecular dynamics simulation

The present recursive algorithm for solving molecular systems' dynamical equations of motion employs internal variable models that reduce such simulations' computation time by an order of magnitude, relative to Cartesian models. Extensive use is made of spatial operator methods recently developed for analysis and simulation of the dynamics of multibody systems. A factor-of-450 speedup over the conventional O(N-cubed) algorithm is demonstrated for the case of a polypeptide molecule with 400 residues.

Jain, A.↗
DOE OSTI2025

Machine learning-accelerated path integral molecular dynamics simulations of reactive organic electrolytes

Hydrogen bonded electrolytes that exhibit accelerated proton transport via sequential reactive hops have drawn interest for their promise in clean energy applications. Molecular dynamics simulations of these electrolytes offer the opportunity to uncover microscopic mechanistic details that could be used to design and tune the properties of candidate electrolyte technologies. However, accurately modeling the proton transfer reactions and transport properties that give rise to high charge conductivites in these electrolytes proves computationally challenging because of the need to perform lengthy condensed phase simulations, treating both the electronic and nuclear degrees of freedom quantum mechanically. In this paper, we demonstrate that such a modeling task can be efficiently achieved with the use of density functional theory (DFT)-trained machine learning potentials (MLP) to accelerate path integral molecular dynamics (PIMD) simulations. We highlight the practical utility of this approach by using it to benchmark how closely PIMD simulations employing different DFT exchange–correlation functionals reproduce the composition-dependent densities, diffusion coefficients, and electrical conductivities of mixtures consisting of imidazole and levulinic acid. Even with the speedup afforded by our MLPs, PIMD simulations remain quite expensive. Furthermore, in order to render PIMD more computationally tractable, we introduce and benchmark the accuracy of a ring polymer contraction approach that leverages a computationally efficient short-range MLP to accelerate our PIMD simulations by an additional factor of four.

Chemical bonding↗
NASA NTRS2000

Molecular Dynamics Simulations of Carbon Nanotubes in Water

We study the hydrophobic/hydrophilic behavior of carbon nanotubes using molecular dynamics simulations. The energetics of the carbon-water interface are mainly dispersive but in the present study augmented with a carbon quadrupole term acting on the charge sites of the water. The simulations indicate that this contribution is negligible in terms of modifying the structural properties of water at the interface. Simulations of two carbon nanotubes in water display a wetting and drying of the interface between the nanotubes depending on their initial spacing. Thus, initial tube spacings of 7 and 8 A resulted in a drying of the interface whereas spacing of > 9 A remain wet during the course of the simulation. Finally, we present a novel particle-particle-particle-mesh algorithm for long range potentials which allows for general (curvilinear) meshes and "black-box" fast solvers by adopting an influence matrix technique.

Walther, J. H.↗
NASA NTRS1994

Dielectric Properties of Poly(ethylene oxide) from Molecular Dynamics Simulations

The order, conformations and dynamics of poly(oxyethylene) (POE) melts have been investigated through molecular dynamics simulations. The potential energy functions were determined from detailed ab initio electronic structure calculations of the conformational energies of the model molecules 1,2-dimethoxyethane (DME) and diethylether. The x-ray structure factor for POE from simulation will be compared to experiment. In terms of conformation, simulations reveal that chains are extended in the melt relative to isolated chains due to the presence of strong intermolecular O...H interactions, which occur at the expense of intramolecular O...H interactions. Conformational dynamics about the C-C bond were found to be significantly faster than in polymethylene, while conformational dynamics about the C-O bond even faster than the C-C dynamics. The faster local dynamics in POE relative to polymethylene is consistent with C-13 NMR spin-lattice relaxation experiments. Conformational transitions showed significant second-neighbor correlation, as was found for polymethylene. This correlation of transitions with C-C neighbors was found to be reduced relative to C-O neighbors. Dielectric relaxation from simulation will also be compared with experiment.

Smith, Grant D.↗
DOE OSTI2020

Efficient perturbation-tracking method for directly probing the spectral phonon properties from molecular dynamics simulations

Existing methods for directly extracting the spectral phonon properties from molecular dynamics (MD) simulations, like the normal mode analysis (NMA) and spectral energy density analysis, all require a very long simulation time to produce reliable results with good convergence. So far, these methods are mainly applied in studies using small systems and with empirical potentials, as the heavy computational load has greatly hindered their further applications. Here we propose a perturbation-tracking (PT) method for directly probing the mode-wise phonon anharmonic frequencies and lifetimes. We show that results obtained from our method are in excellent agreement with those from the conventional NMA approach, using Si as the model material system. Comparing with the NMA approach, the PT method offers a greater accuracy and significant improvement of efficiency. It takes an average of two orders of magnitude and up to three orders of magnitude less simulation time to obtain the same lifetime result of a phonon mode with intermediate to high accuracy. Meanwhile, our method preserves all the dynamics of probed phonon mode from a particular state, which means it is capable of studying the transient thermal transport processes in a nonequilibrium system. Besides the exceptional efficiency, our method also comes with freedom to choose to probe only those modes of interest. This makes it ideal for use with large systems and in computationally demanding applications, such as ab initio MD simulations. Moreover, the PT method we propose here is very straightforward and easy to implement.

74 ATOMIC AND MOLECULAR PHYSICS↗