Search NASA⌕ Search

SEARCH · Search NASA

Results for “Machine learning molecular dynamics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Integrating machine learning interatomic potentials with hybrid reverse Monte Carlo structure refinements in RMCProfile

Structure refinement with reverse Monte Carlo (RMC) is a powerful tool for interpreting experimental diffraction data. To ensure that the under-constrained RMC algorithm yields reasonable results, the hybrid RMC approach applies interatomic potentials to obtain solutions that are both physically sensible and in agreement with experiment. To expand the range of materials that can be studied with hybrid RMC, we have implemented a new interatomic potential constraint in RMCProfile that grants flexibility to apply potentials supported by the Large-scale Atomic/Molecular Massively Parallel Simulator ( LAMMPS ) molecular dynamics code. This includes machine learning interatomic potentials, which provide a pathway to applying hybrid RMC to materials without currently available interatomic potentials. To this end, we present a methodology to use RMC to train machine learning interatomic potentials for hybrid RMC applications.

Cuillier, Paul↗
DOE OSTI2025

Data from "Deep Potential Molecular Dynamics Simulations of Low-Temperature Plasma-Surface Interactions"

Data and input files related to the paper "Deep Potential Molecular Dynamics Simulations of Low-Temperature Plasma-Surface Interactions" (https://doi.org/10.1116/6.0004027). This includes the final DP model used in all simulations, training data set, example input files to run DeepMD (with LAMMPS), and data tables summarizing the results obtained from the simulations.

machine learning models↗
DOE OSTI2021

Atomic understanding of structural deformations upon ablation of graphene

We investigate the atomic rearrangement in graphene under femtosecond pulse illumination with reactive molecular dynamics simulations and compare with ultra-fast laser ablation experiments. To model the impact of the laser pulse irradiation, heat is locally applied to a selected area of the graphene layer and the resulting structural deformation is simulated as a function of time, providing a detailed understanding of the bond breaking process under laser illumination and subsequent re-equilibration afterAQ3 the pulse is turned off. Analysis of the atomic dynamics indicates that the types of defects formed depend on the pulse energy and exposure duration. By varying the exposed area, we determine that the shape of the ablated area is not only a function of the pulse energy, but also of the beam spot size and pulse repetition. Furthermore, we apply a machine learning approach to extrapolate our simulated data to experimental length scales and reproduce the trends in ablated area as a function of temperature. Furthermore, our study provides a first step towards understanding the design parameters for graphene nano-patterning.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

Simulating Crystallization in a Colloidal System Using State Predictive Information Bottleneck Based Enhanced Sampling

Here, we investigate crystal nucleation in supersaturated colloid suspensions using enhanced molecular dynamics simulations augmented with machine learning techniques. The simulations reveal that crystallization in the model colloidal system studied here, with particles interacting through a repulsive screened Coulomb Yukawa potential, proceeds from vapor to dense liquid droplet to crystalline phases across multiple high barriers. Employing a one-dimensional reaction coordinate derived from the State Predictive Information Bottleneck framework, our simulations capture back-and-forth phase transitions across multiple barriers effectively in biased metadynamics simulations. We obtain relative free energy differences between different phases and also quantify the roles of different molecular level features in driving the phase changes.

Chemistry↗
DOE OSTI2019

Interatomic Potential Model Development: Finite-Temperature Dynamics Machine Learning

Developing an accurate interatomic potential model is a prerequisite for achieving reliable results from classical molecular dynamics (CMD) simulations; however, most of the potentials are biased as specific simulation purposes or conditions are considered in the parameterization. Here, for developing an unbiased potential, a finite-temperature dynamics machine learning (FTD-ML) approach is proposed, and its processes and feasibility are demonstrated using the Buckingham potential model and aluminum (Al) as an example. Compared with conventional machine learning approaches, FTD-ML exhibits three distinguished features: 1) FTD-ML intrinsically incorporates more extensive configurational and conditional space for enhancing the transferability of developed potentials; 2) FTD-ML employs various properties calculated directly from CMD, for ML model training and prediction validation against experimental data instead of first-principles data; 3) FTD-ML is much more computationally cost effective than first-principles simulations, especially when the system size increases over 10 3 atoms as employed in this research for ensuring reliable training data. The Al Buckingham potential developed by the FTD-ML approach exhibits good performance for general simulation purposes. Thus, the FTD-ML approach is expected to contribute to a fast development of interatomic potential model suitable for various simulation purposes and conditions, without limitation of model type, while maintaining experimental-level accuracy.

