Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “oxygen vacancy”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Defect Thermodynamics and Transport Properties of Proton Conducting Oxide BaZr 1–x Y x O 3–δ (x ≤ 0.1) Guided by Density Functional Theory Modeling

Density functional theory-based thermodynamic modeling was performed to determine the effect of humidity and H 2 /O 2 gas pressure on the defect chemistry and transport properties of proton conducting oxide BaZr 1-x Y x O 3-δ (x ≤ 0.1) in the temperature range of 800–1200K, relevant for solid oxide fuel/electrolysis cell applications. The first-principles charge defect analysis was carried out to obtain the defect energetic ensembles as a function of Fermi level as well as the dependence of the hydration and oxidation reaction energies on complex defect-dopant configurations. It is shown that oxygen vacancies introduced to compensate for the extra Ba vacancies may cause a slight increase in the concentration of hydroxyl and proton species upon hydration of Ba deficient BaZr 1-x Y x O 3-δ , while the overall proton diffusivity is predicted to decrease with increased proton hopping barriers upon trapping near the Ba vacancies. The developed defect model is further demonstrated to be able to describe the bulk defect chemistry and transport properties of BaZr 0.9 Y 0.1 O 3-δ under the solid oxide cell operating conditions. Finally, the theoretical framework developed in this work further allows the inclusion of configurational, energetic and electronic characteristics of the defects that can be used as a complementary tool to experimental measurements for testing various mechanistic pathways or to provide essential mechanistic data.

36 MATERIALS SCIENCE↗
DOE OSTI2022

MOF-Transformed In 2 O 3-x @C Nanocorn Electrocatalyst for Efficient CO 2 Reduction to HCOOH

For electrochemical CO 2 reduction to HCOOH, an ongoing challenge is to design energy efficient electrocatalysts that can deliver a high HCOOH current density (J HCOOH ) at a low overpotential. Indium oxide is good HCOOH production catalyst but with low conductivity. In this work, we report a unique corn design of In 2 O 3-x @C nanocatalyst, wherein In 2 O 3-x nanocube as the fine grains dispersed uniformly on the carbon nanorod cob, resulting in the enhanced conductivity. Excellent performance is achieved with 84% Faradaic efficiency (FE) and 11 mA cm –2 J HCOOH at a low potential of – 0.4 V versus RHE. At the current density of 100 mA cm –2 , the applied potential remained stable for more than 120 h with the FE above 90%. Density functional theory calculations reveal that the abundant oxygen vacancy in In 2 O 3-x has exposed more In 3+ sites with activated electroactivity, which facilitates the formation of HCOO* intermediate. Operando X-ray absorption spectroscopy also confirms In 3+ as the active site and the key intermediate of HCOO* during the process of CO 2 reduction to HCOOH.

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

Ink-based laser powder bed fusion of barium titanate

Barium titanate (BTO) is a lead-free functional ceramic that is widely used in sensors, transducers, and actuators. The increasing demand for complex geometries and quick design iterations motivates the use of additive manufacturing techniques. Laser powder bed fusion (L-PBF) additive manufacturing technique enables the net-shape or near-net-shape fabrication of metals and alloys. Ceramic L-PBF remains challenging due to their high melting temperature. Here, to address this challenge, a nanoparticle ink feedstock was incorporated into the L-PBF process to allow low power melting of BTO. Key print parameters (i.e. laser power, scanning speed, and hatching spacing) were optimized to fabricate millimeter-scale BTO parts. The printed BTO samples possessed an inhomogeneous equiaxed-columnar-cellular microstructure. Semiconducting current-voltage behavior was observed in as-printed samples, which was attributed to the creation of oxygen vacancies. The piezoelectric d 33 coefficient was measured to be 4.6 pC/N.

Dielectrics↗
DOE OSTI2023

Silica modulated palladium catalyst with superior activity for the selective catalytic reduction of nitrogen oxides with hydrogen

