Search NASA⌕ Search

SEARCH · Search NASA

Results for “Fuel property prediction”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Properties, Theory, and Measurements for Understanding the Function of Heavy Elements (Project 20180474CR, Final Technical Report)

The FY18-20 Laboratory Directed Research and Development (LDRD) funded Properties, Theory, and Measurements for Understanding the Function of Heavy Elements Project aim has been to advance heavy element science in a comprehensive project connecting targeted research with Los Alamos mission imperatives. The mechanisms for project implementation were (1) support of a multi-year broad-based postdoctoral fellows project, (2) a summer student fellows research project and two seminar series, and (3) organization of structured workshops and integrated educational and career development opportunities. The project was successfully administered by the Glenn T. Seaborg Institute (GTSI) for Transactinium Science and has shown to be an exceptional investment in supporting the Los Alamos National Laboratory (LANL) mission areas requiring heavy element and actinide science by attracting and funding the future generation of scientists and engineers. This project exhibited its value across the laboratory in heavy element and actinide science areas include materials, material properties, signatures, modeling, predictions, fabrication, detection, disposal, global security implications, forensics, and the specialized science surrounding plutonium, uranium, and their surrogates as fuels for energy and nuclear weapons. In this report, the GTSI project staff documents the “high return on investment” implementation as measured against the successes of previous actinide research and science-based projects and as evidenced by the retention of 30% of G.T. Seaborg Postdoctoral Fellows as permanent employees. With additional funding through LDRD Reserve funding during this project period, the GTSI expanded its support for heavy element and actinide science research with a Visiting Research Seaborg Scholar and a Rapid Response Small Projects call. These innovations have served to advance GTSI pursuits as communicated in our FY21-23 LDRD Proposal Project#: 20210527CR Seaborg Institute: Center for Advancing Actinide Science and Technology at LANL

36 MATERIALS SCIENCE↗
DOE OSTI2021

Irradiation-enhanced diffusion and diffusion-limited creep in U 3 Si 2

U 3 Si 2 is an advanced fuel candidate due to its relatively high fissile density and attractive thermal properties. Compared to standard UO 2 fuel, there are significant data gaps for the thermophysical and thermomechanical properties of U 3 Si 2 . Point defect concentrations and mobilities under irradiation govern a number of important fuel performance properties, such as creep and fission gas release. In this work, we utilized density functional theory (DFT) data to inform a cluster dynamics framework to predict point defect concentrations in U 3 Si 2 under irradiation. Molecular dynamics (MD) simulations were used to examine the contribution of atomic mixing during ballistic cascades to diffusion, as well as the diffusivity of U and Si at grain boundaries. These atomic scale models for diffusivity were then used to inform a creep model based on bulk (Nabarro-Herring) and grain boundary (Coble) diffusional creep, and climb-limited dislocation creep. The model compares well against available experimental data and has been implemented in the BISON fuel performance code. Finally, a demonstration case using simple power profiles has been carried out, showing that negligible creep occurs due to the low temperatures experienced by U 3 Si 2 in-reactor, a consequence of its high thermal conductivity.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Are we modeling the properties of unconventional shales correctly?

Technological advancements have led to impressive growth in hydrocarbon (HC) production from unconventional shale reservoirs. However, major limiting factors in shale gas development are i) low recovery of gas in place (<20%), (ii) a rapid decline in well productivity, iii) release of contaminants, despite the use of advanced hydraulic fracturing fluids and multiple fracturing stages. This is mainly due to a lack of understanding of the nanoscale physicochemical properties of shales, especially that of kerogen, the macromolecule that is not only the source of the majority of the HC’s in shales but also holds most of the HC’s in adsorbed state. Over the last few years, a large number of studies have been published on simulating the physicochemical properties of shale utilizing molecular models of kerogen. However, the molecular models of kerogen input into these simulations are based on the kerogen “type” derived from a very limited number of shale samples. In this paper, we examine the variations that can exist in kerogen structure within a particular kerogen “type” across different shale basins and single shale basin at similar thermal maturity levels. We propose that using kerogen “type” based structural models for molecular simulations could lead to inaccurate estimation of HC reserves, HC recovery, HC production, frackability, and quantity and quality of produced water. Our study highlights the need for developing a better classification of kerogen based on its “molecular structure” instead of “type” for a more accurate prediction of physicochemical properties of shales.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2024

