Search NASA⌕ Search

SEARCH · Search NASA

Results for “MOLECULAR DYNAMICS METHOD”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Real-time GW -BSE investigations on spin-valley exciton dynamics in monolayer transition metal dichalcogenide

We develop an ab initio nonadiabatic molecular dynamics (NAMD) method based on GW plus real-time Bethe-Salpeter equation (GW + rtBSE-NAMD) for the spin-resolved exciton dynamics. From investigations on MoS 2 , we provide a comprehensive picture of spin-valley exciton dynamics where the electron-phonon (e-ph) scattering, spin-orbit interaction (SOI), and electron-hole (e-h) interactions come into play collectively. In particular, we provide a direct evidence that e-h exchange interaction plays a dominant role in the fast valley depolarization within a few picoseconds, which is in excellent agreement with experiments. Moreover, there are bright-to-dark exciton transitions induced by e-ph scattering and SOI. Our study proves that e-h many-body effects are essential to understand the spin-valley exciton dynamics in transition metal dichalcogenides and the newly developed GW + rtBSE-NAMD method provides a powerful tool for exciton dynamics in extended systems with time, space, momentum, energy, and spin resolution.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Thermal neutron scattering cross sections for amorphous carbon

Carbon materials are commonly found in both nuclear reactors and experimental systems. Various carbon structures occur in nuclear applications ranging from crystalline and nuclear graphite to the amorphous carbon seen in next-generation advanced reactor designs. Amorphous carbon is based on a randomized graphite-like structure and offers the unique ability to disperse impurities throughout the bulk composition. A graphite-like amorphous carbon system was modeled using the classical molecular dynamics (MD) code LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator). An improved version of the temperature-dependent Adaptive Intermolecular Reactive Empirical Bond Order (AIREBO) potential was used to model the carbon-carbon atomic interactions for the temperature at 300 K along with densities 1.60, 1.70, 1.85, and 2.23 g/cm{sup 3}. From the normalized velocity autocorrelation function (VACF), the phonon density of state (DOS) was then calculated as the Fourier transform of the normalized VACF. This DOS was then used as the primary input for the evaluation of the thermal scattering law (TSL, i.e. S(α,β)) and associated neutron thermal scattering cross sections. The TSL was analyzed using the Full Law Analysis Scattering System Hub (FLASSH). The amorphous structure results in shifts of the phonon DOS to lower energy modes than typically displayed for ideal crystalline graphite. This impact on the DOS is directly reflected in the TSL. Furthermore, the typical optical peak at 0.25 eV for the ideal graphite disappears for amorphous carbon, in good agreement with the expected structure. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Processes at lithium-hydride/deuteride surfaces upon low energy impact of H/D

Sputtering, reflection, and retention processes at amorphous and crystalline lithium hydride surfaces due to impact of low energy (1–100 eV) hydrogen and deuterium atoms over the range of 0 o −85 o angle of incidence at 300 K surface temperature were investigated by atomistic computational methods. Classical molecular dynamics simulations were performed with improved reactive bond-order force field (ReaxFF) potentials that include long-range polarization effects. In addition to probabilities of surface processes, the energy and angular spectra of ejected particles were obtained. Comparison of these results with those previously obtained on pristine lithium surfaces indicates the importance of saturation of the Li surface and near-surface region with hydrogen. We show that such saturation, which is typical in both laboratory and fusion device experiments with lithium coating of the plasma-facing surfaces, significantly changes the surface processes with hydrogen irradiation in the understudied low-energy region of impact energies.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Molecular Dynamics Simulation of Complex Reactivity with the Rapid Approach for Proton Transport and Other Reactions (RAPTOR) Software Package

Simulating chemically reactive phenomena such as proton transport on nanosecond to microsecond and beyond time scales is a challenging task. Ab initio methods are unable to currently access these time scales routinely, and traditional molecular dynamics methods feature fixed bonding arrangements that cannot account for changes in the system’s bonding topology. The Multiscale Reactive Molecular Dynamics (MS-RMD) method, as implemented in the Rapid Approach for Proton Transport and Other Reactions (RAPTOR) software package for the LAMMPS molecular dynamics code, offers a method to routinely sample longer time scale reactive simulation data with statistical precision. RAPTOR may also be interfaced with enhanced sampling methods to drive simulations toward the analysis of reactive rare events, and a number of collective variables (CVs) have been developed to facilitate this. Key advances to this methodology, including GPU acceleration efforts and novel CVs to model water wire formation are reviewed, along with recent applications of the method which demonstrate its versatility and robustness.

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

A Method for Obtaining Liquid–Solid Adsorption Rates from Molecular Dynamics Simulations: Applied to Methanol on Pt(111) in H 2 O

