SEARCH · Engineering Papers
Results for “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.
Insights into Spontaneous Solid Electrolyte Interphase Formation at Magnesium Metal Anode Surface from Ab Initio Molecular Dynamics Simulations
Spontaneous chemical reactivity at multivalent (Mg, Ca, Zn, Al) electrode surfaces is critical to solid electrolyte interphase (SEI) formation, and hence, directly affects the longevity of batteries. Here, we report an investigation of the reactivity of 0.5 M Mg(TFSI)2 in 1,2-dimethoxyethane (DME) solvent at a Mg(0001) surface using ab initio molecular dynamics (AIMD) simulations and detailed Bader charge analysis. Based on the simulations, the initial degradation reactions of the electrolyte strongly depend on the structure of the Mg(TFSI)2 species near the anode surface. At the surface, the dissociation of Mg(TFSI)2 species occurs via cleavage of the N-S bond for the solvent separated ion pair (SSIP) and via cleavage of the C-S bond for the contact ion pair (CIP) configuration. In the case of the CIP, both TFSI anions undergo spontaneous bond dissociation reactions to form atomic O, C, S, F, and N species adsorbed on the surface of the Mg anode. These products indicate that the initial SEI layer formed on the surface of the pristine Mg anode consists of a complex mixture of multiple components such as oxides, carbides, sulfides, fluorides, and nitrides. We believe that the atomic level insights gained from these simulations will lay the groundwork for the rational design of tailored and functional interphases that are critical for the success of multivalent battery technology.
Unsupervised learning of representative local atomic arrangements in molecular dynamics data
Molecular dynamics (MD) simulations present a data-mining challenge, given that they can generate a considerable amount of data but often rely on limited or biased human interpretation to examine their information content. By not asking the right questions of MD data we may miss critical information hidden within it. Here we combine dimensionality reduction (UMAP) and unsupervised hierarchical clustering (HDBSCAN) to quantitatively characterize prevalent coordination environments of chemical species within MD data. By focusing on local coordination, we significantly reduce the amount of data to be analyzed by extracting all distinct molecular formulas within a given coordination sphere. We then efficiently combine UMAP and HDBSCAN with alignment or shape-matching algorithms to partition these formulas into structural isomer families indicating their relative populations. The method was employed to reveal details of cation coordination in electrolytes based on molecular liquids.
New insights into the mechanism governing the elasticity of calcium silicate hydrate gels exposed to high temperature: A molecular dynamics study
When exposed to fire, the integrity of cement-based materials is governed by thermally-induced changes in the mechanical properties of their binding phase, i.e., the calcium–silicate–hydrate (C–S–H) gel. However, the effect of temperature on the structure, density, and mechanical properties of C–S–H remains only partially known. Here, based on reactive molecular dynamics simulations, we reveal the nature of thermally-induced damage in C–S–H gels. In general, we show that, at the atomic scale, exposure to high temperature results in partial dehydration, volumetric shrinkage, disordering, and stiffening in the C–S–H grains. However, we show that the thermal response of C–S–H strongly depends on its chemical composition, wherein C–S–H systems associated with lower Ca/Si molar ratios are able to undergo higher temperatures before amorphization. Based on these results, we demonstrate that the stiffness of C–S–H gels (i.e., including porosity—as probed by nanoindentation) is governed by a competition between the stiffening of the grains and the decrease in packing density—wherein the latter eventually become predominant.
Structure/property relationship of semicrystalline polymers during tensile deformation: a molecular dynamics approach
Molecular dynamics (MD) simulation with an ab initio numerical model is conducted to study the structure–property relationship of semicrystalline polymers during uniaxial deformation. The effects of chain length, temperature, and strain rate on mechanical properties are discussed. The influences of microstructural evolution such as bond length, entanglement density, and chain orientation are also studied quantitatively. The temperature will greatly affect the chain conformation in amorphous domains, and the results revealed that the interaction of amorphous and crystalline domains played a crucial role during stretching. Accordingly, the yielding of semicrystalline polymers follows different mechanisms at temperatures above and below T g . A melt-recrystallization scheme is observed during yielding at higher temperatures, while destruction of crystal structures is observed at lower temperatures at the yield point. The correlated effects of different temperatures and strain rates on mechanical properties are examined. This research is part of our efforts to develop a digital twin of a real experiment for efficient optimization of polymer material properties.
Oriented attachment mechanism of triangular Ag nanoplates: a molecular dynamics study
Molecular dynamics simulations predict that PVP-capped silver nanoplates exhibit lateral oriented attachment.
Pathways for Electron Transfer at MgO–Water Interfaces from Ab Initio Molecular Dynamics
The nature of electron transfer across metal oxide–water interfaces depends significantly on the band gap of the oxide and its band edge energies relative to the potentials of relevant aqueous redox couples. Here we focus on the water interface with MgO, a prototypical wide band gap oxide whose conduction band edge is close in energy to that of water. We investigate the behavior of an excess electron at and out of equilibrium near the interface using ab initio molecular dynamics based on hybrid density functional theory. Our simulations show that under equilibrium conditions the excess electron (donated by an Al impurity in MgO) localizes to a midgap defect state comparable in energy and shape to a hydrated electron in bulk water. To characterize the electron transfer from the conduction band of MgO to interfacial product states, we dope near-equilibrium configurations of the pristine MgO–water system with Al and run short trajectories of these instantaneously out-of-equilibrium systems. Here we observe two distinct products associated with the excess electron: a surface-localized electron (e surf –) and an aqueous hydrogen radical (H • ). The H • pathway exhibits a much higher activation barrier despite being more exoergic, making e surf – the kinetic product. Our characterization of the pathways on the basis of Marcus theory is consistent with the poor observed utility of MgO for water radiolysis. Moreover, we anticipate that the computational framework employed here will be broadly applicable to assessing electron transfer mechanisms at aqueous, photocatalytic interfaces.
Quantitative structure-property relationship (QSPR) analysis of calcium aluminosilicate glasses based on molecular dynamics simulations
Not provided.
Resolving the organization of CO 2 molecules confined in silica nanopores using in situ small-angle neutron scattering and molecular dynamics simulations
This work quantifies the core–shell structure of nano-confined CO 2 molecules under pressure in silica pores using in operando small-angle neutron scattering (SANS) measurements and molecular dynamics simulations.
Molecular dynamics simulations of a hydrophilic MIL-160-based membrane demonstrate pressure-dependent selective uptake of industrially relevant greenhouse
Oral presentation for ACS Spring Meeting
Enhancing the formation of ionic defects to study the ice Ih/XI transition with molecular dynamics simulations
Ice Ih, the common form of ice in the biosphere, contains proton disorder. Its proton-ordered counterpart, ice XI, is thermodynamically stable below 72 K. However, the formation of ice XI is kinetically hindered, and experimentally it is obtained by doping with KOH. Doping creates ionic defects that promote the migration of protons and the associated change in proton configuration. In this article, we mimic the effect of doping with a bias potential that enhances the formation of ionic defects in molecular dynamics simulations. The recombination of the ions thus formed proceeds through fast migration of the hydroxide along hydrogen bond loops, providing a physical and expedite way to change the proton configuration. Here, a key ingredient of this approach is a machine learning potential trained with density functional theory data and capable of modelling molecular dissociation. We exemplify the usefulness of this idea by studying the order-disorder transition using an appropriate order parameter that distinguishes the proton environments in ice Ih and XI. We calculate the changes in free energy, enthalpy, and entropy associated with the transition. Our estimated entropy agrees with experiment within the error bars of the calculation.
Nuclear–Electronic Orbital QM/MM Approach: Geometry Optimizations and Molecular Dynamics
Hybrid quantum mechanical/molecular mechanical (QM/MM) methods allow simulations of chemical reactions in atomistic solvent and heterogeneous environments such as proteins. Herein, the nuclear–electronic orbital (NEO) QM/MM approach is introduced to enable the quantization of specified nuclei, typically protons, in the QM region using a method such as NEO-density functional theory (NEO-DFT). This approach includes proton delocalization, polarization, anharmonicity, and zero-point energy in geometry optimizations and dynamics. Expressions for the energies and analytical gradients associated with the NEO-QM/MM method, as well as the previously developed polarizable continuum model (NEO-PCM), are provided. Geometry optimizations of small organic molecules hydrogen bonded to water in either dielectric continuum solvent or explicit atomistic solvent illustrate that aqueous solvation can strengthen hydrogen-bonding interactions for the systems studied, as indicated by shorter intermolecular distances at the hydrogen-bond interface. We then performed a real-time direct dynamics simulation of a phenol molecule in explicit water using the NEO-QM/MM method. Furthermore, these developments and initial examples provide the foundation for future studies of nuclear–electronic quantum dynamics in complex chemical and biological environments.
Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory Based Ab Initio Molecular Dynamics. 1. Theory, Algorithm, and Performance
By including a fraction of exact exchange (EXX), hybrid functionals reduce the self-interaction error in semilocal density functional theory (DFT) and thereby furnish a more accurate and reliable description of the underlying electronic structure in systems throughout biology, chemistry, physics, and materials science. However, the high computational cost associated with the evaluation of all required EXX quantities has limited the applicability of hybrid DFT in the treatment of large molecules and complex condensed-phase materials. To overcome this limitation, we describe a linear-scaling approach that utilizes a local representation of the occupied orbitals (e.g., maximally localized Wannier functions (MLWFs)) to exploit the sparsity in the real-space evaluation of the quantum mechanical exchange interaction in finite-gap systems. In this work, we present a detailed description of the theoretical and algorithmic advances required to perform MLWF-based ab initio molecular dynamics (AIMD) simulations of large-scale condensed-phase systems of interest at the hybrid DFT level. We focus our theoretical discussion on the integration of this approach into the framework of Car–Parrinello AIMD, and highlight the central role played by the MLWF-product potential (i.e., the solution of Poisson’s equation for each corresponding MLWF-product density) in the evaluation of the EXX energy and wave function forces. We then provide a comprehensive description of the exx algorithm implemented in the open-source Quantum ESPRESSO program, which employs a hybrid MPI/OpenMP parallelization scheme to efficiently utilize the high-performance computing (HPC) resources available on current- and next-generation supercomputer architectures. Furthermore, this is followed by a critical assessment of the accuracy and parallel performance (e.g., strong and weak scaling) of this approach when AIMD simulations of liquid water are performed in the canonical (NVT) ensemble. With access to HPC resources, we demonstrate that exx enables hybrid DFT-based AIMD simulations of condensed-phase systems containing 500–1000 atoms (e.g., (H₂O)₂₅₆) with a wall time cost that is comparable to that of semilocal DFT. In doing so, exx takes us one step closer to routinely performing AIMD simulations of complex and large-scale condensed-phase systems for sufficiently long time scales at the hybrid DFT level of theory.
Molecular Dynamics Free Energy Simulations Reveal the Mechanism for the Antiviral Resistance of the M66I HIV-1 Capsid Mutation
While drug resistance mutations can often be attributed to the loss of direct or solvent-mediated protein−ligand interactions in the drug-mutant complex, in this study we show that a resistance mutation for the picomolar HIV-1 capsid (CA)-targeting antiviral (GS-6207) is mainly due to the free energy cost of the drug-induced protein side chain reorganization in the mutant protein. Among several mutations, M66I causes the most suppression of the GS-6207 antiviral activity (up to ~84,000-fold), and only 83- and 68-fold reductions for PF74 and ZW-1261, respectively. To understand the molecular basis of this drug resistance, we conducted molecular dynamics free energy simulations to study the structures, energetics, and conformational free energy landscapes involved in the inhibitors binding at the interface of two CA monomers. To minimize the protein−ligand steric clash, the I66 side chain in the M66I−GS-6207 complex switches to a higher free energy conformation from the one adopted in the apo M66I. In contrast, the binding of GS-6207 to the wild-type CA does not lead to any significant M66 conformational change. Based on an analysis that decomposes the absolute binding free energy into contributions from two receptor conformational states, it appears that it is the free energy cost of side chain reorganization rather than the reduced protein−ligand interaction that is largely responsible for the drug resistance against GS-6207.
Impact of Cr and Co on 99Tc retention in magnetite: A combined study of ab initio molecular dynamics and experiments
This work explores the effect of co-mingled dopants, Co(II) and Cr(III), on Tc(IV) incorporation and retention in magnetite when heat treated to 625 or 700 °C. Key trends in Tc retention in the high temperature regime were identified using a combination of density-functional-theory based ab initio molecular dynamics (AIMD) simulations, and batch experiments including solid phase characterization techniques, e.g. X-ray absorption spectroscopy. A stabilizing effect on Tc(IV) was observed when the number of Tc and Cr atoms are equal or when the magnetite surface is oversaturated with Tc and Cr inclusions. Here, oversaturation is hypothesized to force Cr from the magnetite surface to form a Cr2O3 phase, which may act as a protective layer that prevents Tc release. With the addition of Co, Tc(IV) is stabilized via redox processes. The presence of Cr in low concentrations interferes with this redox stabilization and Cr is preferentially stabilized as opposed to Tc. As a result, using Co as a stabilizing dopant for Tc may be compromised in the presence of Cr. When the relative concentration of Tc, Cr and Co is the same, or more Co atoms are added to high Cr incorporated systems, the formation of the Cr2O3 phase may be suppressed. Although waste streams with co-mingled Tc and Cr potentially may benefit from a Co dopant, since the formation of a Cr2O3 passivation layer may protect incorporated Tc from being released, the relative concentration of the three elements will be a critical parameter for maximizing effectiveness of this strategy.
Ab initio Molecular Dynamics Assessment of Thermodynamic and Transport Properties in (K,Li)Cl and (K, Na)Cl Molten Salt Mixtures
Molten salt mixtures are integral part of highly important technological applications such as nuclear reactors. However, due to inherent difficulties associated with experiment at high temperatures and the intrinsic complexity of liquid-phase multi-component systems, understanding their properties at a molecular level remains a challenge. Here, we report on an ab initio molecular dynamics investigation on structural, electronic, transport, and thermal properties of two common molten salt mixtures, (K, Li)Cl and (K,Na)Cl, at five different compositions and three temperatures. Most of the properties were found to depend on both composition and temperature. While properties, like atomic charges, show additive behaviors, other properties, such as electrical conductivity, show considerable deviations from additivity. We shall show that the mixing of the molten salt mixtures is mainly driven by entropy, and that the KCl and LiCl mix better than KCl and NaCl. Our computational results are in general consistent with available experimental data. Comparison with available theoretical data is also provided.