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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

First principles Peierls-Boltzmann phonon thermal transport: A topical review

The advent of coupled thermal transport calculations with interatomic forces derived from density functional theory has ushered in a new era of fundamental microscopic insight into lattice thermal conductivity. Subsequently, significant new understanding of phonon transport behavior has been developed with these methods, and because they are parameter free and successfully benchmarked against a variety of systems, they also provide reliable predictions of thermal transport in systems for which little is known. This topical review will describe the foundation from which first principles Peierls-Boltzmann transport equation methods have been developed, and briefly describe important necessary ingredients for accurate calculations. Sample highlights of reported work will be presented to illustrate the capabilities and challenges of these techniques, and to demonstrate the suite of tools available, with an emphasis on thermal transport in micro- and nano-scale systems. In conclusion, future challenges and opportunities will be discussed, drawing attention to prospects for methods development and applications.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

Fast and accurate calculation of EXAFS Debye-Waller factors in U⁢O2 using the dynamical matrix method

Theoretical modeling of bonding dynamics in metal oxides is required for predicting their thermal conductivity, catalytic activity, and mechanical properties. A primary challenge is the scarcity of experimental methods for validating theoretical predictions of these atomic-scale dynamics. This work presents a workflow that uses experimental extended x-ray absorption fine structure (EXAFS) data collected at high temperatures to validate an interatomic force field for uranium dioxide (UO2), an important model material. The validated force field is then used to drive computationally intensive molecular dynamics (MD) simulations and as input for the much faster dynamical matrix Debye-Waller (DMDW) method. The predicted values of the Debye-Waller factors from the DMDW calculations are in good agreement with those obtained from the MD simulations, with residual pair-specific differences attributable to quantum zero-point motion at low temperatures and lattice anharmonicity at high temperatures. We further show that theoretical EXAFS spectra constructed directly from DMDW-derived Debye-Waller factors reproduce the experimental data (at relatively low temperatures) with accuracy comparable to full MD-EXAFS, providing an additional validation of the choice of the potential. This study establishes a validated, rapid computational pathway for modeling bond dynamics, naturally incorporating quantum nuclear\\\\r\\\\nstatistics absent in classical simulations, which are essential for the mechanistic understanding of complex oxide materials.

58 GEOSCIENCES↗
DOE OSTI2019

Crystalline phase transitions and vibrational spectra of silicon up to multiterapascal pressures

A composite high-pressure phase diagram for silicon has been predicted for up to 4 TPa. Here, this diagram has been built using a combination of evolutionary algorithm-based structure searches, electronic density functional theory, lattice dynamics of perfect crystals, and anharmonic corrections for the solid state in conjunction with molecular dynamics for evaluating the melt curve. The anharmonic corrections to free energy, arising from both finite-temperature three-phonon interactions and temperature dependence of the lattice axial ratios, play a critical role in properly identifying the solid–solid transition boundaries of the orthorhombic structures. A double hexagonal close-packed structure has been found to be thermodynamically and dynamically stable, sandwiched between the experimentally observed base-centered orthorhombic and hexagonal close-packed structures. In addition, beyond 2.8 TPa, there exists a sequence of face-centered cubic-to-body-centered cubic-to-simple cubic transitions that are accompanied by the localization of electrons in the interstitial spaces between the ions. Supplementing the structural calculations, second- and third-order interatomic force constants were evaluated to compute the phonon vibration modes and linewidths, respectively. This allowed an elaborate analysis of the Raman and infrared spectra for all of the structures of silicon identified so far.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Dominant effect of anharmonicity on the equation of state and thermal conductivity of MgO under extreme conditions

