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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Nuclear energy density functionals grounded in ab initio calculations

Here, we discuss the construction of a nuclear energy density functional (EDF) from ab initio computations and advocate the need for a methodical approach that is free from ad hoc assumptions. The equations of state (EoSs) of symmetric nuclear and pure neutron matter are computed using the chiral NNLO sat and the phenomenological AV4' + UIX c Hamiltonians as inputs to self-consistent Green's function (SCGF) and auxiliary field diffusion Monte Carlo (AFDMC) methods. We propose a convenient parametrization of the EoS as a function of the Fermi momentum and fit it on the SCGF and AFDMC calculations. We apply the ab initio based EDF to carry out an analysis of the binding energies and charge radii of different nuclei in the local density approximation. The NNLO sat -based EDF produces encouraging results, whereas the AV4' + UIX c -based one is farther from experiment. Possible explanations of these different behaviors are suggested, and the importance of gradient and spin-orbit terms is analyzed. Our paper paves the way for a practical and systematic way to merge ab initio nuclear theory and density functional theory, while shedding light on some critical aspects of this procedure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Automatic Prediction of Surface Phase Diagrams Using Ab Initio Grand Canonical Monte Carlo

The properties of a material are often strongly influenced by its surfaces. Depending on the nature of the chemical bonding in a material, its surface can undergo a variety of stabilizing reconstructions that dramatically alter the chemical reactivity, light absorption, and electronic band offsets. For decades, ab initio thermodynamics has been the method of choice for computationally determining the surface phase diagram of a material under different conditions. The surfaces considered for these studies, however, are often human-selected and too few in number, leading both to insufficient exploration of all possible surfaces and to biases toward portions of the composition–structure phase space that often do not encompass the most stable surfaces. To overcome these limitations and automate the discovery of realistic surfaces, we combine density functional theory and grand canonical Monte Carlo (GCMC) into “ab initio GCMC.” This paper presents the successful application of ab initio GCMC to the study of oxide overlayers on Ag(111), which, for many years, mystified experts in surface science and catalysis. Specifically, we report that ab initio GCMC is able to reproduce the surface phase diagram of Ag(111) with no preconceived notions about the system. Using nonlinear, random forest regression, we discover that the Ag coordination number with O and the surface O–Ag–O bond angles are good descriptors of the surface energy. Additionally, using the composition–structure evolution histories produced by ab initio GCMC, we deduce a mechanism for the formation of oxide overlayers based on the Ag 3 O 4 pyramid motif that is common to many reconstructions of Ag(111). Lastly, ab initio GCMC is a promising tool for the discovery of realistic surfaces that can then be used to study phenomena on complex surfaces such as heterogeneous catalysis and material growth, enabling reliable and insightful interpretations of experiments.

36 MATERIALS SCIENCE↗
NASA NTRS1996

Ab Initio Calculation of Accurate Vibrational Frequencies for Molecules of Interest in Atmospheric Chemistry

Due to advances in quantum mechanical methods over the last few years, it is now possible to determine ab initio potential energy surfaces in which fundamental vibrational frequencies are accurate to within +/- 8 cm(sup -1) on average, and molecular bond distances are accurate to within +/- 0.001-0.003 A, depending on the nature of the bond. That is, the potential energy surfaces have not been scaled or empirically adjusted in any way, showing that theoretical methods have progressed to the point of being useful in analyzing spectra that are not from a tightly controlled laboratory environment, such as rovibrational spectra from the interstellar medium. Some recent examples demonstrating this accuracy win be presented and discussed. These include the HNO, CH4, C2H4, and ClCN molecules. The HNO molecule is interesting due to the very large H-N anharmonicity, while ClCN has a very large Fermi resonance. The ab initio studies for the CH4 and C2H4 molecules present the first accurate full quartic force fields of any kind (i.e., whether theoretical or empirical) for a five-atom and six-atom system, respectively.

Lee, Timothy J.↗
DOE OSTI2016

Accurate determination of the valence band edge in hard x-ray photoemission spectra using GW theory