A significant promotion effect of colloidal SiO 2 on Pd/TiO 2 catalyst for improving the activity and N 2 selectivity in the selective catalytic reduction of NO with H 2 (H 2 -SCR) was reported, both before and after hydrothermal aging. The SiO 2 addition not only benefited the formation of more oxygen vacancies but also increased Pd dispersion and created rich Pd-SiO 2 interfaces on Pd-SiO 2 /TiO 2 catalyst. The H 2 -SCR of NO on both Pd/TiO 2 and Pd-SiO 2 /TiO 2 followed a Langmuir-Hinshelwood mechanism, in which the adsorbed bridging bidentate nitrates on TiO 2 could react with dissociated H 2 species on Pd sites at Pd-TiO 2 interface. Particularly, on the promoted Pd-SiO 2 /TiO 2 , the Pd-TiO 2 interface could facilitate NO adsorption and activation, and the created Pd-SiO 2 interface could benefit H 2 adsorption and dissociation, thus contributing to its significantly improved H 2 -SCR activity. Additionally, the active Pd species (which should be Pd0 under reaction condition with higher surface concentration) within Pd-SiO 2 /TiO 2 could be well stabilized during long-term hydrothermal aging process through the formation of rich Pd-SiO 2 interfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2023

LaCeO x coupled N-doped graphene/Ru single-atoms as a binary-site catalyst for efficient hydrogen evolution based on hydrogen spillover

Based on the hydrogen spillover effect, versatile binary-site catalysts are promising for the alkaline hydrogen evolution reaction (HER). In which one site can adsorb and dissociate water, while the neighboring site is favorable for liberating hydrogen. Inspired by these possibilities, oxygen vacancy (O v )-rich LaCeO x coupled N-doped graphene/Ru single-atoms (LaCeO x @NGr/Ru 1 ) was synthesized as an efficient dual-site HER catalyst. It delivered an impressive low overpotential of 22 mV at a current density of 10mAcm –2 and a slight Tafel slope of 40 mV dec –1 in an alkaline medium, outstanding the advanced Ru-based catalysts. Moreover, this promising binary-component catalyst exhibited higher mass activity and longer lasting durability than commercial Pt/C catalyst (20 wt%). Finally, experimental and theoretical investigations provided insights into the HER mechanisms of LaCeO x @NGr/Ru 1 based on three indispensable steps: water adsorption and dissociation on O v -rich LaCeO x , diffusion of generated H* species towards Ce 3+ -N-Ru 1 bridges, and hydrogen evolution on Ru 1 sites.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Co-existence of atomically dispersed Ru and Ce 3+ sites is responsible for excellent low temperature N 2 O reduction activity of Ru/CeO 2

Nitrous oxide N 2 O reduction is a big challenge due to high global warming potential of N 2 O. (~300 times higher compared with CO 2 ). The best known catalysts, such as Rh/ceria, require relatively high temperatures for N 2 O decomposition. Herein, we report that Ru/ceria catalysts with low Ru loading of ~0.25 wt% efficiently catalyze low temperature N 2 O reduction by CO starting at 100 °C (full N 2 O conversion below 200 °C) under industrially relevant flow rates and gas concentrations. Further, this remarkable performance stems from maintaining isolated Ru cations even on reduced ceria surface and, simultaneously, the propensity of Ru to affect ceria surface to form labile surface oxygen thereby creating large number of oxygen vacancies (Ce +3 cations) in the presence of CO. In contrast, for Rh/CeO 2 catalysts with equivalent metal loading, the activity is much lower because atomically dispersed Rh sinters into metallic clusters at the onset reaction temperature (~200 °C): these clusters are much less effective than isolated single Ru ions, with lower Ce +3 concentration maintained on reduced Rh/CeO 2 catalyst. Our study highlights the benefits of gaining molecular-level insight into the dynamic nature of catalytically active sites under reaction conditions for preparing catalysts containing low loading of precious metals with unsurpassed low temperature activity.

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

Exsolution of NiCo alloys over Ruddlesden-Popper perovskite for mild electrochemical synthesis of ammonia on protonic ceramic electrochemical cells

