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 91 records · Page 5

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↗

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↗

Vacancy formation and extraction energies in semiconductor compounds and alloys

A model for calculating the extraction energies and vacancy-formation energies in diamond-cubic and zinc-blende semiconductors is developed on the basis of Harrison's (1980, 1983) tight-binding theory. The extraction energies provide a reference from which other final states of the removed atoms can be calculated. The results of calculations show that, in a given compound, the calculated extraction energies are larger for the anion than for the cation, with the difference between the cation and the anion being larger in the II-VI than in the III-V compounds. This is in agreement with experimental results.

Berding, M. A.↗

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↗

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↗

First-principles study of intrinsic point defects and Xe impurities in uranium monocarbide

Based on density functional theory (DFT) calculations, we perform an extensive investigation of intrinsic point defects and Xe impurities in uranium monocarbide (UC). The DFT calculations involve both the conventional generalized gradient approximation (GGA) and the GGA+U approach with the Hubbard parametric term (U), using up to 5×5×5 supercells. GGA calculations for the formation energy of intrinsic defects demonstrate the significant effect of using larger supercells than in previous studies. Results confirm that the ⟨111⟩ and ⟨100⟩ dumbbell interstitials are the most stable interstitial configurations for U and C, respectively. The interstitial mechanisms are favored for self-diffusion of both uranium and carbon and diffusion of Xe under equilibrium conditions. Calculations also reveal that the Xe substitutional defect at the C lattice site tends to adopt an off-site configuration, which can be interpreted as a Xe interstitial–C vacancy complex. We also utilize GGA+U to assess the impact of effective U parameter (Ueff) on the results. Moreover, we introduce a method to estimate the carbon chemical potential by fitting the phase diagram composition data and propose a selection Ueff=1.25 eV based on the experimental Xe diffusion activation energy. With this approach, GGA+U calculations reproduce the available experimental data for the formation energy of the carbon Frenkel pair and can explain the stepwise recovery of intrinsic properties and burst Xe release behavior in UC observed in annealing experiments.

Huang, Gui-Yang (ORCID:0000000301447376)↗

Influence of point defects on the electronic and topological properties of monolayer WTe 2

In some topological insulators such as monolayers of WTe 2 or graphene with spin-orbit coupling, band inversion originates from chemical bonding and space group symmetry, in contrast to materials such as Bi 2 Se 3 , where the band inversion derives from relativistic effects in the atoms. In the former, band inversion is susceptible to changes of the chemical environment, e.g., by defects, while the latter are less affected by defects due to the larger energy scale associated with atomic relativistic effects. Motivated by recent experiments, we study the effect of Te vacancies and Te adatoms on the electronic properties of WTe 2 . We find that the Te vacancies have a formation energy of 2.21 eV, while the formation energy of the Te adatoms is much lower with 0.72 eV. Here, the vacancies strongly influence the band structure and we present evidence that band inversion is already reversed at the nominal composition of WTe 1.97 . In contrast, we show that the adatoms do not change the electronic structure in the vicinity of the Fermi level and thus the topological properties.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

First principles investigation of dopants and defect complexes in CdSe$_x$Te$_{1-x}$

Se alloying is a common approach to improve the performance of CdTe solar cells by tuning the bandgap, defect levels, and carrier density. A fundamental understanding of these improvements, specifically the effect of Se alloying on the behavior of defects and dopants in CdTe, remains unclear. Here, in this work, we present a density functional theory (DFT) study of point defect energetics in CdTe and CdSe x Te 1-x with x = 0.25, leading to a comparison of how native defects, dopants (As and Cu), impurities (Cl and O), and related defect complexes behave in CdTe vs CdSe x Te 1-x . Our calculations, performed by combining semi-local and nonlocal hybrid functionals, show a general lowering of the formation energies of native defects as well as substitutional defects formed by As and Cl upon Se addition. For successful p-type doping with As, destabilizing Cl-based defects in the CdSeTe lattice would be essential. We find evidence for some low-energy defect complexes of As, Cl, and O in CdSe 0.25 Te 0.75 . The computed defect formation energies further enable estimates of temperature-dependent defect concentrations and self-consistent Fermi levels. A comparison of defect energetics with the energies of impurity phases reveals that As, Cu, Cl, and O overwhelmingly prefer being segregated to unwanted As 2 O 5 , AsCl 3 , Cd 2 AsCl 2 , and CuO x phases rather than remain at defect sites, but such segregation is less likely to happen in CdSe 0.25 Te 0.75 than in CdTe. Overall, our work presents a list of likely defects and complexes in CdTe and Se-incorporated CdTe, paving the way to explain and mitigate limited dopant activation in experimental observations.

