Search NASA⌕ Search

SEARCH · Search NASA

Results for “molecular dynamic 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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Testing thermal gradient driving force for grain boundary migration using molecular dynamics simulations

Strong thermal gradients in low-thermal-conductivity ceramics may drive extended defects, such as grain boundaries and voids, to migrate in preferential directions. In this work, molecular dynamics simulations are conducted to study thermal gradient driven grain boundary migration and to verify a previously proposed thermal gradient driving force equation, using uranium dioxide as a model system. It is found that a thermal gradient drives grain boundaries to migrate up the gradient and the migration velocity increases under a constant gradient owing to the increase in mobility with temperature. Different grain boundaries migrate at very different rates due to their different intrinsic mobilities. The extracted mobilities from the thermal gradient driven simulations are compared with those calculated from two other well-established methods and good agreement between the three different methods is found, demonstrating that the theoretical equation of the thermal gradient driving force is valid, although a correction of one input parameter should be made. The discrepancy in the grain boundary mobilities between modeling and experiments is also discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Mean ionic activity coefficients in aqueous NaCl solutions from molecular dynamics simulations

The mean ionic activity coefficients of aqueous NaCl solutions of varying concentrations at 298.15 K and 1 bar have been obtained from molecular dynamics simulations by gradually turning on the interactions of an ion pair inserted into the solution. Several common non-polarizable water and ion models have been used in the simulations. Gibbs-Duhem equation calculations of the thermodynamic activity of water are used to confirm the thermodynamic consistency of the mean ionic activity coefficients. While the majority of model combinations predict the correct trends in mean ionic activity coefficients, they overestimate their values at high salt concentrations. The solubility predictions also suffer from inaccuracies, with all models underpredicting the experimental values, some by large factors. These results point to the need for further ion and water model development.

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

Theoretical design of a novel copper doped gold cluster supported on graphene utilizing ab initio molecular dynamics simulations

Ab initio molecular dynamics simulations have been used to inspect the adsorption of O{sub 2} to a small gold-copper alloy cluster supported on graphene. The exposed Cu atom in this cluster acts as a crucial attractive site for the approaching of O{sub 2} and consequently widens the reaction channel for the adsorption process. Conversely, a pure Au cluster on the same graphene support is inactive for the O{sub 2} adsorption because the corresponding reaction channel for the adsorption is very narrow. These results clearly indicate that doping a different metal to the Au cluster is a way to enhance the oxygen adsorption and to promote catalytic reactions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Protein–Mineral Interactions: Molecular Dynamics Simulations Capture Importance of Variations in Mineral Surface Composition and Structure

Molecular dynamics simulations, conventional and metadynamics, were performed to determine the interaction of model protein Gb1 over kaolinite (001), Na+-montmorillonite (001), Ca2+-montmorillonite (001), goethite (100), and Na+-birnessite (001) mineral surfaces. Gb1, a small (56 residue) protein with a well-characterized solution-state nuclear magnetic resonance (NMR) structure and having α-helix, four-fold β-sheet, and hydrophobic core features, is used as a model protein to study protein soil mineral interactions and gain insights on structural changes and potential degradation of protein. From our simulations, we observe little change to the hydrated Gb1 structure over the kaolinite, montmorillonite, and goethite surfaces relative to its solvated structure without these mineral surfaces present. Over the Na+-birnessite basal surface, however, the Gb1 structure is highly disturbed as a result of interaction with this birnessite surface. Unraveling of the Gb1 β-sheet at specific turns and a partial unraveling of the α-helix is observed over birnessite, which suggests specific vulnerable residue sites for oxidation or hydrolysis possibly leading to fragmentation.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2020

Lithiation of Sulfur-Graphene Compounds Using Reactive Force-Field Molecular Dynamics Simulations

We performed molecular dynamics simulations of lithium-sulfur-graphene compounds using reactive force fields, providing a time scale to observe atomistic features relevant to the microscopic behavior of the of the bulk of sulfur-based cathodes to be used beyond our present Li-ion batteries. The samples we used were set to realistic geometries through sophisticated protocols to simulate ultrafast reactions that occur within the picosecond range, thus allowing us to get some insights into the characteristics of the bulk material in working cathodes of Li–S batteries, which are mixed with carbon to increase the poor electronic conductivity of S. We report chemical speciation and geometrical data at atomistic levels. We observed that slowly lithiated cathodes were more stable and with higher density than those that were suddenly fully-lithiated. We did not observe molecular Li 2 S formation; however, we observed an amorphous solid arrangement with the same stoichiometry of Li and S, with S–Li–S angles of ~111° and smaller ones due to the interaction between polysulfides that did not reacted totally. In addition, graphene keeps its planar shape; however, S8 changes its shape from rings to chains. Lithiated structures are more stable with lower energies, and more close-packed structures than structures with Li already inserted.

