Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

User Manual - HydraGNN v5.0: Distributed Implementation of Multi-Tasking Graph Neural Networks

This document serves as the user manual for HydraGNN v5.0, a scalable graph neural network (GNN) architecture for simultaneous prediction of multiple target properties using multi-task learning (MTL). This version of HydraGNN has been developed primarily to support the development, training, and deployment of predictive graph-based deep learning (DL) models for atomistic materials modeling. HydraGNN is templated over 13 message-passing policies, including invariant models (GIN, PNA, PNAPlus, GAT, MFC, CGCNN, SAGE, SchNet, DimeNet) and equivariant models (EGNN, PNAEq, PAINN, MACE), and supports distributed training via distributed data parallelism (DDP), DeepSpeed, and Fully Sharded Data Parallelism (FSDP) on leadership-class supercomputers. Although HydraGNN can be applied to problems beyond atomistic materials modeling, its current use is confined to homogeneous graphs. Additional capabilities include machine-learned interatomic potentials with energy-conserving forces, General, Powerful, and Scalable Graph Transformer (GraphGPS) global attention, periodic boundary conditions, hyperparameter optimization, mixed-precision training, and uncertainty quantification.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Structural Principles in Liquids and Glasses: Bottom-Up or Top-Down

The conventional approach to elucidate the atomic structure of liquid and glass is to start with local structural units made of several atoms, and to use them as building blocks to form a global structure, the bottom-up approach. We propose to add an alternative top-down approach in which we start with a global high-temperature gas state and then apply interatomic potentials to all atoms at once. This causes collective density wave instability in all directions with the same wavelength. These two driving forces, local and global, are in competition and are mutually frustrated. The final structure is determined through the compromise of frustration between these two, which creates the medium-range-order. This even-handed approach on global and local potential energy landscapes explains the distinct natures of short-range order and medium-range order, and strong temperature dependence of various properties of liquid.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Atomistic Assessment of Solute-Solute Interactions during Grain Boundary Segregation

Grain boundary solute segregation is becoming increasingly common as a means of stabilizing nanocrystalline alloys. Thermodynamic models for grain boundary segregation have recently revealed the need for spectral information, i.e., the full distribution of environments available at the grain boundary during segregation, in order to capture the essential physics of the problem for complex systems like nanocrystalline materials. However, there has been only one proposed method of extending spectral segregation models beyond the dilute limit, and it is based on simple, fitted parameters that are not atomistically informed. In this work, we present a physically motived atomistic method to measure the full distribution of solute-solute interaction energies at the grain boundaries in a polycrystalline environment. We then cast the results into a simple thermodynamic model, analyze the Al(Mg) system as a case study, and demonstrate strong agreement with physically rigorous hybrid Monte Carlo/molecular statics simulations. This approach provides a means of rapidly measuring key interactions for non-dilute grain boundary segregation for any system with an interatomic potential.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Data and Code for Atomic Scale Etching of Diamond: Insights from Molecular Dynamics Simulations

This work investigates the effects of argon ions, hydrogen atoms, and hydrogen ions on the diamond (100) surface using classical molecular dynamics simulations. The purpose of this investigation was to asses plasma processing techniques for applications in quantum device manufacturing. The simulations suggest that combining argon ion smoothing with selective, near threshold energy H removal of amorphous C could be an effective strategy for diamond surface engineering, leading to more reliable and sensitive diamond color center devices. Results were found to differ significantly with interatomic potential, and an analysis of these differences was also carried out. Included in this repository are LAMMPS source files, input scripts, and plotting scripts required to reproduce the data. Also included are the output data required to make all the plots included in the associated publication.

Brenner↗
NASA NTRS1975

Lattice model theory of the equation of state covering the gas, liquid, and solid phases

