Search NASA⌕ Search

SEARCH · Search NASA

Results for “activation energies”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Explainable machine learning for hydrogen diffusion in metals and random binary alloys

Hydrogen diffusion in metals and alloys plays an important role in the discovery of new materials for fuel cell and energy storage technology. While analytic models use hand-selected features that have clear physical ties to hydrogen diffusion, they often lack accuracy when making quantitative predictions. Machine learning models are capable of making accurate predictions, but their inner workings are obscured, rendering it unclear which physical features are truly important. To develop interpretable machine learning models to predict the activation energies of hydrogen diffusion in metals and random binary alloys, we create a database for physical and chemical properties of the species and use it to fit six machine learning models. Our models achieve root-mean-squared errors between 98–119 meV on the testing data and accurately predict that elemental Ru has a large activation energy, while elemental Cr and Fe have small activation energies. By analyzing the feature importances of these fitted models, we identify relevant physical properties for predicting hydrogen diffusivity. While metrics for measuring the individual feature importances for machine learning models exist, correlations between the features lead to disagreement between models and limit the conclusions that can be drawn. Instead grouped feature importance, formed by combining the features via their correlations, agree across the six models and reveal that the two groups containing the packing factor and electronic specific heat are particularly significant for predicting hydrogen diffusion in metals and random binary alloys. In conclusion, this framework allows us to interpret machine learning models and enables rapid screening of new materials with the desired rates of hydrogen diffusion.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Mechanistic study of chemical looping reactions between solid carbon fuels and CuO

Copper (Cu) based chemical-looping combustion (CLC) is a promising process that utilizes solid carbon fuel such as coal and biomass. Understanding the reaction kinetics in this process can facilitate the industrial design of CLC units. Here, in this work, molecular dynamics (MD) simulations were performed to investigate the reaction kinetics of n-butane and two simplified solid carbon fuels (lignite and anthracite) with/without copper oxide (CuO) nanoparticle. In addition, experiments were conducted on the thermal characteristics, flammability, and flame speeds of CuO and solid carbon mixtures. For n-butane, oxidation on the CuO surface is the primary reaction since the activation energy of the surface reaction is much lower than that of oxygen (O 2 ) combustion (9.2 vs. 53.3 kcal/mol). On the other hand, for lignite, there is a smaller difference in the activation energy of O 2 combustion and surface reaction (23.04 vs. 6.34 kcal/mol). We hypothesize that some solid carbon fuels can have different reaction kinetics dependent upon temperature. This is proven by the case of a simplified anthracite coal: an Arrhenius plot shows that this solid carbon fuel has two different reaction kinetic regimes and the critical temperature for the change in kinetics is related to the oxygen uncoupling of the CuO nanoparticle. Like the modeling simulations, a change in activation energy is observed in the experimental results, where desorption of molecular oxygen from CuO becomes important.

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

An energetic link between order and strength in metals: A nanocrystalline strength limit in high-entropy alloys and intermetallic compounds

The metallurgy and materials communities have long understood and exploited fundamental links between chemical and structural ordering in metallic solids to tailor their mechanical properties. We extend these ideas to include prediction of the nanocrystalline strength limit in high-entropy alloys and intermetallic compounds, where a breakdown occurs in the classical Hall-Petch strengthening behavior. The highest reported strength achievable through alloying has rapidly climbed and given rise to new classifications of materials with extraordinary properties, with a notable case being nanocrystalline metals. High-entropy alloys (chemically disordered, concentrated solid solutions) and intermetallic compounds are two boundary cases of how tailored order can be used to manipulate mechanical behavior. Here, we show that the complex electronic-structure mechanisms governing the peak strength of alloys and pure metals can be reduced to a few physically meaningful parameters based on their atomic arrangements and used – with no fitting parameters – to predict the maximum strength of these materials. This includes a generalized energy-based accounting for the degree of structural and chemical ordering that allows for rapid and reasonably accurate prediction of peak strength (validated in the nanocrystalline limit) as a function of temperature. Predictions of maximum strength based on the activation energy (with all materials properties derived from DFT calculations or experiments) for a stress-driven transition to an amorphous state is shown to accurately describe the breakdown in Hall-Petch behavior at the smallest crystallite sizes for pure metals, intermetallic compounds, high-entropy alloys, and metallic glasses. Further, this activation energy is also shown to be directly proportional to interstitial electronic charge density, which is a good predictor of ductility, stiffness (moduli), and phase stability in high-entropy alloys and solid metals generally. The proposed framework suggests the possibility of coupling ordering and intrinsic strength to mechanisms like dislocation nucleation, hydrogen embrittlement, and transport properties, such as through correlations between the activation energies for amorphization with stacking-fault and grain boundary energies. It additionally opens the prospect for greatly accelerated structural materials design and development to address materials challenges limiting more sustainable and efficient use of energy.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Stability of Cu(In x Ga 1− x )Se 2 Solar Cells Utilizing RbF Postdeposition Treatment under a Sulfur Atmosphere

