Search NASA⌕ Search

SEARCH · Search NASA

Results for “Molecular dynamics computer 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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Atomic Diffusivities of Yttrium, Titanium and Oxygen Calculated by Ab Initio Molecular Dynamics in Molten 316L Oxide-Dispersion-Strengthened Steel Fabricated via Additive Manufacturing

Oxide-dispersion-strengthened (ODS) steels have long been viewed as a prime solution for harsh environments. However, conventional manufacturing of ODS steels limits the final product geometry, is difficult to scale up to large components, and is expensive due to multiple highly involved, solid-state processing steps required. Additive manufacturing (AM) can directly incorporate dispersion elements (e.g., Y, Ti and O) during component fabrication, thus bypassing the need for an ODS steel supply chain, the scale-up challenges of powder processing routes, the buoyancy challenges associated with casting ODS steels, and the joining issues for net-shape component fabrication. In the AM process, the diffusion of the dispersion elements in the molten steel plays a key role in the precipitation of the oxide particles, thereby influencing the microstructure, thermal stability and high-temperature mechanical properties of the resulting ODS steels. In this work, the atomic diffusivities of Y, Ti, and O in molten 316L stainless steel (SS) as functions of temperature are determined by ab initio molecular dynamics simulations. The latest Vienna Ab initio Simulation Package (VASP) package that incorporates an on-the-fly machine learning force field for accelerated computation is used. At a constant temperature, the time-dependent coordinates of the target atoms in the molten 316L SS were analyzed in the form of mean square displacement in order to obtain diffusivity. The values of the diffusivity at multiple temperatures are then fitted to the Arrhenius form to determine the activation energy and the pre-exponential factor. Given the challenges in experimental measurement of atomic diffusivity at such high temperatures and correspondingly the lack of experimental data, this study provides important physical parameters for future modeling of the oxide precipitation kinetics during AM process.

Chemistry↗
DOE OSTI2023

Non-Empirical and Self-Interaction Corrections for DFTB: Towards Accurate Quantum Simulations for Large Mesoscale Systems (Final Report)

This project was comprised of two complementary (but parallel) thrusts: (1) implementing massively-parallelized computing hardware (with new computational hardware that may replace GPUs) in the density functional tight binding (DFTB) approach and (2) developing new capabilities in DFTB to calculate the electronic structure and dynamics of large chemical systems. While classical molecular dynamics can handle hundreds of thousands of atoms, it cannot provide a first-principles based description of chemical systems at the quantum level. At the other extreme, conventional Kohn-Sham DFT methods can probe the true quantum mechanical nature of chemical systems; however, these methods cannot tackle the large sizes and length scales relevant to dynamics simulations of realistic systems. The DFTB formalism utilized in this project provides a viable approach for probing these large systems at a quantum mechanical level of detail. However, to utilize the DFTB approach for accurate calculations of electronic properties, it is crucial to incorporate quantum-based corrections in DFTB since exchange-correlation effects can still remain very strong in these large systems. At the same time, enhancing the computational efficiency of DFTB is also essential since optimal computational performance is required for addressing the large size scales associated with realistic chemical systems. As such, the new non-empirical corrections and computing hardware enhancements implemented in this project will enable accurate and computationally efficient approaches to directly probe electronic properties in these large, complex systems.

36 MATERIALS SCIENCE↗
DOE OSTI2022

GLUE Code: A framework handling communication and interfaces between scales

