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
DOE OSTI2020

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)↗
DOE OSTI2020

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↗
DOE OSTI2025

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↗
DOE OSTI2022

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↗
DOE OSTI2023

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↗
DOE OSTI2022

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↗
DOE OSTI2023

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↗
DOE OSTI2024

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↗
DOE OSTI2023

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↗
DOE OSTI2020

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↗
DOE OSTI2024

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↗
DOE OSTI2024

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↗
DOE OSTI2022

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↗
DOE OSTI2020

The role of water vapor during the synthesis of hydrogen doped In 2 O 3

Hydrogen-doped In 2 O 3 synthesized using magnetron sputtering or atomic layer deposition in the presence of water vapor has high transparency in the near infra-red region and mobility values that are almost three to four times of those in commercially used In 2 O 3 :Sn at a much lower carrier concentration. However, simple questions like what happens when water molecules enter In 2 O 3 or why the dominant charge carrier changes from an oxygen vacancy in In 2 O 3 :Sn to H + in In 2 O 3 are not clear. Using hybrid functional based density functional theory calculations, in this study we show that water molecules spontaneously split into H + and OH – . The H + interstitial bonds with an oxygen in In 2 O 3 , and the OH – interstitial occupies a vacant lattice site or is trapped by an oxygen vacancy. Thus, in agreement with experimental findings, our results suggest that the concentration of oxygen vacancies decreases as the partial pressure of water increases. The charge transition level, ϵ ( – / + ) , of a hydrogen interstitial and the Fermi level at which the formation energies of OH – and H + are equal lie above the conduction band minimum, suggesting that the H + interstitial has a lower formation energy. In addition, water molecules do not exhibit any propensity to form clusters, meaning that H + ions can diffuse and remain spatially separated from each other.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Moving closer to experimental level materials property prediction using AI

Abstract While experiments and DFT-computations have been the primary means for understanding the chemical and physical properties of crystalline materials, experiments are expensive and DFT-computations are time-consuming and have significant discrepancies against experiments. Currently, predictive modeling based on DFT-computations have provided a rapid screening method for materials candidates for further DFT-computations and experiments; however, such models inherit the large discrepancies from the DFT-based training data. Here, we demonstrate how AI can be leveraged together with DFT to compute materials properties more accurately than DFT itself by focusing on the critical materials science task of predicting “formation energy of a material given its structure and composition”. On an experimental hold-out test set containing 137 entries, AI can predict formation energy from materials structure and composition with a mean absolute error (MAE) of 0.064 eV/atom; comparing this against DFT-computations, we find that AI can significantly outperform DFT computations for the same task (discrepancies of $$>0.076$$ > 0.076 eV/atom) for the first time.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Ligand removal energetics control CO 2 electroreduction selectivity on atomically precise, ligated alloy nanoclusters

Atomically precise, thiolate-protected gold nanoclusters (TPNCs) exhibit remarkable catalytic performance for the electrochemical reduction of carbon dioxide (CO 2 R) to CO. The origin of their high CO 2 R activity and selectivity has been attributed to partial ligand removal from the thiolate-covered surfaces of TPNCs to expose catalytically active sulfur atoms. Recently, heterometal doped (alloy) TPNCs have been shown to exhibit enhanced CO 2 R activity and selectivity compared to their monometallic counterparts. However, systematic studies on the effect of doping (metal type and location on TPNC) on active site exposure and CO 2 R activity are missing in literature. Herein, we apply Density Functional Theory calculations to investigate the effect of heterometal (Pt, Pd, Hg and Cd) doping of Au 25 (SR) 18 TPNC on the active site exposure and CO 2 R activity and selectivity. In this study, we reveal that doping significantly modifies relevant TPNC electronic properties, such as electron affinity, while also altering partial ligand removal and carboxyl (*COOH) intermediate formation energies. Furthermore, we demonstrate that changing the dopant (e.g. Hg) position can change the selectivity of the TPNC towards CO (g) or H 2(g) formation, highlighting the importance of dopant locations in TPNC-based CO 2 R. Most notably, we report a universal (i.e. capturing different dopant types and positions) linear trend between the ligand removal energy and i) the *COOH formation energy, as well as, ii) the hydrogen (*H) formation energy on the different alloy TPNCs. Thus, utilizing the ligand removal energy as a descriptor for CO 2 RR activity and selectivity, our work opens new avenues for accelerated computational screening of different alloy TPNCs for electrocatalytic CO 2 R applications.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

ORNL_AISD_NiNb