Alkali halide postdeposition treatments (PDTs) have become a key tool to maximize efficiency in Cu(In x Ga 1− x )Se 2 (CIGS) photovoltaics. RbF PDTs have emerged as an alternative to the more common Na‐ and K‐based techniques. This study utilizes temperature‐dependent current–voltage ( JVT ) measurements to study a unique RbF PDT performed in a S atmosphere. The samples are measured before and after 6 months in a desiccator to study device stability. Both samples contain Na and K which diffuse from the soda–lime glass substrate. A reference sample and a RbF + S PDT sample both show the development of a rear contact barrier after aging. The contact barrier is higher for the RbF + S PDT sample, leading to decreased current in forward bias. Series resistance is also higher in the RbF + S PDT device which leads to lower fill factor. However, after aging the reference sample has a larger decrease in open‐circuit voltage ( V OC ). Ideality factor measurements suggest Shockley–Read–Hall recombination dominates both samples. V OC versus temperature and a temperature‐dependent activation energy model are used to calculate diode activation energies for each sample condition. Both techniques produce similar values that indicate recombination primarily occurs within the bulk absorber.

14 SOLAR ENERGY↗
DOE OSTI2022

Phase transformation and growth mechanism of RF sputtered ferroelectric lead scandium tantalate (PbSc 0.5 Ta 0.5 O 3 ) films

Lead scandium tantalate (PbSc 0.5 Ta 0.5 O 3 , PST), an order/disorder ferroelectric, is a potential candidate for electrocaloric cooling and pyroelectric infrared (IR) detector. In this work, we report the phase transformation kinetics from two series of samples containing pure amorphous and mixture of amorphous and pyrochlore to desired perovskite phase using postdeposition rapid thermal processing (RTP) as well as growth mechanism of RF sputtered PST thin films using excess lead target on platinized silicon (Pt/Ti/SiO 2 /Si) substrates. We find that small changes in the temperature ramp have a large effect on the degree of perovskite conversion (ferroelectric phase), orientation (crystallographic texture), and long-range order parameter (< S 111 >). Through isothermal annealing, we obtained optimal perovskite phase at ≥700°C temperature. The phase transformation is characterized by spontaneous formation of center-type in-plane radial rosette-like structures revealed by scanning electron microscopy. The PST perovskite crystallites were found to coexist with pyrochlore in RTP annealed films. The volume fractions for perovskite and pyrochlore phase were obtained from the analysis of “rosettes” and respective X-ray diffraction intensities which helped to determine various parameters associated with phase kinetics (n, k, and activation energy, Ea) and accompanying growth. The effective activation energies of perovskite transition and growth were found to be 332 ± 11 kJ/mol (345 ± 11 kJ/mol) and 114 ± 10 kJ/mol (122 ± 10 kJ/mol), respectively, for pure amorphous only (and mixed amorphous and pyrochlore) phase following nucleation-growth controlled Avrami's equation. A linear growth rate (n~1) for the perovskite phase indicates predominant interface-controlled process and diffusion-limited phenomena thus inhibiting rosette size owing to reactant depletion and soft impingement at the grain boundary. However, the growth behavior is isotropic in two-dimension parallel to the plane of the substrates for both sample series. Furthermore, lead loss was severe for in-situ growth and RTP combined with conventional furnace annealing than those of RTP only films, which were closer to stoichiometric albeit with excess lead and marginal oxygen vacancies (V o ).

