Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Inter-Atomic Potential”

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 OSTI2019

Performance Modeling of Vectorized SNAP Inter-Atomic Potentials on CPU Architectures.

SNAP potentials are inter-atomic potentials for molecular dynamics that enable simulations at accuracy levels comparable to density functional theory(DFT) at a fraction of the cost. As such, SNAP scales to on the order of 10 4 — 10 6 atoms. In this work, we explore CPU optimization of potentials computation using SIMD. We note that efficient use of SIMD is non-obvious as the application features an irregular iteration space for various potential terms, necessitating use of SIMD across atoms in a cross matrix, batched fashion. We present a preliminary analytic model to determine the correct batch size for several CPU architectures across several vendors, and show end-to-end speedups between 1.66x and 3.22x compared to the original.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Multi-reward reinforcement learning based development of inter-atomic potential models for silica

Abstract Silica is an abundant and technologically attractive material. Due to the structural complexities of silica polymorphs coupled with subtle differences in Si–O bonding characteristics, the development of accurate models to predict the structure, energetics and properties of silica polymorphs remain challenging. Current models for silica range from computationally efficient Buckingham formalisms (BKS, CHIK, Soules) to reactive (ReaxFF) and more recent machine-learned potentials that are flexible but computationally costly. Here, we introduce an improved formalism and parameterization of BKS model via a multireward reinforcement learning (RL) using an experimental training dataset. Our model concurrently captures the structure, energetics, density, equation of state, and elastic constants of quartz (equilibrium) as well as 20 other metastable silica polymorphs. We also assess its ability in capturing amorphous properties and highlight the limitations of the BKS-type functional forms in simultaneously capturing crystal and amorphous properties. We demonstrate ways to improve model flexibility and introduce a flexible formalism, machine-learned ML-BKS, that outperforms existing empirical models and is on-par with the recently developed 50 to 100 times more expensive Gaussian approximation potential (GAP) in capturing the experimental structure and properties of silica polymorphs and amorphous silica.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Atomistic measurement and modeling of intrinsic fracture toughness of two-dimensional materials

Quantifying the intrinsic mechanical properties of two-dimensional (2D) materials is essential to predict the long-term reliability of materials and systems in emerging applications ranging from energy to health to next-generation sensors and electronics. Currently, measurements of fracture toughness and identification of associated atomistic mechanisms remain challenging. Herein, we report an integrated experimental–computational framework in which in-situ high-resolution transmission electron microscopy (HRTEM) measurements of the intrinsic fracture energy of monolayer MoS 2 and MoSe 2 are in good agreement with atomistic model predictions based on an accurately parameterized interatomic potential. Changes in crystalline structures at the crack tip and crack edges, as observed in in-situ HRTEM crack extension tests, are properly predicted. Such a good agreement is the result of including large deformation pathways and phase transitions in the parameterization of the inter-atomic potential. The established framework emerges as a robust approach to determine the predictive capabilities of molecular dynamics models employed in the screening of 2D materials, in the spirit of the materials genome initiative. Moreover, it enables device-level predictions with superior accuracy (e.g., fatigue lifetime predictions of electro- and opto-electronic nanodevices).

Zhang, Xu↗
DOE OSTI2016

Dynamic stabilities of icosahedral-like clusters and their ability to form quasicrystals

The dynamic stabilities of the icosahedral-like clusters containing up to 2200 atoms are investigated for 15 metal elements. The clusters originate from five different initial structures (icosahedron, truncated decahedron, octahedron, closed-shell fragment of an HCP structure, and non-closed-shell fragment of an HCP structure). The obtained order of the dynamic stabilities of the icosahedral-like clusters can be assigned to three groups, from stronger to weaker, according to the size ranges involved: (Zr, Al, Ti) > (Cu, Fe, Co, Ni, Mg, Ag) > (Pb, Au, Pd, Pt, Rh, Ir), which correspond to the predicted formation ability of the quasicrystals. The differences of the sequences can be explained by analyzing the parameters of the Gupta-type many-body inter-atomic potentials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

