Search NASA⌕ Search

SEARCH · Search NASA

Results for “INTERATOMIC FORCES”

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 OSTI2022

Structural and vibrational behavior of 2 H -NbSe 2 at high pressures

In this work, the high-pressure structural and vibrational properties of a layered transition metal dichalcogenide 2 H -NbSe 2 were investigated using single-crystal x-ray diffraction and Raman spectroscopy, demonstrating its structural stability up to 35 GPa. The lattice compressibility changes character from being highly anisotropic at low pressures to largely isotropic at high pressures. Concomitantly, the interatomic bonds demonstrate highly anisotropic compression behavior with the Se-Se interlayer bonds compressing by >20%, while the intramolecular Se-Se distance shows a nonmonotonic pressure dependence with a maximum at-12 GPa. The nearest-neighbor central force lattice vibrational model yields pressure dependencies of the interatomic forces in qualitative agreement with bond length compression, providing insight into the vibrational properties of 2 H -NbSe 2 at high pressures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Liquid Sorption-Enhanced Haber–Bosch Process

The use of a liquid sorbent in a traditional Haber-Bosch process enables significant improvements in energy efficiency and potential cost savings for arguably the most important chemical process on the planet. The approach presented in this report employs an incompressible liquid sorbent that absorbs and releases ammonia (NH 3 ) under specific conditions. To achieve this, we investigate reactions of ammonia and pure phosphoric acid (H 3 PO 4 , PA), which rapidly neutralize to form an equilibrated solution of monoammonium phosphate (MAP) and diammonium phosphate (DAP) that functions as a reversible and regenerable sorbent. Through intimate contact of the gas-phase Haber-Bosch reaction mixture with this liquid absorbent, complete equilibrium uptake may be achieved in an appropriately sized separator, and facile separation occurs through the use of independent liquid and gas phases. Following depressurization and release of the ammonia product, only the incompressible fluid needs to be repressurized and returned to the reactor. This study documents proof-of-concept absorption and desorption experiments carried out in 75 mL batch reactors, predominantly charged with precise MAP and DAP mixtures that equilibrate at process-relevant temperatures and pressures. We then assemble the first thermodynamic relationships that underlie this advantaged separation strategy, validated by reactive force field (ReaxFF) interatomic potential simulations, and benchmarked with traditional separation routes via process modeling and technoeconomic analysis. The scale of energy consumption in the century-old Haber-Bosch process is massive, and the elegant liquid sorption approach reported here offers opportunities to enhance its energy efficiency for the next frontier of ammonia synthesis.

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

Process Image Analysis using Big Data, Machine Learning, and Computer Vision

The development of algorithms for machine learning and data analysis for the 3013 MIS corrosion surveillance program is a collaborative effort by SRNL, USC and GT. For corrosion detection, LCM image data is extracted from large binary files, with software written to convert the data to physical attributes (i.e. height, color and grayscale values; all as functions of a location in a plane projection). The user interface for the software permits selective downloading of binary data and interrogation of attributes. User input thresholds are used to flag attributes of interest. Machine learning algorithms, developed for this application, are used to determine whether the features are the result of corrosion. To address the fundamental mechanisms of corrosion, machine learning algorithms are being developed to derive interatomic potential force-fields from ab-initio DFT calculations. The goal is to apply molecular modeling on a large enough scale to guide the design of resistant materials.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Self-Healing Ceramics and Mechanical Metamaterials

We performed large scale molecular dynamics (MD) simulations and reactive & quantum molecular dynamics simulations to investigate: (1) faceting, grain growth, and crack healing in ceramics; (2) stress-release mechanism and healing in fibers; and (3) mechanical properties of two-dimensional metamaterials. We also designed highly accurate interatomic potentials and forces using artificial neural networks (ANN) and first-principles molecular dynamics simulation data for use in our MD and reactive MD simulations. We demonstrated the efficacy of modeling quantum materials atomistic simulation data using unsupervised deep learning approaches on a quantum computer. These achievements are elaborated in the following sections.

36 MATERIALS SCIENCE↗
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 OSTI2025

Development and transferability of neural-network models for plasma-surface interactions

