Search NASA⌕ Search

SEARCH · Search NASA

Results for “Formation energy”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1980

Meteor Crater: Energy of formation - Implications of centrifuge scaling

Recent work on explosive cratering has demonstrated the utility of performing subscale experiments on a geotechnic centrifuge to develop scaling rules for very large energy events. The present investigation is concerned with an extension of this technique to impact cratering. Experiments have been performed using a projectile gun mounted directly on the centrifuge rotor to launch projectiles into a suitable soil container undergoing centripetal accelerations in excess of 500 G. The pump tube of a two-stage light-gas gun was used to attain impact velocities of approximately 2 km/sec. The results of the experiments indicate that the energy of formation of any large impact crater depends upon the impact velocity. This dependence, shown for the case of Meteor Crater, is consistent with analogous results for the specific energy dependence of explosives and is expected to persist to impact velocities in excess of 25 km/sec.

Schmidt, R. M.↗
DOE OSTI2023

Development of an interatomic potential for the W–Ta system

A physics-inspired, data-driven interatomic potential framework for multi-element system is presented. This potential is based on the generalization of the embedded atom method potential and systematically incorporates two-, three- and many-body effects. Different atomic environments are described by atom-centered Gaussian basis sets which provides sufficient flexibility to capture diverse atomic environments and allows for easy optimization of the free parameters. An interatomic potential model for the tungsten–tantalum (W–Ta) system is developed using this framework by training on data from ab initio density functional theory (DFT) calculations. The potential is thoroughly tested at various compositions by comparing bulk and defect properties from experimental and DFT data. It is shown that the potential predicts the elastic constants, defect properties, such as vacancy and interstitial formation energies, core structures of dislocations, and melting points for both pure and tantalum and tungsten with an accuracy comparable to other single element machine learning based potentials. In conclusion, the potential model is used to investigate the formation energies of ordered alloys and vacancy formation energies for chemically random W–Ta alloys.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

First-principles calculations of metal surfaces. I. Slab-consistent bulk reference for convergent surface properties

The first-principles computation of the surfaces of metals is typically accomplished through slab calculations of finite thickness. The extraction of a convergent surface formation energy from slab calculations is dependent upon defining an appropriate bulk reference energy. I describe a method for an independently computed, slab-consistent bulk reference that leads to convergent surface formation energies from slab calculations that also provides realistic uncertainties for the magnitude of unavoidable nonlinear divergence in the surface formation energy with slab thickness. In conclusion, the accuracy is demonstrated on relaxed, unreconstructed low-index aluminum surfaces with slabs with up to 35 layers.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Analysis of defects in In 2 O 3 :H synthesized in presence of water vapor and hydrogen gas mixture

Using hybrid functional-based density functional theory calculations, we analyze the structure and kinetics of defects formed in two competing synthesis routes to prepare hydrogen-doped In 2 O 3 films, using a hydrogen and oxygen gas mixture vs synthesis in the presence of water vapor. For both of these synthesis routes, we find that H + is the dominant defect species: when the Fermi level is close to the conduction band, H + has a lower formation energy than other intrinsic or extrinsic defects. Our results also suggest that water molecules spontaneously split into H + (which occupies octahedral voids) and OH – interstitials (which occupies vacant oxygen lattice sites or oxygen vacancies). From the analysis of the binding energies between these different defects, we conclude that these defects do not cluster and are most likely to stay spatially distributed throughout the films. In addition, the sum of formation energies of an oxygen (i.e., O$^{2–}_i$) and a H + interstitial is close to the formation energy of a OH – interstitial, meaning that water molecules are completely split into 2H + and O$^{2–}_i$ at the synthesis conditions. Further, in the presence of H 2 + O 2 gas mixture, oxygen interstitials occupy oxygen vacancies while hydrogen interstitials occupy vacant oxygen lattice sites and form bonds with lattice oxygens. Our analysis of the defect equilibria suggests that the hydrogen content in films synthesized in the presence of water vapor is higher than films synthesized in the presence of a hydrogen gas mixture. At high dopant concentrations, a hydrogen bond network is formed in the system and this leads to large distortions in the lattice.

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