ORNL’s Continuing Efforts in Evaluating and Validating TSLs [Slides]

This presentation finds that INS measurements are not as variable based on sample composition, because phonon spectra is a properly of the bulk material. Additionally, INS and transmission measurements are and should be the most significant part of TSL validation, as they are the fundamental quantities used by Monte Carlo neutronics codes. This presentation reasons that relaying solely on theoretical atomistic calculations for TSL evaluation is not enough, they are highly dependent on inter-atomic potentials (which are not the most accurate depending on the quantity trying to be reproduced) and evaluators knowledge and understanding of the material being studied. Finally on Critical benchmarks, while extremely useful for validation of nuclear data at all energies, critical benchmarks are not the best tool to provide a definitive answer on conflicting TSLs. The presentation concludes that INS and transmission measurements need to be the basis of validation of TSLs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Advancements in Validation of TSLs through Inelastic Neutron Scattering and Transmission Measurements [Slides]

This presentation finds that INS measurements are not as variable based on sample composition, because phonon spectra is a properly of the bulk material. Additionally, INS and transmission measurements are and should be the most significant part of TSL validation, as they are the fundamental quantities used by Monte Carlo neutronics codes. Relaying solely on theoretical atomistic calculations for TSL evaluation is not enough, they are highly dependent on inter-atomic potentials (which are not the most accurate depending on the quantity trying to be reproduced) and evaluators knowledge and understanding of the material being studied. When dealing with Crit. benchmarks - While extremely useful for validation of nuclear data at all energies, critical benchmarks are not the best tool to provide a definitive answer on conflicting TSLs. Finally, INS and transmission measurements need to be the basis of validation of TSLs.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Modeling graphene sheet growth and dynamical matrix calculations using molecular dynamics

Molecular dynamics (MD) has been an incredibly useful tool to model physical processes that were synthesized experimentally but not fully understood. MD, through the use of semi-empirical inter-atomic potentials, has allowed understanding of different physical processes in materials science. Yet as well as providing useful insights into materials science, molecular dynamics has a wider range of usability. In this report, I will be detailing how MD can be used to study graphene formation from a carbon liquid which requires high temperatures and pressures. Beyond this, I will describe the usefulness of MD for understanding the physics for phonon transport quantum sensors. To do this, MD was employed to determine the dynamical matrix by treating atoms as coupled oscillators. An accurate understanding of the dynamical matrix of a system is required to calculate the non-equilibrium Green’s function used to describe the phonon transport within phonon wave-guides. I found that, across multiple pressures and temperatures, randomly placed carbon atoms will show evidence of pent-first formation with semi-empirical models. Density functional theory (DFT), on the other hand, was too computationally expensive to use for full scale MD simulations, but we have the possibility of training a machine learned interatomic potential to approximate DFT for carbon in the environments being studied for pent-first graphene sheet formation.

36 MATERIALS SCIENCE↗
NASA NTRS1972

Lattice model calculation of elastic and thermodynamic properties at high pressure and temperature

The elastic constants and the entire frequency spectrum were calculated up to high pressure for the alkali halides in the NaCl lattice, based on an assumed functional form of the inter-atomic potential. The quasiharmonic approximation is used to calculate the vibrational contribution to the pressure and the elastic constants at arbitrary temperature. By explicitly accounting for the effect of thermal and zero point motion, the adjustable parameters in the potential are determined to a high degree of accuracy from the elastic constants and their pressure derivatives measured at zero pressure. The calculated Gruneisen parameter, the elastic constants and their pressure derivatives are in good agreement with experimental results up to about 600 K. The model predicts that for some alkali halides the Grunesen parameter may decrease monotonically with pressure, while for others it may increase with pressure, after an initial decrease.

Demarest, H. H., Jr.↗
DOE OSTI2016