Many scientific applications are inherently multiscale in nature. Such complex physical phenomena often require simultaneous execution and coordination of simulations spanning multiple time and length scales. This is possible by combining expensive small-scale simulations (such as molecular dynamics simulations) with larger scale simulations (such continuum limit/hydro solvers) to allow for considerably larger systems using task and data parallelism. However, the granularity of the tasks can be very large and often leads to load imbalance. Traditionally, we use approximations to streamline the computation of the more costly interactions and this introduces trade-offs between simulation cost and accuracy. In recent years, the available computational power and the advances in machine learning have made computing these scale-bridging interactions and multiscale simulations more feasible. One driving application has been plasma modeling in inertial confinement fusion (ICF), which is fundamentally multiscale in nature. This requires deep understanding of how to extrapolate microscopic information into macroscopically relevant scales. For example, in ICF one needs an accurate understanding of the connection between experimental observables and the underlying microphysics. The properties of the larger scales are often affected by the microscale behavior incorporated usually into the equations of state and ionic and electronic transport coefficients (Liboff, 1959; Rinderknecht et al., 2014; Rosenberg et al., 2015; Ross et al., 2017). Instead of incorporating this information using reliable molecular dynamics (MD) simulations, one often needs to use theoretical models, due to the inability of MD to reach engineering scales (Glosli et al., 2007; Marinak et al., 1998). One approach to resolve this issue is by coupling two MD simulations of different scales via force interpolation, e.g., the AdResS method (Krekeler et al., 2018; Nagarajan et al., 2013). Another approach, which we will pursue in the scope of this work, is by enabling scale bridging between MD simulations and meso/macro-scale models through the development and support of application programming interfaces that these different applications can interact through.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2004

Unfolding of Proteins: Thermal and Mechanical Unfolding

We have employed a Hamiltonian model based on a self-consistent Gaussian appoximation to examine the unfolding process of proteins in external - both mechanical and thermal - force elds. The motivation was to investigate the unfolding pathways of proteins by including only the essence of the important interactions of the native-state topology. Furthermore, if such a model can indeed correctly predict the physics of protein unfolding, it can complement more computationally expensive simulations and theoretical work. The self-consistent Gaussian approximation by Micheletti et al. has been incorporated in our model to make the model mathematically tractable by signi cantly reducing the computational cost. All thermodynamic properties and pair contact probabilities are calculated by simply evaluating the values of a series of Incomplete Gamma functions in an iterative manner. We have compared our results to previous molecular dynamics simulation and experimental data for the mechanical unfolding of the giant muscle protein Titin (1TIT). Our model, especially in light of its simplicity and excellent agreement with experiment and simulation, demonstrates the basic physical elements necessary to capture the mechanism of protein unfolding in an external force field.

Hur, Joe S.↗
DOE OSTI2021

Thermal gradient effect on helium and self-interstitial transport in tungsten

First-wall materials in a fusion reactor are expected to withstand harsh conditions, with high heat and particle fluxes that modify the materials microstructure. These fluxes will create strong gradients of temperature and concentration of diverse species. Besides the He ash and the hydrogenic species, neutron particles generated in the fusion reaction will collide with the material creating intrinsic defects, such as vacancies, self-interstitials atoms (SIAs), and clusters of such point defects. These defects and the He atoms will then migrate in the presence of the aforementioned gradients. In this study, we use nonequilibrium molecular dynamics to analyze the transport of He and SIAs in the presence of a thermal gradient in tungsten. We observe that, in all cases, the defects and impurity atoms tend to migrate toward the hot regions of the tungsten sample. The resulting species concentration profiles are exponential distributions, rising toward the hot regions of the sample, in agreement with irreversible thermodynamics analysis. For both He atoms and SIAs, we find that the resulting species flux is directed opposite to the heat flux, indicating that species transport is governed by a Soret effect (thermal-gradient-driven diffusion) characterized by a negative heat of transport that drives species diffusion uphill (from the cooler to the hot regions of the sample). Here, we demonstrate that the steady-state species profiles obtained accounting for the Soret effect vary significantly from those where temperature-gradient-driven transport is not considered and discuss the implications of such a Soret effect on the response to plasma exposure of plasma-facing tungsten.

36 MATERIALS SCIENCE↗
DOE OSTI2024

PigmentHunter: A point-and-click application for automated chlorophyll-protein simulations