Free energies of formation of WC and WzC and the thermodynamic properties of carbon in solid tungsten

The activity of carbon in the two-phase regions - W + WC and W + W2C was obtained from the carbon content of iron rods equilibrated with mixtures of metal plus carbide powders. From this activity data the standard free energies of formation of WC and W2C were calculated. The temperature of the invariant reaction W2C = W + WC was fixed at 1570 + or - 5K. Using available solubility data for C in solid W, the partial molar free energy of C in the dilute solid solution was also calculated. The heat of solution of C in W, and the excess entropy for the interstitial solid solution, were computed, assuming that the carbon atoms reside in the octahedral interstices of bcc W.

Gupta, D. K.↗
NASA NTRS1977

Energies of formation for ejecta blankets of giant impacts

Ejection energy-diameter scaling relations for impact craters are studied by assuming the energy is equal to a constant times the diameter raised to a (different) constant power. These scaling relations are derived from least-squares fits to calculated ejection energies over a range of crater diameters and vary according to ejection angle, ejecta blanket distribution, planetary radius, and the depth-diameter ratio. Extremely broad changes in ejection angle and ejecta blanket distribution result in relatively small changes in the total energy. The way in which the exponent in the scaling relation is dependent on the planetary radius of curvature as well as on the exponent of the power-law depth-diameter scaling relation is described. Two ejection energy-diameter scaling relations, one for simple craters with depth-diameter ratios of approximately 0.2 and the other for large complex craters, are presented. The scaling relations explicitly include target density, surface gravity, and planetary radius.

Croft, S. K.↗
DOE OSTI2023

Machine learning guided discovery of superconducting calcium borocarbides

Pursuit of superconductivity in light-element systems at ambient pressure is of great experimental and theoretical interest. In this work, we combine a machine learning (ML) method with first-principles calculations to efficiently search for the energetically favorable ternary Ca-B-C compounds. Three layered borocarbides (stable CaBC 5 and metastable Ca 2 BC 11 and CaB 3 C 3 ) are predicted to be phonon-mediated superconductors at ambient pressure. The stable CaBC 5 and the low-energy metastable Ca 2 BC 11 (with formation energy only 9.5 meV/atom above the convex hull) have a superconducting T c of 5.2 and 8.9 K, respectively. While the hexagonal CaB 3 C 3 possesses a T c of 26.1 K, it is metastable with formation energy of 153 meV/atom above the convex hull. The ML-guided approach opens up a way for greatly accelerating the discovery of new high-T c superconductors.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Exchange-correlation functional challenges in modeling quaternary chalcogenides

The development of next-generation quaternary chalcogenides, such as Cu2ZnSnS4 (CZTS) and Cu2ZnGeS4 (CZGS), for solar energy and thermoelectric applications hinges upon both careful experimentation and accurate quantum mechanical modeling. To address the latter, many have turned to density functional theory (DFT), which offers several choices for the approximate treatment of electron exchange and correlation (XC). Popular XC functionals include the Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation (GGA) and the recently developed strongly constrained and appropriately normed (SCAN) meta-GGA. Extensions of DFT functionals, such as adding a Hubbard U correction and introducing a fraction of the Fock exchange (hybrid functionals), have been used widely to model systems containing 3d metal ions. However, no studies yet have compared comprehensively PBE(+U) and SCAN(+U) in the quality of their predictions of the bulk and defect thermodynamics of quaternary chalcogenides, which play a critical role in device fabrication and performance. Hence, here we calculate the (i) 0 K formation energies of bulk Ge compounds and (ii) neutral defect formation energies including charge-balanced (e.g., CuZn + ZnCu) and charge-imbalanced (e.g., CuSn) combinations of antisites and vacancies in CZTS and CZGS using the PBE, PBE +U, SCAN, SCAN +U, and the hybrid Heyd-Scuseria-Ernzerhof XC frameworks. We find that the formation energies of charge-imbalanced defects are more sensitive to the choice of the XC functional than those of charge-balanced defects, which can be explained by the differences in the extent of penalization of defect-generated delocalized electrons/holes by PBE, PBE +U, SCAN, and SCAN +U. Additionally, our results show that SCAN systematically underbinds Ge-containing compounds, thus highlighting the need for even further improvement of XC functionals. Based on our findings, we recommend the use of SCAN for modeling quaternary chalcogenides because its errors are systematic, and it has the firmest theoretical underpinning. Our work provides guidance for future modeling of quaternary chalcogenides.