36 MATERIALS SCIENCE↗
DOE OSTI2024

Development of Leaching Processes for the Recovery of Rare Earth Elements from Acid Mine Drainage Precipitates

Acid mine drainage (AMD) has been a challenge for mine operators to address and has had significant environmental impacts when there is a failure to address it. Fortunately, there is a regulation in the US that requires the treatment of AMD before water can be discharged into the environment. AMD is generated when sulfide minerals are exposed to water and air from mining. The acidic water then leaches metals from the surrounding rock, creating the potential for environmental contamination of this acidic water containing dissolved metals. AMD can contain high concentrations of metals like iron, aluminum, and manganese and also have been shown to contain trace concentrations of critical minerals, including rare earth elements (REEs). This environmental waste stream is now being researched as a potential source for REEs. There is a patented process for processing and extracting REEs from AMD which produces rare earth oxide preconcentrate (REOP) from AMD treatment precipitates and a final stage mixed rare earth oxide (MREO) product. This body of research examines a selective leaching process for each of these products and determines the activation energy associated with the leaching processes. Acid leaching with HCl was examined to selectively leach REEs from REOP while contaminant metals remained in the solid residue. The maximum leaching recovery of the REEs from the REOP was approximately 90% and an activation energy of 6.4 kJ/mol at pH 3.0. Ammonium chloride leaching was examined to selectively remove contaminant metals from a MREO product. The ammonium chloride process successfully leached major contaminant metals in excess of 85% recovery and had a maximum activation energy of 40.0 kJ/mol.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2025

Low-energy pathways lead to self-healing defects in CsPbBr 3

Self-regulation of free charge carriers in perovskites via Schottky defect formation has been posited as the origin of the well-known defect-tolerance of metal halide perovskite materials. Understanding the mechanisms of self-regulation promises to lead to the fabrication of better performing solar cell materials with higher efficiencies. We investigate many such mechanisms here for CsPbBr 3 , a popular representative of a more commercially viable all-inorganic metal halide perovskite. We investigate different atomic-level mechanisms and pathways of the diffusion and recombination of neutral and charged interstitials and vacancies (Schottky pairs) in CsPbBr 3 . We use nudged elastic band calculations and ab initio-derived pseudopotentials within quantum ESPRESSO to determine energies of formation and migration and hence the activation energies for these defects. While halide vacancies are known to exhibit low formation energies, the migration of interstitials is less studied. Our calculations uncover interstitial defect pathways capable of producing an activation energy at, or below, the single experimental value of 0.53 eV observed for the slow, temperature-dependent recovery of light-induced conductivity in bulk CsPbBr 3 . Our work reveals the existence of a low-energy diffusion pathway involving a concerted “domino effect” of interstitials, with the net result that interstitials can diffuse more readily over longer distances than expected. This observation suggests that defect self-healing can be promoted if the “domino effect” strategy can be engaged.

14 SOLAR ENERGY↗
DOE OSTI2020

Kinetic Parameters and Diffusivity of Uranium in FLiNaK and ClLiK

Anodic stripping voltammetry (ASV) and cyclic voltammetry (CV) measurements at 773, 823, and 873 K were made of uranium trifluoride (UF 3 ) in lithium fluoride-sodium fluoride-potassium fluoride eutectic (FLiNaK) and uranium trichloride (UCl 3 ) in lithium chloride-potassium chloride eutectic (ClLiK). ASV data were used to estimate the charge transfer coefficients, exchange current densities, and activation energies of the uranium reactions. Charge transfer coefficients of both salt systems were within the range of 0.17 to 0.38. Exchange current densities in the fluoride and chloride salts were estimated within the range of 0.060 to 0.12 A cm -2 . Activation energy of uranium exchange current was 16.0 kJ mol -1 in the fluoride salt and 28.9 kJ mol -1 in the chloride salt. Kinetics of charge transfer were found to be faster in FLiNaK. Analyses of the CV data suggest the electrochemical system was diffusion controlled and irreversible in the chloride and fluoride salt mixtures. Diffusion coefficients of uranium in the range of temperatures were on the order of 10-5 cm 2 s -1 in both systems. Greater values of diffusivity in ClLiK are attributed to its lower density compared with FLiNaK. Activation energy of uranium diffusion in the fluoride and chloride salt mixtures were 53.4 and 89.4 kJ mol -1 , respectively.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Particle Markov Chain Monte Carlo Approach to Inference in Transient Surface Kinetics

