Search NASA⌕ Search

SEARCH · Search NASA

Results for “Ab initio calculations”

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

A Brief Account of Steven Weinberg’s Legacy in ab initio Many-Body Theory

In this contribution to the special issue “Celebrating 30 years of Steven Weinberg’s papers on Nuclear Forces from Chiral Lagrangians,” we emphasize the important role chiral effective field theory has played in leading nuclear physics into a precision era. To this end, we share our perspective on a few of the recent advances made in ab initio calculations of nuclear structure and nuclear matter observables, as well as Bayesian uncertainty quantification of effective field theory truncation errors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Phase field modeling of shearing processes of a dual-lobed γ"|γ'|γ" coprecipitate

The Ni-based superalloy IN718 is one of the most widely used commercial alloy in the aerospace industry since its development in 1960s. The excellent mechanical properties have been attributed in a great deal to strengthening by coherent γ' and γ" precipitates. The deformation mechanisms of these two phases have been well studied individually. Recent experimental characterization has shown coprecipitates of these two phases in a variety of morphologies and it was argued that these coprecipitates may lead to improved strengthening as compared to their monolithic counterparts. However, the deformation mechanisms of these coprecipitates are still not well understood. In this study, we performed microscopic phase field simulations, with generalized-stacking-fault (GSF) energy surfaces from ab initio calculations as inputs, to systematically study the shearing processes of a periodical array of dual-lobed coprecipitates as well as monolithic precipitates. We found that the coupling between the γ' and γ" phases in the coprecipitates forces dislocations to take high energy shearing pathways in both phases that would not occur if they were in monolithic forms. The coupling also creates stacking fault configurations in the coprecipitates that require high stress to form. Thus, the presence of coprecipitates in general should increase the resistance to dislocation shearing and lead to higher strength levels. Various fault configurations observed during the shearing process are documented as a reference for future comparison with experimental observations. The link between stacking fault shearing and microtwinning is also discussed. Here, the mechanisms analyzed in this study deepens our understanding of coprecipitation effects on alloy strength and may form a cornerstone for multi-precipitate strengthened alloy design strategies.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Structural complexity of γ-Al 2 O 3 : The nature of vacancy ordering and the structure of complex antiphase boundaries

The structure of γ-Al2O3 remains largely undetermined despite decades of research. This is due to the high degree of disorder, which poses significant challenges for structural analysis using conventional crystallographic approaches. Herein, we study the structure of γ-Al2O3 with Scanning Transmission Electron Microscopy (STEM) and ab-initio calculations to provide a complete structural description. We show that the microstructure can be understood in terms of two key structural features of nanoscale spinel domains and finite thickness segments termed as complex antiphase boundaries (cAPB) that provide the domain interconnectivity. The spinel domains have a distinctive preference for vacancy ordering, which can be rationalized in terms of a structure with a stacking disorder. Tetragonal P4 1 2 1 2 or monoclinic P2 1 models, all based on the identical motif, can be considered as representative ordered forms. Individual spinel domains are interconnected via cAPBs, which adopt a distinct non-spinel bonding environment of δ-Al 2 O 3 . The most common cAPB consists of a single delta motif with thickness of just 0.6 nm on (001), with the resulting displacement a/4 [101]. Remarkably, the cAPBs are shown to energetically stabilize the spinel domains of γ-Al 2 O 3 explaining their high abundance. We demonstrate how the tetragonal distortions naturally arise in this intricate microstructure and place the proposed model in the context of phase transformations to high temperature transition aluminas.

36 MATERIALS SCIENCE↗
DOE OSTI2023

The effect of pore structure in ethane-selective metal-organic frameworks for ethylene purification

Separation of ethylene from ethane is of great importance in petrochemical industry. Adsorptive separation making use of porous materials such as metal-organic frameworks (MOFs) is an energy efficient technology compared to the conventional cryogenic distillation. However, preferential adsorption of ethylene over ethane represents the major current challenge. Here we demonstrate how pore structure can significantly alter the adsorption selectivity in two isoreticular Zr-MOFs. A simple ligand functionalization in UiO-66 yields smaller/better matching pore for ethane, and weaker interaction with ethylene, greatly enhancing ethane/ethylene selectivity in MOF-801. Ab initio calculations combined with in situ infrared spectroscopic analysis unveil the nature of host-guest interactions and differences in the overall binding energies. Polymer-grade ethylene (99.9%) can be produced directly from ethane/ethylene mixture in a single step. Having well-balanced adsorption capacity and selectivity, low-cost synthesis and easy scalability MOF-801 stands out as a promising ethane-selective adsorbent for one-step ethylene purification under ambient conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Zero-strain high entropy spinel oxide (FeNiCuCrMn) 3 O 4 @rGO as high-performance anode for sodium ion battery

