Search NASA⌕ Search

SEARCH · Search NASA

Results for “Empirical potentials”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A comparative study of the effectiveness of empirical potentials for molecular dynamics simulations of borosilicate glasses

Despite their practical importance, atomistic modeling of B 2 O 3 containing glasses have been challenging due to the lack of reliable empirical potentials. Fortunately, a few recent developments have shown promises to simulate these glasses where the boron coordination has complex and non-linear dependence on glass composition. This work aims to provide an evaluation of the effectiveness of three recently developed potentials by a systematic study of a series (~20) of sodium borosilicate glasses with constant K (ratio of [SiO 2 ]/[B 2 O 3 ]) =2 and varying R (ratio of [Na 2 O]/[B 2 O 3 ] ranging from 0.1 to 4) values and several sodium boroaluminosilicate glasses. A comparative assessment was established on the basis of the short- and medium-range structure features, such as boron N 4 values, total correlation functions, bond angle distribution, oxygen speciation, and mechanical properties using experimental or well-established models as criteria. Furthermore, this work provides insights on the choices of empirical potentials for MD simulations of borosilicate glasses and gives directions of future potential development and refining.

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

Development of a semi-empirical potential for simulation of Ni solute segregation into grain boundaries in Ag

An Ag–Ni semi-empirical potential was developed to simulate the segregation of Ni solutes at Ag grain boundaries (GBs). The potential combines a new Ag potential fitted to correctly reproduce the stable and unstable stacking fault energies in this metal and the existing Ni potential from Mendelev et al (2012 Phil. Mag. 92 4454–69). The Ag–Ni cross potential functions were fitted to ab initio data on the liquid structure of the Ag 80 Ni 20 alloy to properly incorporate the Ag–Ni interaction at small atomic separations, and to the Ni segregation energies at different sites within a high-energy Σ9 (221) symmetric tilt GB. By deploying this potential with hybrid Monte Carlo/molecular dynamics simulations, it was found that heterogeneous segregation and clustering of Ni atoms at GBs and twin boundary defects occur at low Ni concentrations, 1 and 2 at%. This behavior is profoundly different from the homogeneous interfacial dispersion generally observed for the Cu segregation in Ag. A GB transformation to amorphous intergranular films was found to prevail at higher Ni concentrations (10 at%). In conclusion, the developed potential opens new opportunities for studying the selective segregation behavior of Ni solutes in interface-hardened Ag metals and its effect on plasticity.

36 MATERIALS SCIENCE↗
NASA NTRS2001

Ab Initio and Improved Empirical Potentials for the Calculation of the Anharmonic Vibrational States and Intramolecular Mode Coupling of N-Methylacetamide

The second-order Moller-Plesset ab initio electronic structure method is used to compute points for the anharmonic mode-coupled potential energy surface of N-methylacetamide (NMA) in the trans(sub ct) configuration, including all degrees of freedom. The vibrational states and the spectroscopy are directly computed from this potential surface using the Correlation Corrected Vibrational Self-Consistent Field (CC-VSCF) method. The results are compared with CC-VSCF calculations using both the standard and improved empirical Amber-like force fields and available low temperature experimental matrix data. Analysis of our calculated spectroscopic results show that: (1) The excellent agreement between the ab initio CC-VSCF calculated frequencies and the experimental data suggest that the computed anharmonic potentials for N-methylacetamide are of a very high quality; (2) For most transitions, the vibrational frequencies obtained from the ab initio CC-VSCF method are superior to those obtained using the empirical CC-VSCF methods, when compared with experimental data. However, the improved empirical force field yields better agreement with the experimental frequencies as compared with a standard AMBER-type force field; (3) The empirical force field in particular overestimates anharmonic couplings for the amide-2 mode, the methyl asymmetric bending modes, the out-of-plane methyl bending modes, and the methyl distortions; (4) Disagreement between the ab initio and empirical anharmonic couplings is greater than the disagreement between the frequencies, and thus the anharmonic part of the empirical potential seems to be less accurate than the harmonic contribution;and (5) Both the empirical and ab initio CC-VSCF calculations predict a negligible anharmonic coupling between the amide-1 and other internal modes. The implication of this is that the intramolecular energy flow between the amide-1 and the other internal modes may be smaller than anticipated. These results may have important implications for the anharmonic force fields of peptides, for which N-methylacetamide is a model.

Gregurick, Susan K.↗
NASA NTRS2016

Comparison of Quantum Mechanical and Empirical Potential Energy Surfaces and Computed Rate Coefficients for N2 Dissociation