BISON Simulated and Experimental Fission Product Release Comparisons from Reradiated AGR-3/4 Compacts During High Temperature Heating Tests

The fuel performance modeling code BISON was used to predict the release of fission products iodine-131 (131I), xenon-133 (133Xe), and krypton-85 (85Kr) from four re-irradiated AGR-3/4 fuel compacts containing tristructural isotropic (TRISO) coated particles during high-temperature isothermal heating tests. The AGR-3/4 fuel compacts were irradiated in the Advanced Test Reactor (ATR) as part of the third and fourth series of planned experiments to support the Advanced Gas Reactor (AGR) Program. They were subsequently stored and re-irradiated in the Neutron Radiography (NRAD) reactor for approximately five days and then stored for another five to eight days before being subjected to isothermal heating tests in the Fuel Accident Condition Simulation (FACS furnace) for 200 to 300 hours at temperatures between 1000°C and 1600°C to evaluate fission product release at elevated temperatures. New nuclide-specific fission product source term models for the three nuclides of interest were developed using the reactor multiphysics code Griffin and implemented into BISON to support this work. The new source term models were incorporated into coupled compact- and particle-scale BISON simulations, which predict spatially- and temporally-resolved radionuclide generation, radioactive decay, transport, and release throughout the entire irradiation history, including the initial ATR irradiation, NRAD re-irradiations, FACS heating tests, and intermediate periods spent in storage. The experimentally measured fission product release from the heating tests were compared to modeling release predictions calculated by BISON to evaluate how well the code compares to experimental results. Overall, the experimental measured and BISON predicted comparative release results varied but generally agreed to within 5 particle equivalents. Comparative release results identified general observations to take into consideration to help refine future models and reduce uncertainties associated with both the measurement results and predictive results. This includes developing new uranium oxycarbide (UCO) specific kernel diffusivities for the three isotopes examined to more accurately reflect the material properties of the fuel form. Deriving new diffusivities will aid in producing a more informed BISON model

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

3D TRISO particle-explicit compact meshing

The TRI-structural ISOtropic (TRISO) layered fuel particle is a robust nuclear fuel form offering enhanced safety and performance for advanced reactor concepts, including high-temperature gas-cooled reactors and other Generation IV designs. These poppy-seed-sized particles are embedded in a graphite matrix to form fuel elements that must withstand elevated temperatures and high burn-up levels. The heterogeneous nature of these fuel elements — comprising thousands of randomly distributed TRISO particles — produces complex stress fields and thermal gradients that one- and two-dimensional models cannot accurately capture. While three-dimensional modeling has improved predictions of dimensional changes, internal pressure buildup, and fission product transport under irradiation, current approaches rely on homogenized material properties that are known to have considerable divergence from experimental observations. This work presents a methodology for optimized random packing of TRISO fuel compacts and full three-dimensional mesh generation within the BISON fuel performance code, with each particle coating layer individually discretized. The resulting mesh was demonstrated through heat conduction simulations under representative in-reactor operating conditions, showing strong agreement with expected behavior. This capability enables detailed analysis of particle-to-particle interactions, matrix cracking mechanisms, and the statistical distribution of coating layer failures — all of which directly govern fuel performance and safety margins.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Actinide Radiation Chemistry and Used Nuclear Fuel Reprocessing

Although the actinides boast many unique physical and chemical properties, their inherent susceptibility to radioactive decay—and subsequent consequences of radiation-induced chemistry—are what make them truly interesting and challenging elements to understand. From a closed nuclear fuel cycle perspective, the ability to predict and control the effects of actinide-driven radiolysis is critical for the design, development, and deployment of advanced used nuclear fuel reprocessing strategies and technologies. The absorption of ionizing radiation from actinide decay leads to the formation of a variety of transient and steady-state radicals, ions, and molecular radiolysis products that can lead to significant changes in a reprocessing solvent system’s physical and chemical properties, which ultimately limits that process’ efficiency and longevity. Presented here is an overview of recent advances in actinide radiation chemistry as it applied to used nuclear fuel reprocessing.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2021

Model predictive combustion control of a Gasoline Compression Ignition engine