Ammonia synthesis from renewable energies on protonic ceramic electrochemical cells (PCECs) shows great potential. The primary challenges in ammonia synthesis on PCECs include sluggish catalytic activity, competition from the hydrogen evolution reaction, and unsatisfactory durability of the cathode, which is the active site of ammonia generation. Here, in this study, we report an elaborate design of the cathode with an intended formula of Pr 4 Ni 1.79 Co 1.2R u 0.01 O 10-δ , where NiCo alloy nanoparticles are exsolved from Ruddlesden-Popper perovskite substrates after the reduction in 5 % H 2 /Ar at 400 °C for 1 h, for electrocatalysis of the nitrogen reduction reaction to ammonia. The host material Pr 4 Ni 1.8 Co 1.2 O 10-δ with embeddable layered structure and stability was deliberately chosen to fix the Ru cation and maximize the catalytic activity of NiCo. Also, density functional theory calculations suggest that Ru doping provides an optimal balance between structural stability and redox activity, facilitating the controlled exsolution of NiCo nanoparticles and enhancing catalytic performance. As a result, the composite electrode with exsolved NiCo alloy and abundant oxygen vacancies on fuel-electrode-supported PCECs achieves a superior electrochemical activity towards ammonia synthesis: a peak ammonia formation rate of 27.84 μg h −1 cm −2 and excellent Faradaic efficiencies of 62.6 % at 350 °C.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2022

Transition-metal doped titanate nanowire photocatalysts boosted by selective ion-exchange induced defect engineering

Defect engineering through elemental doping is an efficient way to boost the performance of semiconductor photocatalysts. For this work, transition-metal (TM) doped titanate nanowires (TNWs) were prepared via ion-exchange over the titanate precursors and demonstrated for the Rhodamine B (RhB) degradation under ultraviolet (UV) light irradiation. The ion-exchange of selective ions (V 5+ , Cr 3+ , Ni 2+ , and Zn 2+ ) with protons from pristine TNWs resulted in the hierarchical meso-porosity of nanowires with large pores of ~5–20 nm by TM doping and small pores of ~3.6–4.5 nm inherited from pristine TNWs, which facilitates the mass transfer while maintaining high surface area and active sites. Meanwhile, the TM intercalation partially reduces the Ti 4+ to Ti 3+ and narrows the optical bandgap, which, together with oxygen vacancies and superoxide radicals from pristine TNWs, enhance the adsorption and photocatalytic degradation performance of RhB. This work helps to elucidate the effects of transition-metal doping and provides a rational strategy towards high performance titanate-based photocatalysts for efficient and sustainable wastewater treatment.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Multiple redox mechanisms for water-gas shift reaction on Fe 3 O 4 (1 1 1) surface: A density functional theory and mean-field microkinetic modeling study

In this work, we investigated redox mechanisms for the high temperature (HT) water–gas shift reaction (WGSR) on the Fe 3 O 4 (1 1 1) surface using density functional theory (DFT) calculations and mean-field microkinetic modeling (MKM). The redox pathways branched into three Langmuir-Hinshelwood processes (LH1, LH2 and LH3) and one Mars-van-Krevelen (MVK) process (in the presence of oxygen vacancy) from multidentate binding of CO and CO 2 over four available reactive sites (Fe oct2 , Fe tet1 , Fe bridge, and O1). We found that the LH1 and LH2 processes have CO adsorption at a single iron site (Fe oct2 or Fe tet1 ), while the LH3 and MVK processes have stronger chemisorption of CO or CO 2 by both the Fe oct2 site and the O 1 site. From the mechanistic study of these reaction path, we recognized that availability of O 1 sites was key to proceed to either CO oxidation by the single site (LH1, LH2) or by the dual site (LH3, MVK). We observed that the reaction energetics were significantly different in the CO oxidation steps, where the single or dual site results in markedly different apparent activation energy (73–281 kJ/mol) and reaction rates (10 -3 –10 -10 mol∙m -2 ∙s -1 ) among the four reaction mechanisms. The utilization of mean field MKM with DFT reaction energetics helps explain the experimental debates for the catalytic reaction details as well as providing a possible direction to engineer the catalyst with higher activity.

42 ENGINEERING↗
DOE OSTI2025

Understanding how defects and dopant atoms in copper surface oxides affect reactivity

Copper and its oxides are key catalytic materials, on which reactions often occur at the metal/oxide interface. Here, in this work, we directly connect the induction period observed during methanol-driven reduction of thin-film copper oxides to their atomic-scale structural order. Using temperature-programmed desorption (TPD) methanol titrations combined with scanning tunneling microscopy, we show that highly ordered oxide phases – particularly the “29” structure with its low defect density – exhibit long induction periods and initially low reactivity. The induction period, defined as the number of methanol TPD cycles required to reach half of the maximum formaldehyde yield, scales with oxide order and oxygen coverage. Enhanced reactivity of well-ordered oxides emerges only after repeated methanol adsorption/desorption cycles generate oxygen vacancies and new Cu(111)/Cu x O interfacial sites. In contrast, disordered or sub-stoichiometric oxides, which contain more intrinsic defects and interfaces, are active from the first TPD cycle. We further examine how dilute Pt and Rh dopants influence oxide order and reactivity: 1% Pt increases defect density and catalytic activity, while 1% Rh promotes oxide ordering and longer induction periods. These findings demonstrate that dilute alloying provides a potential method for tuning the structure and reactivity of Cu(111)/Cu x O interfaces.

