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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Active learning of reactive Bayesian force fields applied to heterogeneous catalysis dynamics of H/Pt

Abstract Atomistic modeling of chemically reactive systems has so far relied on either expensive ab initio methods or bond-order force fields requiring arduous parametrization. Here, we describe a Bayesian active learning framework for autonomous “on-the-fly” training of fast and accurate reactive many-body force fields during molecular dynamics simulations. At each time-step, predictive uncertainties of a sparse Gaussian process are evaluated to automatically determine whether additional ab initio training data are needed. We introduce a general method for mapping trained kernel models onto equivalent polynomial models whose prediction cost is much lower and independent of the training set size. As a demonstration, we perform direct two-phase simulations of heterogeneous H 2 turnover on the Pt(111) catalyst surface at chemical accuracy. The model trains itself in three days and performs at twice the speed of a ReaxFF model, while maintaining much higher fidelity to DFT and excellent agreement with experiment.

42 ENGINEERING↗
DOE OSTI2025

Mapping structures and dynamics with frequency-correlated diffusion exchange

Understanding molecular motion in diffusion-driven complex environments is critical for designing sustainable materials and improving chemical processes. Here, we introduce a multidimensional nuclear magnetic resonance (NMR) method that captures how molecular populations exchange across different dynamic regimes. By extending the modulated gradient spin-echo technique to include frequency-frequency correlations, our approach reveals diffusion pathways that are otherwise obscured in heterogeneous systems. Implemented on a unilateral NMR magnet, the method eliminates gradient pulsing constraints and accesses dynamics in the kilohertz regime. We apply this technique to swelling and acid-catalyzed deconstruction of cross-linked and linear polymers to observe how structural heterogeneity evolves over time. By linking molecular motion to topology and chemical state, we extract physical metrics such as fractal surface dimensionality and reaction wavefront velocity, properties inaccessible with standard diffusion measurements. This work expands the capabilities of NMR for probing soft matter, with implications for polymer recycling and materials design.