High entropy spinel oxides (HEOs), as a new type of anode material with at least five transition metal elements and high theoretical specific capacity, have enormous promise in sodium ion batteries (SIBs). However, the influence of different transition metal elements on the structure and electrochemical performances of HEOs is not clear. Here, in this study, the addition of Mn to HEO can increase the surface oxygen vacancy (O v ) concentration, shorten the band gap, and enhance the electrochemical performance of (FeNiCuCrMn) 3 O 4 @rGO. XPS results show (FeNiCuCrMn)3O4@rGO has the highest surface O v concentration. Ab initio calculations demonstrate that the indirect band gap of (FeNiCuCrMn)3O4@rGO is 0.068 eV. In situ X-ray diffraction and synchrotron X-ray absorption spectroscopy reveal that (FeNiCuCrMn) 3 O 4 @rGO exhibits a zero-strain characteristic with an inconsequential unit cell volume change of 0.2 %. Specifically, (FeNiCuCrMn) 3 O 4 @rGO as anode materials for SIBs show good cyclic stability (with a capacity retention of 84 % at 500 mA g −1 after 3000 cycles). Overall, this work provides a new direction for optimizing transition metal elements in HEOs.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Theoretical and kinetic modeling study of H 2 S pyrolysis

Hydrogen sulfide pyrolysis was investigated theoretically and through chemical kinetic modeling. Reactions on the SHH potential energy surface, primarily S + H 2 (+Ar) ⇌ H 2 S (+ Ar) (R1) and S + H 2 ⇌ SH + H (R6b) were characterized by ab initio calculations. Results for k 1 were in good agreement with experiment, but deviated strongly below 2000 K from values previously used in modeling. Collider efficiencies for H 2 S, S 2 , and N 2 compared to Ar were calculated for R1. Hydrogen sulfide decomposition experiments reported in literature were re-examined in terms of an updated detailed chemical kinetic model. Concentration profiles for the atomic S at high temperature in shock tubes supported the present value of k 1 and served to constrain the rate constants for reaction of S with SH and H 2 S. To explain results from batch and flow reactors, conducted at high H 2 S concentrations in the 900–1400 K range, a very fast rate constant was required for HSS + H ⇌ SH + SH. Under dilute conditions, the gas-phase chemistry was too slow to compete and the decomposition of H 2 S was controlled by loss on the reactor surface.

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

Conversion of acenaphthalene to phenalene via methylation: A theoretical study

Ab initio calculations of the C 13 H 10 and C 13 H 9 potential energy surfaces related to the reaction of 1-acenaphthyl and methyl radicals and secondary isomerization of C 13 H 9 primary radical products have been performed at the chemically accurate G3(MP2,CC)//B3LYP/6-311G** level of theory to unravel the mechanism of conversion of acenaphthalene to phenalene or phenalenyl radical + H. The computed energetics and molecular parameters were utilized in Rice-Ramsperger-Kassel-Marcus Master Equation (RRKM-ME) calculations of reaction rate constants and relative product yields. The 1-acenaphthyl + CH 3 reaction is predicted to proceed by a fast radical-radical recombination mechanism and to predominantly produce collisionally stabilized 1-methylacenaphthalene B1 or a C 13 H 9 benzylic radical A1 + H with exothermicities of 106.8 and 25.2 kcal/mol, respectively. The A1 + H channel is preferable at higher temperatures, whereas the stabilization of 1-methylacenaphthalene is favored at higher pressures. The radical A1 can nearly irreversibly interconvert to phenalenyl radical AP via a 22.7 kcal/mol exothermic isomerization process involving formal insertion of the CH 2 group into a C-C bond of the five-member ring leading to the expansion of this ring to a six-member ring. The rate constants for the H addition reaction to phenalenyl radical to form phenalene and the reverse H loss from phenalene were also computed and the strength of the weakest C-H bond in the CH 2 group of phenalene is evaluated as 62.2 kcal/mol. The analysis of the reaction kinetics allowed us to deduce a mechanism for the conversion of 1-acenaphthyl radical to phenalene or phenalenyl radical + H via methylation involving the formation of A1 via a well-skipping channel or stabilization and dissociation of (or H abstraction from) B1 followed by isomerization of A1 to phenalenyl AP, which can add an H atom producing phenalene. As a result, rate constants for the significant elementary reactions are fitted to modified Arrhenius expressions and are proposed for kinetic modeling of the expansion of a five-member ring on a zigzag edge of PAH to a six-member ring by methylation.

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

