Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interatomic potential energy”

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
NASA NTRS1974

Investigation of matter-antimatter interaction for possible propulsion applications

Matter-antimatter annihilation is discussed as a means of rocket propulsion. The feasibility of different means of antimatter storage is shown to depend on how annihilation rates are affected by various circumstances. The annihilation processes are described, with emphasis on important features of atom-antiatom interatomic potential energies. A model is developed that allows approximate calculation of upper and lower bounds to the interatomic potential energy for any atom-antiatom pair. Formulae for the upper and lower bounds for atom-antiatom annihilation cross-sections are obtained and applied to the annihilation rates for each means of antimatter storage under consideration. Recommendations for further studies are presented.

Morgan, D. L., Jr.↗
NASA NTRS1993

Condensed-matter energetics from diatomic molecular spectra

Analyses of molecular spectra and compression data from crystals show that a single function successfully describes the dependence on interatomic separation of both the potential energy of diatomic molecules and the cohesive binding energy of condensed matter. The empirical finding that one function describes interatomic energies for such diverse forms of matter and over a wide range of conditions can be used to extend condensed-matter equations of state but warrants further theoretical study.

Kim, In H.↗
NASA NTRS1993

Ab Initio Potential Energy Surface for H-H2

Ab initio calculations employing large basis sets are performed to determine an accurate potential energy surface for H-H2 interactions for a broad range of separation distances. At large distances, the spherically averaged potential determined from the calculated energies agrees well with the corresponding results determined from dispersion coefficients; the van der Waals well depth is predicted to be 75 +/- 3 micro E(h). Large basis sets have also been applied to reexamine the accuracy of theoretical repulsive potential energy surfaces (25-70 kcal/mol above the H-H2 asymptote) at small interatomic separations; the Boothroyd, Keogh, Martin, and Peterson (BKMP) potential energy surface is found to agree with results of the present calculations within the expected uncertainty (+/- 1 kcal/mol) of the fit. Multipolar expansions of the computed H-H2 potential energy surface are reported for four internuclear separation distances (1.2, 1.401, 1.449, and 1.7a(0)) of the hydrogen molecule. The differential elastic scattering cross section calculated from the present results is compared with the measurements from a crossed beam experiment.

Patridge, Harry↗
NASA NTRS1974

Kinetics of homogeneous nucleation in many component systems

Reiss's classical treatment of the kinetics of homogeneous nucleation in a system containing two chemical components is extended to many-component systems. The formulation is analogous to the pseudo-stationary state theory of chemical reaction rates with the free energy as a function of the composition of the embryo taking the place of the potential energy as a function of interatomic distances.

Hirschfelder, J. O.↗
NASA NTRS1974

Kinetics of homogeneous nucleation on many-component systems

Reiss's (1950) classical treatment of the kinetics of homogeneous nucleation in a system containing two chemical components is extended to many-component systems. The formulation is analogous to the pseudostationary-state theory of chemical reaction rates, with the free energy as a function of the composition of the embryo taking the place of the potential energy as a function of interatomic distances.

Hirschfelder, J. O.↗
NASA NTRS2001

Potential Energy Curves and Transport Properties for the Interaction of He with Other Ground-state Atoms

The interactions of a He atom with a heavier atom are examined for 26 different elements, which are consecutive members selected from three rows (Li - Ne, Na - Ar, and K,Ca, Ga - Kr) and column 12 (Zn,Cd) of the periodic table. Interaction energies are determined wing high-quality ab initio calculations for the states of the molecule that would be formed from each pair of atoms in their ground states. Potential energies are tabulated for a broad range of Interatomic separation distances. The results show, for example, that the energy of an alkali interaction at small separations is nearly the same as that of a rare-gas interaction with the same electron configuration for the dosed shells. Furthermore, the repulsive-range parameter for this region is very short compared to its length for the repulsion dominated by the alkali-valence electron at large separations (beyond about 3-4 a(sub 0)). The potential energies in the region of the van der Waals minimum agree well with the most accurate results available. The ab initio energies are applied to calculate scattering cross sections and obtain the collision integrals that are needed to determine transport properties to second order. The theoretical values of Li-He total scattering cross sections and the rare-gas atom-He transport properties agree well (to within about 1%) with the corresponding measured data. Effective potential energies are constructed from the ab initio energies; the results have been shown to reproduce known transport data and can be readily applied to predict unknown transport properties for like-atom interactions.

Partridge, Harry↗
NASA NTRS2004

Effect of anharmonicity of interatomic potential on strain distribution in semiconductor nanostructures

Experiments and theory have shown that the energy spectrum of nanostructures is extremely sensitive to the built-in strain. Knowledge of the strain distribution is therefore Experiments and theory have shown that the energy spectrum of nanostructures is extremely sensitive to the built-in strain. Knowledge of the strain distribution is therefore of utmost importance for the design of optical devices with prescribed light emission spectrum.

nanostructures↗
NASA NTRS1985

Statistical mechanics of light elements at high pressure. VIII - Thomas-Fermi-Dirac theory for binary mixtures of H with He, C, and O