Understanding heat transfer in the earth mantle region is of great scientific interest but highly challenging due to dissimilar conditions from our daily lives. Thermal conductivity of highly pressured and thermally activated materials can be completely different from that at ambient condition. In this study we calculate the equation of state and thermal conductivity of cubic crystalline MgO, a major component of the mantle, in such an environment. It is shown that the material properties are not accurately captured unless rigorous treatment of anharmonicity of phonon vibrational motion is explicitly applied. To predict renormalized phonon dispersion relations, interatomic force constants of high-order expansion are required, especially at elevated temperature. The anharmonicity plays an increasingly substantial role for regulating volume expansion as temperature rises, because the Gibbs free energy of MgO considerably deviates from what was predicted by the quasiharmonic approximation. We find that the lattice thermal conductivity of MgO and the heat transfer of the earth's mantle are reliably estimated only when both the temperature-dependent renormalization of phonon frequencies and four-phonon scattering effects are taken into consideration. We map the characteristics at Geotherm into contour diagrams in which equilibrium volumes and the thermal conductivities are described as functions of temperature and pressure. Furthermore, the thermal conductivity values of MgO in this study are smaller than those in previous theoretical studies due to the additional incorporation of four-phonon scatterings in the calculations. Specifically, deviation between the values in previous studies and in this study becomes larger as thermodynamic variables reach those near the core-mantle boundary. For instance, thermal conductivities of 52.15 and 23.67 Wm -1 K -1 are estimated at ambient and extreme conditions (T = 4000K and P = 135 GPa), respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Ab initio electronic density in solids by many-body plane-wave auxiliary-field quantum Monte Carlo calculations

We present accurate many-body results of the electronic densities in several solid materials, including Si, NaCl, and Cu. These results are obtained using the ab initio auxiliary-field quantum Monte Carlo (AFQMC) method working in a plane-wave basis with norm-conserving, multiple-projector pseudopotentials. AFQMC has been shown to be an excellent many-body total energy method. Computation of observables and correlation functions other than the ground-state energy requires back-propagation, whose adaption and implementation in the plane-wave basis AFQMC framework are discussed in the present paper. This development allows us to compute correlation functions, electronic densities, and interatomic forces, paving the way for geometry optimizations and calculations of thermodynamic properties in solids. Finite supercell size effects are considerably more subtle in the many-body framework than in independent-electron calculations. We analyze the convergence of the electronic density, and obtain best estimates for the thermodynamic limit. The densities from several typical density functionals are benchmarked against our near-exact results. The electronic densities we obtained can also be used to help construct improved density functionals.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Hybrid Kohn-Sham+Thomas-Fermi scheme for high-temperature density functional theory

We report density functional theory has revolutionised our understanding of dense, strongly coupled systems in the condensed-matter limit. At elevated temperatures, such as those found in warm dense matter and dense plasma, the computational cost of calculating enough states to accurately describe a thermal ensemble of electrons is a major obstacle to practical Kohn-Sham calculations. This problem has recently been tackled from a number of approaches, including a plane-wave approximation for electrons at high energy. In this work, we demonstrate how the Thomas-Fermi density of states, equivalent to a local plane-wave approximation, can yield improved results for thermodynamic quantities, electron density profiles, and interatomic forces while requiring fewer states than existing methods.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Low-energy moiré phonons in twisted bilayer van der Waals heterostructures

We develop a low-energy continuum model for phonons in twisted moiré bilayers, based on a configuration-space approach. In this approach, interatomic force constants are obtained from density functional theory (DFT) calculations of untwisted bilayers with various in-plane shifts. This allows for efficient computation of phonon properties for any small twist angle, while maintaining DFT-level accuracy. Based on this framework, we show how the low-energy phonon modes, including interlayer shearing and layer-breathing modes, vary with the twist angle. As the twist angle decreases, the frequencies of the low-energy modes are reordered and the atomic displacement fields corresponding to phonon eigenmodes break translational symmetry, developing periodicity on the moiré length scale. Here, we demonstrate the capabilities of our model by calculating the phonon properties of three specific structures: bilayer graphene, bilayer molybdenum disulfide (MoS 2 ), and molybdenum diselenide-tungsten diselenide (MoSe 2 -WSe 2 ).

2-dimensional systems↗
DOE OSTI2015

Reflection thermal diffuse x-ray scattering for quantitative determination of phonon dispersion relations