Comparisons are made between potential energy surfaces (PES) for N2 + N and N2 + N2 collisions and between rate coefficients for N2 dissociation that were computed using the quasiclassical trajectory method (QCT) on these PESs. For N2 + N we compare the Laganas empirical LEPS surface with one from NASA Ames Research Center based on ab initio quantum chemistry calculations. For N2 + N2 we compare two ab initio PESs (from NASA Ames and from the University of Minnesota). These use different methods for computing the ground state electronic energy for N4, but give similar results. Thermal N2 dissociation rate coefficients, for the 10,000K-30,000K temperature range, have been computed using each PES and the results are in excellent agreement. Quasi-stationary state (QSS) rate coefficients using both PESs have been computed at these temperatures using the Direct Molecular Simulation of Schwartentruber and coworkers. The QSS rate coefficients are up to a factor of 5 lower than the thermal ones and the thermal and QSS values bracket the results of shock-tube experiments. We conclude that the combination of ab initio quantum chemistry PESs and QCT calculations provides an attractive approach for the determination of accurate high-temperature rate coefficients for use in aerothermodynamics modeling.

quasiclassical trajectory method↗
DOE OSTI2018

Implementation of metal-friendly EAM/FS-type semi-empirical potentials in HOOMD-blue: A GPU-accelerated molecular dynamics software

We present an implementation of EAM and FS interatomic potentials, which are widely used in simulating metallic systems, in HOOMD-blue, a software designed to perform classical molecular dynamics simulations using GPU accelerations. We first discuss the details of our implementation and then report extensive benchmark tests. We demonstrate that single-precision floating point operations efficiently implemented on GPUs can produce sufficient accuracy when compared against double-precision codes, as demonstrated in test simulations of calculations of the glass-transition temperature of Cu 64.5 Zr 35.5 , and pair correlation function of liquid Ni 3 Al. Our code scales well with the size of the simulating system on NVIDIA Tesla M40 and P100 GPUs. Compared with another popular software LAMMPS running on 32 cores of AMD Opteron 6220 processors, the GPU/CPU performance ratio can reach as high as 4.6. In conclusion, the source code can be accessed through the HOOMD-blue web page for free by any interested user.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Small uranium and oxygen interstitial clusters in UO 2 : An empirical potential study

The aggregation of defect clusters in UO 2 under irradiations creates dislocation loops and cavities, which in turn contribute to swelling of the fuels, leading to detrimental pellet-clad mechanical interactions. The kinetics of the defect clusters, however, have received less attention and are not well understood. As such, using two interatomic potentials, the Morelon potential and a many-body model by Cooper, Rushton, and Grimes (CRG), the energetics and kinetics of small uranium interstitial clusters, with and without bound oxygen interstitials, are investigated with molecular dynamics and accelerated molecular dynamics simulations. For the single bound anti-Schottky interstitial cluster (U i 2O i ), both potentials find the interstitialcy migration mechanism to be a low barrier migration mechanism. However, for the CRG potential, an oxygen Frenkel pair assisted interstitial cluster migration mechanism with a lower migration barrier is identified during accelerated molecular dynamics simulations. For the Morelon potential, the interstitial clusters have geometrically compact structures, and small interstitial clusters (two and three anti-Schottky clusters) are more mobile than the single bound anti-Schottky interstitial cluster; for the CRG potential, the interstitial clusters have geometrically non-compact structures, yielding a more complex energy landscape. The rapid migration of interstitials and interstitial clusters are important for the irradiation response of UO 2 and correlates well with the very low barriers reported in a few historic experimental studies.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2022

Molecular Dynamics of High Pressure Tin Phases I: Strength and deformation evaluations of empirical potentials [Slides]

Multi-phase problems have so many more unknowns, we’d like to have a tool to constrain some open questions related to microstructure and twin & dislocation behavior. We want an atomistic scale perspective on aspects of strength. Some multi-scale questions accessible to atomistic study: What lattice-specific behavior influences dislocation production/mobility and/or twinning? Do the phase transformations wipe-out, modify or preserve grain size and orientation? Does plastic strain reset at phase transition? If so under what conditions? Tin is the material chosen for the effort because it is non-hazardous and has multiple accessible solid phases at relatively low pressures.

36 MATERIALS SCIENCE↗
NASA NTRS2016

Comparison of Quantum Mechanical and Empirical Potential Energy Surfaces and Computed Rate Coefficients for N2 Dissociation