Gasoline Compression Ignition is a novel combustion concept that derives its superiority from the high compression ratio of a compression ignition engine as well as the properties of gasoline fuel, such as longer ignition delay and higher volatility compared to diesel. Here, this combustion concept was experimentally tested on a 12.4L Class 8 truck engine. Based on these experimental data, prior efforts by the authors focused on the development of an engine model for a heavy-duty engine operating on a low-reactivity fuel. This engine model was leveraged within this study to investigate a combustion control strategy at different engine conditions and injection methods and was augmented to incorporate cycle-to-cycle combustion variations. State estimation is performed by means of a Kalman filter which feeds into a model predictive controller. The model predictive controller chooses control actions based on a predefined cost function under consideration of bounds reflecting physical constraints. The engine model was utilized to establish a state-space model that serves the Kalman filter and model predictive controller for estimation and prediction. A comparative study investigating control actions and engine behavior was performed with and without limiting in-cylinder peak pressure as well as combustion noise, which is of particular interest for early pilot injection strategies. In addition, the proposed control architecture was investigated at two different levels of cycle-to-cycle variations and compared to the performance of a control structure with input disturbance rejection. For increased cycle-to-cycle variations, disturbance estimation reduces state fluctuations and control effort. In general, this investigation highlights control aspects specific to a compression-ignited combustion regime with low-reactivity fuel. The control algorithm is able to maintain the desired references for brake mean effective pressure and combustion phasing while controlling peak in-cylinder pressure and combustion noise.

42 ENGINEERING↗
DOE OSTI2025

Development and preliminary validation of a mechanistic multiscale model for fuel-cladding chemical interaction in metallic nuclear fuels

Despite decades of fuel rod material and design improvements, fuel-cladding chemical interaction (FCCI) remains the single-most lifetime-limiting behavior for modern metallic fuel rods. Constraining fuel lifetime increases operating costs, limiting the economic viability of commercializing metallic nuclear fuel technology. A mechanistic multiscale model utilizing the finite element method-based MARMOT and BISON codes was developed to more confidently predict cladding-side FCCI and its impact on fuel performance. The new BISON model incorporates mesoscale models for the effects of fuel microstructure evolution on the transport of wastage-inducing lanthanides through the fuel and for the kinetics of cladding wastage layer growth. The mesoscale models, in turn, build on lanthanide transport property data obtained from the atomistic scale. Preliminary validation studies using wastage thickness and cladding profilometry data from four fuel rods irradiated in Experimental Breeder Reactor II experiment X447 and one fuel rod from Fast Flux Test Facility experiment IFR1 show that the new model predicts cladding wastage and its effects on cladding deformation as well as existing empirical FCCI correlations. The new model is expected to aid in the design of new metallic fuel concepts, including fuel additives, cladding liners, and sodium-free annular fuel geometries. In conclusion, future work will focus on broader validation and refinement of the model’s treatment of different fuel alloys and cladding materials.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Status Report on Design of In-situ Thermomechanical Testing at LANSCE

Nuclear fuel encounters severe thermomechanical environments in which its mechanical response is determined by its microstructure, temperature and stress level histories. Simulating the response of such microstructures is crucial for predicting both performance and transient fuel mechanical responses and experimental verification of such predictions is therefore of great interest. While most of the deformation in a nuclear fuel rod occurs in the cladding, deformation of the fuel itself is still of interest with deformation mechanisms at operating temperature and above including creep, swelling, cracking as well as pellet-clad interaction. Characterization of these properties and understanding of the underlying deformation phenomena at operating or excursion temperatures is therefore of great importance for development and ultimately licensing of improved and novel nuclear fuel forms. Diffraction techniques offer unique insight on the atomistic (e.g. crystal structure) and microstructure (e.g. phase transformations, texture, defects) length scales and have a long history of providing unique data to inform relevant deformation models that enable the required predictive capabilities. For example, dislocations lead to diffraction peak broadening that can be characterized to estimate the dislocation density and study the role of dislocations on the deformation while measuring lattice strains allows to studie load sharing in two phase materials. In this report the requirements for a sample environment for high temperature deformation of nuclear fuels are defined. The HIPPO neutron time-of-flight diffractometer at LANSCE will host this sample environment and is also described. This instrument covers diffraction angles from 140° to 40° and is also equipped with an event-mode neutron imaging detector system, enabling energy-resolved neutron imaging in parallel with the diffraction that could measure sample temperature from Doppler broadening of neutron absorption resonances or measure pore densities from changes in the attenuation. Designs of devices to characterize thermomechanical properties of nuclear fuel without diffraction are also considered to guide the design. While this report is focused on applications for nuclear fuels, the device can also characterize cladding, moderator or structural materials and therefore contribute to other fields of research and development for advanced reactors. The temperatures planned to be reached are above 2000℃, thus enabling characterization of LWR reactor fuels under accident scenarios but also reaching temperatures of fuels developed for nuclear thermal propulsion and providing opportunities to characterize those. In conjunction with the energy-resolved neutron imaging detector, this setup would allow to measure neutron cross-sections at high temperatures, filling a gap towards development of reactors operating at high temperatures.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Evaluating Diesel/Biofuel Blends Using Artificial Neural Networks and Linear/Nonlinear Equations