We introduce a new method for determining accurate values of the valence-band maximum in x-ray photoemission spectra. Specifically, we align the sharpest peak in the valence-band region of the experimental spectrum with the corresponding feature of a theoretical valence-band density of states curve from ab initio GW theory calculations. This method is particularly useful for soft and hard x-ray photoemission studies of materials with a mixture of valence-band characters, where strong matrix element effects can render standard methods for extracting the valence-band maximum unreliable. We apply our method to hydrogen-terminated boron-doped diamond, which is a promising substrate material for novel solar cell devices. By carrying out photoemission experiments with variable light polarizations, we verify the accuracy of our analysis and the general validity of the method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Bayesian Model Averaging for Ensemble-Based Estimates of Solvation Free Energies

This paper applies the Bayesian Model Averaging (BMA) statistical ensemble technique to estimate small molecule solvation free energies. There is a wide range methods for predicting solvation free energies, ranging from empirical statistical models to {\it ab initio} quantum mechanical approaches. Each of these methods are based on a set of conceptual assumptions that can affect a method's predictive accuracy and transferability. Using an iterative statistical process, we have selected and combined solvation energy estimates using an ensemble of 17 diverse methods from the SAMPL4 blind prediction study to form a single, aggregated solvation energy estimate. The ensemble design process evaluates the statistical information in each individual method as well as the performance of the aggregate estimate obtained from the ensemble as a whole. Methods that possess minimal or redundant information are pruned from the ensemble and the evaluation process repeats until aggregate predictive performance can no longer be improved. We show that this process results in a final aggregate estimate that outperforms all individual methods by reducing estimate errors by as much as 91% to 1.2 kcal/mol accuracy. We also compare our iterative refinement approach to other statistical ensemble approaches and demonstrate that this iterative process reduces estimate errors by as much as 61%. This work provides a new approach for accurate solvation free energy prediction and lays the foundation for future work on aggregate models that can balance computational cost with predictive accuracy.

solvation, Bayesian model averaging, Electrostatic↗
DOE OSTI2015

Structural stability and defect energetics of ZnO from diffusion quantum Monte Carlo

We have applied the many-body ab initio diffusion quantum Monte Carlo (DMC) method to study Zn and ZnO crystals under pressure and the energetics of the oxygen vacancy, zinc interstitial, and hydrogen impurities in ZnO. We show that DMC is an accurate and practical method that can be used to characterize multiple properties of materials that are challenging for density functional theory (DFT) approximations. DMC agrees with experimental measurements to within 0.3 eV, including the band-gap of ZnO, the ionization potential of O and Zn, and the atomization energy of O{sub 2}, ZnO dimer, and wurtzite ZnO. DMC predicts the oxygen vacancy as a deep donor with a formation energy of 5.0(2) eV under O-rich conditions and thermodynamic transition levels located between 1.8 and 2.5 eV from the valence band maximum. Our DMC results indicate that the concentration of zinc interstitial and hydrogen impurities in ZnO should be low under n-type and Zn- and H-rich conditions because these defects have formation energies above 1.4 eV under these conditions. Comparison of DMC and hybrid functionals shows that these DFT approximations can be parameterized to yield a general correct qualitative description of ZnO. However, the formation energy of defects in ZnO evaluated with DMC and hybrid functionals can differ by more than 0.5 eV.

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

Probing the Origin of Interfacial Carriers in SrTiO 3 –LaCrO 3 Superlattices

