Search NASA⌕ Search

SEARCH · Search NASA

Results for “ab initio methods”

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 1,351 records · Page 75
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Exploring structural transitions at grain boundaries in Nb using a generalized embedded atom interatomic potential

The advancement in experimental techniques, like the atom probe tomography and high resolution electron microscopy, is fueling interest in studying structural transformations of grain boundaries in metal and alloys to uncover correlations between mechanical properties and solute or impurity segregation to grain boundaries. Atomistic modeling is an important tool that can pinpoint the intricate dynamics of grain boundary phase transitions, but the lack of accurate interatomic potentials needed to simulate the complex dynamics of grain boundary structural transitions and identify different metastable phases has been the bottleneck. To this end, we use niobium as a model body centered cubic (BCC) metal and develop an interatomic potential to study grain boundary phase transitions. The potential for Nb is based on a generalization of the embedded atomic method potential and has sufficient flexibility to learn complex energy landscapes using a small set of training structures. We systematically test and validate the using data from ab initio density functional theory calculations and experiments. Using this potential, we calculate energies of multiple symmetric-tilt grain boundaries spanning a wide range of misorientation angles. Additionally, we explore different metastable structures of the Σ 27(552) [$1\overline{1}0$] grain boundary and use molecular dynamic simulations to study the coexistence of metastable phases and grain boundary transitions at finite temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Ab initio study and assignment of electronic states in molecular RaCl

Radium compounds have attracted recently considerable attention due to both development of experimental techniques for high-precision laser spectroscopy of molecules with short-lived nuclei and amenability of certain radium compounds for direct cooling with lasers. Currently, radium monofluoride (RaF) is one of the most studied molecules among the radium compounds, both theoretically and recently also experimentally. Complementary studies of further diatomic radium derivatives are highly desired to assess the influence of chemical substitution on diverse molecular parameters, especially on those connected with laser cooling, such as vibronic transition probabilities, and those related to violations of fundamental symmetries. Here in this article high-precision ab initio studies of electronic and vibronic levels of diatomic radium monochloride (RaCl) are presented. Recently developed approaches for treating electronic correlation with Fock-space coupled cluster methods are applied for this purpose. Theoretical results are compared to an early experimental investigation by Lagerqvist and used to partially reassign the experimentally observed transitions and molecular electronic levels of RaCl. Effective constants of $\mathscr{P}$-odd hyperfine interaction W a and $\mathscr{P, T}$-odd scalar-pseudoscalar nucleus-electron interaction Ws in the ground electronic state of RaCl are estimated within the framework of a quasirelativistic Zeroth-Order Regular Approximation approach and compared to parameters in RaF and RaOH.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2014

Precise Lifetime Measurements in Light Nuclei for Benchmarking Modern Ab-initio Nuclear Structure Models

A new generation of ab-initio calculations, based on realistic two- and three-body forces, is having a profound impact on our view of how nuclei work. To improve the numerical methods, and the parameterization of 3-body forces, new precise data are needed. Electromagnetic transitions are very sensitive to the dynamics which drive mixing between configurations. We have made a series of precise (< 3%) measurements of electromagnetic transitions in the A=10 nuclei {sup 10}C and {sup 10}Be by using the Doppler Shift Attenuation method carefully. Many interesting features can be reproduced including the strong α clustering. New measurements on {sup 8}Be and {sup 12}Be highlight the interplay between the alpha clusters and their valence neutrons.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Revisiting the K-edge X-ray absorption fine structure of Si, Ge–Si alloys, and the isoelectronic series: CuBr, ZnSe, GaAs, and Ge