14 SOLAR ENERGY↗
NASA NTRS1993

The Manson impact crater: Estimation of the energy of formation, possible size of the impacting asteroid or comet, and ejecta volume and mass

A research program on the Manson impact structure has substantially improved our knowledge of the detailed features of this eroded crater. As part of our structural studies, we have derived a value of 21 km for the diameter of the final transient cavity formed during crater excavation. With this information, we can estimate the energy of formation of the Manson crater and the possible size of the impacting asteroid or comet. In addition, we have estimated the near- and far-field ejecta volumes and masses.

Roddy, D. J.↗
DOE OSTI2022

Role of Sr doping and external strain on relieving bottleneck of oxygen diffusion in La 2- x Sr x CuO 4-δ

In many complex oxides, the oxygen vacancy formation is a promising route to modify the material properties such as a superconductivity and an oxygen diffusivity. Cation substitutions and external strain have been utilized to control the concentration and diffusion of oxygen vacancies, but the mechanisms behind the controls are not fully understood. Using first-principles calculations, we find how Sr doping and external strain greatly enhances the diffusivity of oxygen vacancies in La 2-x Sr x CuO 4-δ (LSCO) in the atomic level. In hole-doped case (2x > δ), the formation energy of an apical vacancy in the LaO layer is larger than its equatorial counterpart by 0.2 eV that the bottleneck of diffusion process is for oxygen vacancies to escape equatorial sites. Such an energy difference can be reduced and even reversed by either small strain (< 1.5%) or short-range attraction between Sr and oxygen vacancy, and in turn, the oxygen diffusivity is greatly enhanced. For fully compensated hole case (2x ≦ δ), the formation energy of an apical vacancy becomes too high that most oxygen vacancies cannot move but would be trapped at equatorial sites. From our electronic structure analysis, we found that the contrasting change in the formation energy by Sr doping and external strain is originated from the different localization natures of electron carrier from both types of oxygen vacancies.

25 ENERGY STORAGE↗
DOE OSTI2025

A first principles study on the adsorbate-adsorbate interactions on the CdTe(111) surface with Cd, Te, Zn, and Se adatoms

The study of adsorbate-adsorbate interactions is essential to understanding early crystal growth dynamics. Here, we employ planewave density functional theory to study the binary adatom pair interactions between Cd-Cd, Te-Te, Zn-Zn, Se-Se, Cd-Te, Cd-Se, Cd-Zn, Te-Se, Te-Zn, and Se-Zn adatom pairs on two CdTe(111) surfaces. An analysis of the interaction energies between binary adatom pairs suggests repulsive interactions are common regardless of the relative distance between adatoms. For the CdTe(111)A surface, attractive interactions occur between neighboring chalcogen (i.e., Te and Se) and Group 12 (i.e., Cd and Zn) adatom pairs. For the CdTe(111)B surface, attractive interactions occur between neighboring Group 12 adatoms forming a surface dimer configuration. Furthermore, the formation energy of an adatom pair is decomposed in terms of the electronic, elastic, and adatom binding contributions. For smaller interatomic distances between the adatoms, the formation energy is primarily a function of the electronic interactions, with null contributions from the elastic and adatom binding interactions for Group 12-containing pairs. Because of the less favorable electronic interactions for larger interatomic distances between the adatoms, the formation energies are typically more positive. Lastly, neighboring adatoms significantly increase the barriers of migration on the CdTe(111)A surface relative to unary adatoms for the top-to-fcc and fcc-to-fcc sites, while the migration barriers on the CdTe(111)B surface only increases for the fcc-to-fcc migration of chalcogen species. From this analysis, we illustrate the role of adatom interactions during the early stages of the surface nucleation processes on CdTe(111) thin films.