Emergent phenomena at complex oxide interfaces could provide the basis for a wide variety of next-generation devices, including photovoltaics and spintronics. To date, detailed characterization and computational modeling of interfacial defects, cation intermixing, and film stoichiometry have helped to explain many of the novel behaviors observed at a single heterojunction. Unfortunately, many of the techniques employed to characterize a single heterojunction are less effective for a superlattice made up of a repeating series of interfaces that induce collective interfacial phenomena throughout a film. However, these repeating interfaces present an untapped opportunity to introduce an additional degree of complexity, such as confined electric fields, that cannot be studied in a single heterojunction. In this work, we explore the properties of SrTiO3-LaCrO3 superlattices to understand the role of defects, including variations in cation stoichiometry of individual layers of the superlattice, intermixing across interfaces, and interfacial oxygen vacancies. Using x-ray photoelectron spectroscopy (XPS) and scanning transmission electron microscopy electron energy-loss spectroscopy (STEM-EELS), we quantify the stoichiometry of individual layers of the superlattice and determine the degree of intermixing in these materials. By comparing these results to both density functional theory (DFT) models and STEM-EELS measurements of the Ti and Cr valence in each layer of the superlattice, we correlate different types of defects with the associated materials properties of the superlattice. We show that a combination of ab initio modeling and complementary structural characterization methods can offer unique insight into structure-property relationships in many oxide superlattice systems.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Absorption Spectra of Acetylene, Vinylidene, and Their Deuterated Isotopologues on Ab Initio Potential Energy and Dipole Moment Surfaces

Here, the absorption spectra of acetylene (HCCH) and vinylidene (H 2 CC) as well as their deuterated isotopologues are investigated theoretically on a near spectroscopically accurate full-dimensional potential energy surface reported in an earlier publication, using dipole moment surfaces reported in this work, which are constructed with a neural network method from a large number of ab initio data points. These global surfaces cover not only the deep acetylene well, but also the vinylidene well, as well as the transition region between the two isomers. The agreement with available experimental data for acetylene is excellent, validating both the potential energy surface and the dipole moment surfaces. The infrared spectra of vinylidene and its deuterated isotopologues are predicted. The potential and dipole moment surfaces lay the foundation for future spectroscopic studies of the acetylene-vinylidene isomerization involving large amplitude motions.

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

Communication: Hole localization in Al-doped quartz SiO{sub 2} within ab initio hybrid-functional DFT

We investigate the long-standing problem of hole localization at the Al impurity in quartz SiO{sub 2}, using a relatively recent DFT hybrid-functional method in which the exchange fraction is obtained ab initio, based on an analogy with the static many-body COHSEX approximation to the electron self-energy. As the amount of the admixed exact exchange in hybrid functionals has been shown to be determinant for properly capturing the hole localization, this problem constitutes a prototypical benchmark for the accuracy of the method, allowing one to assess to what extent self-interaction effects are avoided. We obtain good results in terms of description of the charge localization and structural distortion around the Al center, improving with respect to the more popular B3LYP hybrid-functional approach. We also discuss the accuracy of computed hyperfine parameters, by comparison with previous calculations based on other self-interaction-free methods, as well as experimental values. We discuss and rationalize the limitations of our approach in computing defect-related excitation energies in low-dielectric-constant insulators.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1990

Ab initio R-matrix calculations of e+-molecule scattering

The adaptation of the molecular R-matrix method, originally developed for electron-molecule collision studies, to positron scattering is discussed. Ab initio R-matrix calculations are presented for collisions of low energy positrons with a number of diatomic systems including H2, HF and N2. Differential elastic cross sections for positron-H2 show a minimum at about 45 deg for collision energies between 0.3 and 0.5 Ryd. The calculations predict a bound state of positronHF. Calculations on inelastic processes in N2 and O2 are also discussed.

Danby, Grahame↗
DOE OSTI2023

Microstructure modeling of nuclear structural materials: Recent progress and future directions

Modeling and simulation of microstructures are essential to understand the complex responses and behaviors of nuclear materials in extreme environments. The needs to assess the extended life operation as well as the growing interest in accelerating nuclear materials development and qualification have stimulated the use of high-fidelity multiscale models aided by empirical and ab initio data. This paper reviews the role of various models across different length and time scales in investigating irradiation effects on microstructure evolution and degradation, in particular the embrittlement caused by radiation induced or enhanced formation of nanoscale chemical heterogeneities. The strength and limitations of these models, including classical rate theories, cluster dynamics, phase-field methods, and atomistic models informed by ab initio energies, are discussed with seminal examples. Challenges regarding the lack of thermo-kinetic data and theoretical treatments considering chemical complexities and magnetic excitations, as well as the stabilizing effect by excess point defects in nuclear structural materials are presented, along with potential solutions based on ab initio informed surrogate energy models and statistical sampling by Monte Carlo simulations. Further, the review then highlights the opportunities to leverage the advantages of different methods by establishing hybrid models by shared variables or coupled codes and applications. Finally, the review concludes with forward-looking remarks on how the use of physics-based models can aid the improvement of machine-learning models of property degradation and vice versa.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Spin-orbit coupling induced demagnetization in Ni: Ab initio nonadiabatic molecular dynamics perspective

