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

Dynamic Instability of Solid Surfaces under Load

The mechanism of dynamic instability of a solid surface under a load, determined by perturbations of the electron density and the change in the interatomic interaction, is considered. This instability is manifested upon the excitation of dynamic displacements of atoms in the surface layer, which can lead to formation of finite-lifetime structures in the form of waves with a large amplitude. Such structures were earlier observed experimentally on a stretched germanium (111) surface.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Structural, Elastic, Electronic Properties and Interatomic Interactions in Metallic Tetraboride Series MB{sub 4} (M = Ru, Rh, Pd, Os, Ir, and Pt) Obtained from FLAPW–GGA Calculations

The report presents the results of systematic first-principle FLAPW–GGA calculations of a series of metallic tetraborides MB{sub 4} (where M = Ru, Rh, Pd, Os, Ir, and Pt) with space groups Pmmn, Immm, R-3m, and Pnnm. Their equilibrium structural parameters, cohesion and formation energies, elastic constants, anisotropies, Vickers microhardnesses, Debye temperatures, sound velocities, as well as electronic properties and interatomic interactions are estimated and analyzed as compared to each other and the available theoretical data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Simulation of the NaGd(MoO{sub 4}){sub 2}–NaEu(MoO{sub 4}){sub 2} and Na{sub 2}Gd{sub 4}(MoO{sub 4}){sub 7}–Na{sub 2}Eu{sub 4}(MoO{sub 4}){sub 7} Solid Solutions by the Interatomic Potential Method

Simulation of the solid solutions in the system of double sodium–gadolinium and sodium–europium molybdates, which are promising matrices for solid state lasers and phosphors has been carried out by the method of interatomic potentials. Two types of solid solutions have been studied, one of which contains finite components corresponding to the stoichiometric NaGd(MoO{sub 4}){sub 2}–NaEu(MoO{sub 4}){sub 2} compositions with statistical distribution of cations in the crystal lattice. Another object is a cation-deficient Na{sub 2}Gd{sub 4}(MoO{sub 4}){sub 7}–Na{sub 2}Eu{sub 4}(MoO{sub 4}){sub 7} system, in which we have examined the variants of statistical distribution and partial ordering of cations over structural positions. Atomistic simulation has been performed using the GULP 4.0.1 software package (General Utility Lattice Program). It is shown that when we pass from sodium-gadolinium molybdate to sodium-europium molybdate, both of stoichiometric and cation-deficient compositions, an increase in the unit cell volume is observed, while the density of the crystal, the energy of interatomic interactions in the structure, the vibrational entropy and the heat capacity decrease along with increasing europium content. The energy of interatomic interactions in the structure for cation-deficient solid solutions is less than for stoichiometric ones. Other aforementioned characteristics for cation-deficient solid solutions have greater values than for stoichiometric ones. The role of cluster europium centers in concentration quenching in NaGd(MoO{sub 4}){sub 2}–NaEu(MoO{sub 4}){sub 2} solid solutions has been examined.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Dynamic Spatial and Temporal Structures on the Surface of Solids under Load

A model and a method are proposed to explain the origin and main specific features of the dynamic patterns of different types, which were previously observed on the free surface of solids under load in experiments. It is shown that the pattern is formed due to the dynamic instability of the flat surface of the solid under load. This instability develops at the relaxation mechanism caused by the dynamic atomic displacements due to a change in the interatomic interaction during the electron-density redistribution.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Changes of the Thermodynamic Properties at Isochoric and Isobaric Decrease of the Silicon Nanocrystal Size

Equation of state P(ν/ν{sub o}) and the baric dependences of the lattice and surface properties of silicon macro- and nanocrystals have been calculated using the method of calculation of crystal properties from the pair Mie–Lennard-Jones interatomic potential and the RP-model of nanocrystal. The isothermal dependences of P(ν/ν{sub o}) for the macro- and the nanocrystal are shown to be intersected at a certain value of relative volume (ν/ν{sub o}){sub 0}. The surface pressure becomes zero at the intersection point (at (ν/ν{sub o}){sub 0}). The value of (ν/ν{sub o}){sub 0} decreases upon isomorphic–isomeric increase in temperature and also at isomorphic–isothermic decrease in the number of atoms N in the nanocrystal, or at isomeric–isothermic deviation of the nanocrystal shape from the most energetically optimal shape (in the RP-model, this shape is a cube). The obtained equation of state is used to study the changes of the silicon properties at isochoric (ν/ν{sub o} = 1) and also isobaric (P = 0) decrease in N at temperatures 300 and 1000 K.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

The Model of Decomposition of a Fe–Cu Alloy with Concentration-Depending Interatomic Interactions

The thermodynamics and kinetics of decomposition of a Fe–Cu alloy are investigated in the context of a simple ab initio parameterization model taking into account the concentration dependence of the magnetic contribution to the free energy. It is shown that taking into account a difference between the short-range and long-range magnetic orders near the Curie temperature is very important for calculating the solubility of copper in iron. The solubility of copper in bcc iron and the stability limit of a homogeneous state (physical spinodal) are evaluated via the kinetic Monte Carlo method. The impact of structural defects (dislocations and grain boundaries) on these curves, as well as on the kinetic time–temperature—transformation (TTT) diagram, is also discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Excitation of Soliton-Type Waves in Crystals of the A{sub 3}B Stoichiometry