25 ENERGY STORAGE↗
DOE OSTI2018

Guest–Host Interactions in Mixed CH 4 –CO 2 Hydrates: Insights from Molecular Dynamics Simulations

Here, classical molecular dynamics simulations of mixed CH 4 –CO 2 gas hydrates provide a thorough analysis of the energetics of both pure CH 4 and CO 2 hydrates and three intermediate compositions along the (CH 4 ) 1–x (CO 2 ) x ·5.75(H 2 O) solid solution. The energy is broken into guest–guest, guest–host, and host–host contributions. Radial distribution functions and three-dimensional density distributions provide insight into the changes in guest orientation and interactions with the host framework as a function of guest composition. Both energetic and structural analyses provide complementary information to previous experimental studies of the system. The experimentally observed isotropic orientation of CH 4 molecules in both small and large cages is confirmed and the description of the anisotropic orientation of CO 2 molecules in the large cage is confirmed and further enhanced. In mixed hydrates, the presence of CH 4 alters the orientation of CO 2 , an indication that the interactions between guests, either direct (guest–guest) or mediated through interactions with the host (host–guest), are an important phenomenon in these systems.

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

Ab Initio Molecular Dynamics Simulations of Amorphous Calcium Carbonate: Interpretation of Pair Distribution Function and X-Ray Absorption Spectroscopy Data

The structure and transformation of hydrous amorphous calcium carbonate (ACC) are key to understanding biomineralization pathways and their relationship with the properties of the resulting material. Quantitative interpretation of scattering experiments aimed at elucidating the structure of ACC is challenging, due to the amorphous nature of this material, and, therefore, requires models for the structure and the scattering physics. Here, we generate physically realistic ensembles of hydrated ACC structures and their vibrational disorder from ab initio molecular dynamics (AIMD) simulations with an emphasis on enabling the consistent interpretation of the finer details of three complementary structural probes: neutron and x-ray pair distribution experiments and x-ray absorption spectroscopy (XAS). In each case, we simulate the signal directly in reciprocal space and then manipulate it into the real-space pair distribution function (PDF) or spectrum using the same procedures for the experimental and theoretical data. Good agreement with experiment was obtained across the three techniques with the simulations accounting well for all features in the spectra. Remaining small discrepancies pointed to differences between real samples and the idealized simulated systems such as deviations from the nominal CaCO3·nH2O stoichiometry. Additionally, the simulations offered a more accurate description of the local coordination environment of calcium than previous shell-by-shell fits to spectra of synthetic ACC and classical molecular dynamics simulations. This work demonstrates that AIMD is a powerful approach for extracting detailed structural information from neutron PDF, x-ray PDF, and XAS of amorphous carbonate phases.

Prange, Micah P.↗
DOE OSTI2015

Extended Lagrangian Born-Oppenheimer molecular dynamics simulations of the shock-induced chemistry of phenylacetylene

The initial chemical events that occur during the shock compression of liquid phenylacetylene have been investigated using self-consistent tight binding molecular dynamics simulations. The extended Lagrangian Born-Oppenheimer molecular dynamics formalism enabled us to compute microcanonical trajectories with precise conservation of the total energy. Our simulations revealed that the first density-increasing step under shock compression arises from the polymerization of phenylacetylene molecules at the acetylene moiety. The application of electronic structure-based molecular dynamics with long-term conservation of the total energy enabled us to identify electronic signatures of reactivity via monitoring changes in the HOMO-LUMO gap, and to capture directly adiabatic shock heating, transient non-equilibrium states, and changes in temperature arising from exothermic chemistry in classical molecular dynamics trajectories.

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

GPU-accelerated Tersoff potentials for massively parallel Molecular Dynamics simulations

The Tersoff potential is one of the empirical many-body potentials that has been widely used in simulation studies at atomic scales. Unlike pair-wise potentials, the Tersoff potential involves three-body terms, which require much more arithmetic operations and data dependency. In this contribution, we have implemented the GPU-accelerated version of several variants of the Tersoff potential for LAMMPS, an open-source massively parallel Molecular Dynamics code. Compared to the existing MPI implementation in LAMMPS, the GPU implementation exhibits a better scalability and offers a speedup of 2.2X when run on 1000 compute nodes on the Titan supercomputer. On a single node, the speedup ranges from 2.0 to 8.0 times, depending on the number of atoms per GPU and hardware configurations. The most notable features of our GPU-accelerated version include its design for MPI/accelerator heterogeneous parallelism, its compatibility with other functionalities in LAMMPS, its ability to give deterministic results and to support both NVIDIA CUDA- and OpenCL-enabled accelerators. Our implementation is now part of the GPU package in LAMMPS and accessible for public use.