Ab Initio Molecular Dynamics Analysis of Xenon Incorporation in Uranium Dioxide

The fission process in nuclear fuel produces volatile fission product gases, such as xenon, that migrate out of the fuel and into surrounding structures. Consequently, these fission product gases may affect thermal performance and integrity of the fuel. Fission gas release models, used in modern fuel performance codes such as BISON, are currently under investigation to improve their accuracy and ability to simulate the fission gas release phenomenon and its consequences. For uranium dioxide (UO{sub 2}), the most commonly used fuel in commercial power reactors, fission gas release models are known to have significant uncertainties. In the interest of improving such models, several complex diffusional mechanisms have been examined using a variety of techniques such as Nudge Elastic Band (NEB), utilizing both Density Functional Theory (DFT) and semi-empirical inter-atomic potential functions. One proposed diffusion mechanism in UO{sub 2} is through tetra-vacancy (i.e., 2U{sub V} + 2O{sub V}) structures, which is expected to provide a suitable means for the incorporation and diffusion of xenon in the fuel crystal. In this work, a model was developed to study the incorporation of xenon in UO{sub 2} via a tetra-vacancy using ab initio molecular dynamics (AIMD). This is a first principles technique that abandons the use empirical force fields as usually is the case of classical MD analysis. Furthermore, in contrast with static transition analysis tools such as NEB, AIMD is able to capture the dynamics of a given process. The translation of xenon within a tetra-vacancy has been investigated using ab-initio molecular dynamics. The transition frequency, barrier energy, and incorporation energy were extracted from the model, demonstrating its ability to capture a variety of parameters and its potential for use in studying fission gas diffusive behavior. Future refinements to the model will seek to develop additional insight into Xe and Kr mobility and will likely involve executing more elaborate AIMD simulations to investigate this phenomenon. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Final Report (October 2024): University of Tennessee, Knoxville (UTK) contribution to: FusMatML: Machine Learning Atomistic Modeling for Fusion Materials Collaborative Project led by Dr. Aidan Thompson, Sandia National Laboratory

The rapid growth of the field of Machine Learning Inter-Atomic Potentials (MLIAP) has lead to a profusion of methods, all of which have some similarity to each other, but each also restricted to particular design choices, often arrived at in a rather ad hoc fashion. Beyond anecdotal evidence, and some benchmarking studies on specific problems, little progress has been made in developing design principles for MLIAPs. The goal of this project is to use machine learning, data science, and uncertainty quantification methods to optimize the design choices for MLIAP.

Density functional theory, Helium and Hydrogen↗
DOE OSTI2017

Modeling material interfaces with hybrid adhesion method

A molecular dynamics simulation approach is presented to approximate layered material structures using discrete interatomic potentials through classical mechanics and the underlying principles of quantum mechanics. This method isolates the energetic contributions of the system into two pure material layers and an interfacial region used to simulate the adhesive properties of the diffused interface. The strength relationship of the adhesion contribution is calculated through small-scale separation calculations and applied to the molecular surfaces through an inter-layer bond criterion. By segregating the contributions into three regions and accounting for the interfacial excess energies through the adhesive surface bonds, it is possible to model each material with an independent potential while maintaining an acceptable level of accuracy in the calculation of mechanical properties. This method is intended for the atomistic study of the delamination mechanics, typically observed in thin-film applications. Therefore, the work presented in this paper focuses on mechanical tensile behaviors, with observations in the elastic modulus and the delamination failure mode. To introduce the hybrid adhesion method, we apply the approach to an ideal bulk copper sample, where an interface is created by disassociating the force potential in the middle of the structure. Various mechanical behaviors are compared to a standard EAM control model to demonstrate the adequacy of this approach in a simple setting. In addition, we demonstrate the robustness of this approach by applying it on (1) a Cu-Cu 2 O interface with interactions between two atom types, and (2) an Al-Cu interface with two dissimilar FCC lattices. These additional examples are verified against EAM and COMB control models to demonstrate the accurate simulation of failure through delamination, and the formation and propagation of dislocations under loads. Finally, the results conclude that by modeling the energy contributions of an interface using hybrid adhesion bonds, we can provide an accurate approximation method for studies of large-scale mechanical properties, as well as the representation of various delamination phenomena at the atomic scale.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Growth of bi- and tri-layered graphene on silicon carbide substrate via molecular dynamics simulation

