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 91 records · Page 5

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↗

Automating Anharmonic Lattice Thermal Conductivity Calculations using Compressed Sensing Lattice Dynamics [Poster]

Engineering heat transfer is critical for applications in heat exchangers, semiconductor devices, thermoelectrics and more. This demand motivates a high-throughput computational methodology for systematically designing materials with improved thermal properties. One emerging method which can rapidly and explicitly consider phonon-phonon interactions even for highly anharmonic materials is Compressive Sensing Lattice Dynamics (CSLD). In fact, CSLD accurately predicts 2nd (harmonic) and 3rd+ order (anharmonic) interatomic force constants (IFCs) with orders of magnitude fewer Density Functional Theory (DFT) calculations than conventional methods. However, doing CSLD can be an exhaustive, many-step process requiring an intimate knowledge of its sensitivity to parameters. Consequently, this work implements an automatic CSLD workflow capable of obtaining the thermal conductivity of potentially thousands of materials, benchmarks the stability of the workflow against materials with a range of anharmonicity, and begins constructing a dataset of thermal conductivity values and phonon dispersion curves to be stored in the Materials Project for public use.

Chang, Rees↗

Phonon-Spin Scattering from Unpaired f-electrons in U atoms and It’s Influence on Thermal Transport in Uranium-doped Thorium Dioxide Single Crystals

In this work, impact of low level of uranium (U) atom substitution on thermal conductivity of thorium dioxide (ThO2) is investigated. ThO2 is an electronic insulator with a wide optical band-gap and no unpaired electrons whose thermal transport is governed by phonons. U-substitution introduces unpaired f-electrons resulting in paramagnetic behavior of U-ThO2 at room temperature, which significantly suppresses its thermal conductivity. A single crystal of U-ThO2 with graded composition of U is grown using a hydrothermal synthesis method, and thermal conductivity measurements are performed in regions with uniform composition of U at levels of 0%, 6%, 9% and 16%. Measured thermal conductivity profiles over 77–300 K temperature range are analyzed using an analytical expression for phonon-mediated thermal transport based on Klemens-Callaway model. Temperature dependent thermal conductivity is found to deviate significantly from the Rayleigh scattering trend expected for a simple substitutional point defect with a small perturbation to mass and interatomic forces. With the resonant scattering term, observed large suppression of thermal conductivity at low temperatures can be closely reproduced. Additionally, the extracted phonon-spin coupling constants imply a nonlinear relation of phonon-spin interaction intensity with respect to U doping percentage. Our study reveals how phonon-spin scattering contributed by unpaired f-electrons in U atoms influences thermal transport in the U-ThO2 system.

36 - MATERIALS SCIENCE↗

Multi-fidelity learning for interatomic potentials: low-level forces and high-level energies are all you need

The promise of machine learning interatomic potentials (MLIPs) has led to an abundance of public quantum mechanical (QM) training datasets. The quality of an MLIP is directly limited by the accuracy of the energies and atomic forces in the training dataset. Unfortunately, most of these datasets are computed with relatively low-accuracy QM methods, e.g. density functional theory with a moderate basis set. Due to the increased computational cost of more accurate QM methods, e.g. coupled-cluster theory with a complete basis set (CBS) extrapolation, most high-accuracy datasets are much smaller and often do not contain atomic forces. The lack of high-accuracy atomic forces is quite troubling, as training with force data greatly improves the stability and quality of the MLIP compared to training to energy alone. Because most datasets are computed with a unique level of theory, traditional single-fidelity (SF) learning is not capable of leveraging the vast amounts of published QM data. In this study, we apply multi-fidelity learning (MFL) to train an MLIP to multiple QM datasets of different levels of accuracy, i.e. levels of fidelity. Specifically, we perform three test cases to demonstrate that MFL with both low-level forces and high-level energies yields an extremely accurate MLIP—far more accurate than a SF MLIP trained solely to high-level energies and almost as accurate as a SF MLIP trained directly to high-level energies and forces. Therefore, MFL greatly alleviates the need for generating large and expensive datasets containing high-accuracy atomic forces and allows for more effective training to existing high-accuracy energy-only datasets. Indeed, low-accuracy atomic forces and high-accuracy energies are all that are needed to achieve a high-accuracy MLIP with MFL.

36 MATERIALS SCIENCE↗

Simulation of Simple and Complex Gadolinium Molybdates by the Interatomic Potential Method

Crystals of ferroelectric‒ferroelastic gadolinium molybdate Gd{sub 2}(MoO{sub 4}){sub 3}, calcium molybdate CaMoO{sub 4}, and double sodium‒gadolinium molybdates of stoichiometric (Na{sub 1/2}Gd{sub 1/2}MoO{sub 4}) and cationdeficient (Na{sub 2/7}Gd{sub 4/7}MoO{sub 4}) compositions, which are used to design solid-state lasers, phosphors, and white LEDs, have been simulated by the interatomic potential method. Their structural, mechanical, and thermodynamic properties are calculated using a unified system of interatomic potentials and effective ion charges, which demonstrated transferability and made it possible not only to describe the existing experimental data but also to predict some important physical and thermodynamic properties of molybdate crystals. The influence of the deviation from stoichiometry and partial ordering of cations on sites in nonstoichiometric crystals on the properties and local structure of sodium‒gadolinium molybdates is discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Anisotropy of Optical Activity in Langasite Family Crystals

The rotation of the plane of polarization of light transmitted in the directions parallel (ρ{sub ||}) and perpendicular (ρ{sub ⊥}) to the optical axis has been calculated for langasite family crystals. It is shown that the ρ{sub ||} and ρ{sub ⊥} values are, respectively, negative and positive for all the right-handed crystals under study. An analysis of the interatomic potential constructed based on X-ray diffraction data has suggested that this anisotropy of optical activity is structurally based of electron density regions shaped as helices. The atomic helix with an axis parallel to the optical axis is right-handed and three-start, whereas in the direction perpendicular to the optical axis there are three left-handed two-start helices, which intersect in the origin of the unit-cell coordinate system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

Diffusion bonding and its application to manufacturing

In its simplest form diffusion bonding is accomplished by placing clean metal surfaces together under a sufficient load and heating. The natural interatomic attractive force between atoms transforms the interface into a natural grain boundary. Therefore, in principle, the properties of the bond area are identical to those of the parent metal. Other advantages of diffusion bonding over conventional methods of bonding include freedom from residual stresses, excessive deformation, foreign metals, or changed crystal structures. Stainless steels, nickel-base superalloys, and aluminum alloys have all been successfully joined. Complex hardware, including integrated flueric devices, jet engine servovalves, and porous woven structures have been fabricated. The processing involved is discussed, along with such theoretical considerations as the role of metal surfaces, the formation of metal contact junctions, and the mechanisms of material transport in diffusion bonding.

Spurgeon, W. M.↗

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↗