Here, in this work, we develop a novel Bayesian approach to study the adsorption and desorption of CO onto a Pd(111) surface, a process of great importance in natural sciences. The motivation for this work comes from the recent availability of time-resolved infrared spectroscopy data and the need for model interpretability and uncertainty quantification in chemical processes. The objective is to learn the relevant parameters that characterize the process: coverage with time, rate constants, activation energies, and pre-exponential factors. Our approach consists of three main schemes: (i) a problem design and probabilistic model for the whole system, (ii) a particle Markov chain Monte Carlo sampler to learn the hidden coverages and rate constant parameters, and (iii) two Bayesian formulations to infer the activation energies and pre-exponential factors. The flexibility of the Bayesian framework allows for uncertainty quantification where possible and integration of mathematical constraints in the model to reflect the system physically. We found that our results for the activation energies and pre-exponential factor are in agreement with those reported in the experimental literature, independently, and we provide discussions on the advantages and disadvantages as well as applicability to other systems.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Temperature Effects on Interdiffusion of Al and U-Mo under Irradiation

A high-energy Xe ion irradiation experiment was conducted to investigate the temperature dependence of interdiffusion in bilayer Al-UMo samples under irradiation. The amount of interdiffusion achieved at a fixed dose with the increase of temperature showed a clear transition at 175 degrees C (with an estimated error in the range of ± 10 degrees C) from temperature-independent to temperature-dependent behavior. The activation energy derived from the curve of interdiffusion quantity vs. irradiation temperature is 0.77 ± 0.16 eV. This information has been utilized to understand the temperature effect on the interdiffusion process that occurred at the interfaces of U-Mo particles and the Al matrix in U-Mo/Al dispersion fuels, whose magnitude significantly impacts the fuel's performance. Although this temperature effect was deemed important, it cannot be examined directly using in-pile irradiation data, as fuel temperatures cannot be measured in reactor irradiation and are highly correlated with fission rate and thermal conductivity evolution. To connect the knowledge accumulated from ion irradiation with in-pile irradiation data, simulation of a full-sized U-Mo/Al dispersion fuel plate irradiated in the FUTURE test in the BR2 reactor was performed with the Dispersion Analysis Research Tool (DART), a dispersion fuel performance code. DART is equipped with an interaction or interdiffusion layer (IL) growth correlation formulated to describe the temperature dependence of ion mixing results. The agreement between calculated and measured fuel meat constituent volume fractions and swelling data demonstrated that the temperature effect on in-pile Al-UMo interdiffusion is well captured with the correlation. In this case, the fitted activation energy is 0.70 eV. Considering the uncertainties associated with the ion irradiation data, the activation energy obtained from in-pile data fitting is in accord with that from ion irradiation results.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Reduction of Iron (III) Oxide in Microwaves Toward Gasification Studies

The purpose of this study is to reduce hematite (Fe2O3) using hydrogen (H2), carbon monoxide (CO), and methane (CH4) gases using microwave irradiation. Conventional temperature programmed reduction of iron phase transformation under H2, CO, and CH4 agree with the literature. Furthermore, the activation energy for the iron phase transformation in a conventional reactor was in the order of H2>CO>CH4, suggesting that CH4 was the best reductant to initiate the reduction at a much lower temperature than CO, and CO was significantly better than H2 atmosphere. Due to the uneven temperature distribution, the apparent activation energy of microwave reduction under H2 was approximately 1/4th when compared to the conventional reduction in H2. X-ray diffraction (XRD) analysis showed a mixed oxide phase in microwave reduction as opposed to the clear phase transformation in conventional studies. To combat the temperature measurement issues in microwave studies, fiber optic sensors and Forward Looking InfraRed (FLIR) sensors were used to acquire the actual temperature distribution in the catalyst bed. Therefore, a reasonably accurate activation energy estimation was achieved in the catalyst bed. Combining these with the characterization and quantification of the different iron phases in microwave reduction will provide a blueprint for future microwave-assisted gasification studies.