We present three-dimensional Thomas-Fermi-Dirac calculations of lattice mixing energies of hydrogen with carbon and oxygen atoms, respectively. The results are used to derive effective interatomic potentials for use in liquid-state mixture calculations. We then use the potentials to derive analytic expressions for binary mixture-free energies and to map out the phase diagrams of mixtures of hydrogen with, respectively, helium, carbon, and oxygen, over a pressure range of about 5 to about 10 to the 3rd Mbar. Within this pressure range, all three of the latter elements are found to have unlimited solubility in metallic hydrogen over a temperature range which lies above their pure-element melting temperatures, and which includes likely interior temperatures in the Jovian planets.

Hubbard, W. B.↗
NASA NTRS2019

AladynPi – Adaptive Neural Network Molecular Dynamics Simulation Code with Physically Informed Potential: Computational Materials Mini-Application

This report provides an overview and description of commands used in the Computational Materials mini-application, AladynPi. AladynPi is an extension of a previously released mini-application, Aladyn (https://github.com/nasa/aladyn; Yamakov, V.I., and Glaessgen, E.H., NASA/TM-2018-220104). Aladyn and AladynPi are basic molecular dynamics codes written in FORTRAN 2003, which are designed to demonstrate the use of adaptive neural networks (ANNs) in atomistic simulations. The role of ANNs is to efficiently reproduce the very complex energy landscape resulting from the atomic interactions in materials with the accuracy of the more expensive quantum mechanics-based calculations. The ANN is trained on a large set of atomic structures calculated using the density functional theory method. An input for the ANN is a set of structure coefficients, characterizing the local atomic environment of each atom, for which the atomic energy is obtained in the ANN inference process. In Aladyn, the ANN gives directly the energy of interatomic interactions. In AladynPi, the ANN gives optimized parameters for a predefined empirical function, known as bond-order-potential (BOP). The parameterized BOP function is then used to calculate the energy. AladynPi code is being released to serve as a training testbed for students and professors in academia to explore possible optimization algorithms for parallel computing on multicore central processing unit (CPU) computers or computers utilizing manycore architectures based on graphic processing units (GPUs). The effort is supported by the High Performance Computing incubator (HPCi) project at NASA Langley Research Center.

Yamakov, Vesselin I.↗
NASA NTRS2021

Machine Learning the COSMO Model for Predicting Thermodynamics of Electrolyte Mixtures

Bottom-up design of electrolyte mixtures for battery systems requires predicting macro thermodynamic properties from molecular constituents. For instance, molten salt electrolyte batteries require conditions far above room temperature to operate. Therefore, discovering mixtures with increasingly lower eutectic melting points is desirable. A model that can approximate chemical activity is a valuable tool to search through the vast compositional design space. Machine learning can predict properties of materials such as vibrational free energies, electronic energy gaps, and thermal conductivities. Moreover, they can learn physical models such as interatomic potentials. The COSMO-SAC model uses theory and empirical parameterization to predict liquid-vapor and liquid-solid properties using first-principles calculations. However, obtaining activity coefficients required for parameterizing the COSMO-SAC model is costly and limited to a select chemical space. In this work, we explored if machine learning methods could improve the COSMO-SAC model and bridge density functional theory calculations to liquid phase thermodynamic properties. Our data-driven approach uses existing databases for sigma-profiles of organic solvents and reconciles their methodological differences via ensemble averaging. First, an optimal machine learning model is constructed for each dataset. Our machine learning algorithms use the sigma-profile as an input feature to predict binary mixtures' activity coefficients using multi-output regression. Each dataset uses different choices of functionals, methods, and basis sets. Therefore, our ensemble model attempts to predict corrected activity coefficients given the combination of all the model outputs. The activity coefficients used for training are generated using the COSMO-SAC model. This approach enables the extraction of meaningful information from the existing datasets to improve the COSMO-SAC model for obtaining thermodynamic properties of electrolyte mixtures. With the liquid phase activities, we can identify electrolyte mixtures that meet desired phase equilibria conditions.

Thermodynamics↗
NASA NTRS2017

The Statistical Mechanics of Solar Wind Hydroxylation at the Moon, Within Lunar Magnetic Anomalies, and at Phobos

We present a new formalism to describe the outgassing of hydrogen initially implanted by the solar wind protons into exposed soils on airless bodies. The formalism applies a statistical mechanics approach similar to that applied recently to molecular adsorption onto activated surfaces. The key element enabling this formalism is the recognition that the interatomic potential between the implanted H and regolith-residing oxides is not of singular value but possess a distribution of trapped energy values at a given temperature, F(U,T). All subsequent derivations of the outward diffusion and H retention rely on the specific properties of this distribution. We find that solar wind hydrogen can be retained if there are sites in the implantation layer with activation energy values exceeding 0.5eV. We especially examine the dependence of H retention applying characteristic energy values found previously for irradiated silica and mature lunar samples. We also apply the formalism to two cases that differ from the typical solar wind implantation at the Moon. First, we test for a case of implantation in magnetic anomaly regions where significantly lower-energy ions of solar wind origin are expected to be incident with the surface. In magnetic anomalies, H retention is found to be reduced due to the reduced ion flux and shallower depth of implantation. Second, we also apply the model to Phobos where the surface temperature range is not as extreme as the Moon. We find the H atom retention in this second case is higher than the lunar case due to the reduced thermal extremes (that reduces outgassing).

Solar Wind↗
NASA NTRS1978

Modified lattice-statics approach to dislocation calculations. II - Application

The atomic structure of a screw dislocation core of the 110 line type in aluminum is calculated by the modified lattice-statics method developed in the preceding paper. The method includes anharmonic as well as harmonic forces and permits relaxation of the atoms in all three dimensions. All forces used in the present calculations were derived from a first-principles interatomic pair potential obtained via pseudopotential theory. Several significant differences from the ordinary lattice statics results are noted, including the displacement field, Peierl's energy barrier, and the equilibrium core-center location.

Esterling, D. M.↗
NASA NTRS1990

A study of microkinetic adjustments required to match shock wave experiments and Monte Carlo Direct Simulation for a wide Mach number range

Shock wave thickness predictions from Monte Carlo Direct Simulations, using differential scattering and the Maitland-Smith-Aziz interatomic potential, underpredict experiments as shock Mach numbers increase above about 4. Examination of several sources of data has indicated that at relatively high energies the repulsive portion of accepted potentials such as the Maitland-Smith-Aziz may be too steep. An Exponential-6 potential due to Ross, based on high energy molecular beam scattering data and shock velocity measurements in liquid argon, has been combined with the lower energy portion of the Maitland-Smith-Aziz potential. When this hybrid potential is used in Monte Carlo Direct Simulations, agreement with experiments is improved over the previous predictions using the pure Maitland-Smith-Aziz form.

Pham-Van-diep, Gerald C.↗
NASA NTRS2002

Equivalent-Continuum Modeling With Application to Carbon Nanotubes

A method has been proposed for developing structure-property relationships of nano-structured materials. This method serves as a link between computational chemistry and solid mechanics by substituting discrete molecular structures with equivalent-continuum models. It has been shown that this substitution may be accomplished by equating the vibrational potential energy of a nano-structured material with the strain energy of representative truss and continuum models. As important examples with direct application to the development and characterization of single-walled carbon nanotubes and the design of nanotube-based devices, the modeling technique has been applied to determine the effective-continuum geometry and bending rigidity of a graphene sheet. A representative volume element of the chemical structure of graphene has been substituted with equivalent-truss and equivalent continuum models. As a result, an effective thickness of the continuum model has been determined. This effective thickness has been shown to be significantly larger than the interatomic spacing of graphite. The effective thickness has been shown to be significantly larger than the inter-planar spacing of graphite. The effective bending rigidity of the equivalent-continuum model of a graphene sheet was determined by equating the vibrational potential energy of the molecular model of a graphene sheet subjected to cylindrical bending with the strain energy of an equivalent continuum plate subjected to cylindrical bending.

Odegard, Gregory M.↗
NASA NTRS2000

Sputtering Erosion in the Ion Thruster

During the first phase of this research, the sputtering yields of molybdenum by low energy (100 eV and higher) xenon ions were measured by using the methods of secondary neutral mass spectrometry (SNMS) and Rutherford backscattering spectrometry (RBS). However, the measured sputtering yields were found to be far too low to explain the sputtering erosions observed in the long-duration tests of ion thrusters. The only difference between the sputtering yield measurement experiments and the ion thruster tests was that the later are conducted at high ion fluences. Hence, a study was initiated to investigate if any linkage exists between high ion fluence and an enhanced sputtering yield. The objective of this research is to gain an understanding of the causes of the discrepancies between the sputtering rates of molybdenum grids in an ion thruster and those measured from our experiments. We are developing a molecular dynamics simulation technique for studying low-energy xenon ion interactions with molybdenum. It is difficult to determine collision sequences analytically for primary ions below the 200 eV energy range where the ion energy is too low to be able to employ a random cascade model with confidence and it is too high to have to consider only single collision at or near the surface. At these low energies, the range of primary ions is about 1 to 2 nm from the surface and it takes less than 4 collisions on the average to get an ion to degrade to such an energy that it can no longer migrate. The fine details of atomic motion during the sputtering process are revealed through computer simulation schemes. By using an appropriate interatomic potential, the positions and velocities of the incident ion together with a sufficient number of target atoms are determined in small time steps. Hence, it allows one to study the evolution of damages in the target and its effect on the sputtering yield. We are at the preliminary stages of setting up the simulation program.

Ray, Pradosh K.↗
NASA NTRS1972

Optimal utilization of total elastic scattering cross section data for the determination of interatomic potentials

The problem of inversion is considered in relation to absolute total cross sections Q(v) for atom-atom collisions and their velocity dependence, and the glory undulations and the transition to high velocity behavior. There is a limit to the amount of information available from Q(v) even when observations of good accuracy (e.g., + or - 0.25%) are in hand over an extended energy range (from thermal energies upward by a factor of greater than 1000 in relative kinetic energy). Methods were developed for data utilization, which take full advantage of the accuracy of the experimental Q(v) measurements.

Bernstein, R. B.↗