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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Heat flux expressions that satisfy the conservation laws in atomistic system involving multibody potentials

Heat flux expressions are derived for multibody potential systems by extending the original Hardy's methodology and modifying Admal & Tadmor's formulas. The continuum thermomechanical quantities obtained from these two approaches are easy to compute from molecular dynamics (MD) results, and have been tested for a constant heat flux model in two distinctive systems: crystalline iron and polyethylene (PE) polymer. The convergence criteria and affecting parameters, i.e. spatial and temporal window size, and specific forms of localization function are found to be different between the two systems. The conservation of mass, momentum, and energy are discussed and validated within this atomistic–continuum bridging.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Formation of 2D nanoparticles with block structure in simultaneous electric explosion of conductors

A molecular dynamics simulation of nanoparticle formation in simultaneous electric explosion of conductors is performed. Interatomic interaction is described using potentials calculated in the framework of the embedded atom method. High-rate heating results in failure of the conductors with the formation of nanoparticles. The influence of the heating rate, temperature distribution over the specimen cross-section and the distance between simultaneously exploded conductors on the structure of formed nanoparticles is studied. The calculation results show that the electric explosion of conductors allows the formation of nanoparticles with block structure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Hyperelastic tension of graphene

In this paper, we investigate the hyperelastic tensile behaviour of single layer graphene sheets (SLGSs). A one-term incompressible Ogden-type hyperelastic model is chosen to describe the mechanical response of C-C bonds. By establishing equality between the Ogden strain-energy and the variation of the Tersoff-Brenner interatomic potential, three different geometries of SLGSs are studied under tensile loading. We compute the Young's modulus, the finite-deformation Poisson's ratio, ultimate strains, total reactions, and the variation of the potential energy per carbon atom for large strains. Numerical simulations are compared with results obtained by molecular mechanics and molecular dynamics simulations, finite elements, continuum mechanics theory, and experiments. Our predictions are validated, revealing the potential predictive capabilities of the present hyperelastic framework for the analysis of graphene in the context of infinitesimal and large deformations. The good agreement found between our calculations and the published data suggests that graphene may be described as a hyperelastic material.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Growth of bi- and tri-layered graphene on silicon carbide substrate via molecular dynamics simulation

Molecular dynamics (MD) simulation with simulated annealing method is used to study the growth process of bi- and tri-layered graphene on a 6H-SiC (0001) substrate via molecular dynamics simulation. Tersoff-Albe-Erhart (TEA) potential is used to describe the inter-atomic interactions among the atoms in the system. The formation temperature, averaged carbon-carbon bond length, pair correlation function, binding energy and the distance between the graphene formed and the SiC substrate are quantified. The growth mechanism, graphitization of graphene on the SiC substrate and characteristics of the surface morphology of the graphene sheet obtained in our MD simulation compare well to that observed in epitaxially grown graphene experiments and other simulation works.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Structures of 38-atom gold-platinum nanoalloy clusters

Bimetallic nanoclusters, such as gold-platinum nanoclusters, are nanomaterials promising wide range of applications. We perform a numerical study of 38-atom gold-platinum nanoalloy clusters, Au{sub n}Pt{sub 38−n} (0 ≤ n ≤ 38), to elucidate the geometrical structures of these clusters. The lowest-energy structures of these bimetallic nanoclusters at the semi-empirical level are obtained via a global-minimum search algorithm known as parallel tempering multi-canonical basin hopping plus genetic algorithm (PTMBHGA), in which empirical Gupta many-body potential is used to describe the inter-atomic interactions among the constituent atoms. The structures of gold-platinum nanoalloy clusters are predicted to be core-shell segregated nanoclusters. Gold atoms are observed to preferentially occupy the surface of the clusters, while platinum atoms tend to occupy the core due to the slightly smaller atomic radius of platinum as compared to gold’s. The evolution of the geometrical structure of 38-atom Au-Pt clusters displays striking similarity with that of 38-atom Au-Cu nanoalloy clusters as reported in the literature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Thermal transport across graphene and single layer hexagonal boron nitride