Aireddy, Divakar Reddy↗
DOE OSTI2024

In-situ hydrogen microstructural characterization of Si heterojunction passivation: Addressing V OC degradation and mitigation pathways

Si heterojunction (SHJ) solar cells have demonstrated record efficiency >27%, approaching the theoretical limit of ≈ 29%, primarily due to best surface/interface defect passivation provided by deposited thin layers of hydrogenated amorphous silicon (a-Si:H). Such excellent surface/interface passivation reduces recombination loss and result in >100 mV improvement of cell open circuit voltage (V OC ) to ≈ 750 mV, thus the cell efficiency. However, fielded SHJ modules exhibit loss of V OC and hence efficiency over time in years, presumably due to degradation related to a-Si:H layers. This adversely affects the technology’s market acceptance, and levelized cost of energy (LCOE). It is hypothesized that the origin of a-Si:H degradation is somehow related to the presence of weak Si–Si bonds and hydrogen in a-Si:H films. The objective of this project is to test this hypothesis by directly measuring chemical and structural changes occurring within SHJ component layers and solar cells. This is achieved by developing an innovative in-situ Fourier transform infrared (FTIR) spectrometry apparatus to monitor hydrogen microstructural changes occurring within amorphous silicon and decipher hydrogen evolution kinetics over time when samples are exposed to heat and/or light stress. These in-situ measured hydrogen microstructural changes are correlated to the changes in effective minority carrier lifetime (τ eff ), implied V OC (iV OC ), surface recombination velocity (S), and cell V OC . These mechanistic understandings will provide critical guidance to mitigate the V OC -driven degradation of SHJ solar cell performance. Passivation optimization and degradation analysis of individual SHJ component structures were achieved through systematic deposition of three symmetric structures and the completed SHJ solar cell structure. The three symmetric structures used were intrinsic a-Si:H [(i)a-Si:H] layers in a bilayer structure, intrinsic and p-type doped stacked layers [(i-p)a-Si:H] representing the front heterojunction in the SHJ cell, and intrinsic and n-typed doped stacked layers [(i-n)a-Si:H] representing the back-side back surface field (BSF) in the SHJ cell. State-of-the-art passivation qualities are demonstrated by a champion iV OC of 740 mV for the (i)a-Si:H layers, and the (i-n)a-Si:H symmetric structure. A 725 mV iV OC is observed for the (i-p)a-Si:H symmetric structure. These symmetric passivated SHJ component structures were subsequently subjected to different accelerated lifetime (ALT) stressors to identify which conditions contribute the most to iV OC degradation. Degradation of the thin (10 nm) (i)a-Si:H passivation layers without any additional overlying layers is minimal; complexity of this study arises due to unavoidable surface oxidation of (i)a-Si:H layer during most of the stress application, which is likely irrelevant for a full SHJ cell configuration with overlying protective layers. The iV OC degradation of symmetric structures is found to occur primarily at the (i-p)a-Si:H passivation stack under dark heat stress with associated hydrogen loss from the (p)a-Si:H layer. An activation energy for increase in S (defect creation) of 0.65 eV can be correlated to the activation energy of ≈ 0.4 eV for hydrogen loss from the (i-p)a-Si:H stack. This also suggests the presence of weakly bonded hydrogen in the (p)a-Si:H films, which effuses out of the film stack at such low activation energy. When light and heat stress are applied together, similar hydrogen loss from (i-p)a-Si:H stack is observed, however, does not appreciably degrade iV OC or increase S. This is an important result and departure from direct correlation between hydrogen loss and defect creation. This perhaps indicates additional defect chemistries or annealing that might be occurring in the presence of light requiring further detailed defect measurements. The full SHJ cell structure used for this project is depicted in Fig.1(d). SHJ cells with an initial V OC ≈ 700 mV were fabricated and subjected to similar ALT stress conditions. Cell V OC is found to degrade the most under dark heat stress and is confirmed by observed hydrogen migration out of the (i-p)a-Si:H stack. However, hydrogen cannot escape from the cell stack, it accumulates near the (p)a-Si:H/ITO contact interface, where ITO acts as a barrier preventing hydrogen loss. Furthermore, light-heat combined stress does not degrade V OC appreciably, confirming the occurrence of a defect annealing process.