Physics-based modeling of hypersonic flows is predicated on the availability of chemical reaction rate coefficients and cross sections for the collisional processes. This approach has been built around the use of quantum mechanical calculations to describe the interaction between the colliding particles. In this approach a potential energy surface (PES) is computed by solving the electronic Schrödinger equation and collision cross sections are determined for that PES using classical, semiclassical or quantum mechanical scattering methods. The rate coefficients are computed by integrating the thermally weighted cross sections. State-to-state rate coefficients are determined by only integrating over a thermal distribution of collisional energies. Finally, thermal rate coefficients are determined by summation of the state-to-state rate coefficients for reactions of molecules in all relevant ro-vibrational energy levels. If the flow is in thermal non-equilibrium, the translational, vibrational and rotational energy modes can be represented in different ways: three unique temperatures can be used to describe the distributions, the populations of individual ro-vibrational energy levels can be determined by solving the Master Equation, or through the use of direct simulation in particle-based Monte Carlo sampling. The PES-to-rate coefficient approach had been proposed and attempted in the early days of digital computing, but it is only in the last 15 years that computer hardware and software have been up to the task of calculating accurate interatomic and intermolecular potentials.

Jaffe, Richard L.↗
DOE OSTI2021

Assessment of empirical interatomic potential to predict thermal conductivity in ThO 2 and UO 2

Computing vibrational properties of crystals in the presence of complex defects often necessitates the use of (semi-)empirical potentials, which are typically not well characterized for perfect crystals. In this study we explore the efficacy of a commonly used embedded-atomempirical interatomic potential for the U x Th 1- x O 2 system, to compute phonon dispersion, lifetime, and branch specific thermal conductivity. Our approach for ThO 2 involves using lattice dynamics and the linearized Boltzmann transport equation to calculate phonon transport properties based on second and third order force constants derived from the empirical potential and from first-principles calculations. For UO 2 , to circumvent the accuracy issues associated with first-principles treatments of strong electronic correlations, we compare results derived from the empirical interatomic potential to previous experimental results. It is found that the empirical potential can reasonably capture the dispersion of acoustic branches, but exhibits significant discrepancies for the optical branches, leading to overestimation of phonon lifetime and thermal conductivity. The branch specific conductivity also differs significantly with either first-principles based results (ThO 2 ) or experimental measurements (UO 2 ). These findings suggest that the empirical potential needs to be further optimized for robust prediction of thermal conductivity both in perfect crystals and in the presence of complex defects.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1979

Close-coupling calculations with an INDOX/1s static potential, semiclassical exchange, and a semi-empirical polarisation potential for electron-CO2 elastic scattering and rotational excitation

A semiempirical molecular-orbital method for modelling the effective potential for electron-molecule scattering is applied to elastic scattering and rotational excitation of CO2 at 20 eV impact energy. Agreement with experiment is reasonably good. The calculated rotationally summed integral cross section is 67.8 a sub 0 squared.

Onda, K.↗
NASA NTRSDate not supplied

Semi-Empirical Interatomic Potential for Large-Scale Molecular Dynamics Simulation of Metal-Oxide Systems

Classical molecular dynamics (MD) simulation can be applied to systems containing billions of atoms during times up to microseconds. However, utilization of a large-scale MD simulation requires reliable but computationally cheap interatomic potentials. In the case of metallic systems, embedded atom method (EAM) and Finnis-Sinclair (FS) potentials are the reasonable choices, but their development for multicomponent alloys is a challenge. Addition of oxygen atoms makes the problem of developing semi-empirical potentials even more difficult because the long-range Coulomb interaction cannot be neglected, and different atoms have different charges which vary during MD relaxation such that a charge equilibration procedure should be applied every MD step. As results researchers tend to avoid simulating metal-oxide systems. In this talk, a new Finnis-Sinclair (FS) Ni-Co-Cr potential will be presented. Special attention will be paid to reproducing of the element melting temperatures and clustering. Next, a charge transfer ionic potential (CTIP) developed to simulate the NiO properties will be presented. This potential utilizes the Ni FS potential and accounts for the Coulomb interaction in the vicinity of oxygen atoms. Fitting of the CTIP parameters and FS potential functions describing the interaction with oxygen atoms to the available experimental and ab initio data will be discussed. Results of MD simulation of interaction of NiO particles with dislocations in Ni will be shown.

molecular dynamics simulation↗
DOE OSTI2025

Parameterizing empirical interatomic potentials for predicting thermophysical properties via an irreducible derivative approach: the case of ThO 2 and UO 2