Molecular dynamics (MD) simulation with simulated annealing method is used to study the growth process of bi- and tri-layered graphene on a 6H-SiC (0001) substrate via molecular dynamics simulation. Tersoff-Albe-Erhart (TEA) potential is used to describe the inter-atomic interactions among the atoms in the system. The formation temperature, averaged carbon-carbon bond length, pair correlation function, binding energy and the distance between the graphene formed and the SiC substrate are quantified. The growth mechanism, graphitization of graphene on the SiC substrate and characteristics of the surface morphology of the graphene sheet obtained in our MD simulation compare well to that observed in epitaxially grown graphene experiments and other simulation works.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Structures of 38-atom gold-platinum nanoalloy clusters

Bimetallic nanoclusters, such as gold-platinum nanoclusters, are nanomaterials promising wide range of applications. We perform a numerical study of 38-atom gold-platinum nanoalloy clusters, Au{sub n}Pt{sub 38−n} (0 ≤ n ≤ 38), to elucidate the geometrical structures of these clusters. The lowest-energy structures of these bimetallic nanoclusters at the semi-empirical level are obtained via a global-minimum search algorithm known as parallel tempering multi-canonical basin hopping plus genetic algorithm (PTMBHGA), in which empirical Gupta many-body potential is used to describe the inter-atomic interactions among the constituent atoms. The structures of gold-platinum nanoalloy clusters are predicted to be core-shell segregated nanoclusters. Gold atoms are observed to preferentially occupy the surface of the clusters, while platinum atoms tend to occupy the core due to the slightly smaller atomic radius of platinum as compared to gold’s. The evolution of the geometrical structure of 38-atom Au-Pt clusters displays striking similarity with that of 38-atom Au-Cu nanoalloy clusters as reported in the literature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Bose–Einstein condensation under the cubic–quintic Gross–Pitaevskii equation in radial domains

We study stationary and quasi-stationary solutions for the cubic–quintic Gross–Pitaevskii equation modeling Bose–Einstein condensates (BECs) in one, two, and three spatial dimensions under the assumption of radial symmetry with the BEC dynamics influenced by a confining potential. We consider both repulsive and attractive cubic interactions – corresponding respectively to repulsive and attractive two-body interactions – under similar frameworks in order to deduce the effects of the potentials in each case. We also carefully consider the role played by the quintic nonlinearity (modeling the strength of inter-atomic coupling) in modifying the solutions arising due to a purely cubic interaction term. In one spatial dimension, we obtain a variety of exact solutions in the zero-potential limit (including new periodic solutions which generalize known soliton solutions) as well as perturbation solutions for small amplitude confining potentials. For more general forms of the confining potential, we rely on numerical simulations, but these agree with the analytical results when the latter are valid. We also consider the limit where the quintic term dominates the cubic term (with such a limit relevant in the study of a Tonks–Girardeau gas). Under the assumption of radial symmetry, we also consider cylindrical (or, cigar-shaped) and spherical BECs. We consider the nonperturbative regime where either the potential or the amplitude of the solutions is large, obtaining various qualitative analytical results. When the kinetic energy term is small (relative to the nonlinearity and the confining potential), we recover the expected Thomas–Fermi approximation for the stationary solutions. Numerical simulations, under a variety of external confining potentials, are then used to understand the role these potentials play on the BEC solution structure for both the attractive and repulsive regimes. This assortment of analytical and numerical results allows us to better understand the structure of BECs modeled via cubic–quintic Gross–Pitaevskii equation under confining potentials with radial symmetry. Some of the results are natural extensions of those for the cubic Gross–Pitaevskii equation, while others appear new.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS1984