Adsorption is an important step in heterogeneous catalysis as it predetermines how many reactant molecules can participate in a surface reaction per unit time. While the rate of adsorption processes is well studied in gas–solid adsorption in both theory and experiment, such rates are still not well studied for liquid–solid adsorption. This is partly because the ever-changing configurations of liquid-phase solvent molecules impede the ability to study a molecule approaching a surface from a liquid phase by either experiment or theory. In this work, we develop a method using molecular dynamics (MD) simulations to study the rate of adsorption in liquid–solid adsorption processes. Specifically, we use MD to model the diffusion of a methanol molecule in aqueous solvent and its adsorption to a Pt(111) surface. We find that by approximating the solute motion as following the same displacement rates as a random walk model, the adsorbed and non-adsorbed states of the methanol molecule near the Pt(111) surface can be discerned and quantified. In particular, this methodology enables extracting a sticking coefficient and a macroscopically relatable adsorption rate. This method can be applied to arbitrary types of reactants and surfaces, as well as different liquid environments, thus providing a general tool for predicting quantitative adsorption rates of liquid–solid adsorption systems.

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

Fingerprinting Brownian Motions of Polymers under Flow

Here, we present a quantitative approach to the self-dynamics of polymers under steady flow by employing a set of complementary reference frames and extending the spherical harmonic expansion technique to dynamic density correlations. Application of this method to nonequilibrium molecular dynamics simulations of polymer melts reveals a number of universal features. For both unentangled and entangled melts, the center-of-mass motions in the flow frame are described by superdiffusive, anisotropic Gaussian distributions, whereas the isotropic component of monomer self-dynamics in the center-of-mass frame is strongly suppressed. Spatial correlation analysis shows that the heterogeneity of monomer self-dynamics increases significantly under flow.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Path Integrals for Nonadiabatic Dynamics: Multistate Ring Polymer Molecular Dynamics

This review focuses on a recent class of path-integral-based methods for the simulation of nonadiabatic dynamics in the condensed phase using only classical molecular dynamics trajectories in an extended phase space. Specifically, a semiclassical mapping protocol is used to derive an exact, continuous, Cartesian variable path-integral representation for the canonical partition function of a system in which multiple electronic states are coupled to nuclear degrees of freedom. Building on this exact statistical foundation, multistate ring polymer molecular dynamics methods are developed for the approximate calculation of real-time thermal correlation functions. As a result, the remarkable promise of these multistate ring polymer methods, their successful applications, and their limitations are discussed in detail.

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

An Implementation of Replica Exchange with Dynamical Scaling for Efficient Large-Scale Simulations

An implementation of the replica exchange with dynamical scaling (REDS) method in the commonly used molecular dynamics program GROMACS is presented. REDS is a replica exchange method that requires fewer replicas than conventional replica exchange while still providing data over a range of temperatures and can be used in either constant volume or constant pressure ensembles. Details for running REDS simulations are given, and an application to the human islet amyloid polypeptide (hIAPP) 11-25 fragment shows that the model efficiently samples conformational space.

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

Inertial dynamics of an interface with interfacial mass flux: Stability and flow fields’ structure, inertial stabilization mechanism, degeneracy of Landau’s solution, effect of energy fluctuations, and chemistry-induced instabilities

This work focuses on the long-standing problem of inertial dynamics of an interface with interfacial mass flux and reports new mechanisms for the interface stabilization and destabilization. The interface is a phase boundary separating fluids of different densities and having interfacial mass flux. To analyze the interface dynamics from a far field, we develop and apply the general matrix method to rigorously solve the boundary value problem involving the governing equations in the fluid bulk and the boundary conditions at the interface and at the outside boundaries of the domain. We find the fundamental solutions for the linearized system of equations and analyze the interplay of interface stability with flow fields’ structure by directly linking rigorous mathematical attributes to physical observables. We find that the interface is stable when the dynamics conserves the fluxes of mass, momentum, and energy; the stabilization is due to an inertial mechanism causing small oscillations of the interface velocity. In the classic Landau’s dynamics, the postulate of perfect constancy of the interface velocity leads to the development of Landau–Darrieus instability. This destabilization is also linked to the imbalance of the perturbed energy at the interface. The classic Landau’s solution is found to have degeneracy; lifting of the degeneracy may lead to singularity and self-similar dynamics. Our results compare well with traditional theories of combustion and propose new experiments to study the dynamics of the interface and the flow fields in combustible systems. We further conduct reactive molecular dynamics simulations to elucidate the complexity of chemical processes, to study the destabilizing effect of energy fluctuations on the interface stability, and to illustrate the chemistry-induced instabilities. In summary, we identify the extreme sensitivity of the interface dynamics to the interfacial boundary conditions, including the formal properties of fundamental solutions and the qualitative and quantitative properties of the flow fields. This provides new opportunities for studies, diagnostics, and control of multiphase flows in a broad range of processes in nature and technology.