Extended X-ray absorption fine structure (EXAFS) has evolved into an unprecedented local-structure technique that is routinely used to study materials’ problems in the biological, chemical, and physical sciences. Like many other experimental techniques, EXAFS also requires that several key atomic parameters must be known a priori before structural information can be quantitatively determined. Utilizing current analytical methods, we revisit the isoelectronic series CuBr, ZnSe, GaAs, and Ge originally studied by Stern et al. during the early development of EXAFS. We demonstrate that the ab initio EXAFS code FEFF accurately predicts the atomic phase shifts and backscattering amplitudes that are primarily functions of the sum of atomic numbers Z along an EXAFS scattering path. We also investigate quantitative fitting and first- and second-shell phase transferability together with problems that arise if a backscattering atom is identified incorrectly in an EXAFS fitting model. Features in the near-edge region, on the other hand, are shown to require a comprehensive treatment of the band structure and density-of-states, including effects of the screened Coulomb interaction between the photoelectron and core hole. Here, we demonstrate that the Bethe–Salpeter equation (BSE) accurately captures the NEXAFS (or XANES) portion of the spectrum for the isoelectronic series in addition to Si and Ge–Si alloys, including within a few eV of the absorption edge, where band structure and excitonic effects are most important.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2011

Relativistic Jahn-Teller effects in the photoelectron spectra of tetrahedral P{sub 4}, As{sub 4}, Sb{sub 4}, and Bi{sub 4}

The group-V tetrahedral cluster cations P{sub 4}{sup +}, As{sub 4}{sup +}, Sb{sub 4}{sup +}, and Bi{sub 4}{sup +} are known to exhibit exceptionally strong Jahn-Teller (JT) effects of electrostatic origin in their {sup 2}E ground states and {sup 2}T{sub 2} excited states. It has been predicted that there exist, in addition, JT couplings of relativistic origin (arising from the spin-orbit (SO) operator) in {sup 2}E and {sup 2}T{sub 2} states of tetrahedral systems, which should become relevant for the heavier elements. In the present work, the JT and SO couplings in the group-V tetramer cations have been analyzed with ab initio relativistic electronic structure calculations. The vibronic line spectra and the band shapes of the photoelectron spectra were simulated with time-dependent quantum wave-packet methods. The results provide insight into the interplay of electrostatic and relativistic JT couplings and SO splittings in the complex photoelectron spectra of these systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Long-range quadrupole electron-phonon interaction from first principles

Lattice vibrations in materials induce perturbations on the electron dynamics in the form of long-range (dipole and quadrupole) and short-range (octopole and higher) potentials. The dipole Fröhlich term can be included in current first-principles electron-phonon (e-ph) calculations and is present only in polar materials. The quadrupole e-ph interaction is present in both polar and nonpolar materials, but currently it cannot be computed from first principles. In this work, we show an approach to compute the quadrupole e-ph interaction and include it in ab initio calculations of e-ph matrix elements. The accuracy of the approach is demonstrated by comparing with direct density functional perturbation theory calculations. We apply our method to silicon as a case of a nonpolar semiconductor and tetragonal PbTiO 3 as a case of a polar piezoelectric material. In both materials we find that the quadrupole term strongly impacts the e-ph matrix elements. Analysis of e-ph interactions for different phonon modes reveals that the quadrupole term mainly affects optical modes in silicon and acoustic modes in PbTiO 3 , although the quadrupole term is needed for all modes to achieve quantitative accuracy. The effect of the quadrupole e-ph interaction on electron scattering processes and transport is shown to be important. Our approach enables accurate studies of e-ph interactions in broad classes of nonpolar, polar, and piezoelectric materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2011

Influence of CH{sub 2} content and network defects on the elastic properties of organosilicate glasses

Organosilicate glasses (OSGs) are used as low-k intermetal dielectrics for advanced integrated circuits. In this application, the material must fulfill two conflicting requirements: It has to have low density to reduce the dielectric constant while being sufficiently mechanically stable to withstand thermomechanical and other stresses during subsequent steps of integrated circuit manufacture. Recent experimental advances in improving the mechanical and electrical properties of these materials have not yet been systematically studied theoretically at the ab initio level due to the large model sizes necessary to realistically describe amorphous materials. In this paper we employ the density-functional-based tight-binding method to achieve an accurate description of OSG properties at different compositions. We analyze the influence of composition and local network defects on the density and bulk modulus of nonporous OSG. We find that the dependence of density and that of mechanical stiffness on chemical composition are of different natures. This difference is traced to a transition between mechanisms of elastic deformation in silica glass and in silicon hydrocarbide, which is also the reason for the two materials' different sensitivities to network defects.

