Synchrotron-based techniques for characterizing STCH water-splitting materials

Understanding the role of oxygen vacancy–induced atomic and electronic structural changes to complex metal oxides during water-splitting processes is paramount to advancing the field of solar thermochemical hydrogen production (STCH). The formulation and confirmation of a mechanism for these types of chemical reactions necessitate a multifaceted experimental approach, featuring advanced structural characterization methods. Synchrotron X-ray techniques are essential to the rapidly advancing field of STCH in part due to properties such as high brilliance, high coherence, and variable energy that provide sensitivity, resolution, and rapid data acquisition times required for the characterization of complex metal oxides during water-splitting cycles. X-ray diffraction (XRD) is commonly used for determining the structures and phase purity of new materials synthesized by solid-state techniques and monitoring the structural integrity of oxides during water-splitting processes (e.g., oxygen vacancy–induced lattice expansion). X-ray absorption spectroscopy (XAS) is an element-specific technique and is sensitive to local atomic and electronic changes encountered around metal coordination centers during redox. While in operando measurements are desirable, the experimental conditions required for such measurements (high temperatures, controlled oxygen partial pressures, and H2O) practically necessitate in situ measurements that do not meet all operating conditions or ex situ measurements. Here, we highlight the application of synchrotron X-ray scattering and spectroscopic techniques using both in situ and ex situ measurements, emphasizing the advantages and limitations of each method as they relate to water-splitting processes. The best practices are discussed for preparing quenched states of reduction and performing synchrotron measurements, which focus on XRD and XAS at soft (e.g., oxygen K-edge, transition metal L-edges, and lanthanide M-edges) and hard (e.g., transition metal K-edges and lanthanide L-edges) X-ray energies. The X-ray absorption spectra of these complex oxides are a convolution of multiple contributions with accurate interpretation being contingent on computational methods. The state-of-the-art methods are discussed that enable peak positions and intensities to be related to material electronic and structural properties. Through careful experimental design, these studies can elucidate complex structure–property relationships as they pertain to nonstoichiometric water splitting. A survey of modern approaches for the evaluation of water-splitting materials at synchrotron sources under various experimental conditions is provided, and available software for data analysis is discussed.


Introduction
Solar thermochemical hydrogen production, termed STCH, is an emerging process for thermochemical water-splitting technologies with little or no greenhouse gas emissions (Steinfeld, 2005;Rao and Dey, 2017;Lu et al., 2019).Out of the numerous STCH cycles described in the literature (Steinfeld, 2005;Perret, 2011;Muhich et al., 2015;Rao and Dey, 2017;Lu et al., 2019), the two-step redox-active metal oxide cycle is one of the more promising candidate approaches.The fundamental basis for the oxide cycle is simple-a metal oxide is reduced at a high temperature to create oxygen vacancies that are subsequently filled when the defected material is exposed to steam, splitting the H 2 O molecule, generating H 2 , and completing the cycle (Figure 1).Efficacious STCH water splitting necessitates a metal oxide that can undergo reduction via non-stoichiometry at sufficiently low temperatures, demonstrate fast reoxidation kinetics with H 2 present in the steam feed, and exhibit structural stability over the course of up to hundreds of thousands of redox cycles.Identifying an oxide material with the necessary chemical properties for robust water splitting has thus far eluded researchers and is a key challenge that must be overcome for STCH to be considered a viable technology for large-scale hydrogen production.The development of novel STCH materials is contingent on unraveling the relationship between a material's structural and electronic response to the enthalpy and entropy of vacancy formation, requiring empirical and modeling efforts working in tandem (Zinkevich et al., 2006;Lany, 2008;Deml et al., 2014;Hao et al., 2014;Bork et al., 2015;Deml et al., 2015;Muhich et al., 2015;Zhang et al., 2015;Emery et al., 2016;Tsvetkov et al., 2016;Naghavi et al., 2017;Sai Gautam et al., 2020).
Currently, ceria is the benchmark material for two-step metal-oxide-based STCH (Hao et al., 2014;Rothensteiner et al., 2015;Tanwar et al., 2016;Lu et al., 2019).Ceria demonstrates phase stability in the reduced state (Zinkevich et al., 2006;Chueh et al., 2012;Hao et al., 2014;Rothensteiner et al., 2015;Naghavi et al., 2017;Lu et al., 2019), fast reoxidation kinetics (McDaniel et al., 2013;Hao et al., 2014;Ackermann et al., 2015), and importantly, a high tolerance to hydrogen in the steam stream during reoxidation (Chueh et al., 2012;Hao et al., 2014;Barcellos et al., 2018).Unfortunately, temperatures in excess of 1,500 °C are required to reduce ceria appreciably, limiting practical viability (Barcellos et al., 2018;Cheng et al., 2021).Increasingly, researchers are turning their attention away from simple oxide materials to more complex material systems, which include redox-active nonstoichiometric oxides such as Sr x La 1−x Mn y Al 1−y O 3−δ (SLMA) (McDaniel et al., 2013;Suntivich et al., 2014) and BaCe 0.25 Mn 0.75 O 3−δ (BCM) (Barcellos et al., 2018).Here, δ is used to indicate the deviation in oxygen stoichiometry from the fully oxidized state.However, there are multiple ways to define or discuss δ such as normalizations by formula units of oxygen, cations, or volume (e.g., concentration and vacancies•cm −3 ).The perovskite-based oxides have demonstrated significant promise through an increased per-cycle yield at lower temperatures, governed by changes in oxygen vacancy concentration between the reduced and reoxidized states (Δδ), at given H 2 to H 2 O ratios.Continued progress in STCH materials' development requires a more complete understanding of the water-splitting mechanism.The impact of oxygen vacancies on the structural response of materials can be subtle and is often difficult to observe with conventional laboratory instruments, thus requiring synchrotron radiation-enabled techniques such as X-ray diffraction (XRD) for precise determination of complex structures and element-specific localized electronic responses to redox through X-ray absorption spectroscopy (XAS).
Typical changes in oxygen vacancy concentration, Δδ, during redox cycling between reduced and reoxidized states are in the order of Δδ ≥ 0.1 mol O mol −1 cation.The high brilliance of photons and broad range of energies accessible at a synchrotron X-ray source furnish exceptional signal-to-noise ratio for probing structural and electronic changes arising from the small concentration of oxygen defects and enable measurements under in situ or in operando experimental conditions that align with realistic STCH operating parameters.Herein, in operando refers to an experiment carried out under conditions that fully align with all conditions of a representative redox water-splitting cycle (i.e., an experiment to study reoxidation behavior at realistically high temperatures and representative steam:H 2 ratio), while in situ refers to an experiment in which a single condition, or subset of conditions, is being probed (i.e., an experiment at elevated temperature but under ultrahigh vacuum).However, although possible, the full range of operating conditions of STCH materials are for the most part prohibitively difficult to reach for fully in operando experiments and challenging at best with typical in situ experimental apparatuses, with temperatures ranging from 700 °C to 1,600 °C, and oxygen partial pressures ranging from (pO 2 ) = 1 bar to ≤1 × 10 −20 bar.While certain ambitious in situ experiments can capture subranges of the full operating range of STCH materials and are sometimes necessary for unraveling structural behavior at high temperatures, careful sample preparation and experimental design allow ex situ experiments to reveal a wealth of structural and electronic responses of materials to reduction.
Developing next-generation STCH materials will be dependent on understanding the structural and electronic changes that enable reversible oxygen defect formation in complex metal oxides.The high brilliance and variable energy of synchrotron sources provide the resolution and sensitivity necessary for probing the impact of a relatively small percentage of oxygen defects in a bulk material, unobtainable with conventional laboratory-scale techniques.This informs an increasingly useful computational model of these high defect oxides.Despite this, the use of synchrotron characterization has been limited in the STCH field.Herein, we strive to inform researchers of the mechanistic insight that synchrotron experiments can provide into STCH material behavior and provide the best practices for carrying out experiments and interpreting the resultant data.First, we demonstrate a methodology of quenching complex metal oxides in known reduced states such that oxygen vacancies are preserved upon cooling to room temperature.Quenching of samples into a range of known reduced states is achieved using a combination of variable oxygen partial pressures and determination of mass loss during quenching using thermogravimetric analysis (TGA) to determine, and ultimately control, the extent of reduction.The ability to quench oxygen deficient states enables the ex situ probing of a material's structural and electronic response to reduction.We therefore describe the applicability of XAS and XRD techniques for characterizing STCH materials which includes the information that each technique provides, best practices for carrying out the experiments, and advanced data analysis methods for interpreting complex data.Finally, recently developed tools for in situ characterization of water-splitting materials under STCH operating parameters are summarized.