42 ENGINEERING↗
DOE OSTI2024

Modern chemical graph theory

Abstract Graph theory has a long history in chemistry. Yet as the breadth and variety of chemical data is rapidly changing, so too do graph encoding methods and analyses that yield qualitative and quantitative insights. Using illustrative cases within a basic mathematical framework, we showcase modern chemical graph theory's utility in Chemists' analysis and model development toolkit. The encoding of both experimental and simulation data is discussed at various levels of granularity of information. This is followed by a discussion of the two major classes of graph theoretical analyses: identifying connectivity patterns and partitioning methods. Measures, metrics, descriptors, and topological indices are then introduced with an emphasis upon enhancing interpretability and incorporation into physical models. Challenging data cases are described that include strategies for studying time dependence. Throughout, we incorporate recent advancements in computer science and applied mathematics that are propelling chemical graph theory into new domains of chemical study. This article is categorized under: Molecular and Statistical Mechanics > Molecular Dynamics and Monte‐Carlo Methods Structure and Mechanism > Computational Materials Science Structure and Mechanism > Molecular Structures

Leite, Leonardo S. G.↗
DOE OSTI2020

STRUCTURE AND PROPERTIES STUDY ON ENERGY MATERIALS: THERMOELECTRIC MATERIAL TETRAHEDRITE AND LITHIUM ION CONDUCTOR LiPON

Development of efficient energy materials is critical in order to ease the energy demand and reduce our dependence of fossil fuel. Thermoelectric materials are promising due to their capability of generating electrical power by recovering waste heat. The performance of thermoelectric materials is quantified by a dimensionless figure of merit zT, which depends on their properties such as electrical conductivity, Seebeck coefficient and thermal conductivity. Tetrahedrites, a copper antimony sulfosalt mineral, typified by Cu12-xMxSb4S13, where M is a transition metal element such as Ni, Zn, Fe or Mn, have great potential for thermoelectric application due to their relatively high zT (close to 1 at 700 K), earth-abundance, environmental friendliness, favorable electrical properties, and most importantly intrinsic low lattice thermal conductivity (less than 1 W m-1 K-1) in wide temperature. In addition to energy recovery, reliable energy storage devices are also emerging to relieve the energy demand and improve the efficiency of consuming energy resources. Lithium-ion batteries are known to be reliable and successful electrochemical energy storage devices and appliable in various aspects, including laptops, smartphones and electrical vehicles. Lithium phosphorous oxynitride (LiPON) are widely used as thin-film solid-state electrolytes in Li-ion battery, which is the only demonstrated solid-state electrolyte that is quite stable in direct contact with Li metal at potentials from 0-5 V. However, the structure of LiPON, the effects of N doping, and the origin of its good electrochemical stability remains inconclusive.In this thesis, reliable modeling techniques accompanied with experimental tools, are applied to study the thermoelectric material tetrahedrite and the ionic conductor LiPON, in order to study their structural and dynamical properties. Accurate and efficient density-functional theory (DFT) and density-functional tight-binding (DFTB) methods, combined with molecular dynamics (MD) simulations are utilized in order to investigate the structures and properties of these energy materials. The incoherent and coherent atomic dynamics study of tetrahedrite Cu10.5NiZn0.5Sb4S13 provides the origin of softening upon cooling by investigate the motion of Cu12e at different temperatures. The dynamic structure factors in the longitudinal and transverse direction will also be discussed. The Cu movement of Cu-rich tetrahedrite Cu14Sb4S13 is revealed by Cu self-diffusivity, nuclear density map and "nudged elastic band" (NEB). Moreover, we investigate the effect of simulation cell size and basis sets on the DFT-based MD simulation results using tetrahedrite Cu10Zn2Sb4S13 thermoelectric as a model material, showing the advantage of larger cell by accessing smaller Q range. In addition, the low-temperature structural properties of Cu12Sb4S13 is measured by neutron diffraction, which indicates that no cubic to tetragonal transition occurs at metal-semiconductor transition (MST) temperature. Thermoelectric properties such as Seebeck coefficient, electrical resistivity and electrical thermal conductivity will also be investigated. DFTB method is implemented to study the structure and transport properties of Li3PO4 and LiPON, while the exploration of N doping effect is included. Lastly, the LiPON/Li interphase will be revealed in order to study the origin of electrochemical stability

