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,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Thermophysical properties of liquid Ni around the melting temperature from molecular dynamics simulation

Thermophysical properties of liquid nickel (Ni) around the melting temperature are investigated by means of classical molecular dynamics (MD) simulation, using three different embedded atom method potentials to model the interactions between the Ni atoms. Melting temperature, enthalpy, static structure factor, self-diffusion coefficient, shear viscosity, and thermal diffusivity are compared to recent experimental results. Using ab initio MD simulation, we also determine the static structure factor and the mean-squared displacement at the experimental melting point. For most of the properties, excellent agreement is found between experiment and simulation, provided the comparison relative to the corresponding melting temperature. We discuss the validity of the Hansen-Verlet criterion for the static structure factor as well as the Stokes-Einstein relation between self-diffusion coefficient and shear viscosity. The thermal diffusivity is extracted from the autocorrelation function of a wavenumber-dependent temperature fluctuation variable.

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

Charge transport in organic molecular semiconductors from first principles: The bandlike hole mobility in a naphthalene crystal

Predicting charge transport in organic molecular crystals is notoriously challenging. Carrier mobility calculations in organic semiconductors are dominated by quantum chemistry methods based on charge hopping, which are laborious and only moderately accurate. We compute from first principles the electron-phonon scattering and the phonon-limited hole mobility of naphthalene crystal in the framework of ab initio band theory. Our calculations combine GW electronic bandstructures, ab initio electron-phonon scattering, and the Boltzmann transport equation. The calculated hole mobility is in very good agreement with experiment between 100–300K, and we can predict its temperature dependence with high accuracy. Here we show that scattering between intermolecular phonons and holes regulates the mobility, though intramolecular phonons possess the strongest coupling with holes. We revisit the common belief that only rigid molecular motions affect carrier dynamics in organic molecular crystals. Our paper provides a quantitative and rigorous framework to compute charge transport in organic crystals and is a first step toward reconciling band theory and carrier hopping computational methods.

36 MATERIALS SCIENCE↗
DOE OSTI2021

MULTICOM2 open-source protein structure prediction system powered by deep learning and distance prediction

Protein structure prediction is an important problem in bioinformatics and has been studied for decades. However, there are still few open-source comprehensive protein structure prediction packages publicly available in the field. In this paper, we present our latest open-source protein tertiary structure prediction system—MULTICOM2, an integration of template-based modeling (TBM) and template-free modeling (FM) methods. The template-based modeling uses sequence alignment tools with deep multiple sequence alignments to search for structural templates, which are much faster and more accurate than MULTICOM1. The template-free (ab initio or de novo) modeling uses the inter-residue distances predicted by DeepDist to reconstruct tertiary structure models without using any known structure as template. In the blind CASP14 experiment, the average TM-score of the models predicted by our server predictor based on the MULTICOM2 system is 0.720 for 58 TBM (regular) domains and 0.514 for 38 FM and FM/TBM (hard) domains, indicating that MULTICOM2 is capable of predicting good tertiary structures across the board. It can predict the correct fold for 76 CASP14 domains (95% regular domains and 55% hard domains) if only one prediction is made for a domain. The success rate is increased to 3% for both regular and hard domains if five predictions are made per domain. Moreover, the prediction accuracy of the pure template-free structure modeling method on both TBM and FM targets is very close to the combination of template-based and template-free modeling methods. This demonstrates that the distance-based template-free modeling method powered by deep learning can largely replace the traditional template-based modeling method even on TBM targets that TBM methods used to dominate and therefore provides a uniform structure modeling approach to any protein. Finally, on the 38 CASP14 FM and FM/TBM hard domains, MULTICOM2 server predictors (MULTICOM-HYBRID, MULTICOM-DEEP, MULTICOM-DIST) were ranked among the top 20 automated server predictors in the CASP14 experiment. After combining multiple predictors from the same research group as one entry, MULTICOM-HYBRID was ranked no. 5. The source code of MULTICOM2 is freely available at https://github.com/multicom-toolbox/multicom/tree/multicom_v2.0 .

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Advancements in Validation of TSLs through Inelastic Neutron Scattering and Transmission Measurements [Abstract]

Historically, the free gas approximation has been used to treat the thermal scattering of neutrons with energies below a few electron-volts (eV) in unevaluated materials. However, this method inadequately reproduces neutron scattering at these energies. Until recently, only a limited number of materials had available thermal scattering law (TSL) files/libraries in the ENDF nuclear data libraries in this energy range. With advancements in atomistic modeling techniques, such as molecular dynamics, ab-initio molecular dynamics, and density functional theory, TSL libraries have become available for many more materials. This is particularly relevant due to the rising interest in several advanced reactor systems that require novel moderator and reflector materials. While quasi-integral and integral benchmarks have been designed to validate historically important moderator materials (such as light water and polyethylene), there is currently a lack of standard validation methods for TSLs, especially when multiple conflicting TSLs exist. To address this issue, the Oak Ridge National Lab Nuclear Data group has been working on utilizing inelastic neutron scattering (INS) measurements combined with transmission (i.e., total cross section) measurements to evaluate and validate TSLs for different materials. We plan to demonstrate how this method has worked on materials such as polyethylene, lucite, and polystyrene. In addition, we will compare the newly created libraries to ENDF libraries for these materials and explain why integral benchmarks should not be used for validation when multiple TSLs exist.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Towards Ab-Initio Simulations of Crystalline Defects at the Exascale Using Spectral Quadrature Density Functional Theory

Defects in crystalline solids play a crucial role in determining properties of materials at the nano, meso- and macroscales, such as the coalescence of vacancies at the nanoscale to form voids and prismatic dislocation loops or diffusion and segregation of solutes to nucleate precipitates, phase transitions in magnetic materials via disorder and doping. First principles Density Functional Theory (DFT) simulations can provide a detailed understanding of these phenomena. However, the number of atoms needed to correctly simulate these systems is often beyond the reach of many widely used DFT codes. The aim of this article is to discuss recent advances in first principles modeling of crystal defects using the spectral quadrature method. The spectral quadrature method is linear scaling with respect to the number of atoms, permits spatial coarse-graining, and is capable of simulating non-periodic systems embedded in a bulk environment, which allows the application of appropriate boundary conditions for simulations of crystalline defects. In this article, we discuss the state-of-the-art in ab-initio modeling of large metallic systems of the order of several thousand atoms that are suitable for utilizing exascale computing resourses.

36 MATERIALS SCIENCE↗
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 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 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↗
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 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 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↗
DOE OSTI2016

Electron transmission through a class of anthracene aldehyde molecules

Transmission of electrons via metal-molecule-metal junctions, involving rotor-stator anthracene aldehyde molecules is investigated. Two model barriers having input parameters evaluated from accurate ab initio calculations are proposed and the transmission coefficients are obtained by using the quasiclassical approximation. Transmission coefficients further enter in the integral for the net current, utilizing Simmons’ method. Conformational dependence of the tunneling processes is evident and the presence of the side groups enhances the functionality of the future single-molecule based electronic devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