CdTe↗
NASA NTRS1986

Analysis of boron carbides' electronic structure

The electronic properties of boron-rich icosahedral clusters were studied as a means of understanding the electronic structure of the icosahedral borides such as boron carbide. A lower bound was estimated on bipolaron formation energies in B12 and B11C icosahedra, and the associated distortions. While the magnitude of the distortion associated with bipolaron formation is similar in both cases, the calculated formation energies differ greatly, formation being much more favorable on B11C icosahedra. The stable positions of a divalent atom relative to an icosahedral borane was also investigated, with the result that a stable energy minimum was found when the atom is at the center of the borane, internal to the B12 cage. If incorporation of dopant atoms into B12 cages in icosahedral boride solids is feasible, novel materials might result. In addition, the normal modes of a B12H12 cluster, of the C2B10 cage in para-carborane, and of a B12 icosahedron of reduced (D sub 3d) symmetry, such as is found in the icosahedral borides, were calculated. The nature of these vibrational modes will be important in determining, for instance, the character of the electron-lattice coupling in the borides, and in analyzing the lattice contribution to the thermal conductivity.

Howard, Iris A.↗
DOE OSTI2023

Graph neural networks predict energetic and mechanical properties for models of solid solution metal alloy phases

Here, we developed a PyTorch-based architecture called HydraGNN that implements graph convolutional neural networks (GCNNs) to predict the formation energy and the bulk modulus for models of solid solution alloys for various atomic crystal structures and relaxed volumes. We trained the GCNN surrogate model on a dataset for nickel–niobium (NiNb) generated by the embedded atom model (EAM) empirical interatomic potential for demonstration purposes. The dataset was generated by calculating the formation energy and the bulk modulus as a prototypical elastic property for optimized geometries starting from initial body-centered cubic (BCC), face-centered cubic (FCC), and hexagonal compact packed (HCP) crystal structures, with configurations spanning the possible compositional range for each of the three types of initial crystal structures. Numerical results show that the GCNN model effectively predicts both the formation energy and the bulk modulus as function of the optimized crystal structure, relaxed volume, and configurational entropy of the model structures for solid solution alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A comprehensive first-principles study of the effects of the exchange-correlation functional and magnetism on defect and diffusion properties of the CoCrNi medium-entropy alloy

The present work is a novel, systematic study of the effect of density functional theory input parameters on the vacancy formation energy (VFE), migration barrier for diffusion, and electronic structure for each element in the CoCrNi medium-entropy alloy (MEA). In particular, the novelties include: (1) calculating the aforementioned properties of Co, Cr, or Ni, in the CoCrNi MEA using magnetic and non-magnetic states, and two versions of the generalized gradient approximation: Perdew, Burke, and Ernzerhof (PBE) and the PBE version for solids (PBEsol), and (2) a detailed comparison of 0 K activation energy to experimental creep activation energies. First-principles calculations at 0 K are performed using the Vienna ab-initio simulation package. Special quasirandom structures (SQS) and Widom-type substitution are employed. For each element, Co, Cr, or Ni, non-magnetic calculations result in a higher VFE and larger range of calculated values for the configurations studied. The averaged migration barrier is the highest for Co in the CoCrNi for three of four sets of calculation parameters in the configurations studied. Finally, the results indicate that the average 0 K activation energy for diffusion makes up 70–80% of the experimental creep activation energy, depending on the exchange-correlation functional employed.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Reversible Methanolation of Metal Halide Perovskites

