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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Learning interpretable surface elasticity properties from bulk properties via neural network equation learners

Surface elasticity is central to understanding the mechanics and stability of surfaces and interfaces. It is characterized by quantities such as surface tension, residual surface stress, and surface stiffness. However their analytical expressions are typically difficult to derive from atomistic data, and depend strongly on modeling choices. This work presents a neural network-based equation learner which combines customized activation functions and connection-based pruning to discover parsimonious, closed-form equations for surface elasticity from atomistic simulations. Applying the method to seven face-centered cubic (FCC) metals, our equation learner uncovers interpretable equations that describe both low-Miller index and high-Miller index surface properties, capturing long-tail property distributions accurately. The discovered expressions are decoupled into two components: a universal, geometry-driven orientation function, and material-specific baseline coefficients. We find that lower-order properties such as surface tension are fundamentally geometry dependent, while higher-order properties such as surface stress and elasticity show more complex geometry and material dependence. We also relate material dependent coefficients to bulk properties, forming a clear map from bulk material properties to surface elasticity. Overall, this approach demonstrates that interpretable neurosymbolic machine learning can bridge the gap between atomistic simulations and physical laws, enabling the discovery of generalizable structure–property relationships for materials science phenomena such as surface elasticity.

Equation learning↗
DOE OSTI2024

Neural network interatomic potential-driven analysis of phase stability in Ti–V alloys at the atomistic scale

The evolution of the ω phase in titanium–vanadium (Ti–V) alloys is critical for their mechanical properties, particularly in aerospace and biomedical applications. Here, this study employs a Rapid Artificial Neural Network (RANN) potential to model the ω phase evolution at the atomistic level, demonstrating a high degree of consistency with experimental observations, unlike the Modified Embedded Atom Method (MEAM), which fails to capture this phase transformation accurately. RANN simulations replicate key phenomena such as the nucleation of α precipitates at ω/β interfaces and accurate lattice orientations, enhancing our understanding of phase stability and transformation kinetics. The findings affirm that RANN potentials can significantly improve the prediction accuracy of complex material behaviors, offering a powerful tool for designing advanced materials with tailored properties such as solute effect in various stacking fault energies. This approach not only bridges the gap between theoretical predictions and empirical data but also sets a new direction for future research in materials science, emphasizing the integration of machine learning techniques in the development and optimization of new alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Machine learning maximized Anderson localization of phonons in aperiodic superlattices

Nanostructuring materials to achieve ultra-low lattice thermal conductivity has proven to be extremely attractive for numerous applications such as thermoelectric energy conversion. Anderson localization of phonons due to aperiodicity can reduce thermal conductivity in superlattices, but the lower limit of thermal conductivity remains elusive due to the prohibitively large design space. In this work, we demonstrate that an intuition-based manual search for aperiodic superlattice structures (random multilayers or RMLs) with the lowest thermal conductivity yields only a local minimum, while a genetic algorithm (GA) based approach can efficiently identify the globally minimum thermal conductivity by only exploring a small fraction of the design space. Our results show that this minimum value occurs at an average RML period that is, surprisingly, smaller than the period corresponding to the minimum SL thermal conductivity. Interestingly, above this critical period, scattering of incoherent phonons at interfaces is less, whereas below this period, the room for randomization becomes less, thus putting more coherent phonons out of Anderson localization and causing increased thermal conductivity. Moreover, the lower limit of the thermal conductivity occurs at a moderate rather than maximum randomness of the layer thickness. Our machine learning approach demonstrates a general process of exploring an otherwise prohibitively large design space to gain non-intuitive physical insights.

42 ENGINEERING↗
DOE OSTI2021

Microscopic and Macroscopic Characterization of Grain Boundary Energy and Strength in Silicon Carbide via Machine-Learning Techniques

Predicting the properties of grain boundaries poses a challenge because of the complex relationships between structural and chemical attributes both at the atomic and continuum scales. Grain boundary systems are typically characterized by parameters used to classify local atomic arrangements in order to extract features such as grain boundary energy or grain boundary strength. The present work utilizes a combination of high-throughput atomistic simulations, macroscopic and microscopic descriptors, and machine-learning techniques to characterize the energy and strength of silicon carbide grain boundaries. Additionally, a diverse data set of symmetric tilt and twist grain boundaries are described using macroscopic metrics such as misorientation, the alignment of critical low-index planes, and the Schmid factor, but also in terms of microscopic metrics, by quantifying the local atomic structure and chemistry at the interface. These descriptors are used to create random-forest regression models, allowing for their relative importance to the grain boundary energy and decohesion stress to be better understood. Results show that while the energetics of the grain boundary were best described using the microscopic descriptors, the ability of the macroscopic descriptors to reasonably predict grain boundaries with low energy suggests a link between the crystallographic orientation and the resultant atomic structure that forms at the grain boundary within this regime. For grain boundary strength, neither microscopic nor macroscopic descriptors were able to fully capture the response individually. However, when both descriptor sets were utilized, the decohesion stress of the grain boundary could be accurately predicted. These results highlight the importance of considering both macroscopic and microscopic factors when constructing constitutive models for grain boundary systems, which has significant implications for both understanding the fundamental mechanisms at work and the ability to bridge length scales.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Dissolving salt is not equivalent to applying a pressure on water