14 SOLAR ENERGY↗
DOE OSTI2025

Investigations on the Thermal Stability and Kinetics of Biolubricants Synthesized from Different Types of Vegetable Oils

Petroleum-based lubricants raise environmental concerns due to their non-biodegradability and toxicity, whereas biobased lubricants underperform owing to low thermal stability. This study examined and compared three vegetable oils, along with their chemically modified versions, to better understand their suitability as biolubricants. High oleic soybean oil (HOSOY), regular soybean oil (RSOY), and waste cooking oil (WCO) were subjected to chemical modification, where isopropyl groups were attached to the fatty acid chains of the oils to produce branched oils, i.e., b-HOSOY, b-RSOY, and b-WCO. The detailed kinetic study of each regular and modified sample was investigated using thermogravimetric analysis. The kinetic parameters, such as the activation energies, reaction rate, and pre-exponential factor, were generated via Friedman methods. The differential thermal gravimetric (DTG) analysis showed low volatilization at the onset temperature in each modified oil as compared with the unmodified samples under an oxidative environment. Furthermore, the comparative kinetic studies demonstrated the enhanced thermoxidative stability of the modified products relative to their unaltered counterparts. Among the tested oils, the b-RSOY showed an average activation energy of 325 kJ/mol, followed by the b-WCO: 300 kJ/mol and the b-HOSOY: 251 kJ/mol, indicating the most stable modified product under an oxidative environment. For all the samples, the pre-exponential factors were in good agreement with the activation energies, which validates that finding the pre-exponential components is crucial to the kinetic analysis.

Sarker, Majher I. (ORCID:0000000299509274)↗
DOE OSTI2026

A Unified Workflow for Sensitivity-Based Kinetic Analysis in Microkinetic Models

Degrees of rate control (DRC), apparent activation energies, and apparent reaction orders are established local sensitivity diagnostics for interpreting microkinetic models, but applying them routinely to large mechanisms often requires substantial reaction-specific bookkeeping, perturbation design, and postprocessing. Here, in this study, we present a unified derivative-based workflow that evaluates these quantities from a single compiled reaction-network model and target-rate definition. For any user-provided microkinetic model, the workflow compiles the mechanism into stoichiometrically consistent mass-action rate equations, solves the surface dynamics, and uses automatic differentiation to compute sensitivities with respect to rate constants, temperature, and gas partial pressures. By combining their calculations in the same framework, the workflow clearly demonstrates the relationships between different DRCs and the apparent activation energy. Using existing examples of propylene partial oxidation and methane oxidation on Pd(100), we verify expected transient redistribution of rate control, distinguish net Campbell DRCs from one-sided directional sensitivities, and show how apparent activation energy can be reconstructed either from one-sided DRCs or from state-based DRCs while critical mechanistic insights are obtained consistently. In the methane oxidation case, a pathway-subset test further illustrates how a simplified mechanism preserves key kinetic signatures of a full model, showing the potential of our user-friendly tool for model construction beyond kinetic analysis.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Development and application of a uranium mononitride (UN) potential: Thermomechanical properties and Xe diffusion

We report that atomic-scale modeling of thermophysical and defect properties of uranium mononitride (UN) plays an important role in establishing a better understanding and improved models of UN fuel performance. Having an accurate interatomic potential is crucial for generating reliable data at finite temperatures using molecular dynamic simulations. We report a new interatomic potential for UN, based on a combination of many-body and pairwise interactions, a simple form that we later show could be easily adapted to include Xe-U and Xe-N interactions, i.e., generating a U-N-Xe interatomic potential. The potential was fitted to experimental thermal expansion and single crystal elastic constants, as well as Frenkel, Schottky, anti-Schottky, and antisite pair reaction energies from density functional theory (DFT) calculations. Using the potential, we successfully reproduced experimental lattice parameters, thermal expansion, single crystal elastic constants, and temperature dependent heat capacity. The potential also performs reasonably well in reproducing the energy of the aforementioned stoichiometric defect reactions and defect migration barriers calculated using DFT. However, the potential underestimates the energy difference between the tetrahedral and dumbbell uranium interstitials, and a more complex potential form might be needed to overcome this issue. The potential was also used to predict UN single crystal elastic constants and elastic properties at different temperatures, showing that UN becomes softer and more compressible with increasing temperature. We also compare our potential against literature data from previous empirical potentials, demonstrating similar or better behavior depending on the property of interest. To enable the simulation of Xe in UN a Buckingham potential has been fitted to DFT-derived Xe incorporation energies. The potential was then used to determine the activation energy for Xe diffusion due to various Xe-containing defects, with {XeU : VU} exhibiting the lowest activation energy.

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

