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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Effects of Aluminum and Molybdenum on the Phase Stability of Iron-Chromium Alloys: A First-Principles Study

The interaction between solute atoms is critical to the thermodynamic behavior of Fe-Cr alloys, but the effects of non-dilute Al and Mo on the Fe-Cr phase stability and vacancy formation energy are not clearly understood. In this study, density functional theory, cluster expansion, and Monte Carlo simulation are used to predict the effects of ternary solute elements on the thermodynamic properties in multicomponent Fe-Cr alloys. The machine learning regression approach is applied to train and construct energy models that accurately describe the complex chemical interactions. The computational outcomes include the prediction of the partial ternary phase diagram, mixing enthalpy, and vacancy formation energy for different compositions. The phase boundary calculation predicts a pronounced change of Cr solubility in bcc Fe by the addition of Al and the rejection of Al atoms from a’ precipitates. The mixing enthalpy calculation shows that Mo may also reduce the Cr solubility in bcc Fe. Additionally, the simulations show that non-dilute Cr decreases the vacancy formation energy in bcc Fe, while adding Al results in a less significant effect. Finally, the results demonstrate important applications of using machine learning energy models to study model or commercial alloys with multicomponent solute species and point defects.

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

Helium impurities and interactions in lithium

Here we investigate helium interactions with lithium using density functional theory. Like other body-centered cubic (bcc) metals, the lowest-energy site for interstitial helium is a tetrahedral site. However, helium in lithium shows a higher substitutional formation energy than either the tetrahedral or octahedral interstitial formation energies, which is unique. The calculated migration energies of helium in lithium are also very low ( ~ 1 meV), an order of magnitude lower than those in other bcc metals. We find an increase in the binding energy to a vacancy as the number of helium atoms bound to it increases. The extremely low migration energies suggest that helium transport in lithium will be very fast.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Effect of Linker Structure and Functionalization on Secondary Gas Formation in Metal–Organic Frameworks

Rare-earth terephthalic acid (BDC)-based metal–organic frameworks (MOFs) are promising candidate materials for acid gas separation and adsorption from flue gas streams. However, previous simulations have shown that acid gases (H 2 O, NO 2 , and SO 2 ) react with the hydroxyl on the BDC linkers to form protonated acid gases as a potential degradation mechanism. Herein, gas-phase computational approaches were used to identify the formation energies of these secondary protonated acid gases across multiple BDC linker molecules. Formation energies for secondary protonated acid gases were evaluated using both density functional theory (DFT) and correlated wave function methods for varying BDC–gas reaction mechanisms. Upon validation of DFT to reproduce wave function calculation results, rotated conformational linkers and chemically functionalized BDC linkers with -OH, -NH 2 , and -SH were investigated. Furthermore, the calculations show that the rotational conformation affects the molecule stability. Double-functionalized BDC linkers, where two functional groups are substituted onto BDC, showed varied reaction energies depending on whether the functional groups donate or withdraw electrons from the aromatic system. Based on these results, BDC linker design must balance adsorption performance with degradation via linker dehydrogenation for the design of stable MOFs for acid gas separations.

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

GeSn-on-GaAs with photoconductive carrier lifetime >400 ns: role of substrate orientation and atomistic simulation

Group IV GeSn quantum material finds application in electronics and silicon-compatible photonics. Synthesizing these materials with low defect density and high carrier lifetime is a potential challenge due to lattice mismatch induced defects and tin segregation at higher growth temperature. Recent advancements in the growth of these GeSn materials on Si, Ge, GaAs, and with substrate orientations, demonstrated different properties using epitaxial and chemical deposition methods. This article addresses the effect of GaAs substrate orientation and misorientation on the materials’ properties and carrier lifetimes in epitaxial Ge 0.94 Sn 0.06 layers. With starting GaAs substrates of (100)/2°, (100)/6°, (110) and (111)A orientations, Ge 0.94 Sn 0.06 epitaxial layers were grown with an intermediate Ge buffer layer by molecular beam epitaxy and analyzed by several analytical tools. X-ray analysis displayed good crystalline quality, and Raman spectroscopy measurements showed blue shifts in phonon wavenumber due to biaxial compressive strain in Ge 0.94 Sn 0.06 epilayers. Cross-sectional transmission electron microscopy analysis confirmed the defect-free heterointerface of Ge 0.94 Sn 0.06 /Ge/GaAs heterostructure. Minority carrier lifetimes of the unintentionally doped n-type Ge 0.94 Sn 0.06 epilayers displayed photoconductive carrier lifetimes of >400 ns on (100)/6°, 319 ns on (100)/2°, and 434 ns on (110) GaAs substrate at 1500 nm excitation wavelength. On the other hand, Ge 0.94 Sn 0.06 layer showed poor carrier lifetime on (111)A GaAs substrate. The observed differences in carrier lifetimes were correlated with the formation energy of the Ge on (100)/6° and (100)/2° GaAs heterointerface using Stillinger–Weber interatomic potential model-based atomistic simulation with different heterointerfacial bonding by Synopsys QuantumATK tool. Total energy computation of 6280-atom Ge/GaAs supercell on (100)/6° leads to lower formation energy than (100)/2°, making it more thermodynamically stable. Hence, the growth of the GeSn/III–V material system using misoriented (100) substrates that are more thermodynamically stable will enhance the performances of optoelectronic devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Efficient Modeling of Structural, Electronic, and Optical Properties of Silver and Gold Metal Nanoclusters and Alloys Using Optimized SCC-DFTB Parameters

Computation of optical properties using conventional time-dependent density functional theory (TD-DFT) is time-consuming and memory-intensive. In this study, we investigate the accuracy and efficiency of the density functional tight binding (DFTB) framework with newly optimized Slater–Koster (SK) parameters for modeling the structural, electronic properties, and absorption spectra of silver and gold nanoclusters and their alloys. Our investigation of the ground state (GS) properties demonstrates that the newly developed GS-SK parameters enable DFTB to closely approximate DFT-calculated bond lengths for octahedron, tetrahedron, icosahedra, and truncated octahedron with sizes Ag n /Au n (n = 19, 20, 38, 55), nanoclusters and Ag 20 /Au 20 nanoalloys, with a maximum deviation of approximately 0.15 Å. Formation energy results indicate that the GS-SK parameters can closely estimate changes in formation energies with alloy composition, and the comparison of electronic structures for Ag 20 , Au 20 , and AgAu alloy nanoclusters using the DFTB approximation reveals good agreement in the projected density of states (DOS) profiles and energy levels. A second set of SK parameters, ES-SK, has been developed to describe excited state (ES) properties, including the absorption spectra of silver octahedron Ag 19 , tetrahedral Ag n (n = 20, 56, 84), truncated octahedron Ag 38 , and icosahedra Ag 55 closed-shell clusters and their gold and alloy counterparts over a broad range of alloy compositions. This parametrization uses TD-DFTB calculations and fine-tunes the d and p eigenvalues by comparing them to reference absorption spectra from first-principles TD-DFT. This enables the generation of absorption spectra that closely match the reference spectra when plasmon excitation is dominant, as demonstrated by studying the plasmonic properties of icosahedral Ag n and Au n (n = 309 and 561) nanoparticles. This includes the rapid loss in plasmon quality when Au partially replaces Ag in alloy clusters. Furthermore, these results provide a foundation for addressing computational bottlenecks in plasmonics and with new prospects for applications in the quantum plasmonics for bimetallic alloys.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