36 MATERIALS SCIENCE↗
NASA NTRS1995

A Correlated Ab Initio Study of Linear Carbon-Chain Radicals C(sub n)H (n=2-7)

Linear carbon-chain radicals C(sub n) H for n = 2-7 have been studied with correlation consistent valence and core-valence basis sets and the coupled cluster method RCCSD(T). Equilibrium structures, rotational constants, and dipole moments are reported and compared with available experimental data. The ground state of the even-n series changes from 2Sigma(+) to 2Pi as the chain is extended. For C4H, the 2Sigma(+) state was found to lie only 72 cm(exp -1) below the 2Pi state in the estimated complete basis set limit for valence correlation. The C2H(-) and C3H(-) anions have also been characterized.

Woon, David E.↗
DOE OSTI2024

Wavefunction matching for solving quantum many-body problems

Ab initio calculations have an essential role in our fundamental understanding of quantum many-body systems across many subfields, from strongly correlated fermions to quantum chemistry and from atomic and molecular systems to nuclear physics. One of the primary challenges is to perform accurate calculations for systems where the interactions may be complicated and difficult for the chosen computational method to handle. Here we address the problem by introducing an approach called wavefunction matching. Wavefunction matching transforms the interaction between particles so that the wavefunctions up to some finite range match that of an easily computable interaction. This allows for calculations of systems that would otherwise be impossible owing to problems such as Monte Carlo sign cancellations. We apply the method to lattice Monte Carlo simulations of light nuclei, medium-mass nuclei, neutron matter and nuclear matter. We use high-fidelity chiral effective field theory interactions and find good agreement with empirical data. These results are accompanied by insights on the nuclear interactions that may help to resolve long-standing challenges in accurately reproducing nuclear binding energies, charge radii and nuclear-matter saturation in ab initio calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Computation of the Thermal Expansion Coefficient of Graphene with Gaussian Approximation Potentials

Direct experimental measurement of thermal expansion coefficient without substrate effects is a challenging task for two-dimensional (2D) materials, and its accurate estimation with large-scale ab initio molecular dynamics is computationally very expensive. Machine learning-based interatomic potentials trained with ab initio data have been successfully used in molecular dynamics simulations to decrease the computational cost without compromising the accuracy. In this work, we investigated using Gaussian approximation potentials to reproduce the density functional theory-level accuracy for graphene within both lattice dynamical and molecular dynamical methods, and to extend their applicability to larger length and time scales. Two such potentials are considered, GAP17 and GAP20. GAP17, which was trained with pristine graphene structures, is found to give closer results to density functional theory calculations at different scales. Further vibrational and structural analyses verify that the same conclusions can be deduced with density functional theory level in terms of the reasoning of the thermal expansion behavior, and the negative thermal expansion behavior is associated with long-range out-of-plane phonon vibrations. Thus, it is argued that the enabled larger system sizes by machine learning potentials may even enhance the accuracy compared to small-size-limited ab initio molecular dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2012

Icosahedral short-range order in amorphous Cu80Si20 by ab initio molecular dynamics simulation study

Short-range order in liquid and amorphous structures of Cu80Si20 is studied by ab initio molecular dynamics simulations. We performed the simulations at 1140 and 300 K respectively to investigate the local structure change from liquid to amorphous. The result of structure factor in comparison with experimental data indicates that our simulation of amorphous Cu80Si20 is reliable. By using the bond-angle distribution function, Honeycutt–Andersen index, Voronoi tessellation method, and the atomistic cluster alignment method, the icosahedral short-range order in the system is revealed. Strong Cu–Si interaction was also observed.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Determination of the electron–phonon coupling constant in tungsten