Preparation of controllably reduced samples for ex situ experiments
Although in operando or in situ experiments are ideal, the conditions required for STCH cycling are typically prohibitive and well-designed ex situ experiments can provide valuable insights.In this section, we detail the best practices for quenching in defects for ex situ characterization.
Generating fully oxidized materials for ex situ studies is relatively straightforward.Fully oxidized samples may be prepared by heating and cooling in atmosphere or at 1 bar O 2 .Even under 1 bar of O 2 , most samples develop small levels of oxygen vacancy concentration at or near the surface during cooling, with the depth-dependent concentration depending on bulk diffusivity, surface reaction kinetics, and morphology of samples.However, for fully oxidized samples, the difference in oxygen vacancy concentration between the near surface and bulk should be small when compared with the oxygen vacancy concentration between operational reduced and reoxidized states.In other words, care should be taken to not overanalyze small differences between fully oxidized samples, but when comparing fully oxidized and reduced samples, the results should be relatively insensitive to the exact mechanism of creating the fully oxidized samples.
By contrast, making representative reduced STCH samples for ex situ testing requires careful preparation.Essentially, almost all reduced STCH samples are metastable at room temperature in atmosphere.Only the slow kinetics of surface reoxidation and bulk oxygen diffusion prevent reduced samples from reoxidizing in air.At low temperatures, bulk oxygen diffusion in oxide ceramics is anticipated to be much slower than surface reactions, that is, some surface reoxidation is possible in the top couple of nanometers while the bulk is expected to remain in a reduced state (Li et al., 2021).Therefore, the samples are reduced under controlled reducing conditions and are then "quenched" to room temperature where slow kinetics traps the desired reduced states.
One method for determining the amount of reduction quenched into the ex situ sample is to use thermogravimetric analysis (thermogravimetric analysis) either to quench reduced samples or replicate as precisely as possible the cooling rates, temperatures, and pO 2 values used for sample reduction.In Figure 2, we give an example plot of thermogravimetric analysis data taken during isobaric (fixed pO 2 ) cooling from high temperature for preparation of reduced samples of BaCe 0.25 Mn 0.75 O 3 (BCM) for ex situ studies.In thermogravimetric analysis, the change in oxygen vacancy concentration can be seen as (and calculated from) the loss of sample mass during the experiment, with reduction appearing as mass loss and reoxidation as mass gain.Using the thermogravimetric analysis data, the reduction state trapped at room temperature can be quantified.Similarly, if bulk samples are prepared in a furnace, the pO 2 and cooling rate can be duplicated using thermogravimetric analysis to approximate the ex situ sample reduction.
Achieving the desired degree of reduction in quenched samples requires controlling the temperature and pO 2 .Ideally, pO 2 and temperature are chosen to align with a condition of interest.However, if the equipment used for reducing samples cannot meet the conditions of interest, alternate temperatures and pO 2 can be used to reduce samples to equivalent extents.Environments with the same pO 2 are more reducing at high temperatures than at low temperatures, and similarly, environments with the same temperature are more reducing at lower pO 2 than higher pO 2 .The pO 2 range of interest can span over 20 orders of magnitude, and pO 2 is often very difficult to control precisely over this full range due to the presence of trace leaks and virtual sources of O 2 in many vacuum systems.Options for controlling pO 2 include using blends of gas with known fractions of pO 2 .
Samples must also equilibrate at the reducing conditions, and equilibration times differ dramatically between bulk, nanoparticle, and thin film samples.One ideal way to measure equilibration times is to examine bulk samples using thermogravimetric analysis (TGA) under reducing environments and determine the length of time necessary for the mass loss to reach an asymptote that aligns with the in operando condition that the ex situ experiment is attempting to mimic (for example, if 1% of the oxygen is removed during reduction in a water-splitting cycle of interest, then the sample being prepared for ex situ analysis should also be monitored such that 1% of the oxygen is removed).Monitoring for asymptotic mass loss is important to prevent samples from developing steep reduction gradients that can be unrepresentative when later probed with surface-sensitive measurements.
Controlling the cooling rate is an equally important consideration when preparing ex situ reduced oxides.Under isobaric pO 2 conditions, the samples tend to reoxidize during cooling.The maximum rate of reoxidation that will occur during cooling depends on a sample's reoxidation kinetics (at a given temperature) and the difference between the current and equilibrium oxygen vacancy concentrations and availability of oxygen.In general, the faster the sample is cooled, the more closely the extent of reduction will be preserved from high temperature into the quenched sample.
A number of precautions should be taken when performing measurements on and analyzing the data from ex situ reduced samples.The formation of surface layers of carbonate, hydroxyl, or other non-oxide coatings is possible and potentially more likely in reduced samples due to the reactivity of oxygen vacancies.Special care should be taken when using structural characterization techniques that are selective to only the sample surface.Another important consideration is that the effects of oxygen vacancy concentration on atomic and electronic structures can exhibit temperature dependence.For example, the concentration of oxygen vacancies causes concentrationdependent structural distortions, often referred to as chemical expansion, where crystal lattices tend to expand when oxygen vacancies are present.This chemical expansion is temperature dependent, so the difference in lattice parameters between reduced and reoxidized samples will depend on the temperature (Marrocchelli et al., 2012).Likewise, the location of oxygen vacancies in the crystal lattice may differ between high and low temperatures, with an increased probability of vacancies occurring at higher energy sites at high temperatures due to contributions from configurational entropy and k B T terms.While the examination of samples ex situ can be representative of materials under operating conditions, the aforementioned warnings stress the importance of not assuming an identical behavior between the quenched oxides and samples maintained under operating conditions.The following two sections focus specifically on XAS and XRD techniques for characterizing quenched samples.
3 X-ray absorption spectroscopy for solar thermochemical hydrogen production material characterization X-ray absorption spectroscopy (XAS) is a powerful technique for characterizing STCH materials.XAS provides invaluable insights into the oxygen binding environment and quantification of cation oxidation states, enabling the redoxactive elements to be identified and unraveling more complex phenomena, such as those where the cations reduce further in multi-cation systems with more than one redox-active element.Comparing the oxidation binding environments and cation oxidation states of STCH materials before and after, or during, reduction can provide significant insights into watersplitting mechanisms.Here, we provide a brief introduction to the fundamentals underlying XAS (Section 3.1) and how XAS is applicable for characterizing STCH materials specifically (Section 3.2), and finally, introducing advanced data analysis methods for interpreting the complex spectra (Section 3.3).