Using the method of molecular dynamics and taking Ni{sub 3}Al and Pt{sub 3}Al as examples, crystals of the A{sub 3}B composition are considered for the possibility of excitation of soliton-type waves in them. The potentials obtained by the embedded-atom method were used to describe interatomic interactions. It is shown that the harmonic external stimulus can excite waves of the soliton type in a Pt{sub 3}Al crystal, but not in Ni{sub 3}Al. Such compression–expansion waves are generated because of excitation of discrete breathers with soft nonlinearity that cannot exist in a Ni{sub 3}Al crystal near the affected region. The detected waves are capable of propagating to thousands of nanometers along the Pt{sub 3}Al crystal without losses of integrity and speed. The shape of the obtained wave corresponds to the kink solution of the sine-Gordon equation. The aggregate amount of energy transferred by a wave is determined by the number of rows of atoms involved in fluctuations; this may involve dozens and hundreds of electron volts.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Thermal transport in UO 2 with defects and fission products by molecular dynamics simulations

The importance of the thermal transport in nuclear fuel has motivated a wide range of experimental and modelling studies. In this report, the reduction of thermal transport in UO 2 due to defects and fission products has been investigated using non-equilibrium MD simulations, with two sets of empirical potentials for studying the degregation of UO 2 thermal conductivity including a Buckingham type interatomic potential and a recently developed EAM type interatomic potential. Additional parameters for U 5+ and Zr 4+ in UO 2 have been developed for the EAM potential. The thermal conductivity results from MD simulations are then corrected for the spin-phonon scattering through Callaway model formulations. To validate the modelling results, comparison was made with experimental measurements on single crystal hyper-stoichiometric UO 2+x samples.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2015

Quantum mechanical studies of carbon structures

Carbon nanostructures, such as nanotubes and graphene, are of considerable interest due to their unique mechanical and electrical properties. The materials exhibit extremely high strength and conductivity when defects created during synthesis are minimized. Atomistic modeling is one technique for high resolution studies of defect formation and mitigation. To enable simulations of the mechanical behavior and growth mechanisms of C nanostructures, a high-fidelity analytical bond-order potential for the C is needed. To generate inputs for developing such a potential, we performed quantum mechanical calculations of various C structures.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2016

Milestone report on MD potential development for uranium silicide

This report summarizes the progress on the interatomic potential development of triuranium-disilicide (U 3 Si 2 ) for molecular dynamics (MD) simulations. The development is based on the Tersoff type potentials for single element U and Si. The Si potential is taken from the literature and a Tersoff type U potential is developed in this project. With the primary focus on the U 3 Si 2 phase, some other U-Si systems such as U 3 Si are also included as a test of the transferability of the potentials for binary U-Si phases. Based on the potentials for unary U and Si, two sets of parameters for the binary U-Si system are developed using the Tersoff mixing rules and the cross-term fitting, respectively. The cross-term potential is found to give better results on the enthalpy of formation, lattice constants and elastic constants than those produced by the Tersoff mixing potential, with the reference data taken from either experiments or density functional theory (DFT) calculations. In particular, the results on the formation enthalpy and lattice constants for the U 3 Si 2 phase and lattice constants for the high temperature U 3 Si (h-U 3 Si) phase generated by the cross-term potential agree well with experimental data. Reasonable agreements are also reached on the elastic constants of U 3 Si 2 , on the formation enthalpy for the low temperature U 3 Si (m-U 3 Si) and h-U 3 Si phases, and on the lattice constants of m-U 3 Si phase. All these phases are predicted to be mechanically stable. The unary U potential is tested for three metallic U phases (α, β, γ). The potential is found capable to predict the cohesive energies well against experimental data for all three phases. It matches reasonably with previous experiments on the lattice constants and elastic constants of αU.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Plasma Interactions with Mixed Materials and Impurity Transport