Cu(111)Methanol oxidation↗
DOE OSTI2020

Enhanced tunability of two-dimensional electron gas on SrTiO 3 through heterostructuring

Two-dimensional electron gases (2DEGs) on the SrTiO 3 (STO) surface or in STO-based heterostructures have exhibited many intriguing phenomena, which are strongly dependent on the 2DEG-carrier density. In this paper, we report that the tunability of the 2DEG-carrier density is significantly enhanced by adding a monolayer LaTiO 3 (LTO) onto the STO. Ultraviolet (UV) irradiation induced maximum carrier density of the 2DEG in LTO/STO is increased by a factor of ~4 times, compared to that of the bare STO. By oxygen gas exposure, it becomes 10 times smaller than that of the bare STO. This enhanced tunability is attributed to the drastic surface property change of a polar LTO layer by UV irradiation and O 2 exposure. This indicates that the 2DEG controllability in LTO/STO is more reliable than that on the bare STO driven by defects, such an oxygen vacancy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

A Graph Dynamical neural network approach for decoding dynamical states in ferroelectrics.

Ferroelectric materials such as BaTiO 3 show tremendous potential for emerging advances in memory devices, particular neuromorphic type devices. High density of memory can be obtained by stabilising polar domain walls at the nanoscale, regions of discontinuity between the well-defined polarization order parameter, but little is known about what controls their structure and dynamics in real nanoscale materials. Indeed, chiral polar domain walls have been observed in heterogeneous ferroelectrics, such as oxygen-deficient BaTiO 3 , but very little is known about how such polar-domains walls interact with defects. Indeed, a critical understanding of how dynamics of domain-walls depend on point-defects is crucial to create engineered ferroelectric memory devices. For this work, we perform large-scale simulations of nansocale domain-wall dynamics in pristine and defective BaTiO 3 using reactive force-field developed by us earlier (Phys. Chem. Chem. Phys., 2019, 21, 18240–18249), and capture their dynamical dependence on point defects using a graph dynamical neural-network approach, which we adapted to interrogate solids with well-defined order-parameters, and implemented using Pytorch based libraries. Our machine learning (ML) approach goes beyond the traditional post-processing methods to capture both spatial and temporal heterogeneities of large-scale molecular dynamics simulations of complex defective ferroelectric oxide materials. We crucially find that isolated oxygen vacancies introduce very localized spatial regions (~1–2 unit-cell in length) that show slow dipole relaxation due to formation of defect-dipoles, and that these defect-dipoles in turn slow the intrinsic dynamics of domain walls. Further, the roughness of domain walls, also influenced by vacancies, introduce dynamic heterogeneity along the domain-wall. As such we find a novel mechanism by which quenched disorder due to defects introduce dynamic heterogeneity thereby influencing response to external fields (particularly time varying fields) in a ferroelectric. Our study also emphasizes the need for creating digital twins of dynamical quantities to achieve autonomous in operando control of nanoscale switching.

42 ENGINEERING↗
DOE OSTI2022

Effect of surface acidity modulation on Pt/Al 2 O 3 single atom catalyst for carbon monoxide oxidation and methanol decomposition