Brief introduction to X-ray absorption spectroscopy
X-ray absorption spectroscopy probes the excitation of electrons from core orbitals to unoccupied or partially occupied orbitals, or to the continuum.Incident photons generated by a synchrotron source over a range of welldefined energies interact with a material, and X-ray absorption (I/I o ) at each energy step is measured.Sharp increases in absorption occur at edges where the incident X-ray energy corresponds to the energy necessary for ejecting a core electron to a partially filled or unoccupied orbital.The term XANES, or X-ray absorption near edge structure, refers to the portion of the spectrum containing the edge step, a loosely defined region just before the edge termed the pre-edge, and the region up to approximately 50 eV after the edge (Figure 3).The XANES region can provide quantitative information on the oxidation state, extent of hybridization, and insights into coordination geometry.At higher energies, the core electron is ejected to the continuum with kinetic energy dependent on the energy of the incident photon.To understand the higher energy region, referred to as the extended X-ray absorption fine structure (EXAFS), the wave nature of an electron must be considered.The ejected electron can interact with electrons from neighboring atoms causing it to backscatter.The outgoing wave can interfere with the backscattered wave either constructively or destructively.At incident photon energies resulting in constructive interference, the probability of absorption increases, while destructive interference results in a decreased probability of adsorption (Calvin, 2013).Thus, the EXAFS region often exhibits periodic oscillations (Figure 3).The Fourier transform of this region can be modeled to provide the identity of neighboring atoms and quantitative information such as coordination number and bond distances.One significant advantage of XAS is that the technique is suitable for both crystalline and amorphous materials.For more detailed information on XANES and EXAFS analyses, interested readers are referred to the following sources: Stöhr (1992), Bunker (2010), Calvin (2013), and Frati et al. (2020).
Most XAS data are collected using transmission mode, total electron yield (TEY), or fluorescence yield (FY), as represented in Figure 4.All three of these techniques measure the absorption of X-ray photons as a function of energy.More specifically, they are a measure of the absorption coefficient as a function of energy, as defined in Eq. 1.
where I o is the intensity of the incident photons, I is the intensity of the photons after the beam passes through the sample (unabsorbed photons), μ is the absorption coefficient (linear attenuation coefficient), and t is the thickness of the material.
In the transmission mode, the linear attenuation of X-rays [μ(E) t = ln(I o /I)] is measured directly by monitoring the intensity of photons before (I o ) and after (I) the sample is positioned, often using ionization chambers.Due to the strong attenuation of soft X-rays by many materials, transmission experiments are typically carried out only with hard energy X-rays capable of penetrating the sample.FY detection is an alternative technique for measuring X-ray absorption, whereby a fluorescent photon is emitted during the relaxation of an electron from a higher energy state to the core hole formed by the initial X-ray absorption event.
The intensity of fluorescent photons is proportional to the fraction of the incident photons absorbed μ(E) ∝ I f /I o , where I f is the intensity of fluoresced electrons.Electron yield methods for XAS detection rely on the measurement of electrical current generated from electrons excited from core levels during absorption, where the current is proportional to the fraction of incident X-rays absorbed.These techniques are particularly advantageous at soft X-ray energies where methods based on photon emission struggle due to low escape depths.Additional measurement techniques include partial fluorescent and photoemission-based methods.For additional information on all of the aforementioned techniques and the best practices for collecting high-quality data with each method, readers are directed to the books 'XAFS for Everyone' by Scott Calvin (Calvin, 2013) and 'Introduction to XAFS: A Practical Guide to X-ray Absorption Fine Structure Spectroscopy' by Grant Bunker (Bunker, 2010).