Chlorophyll proteins (CPs) are the workhorses of biological photosynthesis, working together to absorb solar energy, transfer it to chemically active reaction centers, and control the charge-separation process that drives its storage as chemical energy. Yet predicting CP optical and electronic properties remains a serious challenge, driven by the computational difficulty of treating large, electronically coupled molecular pigments embedded in a dynamically structured protein environment. To address this challenge, we introduce here an analysis tool called PigmentHunter, which automates the process of preparing CP structures for molecular dynamics (MD), running short MD simulations on the nanoHUB.org science gateway, and then using electrostatic and steric analysis routines to predict optical absorption, fluorescence, and circular dichroism spectra within a Frenkel exciton model. Inter-pigment couplings are evaluated using point-dipole or transition-charge coupling models, while site energies can be estimated using both electrostatic and ring-deformation approaches. The package is built in a Jupyter Notebook environment, with a point-and-click interface that can be used either to manually prepare individual structures or to batch-process many structures at once. Here, we illustrate PigmentHunter’s capabilities with example simulations on spectral line shapes in the light harvesting 2 complex, site energies in the Fenna–Matthews–Olson protein, and ring deformation in photosystems I and II.

14 SOLAR ENERGY↗
DOE OSTI2023

Uncertainty-aware molecular dynamics from Bayesian active learning for phase transformations and thermal transport in SiC

Abstract Machine learning interatomic force fields are promising for combining high computational efficiency and accuracy in modeling quantum interactions and simulating atomistic dynamics. Active learning methods have been recently developed to train force fields efficiently and automatically. Among them, Bayesian active learning utilizes principled uncertainty quantification to make data acquisition decisions. In this work, we present a general Bayesian active learning workflow, where the force field is constructed from a sparse Gaussian process regression model based on atomic cluster expansion descriptors. To circumvent the high computational cost of the sparse Gaussian process uncertainty calculation, we formulate a high-performance approximate mapping of the uncertainty and demonstrate a speedup of several orders of magnitude. We demonstrate the autonomous active learning workflow by training a Bayesian force field model for silicon carbide (SiC) polymorphs in only a few days of computer time and show that pressure-induced phase transformations are accurately captured. The resulting model exhibits close agreement with both ab initio calculations and experimental measurements, and outperforms existing empirical models on vibrational and thermal properties. The active learning workflow readily generalizes to a wide range of material systems and accelerates their computational understanding.

36 MATERIALS SCIENCE↗
DOE OSTI2022

E(3)-equivariant graph neural networks for data-efficient and accurate interatomic potentials

Abstract This work presents Neural Equivariant Interatomic Potentials (NequIP), an E(3)-equivariant neural network approach for learning interatomic potentials from ab-initio calculations for molecular dynamics simulations. While most contemporary symmetry-aware models use invariant convolutions and only act on scalars, NequIP employs E(3)-equivariant convolutions for interactions of geometric tensors, resulting in a more information-rich and faithful representation of atomic environments. The method achieves state-of-the-art accuracy on a challenging and diverse set of molecules and materials while exhibiting remarkable data efficiency. NequIP outperforms existing models with up to three orders of magnitude fewer training data, challenging the widely held belief that deep neural networks require massive training sets. The high data efficiency of the method allows for the construction of accurate potentials using high-order quantum chemical level of theory as reference and enables high-fidelity molecular dynamics simulations over long time scales.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRSDate not supplied

Semi-Empirical Interatomic Potential for Large-Scale Molecular Dynamics Simulation of Metal-Oxide Systems

Classical molecular dynamics (MD) simulation can be applied to systems containing billions of atoms during times up to microseconds. However, utilization of a large-scale MD simulation requires reliable but computationally cheap interatomic potentials. In the case of metallic systems, embedded atom method (EAM) and Finnis-Sinclair (FS) potentials are the reasonable choices, but their development for multicomponent alloys is a challenge. Addition of oxygen atoms makes the problem of developing semi-empirical potentials even more difficult because the long-range Coulomb interaction cannot be neglected, and different atoms have different charges which vary during MD relaxation such that a charge equilibration procedure should be applied every MD step. As results researchers tend to avoid simulating metal-oxide systems. In this talk, a new Finnis-Sinclair (FS) Ni-Co-Cr potential will be presented. Special attention will be paid to reproducing of the element melting temperatures and clustering. Next, a charge transfer ionic potential (CTIP) developed to simulate the NiO properties will be presented. This potential utilizes the Ni FS potential and accounts for the Coulomb interaction in the vicinity of oxygen atoms. Fitting of the CTIP parameters and FS potential functions describing the interaction with oxygen atoms to the available experimental and ab initio data will be discussed. Results of MD simulation of interaction of NiO particles with dislocations in Ni will be shown.