Abstract The use of biomass-derived additives in diesel fuel mixtures has the potential to increase the fuel’s efficiency, decrease the formation of particulate matter during its combustion, and retain the fuel’s behavior in cold weather. To this end, identifying compounds that enable these behaviors is paramount. The present work utilizes a series of linear and non-linear equations in series with artificial neural networks to predict the cetane number, yield sooting index, kinematic viscosity, cloud point, and lower heating value of multi-component blends. Property values of pure components are predicted using artificial neural networks trained with existing experimental data, and these predictions and their expected errors are propagated through linear and non-linear equations to obtain property predictions for multi-component blends. Individual component property prediction errors, defined by blind prediction median absolute error, are 4.91 units, 7.84 units, 0.06 cSt, 4.00 °C, and 0.55 MJ/kg for cetane number, yield sooting index, kinematic viscosity, cloud point, and lower heating value respectively. On average, property predictions for blends are shown to be accurate to within 6% of the blends’ experimental values. Further, a multitude of compounds expected to be produced from catalytically upgrading products of fast pyrolysis are evaluated with respect to their behavior in diesel fuel blends.

09 BIOMASS FUELS↗
DOE OSTI2022

Short- to Intermediate-Range Structure, Transport, and Thermophysical Properties of LiF–NaF–ZrF4 Molten Salts

LiF–NaF–ZrF 4 multicomponent molten salts are identified as promising candidates for coolant salts in molten salt reactors and advanced high-temperature reactors. This study focused on low-melting point salt compositions of interest: 38LiF–51NaF–11ZrF 4 , 42LiF–29NaF–29ZrF 4 , and 26LiF–37NaF–37ZrF 4 . Ab-initio molecular dynamics (AIMD) calculations were performed and compared with available experimental data to assess the ability of rigid ion models (RIM) to reproduce short- to intermediate-range structure, transport, and thermophysical properties of the LiF–NaF–ZrF 4 salt mixtures. It is found that as ZrF 4 mol% increases, the average cation–anion coordination number (CN) of monovalent cations (Li + , Na + ) obtained from RIM calculations decreases, while multivalent Zr 4+ CN varied from 15% to 19% in comparison to corresponding AIMD values. In addition, RIM is found to predict the existence of 7, 8, and 9 coordinated fluorozirconate complexes, while AIMD and the available experimental data showed an occurrence of 6, 7, and 8 coordinated complexes in the melt. The intermediate-range structure analysis revealed that while the RIM parameters are able to reproduce a local structure for lower ZrF 4 mol% salts such as in 38LiF–51NaF–11ZrF 4 , an extensive fluorozirconate network formation is observed in RIM simulations for higher ZrF 4 mol% compositions. The network generated by RIM parameters is found to be mainly connected by “corner-sharing” fluorozirconate complexes as opposed to both “edge-sharing” and “corner-sharing” connectively portrayed by AIMD. It is found that a close agreement between AIMD and the RIM salt structure for the 11-mol% ZrF 4 salt resulted in good agreement in the calculated Zr diffusivities and the viscosity values. However, due to the inaccurate short- to intermediate-range structure prediction by RIM for higher ZrF 4 mol% compositions, thermophysical properties such as densities and heat capacity differ by up to 26% and 27%, respectively, upon comparison with AIMD and experimental values. Also, the network-dominated properties such as diffusion coefficients and viscosities differed by up to two and three orders of magnitude, respectively. This study signifies the importance of accurate salt structure generation for an accurate prediction of transport and thermophysical properties of multicomponent molten salts.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Dynamic System Scaling Application to Accelerated Nuclear Fuel Testing