CdTe↗

Machine Learning-Guided Discovery of Ternary Compounds Containing La, P, and Group 14 Elements

In this work, we integrate a deep machine learning (ML) method with first-principles calculations to efficiently search for the energetically favorable ternary compounds. Using La–Si–P as a prototype system, we demonstrate that ML-guided first-principles calculations can efficiently explore crystal structures and their relative energetic stabilities, thus greatly accelerate the pace of material discovery. A number of new La–Si–P ternary compounds with formation energies less than 30 meV/atom above the known ternary convex hull are discovered. Among them, the formation energies of La 5 SiP 3 and La 2 SiP phases are only 2 and 10 meV/atom, respectively, above the convex hull. These two compounds are dynamically stable with no imaginary phonon modes. Moreover, by replacing Si with heavier-group 14 elements in the eight lowest-energy La–Si–P structures from our ML-guided predictions, a number of low-energy La–X–P phases (X = Ge, Sn, Pb) are predicted.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modeling of fission gas diffusion and release for Gd 2 O 3 doped UO 2

Uranium dioxide (UO 2 ) is the primary nuclear fuel in light water reactors, and its excess neutronic reactivity can be controlled by adding burnable absorbers, such as Gd 2 O 3 . This burnable absorber has a large neutron absorption cross-section, lowering the high reactivity of the reactor's initial fuel load. However, there needs to be more understanding of how added Gd 2 O 3 influences the properties of UO 2 under irradiation. To understand the behavior of defects and fission gas in the UO 2 /Gd 2 O 3 system under irradiation, we use cluster dynamics modeling supported by density functional theory calculations. First, we calculate the formation energies of Gd point and cluster defects, and evaluate the temperature-dependent defect concentrations using the defect formation energies and entropies. We show that Gd is soluble in UO 2 , introducing a negative charge in the system. Using this information, we adapted the cluster dynamics code Centipede to model the influence of Gd on U self-diffusion and Xe diffusion in UO 2 with 10 wt% Gd 2 O 3 . Also, we analyzed the Xe diffusion as a function of Gd 2 O 3 concentration, showing that the Xe diffusivity is decreased, which means that the athermal diffusivity due to electronic stopping persists at higher temperatures. In conclusion, the decrease in Xe diffusion means that more Xe stays in the matrix, decreasing the Xe release, and lowering its influence of fission gas release on the thermomechanical properties of UO 2 .

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Unusual aliovalent doping effects on oxygen non-stoichiometry in medium-entropy compositionally complex perovskite oxides

Aliovalent doping of perovskite oxides can tune the oxygen vacancy formation energy. This work discovers normal vs. abnormal aliovalent doping effects on redox behaviors in medium-entropy compositionally complex perovskite oxides (CCPOs) (La 1–x Sr x )(Mn 1/3 Fe 1/3 Ti 1/3 )O 3–δ (LS_MFT) vs. (La 1–x Sr x )(Mn 1/3 Fe 1/3 Cr 1/3 )O 3–δ (LS_MFC). In the LS_MFC series, the oxygen non-stoichiometry range Δδ (= δ red – δ ox ) linearly depends on the Sr molar ratio x, while the LS_MFT series shows a V-shape dependence of Δδ on x. This unusual observation is investigated and explained based on the analysis of the energy loss near edge structure (ELNES) in STEM electron energy loss spectroscopy, along with density functional theory (DFT) calculations. In LS_MFC, Cr–L 2,3 , Mn–L 2,3 and Fe–L 2,3 peaks have a similar linear shift to higher energy with increasing x, which indicates higher oxidation states of Cr, Mn, and Fe with lower oxygen vacancy formation energies. In LS_MFT, the V-shape of the Δδ vs. x curve is caused by the stable Ti 4+ state and a V-shape Mn/Fe valency dependence on x. Furthermore, this study suggests the possible existence of different (including unexpected) coupled aliovalent doping effects in CCPOs with multiple B-site redox active elements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Can a deep-learning model make fast predictions of vacancy formation in diverse materials?