Computer Science↗
DOE OSTI2015

Ensemble Sampling vs. Time Sampling in Molecular Dynamics Simulations of Thermal Conductivity

In this report we compare time sampling and ensemble averaging as two different methods available for phase space sampling. For the comparison, we calculate thermal conductivities of solid argon and silicon structures, using equilibrium molecular dynamics. We introduce two different schemes for the ensemble averaging approach, and show that both can reduce the total simulation time as compared to time averaging. It is also found that velocity rescaling is an efficient mechanism for phase space exploration. Although our methodology is tested using classical molecular dynamics, the ensemble generation approaches may find their greatest utility in computationally expensive simulations such as first principles molecular dynamics. For such simulations, where each time step is costly, time sampling can require long simulation times because each time step must be evaluated sequentially and therefore phase space averaging is achieved through sequential operations. On the other hand, with ensemble averaging, phase space sampling can be achieved through parallel operations, since each ensemble is independent. For this reason, particularly when using massively parallel architectures, ensemble sampling can result in much shorter simulation times and exhibits similar overall computational effort.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2012

Hydration structures of U(III) and U(IV) ions from ab initio molecular dynamics simulations

We apply DFT+U-based ab initio molecular dynamics simulations to study the hydration structures of U(III) and U(IV) ions, pertinent to redox reactions associated with uranium salts in aqueous media. U(III) is predicted to be coordinated to 8 water molecules, while U(IV) has a hydration number between 7 and 8. At least one of the innershell water molecules of the hydrated U(IV) complex becomes spontaneously deprotonated. As a result, the U(IV)-O pair correlation function exhibits a satellite peak at 2.15 A associated with the shorter U(IV)-(OH{sup -}) bond. This feature is not accounted for in analysis of extended x-ray absorption fine structure and x-ray adsorption near edge structure measurements, which yield higher estimates of U(IV) hydration numbers. This suggests that it may be useful to include the effect of possible hydrolysis in future interpretation of experiments, especially when the experimental pH is close to the reported hydrolysis equilibrium constant value.

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

Atomic Diffusion, Segregation, and Grain Boundary Migration in Nickel-Based Alloys from Molecular Dynamics Simulations

Grain boundary diffusion and metal mobility in alloys control material performance in many applications and yet remain poorly understood at a mechanistic level. With advances in accessible time and length scales for computational molecular simulations, and recent force field developments, we now possess tools to help unravel those mechanisms. Using large-scale molecular dynamics simulations, here we examined vacancy-mediated diffusion processes in Ni-5Cr alloy with low and high-energy grain boundaries. We show that atomic diffusion inside the grain boundary plane is about four times higher than bulk diffusion, at any temperature, and exhibits a typical Arrhenius behavior with a very small energy barrier (0~.8 eV for Cr and 0.7 eV for Ni within 1300-1600 K). Additionally, the fastest diffusing species inverts; Cr diffusion was faster than Ni in the bulk but slower in the grain boundaries. This is attributed to the creation of high cohesive energy clusters of Cr at the grain boundary. Grain boundary migration was also observed to be temperature dependent and appears to be two times higher in the 5% Cr alloy than in pure Ni, highlighting the important role of the alloying element on grain boundary motion.

Simonnin, Pauline GN↗
DOE OSTI2017

Error correction in multi-fidelity molecular dynamics simulations using functional uncertainty quantification

We use functional, Fréchet, derivatives to quantify how thermodynamic outputs of a molecular dynamics (MD) simulation depend on the potential used to compute atomic interactions. Our approach quantifies the sensitivity of the quantities of interest with respect to the input functions as opposed to its parameters as is done in typical uncertainty quantification methods. We show that the functional sensitivity of the average potential energy and pressure in isothermal, isochoric MD simulations using Lennard–Jones two-body interactions can be used to accurately predict those properties for other interatomic potentials (with different functional forms) without re-running the simulations. This is demonstrated under three different thermodynamic conditions, namely a crystal at room temperature, a liquid at ambient pressure, and a high pressure liquid. The method provides accurate predictions as long as the change in potential can be reasonably described to first order and does not significantly affect the region in phase space explored by the simulation. The functional uncertainty quantification approach can be used to estimate the uncertainties associated with constitutive models used in the simulation and to correct predictions if a more accurate representation becomes available.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