This dataset describes the nickel-niobium solid solution binary alloy, where the two constituent elements nickel (Ni) and niobium (Nb) are randomly placed on an underlying crystal lattice. This dataset for nickel-niobium (Ni-Nb) alloys available includes the formation energy and bulk modulus for each crystal structure. Each atomic sample has a disordered phase which is obtained starting from an initial regular crystal structure of type body-centered cubic (BCC), face-centered cubic (FCC), or hexagonal compact packed (HCP). The geometry optimization ensures that all the alloy samples reached the equilibrium with negative formation energy. We perform geometry optimizations using the LAMMPS simulation package [1], a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. We utilized the embedded atom model (EAM) potential for Ni and Nb developed in a previous study [2]. The potential could describe behaviors of the liquid and solid phases of Ni-Nb alloy. The structural factors and angular distributions of three atoms are well-matched with X-ray and ab initio-based molecular dynamics data. We prepared the three different crystals with different initial lattice parameters (3.52 Ã… for FCC, 3.32 Ã… for BCC, and 3.5 Ã… for HCP). We performed energy minimization in two steps. Firstly, we minimized the structures with an isotropic unit cell to minimize the side effects from our arbitrary lattice parameters for all other compositions. Then, we applied geometry optimization with a triclinic (non-orthogonal) unit cell to fully minimize the stress components to calculate the elastic constants. In this procedure, we chose 10,000 as the maximum number of allowable steps aimed at obtaining fully relaxed atomic geometries. The dataset consists of three sets of crystal structures. The first set contains 46,086 irregular crystal structures, each of them with 54 atoms, obtained through optimization starting from a regular BCC crystal structure. The second set contains 24,543 irregular crystal structures, each of them with 32 atoms, obtained through optimization starting from a regular FCC crystal structure. The third set contains 39,303 irregular crystal structures, each of them with 48 atoms, obtained through optimization starting from a regular HCP crystal structure. The atomic configurations within each set span the possible compositional range. The three sets have been unified in a global dataset, which is extremely heterogeneous in terms of crystal structures, lattice volumes, and atomic configurations. Organization of files inside the dataset: the dataset contains three subdirectories called • BCC_opt • FCC_opt • HCP_opt based on the type of initial regular structure used to start the geometry optimization. Inside each of these folders, every atomic structure is identified by a string “A_B_Câ€, where A denotes the number of Nb in the system, B denotes index of structure with a given Nb number, and C denotes the total number of structures generated with a given Nb number. For each optimized crystal structure identified by the unique string of characters “A_B_Câ€, three files are provided: • A_B_C_opt.xyz: The optimized geometries in xyz format • A_B_C_opt.cfg: The optimized geometries in cfg format. It includes cell information and atomic energy, and forces calculated from LAMMPS. • A_B_C.elastic: Raw data of 21 elastic constants from LAMMPS output. • A_B_C.bulk: Calculated upper and lower bounds of bulk modulus and averaged one based on Voigt-Reuss-Hill approach from *.elastic. References: [1] A. P. Thompson, H. M. Aktulga, R. Berger, D. S. Bolintineanu, W. M. Brown, P. S. Crozier, P. J. in 't Veld, A. Kohlmeyer, S. G. Moore, T. D. Nguyen, R. Shan, M. J. Stevens, J. Tranchida, C. Trott, and S. J. Plimpton. LAMMPS - a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. Comp. Phys. Comm., 271:108171, 2022. [2] Y Zhang, R Ashcraft, MI Mendelev, CZ Wang, and KF Kelton. Experimental and molecular dynamics simulation study of structure of liquid and amorphous ni62nb38 alloy. The Journal of chemical physics, 145(20):204505, 2016.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Understanding the interface interaction between U3Si2 fuel and SiC cladding

Abstract Triuranium disilicide (U 3 Si 2 ) fuel with silicon carbide (SiC) composite cladding is being considered as an advanced concept/accident tolerant fuel for light water reactors thus, understanding their chemical compatibility under operational and accident conditions is paramount. Here we provide a comprehensive view of the interaction between U 3 Si 2 and SiC by utilizing density functional theory calculations supported by diffusion couple experiments. From the calculated reaction energies, we demonstrate that triuranium pentasilicide (U 3 Si 5 ), uranium carbide (UC), U 20 Si 16 C 3 , and uranium silicide (USi) phases can form at the interface. A detailed study of U 3 Si 2 and SiC defect formation energies of the equilibrated materials yielding the interfacial phases U 20 Si 16 C 3 , U 3 Si 5 and UC reveal a thermodynamic driving force for generating defects in both fuel and cladding. The absence of either the U 3 Si 2 or SiC phase, however, causes the defect formation energies in the other phase to be positive, removing the driving force for additional interfacial reactions. The diffusion couple experiments confirm the conclusion with demonstrated restricted formation of U 3 Si 5 , UC, and U 20 Si 16 C 3 /USi phases at the interface. The resulting lack of continuous interaction between the U 3 Si 2 and SiC, reflects the diminishing driving force for defect formation, demonstrating the substantial stability of this fuel-cladding system.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