As the dimensions of nanocircuits and nanoelectronics shrink, thermal energies are being generated in more confined spaces, making it extremely important and urgent to explore for efficient heat dissipation pathways. In this work, the phonon energy transport across graphene and hexagonal boron-nitride (h-BN) interface is studied using classic molecular dynamics simulations. Effects of temperature, interatomic bond strength, heat flux direction, and functionalization on interfacial thermal transport are investigated. It is found out that by hydrogenating graphene in the hybrid structure, the interfacial thermal resistance (R) between graphene and h-BN can be reduced by 76.3%, indicating an effective approach to manipulate the interfacial thermal transport. Improved in-plane/out-of-plane phonon couplings and broadened phonon channels are observed in the hydrogenated graphene system by analyzing its phonon power spectra. The reported R results monotonically decrease with temperature and interatomic bond strengths. No thermal rectification phenomenon is observed in this interfacial thermal transport. Results reported in this work give the fundamental knowledge on graphene and h-BN thermal transport and provide rational guidelines for next generation thermal interface material designs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Anharmonicity and necessity of phonon eigenvectors in the phonon normal mode analysis

It is well known that phonon frequencies can shift from their harmonic values when elevated to a finite temperature due to the anharmonicity of interatomic potential. Here, we show that phonon eigenvectors also have shifts, but only for compound materials in which each atom has at least two types of anharmonic interactions with other atoms. Using PbTe as the model material, we show that the shifts in some phonon modes may reach as much as 50% at 800 K. Phonon eigenvectors are used in normal mode analysis (NMA) to predict phonon relaxation times and thermal conductivity. We show, from both analytical derivations and numerical simulations, that the eigenvectors are unnecessary in frequency-domain NMA, which gives a critical revision of previous knowledge. This simplification makes the calculation in frequency-domain NMA more convenient since no separate lattice dynamics calculations are needed. On the other hand, we expect our finding of anharmonic eigenvectors may make difference in time-domain NMA and other areas, like wave-packet analysis.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Pressure dependence of electronic properties of BaI{sub 2}

We present Density Functional Theoretical (DFT) calculations of the electronic properties of scintillator material BaI{sub 2} under pressure which were carried out using Perdew-Burke-Ernzerhof genralized gradient approximation. We found that BaI{sub 2} is a direct band gap material with band gap calculated as 3.35 eV. The pressure effect on BaI{sub 2} indicates a linear monotonously decreasing band gap and increasing valence band width with pressure. We have observed the shifting of band extrema from the Γ point with pressure. The pressure coefficient of band gap is found to be −0.047 eV/GPa. The interatomic ionicity factor of BaI{sub 2} is found to be 0.51. Trends in bonding and ionicity under pressure are also discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Computer simulation structure and vibrations of small metal cluster on the Cu (111) surface

Vibrational properties of the small tetrahedral cluster of Co on the Cu (111) surface are studied by using tight-binding second moment approximation interatomic interaction potentials. It was shown that interaction of the clusters with substrate leads to arising of frustrated translation and frustrated rotation in-plane polarized vibrational modes localized on the cluster atoms. The Co{sub 4} cluster on the surface the high frequency modes remain strongly localized and mixed with the nearest neighbor atoms vibrations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Quantum effect on the nucleation of plastic deformation carriers and destruction in crystals

New concepts on the irreversible crystal deformation as a structure transformation caused by a change in interatomic interactions at fluctuations of the electron density under loading are described. The change in interatomic interactions lead to the excitation of dynamical displacements of atoms. A model and a theory of a deformable pristine crystal taking into account the excitation of thermally activated and dynamical displacements of atoms are suggested. New mechanisms of the nucleation of plastic deformation carriers and destruction in pristine crystals at the real value of the deforming stress are studied.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Bond orbital modelling of heterostructures

We propose a non-orthogonal sp{sup 3} hybrid bond orbital model to determine the electronic properties of semiconductor heterostructures. The model considers the non-orthogonality of sp{sup 3} hybrid states of nearest neighboring adjacent atoms using the intra-atomic Coulomb interactions corrected Hartree-Fock atomic energies and metallic contribution to calculate the valence band width energies of group IV elemental and group III-V and II-VI compound semiconductors without any adjustable parameter.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Electronic properties of liquid Hg-In alloys : Ab-initio molecular dynamics study

