Search NASA⌕ Search

SEARCH · Search NASA

Results for “Interatomic 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 271 records · Page 15

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.↗

Molecular Dynamics Simulation of a Multi-Walled Carbon Nanotube Based Gear

We used molecular dynamics to investigate the properties of a multi-walled carbon nanotube based gear. Previous work computationally suggested that molecular gears fashioned from (14,0) single-walled carbon nanotubes operate well at 50-100 gigahertz. The gears were formed from nanotubes with teeth added via a benzyne reaction known to occur with C60. A modified, parallelized version of Brenner's potential was used to model interatomic forces within each molecule. A Leonard-Jones 6-12 potential was used for forces between molecules. The gear in this study was based on the smallest multi-walled nanotube supported by some experimental evidence. Each gear was a (52,0) nanotube surrounding a (37,10) nanotube with approximate 20.4 and 16,8 A radii respectively. These sizes were chosen to be consistent with inter-tube spacing observed by and were slightly larger than graphite inter-layer spacings. The benzyne teeth were attached via 2+4 cycloaddition to exterior of the (52,0) tube. 2+4 bonds were used rather than the 2+2 bonds observed by Hoke since 2+4 bonds are preferred by naphthalene and quantum calculations by Jaffe suggest that 2+4 bonds are preferred on carbon nanotubes of sufficient diameter. One gear was 'powered' by forcing the atoms near the end of the outside buckytube to rotate to simulate a motor. A second gear was allowed to rotate by keeping the atoms near the end of its outside buckytube on a cylinder. The ends of both gears were constrained to stay in an approximately constant position relative to each other, simulating a casing, to insure that the gear teeth meshed. The stiff meshing aromatic gear teeth transferred angular momentum from the powered gear to the driven gear. The simulation was performed in a vacuum and with a software thermostat. Preliminary results suggest that the powered gear had trouble turning the driven gear without slip. The larger radius and greater mass of these gears relative to the (14,0) gears previously studied requires a smaller rotation rate and multiple rows of teeth to avoid excessive force on the gear teeth resulting, in slip and failure of the driven gear to turn. We hope that studies such as these will eventually lead to synthesis of components that can be assembled into atomically precise fullerene machines. These machines, in turn, may someday be used in machine-phase fullerene materials with remarkable properties.

Han, Jie↗

Molecular Dynamics Simulation of Carbon Nanotube Based Gears

We used molecular dynamics to investigate the properties and design space of molecular gears fashioned from carbon nanotubes with teeth added via a benzyne reaction known to occur with C60. A modified, parallelized version of Brenner's potential was used to model interatomic forces within each molecule. A Leonard-Jones 6-12 potential was used for forces between molecules. One gear was powered by forcing the atoms near the end of the buckytube to rotate, and a second gear was allowed.to rotate by keeping the atoms near the end of its buckytube on a cylinder. The meshing aromatic gear teeth transfer angular momentum from the powered gear to the driven gear. A number of gear and gear/shaft configurations were simulated. Cases in vacuum and with an inert atmosphere were examined. In an extension to molecular dynamics technology, some simulations used a thermostat on the atmosphere while the hydrocarbon gear's temperature was allowed to fluctuate. This models cooling the gears with an atmosphere. Results suggest that these gears can operate at up to 50-100 gigahertz in a vacuum or inert atmosphere at room temperature. The failure mode involves tooth slip, not bond breaking, so failed gears can be returned to operation by lowering temperature and/or rotation rate. Videos and atomic trajectory files in xyz format are presented.

Han, Jie↗

Conducting a thermal conductivity survey

A physically transparent approximate theory of phonon decay rates is presented starting from a pair potential model of the interatomic forces in an insulator or semiconductor. The theory applies in the classical regime and relates the 3-phonon decay rate to the third derivative of the pair potential. Phonon dispersion relations do not need to be calculated, as sum rules relate all the needed quantities directly to the pair potential. The Brillouin zone averaged phonon lifetime turns out to involve a dimensionless measure of the anharmonicity multiplied by an effective density of states for 3-phonon decay. Results are given for rare gas and alkali halide crystals. For rare gases, the results are in good agreement with more elaborate perturbation calculations. Comparison to experimental data on phonon linewidths and thermal conductivity are made.

Allen, P. B.↗

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.↗

Temperature and density dependent pair potential for deuterium under shock