Synchrotron reflection x-ray thermal diffuse scattering (TDS) measurements, rather than previously reported transmission TDS, are carried out at room temperature and analyzed using a formalism based upon second-order interatomic force constants and long-range Coulomb interactions to obtain quantitative determinations of MgO phonon dispersion relations (h) over bar omega(j) (q), phonon densities of states g((h) over bar omega), and isochoric temperature-dependent vibrational heat capacities c v (T). We use MgO as a model system for investigating reflection TDS due to its harmonic behavior as well as its mechanical and dynamic stability. Resulting phonon dispersion relations and densities of states are found to be in good agreement with independent reports from inelastic neutron and x-ray scattering experiments. Temperature-dependent isochoric heat capacities c v (T), computed within the harmonic approximation from (h) over bar omega(j) (q) values, increase with temperature from 0.4 x 10 -4 eV/atom K at 100 K to 1.4 x 10 -4 eV/atom K at 200 K and 1.9 x 10 -4 eV/atom K at 300 K, in excellent agreement with isobaric heat capacity values c p (T) between 4 and 300 K. We anticipate that the experimental approach developed here will be valuable for determining vibrational properties of heteroepitaxial thin films since the use of grazing-incidence (θ ≲ θ c where θ c is the density-dependent critical angle) allows selective tuning of x-ray penetration depths to ≲ 10 nm.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Electron-phonon coupling and thermal transport in the thermoelectric compound Mo 3 Sb 7–x Te x

Phonon properties of Mo 3 Sb 7–x Te x (x = 0, 1.5, 1.7), a potential high-temperature thermoelectric material, have been studied with inelastic neutron and x-ray scattering, and with first-principles simulations. The substitution of Te for Sb leads to pronounced changes in the electronic struc- ture, local bonding, phonon density of states (DOS), dispersions, and phonon lifetimes. Alloying with tellurium shifts the Fermi level upward, near the top of the valence band, resulting in a strong suppression of electron-phonon screening, and a large overall stiffening of interatomic force- constants. The suppression in electron-phonon coupling concomitantly increases group velocities and suppresses phonon scattering rates, surpassing the effects of alloy-disorder scattering, and re- sulting in a surprising increased lattice thermal conductivity in the alloy. We also identify that the local bonding environment changes non-uniformly around different atoms, leading to variable perturbation strengths for different optical phonon branches. The respective roles of changes in phonon group velocities and phonon lifetimes on the lattice thermal conductivity are quantified. Lastly, our results highlight the importance of the electron-phonon coupling on phonon mean-free-paths in this compound, and also estimates the contributions from boundary scattering, umklapp scattering, and point-defect scattering.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Electronic band structure effects in the stopping of protons in copper [Electronic band structure non-linear effects in the stopping of protons in copper]

Here, we present an ab initio study of the electronic stopping power of protons in copper over a wide range of proton velocities v = 0.02–10a.u. where we take into account nonlinear effects. Time-dependent density functional theory coupled with molecular dynamics is used to study electronic excitations produced by energetic protons. A plane-wave pseudopotential scheme is employed to solve the time-dependent Kohn-Sham equations for a moving ion in a periodic crystal. The electronic excitations and the band structure determine the stopping power of the material and alter the interatomic forces for both channeling and off-channeling trajectories. Our off-channeling results are in quantitative agreement with experiments, and at low velocity they unveil a crossover region of superlinear velocity dependence (with a power of ~1.5) in the velocity range v = 0.07–0.3a.u., which we associate to the copper crystalline electronic band structure. The results are rationalized by simple band models connecting two separate regimes. We find that the limit of electronic stopping v → 0 is not as simple as phenomenological models suggest and it is plagued by band-structure effects.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Effects of functional group mass variance on vibrational properties and thermal transport in graphene

Intrinsic thermal resistivity critically depends on features of phonon dispersions dictated by harmonic interatomic forces and masses. We present the effects of functional group mass variance on vibrational properties and thermal conductivity (κ ) of functionalized graphene from first principles calculations. We also use graphane, a buckled graphene backbone with covalently bonded Hydrogen atoms on both sides, as the base material and vary the mass of the Hydrogen atoms to simulate the effect of mass variance from other functional groups. We find non-monotonic behavior of κ with increasing mass of the functional group and an unusual cross-over from acoustic-dominated to optic-dominated thermal transport behavior. We connect this cross-over to changes in the phonon dispersion with varying mass which suppress acoustic phonon velocities, but also give unusually high velocity optic modes. Further, we show that out-of-plane acoustic vibrations contribute significantly more to thermal transport than in-plane acoustic modes despite breaking of a reflection symmetry based scattering selection rule responsible for their large contributions in graphene. Our work demonstrates the potential for manipulation and engineering of thermal transport properties in two dimensional materials toward targeted applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Ab initio phonon point defect scattering and thermal transport in graphene