With ideal metal utilization efficiency and homogeneous active sites, single atom catalysts (SACs) have attracted extensive attention in heterogeneous catalysis field. The performance of SACs is highly dependent on the local coordination environment of metal single atoms. However, few works reported the effect of surface acidity modification on the performance of SACs such as Pt SACs which have been widely studied nowadays. In this work, a Pt/Al 2 O 3 SAC has been successfully prepared using a bayerite Al(OH) 3 as support through simple IWI method for Pt loading. Through post-modification of Pt/Al 2 O 3 with acidic WO 3 and basic MgO, the surface acidity of Pt/Al 2 O 3 was tuned. Here, activity of CO oxidation and methanol decomposition on the modified catalysts has been evaluated. It was found that the catalytic performance was highly related to the surface acidity of Pt SACs. Basic MgO modified Pt/Al 2 O 3 performed better in CO oxidation under the reaction condition with H 2 O. Acidic WO 3 modified Pt/Al 2 O 3 showed superior activity in methanol decomposition. Detailed characterization and in situ diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS) study suggested that the positive H 2 O effect on CO oxidation was due to the mitigated formation of bicarbonate species on Pt/Al 2 O 3 and WO 3 /Pt/Al 2 O 3 catalysts and the promoted desorption of carbonate species from MgO/Pt/Al 2 O 3 . The superior methanol decomposition performance on WO 3 /Pt/Al 2 O 3 was owing to its abundant oxygen vacancies and high resistance to carbonate accumulation. This study highlights the importance of surface acidity modulation of Pt SACs in heterogeneous catalysis and suggests that fine-tuning the surface acidity can be an effective strategy in new SACs development.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Mechanistic understanding of support effect on the activity and selectivity of indium oxide catalysts for CO 2 hydrogenation

Herein we present a mechanistic study on the support effect (ZrO 2 and CeO 2 ) of In 2 O 3 catalysts in CO 2 hydrogenation by a combined experimental and computational approach. Kinetic experiments and surface characterization suggested that the activity of In 2 O 3 catalysts cannot be simply correlated with the abundance of surface oxygen vacancies (O v ) formed by either H 2 -reduction or thermal treatment, which has been frequently invoked in previous studies. The support effect should originate from the electronic interactions between In 2 O 3 and the support oxide, rather than geometric factors or the difference in the particle size of In 2 O 3 . Theoretical modelling revealed that surface O v facilitate the formation and stabilization of the formate (HCOO*) intermediate. While a carbonate-like structure is favored for CO 2 adsorption on CeO 2 -supported or unsupported In 2 O 3 catalysts, CO 2 tends to bind strongly in a bent configuration on the O v site at the In 2 O 3 -ZrO 2 interface. The distinct CO 2 adsorption structures on different supported In 2 O 3 catalysts may account for the different reaction energy profiles in the subsequent hydrogenation reactions, especially the rate-limiting step, i.e., hydrogenation of HCOO* to CH 2 O* and methoxy (CH 3 O*). The relatively higher methanol selectivity of In 2 O 3 catalyst supported on ZrO 2 with respect to that on CeO 2 are suggested to stem from the greater energy difference (Δ$E_a$) between the parallel hydrogenation and C-O bond cleavage of HCOO*, which leads to the formation of methanol and CO, respectively. We report this study underlines the important role of metal-oxide-interface in determining the catalytic behavior of oxide-supported In 2 O 3 catalysts in CO 2 conversion.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2023

Optimization and understanding of ZnO nanoarray supported Cu-ZnO-Al 2 O 3 catalyst for enhanced CO 2 -methanol conversion at low temperature and pressure

Cu-ZnO-Al 2 O 3 is the most widely applied catalyst for CO 2 hydrogenation to methanol. However, it is still a challenge to produce methanol using this catalyst under low-temperature (<250 °C) and low-pressure (<10 bar) conditions with desirable yield and selectivity. In this work, by tuning the experimental processing parameters such as solvent, loading amount, and annealing temperature, highly improved ZnO nanoarray supported Cu-ZnO-Al 2 O 3 catalysts have been successfully demonstrated. Here, by using organic solvent (N,N-dimethylformamide (DMF), acetone, or isopropanol) for dip-coating loading process instead of deionized (DI) water, Cu-ZnO-Al 2 O 3 nanocatalysts was comparatively better dispersed on the nanorod array support with populated and abundant active sites, thus enhancing the methanol yield. With the control of the loading amount and annealing temperature, finely distributed Cu nanoparticles were obtained on the ZnO nanorod surfaces to enhance the interactions between Cu and ZnO nanorod surfaces. Further improvement of the catalyst performance is demonstrated by tuning the reaction space velocity. At 200 °C and 10 bar conditions, the optimized catalyst achieved a methanol yield of 6.46 mol h –1 kg –1 with 100 % selectivity. The good stability after prolonged testing of the catalysts demonstrates the potential practical implementation. The in situ diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS) measurements under the 1 bar reveal that the CO 2 hydrogenation to methanol on the ZnO nanoarray supported Cu-ZnO-Al 2 O 3 catalyst follows the CO reaction pathway, due to the surface oxygen vacancies on ZnO nanorods which facilitate CO 2 dissociation.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