The presence of point defects, such as vacancies, plays an important role in materials design. Here, we explore the extrapolative power of a graph neural network (GNN) to predict vacancy formation energies. We show that a model trained only on perfect materials can also be used to predict vacancy formation energies (E vac ) of defect structures without the need for additional training data. Such GNN-based predictions are considerably faster than density functional theory (DFT) calculations and show potential as a quick pre-screening tool for defect systems. To test this strategy, we developed a DFT dataset of 530 E vac consisting of 3D elemental solids, alloys, oxides, semiconductors, and 2D monolayer materials. We analyzed and discussed the applicability of such direct and fast predictions. We applied the model to predict 192 494 E vac for 55 723 materials in the JARVIS-DFT database. Our work demonstrates how a GNN-model performs on unseen data.

2D materials↗

Adsorbate-induced adatom formation on Au-Cu bimetallic alloys and its possible consequences for CO 2 electroreduction

The adsorbate-induced formation of sub-nanometer clusters on transition-metal single crystals observed in previous high-pressure microscopic studies hinted at the in-situ formation of unique active sites even on large nanoparticle catalysts. We propose that the adatom formation energy can be used as an energetic descriptor for the initial step toward the adsorbate-induced metal-cluster formation process. This descriptor can be efficiently computed using density functional theory (DFT) calculations and applied for screening and identification of metal catalysts where this phenomenon may play an important role in generating active sites in-situ. As a proof of concept, here, we construct an adatom formation energy database for three Au x Cu y alloys (x:y = 3:1, 1:1, or 1:3) and eighteen adsorbates (H, C, N, O, F, S, Cl, Br, I, CH x , NH x (x = 1 – 3), CO, NO, and OH) commonly involved in catalytic reactions. The energetics of adatom formation were examined in all cases where the (111) terrace, (211) step-edge, and (874) kink were the sources of the adatom. We demonstrate that the presence of an adsorbate could alter not only the energetics for adatom formation but also the elemental nature of the preferred adatom being formed. Using our database, we identified promising systems which favor adsorbate-induced adatom formation under near-ambient conditions. Specifically, CO-induced adatom formation on all three Au-Cu alloy surfaces could occur under CO 2 electroreduction (CO 2 RR) conditions. This phenomenon offers a qualitative explanation for the experimentally observed CO 2 RR activity on Au-Cu alloy catalysts. As a result, our methodology offers an easily expandable and efficient approach for large-scale catalyst screening with regards to adatom/cluster formation under reaction conditions and provides insight into the possible nature of active sites on alloy catalysts from a novel perspective.

Active site↗

Excitonic Condensate in Flat Valence and Conduction Bands of Opposite Chirality

Excitonic Bose-Einstein condensation (EBEC) has drawn increasing attention recently with the emergence of 2D materials. A general criterion for EBEC, as expected in an excitonic insulator (EI) state, is to have negative exciton formation energies in a semiconductor. Here, using exact diagonalization of a multiexciton Hamiltonian modeled in a diatomic kagome lattice, we demonstrate that the negative exciton formation energies are only a prerequisite but insufficient condition for realizing an EI. By a comparative study between the cases of both conduction and valence flat bands (FBs) versus that of a parabolic conduction band, we further show that the presence and increased FB contribution to exciton formation provide an attractive avenue to stabilize the excitonic condensate, as confirmed by calculations and analyses of multiexciton energies, wave functions, and reduced density matrices. In conclusion, our results warrant a similar many-exciton analysis for other known and/or new candidates of EIs and demonstrate the FBs of opposite parity as a unique platform for studying exciton physics, paving the way to material realization of spinor BEC and spin superfluidity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A quantitative phase-field model for gas bubble evolution in UO2

Due to the large formation energy of vacancies and noble gas atoms in the form of interstitials or substitutional atoms in nuclear fuel (UO2), the thermodynamic equilibrium concentrations of these species are very low in the nuclear fuel matrix even at very high temperatures, which imposes difficulties upon the quantitative study of bubble evolution via the phase-field method. In this study, a quantitative phase-field model is proposed to deal with this problem. The system’s free energy density is derived according to the principles of thermodynamics, with consideration of the elastic effect and with the use of real material parameters from experiments. The model is useful for the study of the kinetics of gas bubble growth with very dilute concentrations of vacancy and gas atoms in the matrix. This model is applied to study single bubble growth and multiple bubble growth under various concentrations of vacancy and gas atoms and at various temperatures. The elastic effect and the effects of the generation rate of vacancies and gas atoms on bubble growth are analyzed.

low defect concentration, nuclear fuel, gas bubble↗

Effect of impurities on hydrogen defect stability and migration barrier in yttrium dihydride crystal