The project brings together three discipline areas at LLNL to develop advanced capability to predict the impact of plasma/material interactions (PMI) on metallic surfaces in magnetic fusion energy (MFE) devices. These areas are (1) modeling transport of wall impurity ions through the edge plasma to the core plasma, (2) construction of a laser blow-off (LBO) system for injecting precise amounts of metallic atoms into a tokamak plasma, and (3) material science analysis of fundamental processes that modify metallic surfaces during plasma bombardment. The focus is on tungsten (W), which is being used for the ITER divertor and in designs of future MFE devices. In area (1), we have worked with the University of California, San Diego (UCSD) on applications of the UEDGE/DUSTT coupled codes to predict the influx of impurity ions from W dust through the edge plasma, including periodic edge-plasma oscillations, and revived a parallel version of UEDGE to speed up these simulations. In addition, the impurity transport model in the 2D UEDGE code has been implemented into the 3D BOUT++ turbulence/transport code to allow fundamental analysis of the impact of strong plasma turbulence on the impurity transport. In area (2), construction and testing of the LBO injection system has been completed. The original plan to install the LBO on the National Spherical Torus Experiment Upgrade (NSTX-U) at Princeton and its use to validate the impurity transport simulations is delayed owing to NSTX-U being offline for substantial magnetic coil repair period. In area (3), an analytic model has been developed to explain the growth of W tendrils (or fuzz) observed for helium-containing plasmas. Molecular dynamics calculations of W sputtering by W and deuterium (D) ions shows that a spatial blending of interatomic potentials is needed to describe the near-surface and deeper regions of the material.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Learning-based Computational Study of the Thermodynamic, Structural, and Dynamic Properties of Molten Salt at the Atomic and Electronic Scale and Experimental Validations (Final Report for Project 19-16909)

The project completed several milestones related to the work scope to determine thermophysical properties of the molten salt FLiNaK. A combined approach to integrate advanced computations with neutron scattering experimental interrogation was used to quantify the structure and mass transport of pure FLiNaK and FLiNaK with additives. The additives investigated included Ce, Cs, Zr, and Cr. The primary outcome of the computational work was the development of a neutral network force field interatomic potential (referred to as the NNFF model in this document). The computational work also involved ab initio calculations to determine the ground truth and molecular dynamics simulations to provide computational analogs to the experimental campaign. The experimental work used enriched F 7 LiNaK prepared by Argonne National Laboratory. The use of enriched F 7 LiNaK reduced parasitic neutron absorption and significantly improved data collection and data quality. Several graduate students were supported during the course of the project, which included a one year no-cost extension. These students were trained in advanced computational methods, advanced neutron scattering methods, neutron scattering data analysis and interpretation, and molten salt sample preparation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Machine-Learned Force Field for Molecular Dynamics Simulations of Nonequilibrium Ammonia Synthesis on Iron Catalysts

Ammonia (NH 3 ) is one of the most important industrial chemicals. The conventional NH 3 synthesis method-the Haber–Bosch process-converts atmospheric nitrogen (N 2 ) into NH 3 using H 2 with an iron (Fe) catalyst. However, this process requires high pressures (100–200 atm) and temperatures (700–800 K) near thermal equilibrium. Recently, Fe-based nanocatalysts have been reported to produce promising NH 3 yields under atmospheric pressures and temperature-modulated nonequilibrium conditions. Understanding the mechanism of nonequilibrium catalysis with programmed temperature variation could help to optimize this fully electrified and less energy-intensive process. Although reactive molecular dynamics (RMD) simulations can be a useful tool to model nonequilibrium catalytic processes, they require the development of accurate force fields (i.e., interatomic potentials). Here, we present a machine-learned (ML) force field within the Deep Potential MD (DPMD) framework, trained using periodic density functional theory (DFT) calculations, to model NH 3 synthesis on Fe catalysts with various surface adsorbates such as *N, *H, *N 2 , *H 2 , *NH, *NH 2 , and *NH 3 . Here, we generated the DFT data from static models of elementary reactions on the most stable (110) surface of body-centered cubic Fe, which then were augmented by data from constant number of particles–volume–temperature (NVT) DFT-MD trajectories at various temperatures. Finally, we utilized the fully optimized ML force field to investigate reaction dynamics at an Fe(110) surface at linearly increasing temperatures using NVT-DPMD simulations. Our simulations indicate that pulsed temperature ramping could prove favorable for NH3 synthesis. For example, we conducted ramping under multiple sets of conditions: (i) from 900 to 1200 K over periods of 0.1–0.3 ns for Fe surfaces precovered with N or NH along with H; and (ii) from 300 to 600 K over 0.1–0.3 ns for Fe surfaces precovered with NH 3 . While our simulations so far are limited to short time scales (very rapid heating), these observations shed light on the mechanism of the high NH 3 synthesis rate achieved in a novel temperature-modulated nonequilibrium catalytic reactor using pulsed heating and cooling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 OSTI2016

The ReaxFF reactive force-field: Development, applications, and future directions

The reactive force-field (ReaxFF) interatomic potential is a powerful computational tool for exploring, developing and optimizing material properties. Methods based on the principles of quantum mechanics (QM), while offering valuable theoretical guidance at the electronic level, are often too computationally intense for simulations that consider the full dynamic evolution of a system. Alternatively, empirical interatomic potentials that are based on classical principles require significantly fewer computational resources, which enables simulations to better describe dynamic processes over longer timeframes and on larger scales. Such methods, however, typically require a predefined connectivity between atoms, precluding simulations that involve reactive events. The ReaxFF method was developed to help bridge this gap. Approaching the gap from the classical side, ReaxFF casts the empirical interatomic potential within a bond-order formalism, thus implicitly describing chemical bonding without expensive QM calculations. As a result, this article provides an overview of the development, application, and future directions of the ReaxFF method.

97 MATHEMATICS AND COMPUTING↗
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↗