Theoretical study of electron correlation effects in transition metal dimers

Introduction of partially localized orbitals is shown to reduce the number of terms needed to describe the bonding in transition metal clusters. Using this formalism, it is possible to compute the various intra- and inter-atomic electron correlation contributions to the bond energy. Calculations demonstrate the relative importance of several kinds of electron correlation terms involving the 3p, 3d, and 4s electrons. Improved interaction potentials are obtained for the dimers V(2) and Cr(2) when additional correlation is added to the CAS SCF results of Walch, Bauschlicher, Roos, and Nelin (1983).

Das, G. P.↗
DOE OSTI2016

FP-LAPW calculations of the elastic, electronic and thermoelectric properties of the filled skutterudite CeRu{sub 4}Sb{sub 12}

We have investigated the electronic structure, elastic and thermoelectric properties of the filled skutterudite CeRu{sub 4}Sb{sub 12} using the density functional theory (DFT). The full potential linearized augmented plane wave (FP-LAPW) method within a framework of the generalized gradient approximation (GGA) approach is used to perform the calculations presented here. The electronic structure calculation suggests an indirect band gap semiconducting nature of the material with energy band gap of 0.08 eV. The analysis of the elastic constants at relaxed positions reveals the ductile nature of the sample material with covalent contribution in the inter-atomic bonding. The narrow band gap semiconducting nature with high value of Seebeck coefficient suggests the possibility of the thermoelectric application of the material. The analysis of the thermal transport properties confirms the result obtained from the energy band structure of the material with high thermopower and dimensionless figure of merit 0.19 at room temperature.

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

Probing periodic potential of crystals via strong-field re-scattering

Strong-field ionization and re-scattering phenomena have been used to image angstrom-scale structures of isolated molecules in the gas phase. These methods typically make use of the anisotropic response of the participating molecular orbital. Recently, an anisotropic strong-field response has also been observed in high-order harmonic generation (HHG) from bulk crystals. In a (100) cut magnesium oxide crystal, extreme ultraviolet high-harmonics are found to depend strongly on the crystal structure and inter-atomic bonding. Here, we extend these measurements to other two important crystal orientations: (111) and (110). We find that HHG from these orientations is also strongly anisotropic. The underlying dynamics is understood using a real-space picture, where high-harmonics are produced via coherent collision of strong-field driven electrons from the atomic sites, including from the nearest neighbor atoms. We find that harmonic efficiency is enhanced when semi-classical electron trajectories connect to the concentrated valence charge distribution regions around the atomic cores. Similarly, the efficiency is suppressed when the trajectories miss the atomic cores. Furthermore, these results further support the real-space picture of HHG with implications for retrieving the periodic potential of the crystal, if not the wavefunctions in three-dimensions.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS2002

Enabling Computational Nanotechnology through JavaGenes in a Cycle Scavenging Environment

A genetic algorithm procedure is developed and implemented for fitting parameters for many-body inter-atomic force field functions for simulating nanotechnology atomistic applications using portable Java on cycle-scavenged heterogeneous workstations. Given a physics based analytic functional form for the force field, correlated parameters in a multi-dimensional environment are typically chosen to fit properties given either by experiments and/or by higher accuracy quantum mechanical simulations. The implementation automates this tedious procedure using an evolutionary computing algorithm operating on hundreds of cycle-scavenged computers. As a proof of concept, we demonstrate the procedure for evaluating the Stillinger-Weber (S-W) potential by (a) reproducing the published parameters for Si using S-W energies in the fitness function, and (b) evolving a "new" set of parameters using semi-empirical tightbinding energies in the fitness function. The "new" parameters are significantly better suited for Si cluster energies and forces as compared to even the published S-W potential.

Globus, Al↗