Ab-initio molecular dynamics simulations are performed to study the structural properties of liquid Hg-In alloys. The interatomic interactions are described by ab-initio pseudopotentials given by Troullier and Martins. Three liquid Hg-In alloys (Hg{sub 10}In{sub 90}, Hg{sub 30}In{sub 70,.} Hg{sub 50}In{sub 50}, Hg{sub 70}In{sub 30}, and Hg{sub 90}Pb{sub 10}) at 299 K are considered. The calculated results for liquid Hg (l-Hg) and lead (l-In) are also drawn. Along with the calculated results of considered five liquid alloys of Hg-In alloy. The results obtained from electronic properties namely total density of state and partial density of states help to find the local arrangement of Hg and In atoms and the presence of liquid state in the considered five alloys.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Dynamic stabilities of icosahedral-like clusters and their ability to form quasicrystals

The dynamic stabilities of the icosahedral-like clusters containing up to 2200 atoms are investigated for 15 metal elements. The clusters originate from five different initial structures (icosahedron, truncated decahedron, octahedron, closed-shell fragment of an HCP structure, and non-closed-shell fragment of an HCP structure). The obtained order of the dynamic stabilities of the icosahedral-like clusters can be assigned to three groups, from stronger to weaker, according to the size ranges involved: (Zr, Al, Ti) > (Cu, Fe, Co, Ni, Mg, Ag) > (Pb, Au, Pd, Pt, Rh, Ir), which correspond to the predicted formation ability of the quasicrystals. The differences of the sequences can be explained by analyzing the parameters of the Gupta-type many-body inter-atomic potentials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Optical control of interatomic interaction in a Bose condensate

We propose to control the interatomic interaction in a Bose condensate by using optical beams forming an atomic trap. It is shown that for a Bose condensate in a double-well optical potential the nonresonant quantized radiation simulates a change in the interatomic interaction in the well. The magnitude and sign of this change can be effectively controlled by varying the frequency of the radiation source forming the trap. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

A ∼32–70 K FORMATION TEMPERATURE RANGE FOR THE ICE GRAINS AGGLOMERATED BY COMET 67 P/CHURYUMOV–GERASIMENKO

Grand Canonical Monte Carlo simulations are used to reproduce the N{sub 2}/CO ratio ranging between 1.7 × 10{sup −3} and 1.6 × 10{sup −2} observed in situ in the Jupiter-family comet 67 P/Churyumov–Gerasimenko (67 P) by the ROSINA mass spectrometer on board the Rosetta spacecraft. By assuming that this body has been agglomerated from clathrates in the protosolar nebula (PSN), simulations are developed using elaborated interatomic potentials for investigating the temperature dependence of the trapping within a multiple-guest clathrate formed from a gas mixture of CO and N{sub 2} in proportions corresponding to those expected for the PSN. By assuming that 67 P agglomerated from clathrates, our calculations suggest the cometary grains must have been formed at temperatures ranging between ∼31.8 and 69.9 K in the PSN to match the N{sub 2}/CO ratio measured by the ROSINA mass spectrometer. The presence of clathrates in Jupiter-family comets could then explain the potential N{sub 2} depletion (factor of up to ∼87 compared to the protosolar value) measured in 67 P.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Heavy-impurity resonance, hybridization, and phonon spectral functions in Fe 1-x M x Si, M=Ir,Os

The vibrational behavior of heavy substitutional impurities (M=Ir,Os) in Fe 1-x M x Si (x = 0, 0.02, 0.04, 0.1) was investigated with a combination of inelastic neutron scattering (INS), transport measurements, and first-principles simulations. In this paper, our INS measurements on single-crystals mapped the four-dimensional dynamical structure factor, S(Q;E), for several compositions and temperatures. Our results show that both Ir and Os impurities lead to the formation of a weakly dispersive resonance vibrational mode, in the energy range of the acoustic phonon dispersions of the FeSi host. We also show that Ir doping, which introduces free carriers and increases electron-phonon coupling, leads to softened interatomic force-constants compared to doping with Os, which is isoelectronic to Fe. We analyze the phonon S(Q,E) from INS through a Green's function model incorporating the phonon self-energy based on first-principles density functional theory (DFT) simulations. Calculations of the quasiparticle spectral functions in the doped system reveal the hybridization between the resonance and the acoustic phonon modes. Finally, our results demonstrate a strong interaction of the host acoustic dispersions with the resonance mode, likely leading to the large observed suppression in lattice thermal conductivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

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↗
DOE OSTI2018

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↗