The impurity or alloying atoms in YH 2 can alter the local electronic structure and so the hydrogen defect stability, as well as the H migration barrier energy. Thus, DFT calculations were employed to determine the effect of foreign elements from alkali and alkaline earth metals to transition metals and one critical impurity element, O, on H vacancy stability and retention characteristics in YH2. Results revealed that alloying elements act as hydrogen vacancy sinks by reducing the vacancy formation energy at neighboring sites. The implantation of non-magnetic foreign elements (s1, s2, and d10 valence electrons) in hydrogen energy landscape was calculated to be minor; while the hydrogen vacancy formation energy was reduced from 1.37 eV to 1.00 eV, the migration energy barrier of hydrogen was increased from 0.87 eV to 1.15 eV for non-magnetic foreign elements. The migration energy barrier monotonically decreased with increasing d-shell occupancy, reaching as low as 0.4 eV for Cr, Mo(d4), and Fe (d4). Alloying with late transition metals (d8 and d9) moderately impacted the hydrogen vacancy formation. Finally, it was found to be O addition into the YH 2- lattice did not alter the energy landscape of hydrogen vacancies. Since alloyed YH 2 has not been studied extensively, this study provides an atomistic understanding how alloying elements and impurities trap vacancies and affects hydrogen mobility YH 2 . Meanwhile, the main findings of this study may serve as guidelines for introducing alloying elements in ZrH 2 as well.

08 HYDROGEN↗

Elucidating the Impact of Cis – Trans Organic Structure Directing Agent Isomer Ratios on the Aluminum Distribution Within SSZ-39

Despite their widespread use, the mechanisms governing the synthesis of zeolite catalysts are still poorly understood. A notable example of this problem is the uncertainty surrounding the influence of synthesis conditions on the placement of Al atoms in the zeolite framework which determines the active sites available for catalytic species. In this work, the role of the cis to trans isomer ratio of the OSDA N,N-dimethyl-3-5-dimethylpiperidinium on the energetics of 26 distinct Al pair distributions in SSZ-39 is examined both in the presence and absence of Na using density functional theory calculations. The initial orientation of the OSDA was found to have a significant impact on the final energies present, necessitating the screening of a large number of initial orientations with force field calculations and single point DFT calculations. Ground state energies were found to vary significantly with the ratio of cis to trans OSDAs with a Boltzmann distribution revealing the most likely Al pair distributions shift from sharing the same 8 membered rings to sharing the same double six membered rings to having no shared subunits as one increases the amount of cis OSDA present within the framework. The presence of Na was found to favor Al pair distributions where both Als occupied the same 6-membered ring. When an implicit solvent model was used to evaluate ground state energies the ideal Na sites shifted from 6-membered rings to empty SSZ-39 cages while OSDA positions and orientations remained largely the same. To provide insight on how kinetic factors may influence Al distributions, formation energies we calculated for connected double six membered rings. Further, these formation energies revealed a preference for Al pairs to occupy the same 4-membered ring which indicates kinetic and thermodynamic control may lead to different Al distributions in SSZ-39.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Predicting Oxygen Off-Stoichiometry and Hydrogen Incorporation in Complex Perovskite Oxides

Chemically and structurally complex solid compounds, including those with significant off-stoichiometry, are rapidly extending new material functionality across a variety of applications. Accelerated development of these compounds requires accurate predictions of material defect properties including effective defect formation energies and equilibrium defect concentrations. Traditional first-principles approaches typically examine dilute defect concentrations and relatively ordered atomic structures to identify the lowest energy defect sites. These approaches are rarely suitable for describing the disorder present in these systems and its influence on defect formation, which can lead to unphysically large predictions for defect concentrations. Here, we demonstrate a new method to accurately predict the temperature and pressure dependence of oxygen vacancy concentrations and proton interstitial concentrations in complex oxides. This method extends standard dilute defect calculations to incorporate atomic and magnetic disorder, employs the ensemble descriptions of defect sites resulting in improved predictions of defect formation energies, and accounts for effects beyond the dilute defect limit. To demonstrate our method, we show that the predicted defect concentrations in perovskites used as ceramic fuel cell cathodes, including Ba0.5Sr0.5Fe0.8Zn0.2O3-d, Ba0.5Sr0.5Co0.8Fe0.2O3-d, and BaCo1-x-y-zFexZryYzO3-d, are in good agreement with experimental values, thereby opening the door for predictive design of complex oxides by these applications.

DFT↗