Application of X-ray absorption spectroscopy to solar thermochemical hydrogen production materials
In an XAS experiment, the energy of the absorption event is dependent on both the element and the specific electronic transition probed.The ability to observe an electronic transition is governed by the electric dipole selection rules (e.g., single electron, Δℓ ± 1) 31 , thus transitions from 1s to 2p orbitals and 2p to 3d orbitals result in high absorption intensity.The full range of X-ray energies available at a synchrotron X-ray source enables absorption measurements on most of the elements across the periodic table, which include low Z elements such as oxygen, making XAS a powerful tool for investigating material structures, chemical properties, and electronic behaviors of STCH materials.The oxygen 1s electron binding energy is 543.1 eV, relative to the vacuum level (Thompson et al., 2009).XAS at the oxygen K-edge excites 1s electrons to 2p orbitals and is essentially a measure of empty 2p states, providing a wealth of information on oxygen-metal hybridization.As an example, the oxygen K-edge spectrum of MnO 2 is provided in Figure 5 and is representative of the oxygen in 3d transition metal oxides where hybridization occurs between the oxygen 2p and transition metal 3d orbitals.The lower energy peaks between approximately 525 eV and 535 eV are within the pre-edge region and correspond to O 2p states hybridized with Mn 3d, with the doublet arising from crystal field splitting.The first peak with maximum intensity at 528 eV arises from the spin down t 2g and spin up e g , which are too close in energy to separately resolve, while the second peak, 2.4 eV higher in energy, is related to the spin down e g .The broader peaks at higher energy are O 2p with Mn 4s and 4p character.The area under the curve related to O 2p with transition metal 3d character is quantitatively correlated to the extent of oxygen-transition metal hybridization (Suntivich et al., 2014).Figure 6 overlays the oxygen K-edge spectrum of MnO 2 with the BCM water-splitting material.Significant differences in the pre-edge region are readily observable arising from the different oxygen binding environments, which include covalency and electronic configuration.Transition metal oxides play a prominent role in materials being pursued for STCH applications.XAS provides a means to compare the oxygen binding environments between materials, and more importantly, between a single material in reduced and reoxidized states necessary for providing insights into watersplitting mechanisms.
Reduction of a redox-active cation is integral to the two-step metal-oxide STCH cycle.Lower energy X-rays (termed soft X-rays when <1 keV and tender X-rays when between 1 and

FIGURE 7
X-ray absorption spectra at the Mn 2p edge demonstrating the shift in energy associated with oxidation state which includes a sample that is predominately Mn 2+ from MnO (blue), Mn 3+ from Mn 2 O 3 , and Mn 4+ from MnO 2 (black).The L 3 edge is at approximately 640 eV, while the L 2 edge is at approximately 652 eV (see above for discussion on the L 2 and L 3 edges).Multiple peaks are present for all three samples in the L 3 and L 2 edges, indicating that the samples are not pure and have some contribution from other Mn oxidation states.
5 keV) probe lower energy electron transitions in metals and transition metals, while higher energy X-rays (>5 keV termed hard X-rays) probe higher energy transitions.For example, at the Mn L-edge, 2p electrons are excited to unoccupied or partially occupied 3d states with soft X-rays, whereas at the higher energy K-edge, 1s electrons are excited to 2p states with hard X-rays.The L-edge is split into two states due to spin orbital coupling, resulting in two separate peaks in the L-edge spectrum.The lower energy L 3 edge involves transitions from the 2p 3/2 , while the higher energy L 2 edge involves transitions from the 2p 1/2 with 3/2 and 1/2, referencing the total angular momentum quantum number, J = l + s (Figure 7).The L-and K-edges shift to lower energies with reduction as a consequence of greater nuclear shielding by the surrounding electron cloud.Thus, XAS at both the L-and K-edges provide a means to determine which cations are redox active and the extent of reduction in the mixed cation systems, enabling elucidation of the redox behavior in STCH materials (Figure 7).It should be noted that despite probing transitions to the 3d orbitals (which are hybridized with the oxygen 2p orbital), the L-edge does not necessarily provide information on oxygen-metal hybridization.For example, with 3d transition metals, the local effects such as the 2p spin orbital coupling and 2p-3d electron interaction contribute strongly to the L-edge spectrum shape such that, unlike the oxygen 1s, it is not a measure of unoccupied states (Frati et al., 2020).In addition to oxidation state information from the XANES region of metals and transition metals, detailed analysis of the EXAFS region can provide structural information such as identification of neighboring atoms and bond distances.This method is particularly useful for weakly crystalline or amorphous materials where XRD provides little assistance.
When choosing an XAS measurement technique, it is imperative that the experimentalist be cognizant of limitations and potential pitfalls.TEY probes only the first few nanometers of the sample (de Groot and Kotani, 2008) and is therefore a surface-sensitive technique, the results from which may not be representative of the bulk material.FY probes the bulk of the sample material from the first hundreds of nanometers to several micrometers depending on the incident energy and sample composition, making it an effective method for studying bulk changes in STCH materials.However, FY is not without its limitations.In thick samples where the element of interest is concentrated, as is the case with metal oxide powders, selfabsorption is a prevalent issue.In fluorescence yield, at an absorption edge, the penetration depth of the incident photons decreases such that they interact with fewer atoms and the relationship between the florescence signal and absorption is no longer linear resulting in artifacts being introduced into the spectral intensity and shape (Bunker, 2010).At energies above the edge, when the element being probed is concentrated, it becomes more likely that the incident photons will be absorbed regardless of the energy suppressing the EXAFS oscillations (Calvin, 2013).There are experimental techniques to reduce or eliminate self-adsorption, such as working with thinner samples and changing the angle of the incident energy or detector angle, as well as data analysis methods that correct for self-absorption, details of which can be found elsewhere (Bunker, 2010;Calvin, 2013).Often, it is possible to collect TEY and FY data simultaneously.Although small differences between the two methods may be expected, significant differences in relative peak intensities is a strong indication that self-absorption is an issue, an example of which is provided in Figure 8.It is recommended that both TEY and FY be collected, when possible, for STCH materials.

Data analysis
The oscillations of the EXAFS region are a well-understood phenomenon that can be accurately modeled as a function of wavenumber with the EXAFS equation (Eq.2) 27 .