Plasma-surface interactions are increasingly critical to modern technologies; yet, accurate molecular dynamics simulations remain limited by the capabilities of interatomic potentials. Deep Potentials (DPs) promise to revolutionize the field by providing a systematic method for producing accurate interatomic potentials. The primary challenge of DP development is selecting a dataset, which efficiently spans the set of atomic environments one expects to encounter in the subsequent molecular dynamics simulations. The computational cost of density functional theory calculations, which are the typical basis for DP development, makes it impossible to directly verify the quality of a given DP. To address this challenge, we explore the development of a deep-learned interatomic potential, “DeepREBO,” trained to reproduce the behavior of the REBO2 empirical potential, enabling direct validation of training methodology and transferability. Using an active learning framework, we begin with a minimal dataset and iteratively expand it to train a Deep Potential-Smooth Edition model that faithfully reproduces REBO2 results for 25 eV hydrogen bombardment of diamond (001), a particularly challenging case. We show that small, carefully curated datasets can outperform large, unguided ones, with effective models requiring fewer than 15 000 snapshots. Subsequent transferability tests demonstrate that while DeepREBO generalizes well to diamond (111) surfaces, performance degrades for amorphous carbon or higher-energy impacts, highlighting the need for use-case-specific training data. We also evaluate methods to improve short-range repulsion. This study outlines best practices for training robust deep potentials and underscores the importance of dataset design for predictive plasma simulations.

Ab-initio molecular dynamics↗
DOE OSTI2020

Re-evaluation of the TSL for Yttrium Hydride

Yttrium hydride (YH x ) is of interest as a high-temperature moderator material in advanced nuclear reactor systems because of its superior ability to retain hydrogen at elevated temperatures. Thermal neutron scattering laws (TSL) for hydrogen bound in yttrium hydride (H-YH 2 ) and yttrium bound in yttrium hydride (Y-YH 2 ) were previously evaluated by Naval Nuclear Laboratory using the ab initio approach and released in ENDF/B-VIII.0. In that work, density functional theory, incorporating the generalized gradient approximation (GGA) for the exchange-correlation energy, was used to simulate the face-centered cubic structure of YH 2 and calculate the interatomic Hellmann-Feynman forces for a 2×2×2 supercell containing 96 atoms. Lattice dynamics calculations using PHONON were used to determine the phonon density of states. The calculated phonon density of states for H and Y in YH 2 were then used to prepare H-YH 2 and Y-YH 2 TSL evaluations, in the incoherent approximation, using the LEAPR module of NJOY. In addition, elastic scattering was assumed to be incoherent for both H and Y. While the incoherent elastic scattering approximation is appropriate for H-YH 2 , it introduces an undesirable approximation for Y-YH 2 . In this work, we re-evaluate the TSL for Y-YH 2 using FLASSH (Full Law Analysis Scattering System Hub). Y-YH 2 is evaluated using the FLASSH generalized coherent elastic scattering capability in order to capture the Bragg peaks associated with the YH 2 crystal structure which were neglected in the prior NJOY-based evaluation due to limitations in LEAPR. An experimental approach to validate the Y-YH 2 TSL using neutron transmission measurements is discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Nanoporous Tio2 Water Training Data

Data and input files used to train a Deep Potential (DP) model for the nanoporous TiO2-water interface. The DeepMD-kit code was used to train the DP. Information about data format, and how to use DeepMD-kit can be found at https://docs.deepmodeling.com/projects/deepmd/en/master/

08 HYDROGEN↗
DOE OSTI2026

Deriving effective electrode–ion interactions from free-energy profiles at electrochemical interfaces

Understanding ion adsorption at electrified metal–electrolyte interfaces is essential for accurate modeling of electrochemical systems. Here, in this study, we systematically investigate the free energy profiles of Na + , Cl − , and F − ions at the Au(111)–water interface using enhanced sampling molecular dynamics with both classical force fields and machine-learned interatomic potentials (MLIPs). Our classical metadynamics results reveal a strong dependence of predicted ion adsorption on the Lennard-Jones parameters, highlighting that—without due care—standard mixing rules can lead to qualitatively incorrect descriptions of ion–metal interactions. We present a systematic methodology for tuning the cross term LJ parameters to control adsorption energetics in agreement with more accurate models. As a surrogate for an ab initio model, we employed the recently released Universal Models for Atoms MLIP, which validates classical trends and displays strong specific adsorption for chloride, weak adsorption for fluoride, and no specific adsorption for sodium, in agreement with experimental and theoretical expectations. By integrating molecular-level adsorption free energies into continuum models of the electric double layer, we show that specific ion adsorption substantially alters the interfacial ion population, the potential of zero charge, and the differential capacitance of the system. Our results underscore the critical importance of force field parameterization and advanced interatomic potentials for the predictive modeling of ion-specific effects at electrified interfaces and provide a robust framework for bridging molecular simulations and continuum electrochemical models.