25 ENERGY STORAGE↗
DOE OSTI2020

Development of neural network force fields for corrosion studies

To fully understand the chemistry and physics of corrosion, novel methods of simulation must be developed. One approach is designing machine learning (ML) algorithms integrated with density functional theory to develop adaptive force fields to gain insight into corrosion behavior namely at the surface of metal oxides. Current methods of modeling corrosion are slow due to the computational cost of resolving both reaction mechanics and mass transport processes. Machine learning methods can be implemented to obtain structure-activity relationships at both the molecular and bulk scale while still retaining the accuracy of density functional theory (DFT) and significantly decreasing the time needed for simulations of complex chemical processes in the various environments of corrosion. Multiscale models are needed for corrosion studies to fully understand its processes not only at the atomic length scale (chemical bonding, energies, and forces), but also at the nano and meso length scales (solid-state physics and material science processes). Current methods of study include DFT, molecular dynamics, and Monte Carlo. The limitation of DFT is that only a small number of atoms or molecules can be simulated at that level of theory. Density functional theory is used to study the electronic structure of atoms and molecules, and calculate the force component of each atom. However, these calculations are limited to about 1000 atoms. Custom periodic boundary conditions (PBC) can be used to describe the various environments and defects that affect the atomic forces to produce a large data set from which a training set can be derived. Machine learning can be utilized to overcome the barrier of modeling macroscopic and multi-scale processes from ab initio calculations through the development of adaptive force fields. Local environments determine the atomic forces of a given system, therefore adaptive force fields must be created to produce reliable quantum mechanical calculations. This can be achieved by developing a learning algorithm that uses the mapped atomic forces or fingerprint as an input to produce energies and magnetic moments as output. A systematic approach was used to begin to build a data set in order to accurately describe the atomic forces in various environments. In Figure 4 below, a simple PBC cell of Fe{sub 2}O{sub 3} was first optimized. A surface optimization was performed next, followed by a hydroxylated surface optimization. Once this calculation has converged, the adsorption of halide species to the hydroxylated surface will be investigated. TensorFlow is an open source platform for machine learning developed by Google. Using a high level application program interface (API) such as Keras allows for building and training ML models easily in a number of different environments and languages. For this project, a neural network was developed within Anaconda in Python. Future Work: Further development of reference data set; Refining neural network and learning algorithm; Fingerprinting atomic environment to enable mapping of atomic force components; Choosing appropriate training set from reference data; Learning from training set and enabling non-linear mapping of training set fingerprints and the atomic forces; Estimation of uncertainty to identify ranges of outside applicability; Testing and analysis of molecular dynamic simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Breaking the curse of dimensionality: Solving configurational integrals for crystalline solids by tensor networks

Accurately evaluating configurational integrals for dense solids remains a central and difficult challenge in the statistical mechanics of condensed systems. Here, we present a tensor network approach that reformulates the high-dimensional configurational integral for identical-particle crystals into a sequence of computationally efficient summations. We represent the integrand as a high-dimensional tensor and apply tensor-train (TT) decomposition together with a custom TT-cross interpolation. This approach circumvents the need to explicitly construct the full tensor. We introduce tailored rank-1 and rank-2 schemes optimized for sharply peaked Boltzmann probability densities, typical for identical-particle crystals. When applied to the calculation of internal energy and pressure-temperature curves for crystalline Cu and Ar at high (GPa) pressures, as well as the alpha-to-beta phase transition diagram of Sn, our method accurately reproduces molecular dynamics simulation results using tight-binding, machine learning, hierarchical interacting particle–neural network, and modified embedded atom method potentials,all within seconds of computation time.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

STRUCTURE AND PROPERTIES STUDY ON ENERGY MATERIALS: THERMOELECTRIC MATERIAL TETRAHEDRITE AND LITHIUM ION CONDUCTOR