Here, we study the scattering of phonons from point defects and their effect on lattice thermal conductivity κ using a parameter-free ab initio Green's function methodology. Specifically, we focus on the scattering of phonons by boron (B), nitrogen (N), and phosphorus substitutions as well as single- and double-carbon vacancies in graphene. We show that changes of the atomic structure and harmonic interatomic force constants locally near defects govern the strength and frequency trends of the scattering of out-of-plane acoustic (ZA) phonons, the dominant heat carriers in graphene. ZA scattering rates due to N substitutions are nearly an order of magnitude smaller than those for B defects despite having similar mass perturbations. Furthermore, ZA phonon scattering rates from N defects decrease with increasing frequency in the lower-frequency spectrum in stark contrast to expected trends from simple models. ZA phonon-vacancy scattering rates are found to have a significantly softer frequency dependence (~ω 0 ) in graphene than typically employed in phenomenological models. The rigorous Green's function calculations demonstrate that typical mass-defect models do not adequately describe ZA phonon-defect scattering rates. Our ab initio calculations capture well the trend of κ vs vacancy density from experiments, though not the magnitudes. In conclusion, this work elucidates important insights into phonon-defect scattering and thermal transport in graphene, and demonstrates the applicability of first-principles methods toward describing these properties in imperfect materials.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Thermal conductivity of InN with point defects from first principles

We present ab initio calculations of thermal conductivity of InN with vacancies and substitutional defects using a full solution of the Peierls-Boltzmann transport equation. Our parameter-free calculations are in good agreement with experimental measurements demonstrating the predictive power of this approach. Phonon-defect scattering rates are computed from a Green's function methodology that is nonperturbative and includes interatomic force constant variance induced near the defects. Restricting calculations to first-order perturbation approaches can overestimate optic phonon scattering rates by nearly three orders of magnitude. On the other hand, neglecting the force variance weakens the scattering rates by about an order of magnitude, mostly in the low-frequency region below 2 THz. Here, this work elucidates important properties of phonon-defect scattering in thermal transport and demonstrates the predictive power of the coupling of Peierls-Boltzmann transport, Green's function methods, and density functional theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Phonon thermal conductance across GaN-AlN interfaces from first principles

The vibrational thermal conductances (G) across GaN-AlN interfaces are computed using a nonequilibrium Green's function formalism in the harmonic limit with bulk and interfacial interatomic force constants (IFCs) fully from density functional theory. Several numerical methods and supercell configurations are employed to examine the sensitivity of G to variances of IFCs. In particular, the effects of supercell size, the enforcement of symmetry constraints, and truncation of IFCs near the interface, and atomic relaxation on phonon transmission and conductance are explored. Our fully first-principles calculations are compared with common approximations and measured G values inferred from thermal conductivity measurements for GaN-AlN superlattices. Here, our calculated value, G~300MWm -2 K -1 , is nearly half that from measurements. This discrepancy is critically analyzed in terms of the physical assumptions of the calculations and the derivation of the experimental values. This work provides guidelines to determine “physically correct” sets of interfacial IFCs from first principles for thermal conductance calculations using minimal computational resources. It also contributes toward developing predictive calculations and a more complete picture of thermal conduction across interfaces, a step toward first-principles multiscale thermal transport.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Deep Potential Molecular Dynamics: A Scalable Model with the Accuracy of Quantum Mechanics

We introduce a scheme for molecular simulations, the deep potential molecular dynamics (DPMD) method, based on a many-body potential and interatomic forces generated by a carefully crafted deep neural network trained with ab initio data. The neural network model preserves all the natural symmetries in the problem. It is first-principles based in the sense that there are no ad hoc components aside from the network model. We show that the proposed scheme provides an efficient and accurate protocol in a variety of systems, including bulk materials and molecules. In all these cases, DPMD gives results that are essentially indistinguishable from the original data, at a cost that scales linearly with system size.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2019

Coupling 2D atomistic information to 3D kink-pair enthalpy models of screw dislocations in bcc metals