We used two methods to determine the effective electron-phonon coupling constant (G{sub 0}) in tungsten. Our first principles calculations predict G{sub 0} = 1.65 × 10{sup 17 }W m{sup −3} K{sup −1}. The temporal decay of the femtosecond-resolution optical reflectivity for a (100) surface of bulk W was measured using a pump-probe scheme and analysed using ab initio parameterised two temperature model, which includes both the effects of the electron-phonon coupling and thermal conduction into bulk. This analysis gives G{sub 0} = 1.4(3) × 10{sup 17 }W m{sup −3} K{sup −1}, in good agreement with the theoretical prediction. The described effective method of calculating and measuring G{sub 0} in bulk materials can be easily extended to other metals.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Atomistic k ⋅ p theory

Pseudopotentials, tight-binding models, and k ⋅ p theory have stood for many years as the standard techniques for computing electronic states in crystalline solids. Here, we present the first new method in decades, which we call atomistic k ⋅ p theory. In its usual formulation, k ⋅ p theory has the advantage of depending on parameters that are directly related to experimentally measured quantities, however, it is insensitive to the locations of individual atoms. We construct an atomistic k ⋅ p theory by defining envelope functions on a grid matching the crystal lattice. The model parameters are matrix elements which are obtained from experimental results or ab initio wave functions in a simple way. This is in contrast to the other atomistic approaches in which parameters are fit to reproduce a desired dispersion and are not expressible in terms of fundamental quantities. This fitting is often very difficult. We illustrate our method by constructing a four-band atomistic model for a diamond/zincblende crystal and show that it is equivalent to the sp{sup 3} tight-binding model. We can thus directly derive the parameters in the sp{sup 3} tight-binding model from experimental data. We then take the atomistic limit of the widely used eight-band Kane model and compute the band structures for all III–V semiconductors not containing nitrogen or boron using parameters fit to experimental data. Our new approach extends k ⋅ p theory to problems in which atomistic precision is required, such as impurities, alloys, polytypes, and interfaces. It also provides a new approach to multiscale modeling by allowing continuum and atomistic k ⋅ p models to be combined in the same system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Ab Initio Modeling of fcc Fe-Co-Cr-Ni High Entropy Alloys with Full Composition Range

The Fe-Co-Cr-Ni quaternary system has been studied extensively in the past decade, not only because of the superior properties achieved like high tensile ductility and fracture toughness, but also it is the foundation for the development of FeCrCoNi-based HEAs. However, most of the investigations are mainly focusing on the equiatomic and semi-equiatomic compositions. The physical properties of non-equiatomic compositions are barely explored. In the current work, the ab initio approach is adopted to predict the lattice parameter, structure stabilities, elastic properties, and enthalpy of formations of fcc Fe-Cr-Co-Ni HEA single crystals with the special quasirandom structure (SQS) method. To expand the design of Fe-Co-Cr-Ni HEAs into a full composition range, the current simulations cover all the binaries, ternaries, and the whole quaternary system. The 2-D composition-property contour diagrams are created to show the physical properties in the ternary and quaternary systems. Furthermore, our investigation shows that the fcc alloys' lattice parameter increases with the addition of Fe and Cr and decreases with Ni and Co. In addition to the enthalpy of formation at o K, the shape distortion rate is proposed to be another essential parameter to evaluate crystal stability. The results show the addition of Cr would destabilize the fcc lattice and cause symmetry breaking. Finally, the result of Pugh’s ratio shows that most of the fcc alloys show ductile behavior (P > 1.75), especially for the alloys around Fe 3 Ni, FeCo, and Fe 2 NiCo (P > 2.5). The brittle behavior (P < 1.75) is located around Fe 3 Cr and Fe 2 CoCr. Meanwhile, Cr 2 NiCo and Cr 2 FeNiCo are considered promising compositions to be tried and verified experimentally because of the higher bulk modulus and moderate shear modulus.