42 ENGINEERING↗
DOE OSTI2025

Efficient machine learning interatomic potentials robust for liquid and multiple solid polymorphs of NaF and KF

Achieving atomic-level understanding of crystallization of molten salts is of importance to a wide range of technological applications. Recent work [Fan et al., Proc. Natl. Acad. Sci. USA 122, e2425702122 (2025)] revealed that crystal nucleation in molten LiF salt is a multistage process according to the molecular-dynamics (MD) simulations based on an atomic cluster expansion (ACE) machine-learning interatomic potential (MLIP). In order to understand the influence of increasing cation size on nucleation pathways and nucleation rates of molten fluoride salts, here we develop two new ACE MLIPs for NaF and KF. The two ACE MLIPs feature DFT-SCAN-level accuracy for liquid and multiple solid polymorphs over a wide temperature (0–2000 K) and pressure (0–100 GPa) range, and also reproduce well a number of experimental data for solid and liquid equilibrium properties. The efficiency of the two ACE MLIPs enable million-atom-scale or microsecond-scale MD simulations. The two general-purpose ACE MLIPs are expected to be useful for atomistic simulations for different purposes, in addition to studying crystallization of molten NaF and KF salts.

Crystal melting↗
DOE OSTI2026

From bulk to surface: Structure and dynamics of amorphous alumina from deep potential molecular dynamics

Understanding the atomic-scale structure and dynamics of amorphous oxide surfaces is essential for interpreting their chemical reactivity, mechanical stability, and interfacial behavior, yet direct experimental characterization remains challenging. We employ Deep Potential (DP) molecular dynamics to generate large-scale, ab initio -quality models of amorphous Al 2 O 3 bulk glasses and melt-quenched free surfaces, enabling a quantitative analysis of both structure and relaxation dynamics with statistical confidence inaccessible to direct ab initio simulation. The trained DP model reproduces experimental liquid and glass structure, captures the cooling-rate dependence of the bulk glass transition, and corrects systematic biases in the polyhedral populations predicted by widely used classical force fields. At the free surface, mass density recovers to bulk values over ~10 Å, while local coordination requires a slightly wider subsurface region to fully converge. The outermost layer is oxygen-enriched, exhibits altered polyhedral connectivity with contracted Al–O bonds, and hosts a broad population of under-coordinated motifs (notably AlO 3 and OAl 2 ) whose abundances are governed by glass stability. These under-coordinated surface motifs exhibit distinct vibrational signatures and occur as locally paired Lewis acid and Brønsted base sites consistent with bond-valence compensation, yet remain spatially dispersed rather than aggregating into extended clusters. Despite this pronounced structural heterogeneity, surface relaxation and the glass-transition temperature remain comparable to their bulk counterparts, suggesting that the disordered surface is kinetically stable once formed. Together, these results establish a molecular-level picture of amorphous alumina surfaces and demonstrate the capability of machine-learned potentials to resolve structure–property relationships in disordered oxide interfaces.

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

Biochemical Process Modeling and Simulation (BPMS)

The Biochemical Process Modeling and Simulation project aims to reduce the cost and time of research by applying theory, modeling, and simulation to the most relevant bottlenecks in the biochemical process. We use molecular modeling, quantum mechanics, metabolic modeling, fluid dynamics, and reaction-diffusion methods in close collaboration with pretreatment, hydrolysis, upgrading, and TEA. The project's outcomes are increased yields and efficiency of the biochemical process, added value to products, and reduced price of fuels by specifically targeting catalytic efficiency, reactor design, enzyme efficiency, and microbial design. We work closely with experimental projects to identify problems and iterate with experiments to find and refine solutions. By working with experimentalists, we decide on problems that can be solved with simulation that could otherwise not be solved or would take too long with experiment alone to reach BETO's targets. Over the years, we have produced solutions that have resulted in determining the most likely fatty-acid derivative for passive transport out of bacteria that upgrade biomass, and we have also designed enzyme mutations for enhanced lignin upgrading. Metabolic models have been developed to tune the activity of 2,3 butanediol production for the 2030 target. A computational method to deliver understanding of how complex omics data can be interpreted in the metabolic pathways of organisms used in the Agile Biofoundry. We have found methods to overcome specific barriers and continue to develop those methods. Our reactor studies have guided the design of both the microbes and reactors for aerobic and micro-aerobic production at all scales and have been instrumental in improving the accuracy of techno-economic analysis models. This project is essential in the process of selecting the final processes for 2030 SAF production targets. More specifically, recently, we have: 1) Predicted the strength of the basic structural interactions in commodity plastics to provide guidance for plastics upcycling strategies. 2) Developed computational tool to improve the characterization of lignin-derived compounds 3) Developed new methodologies to enable Machine Learning-based Directed Evolution for protein engineering. 4) Developed Machine Learning methods to predict protein promiscuity and mutations to further improve microbial and enzymatic driven processes and demonstrated the utility of ML approaches to engineering proteins from sparse experimental datasets. 5) Developed new methods to enable high-fidelity simulation of aerobic fermentation at industrial scale and resolving mismatch of time scales through subcycling/operator splitting 7) Identified the difficulty in preventing local high-oxygen conditions in industrial bubble columns, which leads to less-desirable acetoin production, suggesting future research directions in alternative reactor configurations (e.g loop reactors, shallow-channel reactors).