The kink-pair activation enthalpy is a fundamental parameter in the theory of plasticity of body-centered cubic (bcc) metals. It controls the thermally activated motion of screw dislocation at low and intermediate temperatures. While direct atomistic calculations of kink pairs on screw dislocations have reached a high degree of accuracy, they can only be practically performed using semiempirical interatomic force fields, as electronic structure methods have not yet reached the level of efficiency needed to capture the system sizes required to model kink-pair structures. In this context, an alternative approach based on standard three-dimensional elastic models, which are efficient but lack atomic-level information, coupled to a substrate potential that represents the underlying lattice, has been widely applied over the past few years. This class of methods, known as ‘line-on-substrate’ (LOS) models, uses the substrate potential to calculate the lattice contribution to the kink-pair energies. In this work, we introduce the stress dependence of the substrate potential into LOS models to evaluate its impact on kink-pair energies. In addition, we include asymmetric dislocation core energies as an extra descriptor of the dislocation character. This asymmetry is also elevated to the continuum level by adding core energies to the general LOS formulation and used to explain potential energy differences known to exist between left and right kinks in bcc metals. Finally, more importantly, by matching the total LOS energies to previously calculated atomistic energies of kink-pair configurations, we issue a rule to establish the value of the so-called core width in nonsingular elasticity theories and reduce its arbitrariness as a mathematical construct.

36 MATERIALS SCIENCE↗
DOE OSTI2024

First-principles determination of the phonon-point defect scattering and thermal transport due to fission products in ThO 2

This work presents the first-principles calculations of the lattice thermal conductivity degradation due to point defects in thorium dioxide using an iterative solution of the Pierels-Boltzmann transport equation. We have used the non perturbative Green’s function methodology to compute the phonon- point defect scattering rates that consider the local distortion around the point defect, including the mass difference changes, interatomic force constants, and structural relaxation near the point defects. The point defects considered in this work include the vacancy of thorium ($V$ $TH$ ) and oxygen ($V$ $O$ ), substitutions of helium ($He$ $Th$ ), krypton ($Kr$ $Th$ ), zirconium ($Zr$ $Th$ ), iodine ($I$ $Th$ ) and xenon ($Xe$ $Th$ ) in the thorium site, and the three different configurations of the Schottky defects. The results of the phonon-defect scattering rate reveal that, among all the considered intrinsic defects, the thorium vacancy and helium substitution in the thorium site scatter the phonon most due to the substantial changes in the force constant and structural distortions. The scattering of phonons due to the substitutional defects unveils that the zirconium atom scatters phonons the least, followed by xenon, iodine, krypton, and helium. This is contrary to the intuition that the scattering strength follows $He$ $Th$ > $Kr$ $Th$ > $Zr$ $Th$ > $I$ $Th$ > $Xe$ $Th$ based on the mass difference. This striking difference in the zirconium phonon scattering is due to the local chemical environment changes. Zirconium is an electropositive element with valency similar to thorium and, therefore, can bond with the oxygen atoms, thus creating less force constant variance compared to iodine, an electronegative element noble gas helium, xenon, and krypton. In conclusion, these results can serve as the benchmark for the analytical models and help the engineering-scale modeling effort for nuclear design.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Scalable molecular dynamics with NAMD on the Summit system

According to the article, NAnoscale Molecular Dynamics (NAMD) is a parallel molecular dynamics application that has been used to make breakthroughs in understanding the structure and dynamics of large biomolecular complexes, such as viruses like HIV and various types of influenza. State-of-the-art biomolecular simulations often require integration of billions of timesteps, computing all interatomic forces for each femtosecond timestep. Molecular dynamics simulation of large biomolecular systems and long-timescale biological phenomena requires tremendous computing power. NAMD harnesses the power of thousands of heterogeneous processors to meet this demand. In this paper, we present algorithmic improvements and performance optimizations that enable NAMD to achieve high performance on the IBM Newell platform (with IBM POWER9 processors and NVIDIA Volta V100 GPUs), which underpins the Oak Ridge National Laboratory's Summit and Lawrence Livermore National Laboratory's Sierra supercomputers. The Top-500 supercomputers November 2018 list shows Summit at the number one spot, with 200 petaflop/s peak performance, and Sierra second with 125 petaflop/s. Optimizations for NAMD on Summit include: data layout changes for GPU acceleration and CPU vectorization, improving GPU offload efficiency, increasing performance with Parallel Active Messaging Interface support in Charm++, improving efficiency of fast Fourier transform calculations, improving load balancing, enabling better CPU vectorization and cache performance, and providing an alternative thermostat through stochastic velocity rescaling. We also present performance scaling results on early Newell systems.

Computer Science↗