Spin-orbit coupling (SOC), which can induce spin flip during the relaxation of photoexcited charge carrier, plays a crucial role in spin dynamics. Here in this work, we have used time-domain ab initio nonadiabatic molecular dynamics (NAMD) method to study the SOC induced ultrafast demagnetization in Ni at 300K. The spin-diabatic representation using spin-polarized Kohn-Sham (KS) basis sets and spin-adiabatic representation using spinor basis sets have been applied, and both of them achieve demagnetization in Ni with a timescale around 100fs. The spin-diabatic representation suggests a picture that the electron-phonon coupling (EPC) provides direct energy relaxation channel among the same-spin states, while the SOC can induce spin flip. After photoexcitation, it is found the spin-minority electrons relax to the same-spin states rather than the opposite-spin states, since EPC is larger than SOC by one order of magnitude. By contrast, for the spin-majority electrons, spin flip occurs since there are no empty same-spin states as electron acceptor above the Fermi level. The different relaxation pathways for spin-majority and spin-minority electrons induce the demagnetization. The spin-adiabatic representation provides an Elliott-Yafet spin-phonon scattering picture. The SOC induced reduction of magnetic moment in Ni may induce magnon to drive further demagnetization. The ab initio NAMD simulation provides a critical angle to understand how the SOC and EPC affect demagnetization process in Ni.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Influence of Cation Size on the Electronic and Magnetic Structure of M CoO{sub 3} (M = Bi, La) from Ab Initio Calculations

Using the first-principle DFT method and LSDA+ U calculations, we have performed a detailed investigation on the structural and magnetic ground states of strongly correlated M CoO{sub 3} (M= Bi, La) perovskites. Magneticphase stabilities, structural behaviors and electronic properties of both materials have been investigated. LSDA+ U calculations show that the antiferromagnetic (AFM)-C phase in the tetragonal La/BiCoO{sub 3} is energetically more favorable than the other configurations, and both rhombohedral materials are nonmagnetic. All structural data of tetragonal BiCoO{sub 3} and rhombohedral LaCoO{sub 3} compounds are found in excellent agreement with the experiment for U{sub eff}= 6 eV. We have predicted also band gaps and magnetic moments with a good accuracy via the band structure calculations. From the total and partial densities of states, we have shown the important role of O p Co d and Bi/La s-O p hybridizations in both systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Optically isotropy in scintillator host compounds M{sub 2}LaCl{sub 5} (M=Rb, and Cs): Ab-initio study

Full potential linearized augmented plane wave method (FP-LAPW) has been used to calculate the electronic structure and optical properties of high light output scintillator host compounds M{sub 2}LaCl{sub 5} (M=Rb, and Cs) with in the Tran Blaha modified Becke Johnson (TB-mBJ) potential. These are wide band gap materials with the band gap of 4.75, and 4.72 eV for Rb{sub 2}LaCl{sub 5}, and Cs{sub 2}LaCl{sub 5} respectively. From the calculated optical properties of these compounds, we find these compounds to be optically isotropic, though they are structurally anisotropic, which is an important criteria for the ceramic scintillators.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Defect energetics of concentrated solid-solution alloys from ab initio calculations: Ni 0.5 Co 0.5 , Ni 0.5 Fe 0.5 , Ni 0.8 Fe 0.2 and Ni 0.8 Cr 0.2