Radiation-reduced diffusion of Nd in bcc Fe

Here, this study combines ab initio calculations and kinetic Monte Carlo simulations to investigate the bulk diffusion of Nd, an exemplar lanthanide element, in bcc Fe. Contrary to the conventional belief that irradiation accelerates diffusion processes, our study reveals that irradiation can unexpectedly reduce the diffusion of Nd due to the strong binding between substitutional Nd atoms and vacancies in bcc Fe. Combined with the high concentration of nonequilibrium vacancies, this promotes the abundant formation of Nd+3Va defect complexes containing a single Nd substitutional atom and three vacancies, accompanied by the depletion of Nd-vacancy pairs that mediate the diffusion of Nd.

Density functional theory (DFT)↗
DOE OSTI2022

express: Extensible, high-level workflows for swifter ab initio materials modeling

In this work, we introduce an open-source Julia project, express, an extensible, lightweight, high-throughput, high-level workflow framework that aims to automate ab initio calculations for the materials science community. express is shipped with well-tested workflow templates, including structure optimization, equation of state (EOS) fitting, phonon spectrum (lattice dynamics) calculation, and thermodynamic property calculation in the framework of the quasi-harmonic approximation (QHA). It is designed to be highly modularized so that its components can be reused across various occasions, and customized workflows can be built on top of that. Users can also track the status of workflows in real-time, and rerun failed jobs thanks to the data lineage feature express provides. Finally, two working examples, i.e., all workflows applied to lime and akimotoite, are also presented in the code and this paper.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Ultra-large-scale ab initio quantum chemical computation of bio-molecular systems: The case of spike protein of SARS-CoV-2 virus

The COVID-19 pandemic poses a severe threat to human health with an unprecedented social and economic disruption. Spike (S) glycoprotein of the SARS-CoV-2 virus is pivotal in understanding the virus anatomy, since it initiates the first contact with the ACE2 receptor in the human cell. We report results of ab initio computation of the spike protein, the largest ab initio quantum chemical computation to date on any bio-molecular system, using a divide and conquer strategy by focusing on individual structural domains. In this approach we divided the S-protein into seven structural domains: N-terminal domain (NTD), receptor binding domain (RBD), subdomain 1 (SD1), subdomain 2 (SD2), fusion peptide (FP), heptad repeat 1 with central helix (HR1-CH) and connector domain (CD). The entire Chain A has 14,488 atoms including the hydrogen atoms but excluding the amino acids with missing coordinates based on the PDB data (ID: 6VSB). The results include structural refinement, ab initio calculation of intramolecular bonding mechanism, 3- dimensional non-local inter-amino acid interaction with implications for the inter-domain interaction. Details of the electronic structure, interatomic bonding, partial charge distribution and the role played by hydrogen bond network are discussed. In the interaction among structural domains, we present new insights for crucial hinge-like movement and fusion process. Extension of such calculation to the interface between the S-protein binding domain and ACE2 receptor can provide a pathway for computational understanding of mutations and the design of therapeutic drugs to combat the COVID-19 pandemic.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Identification of Auger mechanisms responsible for low energy electron emission from graphene on copper using Auger-gamma coincidence spectroscopy

We have applied positron annihilation induced Auger-gamma coincidence spectroscopy to identify important mechanisms responsible for the emission of low energy electrons following the sudden creation of holes in bilayer graphene on copper substrate. The novel surface spectroscopic method measures the energy of the Doppler shifted annihilation gamma photon in coincidence with the Auger electron emitted following the relaxation of the hole created by the annihilation of a surface electron with a surface trapped positron. By extracting and theoretically modelling the annihilation gamma spectra coincident with low energy electrons, we associate majority of the intensity in the low energy (7 eV–25 eV) region of the Auger spectrum to electron emission following Auger decays of 2s holes in adsorbed oxygen and deep valence holes in graphene. We provide additional support to this conclusion by showing that most of the Auger electrons with energy less than ~ 25 eV are coincident with gamma photons with small Doppler shift (511 keV–512 keV) indicating that the primary electron whose annihilation resulted in low energy Auger electron emission had small momentum parallel to the gamma emission direction. On the other hand, majority of the Auger electrons measured in coincidence with the photons having maximum Doppler shift (515 keV–521 keV) were emitted following the decay of core holes (C 1s and O 1s). Furthermore, by selecting annihilation gamma photons coincident only with the O KVV (from surface adsorbed oxygen), C KVV (graphene), and Cu MVV (copper substrate) Auger electrons, we have also experimentally derived the energy spectrum of Doppler shifted gamma photons resulting from the annihilation of positrons with 1s electrons of O, 1s electrons of C and 3p electrons of Cu respectively. Model Doppler broadened gamma spectra produced using ab-initio calculations agree well with the experimentally derived line shapes. This demonstrates the ability of Auger-gamma coincidence method to experimentally resolve the Doppler broadened annihilation gamma spectrum from surfaces into its veiled electronic level constituents which has, heretofore, relied solely on theoretical analysis. Our results also demonstrate that the line shape of the Doppler broadened annihilation gamma peak reflect the chemical composition of the topmost atomic layer and can thus be used to characterize both external and inaccessible surfaces of porous materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Iron, magnesium, and titanium isotopic fractionations between garnet, ilmenite, f ayalite, biotite, and tourmaline: Results from NRIXS, ab initio , and study of mineral separates from the Moosilauke metapelite