The accuracy of classical physical property predictions using molecular dynamics simulations is determined by the quality of the interatomic potentials. Here we introduce a training approach for empirical interatomic potentials (EIPs) which is well suited for capturing phonons and phonon-related properties. Our approach is based on direct comparisons of the second- and third-order irreducible derivatives (IDs) between an EIP and the Born–Oppenheimer potential within density functional theory (DFT) calculations. IDs fully exploit space group symmetry and allow for training without redundant information. We demonstrate the fidelity of our approach in the context of ThO 2 and UO 2 , where we optimize parameters of an embedded-atom method potential in addition to core–shell interactions. Our EIPs provide thermophysical properties in good agreement with DFT and outperform widely utilized EIPs for phonon dispersion and thermal conductivity predictions. Reasonable estimates of thermal expansion and formation energies of Frenkel pairs are also obtained.

empirical interatomic potential↗
DOE OSTI2025

Deep potential molecular dynamics simulations of low-temperature plasma-surface interactions

Machine learning approaches to potential generation for molecular dynamics (MD) simulations of low-temperature plasma-surface interactions could greatly extend the range of chemical systems that can be modeled. Empirical potentials are difficult to generalize to complex combinations of multiple elements with interactions that might include covalent, ionic, and metallic bonds. This work demonstrates that a specific machine learning approach, Deep Potential Molecular Dynamics (DeepMD), can generate potentials that provide a good model of plasma etching in the Si-Cl-Ar system. Comparisons are made between MD results using DeepMD models and empirical potentials, as well as experimental measurements. Pure Si properties predicted by the DeepMD model are in reasonable agreement with experimental results. Simulations of Si bombardment by Ar + ions demonstrate the ability of the DeepMD method to predict sputtering yields as well as the depth of the amorphous-crystalline interface. Etch yields as a function of flux ratio and ion energy for simultaneous Cl 2 and Ar + impacts are in good agreement with previous simulation results and experiment. Predictions of etch yields and etch products during plasma-assisted atomic layer etching of Si-Cl 2 -Ar are shown to be in good agreement with MD predictions using empirical potentials and with experiment. Finally, good agreement was also seen with measurements for the spontaneous etching of Si by Cl atoms at 300 K. Further, the demonstration that DeepMD can reproduce results from MD simulations using empirical potentials is a necessary condition to future efforts to extend the method to a much wider range of systems for which empirical potentials may be difficult or impossible to obtain.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Enhancing high-fidelity neural network potentials through low-fidelity sampling

The efficacy of neural network potentials (NNPs) critically depends on the quality of the configurational datasets used for training. Prior research using empirical potentials has shown that well-selected liquid–solid transitional configurations of a metallic system can be translated to other metallic systems. This study demonstrates that such validated configurations can be relabeled using density functional theory (DFT) calculations, thereby enhancing the development of high-fidelity NNPs. Training strategies and sampling approaches are efficiently assessed using empirical potentials and subsequently relabeled via DFT in a highly parallelized fashion for high-fidelity NNP training. Our results reveal that relying solely on energy and force for NNP training is inadequate to prevent overfitting, highlighting the necessity of incorporating stress terms into the loss functions. To optimize training involving force and stress terms, we propose employing transfer learning to fine-tune the weights, ensuring that the potential surface is smooth for these quantities composed of energy derivatives. This approach markedly improves the accuracy of elastic constants derived from simulations in both empirical potential-based NNPs and relabeled DFT-based NNPs. Overall, this study offers significant insights into leveraging empirical potentials to expedite the development of reliable and robust NNPs at the DFT level.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Guidelines for creating artificial neural network empirical interatomic potential from first-principles molecular dynamics data under specific conditions and its application to $α$-Ag 2 Se

First-principles molecular dynamics (FPMD) simulations are highly accurate, but due to their high calculation cost, the computational scale is often limited to hundreds of atoms and few picoseconds under specific temperature and pressure conditions. As such, we present here the guidelines for creating artificial neural network empirical interatomic potential (ANN potential) trained with such a limited FPMD data, which can perform long time scale MD simulations at least under the same conditions. The FPMD data for training are prepared on the basis of the convergence of radial distribution function [$g(r)$]. While training the ANN using total energy and atomic forces of the FPMD data, the error of pressure is also monitored and minimized. To create further robust potential, we add a small amount of FPMD data to reproduce the interaction between two atoms that are close to each other. ANN potentials for $α$-Ag 2 Se were created as an application example, and it has been confirmed that not only $g(r)$ and mean square displacements but also the specific heat requiring a long time scale simulation matched the FPMD and the experimental values. In addition, the MD simulation using the ANN potential achieved over 10 4 acceleration over the FPMD one. The guidelines proposed here mitigate the creation difficulty of the ANN potential, and a lot of FPMD data sleeping on the hard disk after the research may be put on the front stage again.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