FIGURE 8
Mn L-edge of 12R-Ba 4 Ce 1 Mn 3 O 12 collected using FY[(A), black] and TEY [(B), red].Due to spin-orbit coupling, the L-edge is comprised of two edges, with the lower energy edge arising from J = (1 + s) = 3/2 transitions termed L 3 and the higher energy from J = (1 − s) = 1/2 transitions termed L 2 .From degeneracy of the states, it is expected that the L3:L2 intensity ratio is approximately 2:1.Self-absorption is readily apparent in the FY spectrum where the relative intensities are not correct.
The EXAFS equation is a sum of all paths between the absorbing atom and the scattering atoms, termed scattering paths, or i in the equation.For examples, when Ce is the absorbing element in CeO 2 , the predominant path contributing to the first peak in the Fourier transform would be the path between Ce and its nearest oxygen neighbors.S o is an amplitude reduction factor that accounts for the experimental Fourier transform amplitude being less than that predicted with the EXAFS equation, N is the degeneracy of the scattering path, D is the half path length or distance between the absorbing and scattering atom, σ 2 is a measure in the variance of the half path length, k is wavenumber, f is related to the probability of scattering, and δ is a phase shift that acts as a correction for assuming the electron maintains a constant wavenumber throughout the round trip (the wavenumber varies due to its interaction with the absorbing and scattering atoms' potential wells).After some data processing (such as normalization to the incident energy, removing the background contribution, and converting from energy to wavenumber), the Fourier transform of the EXAFS region can be fit using the EXAFS equation, where N, D, and σ 2 are the fit parameters.δ(k) and f(k) are atomic number dependent and are typically calculated within the fitting software.S o is chemically transferable and often determined via running a standard with known structure, therefore with known D and N, making S o readily determined from running a fit on a standard where the coordination number is known for each path.Figure 9 provides an example of an X-ray absorption spectrum at the Mn K-edge of an STCH material and its corresponding Fourier transform of the EXAFS region.From fitting the EXAFS, neighboring species can be identified with bond distances and coordination number quantitatively determined (Alia et al., 2017).For example, in Figure 9, the first peak in the Fourier transform corresponds to the Mn-O paths, thus fitting the first peak would provide the average number of oxygen atoms neighboring each manganese and the average bond distance.Depending on the quality and extent of the EXAFS data collected, the fit can incorporate atoms beyond the nearest neighbors to the absorbing element.There are excellent software packages available free to the research community for data processing and fitting that include detailed tutorials and other resources with extensive information on the best practices for data analysis that target both the novice and experienced experimentalists (Ravel et al., 2005;Webb, 2005).
The EXAFS equation is not applicable to the XANES region of the spectra.Analysis of the XANES region is typically more qualitative, although some quantitative information is retrievable.If multiple molecular species are present in a sample, the observed spectrum becomes a linear combination of contributions from the different species.As such, a linear combination fit to the experimental XANES spectrum can often provide the relative weight fractions of the constituents.Particularly applicable to studying STCH materials, oxidation states can also be quantitatively determined.At higher oxidation states, there is less electron shielding of the positively charged atomic nucleus and a greater energy is required to eject a core electron, moving the edge position to higher energy.When using the edge shift to determine the oxidation state, it is important to ensure that the shift is not an artifact of a change in the experimental configuration, such as a shift in the monochromator angle.It is usually possible to collect XAS data on a known reference material simultaneously with data collected on the sample.The spectrum from the reference material provides a reference spectrum for each sample measured.Each sample spectrum is then essentially tied to its corresponding reference spectrum so that by aligning the reference spectra and shifting the sample spectra equivalently, robust alignment is possible and shifts in edge energy can be associated with differences in the oxidation state.
As described above, pre-edge features in the XANES region are often attributable to more complex material properties and behavior, such as transition metal-oxygen hybridization and electronic configuration, both important for relating the STCH material structure to performance.However, the extraction of physical meaning from the pre-edge region requires correlating the peaks in a XANES spectrum to the specific molecular orbitals into which an excited core electron transitions, often necessitating a computational approach.The density functional theory (DFT) is a valuable tool for unraveling complex contributions to X-ray absorption spectra.Like the Schrodinger wave equation that other quantum chemical methods rely on, the DFT is rooted in an electronic structure (Hohenberg and Kohn, 1964).However, instead of attempting to solve the wave equation, DFT uses exchange-correlation energy functional to define a system's electronic density.These functionals are a mathematical description of how electrons interact with each other.This key element of DFT is how one can approach an interacting problem (solving the wave equation) by mapping it exactly to an easier-to-solve noninteracting problem (electron density) (Burke, 2007).From the electron density, it is then possible to determine the ground state electronic structure of the system from which materials' physical properties arise.Calculated vibrational force constants provide insights into chemical bonding or behavior with other species (i.e., electrical polarizability and relative energies in chemical reactions).In relation to spectroscopic properties one can also find excitation energies to the nth excited state and the probabilities of their occurrence, also known as the density of states (DOS) (Kohn et al., 1996).As XAS is essentially a measure of partially occupied or unoccupied states available for an electron transition, DFT is a useful tool for associating X-ray absorption spectra peaks to specific core level electron transitions to unoccupied molecular orbitals.In other words, calculating the DOS identifies the unoccupied orbitals to which the core electrons can transition.However, the substantial difference in resolution between theory and experiment can readily lead to erroneous assignments; for example, from the DFT, it is known that the two lowest unoccupied states in MnO 2 are the t 2g spin-down and e g spin-up, but these do not correspond to the first two peaks in the X-ray absorption spectrum (Figure 5) as the energy states are only separated by a few tenths of an electron volt and cannot be resolved using XAS.The first peak represents a combination of these two electronic transitions.Although computationally more intensive, simulating the X-ray absorption spectrum, as detailed below, can be critical to resolving multiple contributions to complex spectra for robust peak identifications, such as those resulting from hybridized states.In addition, when a simulated spectrum matches the experimental spectrum, it can be assumed that the DFT structural model that was used to simulate the XANES spectrum is accurate and the subsequent computational analyses are more likely to be reliable.The following two subsections provide further details on the DFT, and calculations are used to simulate XANES spectra.

Density functional theory
The density functional theory (DFT), which has consistently been the leading workhorse in electronic structure calculations with applications across a wide variety of systems and phenomena (Becke, 2014;Tozer and Peach, 2014;Jones, 2015;Besley, 2021), can be extremely useful in prediction, verification, and analysis of the experimental results in STCH research.The DFT is rooted in the two Hohenberg-Kohn theorems (Hohenberg and Kohn, 1964) which state that for any system of interacting electrons: 1) the external potential is uniquely determined by the ground state electron density, and 2) this density and the associated energy can be found variationally by minimizing the total energy as a function of the density.In practice, the DFT is used mostly within the Kohn-Sham (KS) (Kohn and Sham, 1965) framework, which seeks a fictitious noninteracting system of electrons whose densities equal that of the interacting system.Such a system is under the influence of the multiplicative KS potential: v KS = v ext + v H + v XC , where v ext is the external potential on the interacting system, v H is the classical Hartree potential corresponding to the electron density, and v XC is the exchange-correlation (XC) potential, which accounts for the quantum mechanical exchange interaction and the effects of correlation that arise from the interaction between individual electron pairs beyond a mean-field treatment.The exact form of v XC is as yet unknown and depending on the problem at hand, various approximations (Kohn and Sham, 1965;Becke, 1988;Lee et al., 1988;Perdew et al., 1996a;Rappoport et al., 2011) are used for this term.Additionally, for the sake of computational efficiency, many DFT calculations (especially those that use plane waves as the basis function) replace the collective effect of the core electrons with a pseudopotential (Hellmann, 1935;Schwerdtfeger, 2011) term dependent on the atomic species.
Even though KS-DFT does not guarantee an equivalence (Kohn et al., 1996) [the highest occupied level is an exception (Janak, 1978;Perdew et al., 1982;Perdew and Levy, 1997)] between the single-particle wave functions/energies of the KS system and those associated with the addition or removal of electrons in the real system, such an equivalence is often observed in many systems to an appreciable extent, although this can be subject to the use of exchange-correlation functionals pertaining to the generalized KS framework (Seidl et al., 1996;Kümmel and Kronik, 2008;Perdew et al., 2017).Therefore, the KS energies and wave functions are routinely used as approximate quasiparticle (Onida et al., 2002) counterparts in several contexts (e.g., band structure plots).However, it must be noted that commonly used (semi-)local XC potentials, within the local density approximation (LDA) (Kohn and Sham, 1965) or the generalized gradient approximation (GGA) (Perdew et al., 1996a;Perdew et al., 1996b), tend to delocalize the electron density unphysically (Mori-Sánchez et al., 2008) due to inherent self-interaction errors (Perdew, 1985), whereby each electron interacts spuriously with itself via the mean-field of all electrons.For STCH materials, the description of the electrons in the localized d and f orbitals of the transition metal/rare earth atoms can be heavily affected by such an error.Therefore, v KS is usually supplemented with an additional orbital-dependent, nonlocal Hubbard correction (Himmetoglu et al., 2014) term that with a user-specified Hubbard parameter U [often obtained empirically (Wang et al., 2006)] enforces electron localization by penalizing fractional occupation in some predefined atomic d and/or f orbitals.Inclusion of a fraction of nonlocal exact exchange in v KS (analogous to the Hartree-Fock calculations) can also counteract the self-interaction error, although typically at a much higher computational expense.
Regarding STCH research, the DFT can play a crucial role in complementing XRD experiments for structural analysis and XAS experiments for the analysis of electronic structure.The former is typically accomplished with the help of a calculation of geometry relaxation, which entails computing the net force on each ion as the gradient of the total energy and optimizing the structure until the force drops below a certain threshold.Additionally, the thermal motion of the ions can be simulated with the help of molecular dynamics (MD) (Car and Parrinello, 1985;Iftimie et al., 2005) or calculation of phonon (Giannozzi et al., 1991;Parlinski et al., 2005) modes.DFT-based computational tools are also routinely used in prediction and analysis of various spectroscopic experiments.KS-DFT is found to be particularly adept at simulating the K-edge absorption spectra (such as the oxygen K-edge) since multiplet effects typically play a minor role in such excitations and consequently, the corresponding many-body state can be approximated as a single-reference system (de Groot and Kotani, 2008), consistent with the existing common approximations to the exchange-correlation functional within the DFT.