molecular dynamics simulation↗
NASA NTRS2019

High-Performance Computing Optimization for Aladyn – Adaptive Neural Network Molecular Dynamics Mini-Application

This report provides a description and performance evaluation of the optimization techniques for high performance computing (HPC) implementation of the open source Computational Materials mini-application Aladyn (https://github.com/nasa/aladyn). Aladyn is a basic molecular dynamics code written in FORTRAN 2003, which is designed to demonstrate the use of adaptive neural networks (ANNs) in atomistic simulations. The role of ANNs is to efficiently reproduce the very complex energy landscape resulting from the atomic interactions in materials with the accuracy of the more expensive quantum mechanics-based calculations. The ANN is trained on a large set of atomic structures calculated using the density functional theory (DFT) method. While achieving orders of magnitude faster computational performance than DFT, the ANN-based approach was still very computationally demanding compared to the conventional approach of using empirically fitted energy functions. After its initial development, Aladyn was evaluated and optimized by experts at the NASA Advanced Supercomputing (NAS) division to exploit modern supercomputer architectures. The code has been optimized for execution on multicore central processing units (CPUs), including Intel® Skylake microarchitecture, and on graphic accelerators, such as Nvidia® V100 graphic processing units (GPUs), using Open Multi-Processing (OpenMP) and Open Accelerators (OpenACC) programming interfaces. The optimization achieved a speedup of 4.7 times the baseline version on CPU performance and an additional 2.4 times on CPU+GPU performance. Atomistic computer simulations are a fundamental tool in materials research to model material properties form physics-based first principles. Atomic interaction, governed by Quantum Mechanics (QM) require sophisticated and highly computationally demanding mathematical models to calculate [1]. Classical methods use approximate functional forms, empirically fitted through a set of variable parameters to emulate atomic energies as direct functions of atomic coordinates [2]. While empirical potentials are computationally much simpler, allowing simulations of large-scale systems of up to a trillion (1012) atoms [3], they are substantially less accurate compared to quantum calculations and applicable only to very specific atomic configurations or predefined crystallographic phases. A recently suggested approach is to use heuristic machine learning methods [4], such as those based on Adaptive Neural Networks (ANNs) to predict atomic energies, after being trained on a sufficiently large database of QM-calculated structures [5,6]. This approach reduces significantly the computational complexity, allowing for simulations of orders of magnitude larger systems compared to QM-based methods without compromising accuracy. Still, compared to classical methods using empirical energy functions, ANN methods remain two- to three orders of magnitude more computationally demanding. Hence, the computational cost of simulations, together with the need for extensive training of ANNs, still makes the practical implementation of ANN-based methods quite challenging. The purpose of the Aladyn mini-application software [7], available as open source at https://github.com/nasa/aladyn, is to be a testbed for exploring possible optimization strategies to develop highly scalable parallel algorithms for ANN-based atomistic simulations. Aladyn is aimed at utilizing the architecture of the high-end modern highperformance computing (HPC) hardware based on multicore central processing units (CPUs) equipped with graphic processing unit (GPU) accelerators. Specifically, the goal is to optimize the performance on a single HPC compute node, before implementing scaling to multi-node parallelization using message passing interface (MPI). At the same time, the open source code of Aladyn can serve as a training model for students and professors in academia.

Yamakov, Vesselin I.↗
DOE OSTI2022

Phase diagram of the TIP4P/Ice water model by enhanced sampling simulations

In this work, we studied the phase diagram for the TIP4P/Ice water model using enhanced sampling molecular dynamics simulations. Our approach is based on the calculation of ice–liquid free energy differences from biased coexistence simulations that reversibly sample the melting and growth of layers of ice. We computed a total of 19 melting points for five different ice polymorphs, which are in excellent agreement with the melting lines obtained from the integration of the Clausius–Clapeyron equation. For proton-ordered and fully proton-disordered ice phases, the results are in very good agreement with previous calculations based on thermodynamic integration. For the partially proton-disordered ice III, we find a large increase in stability that is in line with previous observations using direct coexistence simulations for the TIP4P/2005 model. This issue highlights the robustness of the approach employed here for ice polymorphs with diverse degrees of proton disorder. Our approach is general and can be applied to the calculation of other complex phase diagrams.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Nonadiabatic Molecular Dynamics Simulations Based on Time-Dependent Density Functional Tight-Binding Method

Nonadiabatic excited-state molecular dynamics underpin many photophysical and photochemical phenomena, such as exciton dynamics, charge separation and transport. In this work, we present an efficient nonadiabatic molecular dynamic (NAMD) simulation method based on time-dependent density functional tight-binding (TDDFTB) theory. Specifically, the adiabatic electronic structure, an essential NAMD input, is described at the TDDFTB level. Here, the nonadiabatic effects originating from the coupled motions of electrons and nuclei are treated by the trajectory surface hopping algorithm. To improve the computational efficiency, nonadiabatic couplings between excited states within the TDDFTB method are derived and implemented using an analytical approach. Further, the time-dependent nonadiabatic coupling scalars are calculated based on the overlap between molecular orbitals rather than the Slater determinants to speed up the simulations. In addition, the electronic decoherence scheme and a state reassigned unavoided crossings algorithm, which has been implemented in the NEXMD software, are used to improve the accuracy of the simulated dynamics and handle trivial unavoided crossings. Finally, the photoinduced nonadiabatic dynamics of a benzene molecule are simulated to demonstrate our implementation. The results for excited state NAMD simulations of benzene molecule based on TDDFTB method compare well that obtained with numerically expensive time-dependent density functional theory. The proposed methodology provides an attractive theoretical simulation tool for predicting the photophysical and photochemical properties of complex materials.

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

Reconciling Work Functions and Adsorption Enthalpies for Implicit Solvent Models: A Pt (111)/Water Interface Case Study

Implicit solvent models are a computationally efficient method of representing solid/liquid interfaces prevalent in electrocatalysis, energy storage, and materials science. However, electronic structure changes induced at the metallic surface by the dielectric continuum are not fully understood. To address this, we perform DFT calculations for the Pt(111)/water interface, in order to compare Poisson–Boltzmann continuum solvation methods with ab initio molecular dynamics (AIMD) simulations of explicit solvent. We show that the implicit solvent cavity can be parametrized in terms of the electric dipole moment change at the equilibrated explicit Pt/water interface to obtain the potential of zero charge (PZC). We also compare the accuracy of aqueous enthalpies of adsorption of phenol on Pt(111) using geometry and charge density based dielectric cavitation methods. The ability to parametrize the cavity according to individual atoms, as afforded in the geometry based approach, is key to obtaining accurate enthalpy changes of adsorption under aqueous conditions. Additionally, we show that the electronic structure changes induced by explicit solvent and our proposed implicit solvent parametrization scheme yield comparable density difference profiles and d-band projected density of states. We therefore demonstrate the capability of implicit solvent approaches to capture both the energetics of adsorption processes and the main electronic effects of aqueous solvent on the metallic surface. Therefore, this work provides a scheme for computationally efficient simulations of interfacial processes for applications in areas such as heterogeneous catalysis and electrochemistry.

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

Insights into the Properties of MXenes and MXene Analogs from Atomistic Simulation

We review the properties that have been predicted for MXenes and MXene analogs from computational simulation with a focus on structural and electronic properties, energy storage, and ion transport. Methods considered range from quantum mechanical approaches to classical molecular dynamics. We conclude by reviewing current limitations and outstanding questions for simulation and MXene properties that have been little explored to-date.

Muraleedharan, Murali Gopal↗
DOE OSTI2026

Elastic strain engineering of lattice thermal conductivity of silicon: An ab-initio study

Silicon (Si) is the most essential material in the semiconductor industry. It is important to manage the thermal properties of crystalline Si. Elastic strain engineering (ESE) has proven to be an effective tool in controlling the electrical conductivity of Si in strained-silicon technology; its effects on the thermal conductivity of silicon, therefore, warrants careful investigation. The ESE effect is much more pronounced for nanostructured materials due to the ultralarge elastic strains (on the order of 10%) achievable at the nanoscale. In this work, the lattice thermal conductivity (κ L ) of Si under hydrostatic, biaxial, and uniaxial strain states is studied with ab-initio simulations, and the values of strain-dependent κ L compare well with experimental results and existing molecular dynamics simulations. To understand the mechanisms of strain-modulated κ L , the phonon bands, scattering rate, and Grüneisen parameters of phonon modes are computed. It is shown that strain can significantly change the anharmonicity of the crystal system, thus changing phonon scattering rates and κ L . Our results demonstrate that ESE can reduce silicon κ L by up to approximately 90%. Furthermore, uniaxial and biaxial strains can induce highly anisotropic thermal conductivity in Si, with relative variations up to 58.5% and 14.5%, respectively.

Anisotropy in Thermal Conductivity↗
DOE OSTI2025

Atomistic Insights into Lithium–Glyme Solvate Ionic Liquids: Effects of Chain Length and Anion Coordination

For this study, mixtures of lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) in diglyme (G2), triglyme (G3) and tetraglyme (G4) at solvate ionic liquid (SIL) concentrations were investigated using classical molecular dynamics (cMD) simulations with a physically motivated force-field specifically developed for modeling these systems. The structural and dynamical properties of the mixtures were computed and analyzed. Lithium solvation shells, radial distribution functions, and X-ray structure factors were studied across the different SIL systems. Translational diffusion and rotational relaxation times were also evaluated, exhibiting similar trends with increasing glyme chain length. The results are consistent with experimental data and in good agreement with previous computational studies on G3 and G4. These findings validate the accuracy of the force field in modeling glyme systems and its use for describing the [Li(G2) 4/3 ][TFSI] mixture. Additionally, the thermal and electrochemical stability of these electrolytes were systematically examined. The thermal stability appears to be governed by cooperative interactions among glyme molecules, while the electrochemical stability is primarily influenced by Li + -anion interactions, which vary significantly with glyme chain length. Overall, the study sheds light on the crucial role of the anion in these glyme-based SILs and offers valuable insights into Li + -glyme systems at SIL concentrations, highlighting their promise as potential Li-ion battery electrolytes.

anions↗
DOE OSTI2023

Atomistic structure of the SARS-CoV-2 pseudoknot in solution from SAXS-driven molecular dynamics

SARS-CoV-2 depends on -1 programmed ribosomal frameshifting (-1 PRF) to express proteins essential for its replication. The RNA pseudoknot stimulating -1 PRF is thus an attractive drug target. However, the structural models of this pseudoknot obtained from cryo-EM and crystallography differ in some important features, leaving the pseudoknot structure unclear. We measured the solution structure of the pseudoknot using small-angle X-ray scattering (SAXS). The measured profile did not agree with profiles computed from the previously solved structures. Beginning with each of these solved structures, we used the SAXS data to direct all atom molecular dynamics (MD) simulations to improve the agreement in profiles. In all cases, this refinement resulted in a bent conformation that more closely resembled the cryo-EM structures than the crystal structure. Applying the same approach to a point mutant abolishing -1 PRF revealed a notably more bent structure with reoriented helices. This work clarifies the dynamic structures of the SARS-CoV-2 pseudoknot in solution.

60 APPLIED LIFE SCIENCES↗