Interpreting isotopic signatures documented in natural rocks requires knowledge of equilibrium isotopic fractionation factors. Here, we determine equilibrium Fe isotope fractionation factors between several common rock-forming minerals using a comparative approach involving three independent methods: (i) isotopic analyses of natural minerals from a metapelite from Mt. Moosilauke, New Hampshire, for which equilibration temperature and pressure are well constrained to be near the aluminosilicate triple point (T ≃ 500 °C, P ≃ 4 kbar), (ii) Nuclear Resonant Inelastic X-ray Scattering (NRIXS) measurements of Fe force constants of minerals, and (iii) Density Functional Theory (DFT) ab initio calculations of Fe force constants of minerals. The minerals studied for Fe isotopes include, in increasing order of their β-factors: garnet < ilmenite ≈ fayalite < biotite < tourmaline < muscovite ≈ plagioclase. Some of this ordering is affected by the presence of Fe 3+ in the minerals, which tends to form stiffer bonds and be associated with heavy Fe isotope enrichments relative to Fe 2+ . We are, however, able to assess the magnitude of the effect of the ratio Fe 3+ /ΣFe on equilibrium fractionation factors, notably on the ilmenite-hematite solid solution. Equilibrium Fe isotopic fractionation factors between garnet, ilmenite, biotite, tourmaline and fayalite are determined. We also report Mg and Ti isotopic compositions of selected Moosilauke minerals that allow us to better constrain the equilibrium fractionation factors for garnet-biotite-tourmaline (Mg isotopes) and biotite-ilmenite (Ti isotopes). We show how the newly determined equilibrium fractionation factors can be used to address diverse problems in Earth and planetary sciences, notably (i) Fe and Mg isotopic fractionation during anatexis, (ii) Fe isotopic fractionation in lunar ilmenite, and (iii) Ti isotopic fractionation during fluvial transport of minerals.

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

Contrasting magnetic properties of polymorphic TbPt 3

Here, in this work, we report the contrasting magnetic behaviour of a bimorphic intermetallic compound TbPt 3 that undergoes a structural phase transformation from cubic AuBe 5 -type structure to the cubic AuCu 3 -type structure on annealing at high temperature for a very short time. The magnetic ground state configuration of TbPt 3 changes from ferromagnetic (FM) in AuBe 5 -type crystalline phase to antiferromagnetic (AFM) in AuCu 3 -type phase. However, the AFM ground state of the AuCu 3 -type compound is unstable, as a moderate magnetic field can induce the FM character through metamagnetic transformation. To understand the said crystallographic as well as magnetic transformations, structural and magnetic ab-initio calculations have been carried out for both phases. From theoretical calculations we have proposed that although the AuBe 5 -type structure of TbPt 3 irreversibly transforms to the AuCu 3 -type TbPt 3 by a very short heat treatment, the application of moderate magnetic field on AuCu 3 -type structure may result in recovering the magnetic properties of AuBe 5 -type structure of TbPt 3 . As the two different polymorphs exist simultaneously at ambient condition, TbPt 3 thus provides us a great opportunity to study directly the effect of crystal structure on its physical properties and vice versa.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Thermoelectric performance in disordered Cu 2 ZnSnSe 4 nanostructures driven by ultra-low thermal conductivity