Development of efficient energy materials is critical in order to ease the energy demand and reduce our dependence of fossil fuel. Thermoelectric materials are promising due to their capability of generating electrical power by recovering waste heat. The performance of thermoelectric materials is quantified by a dimensionless figure of merit zT, which depends on their properties such as electrical conductivity, Seebeck coefficient and thermal conductivity. Tetrahedrites, a copper antimony sulfosalt mineral, typified by Cu 12-x MxSb 4 S 13 , where M is a transition metal element such as Ni, Zn, Fe or Mn, have great potential for thermoelectric application due to their relatively high zT (close to 1 at 700 K), earth-abundance, environmental friendliness, favorable electrical properties, and most importantly intrinsic low lattice thermal conductivity (less than 1 W m -1 K -1 ) in wide temperature. In addition to energy recovery, reliable energy storage devices are also emerging to relieve the energy demand and improve the efficiency of consuming energy resources. Lithium-ion batteries are known to be reliable and successful electrochemical energy storage devices and appliable in various aspects, including laptops, smartphones and electrical vehicles. Lithium phosphorous oxynitride (LiPON) are widely used as thin-film solid-state electrolytes in Li-ion battery, which is the only demonstrated solid-state electrolyte that is quite stable in direct contact with Li metal at potentials from 0-5 V. However, the structure of LiPON, the effects of N doping, and the origin of its good electrochemical stability remains inconclusive. In this thesis, reliable modeling techniques accompanied with experimental tools, are applied to study the thermoelectric material tetrahedrite and the ionic conductor LiPON, in order to study their structural and dynamical properties. Accurate and efficient density-functional theory (DFT) and density-functional tight-binding (DFTB) methods, combined with molecular dynamics (MD) simulations are utilized in order to investigate the structures and properties of these energy materials. The incoherent and coherent atomic dynamics study of tetrahedrite Cu 10.5 NiZn 0.5 Sb 4 S 13 provides the origin of softening upon cooling by investigate the motion of Cu12e at different temperatures. The dynamic structure factors in the longitudinal and transverse direction will also be discussed. The Cu movement of Cu-rich tetrahedrite Cu 14 Sb 4 S 13 is revealed by Cu self-diffusivity, nuclear density map and “nudged elastic band” (NEB). Moreover, we investigate the effect of simulation cell size and basis sets on the DFT-based MD simulation results using tetrahedrite Cu1 0 Zn 2 Sb 4 S 13 thermoelectric as a model material, showing the advantage of larger cell by accessing smaller Q range. In addition, the low-temperature structural properties of Cu 12 Sb 4 S 13 is measured by neutron diffraction, which indicates that no cubic to tetragonal transition occurs at metal-semiconductor transition (MST) temperature. Thermoelectric properties such as Seebeck coefficient, electrical resistivity and electrical thermal conductivity will also be investigated. DFTB method is implemented to study the structure and transport properties of Li 3 PO 4 and LiPON, while the exploration of N doping effect is included. Lastly, the LiPON/Li interphase will be revealed in order to study the origin of electrochemical stability.

Li, Junchao↗
DOE OSTI2024

Generation of Enrichment-Dependent Thermal Neutron Scattering Data

This work details the generation of enrichment-dependent thermal neutron scattering cross sections for several crucial uranium fuel compounds. The evaluations of the thermal scattering law (TSL) and associated cross sections for uranium dioxide (UO 2 ), uranium carbide (UC), and uranium nitride (UN) were performed using standard ab initio lattice dynamics (AILD) methods. The data for uranium metal was produced using a novel hybrid approach of molecular dynamics combined with lattice dynamics methods. 235 U enrichments of 5%, 10% (LEU+), 19.75% (HALEU), 93% (HEU), and 100% were considered, in addition to natural uranium. The enrichment-dependent masses and free atom cross sections were used in the generation of elastic and inelastic thermal neutron scattering cross sections, while the calculation of the phonon density of states (DOS) and resulting TSL considered only the natural isotopic composition of uranium. The use of an identical DOS for all enrichments is expected to have minimal impact on the final data, as the small change in uranium mass should not significantly affect lattice vibrations. The cross sections are shown to exhibit significant dependence on 235 U enrichment. The submission of this data to the National Nuclear Data Center (NNDC) for release in the ENDF/B-VIII.1 database should support the design of advanced reactor concepts.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Weighted active space protocol for multireference machine-learned potentials

Multireference methods such as multiconfiguration pair-density functional theory accurately capture electronic correlation in systems with strong multiconfigurational character, but their cost precludes direct use in molecular dynamics. Combining these methods with machine-learned interatomic potentials (MLPs) can extend their reach. However, the sensitivity of multireference calculations to the choice of the active space complicates the consistent evaluation of energies and gradients across structurally diverse nuclear configurations. To overcome this limitation, we introduce the weighted active space protocol (WASP), a systematic approach to assign a consistent active space for a given system across uncorrelated configurations. By integrating WASP with MLPs and enhanced sampling techniques, we propose a data-efficient active learning cycle that enables the training of an MLP on multireference data. We demonstrated the approach on the TiC + -catalyzed C–H activation of methane, a reaction that poses challenges for Kohn–Sham density functional theory due to its significant multireference character. This framework enables accurate and efficient modeling of catalytic dynamics, establishing a paradigm for simulating complex reactive processes beyond the limits of conventional electronic-structure methods.

enhanced sampling↗