The development of nuclear fuel and materials requires a continuous effort to investigate the acute and prolonged effects of irradiation, thermal-material stress, chemical change, or other conceivable damage mechanics acquired during normal operation or accident scenarios throughout its lifetime. As in-core fuel property measurement techniques advance to support in the real-time, non-invasive, and enhanced accuracy realm, it is the future of fuel development to pursue to higher degrees of control, predictability of integral test behavior via separate effects test (SET), and shorter test time intervals. The fuel development life cycle from initial concept to commercial licensing is approximated to be 20 years and current literature suggests by optimizing fuel performance codes with SETs, the process could possibly be compressed to 5 to 10 years. Recently, a research group in the Idaho National Laboratory (INL) is testing reduced scale fuel rods and increased power density to accelerate evolution of fuel phenomena in metallic fuels. In support of nuclear fuel rod development, compressing fuel test process, and accelerating fuel phenomena, it was the purpose of this study to investigate nuclear fuel performance phenomena via literature review and effectively scale the initial conditions, boundary conditions, and geometric properties to describe the time-dependent response including to fuel burnup, thermal and mechanical stress, transmutations and inter-diffusion, and other relevant observed phenomena. The study was based on BISON simulations of historic EBR-II metallic fuel experiments and Dynamical System Scaling (DSS) method are utilized to assess effects of scaling activity including fuel phenomena acceleration and calculations of time-dependent distortions. The research successfully scaled metallic fuel phenomena, accelerated fuel testing, and assessed the distortions for each scaled case derived.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Assessment of uranium nitride interatomic potentials

Uranium mononitride (UN) is a promising nuclear fuel due to its high fissile density, high thermal conductivity, and suitability for reprocessing. In this study, two uranium nitride interatomic potentials are assessed: Tseplyaev and Starikov's angular-dependent potential and Kocevski et al.'s embedded atom model potential. Predictions of the thermophysical and elastic properties of UN, UN 2 , and α- and β-U 2 N 3 computed using both potentials are assessed and compared to available experimental data. Notably, the Tseplyaev potential performs better with the energetic aspects of UN, e.g., specific heat capacity and point defect formation energies, whereas the Kocevski potential performs better with the structural aspects of UN, e.g., thermal expansion as well as with the elastic properties. The reasons why the Kocevski potential underestimates the UN specific heat are explained by examining the UN phonon properties modeled using both potentials. The Kocevski potential shows better identification of the mechanical stability ranges of UN, UN 2 , and α- and β-U 2 N 3 , reasonably predicting the melting point of UN and predicting stable structures for UN 2 and α- and β-U 2 N 3 . On the other hand, the Tseplyaev potential predicts a premature phase change of both UN and UN 2 and cannot stabilize α- nor β-U 2 N 3 . However, the Kocevski potential cannot predict a stable α-U phase and is thus not suitable for the calculation of formation energies for non-stoichiometric point defects.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

High-Pressure Apparatus for Monitoring Solid–Liquid Phase Transitions

This work presents a new technique for observing the solid–liquid phase transformations in complex diesel fuel blends and diesel surrogates under high-pressure conditions intended to simulate those occurring in vehicle fuel injectors. A high-pressure apparatus based on a visual identification of freezing and thawing has been designed and built to monitor phase behavior and determine the crystallization temperature of complex fuels to predict wax precipitation. The proposed methodology was validated using pure substances—n-hexadecane (C16H34), cyclohexane (C6H12), and a binary cyclohexane/n-hexadecane mixture—all of which have been well-characterized previously. The crystallization temperatures of these compounds were measured from atmospheric pressure to 400 MPa for temperatures varying from 290 to 363 K and compared to those reported in the literature. The standard error of the estimated temperatures, based on a given pressure, between the experimental data obtained in this work was compared to data in the literature from Domanska et al. This methodology is being extended to investigate the properties of more complex fuel mixtures.

High-Pressure, phase change, diesel, solidificatio↗
DOE OSTI2021

High Performance PEFC Electrode Structures

