Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “molecular dynamic (MD) 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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Simulations with machine learning potentials identify the ion conduction mechanism mediating non-Arrhenius behavior in LGPS

Abstract Li 10 Ge(PS 6 ) 2 (LGPS) is a highly concentrated solid electrolyte, in which Coulombic repulsion between neighboring cations is hypothesized as the underlying reason for concerted ion hopping, a mechanism common among superionic conductors such as Li 7 La 3 Zr 2 O 12 (LLZO) and Li 1.3 Al 0.3 Ti 1.7 (PO 4 ) 3 (LATP). While first principles simulations using molecular dynamics (MD) provide insight into the Li + transport mechanism, historically, there has been a gap in the temperature ranges studied in simulations and experiments. Here, we used a neural network potential trained on density functional theory (DFT) simulations, to run up to 40-nanosecond long MD simulations at DFT-like accuracy to characterize the ion conduction mechanisms across a range of temperatures that includes previous simulations and experimental studies. We have confirmed a Li + sublattice phase transition in LGPS around 400 K, below which the ab -plane diffusivity D a b ∗ is drastically reduced. Concomitant with the sublattice phase transition near 400 K, there is less cation-cation (cross) correlation, as characterized by Haven ratios closer to 1, and the vibrations in the system are more harmonic at lower temperature. Intuitively, at high temperature, the collection of vibrational modes may be sufficient to drive concerted ion hops. However, near room temperature, the vibrational modes available may be insufficient to overcome electrostatic repulsion, thus resulting in less correlated ion motion and comparatively slower ion conduction. Such phenomena of a sublattice phase transition, below which concerted hopping plays a less significant role, may be extended to other highly concentrated solid electrolytes such as LLZO and LATP.

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

Molecular dynamics simulation of the effect of cooling rate on the structure and properties of lithium disilicate glass

The effects of cooling rate on the structure and properties of lithium disilicate (LS2) glass are investigated using molecular dynamics (MD) computer simulations. The evolution of structural features such as pair distribution function, bond angle distribution, and Li coordination distribution are determined, and correlated with dynamic and static properties to elucidate the effects of cooling rate. The density, elastic moduli, and diffusion coefficient are found to be highly sensitive to cooling rate, whereas the Si pair distribution function and bond angle distribution are weakly affected by the cooling rate. Additionally, the changes of Si-O-Si bond angle and Li coordination number suggest the formation of Li cluster at lower cooling rates. Furthermore, by comparing results from other simulations and reported experiments, we confirm that the increase of cooling rate leads to an increase of conductivity and a decrease of density. Finally, at very high cooling rates, we find that all the atoms in LS2-glass do not relax simultaneously, but do so in two distinct configurations.

36 MATERIALS SCIENCE↗
NASA NTRS2018

Simple Ion Channels: From Structure to Electrophysiology and Back

A reliable way to establish whether our understanding of a channel is satisfactory is to reproduce its measured ionic conductance over a broad range of applied voltages in computer simulations. In molecular dynamics (MD), this can be done by way of applying an external electric field to the system and counting the number of ions that traverse the channel per unit time. Since this approach is computationally very expensive, we have developed a markedly more efficient alternative in which MD is combined with the electrodiffusion (ED) equation. In this approach, the assumptions of the ED equation can be rigorously tested, and the precision and consistency of the calculated conductance can be determined. We have demonstrated that the full current/voltage dependence and the underlying free energy profile for a simple channel can be reliably calculated from equilibrium or non-equilibrium MD simulations at a single voltage. To carry out MD simulations, a structural model of a channel has to be assumed, which is an important constraint, considering that high-resolution structures are available for only very few simple channels. If the comparison of calculated ionic conductance with electrophysiological data is satisfactory, it greatly increases our confidence that the structure and the function are described sufficiently accurately. We examined the validity of the ED for several channels embedded in phospholipid membranes - four naturally occurring channels: trichotoxin, alamethicin, p7 from hepatitis C virus (HCV) and Vpu from the HIV-1 virus, and a synthetic, hexameric channel, formed by a 21-residue peptide that contains only leucine and serine. All these channels mediate transport of potassium and chloride ions. It was found that the ED equation is satisfactory for these systems. In some of them experimental and calculated electrophysiological properties are in good agreement, whereas in others there are strong indications that the structural models are incorrect.