The three stable states of matter and the corresponding phase transitions were obtained with a single model. Patterned after Lennard-Jones and Devonshires's theory, a simple cubic lattice model containing two fcc sublattices (alpha and beta) is adopted. The interatomic potential is taken to be the Lennard-Jones (6-12) potential. Employing the cluster variation method, the Weiss and the pair approximations on the lattice gas failed to give the correct phase diagrams. Hybrid approximations were devised to describe the lattice term in the free energy. A lattice vibration term corresponding to a free volume correction is included semi-phenomenologically. The combinations of the lattice part and the free volume part yield the three states and the proper phase diagrams. To determine the coexistence regions, the equalities of the pressure and Gibbs free energy per molecule of the coexisting phases were utilized. The ordered branch of the free energy gives rise to the solid phase while the disordered branch yields the gas and liquid phases. It is observed that the triple point and the critical point quantities, the phase diagrams and the coexistence regions plotted are in good agreement with the experimental values and graphs for argon.

Bonavito, N. L.↗
NASA NTRS1976

Computer simulation of screw dislocation in aluminum

The atomic structure in a 110 screw dislocation core for aluminum is obtained by computer simulation. The lattice statics technique is employed since it entails no artificially imposed elastic boundary around the defect. The interatomic potential has no adjustable parameters and was derived from pseudopotential theory. The resulting atomic displacements were allowed to relax in all three dimensions.

Esterling, D. M.↗
NASA NTRS1978

Solid-liquid phase transition in argon

Starting from the Lennard-Jones interatomic potential, a modified cell theory has been used to describe the solid-liquid phase transition in argon. The cell-size variations may be evaluated by a self-consistent condition. With the inclusion of cell-size variations, the transition temperature, the solid and liquid densities, and the liquid-phase radial-distribution functions have been calculated. These ab initio results are in satisfactory agreement with molecular-dynamics calculations as well as experimental data on argon.

Tsang, T.↗
NASA NTRS1980

Dislocation dissociation in some f.c.c. metals

The dissociation of a perfect screw dislocation into a stacking fault in an f.c.c. lattice is modeled by the modified lattice statics. The interatomic potentials are obtained from the work of Esterling and Swaroop and differ substantially from those empirical potentials usually employed in defect simulations. The calculated stacking fault widths for aluminum, copper, and silver are in good agreement with weak beam microscopy results.

Esterling, D. M.↗
NASA NTRS1980

Time correlation functions in classical dense fluids

The velocity autocorrelation functions and memory functions of dense classical fluids may be directly obtained from the static radial distribution function g(r) in an approximate way. Following the Mori projection operator formalism, the memory functions may be related to the fluctuating force correlation. At low densities, these functions may be evaluated by following the trajectories of particle pairs in the interatomic potential. At higher densities, the force correlation functions can be evaluated approximately from particle pair trajectories via the potential of the mean force. The results for argon fluids at various densities and temperatures agree satisfactorily with the molecular dynamics and the Enskog values. The decrease of the diffusion coefficient with density is partly due to the nature of g(r) which reflects the stronger clustering of atoms at higher densities.

Tsang, T.↗
NASA NTRS1981

Far wing depolarization of light - Generalized absorption profiles

An absorption (and/or emission) event which takes place during a strong collision is called a 'correlated event'. It is discussed how correlated events affect the far red wing depolarization of fluorescence. Attention is given to an atomic vapor which is irradiated by linearly polarized light of a frequency on the red side of the resonance line. Two limiting cases are considered, corresponding to excitation in the impact region and in the quasi-static wing. In the quasi-static wing, absorption of a photon followed by fluorescence (rather than Rayleigh scattering), occurs mostly during a collision. Correlated events dominate the scattering process. Expressions derived for the polarization of the fluorescent light are applied to far red wing depolarization. It is found that the polarization of the fluorescent light does not go to zero in the far wing, but depends crucially on the detailed nature of the anisotropy in the long-range part of the interatomic potential.

Thomann, P.↗
NASA NTRS1982

A study pertaining to a very low temperature hydrogen maser feasibility