Raytheon Technologies Research Center (RTRC) in collaboration with Ion Power Inc. and the University of Arkansas at Little Rock and the FC-PAD consortium executed a project that: (1) developed comprehensive models of oxygen and proton transport in cathode catalyst layers that provided new insights about losses associated with reducing platinum loading to 0.1 mg/cm 2 , (2) designed, fabricated, characterized, and tested cells featuring carbon-supported catalyst layers with low platinum loadings that met efficiency and power density targets set forth by the Hydrogen and Fuel Cells Technology Office, and (3) designed, fabricated, characterized, and tested novel catalyst structures. The modeling efforts produced eight scientific publications that shared new descriptions of oxygen and proton transport in catalyst layers with the broader technical community. The key physical insights from the models are expressed in straightforward algebraic expressions, facilitating easy adoption by other researchers. The models developed during the project explain why resistance to oxygen transport is inversely proportional to interfacial area of platinum – an observation that eluded quantitative explanation for a decade. The models quantitatively predict polarization when measured morphological parameters and transport properties are input. Membrane-electrode assemblies made during the project met stated goals for current at high voltage and power density at rated voltage for transportation fuel cells at the end of budget period 2. High performance was achieved by developing a diffusion layer with low resistance to oxygen transport and a PtCo/C cathode catalyst layer capable of supporting rapid proton and oxygen transport. The project was subsequently directed to curtail experimental work and focus on modeling during budget period 3. Ion Power Inc. developed new expertise in fabricating catalyst layers containing platinum-cobalt catalysts, and processing very thin membranes. The University of Arkansas at Little Rock developed novel nanocolumnar self-supported thin film electrocatalysts that can be grown on carbon supports using a simple high pressure sputter deposition method. The nanocolumnar microstructure provides adequate surface-to-volume ratio for efficient platinum utilization, and the conformal platinum shell with larger crystal grain sizes covering the carbon support surface may eliminate durability issues associated with catalyst dissolution, agglomeration, and carbon corrosion.

08 HYDROGEN↗
DOE OSTI2022

Predicting Catalyst Surface Stability Under Reaction Conditions Using Deep Reinforcement Learning and Machine Learning Potentials

Catalysts are critical for most large-scale energy intensive chemical transformation processes, such as energy storage, liquid fuel production, and the formation of chemical building blocks. The catalyst composition, structure, and morphology impact the performance under reaction conditions and influences properties like activity and selectivity. Many industrial catalysts offer imperfect activity/selectivity, or contain expensive metals. Further, catalysts can deactivate over time as the morphology changes or harsh reactive environments alter the surface structure. Methods to automatically model and predict the kinetics of how catalyst surfaces will restructure would enable engineers to design around these challenges, improve performance, increase longevity. This project investigated a specific type of machine learning model, deep reinforcement learning, machine learning models to act as surrogates for the physical system, and compared the results of these approaches with a specialized high-throughput experimental synthesis and measurement platform.

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

Factors Governing Oxygen Vacancy Formation in Oxide Perovskites

The control of oxygen vacancy (V O ) formation is critical to advancing multiple metal-oxide-perovskite-based technologies. In this work, we report the construction of a compact linear model for the neutral V O formation energy in ABO 3 perovskites that reproduces, with reasonable fidelity, Hubbard-U-corrected density functional theory calculations based on the state-of-the-art, strongly constrained and appropriately normed exchange-correlation functional. We obtain a mean absolute error of 0.45 eV for perovskites stable at 298 K, an accuracy that holds across a large, electronically diverse set of ABO 3 perovskites. Our model considers perovskites containing alkaline-earth metals (Ca, Sr, and Ba) and lanthanides (La and Ce) on the A-site and 3d transition metals (Ti, V, Cr, Mn, Fe, Co, and Ni) on the B-site in six different crystal systems (cubic, tetragonal, orthorhombic, hexagonal, rhombohedral, and monoclinic) common to perovskites. Physically intuitive metrics easily extracted from existing experimental thermochemical data or via inexpensive quantum mechanical calculations, including crystal bond dissociation energies and (solid phase) reduction potentials, are key components of the model. Beyond validation of the model against known experimental trends in materials used in solid oxide fuel cells, the model yields new candidate perovskites not contained in our training data set, such as (Bi,Y)(Fe,Co)O 3 , which we predict may have favorable thermochemical water-splitting properties. The confluence of sufficient accuracy, efficiency, and interpretability afforded by our model not only facilitates high-throughput computational screening for any application that requires the precise control of V O concentrations but also provides a clear picture of the dominant physics governing V O formation in metal-oxide perovskites.

08 HYDROGEN↗