Large-scale classical molecular dynamics (CMD) simulations naturally include the microscopic physics necessary for atomistic modeling of shock release at the ablator-fuel interface in an inertial confinement fusion (ICF) capsule. Here, the multi-megabar shocks utilized in ICF experiments can drive the deuterium fuel from ambient to electron volt temperatures (T) and multi-fold compression. Modeling interatomic interactions over such an extreme range of conditions is challenging for empirical bond order potentials. We generate a pair potential for deuterium with explicit temperature and mass density dependence from ab initio density functional theory molecular dynamics using the iterative Boltzmann inversion method. This potential accurately reproduces the radial distribution functions and pressures from DFT in CMD equilibrium simulations across a wide range of thermodynamic conditions, yet fails to return the expected Hugoniot relations when used in direct CMD shock simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A calculation of the diffusion energies for adatoms on surfaces of F.C.C. metals

The activation energies for diffusion were determined for gold, platinum and iridium adatoms on plane and plane PT surfaces and were found to be in good agreement with the measurements reported by Bassett and Webber. The Lennard-Jones pair potentials were used to model the interatomic forces, and relaxation of the substrate atoms in near proximity to the adatom was considered in detail. The present calculations clarify the mechanism of the observed two-dimensional diffusion of platinum and iridium atoms on a plane PT surface. The results are compared with those obtained using Morse potential functions and different relaxation techniques.

Halicioglu, T.↗

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↗

Interpretation of experimental differential elastic scattering cross section for H/+/ + Ne.

The interatomic V(r) found by potential-model calculations is compared to an ab initio calculation of the intermolecular potential for NeH+ due to Peyerimhoff (1965). The reason for a rather significant difference in the two results is discussed, and a new method due to Remler (1971) involving Regge poles is applied. The method can be used to calculate the differential cross section, where the starting point for such calculations is based on intuition in light of the classical deflection function.

Bobbio, S. M.↗

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.↗

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.↗

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↗

Development and assessment of hierarchical multi-reward reinforcement learning based potential for silicene with state-of-the-art models

We develop a new interatomic force field for Silicene, a 2D material with a buckled hexagonal lattice structure with high polymorphism. We introduce new parameterizations of a Tersoff model using a hierarchical multi-reward reinforcement learning (RL) methodology coupled with a continuous Monte Carlo Tree Search optimization. Our model significantly outperforms existing methods by enhancing the accuracy of predictions for the structural and thermodynamic properties of seven silicene polymorphs-including structure, energy, equation of state, elasticity, and phonon dispersion-when compared to established models. We further make a comprehensive comparison of the various models in predicting the mechanical and thermal properties of silicene. We trace the origin of the improved performance to the description of the angular dependence in the bond-order term, suggesting that modifying the angular terms in short-range models is essential to capture the structural diversity in low dimensional systems.

2D materials↗

Predicting thermophysical properties of molten salts in the MgCl 2 -NaCl-KCl-LiCl system with a shell-model potential

Ternary eutectic salts composed of MgCl 2 , NaCl, and KCl, referred to as MNK salts, have recently emerged as promising candidates as high-temperature heat transfer fluids and thermal energy storage media. Here we performed classical molecular dynamics (MD) simulations to predict the densities, specific heat capacities, viscosities, and ionic self-diffusivities for MNK salts over a wide temperature range. The impact of LiCl additive on their thermophysical properties was also investigated. To capture the electronic polarization of Cl anions by neighboring cations, we developed a novel shell-model potential using the force-matching method and a dataset of ab initio calculated interatomic forces. Our extensive MD simulations predict structure and properties for pure salts and binary/ternary salt mixtures in the MgCl 2 -NaCl-KCl-LiCl system in overall good agreement with available experimental and theoretical data, which corroborates the accuracy and reliability of our developed potential.

36 MATERIALS SCIENCE↗

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.↗

Ground state properties of spin-alined atomic hydrogen

Theoretical calculation of the ground state properties of spin-alined atomic hydrogen by the Monte Carlo method. The interatomic interaction, as described by the results of Kolos and Wolniewicz (1965), is made fit to a Morse potential form. An appropriate trail wavefunction is formed from the short-range part of the WKB solution for a pair of atoms interacting through a Morse potential.

Dugan, J. V., Jr.↗

Ab initio study of the electrostatic multipole nature of torsional potentials in CH3SSCH3, CH3SSH, and HOOH

The origin of torsional potentials in H3CSSCH3, H3CSSH, and HOOH and the anisotropy of the local charge distribution has been analyzed in terms of atomic multipoles calculated from the ab initio LCAO-MO-SCF wave function in the 6-31G* basis set. The results indicate that for longer -S-S-bonds the major contribution to these torsional barriers are electrostatic interactions of the atomic multipoles located on two atoms forming the rotated bond. This finding demonstrates the important role of electrostatic 1-2 interatomic interactions, usually neglected in conformational studies. It also opens the possibility to derive directly from accurate ab initio wave functions a simple nonempirical torsional potential involving atomic multipoles of two bonded atoms defining the torsional angle. For shorter -O-O- bonds, use of more precise models and inclusion of 1-3 interactions seems to be necessary.

NASA Program Exobiology↗