Oxidation behaviors of matrix-grade graphite during water vapor ingress accidents for high temperature gas-cooled reactors

Water ingress into cores of high temperature gas-cooled reactors (HTGRs) can cause serious safety problems in graphitic components, especially for matrix graphite, which is more vulnerable than nuclear graphite. Here in this study, oxidation behaviors of a matrix-grade graphite, ARB-B1, were investigated under simulated water ingress accidental conditions. The oxidation tests were conducted in water vapor/helium mixed atmospheres at elevated temperatures up to 1200 °C. Oxidation rates and activation energy were evaluated by mass loss measurements. The activation energy was 121.5 kJ/mol in the considered temperature range, indicating only one oxidation regime-in-pore diffusion assisted by cracks. The cracks originated from the fabrication process, acting as penetration paths for the oxidant, played an important role in the oxidation behaviors. Ungraphitized binder oxidation occurred at all temperatures, while apparent oxidation of filler particles started from 1100 °C. Moreover, oxidation effects on compressive properties and hardness of the matrix graphite were discussed in terms of oxidized microstructures. This first detailed study on matrix graphite oxidation in water vapor provides important kinetics knowledge and microstructural and mechanical data for accidental scenarios of HTGRs as well as for development of new matrix graphite materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Biome-scale temperature sensitivity of ecosystem respiration revealed by atmospheric CO2 observations

Abstract The temperature sensitivity of ecosystem respiration regulates how the terrestrial carbon sink responds to a warming climate but has been difficult to constrain observationally beyond the plot scale. Here we use observations of atmospheric CO 2 concentrations from a network of towers together with carbon flux estimates from state-of-the-art terrestrial biosphere models to characterize the temperature sensitivity of ecosystem respiration, as represented by the Arrhenius activation energy, over various North American biomes. We infer activation energies of 0.43 eV for North America and 0.38 eV to 0.53 eV for major biomes therein, which are substantially below those reported for plot-scale studies (approximately 0.65 eV). This discrepancy suggests that sparse plot-scale observations do not capture the spatial-scale dependence and biome specificity of the temperature sensitivity. We further show that adjusting the apparent temperature sensitivity in model estimates markedly improves their ability to represent observed atmospheric CO 2 variability. This study provides observationally constrained estimates of the temperature sensitivity of ecosystem respiration directly at the biome scale and reveals that temperature sensitivities at this scale are lower than those based on earlier plot-scale studies. These findings call for additional work to assess the resilience of large-scale carbon sinks to warming.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Deep chalcogen donors and electron localization in aluminum nitride

Abstract We show with hybrid density functional theory calculations that chalcogen donors other than oxygen (i.e. S N , Se N , and Te N ) give rise to deep donor states in aluminum nitride. These donors trap a localized electron in their neutral charge state, leading to deep ( + / 0 ) donor levels that are 0.45 eV or more from the conduction-band edge. As such, this behavior is distinct from the DX behavior leads to deep (+/−) levels which affects other donors such as O N and Si Al . We highlight how these results hint at the formation of small electron polarons in AlN, which are found to be unstable in the bulk, but metastable when bound to donor dopants like Si Al and the chalocogens, with activation energies on the order of 0.2–0.3 eV. These results indicate that S, Se, and Te are not shallow donor dopants in aluminum nitride and identify origins of the experimentally observed ∼200–300 meV activation energies for dopant activation in donor-doped samples.

aluminum nitride↗