Results for “Interatomic potentials” — Page 13

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

NASA NTRS2000

The Rotational Spectra of IO X(sub 1) (sup 2)pi(sub 3/2), v <= 13 and X(sub 2) (sup 2)pi(sub 1/2), v <= 9

The rotational spectra of IO in vibrational states up to v = 13 in the X(sub 1) (sup 2)pi(sub 3/2) state and up to v = 9 in the X2 (sup 2)pi(sub 1/2) state have been observed in an O2 discharge over molecular I2. In addition, I(18)O has been observed for both the X(sub 1) and X(sub 2) states up to v = 5. All data have been analyzed simultaneously with fixed isotopic ratios among the constants. This extends the data set for the X(sub 1) state described last year at this meeting and provides the first high resolution data for the X(sub 2) state and for I(18)O. An extensive set of parameters has been derived. These will be interpreted in terms of the electronic structure and the interatomic potential.

Miller, Charles E.↗

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 NTRS2011

Lattice Thermal Conductivity of Ultra High Temperature Ceramics ZrB2 and HfB2 from Atomistic Simulations

Atomistic Green-Kubo simulations are performed to evaluate the lattice thermal conductivity for single crystals of the ultra high temperature ceramics ZrB2 and HfB2 for a range of temperatures. Recently developed interatomic potentials are used for these simulations. Heat current correlation functions show rapid oscillations which can be identified with mixed metal-Boron optical phonon modes. Agreement with available experimental data is good.

Lawson, John W.↗

NASA NTRS2011

Lattice Thermal Conductivity of Ultra High Temperature Ceramics (UHTC) ZrB2 and HfB2 from Atomistic Simulations

Ultra high temperature ceramics (UHTC) including ZrB2 and HfB2 are candidate materials for applications in extreme environments because of their high melting point, good mechanical properties and reasonable oxidation resistance. Unlike many ceramics, these materials have high thermal conductivity which can be advantageous, for example, to reduce thermal shock. Recently, we developed Tersoff style interatomic potentials for both ZrB2 and HfB2 appropriate for atomistic simulations. As an application, Green-Kubo molecular dynamics simulations were performed to evaluate the lattice thermal conductivity for single crystals of ZrB2 and HfB2. The atomic mass difference in these binary compounds leads to oscillations in the time correlation function of the heat current. Results at room temperature and at elevated temperatures will be reported.

Lawson, JOhn W.↗

NASA NTRS2012

Multiscale Modeling of Ultra High Temperature Ceramics (UHTC) ZrB2 and HfB2: Application to Lattice Thermal Conductivity

We are developing a multiscale framework in computational modeling for the ultra high temperature ceramics (UHTC) ZrB2 and HfB2. These materials are characterized by high melting point, good strength, and reasonable oxidation resistance. They are candidate materials for a number of applications in extreme environments including sharp leading edges of hypersonic aircraft. In particular, we used a combination of ab initio methods, atomistic simulations and continuum computations to obtain insights into fundamental properties of these materials. Ab initio methods were used to compute basic structural, mechanical and thermal properties. From these results, a database was constructed to fit a Tersoff style interatomic potential suitable for atomistic simulations. These potentials were used to evaluate the lattice thermal conductivity of single crystals and the thermal resistance of simple grain boundaries. Finite element method (FEM) computations using atomistic results as inputs were performed with meshes constructed on SEM images thereby modeling the realistic microstructure. These continuum computations showed the reduction in thermal conductivity due to the grain boundary network.

Lawson, John W.↗

NASA NTRS2012

Multiscale Modeling of UHTC: Thermal Conductivity

We are developing a multiscale framework in computational modeling for the ultra high temperature ceramics (UHTC) ZrB2 and HfB2. These materials are characterized by high melting point, good strength, and reasonable oxidation resistance. They are candidate materials for a number of applications in extreme environments including sharp leading edges of hypersonic aircraft. In particular, we used a combination of ab initio methods, atomistic simulations and continuum computations to obtain insights into fundamental properties of these materials. Ab initio methods were used to compute basic structural, mechanical and thermal properties. From these results, a database was constructed to fit a Tersoff style interatomic potential suitable for atomistic simulations. These potentials were used to evaluate the lattice thermal conductivity of single crystals and the thermal resistance of simple grain boundaries. Finite element method (FEM) computations using atomistic results as inputs were performed with meshes constructed on SEM images thereby modeling the realistic microstructure. These continuum computations showed the reduction in thermal conductivity due to the grain boundary network.

Lawson, John W.↗

NASA NTRS2012

Lattice Thermal Conductivity from Atomistic Simulations: ZrB2 and HfB2

Ultra high temperature ceramics (UHTC) including ZrB2 and HfB2 have a number of properties that make them attractive for applications in extreme environments. One such property is their high thermal conductivity. Computational modeling of these materials will facilitate understanding of fundamental mechanisms, elucidate structure-property relationships, and ultimately accelerate the materials design cycle. Progress in computational modeling of UHTCs however has been limited in part due to the absence of suitable interatomic potentials. Recently, we developed Tersoff style parameterizations of such potentials for both ZrB2 and HfB2 appropriate for atomistic simulations. As an application, Green-Kubo molecular dynamics simulations were performed to evaluate the lattice thermal conductivity for single crystals of ZrB2 and HfB2. The atomic mass difference in these binary compounds leads to oscillations in the time correlation function of the heat current, in contrast to the more typical monotonic decay seen in monoatomic materials such as Silicon, for example. Results at room temperature and at elevated temperatures will be reported.

Lawson, John W.↗

NASA NTRS2012

Lattice Thermal Conductivity of Ultra High Temperature Ceramics (UHTC) ZrB2 and HfB2 from Atomistic Simulations

Ultra high temperature ceramics (UHTC) including ZrB2 and HfB2 have a number of properties that make them attractive for applications in extreme environments. One such property is their high thermal conductivity. Computational modeling of these materials will facilitate understanding of fundamental mechanisms, elucidate structure-property relationships, and ultimately accelerate the materials design cycle. Progress in computational modeling of UHTCs however has been limited in part due to the absence of suitable interatomic potentials. Recently, we developed Tersoff style parameterizations of such potentials for both ZrB2 and HfB2 appropriate for atomistic simulations. As an application, Green-Kubo molecular dynamics simulations were performed to evaluate the lattice thermal conductivity for single crystals of ZrB2 and HfB2. The atomic mass difference in these binary compounds leads to oscillations in the time correlation function of the heat current, in contrast to the more typical monotonic decay seen in monoatomic materials such as Silicon, for example. Results at room temperature and at elevated temperatures will be reported.

Lawson, John W.↗

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↗