BIOMASS FUELS↗
DOE OSTI2025

Molecular simulation and artificial intelligence for the circular economy of bioenergy and bioproducts

The concept of the circular bioeconomy is a carbon neutral, sustainable system with zero waste. One vision for such an economy is based upon lignocellulosic biomass. This lignocellulosic circular bioeconomy requires CO 2 absorption from biomass growth and the efficient deconstruction of recalcitrant biomass into solubilized and fractionated biopolymers, which are then used as precursors for the sustainable production of high-quality liquid fuels, chemical bioproducts, and bio-based materials. Here, in this study, we summarize the roles that molecular dynamics (MD) simulations and machine learning (ML) are playing in overcoming several fundamental challenges hindering the adoption of a circular bioeconomy. Specifically, we discuss the role of MD and ML/AI in overcoming lignocellulose recalcitrance by designing biomass pretreatment methods to efficiently produce solubilized cellulose/lignin/hemicellulose and of that in improving energy-intensive manufacturing of biomass-based materials and their structural and mechanical properties. Quantum mechanical methods and MD simulations, in addition to offering a mechanistic understanding of biomass deconstruction and biomaterials design, can provide meaningful structural, energetics, and physiochemical properties as inputs to train AI/ML models. The ML models can guide the experimental prioritization of materials/solvents and process parameters that significantly accelerate the development of biofuel and biomaterial components of the circular bioeconomy.