It has been shown that concentrated solid solution alloys possess unusual electronic, magnetic, transport, mechanical and radiation-resistant properties that are directly related to underlying chemical complexity. Because every atom experiences a different local atomic environment, the formation and migration energies of vacancies and interstitials in these alloys exhibit a distribution, rather than a single value as in a pure metal or dilute alloy. In this study, using ab initio calculations based on density functional theory and special quasirandom structure, we have characterized the distribution of defect formation energy and migration barrier in four Ni-based solid-solution alloys: Ni 0.5 Co 0.5 , Ni 0.5 Fe 0.5 , Ni 0.8 Fe 0.2 and Ni 0.8 Cr 0.2 . As defect formation energies in finite-size models depend sensitively on the elemental chemical potential, we have developed a computationally efficient method for determining it which takes into account the global composition and the local short-range order. In addition we have compared the results of our ab initio calculations to those obtained from available embedded atom method (EAM) potentials. Our results indicate that the defect formation and migration energies are closely related to the specific atomic size in the structure, which further determines the elemental diffusion properties. In conclusion, different EAM potentials yield different features of defect energetics in concentrated alloys, pointing to the need for additional potential development efforts in order to allow spatial and temporal scale-up of defect and simulations, beyond those accessible to ab initio methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2001

Ab-Initio Description and Prediction of Properties of Carbon-Based and Other Non-Metallic Materials

We have resolved the long-standing problem consisting of 30%-50% theoretical underestimates of the band gaps of non-metallic materials. We describe the Bagayoko, Zhao, and Williams (BZW) method that rigorously circumvents the basis-set and variational effect presumed to be a cause of these underestimates. We present ab-initio, computational results that are in agreement with experiment for diamond (C), silicon (Si), silicon carbides (3C-SiC and 4H-SiC), and other semiconductors (GaN, BaTiO3, AlN, ZnSe, ZnO). We illustrate the predictive capability of the BZW method in the case of the newly discovered cubic phase of silicon nitride (c-Si3N4) and of selected carbon nanotabes [(10,0), and (8,4)]. Our conclusion underscores the inescapable need for the BZW method in ab-initio calculations that employ a basis set in a variational approach. Current nanoscale trends amplify this need. We estimate that the potential impact of applications of the BZW method in advancing our understanding of nonmetallic materials, in informing experiment, and particularly in guiding device design and fabrication is simply priceless.

Bagayoko, D.↗
DOE OSTI2015

Infrared signature of micro-hydration in the organophosphate sarin: An ab initio study

The infrared (IR) spectra of micro-hydrated Sarin•(H 2 O) n clusters containing between one and four explicit waters have been studied using ab initio density functional theory (DFT) methods. The phosphate group P=O bond vibration region (~1270 to 1290 cm –1 ) revealed the largest frequency variation with hydration, with a frequency red shift reflecting the direct hydrogen bond formation between the P=O of Sarin and water. Small variations to the P-F stretch (~810 to 815 cm –1 ) and the C-O-P vibrational modes (~995 to 1004 cm –1 ) showed that the water interactions with these functional groups were minor, and that the structures of Sarin were not extensively perturbed in the hydrated complexes. Increasing the number of explicit hydration waters produced only small vibrational changes in the lowest free energy complexes. These minor changes were consistent with a single water-phosphate hydrogen bond being the dominant structure, though a second water-phosphate hydrogen bond was observed in some complexes and was identified by an additional red shift of the P=O bond vibration. As a result, the H 2 O•H 2 O vibrational modes (~3450 to 3660 cm –1 ) increased in complexity with higher hydration levels and reflect the extended hydrogen bonding networks formed between the explicit waters in the hydrated Sarin clusters.

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

Ab initio calculations on the magnetic properties of transition metal complexes

We present a protocol for the ab initio determination of the magnetic properties of mono- and polynuclear transition metal compounds. First, we obtain the low lying electronic states by multireference methods. Then, we include spin-orbit coupling and an external magnetic field for the determination of zero-field splitting and g-tensors. For the polynuclear complexes the magnetic exchange coupling constants are determined by a modified complete active space self consistent field method. Based on the results of the ab initio calculations, magnetic data such as magnetic susceptibility or magnetization are simulated and compared to experimental data. The results obtained for the polynuclear complexes are further analysed by calculations on model complexes where part of the magnetic centers are substituted by diamagnetic ions. The methods are applied to different Co and Ni containing transition metal complexes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