Here, high-energy reactive mechanical alloying (ball milling) was used to synthesize tetragonal and cubic polymorphs of Cu 2 ZnSnSe 4 . The ordered tetragonal (I-4) polymorph undergoes a phase transition above 400 K into a Cu-Zn disordered tetragonal (I-42 m) polymorph, while the cubic (F-43 m) polymorph with full cation disorder is stabilized at room temperature. Both polymorphs show ultra-low thermal conductivities, 0.42 W m -1 K -1 at 722 K and 0.21 W m -1 K -1 at 523 K for the disordered tetragonal and cubic phases respectively. The cubic polymorph has a higher zT in the low-temperature range, peaking at 0.26 (523 K), while the disordered tetragonal has a maximum zT of 0.46 at 712 K. The latter is the highest reported zT for stoichiometric Cu 2 ZnSnSe 4 , comparable to the best-performing doped materials in the literature. A combination of experimental results and ab-initio calculations point to a coupling between structural disorder and microstructure as the mechanism behind the reported performance.

36 MATERIALS SCIENCE↗
DOE OSTI2020

The correlation between N deficiency and the mechanical properties of the Ti 2 AlNy MAX phase

The role of X deficiency on the mechanical properties of MAX phases was studied by synthesizing Ti 2 AlN through powder metallurgy in stoichiometric and sub/extra-stoichiometric nitrogen compositions. XRD analyses and ab initio calculations indicate that nitrogen vacancies result in a lattice contraction predominantly along the c-axis. The elastic moduli and intrinsic hardness of substoichiometric Ti 2 AlN 0.9 measured from nanoindentation tests are shown to be slightly smaller than that of Ti 2 AlN. The key mechanical indexes, bulk (B), shear (G) and Young’s (E) moduli as well as the hardness variation are calculated in density functional theory, and show different responses depending on the concentration of N vacancies. Overall, this joint experimental and theoretical study provides a full understanding of the energetics, chemical bonding, electronic structure, and mechanics of the N deficient MAX phases which would increase the application of nitride ceramics.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Phase-field modeling of dislocation–interstitial interactions

The mechanical behavior of body-centered cubic (BCC) materials can be dramatically affected by the presence of interstitial solute atoms. Here we present a new phase-field dislocation dynamics formulation to include the diffusion of interstitials. Short-range interactions are accounted for by a concentration-dependent lattice energy, and long-range interactions are accounted for by modifications to the elastic energy. The interstitial diffusion law introduces gradients that require methods for minimizing Gibbs oscillations, which is done via a modified Green’s function. The formulation is general to any solute-solvent system and is applied here to Nb-O as a model system, whose interstitial parameters are obtained from ab initio calculations. The effect of O on the core structures of Nb edge and screw dislocations is calculated. The diffusion of O to form interstitial atmospheres around dislocation cores is simulated, as well as the critical stresses required for dislocations to break away or cross slip from these atmospheres. Future applications of the method to simulate complex interstitial embrittlement mechanisms are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Role of chemical disorder on radiation-induced defect production and damage evolution in NiFeCoCr

Understanding chemical disorder in many concentrated solid solution alloys (CSAs) at the levels of electrons and atoms has attracted increasing attention as a path forward to reveal and identify underlying mechanisms for extraordinary mechanical properties and improved radiation tolerance. Single-phase NiFeCoCr CSA is a common base for many high-entropy alloys (HEAs) that have shown improved mechanical strength and radiation tolerance. In this study, defect production and damage evolution in NiFeCoCr under ion irradiation at room temperature to dose over 20 dpa are determined using ion channeling technique along both <100> and <110> directions utilizing multiple probing beam energies. The results obtained from the multi-axial and multi-energy channeling analysis are compared with those previously obtained for Ni crystals irradiated under similar conditions. The influence of chemical complexity on defect production and clustering at early-stage under room temperature irradiation up to dose of 1 dpa is discussed based on positron annihilation spectroscopy results. Defect structure evaluation in Ni and NiFeCoCr is also discussed based on transmission electron microscopy results over a prolonged irradiation at both room and elevated temperatures. Compared with chemically complex NiFeCoCr, larger dislocation loops thus less lattice strain are expected to form in pure Ni. Moreover, the role of chemical disorder in this CSA is also investigated based on ab initio calculations using large supercells. Finally, to understand the impact of chemical complexity effect on defect structure evolution, this integrated research effort attempts to link the relatively large charge redistribution due to difference in valence electron counts resulting from alloying different 3d transition metal elements, moderate lattice distortion arising from similar adaptable atomic size, and notable suppressed or delayed damage evolution in NiFeCoCr.

36 MATERIALS SCIENCE↗