36 MATERIALS SCIENCE↗
DOE OSTI2010

Ab-initio calculation of structural, electronic, and optical characterizations of the intermetallic trialuminides ScAl{sub 3} compound

We present first-principles study of the electronic and the optical properties for the intermetallic trialuminides ScAl{sub 3} compound using the full-potential linear augmented plane wave method within density-functional theory. We have employed the generalized gradient approximation (GGA), which is based on exchange-correlation energy optimization to calculate the total energy. Also we have used the Engel-Vosko GGA formalism, which optimizes the corresponding potential for calculating the electronic band structure and optical properties. The electronic specific heat coefficient ({gamma}), which is a function of density of states, can be calculated from the density of states at Fermi energy N(E{sub F}). The N(E{sub F}) of the phase L1{sub 2} is found to be lower than that of D0{sub 22} structure which confirms the stability of L1{sub 2} structure. We found that the dispersion of the band structure of D0{sub 22} is denser than L1{sub 2} phase. The linear optical properties were calculated. The evaluations are based on calculations of the energy band structure. - Graphical abstract: Crystal structure of L1{sub 2} (a) and D0{sub 22} (b) phases of ScAl{sub 3} compound.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Computational understanding of Li-ion batteries

Over the last two decades, computational methods have made tremendous advances, and today many key properties of lithium-ion batteries can be accurately predicted by first principles calculations. For this reason, computations have become a cornerstone of battery-related research by providing insight into fundamental processes that are not otherwise accessible, such as ionic diffusion mechanisms and electronic structure effects, as well as a quantitative comparison with experimental results. The aim of this review is to provide an overview of state-of-the-art ab initio approaches for the modelling of battery materials. Here, we consider techniques for the computation of equilibrium cell voltages, 0-Kelvin and finite-temperature voltage profiles, ionic mobility and thermal and electrolyte stability. The strengths and weaknesses of different electronic structure methods, such as DFT+U and hybrid functionals, are discussed in the context of voltage and phase diagram predictions, and we review the merits of lattice models for the evaluation of finite-temperature thermodynamics and kinetics. With such a complete set of methods at hand, first principles calculations of ordered, crystalline solids, i.e., of most electrode materials and solid electrolytes, have become reliable and quantitative. However, the description of molecular materials and disordered or amorphous phases remains an important challenge. We highlight recent exciting progress in this area, especially regarding the modelling of organic electrolytes and solid–electrolyte interfaces.

25 ENERGY STORAGE↗
DOE OSTI2016

A first-principles study on new high-pressure metastable polymorphs of MoO{sub 2}

The pressure-dependence of the stabilities of several MoO{sub 2} phases has been investigated by density-functional theory (GGA/PBE/PAW). Out of a set of 15 MX{sub 2} structures, the [SnO{sub 2}(II)], [α-PbO{sub 2}], and a modified rutile structure type were identified as possible metastable MoO{sub 2} polymorphs based on the analysis of thermodynamic properties and dynamic stability. High-pressure calculations suggest an orthorhombic TiO{sub 2} structure, dubbed [ortho-TiO{sub 2}], as a high-pressure polymorph at around 25 GPa. Furthermore, we find that the previously reported rutile-type MoO{sub 2} may be understood as a modified rutile type similar to the [VO{sub 2}] structure. - Graphical abstract: First-principles electronic structure and thermochemical calculations reveal four structure candidates as possible metastable polymorphs of MoO{sub 2}. Most promising is a distorted rutile-type similar to the known [VO{sub 2}] structure. An orthorhombic polymorph is proposed as a high-pressure polymorph. Display Omitted - Highlights: • Three possible metastable structure candidates for MoO{sub 2}. • Undistorted rutile type is improbable, a new distorted rutile-type MoO{sub 2} was suggested. • Orthorhombic phase of MoO{sub 2} (ortho-TiO{sub 2} type) should form at 25 GPa. • ab initio thermochemical data provided for MoO{sub 2}.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