Zero-strain high entropy spinel oxide (FeNiCuCrMn) 3 O 4 @rGO as high-performance anode for sodium ion battery

High entropy spinel oxides (HEOs), as a new type of anode material with at least five transition metal elements and high theoretical specific capacity, have enormous promise in sodium ion batteries (SIBs). However, the influence of different transition metal elements on the structure and electrochemical performances of HEOs is not clear. Here, in this study, the addition of Mn to HEO can increase the surface oxygen vacancy (O v ) concentration, shorten the band gap, and enhance the electrochemical performance of (FeNiCuCrMn) 3 O 4 @rGO. XPS results show (FeNiCuCrMn)3O4@rGO has the highest surface O v concentration. Ab initio calculations demonstrate that the indirect band gap of (FeNiCuCrMn)3O4@rGO is 0.068 eV. In situ X-ray diffraction and synchrotron X-ray absorption spectroscopy reveal that (FeNiCuCrMn) 3 O 4 @rGO exhibits a zero-strain characteristic with an inconsequential unit cell volume change of 0.2 %. Specifically, (FeNiCuCrMn) 3 O 4 @rGO as anode materials for SIBs show good cyclic stability (with a capacity retention of 84 % at 500 mA g −1 after 3000 cycles). Overall, this work provides a new direction for optimizing transition metal elements in HEOs.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Facet-dependent strong metal-support interactions control the C–O bond activation

Reducible metal oxides are selective and effective for C–O bond scission of hydrodeoxygenation reactions via the reverse Mars-van Krevelen mechanism. Creating efficient redox centers while minimizing metal surfaces leads to highly selective catalysts. Single noble metal atoms activate the surface M–O bond, but the catalyst activity is limited due to low loading. Here, we report that the strong metal-support interaction between Ir and CeO 2 is facet sensitive, and certain facets regulate the C–O bond cleavage. At 300°C reduction, Ir is mostly encapsulated on an octahedron by (111) facets but remains exposed by (110) facets. The former is selective, whereas the latter is not. In conclusion, density functional theory indicates that Ir encapsulation is favored on (111) under reaction conditions, and oxygen vacancies more readily form on encapsulated Ir than on pristine ceria.

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

Multisublattice cluster expansion study of short-range ordering in iron-substituted strontium titanate

Owing to the challenges in obtaining realistic atomic configurations in large chemical phase spaces, it is not straightforward to describe structure–property relations in materials exhibiting configurational disorder. One example is iron-substituted strontium titanate (SrTi 1–x Fe x O 3–d , STF), a promising perovskite-derivative cathode material in solid oxide fuel cells that exhibits full solid solubility 0 ≤ x ≤ 1 and a tendency to exhibit short-range order. Here we demonstrate a multisublattice cluster expansion (CE) framework and apply it to STF across the full composition range. The CE approach is distinct from more traditional CE formulations in that clusters are defined explicitly by the chemical species distributed among multiple sublattices, rather than via cluster functions of occupation variables with decoration. The modified CE approach makes it easy to distinguish meaningful chemical interactions that are harder to extract from conventional CE, since for the latter chemical identity in a cluster is expressed as a product of site occupations. The least absolute shrinkage and selection operator (LASSO) is implemented as a regression analysis tool to select key clusters and avoid overfitting. We demonstrate this formulation on STF, and show that it can accurately predict configurational energies in comparison to conventional CE. From the key clusters, we identify that short-range ordering between substitutional Fe and oxygen vacancies (V O ) results in the formation of Fe–VO strings. In addition, we consider the stability of STF through CE-based Monte Carlo (MC) simulations and confirm the presence of superstructures that were previously observed in transmission electron microscopy. In this work, analysis of atomic configurations from MC samples reveals variations in the oxidation state of Fe atoms, which can be explained by the ordering tendency of Fe and V O . The cluster description and selection formalism described here may be applied to other disordered multisublattice systems for accurate and efficient material modeling.

36 MATERIALS SCIENCE↗