Smith, Jeremy C. [Oak Ridge National Laboratory (O↗
DOE OSTI2021

Azo(xy) vs Aniline Selectivity in Catalytic Nitroarene Reduction by Intermetallics: Experiments and Simulations

Intermetallic nanoparticles are promising catalysts in hydrogenation and fuel cell technologies. Much is known about the ability of intermetallic nanoparticles to selectively reduce nitro vs alkene, alcohol, or halide functional groups; less is known about their selectivity toward aniline vs azo or azoxy condensation products that result from the reduction of a nitro group alone. Because azo(xy)arenes bear promise as dyes, chemical stabilizers, and building blocks to functional materials but can be difficult to isolate, developing high surface area nanoparticle catalysts that display azo(xy) selectivity is desirable. To address this question, we studied a family of nanocrystalline group 10 metal (Pd, Pt)- and group 14 metal (Ge, Sn, Pb)-containing intermetallics-Pd 2 Ge, Pd 2 Sn, Pd 3 Sn 2 , Pd 3 Pb, and PtSn-in the catalytic reduction of nitroarenes. In contrast to monometallic Au, Pt, and Pd nanoparticles and "random" PdxSn1 – x nanoalloys, which are selective for aniline, nanoparticles of atomically precise intermetallic Pd 2 Ge, Pd 2 Sn, Pd 3 Sn 2 , and PtSn prefer an indirect condensation pathway and have a high selectivity for the azo(xy) products. The only exception is Pd 3 Pb, the most active among the intermetallic nanoparticles studied here, which is instead selective for aniline. Employing a novel application of molecular dynamics-based on machine learned potentials within a DeePMD framework-to heterogeneous catalysis, we are able to identify key reaction species on the different types of catalysts employed, furthering our understanding of the unique selectivity of these materials. By demonstrating how intermetallic nanoparticles can be as active yet more selective than other more traditional catalysts, this work provides new physical insights and opens new opportunities in the use of these materials in other important chemical transformations and applications.

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

Theoretical Investigations of Anharmonic Effects and Phonon Transport in the Cubic Phase of Crystalline Perovskite CsPbCl 3

The role of anharmonic effects on lattice dynamics and thermal transport has been investigated in the cubic phase of the CsPbCl 3 perovskite. Limitations of the harmonic approximation, which lead to phonon instabilities in the Brillouin zone, were addressed based on self-consistent phonon theory in combination with first-principles calculations and by incorporation of frequency renormalization effects from bubble and loop diagrams. This theoretical approach demonstrates significant improvements compared to results obtained using a self-consistent phonon dynamic matrix. Both the cubic-to-tetragonal phase transition temperature as well as the predicted lattice thermal conductivity values show good agreement with selected sets of available experimental data. In conclusion, the accuracy of the predicted thermal conductivity data is also tested against systems significantly larger than those allowed by quantum calculations, by using molecular dynamic simulations with machine-learning force fields.

Lattices↗
DOE OSTI2025

Electric Field’s Dueling Effects through Dehydration and Ion Separation in Driving NaCl Nucleation at Charged Nanoconfined Interfaces

Investigating nucleation in charged nanoconfined environments under electric fields is crucial for many scientific and engineering applications. Here we study the nucleation of NaCl from aqueous solution near charged surfaces using machine-learning-augmented enhanced sampling molecular dynamics simulations. Our simulations successfully drive phase transitions between the liquid and solid phases of NaCl. The solid phase is stabilized under electric fields, particularly at an intermediate surface charge density. We examine which physical characteristics drive the nucleation of NaCl from aqueous solutions and find that the removal of solvent water from Cl– at the solid precursor surface plays a more critical role than the accumulation of ions. Our simulations reveal the competing effects of electric fields on nucleation processes: they facilitate the removal of water, promoting nucleation, but also promote the separation of ion pairs, thereby hindering nucleation. Here, this work provides a framework for studying nucleation processes in nanoconfined environments under electric fields and provides physical insights for the design of electrochemistry materials.

Electric fields↗
DOE OSTI2026

Tandem bulk oxygen diffusion and surface reactions in reducible metal oxides control redox cycle dynamics

The interplay between bulk oxygen diffusion and surface reactions in reducible metal oxides is key in heterogeneous catalysts, but direct measurements of oxygen mobility, transient kinetics, and in situ spectroscopies have been lacking. Here, we reveal complex dynamic behavior of ceria-zirconia by H 2 using transient kinetics via mass spectrometry and in situ Raman and near-ambient pressure x-ray photoelectron spectroscopies. Molecular dynamics simulations with a machine learning potential delineate competitive oxygen diffusion mechanisms, with an optimal mobility at intermediate reductions. We expose a compensation between vacancy availability and lattice distortion at intermediate to high reductions and Frenkel defects at low reductions, underscoring a potential deficiency of 16 O/ 18 O exchange experiments in deducing oxygen mobility. Vacancies in proximity require electron localization on Ce atoms further away. The continuous replenishment of surface oxygen results in a varying reduction rate, with H 2 dissociation being the rate-limiting step. Multiscale transient simulations, consistent with experiments, indicate catalysts of potentially spatially varying oxidation states. The approach is broadly applicable to reducible oxide materials.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Solid-state NMR and theoretical studies illuminate lanthanum borohydride C–H borylation catalysts confined within a zeolite

Zeolite-supported single-site lanthanum borohydride catalyzes C–H borylation of hydrocarbons, while the related silica-supported complex is inactive under comparable conditions. The identical composition of support-La(BH 4 ) 2 (THF) 2 sites in the two materials implies that the support bestows underlying structural features onto La that are required for bond activation catalysis, yet the nature of these effects, which could include confinement effects and/or electronic modulation of the site itself, remains to be identified. We used solid-state nuclear magnetic resonance (SSNMR) spectroscopy and molecular dynamics simulations with machine-learning potentials (ML-MD) to analyze the electronic and steric effects imparted by the faujasite support on the precatalyst structure to correlate with catalytic activity. ML-MD simulations show that THF dissociates from La under the influence of confinement, leading to coordinatively unsaturated sites in the zeolite pores. Then, the La complex grafts on Brønsted acid sites (La BAS ) or isolated silanols (La SiO ) or remains physisorbed in the zeolite pores. Catalytic studies comparing compounds supported on faujasite zeolites containing or lacking BAS and/or silanols show that only the former complexes lead to active sites, ruling out confinement as the sole requirement for catalysis. The DFT calculations and ML-MD simulations also reveal that the surface-lanthanum coordination number is two (bidentate) for La BAS , with the metal center forming long, flexible bonds to two oxygen atoms bridging Si and Al, but only one oxygen atom (monodentate) for La SiO . The structure–activity relationship identifies confined, BAS-grafted species as active sites and provides important guidance for the design of enhanced atom-efficient catalysts.

Cui, Jinlei [Ames Laboratory, and Iowa State Univ.↗
DOE OSTI2023

Melting curves of ice polymorphs in the vicinity of the liquid–liquid critical point

The possible existence of a liquid–liquid critical point in deeply supercooled water has been a subject of debate due to the challenges associated with providing definitive experimental evidence. The pioneering work by Mishima and Stanley [Nature 392, 164–168 (1998)] sought to shed light on this problem by studying the melting curves of different ice polymorphs and their metastable continuation in the vicinity of the expected liquid–liquid transition and its associated critical point. Based on the continuous or discontinuous changes in the slope of the melting curves, Mishima [Phys. Rev. Lett. 85, 334 (2000)] suggested that the liquid–liquid critical point lies between the melting curves of ice III and ice V. Herein we explore this conjecture using molecular dynamics simulations with a machine learning model based on ab initio quantum-mechanical calculations. We study the melting curves of ices III, IV, V, VI, and XIII and find that all of them are supercritical and do not intersect the liquid–liquid transition locus. We also find a pronounced, yet continuous, change in the slope of the melting lines upon crossing of the liquid locus of maximum compressibility. Finally, we analyze the literature in light of our findings and conclude that the scenario in which the melting curves are supercritical is favored by the most recent computational and experimental evidence. Although the preponderance of evidence is consistent with the existence of a second critical point in water, the behavior of ice polymorph melting lines does not provide strong evidence in support of this viewpoint, according to our calculations.

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

Generalized representative structures for atomistic systems

A new method is presented to generate atomic structures that reproduce the essential characteristics of arbitrary material systems, phases, or ensembles. Previous methods allow one to reproduce the essential characteristics (e.g. the chemical disorder) of a large random alloy within a small crystal structure. The ability to generate small representations of random alloys, along with the restriction to crystal systems, results from using the fixed-lattice cluster correlations to describe structural characteristics. A more general description of the structural characteristics of atomic systems is obtained using complete sets of atomic environment descriptors. These are used within for generating representative atomic structures without restriction to fixed lattices. A general data-driven approach is provided here utilizing the atomic cluster expansion (ACE) basis. The N-body ACE descriptors are a complete set of atomic environment descriptors that span both chemical and spatial degrees of freedom and are used within for describing atomic structures. The generalized representative structure (GRS) method presented within generates small atomic structures that reproduce ACE descriptor distributions corresponding to arbitrary structural and chemical complexity. It is shown that systematically improvable representations of crystalline systems on fixed parent lattices, amorphous materials, liquids, and ensembles of atomic structures may be produced efficiently through optimization algorithms. With the GRS method, we highlight reduced representations of atomistic machine-learning training datasets that contain similar amounts of information and small 40–72 atom representations of liquid phases. The ability to use GRS methodology as a driver for informed novel structure generation is also demonstrated. The advantages over other data-driven methods and state-of-the-art methods restricted to high-symmetry systems are highlighted.

atomic cluster expansion↗
DOE OSTI2023

From Optimization to Sampling Through Gradient Flows

Optimization and sampling algorithms play a central role in science and engineering as they enable finding optimal predictions, policies, and recommendations, as well as expected and equilibrium states of complex systems. The notion of “optimality” is formalized by the choice of an objective function, while the notion of an “expected” state is specified by a probabilistic model for the distribution of states. Optimizing rugged objective functions and sampling multimodal distributions is computationally challenging, especially in high-dimensional problems. Here, for this reason, many optimization and sampling methods have been developed by researchers working in disparate fields such as Bayesian statistics, molecular dynamics, genetics, quantum chemistry, machine learning, weather forecasting, econometrics, and medical imaging.

Trillos, N. García↗
DOE OSTI2022

Stability and distortion of fcc LaH 10 with path-integral molecular dynamics

The synthesis of the high-temperature superconductor LaH 10 requires pressures in excess of 100 GPa, wherein it adopts a face-centered cubic structure. Upon decompression, this structure undergoes a distortion, which still supports superconductivity, but with a lower critical temperature. Previous calculations have shown that quantum and anharmonic effects are necessary to stabilize the cubic structure, but have not resolved the low pressure distortion. Using large scale path-integral molecular dynamics enabled by a machine learned potential, we show that a rhombohedral distortion appears at sufficiently low pressures. Here, we also highlight the importance of quantum zero-point motion in stabilizing the cubic structure.

08 HYDROGEN↗