molecular dynamics (MD↗
DOE OSTI2022

Modular Approach for the Synthesis of Bottlebrush Diblock Copolymers from Poly(Glycidyl Methacrylate)-block-Poly(Vinyldimethylazlactone) Backbones

A strategy that enables the facile synthesis of bottlebrush block copolymers with flexible backbones was developed. A demonstration of the strategy’s utility was carried out by grafting end-functionalized polymethylmethacrylate (PMMA) and polystyrene (PS) to the dually reactive block copolymer, poly(glycidyl methacrylate)-block-poly(vinyldimethylazlactone) (PGMA-b-PVDMA). Five different bottlebrush diblock copolymers were investigated by size-exclusion chromatography (SEC), 1H NMR, Fourier transform infrared (FT-IR), differential scanning calorimetry (DSC), X-ray scattering methods, atomic force microscopy (AFM), rheology and computational simulations using molecular dynamics (MD), and self-consistent field theory (SCFT). A relationship between the glass transition temperature and the fraction of chain ends was demonstrated by DSC and highlights the potential of this synthetic method to tailor polymer properties. All five samples were found to be in a disordered phase exhibiting multiscale structures revealed by two broad peaks in small-angle X-ray scattering (SAXS) that can be attributed to graft-to-graft and backbone-to-backbone density correlations using MD simulations. The SCFT-based simulations justify the observation of a disordered phase due to its stabilization by the grafts. Additionally, this modular approach can be easily extended to other grafts, including responsive, conducting, and charged polymers with the prerequisite end groups. The versatility and ease of assembling these functional bottlebrushes constitute a powerful “toolbox” method for the rapid and scalable synthesis of novel bottlebrush block copolymers with desired properties.

36 MATERIALS SCIENCE↗
NASA NTRS2020

Ion Channels – From Structure to Electrophysiology and Back

A reliable way to establish whether our understanding of a channel is satisfactory is to reproduce its measured ionic conductance over a broad range of applied voltages in computer simulations. In molecular dynamics (MD), this can be done by way of applying an external electric field to the system and counting the number of ions that traverse the channel per unit time. Since this approach is computationally very expensive, we have developed a markedly more efficient alternative in which MD is combined with the electrodiffusion (ED) equation. In this approach, the assumptions of the ED equation can be rigorously tested, and the precision and consistency of the calculated conductance can be determined. We have demonstrated that the full current/voltage dependence and the underlying free energy profile for a simple channel can be reliably calculated from equilibrium or non-equilibrium MD simulations at a single voltage. Free energy profiles can be obtained from non-equilibrium simulations without a loss of accuracy even without the knowledge of diffusion coefficient. To carry out MD simulations, a structural model of a channel has to be assumed, which is an important constraint, considering that high-resolution structures are available for only very few simple channels. If the comparison of calculated ionic conductance with electrophysiological data is satisfactory, it greatly increases our confidence that the structure and the function are described sufficiently accurately. We examined the validity of the ED for several channels embedded in phospholipid membranes - four naturally occurring channels: trichotoxin, alamethicin, p7 from hepatitis C virus (HCV) and Vpu from the HIV-1 virus, a synthetic, hexameric channel, formed by a 21-residue peptide that contains only leucine and serine and a bacterial pentameric ligand-gated ion channel, GLIC. All these channels mediate transport of potassium and chloride ions. It was found that the ED equation is satisfactory for these systems. In some of them, experimental and calculated electrophysiological properties are in good agreement, whereas in others there are strong indications that the structural models are incorrect.

Pohorille, Andrew↗
DOE OSTI2020

Role of Capping Agents in the Synthesis of Salicylate-Capped Zinc Oxide Nanoparticles

Capping agents are often used for controlling the size, aggregation, and properties of nanoparticles. To guide the design of improved nanomaterials for targeted performance, one can use mechanistic insights into the interactions between capping agents and nanoparticles. Here, we employ density functional theory (DFT), reactive force-field molecular dynamics (ReaxFF MD) simulations, and optical spectroscopy to study the interactions between salicylate, as a model capping agent, and zinc oxide (ZnO) nanoparticles. In this work, we find that salicylate strongly interacts with the nanoparticle via the formation of O–Zn bonds in a distorted six-membered coordination ring structure. We describe the mechanisms of capping of ZnO nanoparticles by salicylate via three different binding modes. Simulations indicate that salicylate undergoes dissociative adsorption at the highly active surface Zn sites via a hydrogen-transfer process, thereby forming a tridentate configuration. The water-mediated interaction also facilitates the dissociative adsorption, leading to two salicylate O atoms coordinating with a surface Zn atom, while the other salicylate O atom bonds with another surface Zn atom. For molecular adsorption, binding free energies indicate that salicylate binds more strongly to ZnO, often in a bidentate configuration, than water does. The formation of the salicylate–ZnO complex is substantiated by UV–visible and Fourier transform infrared spectra. We find that the C=O stretching mode of salicylate becomes softened when it interacts with the nanoparticle, suggesting chemisorption of salicylate on ZnO. Although DFT predicts strong interaction between salicylate and ZnO, ReaxFF MD simulation indicates the moderate interaction between these two components in aqueous solution. Water molecules in close contact with the nanoparticle surface undergo dissociation, thus resulting in a surface hydroxyl, a reactive oxygen species that may influence the nanoparticle’s catalytic properties. Overall, the atomic-level information provided here can guide the selection process of salicylate as an appropriate agent for ZnO nanoparticle synthesis when strong interactions of particles with capping agents are required.

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

Atomistic Modeling of Irradiation-Induced Defects and Clusters in Additively-Manufactured Austenitic Stainless Steel

Understanding the progression of irradiation-induced defects in real materials and alloys requires insight into the interaction between point defects, defect clusters, and the existing microstructure. This is particularly true when it comes to additively manufactured (AM) materials with complicated microstructures and residual stresses. In this work, we employ molecular dynamics (MD) to simulate the creation of point defects via primary radiation damage modeling and point defect interactions with extended defects such as voids and dislocations using simplified but representative alloy compositions. We analyze the effect of alloying elements on the formation energies of the defect clusters and measure interactions such as sink strengths and binding energies. This study provides insight into the unique behaviors of irradiation damage in AM materials and information necessary to bridge the gap in length and time scales to connect the atomistic behavior of defects to longer length scale simulations such as cluster dynamics.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Accelerated Modeling of Lithium Diffusion in Solid State Electrolytes using Artificial Neural Networks

Abstract Previous efforts to understand structure‐function relationships in high ionic conductivity materials for solid state batteries have predominantly relied on density functional theory (DFT‐) based ab initio molecular dynamics (MD). Such simulations, however, are computationally demanding and cannot be reasonably applied to large systems containing more than a hundred atoms. Here, an artificial neural network (ANN) is trained to accelerate the calculation of high accuracy atomic forces and energies used during such MD simulations. After carefully training a robust ANN for four and five element systems, nearly identical lithium ion diffusivities are obtained for Li 10 GeP 2 S 12 (LGPS) when benchmarking the ANN‐MD results with DFT‐MD. Applying the ANN‐MD approach, the effect of chlorine doping on the lithium diffusivity is calculated in an LGPS‐like structure and it is found that a dopant concentration of 1.3% maximizes ionic conductivity. The optimal concentration balances the competing consequences of effective atomic radii and dielectric constants on lithium diffusion and agrees with the experimental composition. Performing simulations at the resolution necessary to model experimentally relevant and optimal concentrations would be infeasible with traditional DFT‐MD. Systems that require a large number of simulated atoms can be studied more efficiently while maintaining high accuracy with the proposed ANN‐MD framework.

Rao, Karun K.↗
DOE OSTI2021

Mechanistic study of pH effect on organic solvent nanofiltration using carboxylated covalent organic framework as a modeling and experimental platform

Covalent organic framework (COF) membranes for organic solvent nanofiltration (OSN) have attracted growing interest due to their high stability, well-defined pore size, and controllable pore functionality. In aqueous nanofiltration processes, solution pH has a significant effect on membrane permeability, selectivity, and fouling. However, our understanding of pH effects on OSN processes is still limited. In this work, we present experimental and reactive force field molecular dynamics (ReaxFF MD) simulation results on the effect of pH on OSN using a well-characterized carboxylated COF (C-COF) membrane. Filtration experimental and modeling results both demonstrate that pH changes can affect solvent permeance and solute rejection dramatically. Specifically, methanol permeance decreases with the increase of NaOH and HCl concentration, respectively. When Alcian Blue dye was used as the solute in methanol, the selectivity increased from 23% to 98% when changing the feed pH from 2.2 to 10.1, respectively. While ReaxFF MD results show that C-COF solvated pore size is only marginally affected by the pH of solution, methanol diffusion is affected by the size of methanol-solvated Na + and Cl - ions, which increase with increasing Na + and Cl - concentrations. Solution pH change also affected solute solvation layer thicknesses and C-COF pore charges, contributing to the significant AB dye selectivity increase with pH increase. To summarize, we have shown that dramatic changes in C-COF selectivity are not anomalous in OSN, but are instead controllable by varying feed solution pH to modify solvated solute radii and surface charge.

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

Direct Modulators of K-Ras–Membrane Interactions

Protein–membrane interactions (PMIs) are ubiquitous in cellular signaling. Initial steps of signal transduction cascades often rely on transient and dynamic interactions with the inner plasma membrane leaflet to populate and regulate signaling hotspots. Methods to target and modulate these interactions could yield attractive tool compounds and drug candidates. Here, we demonstrate that the conjugation of a medium-chain lipid tail to the covalent K-Ras(G12C) binder MRTX849 at a solvent-exposed site enables such direct modulation of PMIs. The conjugated lipid tail interacts with the tethered membrane and changes the relative membrane orientation and conformation of K-Ras(G12C), as shown by molecular dynamics (MD) simulation-supported NMR studies. In cells, this PMI modulation restricts the lateral mobility of K-Ras(G12C) and disrupts nanoclusters. The described strategy could be broadly applicable to selectively modulate transient PMIs.

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

Integrating Ultra-Coarse-Grained Protein Models into Accessible Workflows for Multiscale Molecular Dynamics

To capture protein conformational transitions using molecular dynamics (MD), several simulation resolutions covering different spatial and temporal scales are typically needed. All-atom (AA) simulations provide fine resolution, but are computationally infeasible for large systems over longer durations. Coarse-grained (CG) and ultra-coarse-grained (UCG) models have a lower resolution and computational cost while still being able to conserve essential protein features. Prior work on a Multiscale Machinelearned Modeling Infrastructure (MuMMI) combined both AA and CG simulations to study RAS-RAF protein interactions, leveraging CG models for longer time scales and using AA to investigate unusual conformations in greater detail. However, MuMMI is still resource-intensive, and this study aims to maximize exploration of the protein conformational space while reducing computational cost. In this paper, we build on prior work that integrates UCG models based on heterogeneous elastic network modeling (hENM) into the MuMMI workflow. We demonstrate that UCG models enable accurate sampling of protein conformations, focusing on simulating RAS-RAF protein interactions. Using higher-resolution CG Martini simulation data, we can automatically refine intramolecular interactions in UCG models. We present a scalable Python package that uses fluctuations observed in higher-resolution CG Martini simulations to estimate bond coefficients of the UCG model. We built novel machine learning-based backmapping methods to recover more detailed CG Martini structures from UCG structures, using diffusion models to learn the mapping between scales. Finally, we present UCG-mini-MuMMI, an accessible and less compute-intensive version of MuMMI as a resource for the scientific community. Incorporating UCG models into MD studies is applicable to a broad range of systems and proteins, and our study offers insights into the advantages and limitations of these methods.

Chemical structure↗
DOE OSTI2020

Radiation Damage in Materials—Helium Effects

Despite its scarcity in terrestrial life, helium effects on microstructure evolution and thermo-mechanical properties can have a significant impact on the operation and lifetime of applications, including: advanced structural steels in fast fission reactors, plasma facing and structural materials in fusion devices, spallation neutron target designs, energetic alpha emissions in actinides, helium precipitation in tritium-containing materials, and nuclear waste materials. The small size of a helium atom combined with its near insolubility in almost every solid makes the helium–solid interaction extremely complex over multiple length and time scales. This Special Issue, “Radiation Damage in Materials—Helium Effects”, contains review articles and full-length papers on new irradiation material research activities and novel material ideas using experimental and/or modeling approaches. These studies elucidate the interactions of helium with various extreme environments and tailored nanostructures, as well as their impact on microstructural evolution and material properties.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Quantum bath augmented stochastic nonequilibrium atomistic simulations for molecular heat conduction

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

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

Parallel Grand-Canonical Monte Carlo (ParaGrandMC) User’s Manual Version 2.0

This manual describes the commands and command line options for the Parallel Grand Canonical Monte Carlo version 2.0 (ParaGrandMC.2.0) simulation code. This is a highly scalable parallel FORTRAN 2003 code for simulating the thermodynamic evolution of materials at the atomic level, and predicting their thermodynamic state, phase diagram, chemical composition and mechanical properties. The code is specifically designed to simulate multi-component alloy systems, predict solid-state phase transformations such as austenite-martensite transformations, precipitate formation, recrystallization, capillary effects at interfaces, surface absorption, etc., which can aid the design of novel metallic alloys. While the software is mainly tailored for modeling metal alloys, it can also be used for other types of solid-state systems, and to some degree for liquid or gaseous systems, including multiphase systems forming solid-liquid-gas interfaces. In addition to performing Monte Carlo (MC) simulations, the code can also perform Molecular Dynamics (MD) and Langevin Dynamics (LD) simulations, which can be combined and interchanged with MC for faster and more efficient system evolution. A detailed description of the MC part of the code is provided in the NASA ParaGrandMC report: NASA/CR–2016-219202; http://www.sti.nasa.gov.

High performance computing↗
DOE OSTI2025

An atomic cluster expansion (ACE) potential for water under extreme conditions

We present a machine learning interatomic potential for water designed to capture its complex multiphase behavior, including both molecular and superionic ice phases. The potential is based on the atomic cluster expansion (ACE) formulation and has been parameterized to enable high-fidelity molecular dynamics simulations of water under extreme conditions, for pressures up to 100 GPa and for temperatures between 500 and 6000 K. A diverse range of configurations was generated through ab initio molecular dynamics (AI-MD) simulations, covering insulating and superionic ice phases, liquid water, and dissociated plasma phase. We demonstrate that the H 2 O ACE potential accurately reproduces experimental and DFT predicted isotherms and Hugoniots. Crucially, the potential is able to capture the intricate phase behavior of water, including the transition from molecular fluid to the appropriate solid ice phases, and the superionic ice phases. This work provides a robust interatomic potential that can be used for large-scale, accurate simulations of water under extreme thermodynamic conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Nonadiabatic molecular dynamics simulation of C2H22+ in a strong laser field*

We investigate the alignment dependence of the strong laser dissociation dynamics of molecule C 2 H 2 2 + in the frame of real-time and real-space time-dependent density function theory coupled with nonadiabatic quantum molecular dynamics (TDDFT-MD) simulation. This work is based on a recent experiment study “ultrafast electron diffraction imaging of bond breaking in di-ionized acetylene” [Wolter et al , Science 354 , 308–312 (2016)]. Our simulations are in excellent agreement with the experimental data and the analysis confirms that the alignment dependence of the proton dissociation dynamics comes from the electron response of the driving laser pulse. Our results validate the ability of the TDDFT-MD method to reveal the underlying mechanism of experimentally observed and control molecular dissociation dynamics.

Physics↗
DOE OSTI2020

Calculation of the detonation state of HN 3 with quantum accuracy

HN 3 is a unique liquid energetic material that exhibits ultrafast detonation chemistry and a transition to metallic states during detonation. We combine the Chebyshev interaction model for efficient simulation (ChIMES) many-body reactive force field and the extended-Lagrangian multiscale shock technique molecular dynamics method to calculate the detonation properties of HN 3 with the accuracy of Kohn–Sham density-functional theory. ChIMES is based on a Chebyshev polynomial expansion and can accurately reproduce density-functional theory molecular dynamics (DFT-MD) simulations for a wide range of unreactive and decomposition conditions of liquid HN 3 . We show that addition of random displacement configurations and the energies of gas-phase equilibrium products in the training set allows ChIMES to efficiently explore the complex potential energy surface. Schemes for selecting force field parameters and the inclusion of stress tensor and energy data in the training set are examined. Structural and dynamical properties and chemistry predictions for the resulting models are benchmarked against DFT-MD. We demonstrate that the inclusion of explicit four-body energy terms is necessary to capture the potential energy surface across a wide range of conditions. Our results generally retain the accuracy of DFT-MD while yielding a high degree of computational efficiency, allowing simulations to approach orders of magnitude larger time and spatial scales. The techniques and recipes for MD model creation we present allow for direct simulation of nanosecond shock compression experiments and calculation of the detonation properties of materials with the accuracy of Kohn–Sham density-functional theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