Roncoroni, Fabrice [Lawrence Berkeley National Lab↗
NASA NTRS1980

Time correlation functions in classical dense fluids

The velocity autocorrelation functions and memory functions of dense classical fluids may be directly obtained from the static radial distribution function g(r) in an approximate way. Following the Mori projection operator formalism, the memory functions may be related to the fluctuating force correlation. At low densities, these functions may be evaluated by following the trajectories of particle pairs in the interatomic potential. At higher densities, the force correlation functions can be evaluated approximately from particle pair trajectories via the potential of the mean force. The results for argon fluids at various densities and temperatures agree satisfactorily with the molecular dynamics and the Enskog values. The decrease of the diffusion coefficient with density is partly due to the nature of g(r) which reflects the stronger clustering of atoms at higher densities.

Tsang, T.↗
DOE OSTI2023

Insights into factors that affect non-Arrhenius migration of a simulated incoherent Σ3 grain boundary

Non-Arrhenius grain boundary migration, sometimes referred to as antithermal migration where temperature and GB velocity values are inversely related to each other, is examined in an incoherent twin Σ3 [111] 60° (11 8 5) nickel grain boundary. Molecular dynamics is used to simulate migration and examine the effect of various factors on the migration, including interatomic potential, system size, driving force, and variation of atomic grain boundary structure. A classical model for grain boundary migration, in its unsimplified form, is used to analyze the results. The grain boundaries exhibit migration mechanisms with very low apparent barrier heights to migration. As a result, the boundaries migrate quickly but exhibit a velocity saturation similar to that of dislocations. Furthermore, the various interatomic potentials exhibit different migration velocities, but their similarities suggest they all predict similar overall behaviors of migration. The variation of atomic structure in the same incoherent twin grain boundary leads to diverse behaviors with barrier heights that vary from non-Arrhenius to Arrhenius migration. Facet nucleation is confirmed not to be a requirement for this boundary based on an examination of simulation cell size; however, the presence and/or interaction between numerous facets does suggest a slowing and increased barrier height to migration for larger boundaries.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Active learning for SNAP interatomic potentials via Bayesian predictive uncertainty

Bayesian inference with a simple Gaussian error model is used to efficiently compute prediction variances for energies, forces, and stresses in the linear SNAP interatomic potential. Here, the prediction variance is shown to have a strong correlation with the absolute error over approximately 24 orders of magnitude. Using this prediction variance, an active learning algorithm is constructed to iteratively train a potential by selecting the structures with the most uncertain properties from a pool of candidate structures. The relative importance of the energy, force, and stress errors in the objective function is shown to have a strong impact upon the trajectory of their respective net error metrics when running the active learning algorithm. Batched training of different batch sizes is also tested against singular structure updates, and it is found that batches can be used to significantly reduce the number of retraining steps required with only minor impact on the active learning trajectory.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Effects of Nonequilibrium Atomic Structure on Ionic Diffusivity in LLZO: A Classical and Machine Learning Molecular Dynamics Study

To improve the performance of electrochemical devices, it is essential to understand the effects of nonequilibrium motifs in solids, such as grain boundaries, amorphous phases, and highly strained regions, on atomic-scale transport and stability. Molecular dynamics simulations are used to explore the combined effect of far-from-equilibrium atomic structures and the choice of interatomic potential on ionic diffusivity predictions for Li 7 La 3 Zr 2 O 12 (LLZO), a promising solid electrolyte for all-solid-state batteries. Amorphization and high strain are considered using both classical Buckingham interatomic potentials and machine learning force fields. Here we find that both crystalline expansion and amorphization tend to slow diffusion, although the different physical encodings in the two potentials impact the properties in different ways. We trace these variations to a combination of structural and transport factors, the contributions of which are deconvoluted computationally. Graph-based analysis reveals that the variations for amorphous LLZO arise from the connectivity of diffusion pathways within the predicted structures, which generally correlates with diffusivity and is notably higher for structures generated by the machine learning force fields. Our study provides additional insight into the relationship between atomic structure and diffusivity in LLZO, while also highlighting the need for care in choosing and validating potentials to simulate far from equilibrium structures.

25 ENERGY STORAGE↗
DOE OSTI2019

Data-driven material models for atomistic simulation

The central approximation made in classical molecular dynamics simulation of materials is the interatomic potential used to calculate the forces on the atoms. Great effort and ingenuity is required to construct viable functional forms and find accurate parametrizations for potentials using traditional approaches. Machine learning has emerged as an effective alternative approach to develop accurate and robust interatomic potentials. Beginning with a very general model form, the potential is learned directly from a database of electronic structure calculations and therefore can be viewed as a multiscale link between quantum and classical atomistic simulations. Risk of inaccurate extrapolation exists outside the narrow range of time and length scales where the two methods can be directly compared. Here, we use the spectral neighbor analysis potential (SNAP) and show how a fit can be produced with minimal interpolation errors which is also robust in extrapolating beyond training. To demonstrate the approach, we have developed a tungsten-beryllium potential suitable for the full range of binary compositions. Subsequently, large-scale molecular dynamics simulations were performed of high energy Be atom implantation onto the (001) surface of solid tungsten. The machine learned W-Be potential generates a population of implantation structures consistent with quantum calculations of defect formation energies. A very shallow ( < 2 nm ) average Be implantation depth is predicted which may explain ITER diverter degradation in the presence of beryllium.

36 MATERIALS SCIENCE↗
NASA NTRS1989

Universal aspects of brittle fracture, adhesion, and atomic force microscopy

This universal relation between binding energy and interatomic separation was originally discovered for adhesion at bimetallic interfaces involving the simple metals Al, Zn, Mg, and Na. It is shown here that the same universal relation extends to adhesion at transition-metal interfaces. Adhesive energies have been computed for the low-index interfaces of Al, Ni, Cu, Ag, Fe, and W, using the equivalent-crystal theory (ECT) and keeping the atoms in each semiinfinite slab fixed rigidly in their equilibrium positions. These adhesive energy curves can be scaled onto each other and onto the universal adhesion curve. The effect of tip shape on the adhesive forces in the atomic-force microscope (AFM) is studied by computing energies and forces using the ECT. While the details of the energy-distance and force-distance curves are sensitive to tip shape, all of these curves can be scaled onto the universal adhesion curve.

Banerjea, Amitava↗
DOE OSTI2026

User Manual - HydraGNN v5.0: Distributed Implementation of Multi-Tasking Graph Neural Networks

This document serves as the user manual for HydraGNN v5.0, a scalable graph neural network (GNN) architecture for simultaneous prediction of multiple target properties using multi-task learning (MTL). This version of HydraGNN has been developed primarily to support the development, training, and deployment of predictive graph-based deep learning (DL) models for atomistic materials modeling. HydraGNN is templated over 13 message-passing policies, including invariant models (GIN, PNA, PNAPlus, GAT, MFC, CGCNN, SAGE, SchNet, DimeNet) and equivariant models (EGNN, PNAEq, PAINN, MACE), and supports distributed training via distributed data parallelism (DDP), DeepSpeed, and Fully Sharded Data Parallelism (FSDP) on leadership-class supercomputers. Although HydraGNN can be applied to problems beyond atomistic materials modeling, its current use is confined to homogeneous graphs. Additional capabilities include machine-learned interatomic potentials with energy-conserving forces, General, Powerful, and Scalable Graph Transformer (GraphGPS) global attention, periodic boundary conditions, hyperparameter optimization, mixed-precision training, and uncertainty quantification.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Proton radiation effects in indium oxide using cascade molecular dynamics simulations

Metal oxide (MO) semiconductors, characterized by their wide band gaps and notable charge transport properties, are promising candidates for electronic applications in extreme environments, including near-Earth space. However, atomistic simulations of radiation–matter interactions in MOs remain challenging due to the limitations of existing interatomic potentials, which often fail to capture both the short-range repulsive forces essential for radiation damage modeling and the long-range electrostatic effects governing defect evolution. In this work, we develop a customized interatomic potential tailored for radiation damage simulations in indium oxide (In 2 O 3 ) as a model system, a representative MO material. Our potential integrates the Ziegler-Biersack-Littmark potential to accurately describe short-range interactions with Buckingham and Coulombic potentials to account for long-range forces. We perform molecular dynamics simulations of low-energy proton irradiation using this custom potential. We employ the primary knock-on atom (PKA) cascade method to study atomic displacements and primary defect formation. Simulations were conducted for 1 keV proton irradiation in four randomly chosen directions, and PKA-driven defect analyses at 5, 10, and 15 keV to examine the effects of direction and energy level on damage generation. Our results provide insight into the impact of irradiation direction and energy level on the cascade evolution and defect formation mechanisms.

Atomistic simulations↗
DOE OSTI2023

Stable Solid Molecular Hydrogen above 900 K from a Machine-Learned Potential Trained with Diffusion Quantum Monte Carlo

Here, we survey the phase diagram of high-pressure molecular hydrogen with path integral molecular dynamics using a machine-learned interatomic potential trained with quantum Monte Carlo forces and energies. Besides the HCP and C2/c–24 phases, we find two new stable phases both with molecular centers in the Fmmm–4 structure, separated by a molecular orientation transition with temperature. The high temperature isotropic Fmmm–4 phase has a reentrant melting line with a maximum at higher temperature (1450 K at 150 GPa) than previously estimated and crosses the liquid-liquid transition line around 1200 K and 200 GPa.

08 HYDROGEN↗