Simulation of X-ray absorption spectra
Within the framework of KS-DFT, one tries to simulate the X-ray absorption spectrum by using Fermi's golden rule #2 to compute the absorption probability: where ω is the angular frequency of the absorbed photon, T is the transition operator, |Ψ i 〉 (|Ψ f 〉) is the initial (final) state associated with the X-ray absorption process, and E i (E f ) is the energy thereof.In the large-wavelength limit, T can be expressed as the dipole operator e • R, where e is a unit vector along the polarization direction and R is the many-body position operator.
In practice, the absorption cross-section probability is found with the help of two separate KS self-consistent field (SCF) calculations: 1) a ground state SCF calculation for simulating the initial state and 2) another DFT calculation run on a positively charged system in which the core of the excited atom is represented by a modified pseudopotential mimicking a full core-hole (FCH) [alternatively, this calculation can be run with a neutral cell producing the excited-electron core-hole (XCH) state].The latter calculation, referred to as the FCH state calculation hereafter, is used for extracting information on the final core-excited states.In the so-called single-particle treatment, the dipole matrix element is approximated as where r is the one-body position operator, ϕ core denotes a core orbital, and φc , which is the orbital of the excited electron after the absorption, is an unoccupied (conduction) KS orbital of the FCH state.The recently developed many-body X-ray absorption spectroscopy (MBXAS) method (Liang et al., 2017;Liang and Prendergast, 2018;Liang and Prendergast, 2019) seeks to improve upon the single-particle treatment by approximating the many-body state |Ψ f 〉 (|Ψ i 〉) by a Slater determinant (SD) constructed by populating the relevant KS orbitals obtained from the FCH (ground) state calculation.Note that each possible final state will have a different set of occupied orbitals, and then, with some algebra, the transition matrix element reduces to where the sum is over all the unoccupied orbitals ϕ c of the ground state system and |Ψ c i 〉, which is a neutral core-excited state with an excited electron in the orbital ϕ c , is represented (non-selfconsistently) by an SD composed of the ground state-occupied valence KS orbitals in combination with ϕ c .The term 〈Ψ f |Ψ c i 〉, which is an inner product between two Slater determinants built from orbitals of two different SCF calculations (i.e., FCH and ground state calculation), can be reexpressed as the complex conjugate of a determinant composed of the overlap matrix elements between the FCH and ground state orbitals.Thus, in a nutshell, MBXAS expresses each transition matrix element as a weighted sum of single-particle transitions 〈ϕ c |e.r|ϕ core 〉 with the weighing factor given by the projection of the corresponding SD expressed in terms of the ground state orbitals onto an SD representing the actual final state of interest.It is important to note here that the abovementioned simulation methods, which use the KS eigenvalues/eigenfunctions in their pristine form, present a computationally cheaper and faster alternative to the techniques (Vinson et al., 2011;Gilmore et al., 2015;Vorwerk et al., 2019) rooted in many-body perturbation theory (Shirley, 1998;Rohlfing and Louie, 2000).They are particularly useful and efficient for systems with large unit cells (such as materials with defects) or for systems where the effect of finite-temperature lattice dynamics needs to be incorporated via molecular dynamics (Prendergast and Galli, 2006;Pascal et al., 2014;Roychoudhury et al., 2021a).
Using pseudopotentials (i.e., instead of including the core electrons explicitly) necessitates the use of an overall empirical shift to align to experimental energy scales.In addition, we must also account for the so-called chemical shifts associated with different chemical or coordination environments of the same excited element, for example, at symmetry inequivalent atomic sites in a given crystal or in entirely different materials.The use of two different sets of pseudopotentials for the ground and excited state calculations prevents us from relying on the raw total energy differences to align spectral contributions of inequivalent atoms.The relative excitation energies for distinct atomic sites are determined with respect to a common theoretical reference, the isolated atom (England et al., 2011;Jiang et al., 2013;Roychoudhury et al., 2021b).To this end, for excitation of atom X, the term (E f − E i ) in Eq. 3 is replaced by an effective formation energy difference, where e X (e ′ X ) denotes the total energy of the isolated, neutral atom X using the pseudopotential employed in the ground state (FCH) calculation.The total energy E f of any final state can now be estimated as where ε denotes a Kohn-Sham orbital eigenenergy of the core hole-excited self-consistent field.The index i runs over all orbitals that are occupied in the final state of interest (f) but unoccupied in the so-called excited-electron core-hole (XCH) state (Prendergast and Galli, 2006), which is the core-excited state with the lowest energy.Conversely, the index j runs over all orbitals that are unoccupied in the final state f but occupied in the XCH state.The total energy E XCH of the XCH state can be found from a separate DFT SCF calculation on the neutral system employing the pseudopotentials used in the FCH calculation.In the simplest case involving the creation of only a single electron-hole pair upon X-ray absorption [denoted as f (1) in Liang and Prendergast (2018)], the realignment of the final state energies is accomplished with a simpler expression (since there is only one term in each sum above) for a system with N valence electrons and a final excited state with a hole in the core orbital and an electron in orbital f > N: This alignment scheme is crucial not only for comparing the onset energy of the spectra of different materials but also for simulating the accurate line shape of the resultant spectrum for a material in which the excited atomic species occupy multiple inequivalent sites.In particular, the aforementioned formalism has proven to be highly effective in simulating the O K-edge spectra, both in terms of line shape and onset energy, of transition metal oxides (Roychoudhury et al., 2021b).
Figure 10 shows the simulated O K-edge X-ray absorption spectra of CeO 2 (ceria) using the single-particle treatment (Eq.4) and the MBXAS method (Eq.5), with the latter being calculated separately using three different values of Hubbard U. It must be noted that even though Ce 4+ is in a 4f 0 configuration in ceria, it is still important to use a Hubbard correction for the f electrons due to the hybridized nature of the electronic orbitals.Hubbard parameters ranging from U = 4.5 eV to U = 6 eV have been suggested (Fabris et al., 2005;Da Silva et al., 2007;Loschen et al., 2007;Ismail et al., 2011;Grieshammer et al., 2014;Grieshammer, 2018) in the existing