Abstract Salt water is ubiquitous, playing crucial roles in geological and physiological processes. Despite centuries of investigations, whether or not water’s structure is drastically changed by dissolved ions is still debated. Based on density functional theory, we employ machine learning based molecular dynamics to model sodium chloride, potassium chloride, and sodium bromide solutions at different concentrations. The resulting reciprocal-space structure factors agree quantitatively with neutron diffraction data. Here we provide clear evidence that the ions in salt water do not distort the structure of water in the same way as neat water responds to elevated pressure. Rather, the computed structural changes are restricted to the ionic first solvation shells intruding into the hydrogen bond network, beyond which the oxygen radial-distribution function does not undergo major change relative to neat water. Our findings suggest that the widely cited pressure-like effect on the solvent in Hofmeister series ionic solutions should be carefully revisited.

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

Exploring the formation of gold/silver nanoalloys with gas-phase synthesis and machine-learning assisted simulations

While nanoalloys are of paramount scientific and practical interest, the main processes leading to their formation are still poorly understood. Key structural features in the alloy systems, including the crystal phase, chemical ordering, and morphology, are challenging to control at the nanoscale, making it difficult to extend their use to industrial applications. Here, in this contribution, we focus on the gold/silver system that has two of the most prevalent noble metals and combine experiments with simulations to uncover the formation mechanisms at the atomic level. Nanoparticles were produced using a state-of-the-art inert-gas aggregation source and analyzed using transmission electron microscopy and energy-dispersive X-ray spectroscopy. Machine-learning-assisted molecular dynamics simulations were employed to model the crystallization process from liquid droplets to nanocrystals. Our study finds a preponderance of nanoparticles with five-fold symmetric morphology, including icosahedra and decahedra which is consistent with previous results on mono-metallic nanoparticles. However, we observed that gold atoms, rather than silver atoms, segregate at the surface of the obtained nanoparticles for all the considered alloy compositions. These segregation tendencies are in contrast to previous studies and have consequences on the crystallization dynamics and the subsequent crystal ordering. We finally showed that the underpinning of this surprising segregation dynamics is due to charge transfer and electrostatic interactions rather than surface energy considerations.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Liquid to crystal Si growth simulation using machine learning force field

Machine learning force field (ML-FF) has emerged as a potential promising approach to simulate various material phenomena for large systems with ab initio accuracy. However, most ML-FFs have been used to study the phenomena relatively close to the equilibrium ground states. In this work, we have studied a far from equilibrium system of liquid to crystal Si growth using ML-FF. Here, we found that our ML-FF based on ab initio decomposed atomic energy can reproduce all the aspects of ab initio simulated growth, from local energy fluctuations to transition temperatures, to diffusion constant, and growth rates. We have also compared the growth simulation with the Stillinger-Weber classical force field and found significant differences. A procedure is also provided to correct a systematic fitting bias in the ML-FF training process, which exists in all training models, otherwise critical results like transition temperature will be wrong.

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

An Entropy-Maximization Approach to Automated Training Set Generation for Interatomic Potentials

Machine learning-based interatomic potentials are currently garnering a lot of attention as they strive to achieve the accuracy of electronic structure methods at the computational cost of empirical potentials. Given their generic functional forms, the transferability of these potentials is highly dependent on the quality of the training set, the generation of which can be highly labor-intensive. Good training sets should at once contain a very diverse set of configurations while avoiding redundancies that incur cost without providing benefits. We formalize these requirements in a local entropy-maximization framework and propose an automated sampling scheme to sample from this objective function. We show that this approach generates much more diverse training sets than unbiased sampling and is competitive with hand-crafted training sets.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

When do short-range atomistic machine-learning models fall short?

In this work, we explore the role of long-range interactions in atomistic machine-learning models by analyzing the effects on fitting accuracy, isolated cluster properties, and bulk thermodynamic properties. Such models have become increasingly popular in molecular simulations given their ability to learn highly complex and multi-dimensional interactions within a local environment; however, many of them fundamentally lack a description of explicit long-range interactions. In order to provide a well-defined benchmark system with precisely known pairwise interactions, we chose as the reference model a flexible version of the Extended Simple Point Charge (SPC/E) water model. Our analysis shows that while local representations are sufficient for predictions of the condensed liquid phase, the short-range nature of machine-learning models falls short in representing cluster and vapor phase properties. These findings provide an improved understanding of the role of long-range interactions in machine learning models and the regimes where they are necessary.

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

Investigating 3,4-bis(3-nitrofurazan-4-yl)furoxan detonation with a rapidly tuned density functional tight binding model