Fricke, Sophia N. [University of California, Berke↗
DOE OSTI2022

Efficient and universal characterization of atomic structures through a topological graph order parameter

Abstract A graph-based order parameter, based on the topology of the graph itself, is introduced for the characterization of atomistic structures. The order parameter is universal to any material/chemical system and is transferable to all structural geometries. Four sets of data are used to validate both the generalizability and accuracy of the algorithm: (1) liquid lithium configurations spanning up to 300 GPa, (2) condensed phases of carbon along with nanotubes and buckyballs at ambient and high temperature, (3) a diverse set of aluminum configurations including surfaces, compressed and expanded lattices, point defects, grain boundaries, liquids, nanoparticles, all at nonzero temperatures, and (4) eleven niobium oxide crystal phases generated with ab initio molecular dynamics. We compare our proposed method to existing, state-of-the-art methods for the cases of aluminum and niobium oxide. Our order parameter uniquely classifies every configuration and outperforms all studied existing methods, opening the door for its use in a multitude of complex application spaces that can require fine structure-level characterization of atomistic graphs.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A new self-adaptive reconstruction method to identify defects through Wigner–Seitz approach

A new self-adaptive reconstruction method based on local atomic structure at any given molecular dynamics (MD) step has been developed in this article. The method can be used in Wigner–Seitz defect analysis approach to correctly and efficiently explore the information of both point defects and complex defect clusters (e.g. dislocation loops and voids) formed after a displacement cascade where the cascade interacts with grain boundaries and/or dislocations. The algorithm and validation are provided in detail. Results for identification of radiation defects during and after cascades interacting with a dislocation network show that the new method can well recognize all simple and complex defects and defect clusters. Thus, this new method provides a totally new way to explore the density and size of radiation defects at atomic scale after complex MD evolution processes, providing correct information to understand and predict radiation damage in materials through atomic simulations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Differences in Sb 2 Te 3 growth by pulsed laser and sputter deposition

High quality van der Waals chalcogenides are important for phase change data storage, thermoelectrics, and spintronics. Using a combination of statistical design of experiments and density functional theory, we clarify how out-of-equilibrium van der Waals epitaxial deposition methods can improve the crystal quality of Sb 2 Te 3 films. We compare films grown by radio frequency sputtering and pulsed laser deposition (PLD). The growth factors that influence the crystal quality for each method are different. For PLD grown films a thin amorphous Sb 2 Te 3 seed layer most significantly influences the crystal quality. In contrast, the crystalline quality of films grown by sputtering is rather sensitive to the deposition temperature and less affected by the presence of a seed layer. This difference is somewhat surprising as both methods are out-of-thermal-equilibrium plasma-based methods. Non-adiabatic quantum molecular dynamics simulations show that this difference originates from the density of excited atoms in the plasma. The PLD plasma is more intense and with higher energy than that used in sputtering, and this increases the electronic temperature of the deposited atoms, which concomitantly increases the adatom diffusion lengths in PLD. In contrast, the adatom diffusivity is dominated by the thermal temperature for sputter grown films. These results explain the wide range of Sb 2 Te 3 and superlattice crystal qualities observed in the literature. Finally, these results indicate that, contrary to popular belief, plasma-based deposition methods are suitable for growing high quality crystalline chalcogenides.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Nuclear quantum effects on the quasiparticle properties of the chloride anion aqueous solution within the GW approximation

Photoelectron spectroscopy experiments in ionic solutions reveal important electronic structure information, in which the interaction between hydrated ions and water solvent can be inferred. Based on many-body perturbation theory with GW approximation, we theoretically compute the quasiparticle electronic structure of chloride anion solution, which is modeled by path-integral abinitio molecular dynamics simulation by taking account the nuclear quantum effects (NQEs). The electronic levels of hydrated anion as well as water are determined and compared with the recent experimental photoelectron spectra. It is found that NQEs improve the agreement between theoretical prediction and experiment because NQEs effectively weaken the hybridization of the between the Cl – anion and water. Furthermore, our study indicates that NQEs plays a small but non-negligible role in predicting the electronic structure of the aqueous solvation of ions of the Hofmeister series.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Molecular Driving Force for Facet Selectivity of Sequence-Defined Amphiphilic Peptoids at Au–Water Interfaces

Shape-controlled colloidal nanocrystal syntheses often require facet-selective solution-phase chemical additives to regulate surface free energy, atom addition/migration fluxes or particle attachment rates. Because of their highly tunable properties and robustness to a wide range of experimental conditions, peptoids represent a very promising class of next generation functional additives for control over nanocrystal growth. However, understanding the origin of facet selectivity at the molecular level is critical to generalizing their design. In this work, we employ molecular dynamics simulations and biased sampling methods and report stronger selectivity to Au(111) than to Au(100) for Nce3Ncp6, a peptoid that has been shown to assist the formation of five-fold twinned Au nanostars. We find that facet-selectivity is achieved through synergistic effects of both peptoid-surface and solvent-surface interactions. Moreover, the amphiphilic nature of Nce3Ncp6 together with the order of peptoid-peptoid and peptoid-surface binding energies, i.e., peptoid-Au(100) < peptoid-peptoid < peptoid-Au(111), further amplifies its distinct collective behavior on different Au surfaces. Our studies provide a fundamental understanding of the molecular origin of facet-selective adsorption and highlight the possibility of future designs and uses of sequence-defined peptoids for predictive syntheses of nanocrystals with designed shapes and properties.

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

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.

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

The behavior of methane–water mixtures under elevated pressures from simulations using many-body potentials

Non-polarizable empirical potentials have been proven to be incapable of capturing the mixing of methane–water mixtures at elevated pressures. Although density functional theory-based ab initio simulations may circumvent this discrepancy, they are limited in terms of the relevant time and length scales associated with mixing phenomena. Here, we show that the many-body MB-nrg potential, designed to reproduce methane–water interactions with coupled cluster accuracy, successfully captures this phenomenon up to 3 GPa and 500 K with varying methane concentrations. Two-phase simulations and long time scales that are required to fully capture the mixing, affordable due to the speed and accuracy of the MBX software, are assessed. Constructing the methane–water equation of state across the phase diagram shows that the stable mixtures are denser than the sum of their parts at a given pressure and temperature. As a result, we find that many-body polarization plays a central role, enhancing the induced dipole moments of methane by 0.20 D during mixing under pressure. Overall, the mixed system adopts a denser state, which involves a significant enthalpic driving force as elucidated by a systematic many-body energy decomposition analysis.

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

Exploring cryo-electron microscopy with molecular dynamics

Single particle analysis cryo-electron microscopy (EM) and molecular dynamics (MD) have been complimentary methods since cryo-EM was first applied to the field of structural biology. The relationship started by biasing structural models to fit low-resolution cryo-EM maps of large macromolecular complexes not amenable to crystallization. The connection between cryo-EM and MD evolved as cryo-EM maps improved in resolution, allowing advanced sampling algorithms to simultaneously refine backbone and sidechains. Moving beyond a single static snapshot, modern inferencing approaches integrate cryo-EM and MD to generate structural ensembles from cryo-EM map data or directly from the particle images themselves. We summarize the recent history of MD innovations in the area of cryo-EM modeling. The merits for the myriad of MD based cryo-EM modeling methods are discussed, as well as, the discoveries that were made possible by the integration of molecular modeling with cryo-EM. Lastly, current challenges and potential opportunities are reviewed.

Biochemistry & Molecular Biology↗
DOE OSTI2020

Using density functional theory to construct multiphase equations of state: tin as an example

We present a general methodology for using density functional theory (DFT) calculations as a basis for constructing multiphase equations of state (EOS). Focusing on tin as an example, we discuss the full process of generating a tabular EOS, starting from first-principles calculations and arriving at a full EOS in agreement with existing experimental data. We begin by describing DFT calculations for the five solid phases and liquid phase of tin, including cold curve, phonon, and DFT-based molecular dynamics calculations. We then discuss methods for incorporating DFT results into materials models used in OpenSesame, the program used to create a full tabular multiphase EOS. Next, we outline a general strategy for adjusting model parameters in OpenSesame to fit to a variety of experimental data, including isobaric data, isothermal data, shock data, solid-solid phase boundary measurements, and measurements of the melt curve. We end with a discussion of the advantages of using DFT data as a foundation for EOS generation and discuss potential strategies for automated EOS generation based on the fitting process highlighted in this work.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Stochastic and mixed density functional theory within the projector augmented wave formalism for simulation of warm dense matter

Stochastic density functional theory (DFT) and mixed stochastic-deterministic DFT are burgeoning approaches for the calculation of the equation of state and transport properties in materials under extreme conditions. In the intermediate warm dense matter regime, a state between correlated condensed matter and kinetic plasma, electrons can range from being highly localized around nuclei to delocalized over the whole simulation cell. The plane-wave basis pseudopotential approach is thus the typical tool of choice for modeling such systems at the DFT level. Unfortunately, stochastic DFT methods scale as the square of the maximum plane-wave energy in this basis. To reduce the effect of this scaling and improve the overall description of the electrons within the pseudopotential approximation, we present stochastic and mixed DFT approaches developed and implemented within the projector augmented wave formalism. In conclusion, we compare results between the different DFT approaches for both single-point and molecular dynamics trajectories and present calculations of self-diffusion coefficients of solid density carbon from 1 to 50 eV.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Neural Phase Simulation

The Neural Phase Simulation (NPS) is a package of codes for simulating microstructure evolution and accelerated molecular dynamics with deep neural-networks based surrogate models. NPS is designed to offer quantitatively accurate and computationally efficient simulation capabilities by leveraging modern machine-learning techniques. The primary intended use cases of NPS are training neural network surrogate models, though performing simulations on a single node is also supported. The NPS surrogate models can be trained from ground truth simulation methods, which are supposed to be accurate but expensive, such as molecular dynamics, phase field methods, kinetic Monte Carlo and discrete dislocation dynamics.

Zhou, Fei↗
DOE OSTI2025

Understanding the structural mechanics of ligated DNA crystals via molecular dynamics simulation

DNA self-assembly is a highly programmable method to construct arbitrary architectures based on sequence complementarity. Among various constructs, DNA crystals are macroscopic crystalline materials formed by assembling motifs via sticky end association. Due to their high structural integrity and size ranging from tens to hundreds of micrometers, DNA crystals offer unique opportunities to study the structural properties and deformation behaviors of DNA assemblies. For example, enzymatic ligation of sticky ends can selectively seal nicks resulting in more robust structures with enhanced mechanical properties. However, the research efforts have been mostly on experiments involving different motif designs, structural optimization, or new synthesis methods, while their mechanics are not yet fully understood. The complex properties of DNA crystals are difficult to study via experiments alone, and numerical simulation can complement and aid the experiments. The coarse-grained molecular dynamics (MD) simulation is a powerful tool that can probe the mechanics of DNA assemblies. Here, we investigate DNA crystals made of four different motif lengths with various ligation patterns (full ligation, major directions, connectors, and in-plane) using oxDNA, an open-source, coarse-grained MD platform. We found that several distinct deformation stages emerge in response to mechanical loading and that the number and the location of ligated nucleotides can significantly modulate structural behaviors. These findings should be useful for predicting crystal properties and thus improving the design.

DNA crystal↗
DOE OSTI2020

Steering Molecular Dynamics Simulations of Membrane-Associated Proteins with Neutron Reflection Results

Herein, we present a novel method to incorporate structural results from surface-sensitive scattering, such as X-ray or neutron reflectometry, into molecular dynamics simulations. While reflectometry techniques generally provide a means to determine the molecular-scale structures of organized interfacial films, they were recently shown to offer the capability to characterize the structures of protein-membrane complexes supported by a solid substrate. One-dimensional information inherent in the experimental results is used in the form of component volume occupancy (CVO) profiles, which describe the distribution of molecular components within an interfacial architecture, to construct real-space constraints in the form of a biasing potential for the simulation that vanishes when the simulated and experimental profiles agree. This approach improves the correspondence between simulation and experiment, as shown in the re-evaluation of an neutron-reflection-derived structure which was approximated by an independent molecular dynamics simulation in earlier work, and it also leads to faster equilibration of ensemble structures. We further show that time averaging the CVO profile that develops in the simulation while biasing with this approach permits fluctuations about the average that are necessary for conformational exploration of the system. This method is particularly valuable for studies of proteins at interfaces that contain disordered regions since the conformation of such regions is difficult to judge from the analysis of one-dimensional experimental profiles and may take prohibitively long to equilibrate in simulations.

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

Combining Experimental and Theoretical Techniques to Gain an Atomic Level Understanding of the Defect Binding Mechanism in Hard Carbon Anodes for Sodium Ion Batteries

Sodium ion batteries (NIBs) are an attractive alternative to lithium–ion batteries in applications that require large–scale energy storage due to sodium's high natural abundance and low cost. Hard carbon (HC) is the most promising anode material for NIBs; however, there is a knowledge gap in the understanding of the sodium binding mechanism that prevents a rational design of HC. This study tunes sucrose–derived HC via synthesis temperature then evaluates the structural, physical, and electrochemical properties. Neutron total scattering is used to generate structural models by fitting pair distribution functions (PDF) with a combination of molecular dynamics and reverse Monte Carlo methods. From this model, the number and type of structural features are identified, quantified, and correlated to the galvanostatic charge/discharge. A method of PDF “fingerprinting” binding sites using Na probe atoms is developed and analyzing these PDFs reveals an atomistic view of ion binding sites responsible for “defect” storage mechanisms. Combining these techniques results in an atomic–level study that provides a big picture of the Na–binding mechanism in NIBs, which allows for more precise tuning of the structure–property relationships in the future. Finally, the methodologies developed will also enable new strategies for the analysis of amorphous functional materials.

25 ENERGY STORAGE↗
DOE OSTI2023

An interoperable implementation of collective‐variable based enhanced sampling methods in extended phase space within the OpenMM package

Collective variable (CV)-based enhanced sampling techniques are widely used today for accelerating barrier-crossing events in molecular simulations. A class of these methods, which includes temperature accelerated molecular dynamics (TAMD)/driven-adiabatic free energy dynamics (d-AFED), unified free energy dynamics (UFED), and temperature accelerated sliced sampling (TASS), uses an extended variable formalism to achieve quick exploration of conformational space. These techniques are powerful, as they enhance the sampling of a large number of CVs simultaneously compared to other techniques. Extended variables are kept at a much higher temperature than the physical temperature by ensuring adiabatic separation between the extended and physical subsystems and employing rigorous thermostatting. Here, in this work, we present a computational platform to perform extended phase space enhanced sampling simulations using the open-source molecular dynamics engine OpenMM. The implementation allows users to have interoperability of sampling techniques, as well as employ state-of-the-art thermostats and multiple time-stepping. This work also presents protocols for determining the critical parameters and procedures for reconstructing high-dimensional free energy surfaces. As a demonstration, we present simulation results on the high dimensional conformational landscapes of the alanine tripeptide in vacuo, tetra-N-methylglycine (tetra-sarcosine) peptoid in implicit solvent, and the Trp-cage mini protein in explicit water.

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

An adaptive granular representative volume element model with an evolutionary periodic boundary for hierarchical multiscale analysis

The hierarchical multiscale analysis normally utilizes a microscopic representative volume element (RVE) model to capture path/history-dependent macroscopic responses instead of using phenomenological constitutive models. However, for problems involving large deformation, the current RVE model used in geomechanics may lose representative properties due to the progressive distortion of the RVE box, unless a proper reinitialization is applied. This work develops an adaptive RVE model in conjunction with an evolutionary periodic boundary (EPB) algorithm for hierarchical multiscale analysis of granular materials undergoing large deformation based on a recent RVE model proposed for coupling molecular dynamics and the material point method. The proposed adaptive RVE model avoids the reinitialization of the RVE box that even undergoes extremely large shear deformation; meanwhile, it accounts for the deformation history of the RVE model and treats the interaction between boundary particles and other image particles in a more efficient way. Numerical examples with extremely large deformation are used to illustrate the adaptive granular RVE model enhanced by the proposed EPB algorithm. Furthermore, some key features of this new methodology are further discussed for clarification.

42 ENGINEERING↗