The low formation energies inherent in metal halide perovskites enable the structure to be easily broken and remade with little energy input. In this work, we leverage low formation energy to demonstrate 0D/3D structural transformations induced by reversible intercalation of methanol to yield dramatic control of visible light absorption. We identify a methanolated structure that features sheets of 0D isolated [PbI6]4- octahedra separated by MAI and methanol. Methanol and water reversibly displace each other in the 0D complex by controlling the chemical potential of the system via Le Chatelier's principle. The weaker H-bonding of methanol compared to water effectively reduces the complex dissociation temperature from 70 to 50 degrees C, rendering the methanol complex more desirable for an array of next-generation applications spanning low-power-consumption memory to switchable photovoltaics.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Node support input to HydroGEN STCH project quarterly report

n this quarter, we further refined our water splitting (WS) reaction calculations on the surface of SrTi 0.5 Mn 0.5 O 3 . Among various surface configurations for WS imposed by surface orientations, terminations, and neighboring chemical species, we hypothesized that (1) the WS reaction energy should be exothermic and (2) the surface oxygen vacancy formation energy should be in a reasonable range compared to the bulk oxygen vacancy formation energy to be accessible during the reduction and WS reaction conditions.

08 HYDROGEN↗
DOE OSTI2022

Anomalous stability of non-van der Waals bonded B4C nanosheets through surface reconstruction

Boron carbide (B4C) has been well studied both theoretically and experimentally in its bulk form due to its exceptional hardness and use as a high-temperature thermoelectric. However, the properties of its two-dimensional nanosheets are not well established. In this paper, using van der Waals-corrected density-functional theory simulations, we show that bulk B4C can be cleaved along different directions to form B4C nanosheets with low formation energies. We find that there is minimal dependence of formation energies on cleavage planes and surface terminations, even though the bulk is not van der Waals layered. This anomalous stability of B4C nanosheets is found to be a result of surface reconstructions that are unique to B-rich systems. While the density of states of the bulk B4C indicate that it is a semiconductor, the B4C nanosheets are found to be predominantly metallic. We attribute this metallic behavior to a redistribution of charges on the surface bonds of the films. The Seebeck coefficients of the B4C films remain comparable to those of the bulk and are nearly constant as a function of temperature. Our results provide guidance for experimental synthesis efforts and future application of B4C nanosheets in nanoelectronic and thermoelectric applications.

Physics↗
DOE OSTI2021

First-principles investigation of uranium mononitride (UN): Effect of magnetic ordering, spin-orbit interactions and exchange correlation functional

Uranium mononitride (UN) is a promising nuclear fuel that combines the advantageous properties of readily used UO 2 and uranium alloys. Various properties of UN have been previously studied using different density functional theory (DFT) methodologies; however, there are still inconsistencies when it comes to the dynamical stability and defect properties of UN. We address these inconsistencies by studying the UN phonons and defect properties using DFT calculations employing two generalized gradient approximation (GGA) exchange-correlation functionals: PBE and AM05, with and without an added on-site Coulomb repulsion term (+U). Furthermore, we investigate the importance of spin-orbit coupling (SOC) when calculating the properties of UN. We use the different methodologies to determine the preference of UN to have antiferromagnetic (AFM) ordering, as seen in experiments, or ferromagnetic (FM) ordering of the uranium spins. We compare the crystallographic properties, density of states, the DFT X-ray photoelectron spectra and phonon dispersions calculated using the different methodologies. We demonstrate that GGA + U reproduces the AFM ordering in UN, but the crystal structure is dynamically unstable. We also show that magnetic ordering is important in finding the lowest energy defective structure, and that SOC has a distinct influence on the energy of the different uranium interstitial defects. Lastly, we discuss the point defect formation energies under U-rich and N-rich conditions, and the stoichiometric formation energies calculated with the different methodologies, providing insight into the observed tendency for forming hypostoichiometric UN.

36 MATERIALS SCIENCE↗