FIGURE 10
Simulated XAS plots for oxygen K-edge of CeO 2 .In the "FCH" plot, the dipole matrix element is obtained with the single-particle treatment presented in Eq. 4. This plot is obtained from KS-DFT calculations using a Hubbard parameter of U = 5.4 eV.The three MBXAS plots, which correspond to different values of Hubbard U (namely, 5.4, 4.5, and 0 eV), use the many-body expression of Eq. 5 to calculate the dipole term.
literature.Comparison of the top three spectra in Figure 10 can attest to the importance of Hubbard U in XAS simulation of STCH materials.The two plots at the bottom of Figure 10, both of which are obtained using the same Hubbard parameter (U = 5.4 eV), show a severe underestimation of the pre-peak height in the single-particle FCH treatment when compared to the MBXAS spectrum, which is in good agreement with the experimental results (Aguiar et al., 2010).This underestimation is reported for oxides of a large number of transition metals and can be explained with the help of a simple tight-binding model (Liang and Prendergast, 2018).Thus, MBXAS, which can be seen to rectify the relative peak heights, is a particularly useful tool for research on STCH materials, which are typically transition metal/rare-earth oxides.Finally, we note that DFT-based simulations provide valuable information regarding the nature and constitution of the single-particle orbitals of the excited electrons.As representative examples, in Figure 11, we show the isovalue plots of three KS orbitals, transitions to which contribute appreciably to the first three peaks in the MBXAS spectrum shown with the green curve in Figure 10.

X-ray diffraction for solar thermochemical hydrogen production material characterization
XRD with a synchrotron X-ray source provides sensitivity and resolution unachievable with conventional laboratory-source diffractometers.Synchrotron XRD enables complex crystal structures to be solved (and resolved), subtle changes in the structure (such as strain or changes in bond distances) to be detected, and low-concentration impurity phases to be identified and quantified, features critical to unraveling STCH material behavior.The following provides a brief description of the fundamentals underlying XRD (Section 4.1) and its applicability for characterizing STCH materials (Section 4.2).

Brief introduction to X-ray diffraction
X-ray diffraction (XRD) is a scattering-based technique suitable for investigating the identity and arrangement of atoms within a given lattice.Here, an elementary background on the theory behind X-ray diffraction will be presented such that the reader is made familiar with how a diffraction experiment can aid the structural characterization and development of STCH materials.
Photon scattering is encountered when the wavelength of radiation is comparable to interatomic spacing.For incident photons having wave vector, k, scattering from two positions separated by r into scattered photons with wave vector, k', the difference in the path length between the scattered wave vectors must be equal to (k-k') • r = Q r, i.e., the normal projection of r onto Q.For elastic scattering, the phase difference between the incident and scattered waves must be equal to 2π/λ times the path difference in order for the photon to constructively interfere with the wave scattered from the adjacent atom.The Laue equations are formed when this one-dimensional example is further expanded into three dimensions.The intensity of scattering from a lattice of ordered atoms is defined by the crystal structure factor: The value of F(Q) is nonzero when Q coincides with a reciprocal lattice vector.f j (Q) in the first summation of the expression above defines the atomic form factor, which is a measure of the X-ray "scattering efficiency" for a given atom.Since photons scatter via interactions with electrons, the form factor naturally increases monotonically as a function of atomic number, Z.The first summation in Eq. 9 describes scattering from atomic sites within a given unit cell, where j is the number of distinct atoms within the unit cell.The second summation extends the scattering contribution to the lattice sum and together these summations make up the crystal structure factor (Als-Nielsen and McMorrow, 2011).

Application of X-ray diffraction to solar thermochemical hydrogen production materials
As it pertains to STCH processes, ex situ XRD measurements provide the ability to solve the crystal structures of new candidate STCH compounds, demonstrate phase purities of as-synthesized and redox-cycled materials, and determine/differentiate lattice changes associated with thermal expansion and oxygen nonstoichiometry in the reduced state (Metcalfe et al., 2019).These properties are critical for evaluating the performance of a new candidate STCH material by identifying/separating contributions from the known impurity compounds/phases.In situ XRD measurements (e.g., high temperature and controlled gaseous environments) provide a valuable route for deriving relationships between atmospheric redox conditions and material structure (i.e., structure-property relationships) (Metcalfe et al., 2019;Mastronardo et al., 2020).The advantage of XRD with a synchrotron source is exemplified by Bell et al. (2022), where impurity phases at less than 0.4 weight percent were identified within the high-purity STCH material BCM, as well as by Strange et al. (2022), where a novel BCM polytype was identified and found to form under high-temperature reducing conditions.The results described in both studies would have not been observed with conventional laboratory-source XRD.The high brilliance offered by a synchrotron X-ray source also enables kinetic studies where structural changes can be monitored as fast as 1 kHz with hybrid photon counting detectors.
Using the most basic definition of a structural refinement, the atomic identities and positions are determined from experimental diffraction patterns in order to describe the crystal structure.Since the atomic form factor exhibits a dependence on atomic number, X-rays are relatively insensitive to oxygen atoms, especially when present in a lattice of elements with high Z (e.g., period 6 and 7 elements), which is the case for the top contending STCH materials.This limitation can be overcome with neutron diffraction, where the coherent scattering cross sections are independent of Z, and scattering from oxygen is appreciably relative to the cations.Additionally, in XRD, atoms with similar atomic numbers are not readily distinguished in a structural refinement since differences in form factor contributions are relatively low.To overcome this limitation, resonant X-ray diffraction is more suitable, whereby the diffraction intensities are recorded as the energy scanned over an atom's absorption edge.The relative changes in diffraction intensities then allow for differentiation of site occupancies.
The two most common XRD measurement geometries are Debye-Scherrer (transmission) and Bragg-Brentano (reflection), both exhibiting advantages and disadvantages.Transmission measurements require small sample volumes and are often performed in capillaries that offer a uniform cylindrical geometry with respect to the scattering angle.Sample absorption is of particular concern in transmission geometry, but a correction to the observed diffraction intensities as a function of the scattering angle can be made with a priori knowledge of μ t (where μ is the linear attenuation coefficient and t is the sample thickness) and an estimated particle packing fraction.The primary advantage of transmission geometry is the ability to perform measurements in extreme environments (e.g., high/low temperatures and pressures and controlled atmospheres) with low attenuation from the sample cell (typically thin-walled borosilicate glass, quartz, or sapphire).Reflection geometries overcome the absorption problem but present additional complications.When performing in situ measurements with a heating stage, the poor thermal conductivity of STCH materials (e.g., thermally insulating ceramics) results in a severe temperature gradient between the heating element and radiated sample volume, particularly when the sample is in the form of a pressed puck.Powders measured in reflection geometry must be uniformly deposited using a suitable nonreactive solvent in order to prevent movement of the sample during gas flow.Using this approach, there are also potential complications with sample roughness, i.e., the change in apparent density as a function of powder depth.Alternatively, a thin film may be used, but there are uncertainties regarding the structure and STCH performance of thin films when compared to their bulk counterparts.In reflection geometries, since the heating element and sample are likely to expand as a function of increasing temperature, an external calibration should be performed such that the degree of expansion is known (and reproducible) and a manual correction of offset sample position can be applied to the data.Knowledge of the temperaturedependent sample displacement is especially important when differentiating the relative contributions of thermal vs chemical expansion of an oxide during STCH.