In this work, we describe a machine learning approach to rapidly tune density functional tight binding models for the description of detonation chemistry in organic molecular materials. Resulting models enable simulations on the several 10s of ps scales characteristic to these processes, with “quantum-accuracy.” We use this approach to investigate early shock chemistry in 3,4-bis(3-nitrofurazan-4-yl)furoxan, a hydrogen-free energetic material known to form onion-like nanocarbon particulates following detonation. We find that the ensuing chemistry is significantly characterized by the formation of large C x N y O z species, which are likely precursors to the experimentally observed carbon condensates. Beyond utility as a means of investigating detonation chemistry, the present approach can be used to generate quantum-based reference data for the development of full machine-learned interatomic potentials capable of simulation on even greater time and length scales, i.e., for applications where characteristic time scales exceed the reach of methods including Kohn–Sham density functional theory, which are commonly used for reference data generation.

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

A deep potential model with long-range electrostatic interactions

Machine learning models for the potential energy of multi-atomic systems, such as the deep potential (DP) model, make molecular simulations with the accuracy of quantum mechanical density functional theory possible at a cost only moderately higher than that of empirical force fields. However, the majority of these models lack explicit long-range interactions and fail to describe properties that derive from the Coulombic tail of the forces. To overcome this limitation, we extend the DP model by approximating the long-range electrostatic interaction between ions (nuclei + core electrons) and valence electrons with that of distributions of spherical Gaussian charges located at ionic and electronic sites. The latter are rigorously defined in terms of the centers of the maximally localized Wannier distributions, whose dependence on the local atomic environment is modeled accurately by a deep neural network. In the DP long-range (DPLR) model, the electrostatic energy of the Gaussian charge system is added to short-range interactions that are represented as in the standard DP model. The resulting potential energy surface is smooth and possesses analytical forces and virial. Missing effects in the standard DP scheme are recovered, improving on accuracy and predictive power. By including long-range electrostatics, DPLR correctly extrapolates to large systems the potential energy surface learned from quantum mechanical calculations on smaller systems. We illustrate the approach with three examples: the potential energy profile of the water dimer, the free energy of interaction of a water molecule with a liquid water slab, and the phonon dispersion curves of the NaCl crystal.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Exploring solvation structure and transport behavior for rational design of advanced electrolytes for next generation of lithium batteries

The efficacy of electrolytes significantly affects battery performance, leading to the development of several strategies to enhance them. Despite this, the understanding of solvation structure remains inadequate. It is imperative to understand the structure–property–performance relationship of electrolytes using diverse techniques. Here, this review explores the recent advancements in electrolyte design strategies for high capacity, high-voltage, wide-temperature, fast-charging, and safe applications. To begin, the current state-of-the-art electrolyte design directions are comprehensively reviewed. Subsequently, advanced techniques and computational methods used to understand the solvation structure are discussed. Additionally, the importance of high-throughput screening and advanced computation of electrolytes with the help of machine learning is emphasized. Finally, future horizons for studying electrolytes are proposed, aimed at improving battery performance and promoting their application in various fields by enhancing the microscopic understanding of electrolytes.

25 ENERGY STORAGE↗
DOE OSTI2024

Interpretation of autoencoder-learned collective variables using Morse–Smale complex and sublevelset persistent homology: An application on molecular trajectories

Dimensionality reduction often serves as the first step toward a minimalist understanding of physical systems as well as the accelerated simulations of them. In particular, neural network-based nonlinear dimensionality reduction methods, such as autoencoders, have shown promising outcomes in uncovering collective variables (CVs). However, the physical meaning of these CVs remains largely elusive. In this work, we constructed a framework that (1) determines the optimal number of CVs needed to capture the essential molecular motions using an ensemble of hierarchical autoencoders and (2) provides topology-based interpretations to the autoencoder-learned CVs with Morse–Smale complex and sublevelset persistent homology. Furthermore, this approach was exemplified using a series of n-alkanes and can be regarded as a general, explainable nonlinear dimensionality reduction method.

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

Interacting dust grains in complex plasmas: Ion wake formation and the electric potential

Dust grains have been used as minimally invasive probes to determine plasma parameters including the plasma density, temperature, and electric field in a plasma discharge. However, the dust grains in a plasma generate local potential disturbances due to the collection of charge and the subsequent electrostatic interactions between the dust and charged plasma particles. Dust grains in close proximity to one another exhibit interesting non-reciprocal interactions and self-organize into structures such as one-dimensional filamentary chains, two-dimensional “zigzags,” and three-dimensional helices, among others. The formation of these structures suggests that although the dust grains may be less invasive than traditional plasma probes, the disturbance to the local plasma environment introduced by dust grains is non-trivial. Commonly used analytic forms of the electric potential describing complex plasmas have failed to resolve the near-dust region, and as a result are insufficient to provide insight about the formation of complex dust structures. Here, we use an N-body simulation to compute the electric potential from ion densities near various dust grain configurations. We provide an alternative description to the standard analytic model for the electric potential of dust and ion wakes based on a Gaussian shaped cloud of ions. The electric potential obtained from simulations is used to identify minimum energy configurations for two and three dust grains. It is further demonstrated that the minimum potential region identified for N dust grains and their associated ion wakes does not predict the minimum-energy configuration of N + 1 dust grains.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