Very low temperature atomic hydrogen masers were developed. The advantages over room temperature hydrogen masers include higher radiated power due to the higher beam intensities possible using low temperature techniques and the much smaller electron spin flip cross section at low temperatures, lower cavity and amplifier noise temperatures, increased stability against mechanical creep at low temperatures, and opportunities to extend H atom collision studies to low temperatures, where quantum effects and details of the interatomic potentials are much more important.

Crampton, S. B.↗
NASA NTRS1982

Anharmonic properties of solids from measurements of the stress acoustic constant

The equations of elastic motion and their solutions are generalized in order to include nonzero homogeneous initial stresses and redefine the stress acoustic constants to include the effect of initial stress. In deriving the relationship between the stress acoustic constants and the strain-generalized Gruneisen parameters, implications in material anharmonicity and nonlinear thermoelasticity are discovered. It is found that the linear change in velocity obtained from these measurements is a nonlinear process resulting from anharmonicity in the interatomic potential. In addition, the thermal expansion coefficient can also be expressed in terms of linear combinations.

Cantrell, J. H., Jr.↗
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 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 NTRS1991

Nonequilibrium gas flows. I - A detailed validation of Monte Carlo direct simulation for monatomic gases

One-dimensional shock wave properties in helium and argon are predicted using Monte Carlo direct simulation. The collision model is based directly on the interatomic potential taking angular scattering into account. The potential is assumed to be of the Maitland-Smith n(r)-6 form. The detailed validity of the simulation is studied by comparing the predicted macroscopic and miroscopic flow properties in shock waves to a wide range of available data.

Erwin, Daniel A.↗
NASA NTRS1993

A laser-cooled cesium fountain frequency standard and a measurement of the frequency shift due to ultra-cold collisions

A frequency standard based on an atomic fountain of cesium atoms may have an accuracy of 10(exp -16) due to longer interaction times and smaller anticipated systematic errors. All of the known systematic effects that now limit the accuracy of the Cs frequency standard increase either linearly or as some higher power of the atom's velocity. The one systematic frequency shift which is dramatically different is the frequency shift due to the collisions between the laser cooled atoms. At a temperature of a few micro-K, the de Broglie wavelength (lambda(sub deB) = h/p, where h is Planck's constant and p is the momentum of the atom) is much larger than the scale of the interatomic potential. Under these conditions the collision cross sections can be as large as (lambda(sub deB)(sup 2))/Pi and the frequency shift due to these collisions was recently calculated. In our Cs atomic fountain, we laser cooled and trapped 10(exp 10) Cs atoms in 0.4 s. By shifting the frequencies of the laser beams, the atoms were launched upwards at 2.5 m/s and a fraction of the atoms were optically pumped into the F=3 ground state. The unwanted atoms in the F=4 ground state were removed from the fountain with radiation pressure from a laser beam tuned to excite only those atoms. The Cs atoms in the F=3 state traveled ballistically upwards, were excited by the microwave cavity, and then returned back through the same cavity in the atomic fountain configuration. By varying the cold atom density, a density dependent shift of -12.9 +/- 0.7 mHz or -1.4 x 10-12 for an average fountain density of (2.7 +/- 1.5) 10(exp 9) atoms/cm(sup 3) was measured.

Gibble, Kurt↗
NASA NTRS1993

Evaporation kinetics of Mg2SiO4 crystals and melts from molecular dynamics simulations

Computer simulations based on the molecular dynamics (MD) technique were used to study the mechanisms and kinetics of free evaporation from crystalline and molten forsterite (i.e., Mg2SiO4) on an atomic level. The interatomic potential employed for these simulations reproduces the energetics of bonding in forsterite and in gas-phase MgO and SiO2 reasonably accurately. Results of the simulation include predicted evaporation rates, diffusion rates, and reaction mechanisms for Mg2SiO4(s or l) yields 2Mg(g) + 20(g) + SiO2(g).

Kubicki, J. D.↗