Survey of solar thermochemical hydrogen production applicable in situ tools
While laboratory in situ measurements of STCH materials are widely performed (such as thermogravimetric analysis mass spectrometry of the effluent gas to correlate mass loss or gain to the species lost), in situ measurements via X-ray methods at a synchrotron are much less common.For redox-active materials used in water-splitting processes, the temperatures required for reduction typically exceed 1,000 °C and require fine control over pO 2 .Some of the primary challenges associated with in situ measurements are access of photons at necessary energies to the sample (typically achieved through the choice of the window/ capillary material), and a well-controlled sample environment exhibiting uniform temperature under simultaneous gas flow with a known composition.Heating samples to over 1,000 °C is particularly challenging due to instrumental and material-based limitations.Conduction heating is highly inefficient for X-ray cells, which cannot be completely enclosed by heating elements and insulation.On the materials side, the number of mechanically stable and inert compounds, which are still permeable to X-rays, is very limited (Rothensteiner et al., 2015).Furthermore, the interpretation of structural data can be especially difficult at high temperatures where lattice motions distort the time-averaged depiction of the crystal structure.However, there are notable successes in the literature where synchrotron experiments have been carried out under realistic operating conditions, or a subset thereof, which are summarized below.
Among the studies that have used in situ synchrotron X-ray measurements for thermochemical cycling conditions, the most common techniques used are XPS, XAS, and XRD.Thermochemical energy storage materials have been investigated at temperatures under 1,000 °C for materials systems such as Mg-H-F (Tortoza et al., 2018) for using H 2 as a hydrogen storage material, SrFeCuO 3 (Vieten et al., 2019) for oxygen storage and air separation, and reactive carbon composites such as BaCO 3 (Møller et al., 2020) and CaMg(CO 3 ) 2 (Humphries et al., 2019).These studies often use a quartz capillary tube with thin walls to hold the sample, allowing a controlled gas flow across the sample during transmission XRD, X-ray photoemission spectroscopy (XPS), and/or XAS measurements.
For ceria-based materials, which in some cases contain other elements such as Hf or Zr, the cerium K-edge has been measured under realistic thermochemical conditions, with the reduction occurring at 1,500 °C and reoxidation at 800 °C in the presence of steam or CO 2 (Rothensteiner et al., 2015;Rothensteiner et al., 2016;Rothensteiner et al., 2017).These studies use synchrotron X-rays to ultimately perform XAS, XRD, and simultaneous mass spectrometry (MS) in transmission geometry to measure full thermochemical cycles, ramping up to 1,500 °C in Ar for reduction of the ceria, and cooling to 800 °C before introducing CO 2 or steam for reoxidation.The heating source used was an infrared focusing furnace to heat an alumina inner tube while maintaining a cooler environment for the outer quartz tube leading to the gas outlet and MS (Rothensteiner et al., 2017).Additional studies have also been done using XANES under conditions up to 1,100 °C with exposure to ambient air, vacuum, and hydrogen to examine surface concentration of Ce 3+ in ceria-zirconia catalyst materials (Yuan et al., 2020).

Conclusion
Developing a robust redox-active metal oxide that performs well under reasonable operating conditions is key for realizing STCH as a commercially viable process for clean hydrogen generation.Progress in material development is contingent on developing structure-property relationships, which requires a mechanistic understanding of reactions, material stability, and failure routes.Synchrotron radiation is a powerful tool for characterizing STCH materials.X-ray absorption spectroscopy identifies those cations that are redox active and the extent to which they are reduced under quenched conditions.By probing the oxygen K-edge, XAS provides quantitative information on the extent of metal-oxygen hybridization and qualitative information on how the oxygen electronic environment responds to defects and reoxidation.X-ray diffraction provides the ability to solve the crystal structure of new materials and quantify purity, as well as identify secondary phases and determine how the crystal structure responds to oxygen defects during a redox cycle.The high-energy X-rays generated by a synchrotron source opens up opportunities for in situ experiments unobtainable with laboratory-scale techniques.While it would be impossible to cover all the information and details necessary to design and carry out the various synchrotron experiments, the goal herein has been to make interested readers aware of what can be learned from the various techniques relevant to STCH materials and provide a wealth of resources that can be further consulted.

FIGURE 2
FIGURE 2 Results of representative BCM thermogravimetric analysis under reducing conditions.(A) Change in relative mass (black), instantaneous mass/ initial mass, with temperature (red), demonstrating mass loss during 1,350 °C isothermal reduction.(B) Relative mass from (A) plotted as a function of temperature (left axis) and with relative mass transformed into oxygen off-stoichiometry, δ (inverted right axis), for BCM.Results demonstrate quenching of a known reduced state for subsequent ex situ investigations.
FIGURE 3X-ray absorption spectrum of Mn K-edge with XANES and EXAFS regions identified.

FIGURE 4
FIGURE 4 Simplified schematic of an XAS experiment, with the sample in green.The orange boxes represent detectors, where I o measures the incident photons, I T measures the transmitted photons, and I F measures the fluoresced photons.The arrow with e − indicates electrons at the sample surface from excitation measured via electron yield experiments, indicated by the gray box labeled TEY.
FIGURE 5 X-ray absorption spectrum of MnO 2 at the O K-edge.
FIGURE 6 X-ray absorption spectra of reoxidized BCM (solid line) and MnO 2 (dashed line) at the O K-edge.

FIGURE 9
FIGURE 9Mn K-edge X-ray absorption spectrum (A) of a calcium cerium-manganese oxide STCH material and its corresponding Fourier transform (B).It is the Fourier transform that would be fit with the EXAFS equation to determine nearest neighbors, bond distances, and coordination numbers.

FIGURE 11
FIGURE 11Isovalue plots of KS orbitals of the core-excited state with significant contribution in the first (A), second (B), and third (C) peaks in the O K-edge MBXAS plot, as shown in Figure10(green line).The red and gray spheres represent oxygen and cerium atoms, respectively.Note that these are all hybrid orbitals containing O p and Ce d/f characters, with different degrees of hybridization.