research papers
GIWAXS using microbeam applied in halide perovskite films for high spectral resolution and mapping capability
aShanghai Institute of Applied Physics, Chinese Academy of Sciences, Shanghai 201800, People's Republic of China, bShanghai Synchrotron Radiation Facility, Shanghai Advanced Research Institute, Chinese Academy of Sciences, Shanghai 201204, People's Republic of China, cUniversity of Chinese Academy of Sciences, Beijing 100049, People's Republic of China, dThe Institute for Advanced Studies, Wuhan University, Wuhan 430072, People's Republic of China, and eSchool of Environmental Science and Engineering, Frontiers Science Center for Transformative Molecules, Shanghai Jiao Tong University, Shanghai 200240, People's Republic of China
*Correspondence e-mail: [email protected], [email protected], [email protected]
Synchrotron-based grazing-incidence wide-angle X-ray scattering (GIWAXS) has become an essential tool for investigating halide perovskite film microstructures. Nevertheless, GIWAXS employing 100 µm X-ray beams inevitably suffers from footprint effects manifesting as diffraction peak broadening and spectral distortions. Herein, we demonstrate that microbeam GIWAXS (µ-GIWAXS) effectively mitigates these footprint effects, achieving enhanced spectral resolution and accuracy. By systematically tuning the grazing incidence angle, we quantitatively illustrated how the footprint length contributes to diffraction peak broadening in µ-GIWAXS. In addition, we optimize µ-GIWAXS data acquisition through beam-direction-perpendicular fly scanning, which improves data quality while preserving resolution and minimizing radiation damage at individual spots. Through extended exposure experiments, the effect of µ-GIWAXS microbeam irradiation on the structural stability of perovskite films was systematically evaluated. Furthermore, µ-GIWAXS enables spatial mapping of perovskite film microstructural characteristics including crystallinity, preferential crystal orientation, PbI2 concentration and lattice strain – critical information for addressing inhomogeneity challenges in large-area perovskite devices. Given the growing demand in scalable perovskite optoelectronic applications, µ-GIWAXS warrants increased attention due to its dual advantages of superior spectral resolution and spatially resolved analysis capabilities.
Keywords: perovskite; µ-GIWAXS; radiation damage; mapping; footprint effect.
1. Introduction
Metal halide perovskites have emerged as a research focus in optoelectronic fields (Zhang et al., 2023
; Song et al., 2025
; Sakhatskyi et al., 2023
) due to their strong light absorption, high carrier mobility, low cost, and low-temperature solution processability (Wei et al., 2024
; Duan et al., 2024
). Single-junction perovskite solar cells (PSCs) have already reached a certified power conversion efficiency (PCE) of 27.3%, rivaling the performance of commercial crystalline silicon photovoltaics (see https://www.nrel.gov/pv/interactive-cell-efficiency.html). In addition, perovskite–silicon tandem devices also have reached a PCE of 35.0%, highlighting their great potential for commercialization (Aydin et al., 2024
). The microstructural characteristics of perovskite films are universally recognized as critical determinants of device performance, a subject of extensive investigation in recent years (Zhou et al., 2022
; Guo et al., 2024
). Among various characterization techniques, synchrotron-based grazing-incidence wide-angle X-ray scattering (GIWAXS)/grazing-incidence X-ray diffraction (GIXRD) has emerged as one of the most widely applied techniques to investigate the microstructures of perovskite films (Liang et al., 2023
). Recent years have witnessed a rapid increase in GIWAXS-based studies of PSCs (Park et al., 2023
; Li, Su et al., 2023
; Huang et al., 2023
).
In GIWAXS measurements, the diffraction signals are collected from the entire X-ray illuminated region on the sample surface, which is significantly expanded by the so-called footprint effect at grazing incident angles. This geometric effect advantageously increases the X-ray–sample interaction volume at the surface, contributing to GIWAXS's widespread application in thin film characterization. While surface sensitivity can be enhanced by reducing the grazing angle, this comes at the cost of further elongating the X-ray footprint. However, at synchrotron facilities dedicated to X-ray diffraction/scattering techniques, GIWAXS is typically conducted with beam sizes ranging from hundreds of micrometres. Consequently, the footprint can exceed common sample dimensions at extreme grazing angles, leading to peak broadening since X-ray diffraction peaks originating from different positions along the elongated footprint project as overlapping signals on the area detector (Smilgies, 2009
). Furthermore, both non-uniform beam intensity profiles and sample microstructural inhomogeneities may introduce distortions in the diffraction patterns.
An effective strategy to mitigate this issue involves replacing conventional millimetre-scale beams with microbeams or even nanobeams in GIWAXS experiments. The elongated footprint thus becomes much shorter, thereby narrowing the diffraction peak distribution projected onto detectors. In fact, microbeam X-ray diffraction (µ-XRD) or even nanobeam XRD (n-XRD) has been proven successful in analyzing small samples or a small area/volume of samples (Goldsmith & Walker, 1983
; Nozaki et al., 2020
). By combining high-brilliance synchrotron sources with advanced X-ray focusing optics, sub-micrometre polychromatic/monochromatic beams can be achieved for spatial mapping of grain orientation, strain/stress, and micro-topography (Li et al., 2022
). For example, scanning n-XRD in transmission geometry was used to simultaneously probe the morphology and the structural properties of spin-coated MAPbI3 perovskite films (Lilliu et al., 2016
). Nevertheless, µ-XRD/n-XRD does not take intrinsic advantages of grazing incidence ideal for detecting thin films. For such purposes, a few GIXRD experiments, using a beam with a diameter of only a few micrometres (µ-GIXRD), have been reported to investigate structurally inhomogeneous samples (Sáez-Martínez et al., 2021
). Recently, the application of µ-GIWAXS has been successfully introduced to halide perovskite films, demonstrating its powerful capability for spatially resolved microstructural mapping (Chan et al., 2025
; Zhang et al., 2025
). However, a systematic investigation of the fundamental aspects of µ-GIWAXS—including the quantitative relationship between beam footprint and spectral broadening, optimized data acquisition strategies, radiation damage assessment under high-flux-density microbeam and comprehensive multi-parameter spatial mapping—remains lacking.
In this study, µ-GIWAXS/µ-GIXRD measurements of halide perovskite films were explored by using both an X-ray beam with vertical size of about 3 µm at the BL03HB beamline and that of about 1 µm at the BL17UM beamline of the Shanghai Synchrotron Radiation Facility (SSRF). Parallel GIWAXS measurements with a beam size of about 300 µm were simultaneously conducted at BL14B1 for direct comparison. Unless otherwise specified, beam sizes are defined as the full width at half-maximum (FWHM). The experimental results demonstrate significantly reduced footprint effects in µ-GIWAXS spectra, yielding enhanced spectral resolution. Additionally, scanning µ-GIWAXS was implemented to spatially resolve microstructural heterogeneities across perovskite films. Given that these local microstructural variations are known to influence the long-term stability and efficiencies of large-area perovskite devices, µ-GIWAXS represents an effective and valuable tool for accelerating their commercial development (He et al., 2023
; Jiang et al., 2025
; Yin et al., 2025
).
2. Material and methods
2.1. Materials
All the raw materials were purchased from commercial suppliers without further purification. Patterned indium tin oxide (ITO) glass was purchased from Advanced Election Technology Co. Ltd. N,N-dimethylformamide (DMF), dimethyl sulfoxide (DMSO), chlorobenzene (CB) and isopropanol (IPA) were purchased from Sigma Aldrich. Lead iodide (PbI2) was purchased from TCI. Formamidinium iodide (FAI), methylammonium chloride (MACl), methylammonium iodide (MAI) and caesium iodide (CsI) were all purchased from Xi'an Polymer Light Technology. Except for ITO glass, all the materials were stored in a nitrogen-filled glove box to avoid water.
2.2. Perovskite precursor and films preparation
ITO glass was ultrasonically cleaned by sequentially washing with detergent, deionized water, ethyl alcohol and isopropanol (IPA). Before use, the cleaned ITO substrates (with lateral dimensions of 15 mm × 20 mm) were dried using a nitrogen gun and then treated with ultraviolet ozone for 20 min to improve the surface wettability of the surfaces. For the pure MAPbI3 film preparation, PbI2 (461 mg) and MAI (159 mg) in DMF:DMSO = 700 µL:78 µL was stirred for 2 h and then filtered with a 0.22 µm polytetrafluoroethylene (PTFE) filter. 40 µL of MAPbI3 precursor was spin-coated at 4000 rpm for 30 s (4 s acceleration to 4000 rpm). 100 µL CB was dropped on the film at 22 s before the end of the spin-coating. The film was annealed by 100 °C hot plate in nitrogen glovebox for 10 min. For the preparation of the mixed hybrid Cs0.05FA0.85MA0.1PbI3 (hereafter denoted as CsFAMAPbI3 for brevity) film, a precursor solution was first prepared by dissolving PbI2 (691.5 mg) and CsI (19.5 mg) in a solvent mixture of DMF and DMSO in a 900 µL:100 µL ratio. The solution was stirred for 1 h and subsequently filtered using a 0.22 µm PTFE filter. A 40 µL of the PbI2 solution was spin-coated onto ITO glass at 1500 rpm for 30 s, followed by annealing at 70 °C for 1 min in a nitrogen atmosphere. After cooling to room temperature, 80 µL of a pre-filtered organic salt solution—comprising FAI:MAI:MACl in a weight ratio of 90 mg:6.39 mg:9 mg dissolved in 1 ml of isopropanol—was spin-coated onto the PbI2 layer at 2000 rpm for 30 s. The resulting film was then transferred from the nitrogen glovebox to a hot plate in ambient air (30–40% relative humidity) and annealed at 150 °C for 15 min.
2.3. GIWAXS characterization
GIWAXS measurements were performed at the BL03HB, BL14B1 and BL17UM beamlines of the SSRF, China (Xing-Yu et al., 2015
; Wang, Wang et al., 2024
; Kong et al., 2025
). The key experimental parameters are summarized in Table 1
and are described below.
| ||||||||||||||||||||||||||||||||||||||||||||||
BL14B1 is a bending magnet beamline (1.27 T) equipped with a collimating mirror, a focused Si(111) double-crystal monochromator (DCM), and a vertically focusing mirror. No micro-focusing optics are employed; the beam is shaped by a slit system to approximately 300 µm (V) × 500 µm (H) (FWHM) at the sample position, with an energy resolution of ΔE/E ≃ 1.91 × 10−4 and a beam divergence of ∼0.25 mrad (V). The intrinsic beamline divergence and the divergence at the sample position are essentially identical, as no focusing optics are used. A Huber 5021 six-circle diffractometer provides multi-axis sample positioning. The 2D GIWAXS patterns were recorded using a MarCCD225 detector (pixel size ≃ 73 µm × 73 µm) at a sample-to-detector distance of ∼250 mm.
BL03HB is a Laue microdiffraction beamline based on a superbend dipole magnet (2.29 T), featuring a two-stage focusing design. The optical system comprises a toroidal mirror for horizontal pre-focusing, a four-bounce Si(111) monochromator that can switch between monochromatic and white-beam modes, and a vertically focusing mirror that creates a secondary source. Two sets of Pt-coated Kirkpatrick–Baez (KB) mirrors are installed downstream: the first set (at 46 m from the source) focuses the beam to ∼4 µm × 4 µm in the protein crystallography sector, and the second set (at 51 m) further demagnifies the beam. The µ-GIWAXS measurements in this work were conducted in the materials science sector using the monochromatic mode at 10 keV (ΔE/E ≃ 0.96 × 10−4), with a beam size of approximately 3 µm (V) × 5 µm (H) at the sample position. Due to the large numerical aperture required for KB micro-focusing, the effective beam divergence at the sample position (∼2.8 mrad) is substantially larger than the intrinsic beamline divergence upstream of the KB system. The 2D µ-GIWAXS patterns were recorded using a Pilatus 2M area detector (pixel size 172 µm × 172 µm) at a sample-to-detector distance of ∼250 mm. A six-degrees-of-freedom sample stage (three translational and three rotational motions) enables precise positioning and alignment.
BL17UM is a high-performance undulator beamline (CPMU20, 20 mm period, 160 periods) designed for micro-crystallography. The beamline features interchangeable optics: a cryo-cooled Si(111) DCM (ΔE/E ≃ 2 × 10−4) for high-resolution applications and a double multilayer monochromator (DMM) for high-flux experiments. The beam is focused by a pair of KB mirrors to approximately 1 µm × 0.7 µm (FWHM) at the sample position, with an effective divergence of ∼2.5 mrad (H) × ∼1.0 mrad (V). In this work, the DCM mode at 10 keV was used for µ-GIWAXS measurements. The 2D patterns were recorded using an Eiger X 16M detector (pixel size 75 µm × 75 µm) at a sample-to-detector distance of ∼300 mm. An MD3-UP diffractometer is used to control precise positioning and alignment.
At all of the beamlines, the diffraction spectra were calibrated using NIST SRM660b (LaB6) sample. The GIWAXS/µ-GIWAXS data were processed by using MATLAB and Fit2d software.
3. Results and discussion
3.1. Comparison of µ-GIWAXS and GIWAXS
Fig. 1
(a) presents the µ-GIWAXS geometry featuring a set of KB mirrors for X-ray micro-focusing. To quantify footprint effects, Fig. 1
(b) models a grazing-incidence (α) X-ray beam (of diameter d) generating an out-of-plane diffraction peak at 2θ onto the area detector with a vertical projection of D. From geometric analysis, the footprint in Fig. 1
(b) is deduced to be Δ = d/sinα with the projection D = . The beam diameters at the BL14B1 and the BL03HB beamlines are approximately 300 µm and 3 µm, respectively. At an incidence angle of 1° for a diffraction angle of 11.33° [corresponding to the perovskite (001) reflection at a photon energy of 10 keV], the Δ values at the BL14B1 and the BL03HB beamlines are 17.19 mm and 171.9 µm with a D of 3143.66 µm and 31.44 µm, respectively. Considering the significant beam divergence angle resulting from KB mirror focusing at the BL03HB beamline, we have incorporated the contribution of beam divergence into our analysis (Fig. S1 of the supporting information). Under identical conditions, we further incorporated the specific beam divergence angles of the two beamlines: β1 = 0.25 mrad for BL14B1 and β2 = 2.8 mrad for BL03HB. A typical sample-to-detector distance of 250 mm was used. Substituting these values into D =
−
, the revised D at the BL14B1 and the BL03HB beamlines are calculated as 3573.68 µm and 1486.76 µm, respectively. These results clearly demonstrate that the use of a large beam size leads to considerable diffraction peak broadening, due to the increased beam projection. It should be noted that this estimation of instrumental broadening arising from footprint effects is simplified and does not account for the intrinsic diffraction peak width. The observed peak width is also determined by the natural diffraction profile, which is governed by the microstructural properties of the films, including grain size, strain, defects and (Holder & Schaak, 2019
).
| Figure 1 (a) Schematic of X-ray microfocusing and the µ-GIWAXS experimental setup. (b) Schematic illustrating the broadening of a diffraction peak in GIWAXS due to the footprint effects. 2D GIWAXS patterns of a CsFAMAPbI3 perovskite film measured with an incidence angle of 1° at the beamline BL14B1 in (c) and BL03HB in (d), respectively. The corresponding 1D GIWAXS diffraction spectra derived from (c) and (d) are shown in (e) and (f), respectively. |
Beyond the footprint contribution, the instrumental broadening of the peak is a convolution of several additional factors. First, the beam divergence at the sample position—arising from the combined contributions of the intrinsic source emittance and the focusing optics—broadens diffraction peaks along the radial direction. This specific contribution has been explicitly incorporated into our projection analysis above. Second, the detector point-spread function (PSF) and the finite pixel sizes of the detectors used in this work impose a fundamental lower limit on the achievable angular resolution at a given sample-to-detector distance (Smilgies, 2009
). Third, the geometric projection of scattering from the curved Ewald sphere onto a flat planar detector introduces a non-linear distortion of diffraction rings that exacerbates at higher scattering angles, contributing to additional apparent broadening (Jiang, 2015
). To mitigate this, the detectors were positioned approximately normal to the direct beam, and rigorous geometric corrections—accounting for detector tilt and sample-to-detector distance uncertainty—were applied during data reduction using LaB6 calibration. Any error in this calibrated sample-to-detector distance would propagate directly into the q-vector scale, introducing systematic shifts and artificial broadening in the integrated profiles. Furthermore, the finite energy bandwidth (ΔE/E ≃ 10−4) of the Si(111) monochromator produces a wavelength spread that broadens peaks proportionally to tanθ (Birkholz, 2006
). Finally, the integration geometry during data reduction (e.g. azimuthal integration) introduces geometric smearing when curved diffraction rings are binned over a finite azimuthal angle onto a Cartesian grid (Jiang, 2015
). Among these contributions, the beam footprint and beam divergence are the dominant instrumental factors in the present GIWAXS/µ-GIWAXS measurements, while the remaining contributions are comparatively minor but are accounted for in the data reduction procedure.
To systematically investigate the influence of footprint effects on GIWAXS measurements in perovskite films, we conducted comparative experiments on CsFAMAPbI3 samples using two different beamline setups: BL14B1 and BL03HB. The GIWAXS diffraction patterns collected at the same photon energy of 10 keV with an incidence angle of 1° are presented in Figs. 1
(c) and 1
(d), respectively. This incidence angle was chosen to ensure that the standard beam footprint remains within the sample boundaries; at smaller grazing angles the footprint would exceed the sample dimensions, convolving footprint broadening with sample-edge truncation effects (Pflüger et al., 2017
). The intense scattering feature near the beamstop arises from specular reflection, Yoneda diffuse scattering at the critical angle (αc), and parasitic scattering, and is confined to the low-q region (q < ∼5 nm−1), well separated from the perovskite and PbI2 reflections of interest (Yoneda, 1963
; Hexemer & Müller-Buschbaum, 2015
). Both datasets display the characteristic diffraction features of the α-phase perovskite with a preferential crystal orientation, along with a minor PbI2 phase. Notably, every diffraction feature recorded with the microbeam is visibly sharper [Fig. 1
(d)] than its counterpart obtained with the larger beam [Fig. 1
(c)], signaling a substantial mitigation of footprint broadening.
To quantitatively assess the impact of footprint reduction on diffraction quality, we analyze the diffraction peaks corresponding to the CsFAMAPbI3 (001) reflection at q ≃ 10.0 nm−1 and the PbI2 impurity phase at q ≃ 9.0 nm−1. As shown in Fig. 1
(f), both peaks appear significantly sharper, more symmetric, and closely approach an ideal Gaussian profile in contrast to their corresponding counterparts in Fig. 1
(e). Critically, the FWHM of the CsFAMAPbI3 (001) peak decreases notably from 0.462 nm−1 [Fig. 1
(e)] to 0.098 nm−1 [Fig. 1
(f)], underscoring the dramatic enhancement in spectral resolution achieved via footprint minimization. Importantly, the broader, asymmetric peak shapes observed at BL14B1 [Fig. 1
(e)] arise from a convolution of two factors: (i) non-Gaussian flux density distribution across the elongated footprint, and (ii) spatial inhomogeneity in the perovskite film's microstructure, as further illustrated in Fig. S2. This coupling of instrumental and sample-related inhomogeneity effects leads to peak distortions, including apparent peak splitting in the PbI2 signal as well as an irregular asymmetry in the CsFAMAPbI3 (001) peak. Such distortions not only complicate the interpretation of GIWAXS data but also risk introducing unwanted errors in and structural analysis (Qin et al., 2021
). For example, the splitting of the PbI2 peak in Fig. 1
(e) could easily be mistaken for a secondary phase. Moreover, footprint broadening often merges adjacent diffraction peaks, making it difficult to distinguish between structurally similar phases, such as tetragonal and cubic perovskites (Yalcinkaya et al., 2022
; Steele et al., 2023
). In particular, studies aiming to resolve depth-dependent microstructural variations via GIWAXS measurements at different incident angles must be approached with caution: changes due to footprint variations with varying incidence can masquerade as spurious depth-sensitive structural transitions. While this phenomenon is also present in µ-GIWAXS, the absolute footprint lengths are approximately two orders of magnitude smaller, making this effect considerably less critical but not entirely negligible. In contrast, the µ-GIWAXS data exhibit no discernible peak splitting artifacts in either the perovskite or PbI2 reflections, and deliver more symmetric peak profiles (asymmetry factor closer to 1.0) that closely follow a Gaussian line shape (adjusted R2 > 0.99; see Table S1). These well defined peak profiles, with spurious modulations largely eliminated, significantly simplify data analysis.
To further elucidate how geometric footprint effects constrain reciprocal-space resolution, we performed µ-GIWAXS measurements on an MAPbI3 film by varying the incident angle at the BL03HB beamline. The FWHM of the MAPbI3 (110) diffraction peak from the 1D integrated GIWAXS was observed to change accordingly [Figs. 2
(a)–2(c)]. With an X-ray beam longitudinal size of approximately ∼2.9 µm, the footprint on the sample ranges from 553.9 µm at an incident angle of 0.3° to 18.5 µm at 9°, with the corresponding FWHM of the (110) diffraction peak decreasing from 0.120 nm−1 to 0.079 nm−1, as shown in Figs. 2
(a) and 2(b). Fig. 2
(c) presents a stitched image of diffraction patterns for perovskite (110) at various grazing incident angles, intuitively showing the reduction of the diffraction ring width with increasing grazing angle. By correlating the grazing incident angle, X-ray footprint length and the FWHM values [shown in Fig. 2
(d)], we can further quantitatively illustrate the footprint effects on diffraction peak broadening. Specifically, as the grazing incident angle increases from 0° to 3°, the FWHM of the diffraction peaks decreases significantly, and, beyond 3°, the reduction in FWHM slows down. There is an almost linear relationship between the footprint length and the FWHM of the diffraction peaks, which is consistent with the established geometrical broadening models (Fig. S3): the projection D on the detector scales proportionally with the footprint Δ when footprint-induced broadening is the dominant instrumental contribution. By extrapolating the curve in Fig. 2
(d) to zero footprint length, the intrinsic peak width of the present µ-GIWAXS setup at BL03HB can be estimated, which represents the combined contribution of microstructural broadening and the residual non-footprint instrumental terms (beam divergence, detector PSF and energy bandwidth). This footprint-corrected width can serve as the input for quantitative line-profile analysis such as Scherrer grain-size estimation and Williamson–Hall strain-size separation (Smilgies, 2009
; Williamson & Hall, 1953
). We emphasize that peak positions and derived lattice parameters remain largely robust against footprint broadening, preserving the quantitative reliability of standard GIWAXS for phase identification and d-spacing determination.
| Figure 2 (a) Waterfall plot showing 1D µ-GIWAXS spectra by the azimuthal integration for 7.5 nm−1 ≤ q ≤ 18 nm−1 at varying incident angles measured from an MAPbI3 film. (b) 1D integrated µ-GIWAXS profiles of the MAPbI3 (110) peak at varying incident angles. (c) Stitched 2D diffraction pattern from the same perovskite film measured at different grazing incident angles. (d) Correlation between the MAPbI3 (110) diffraction peak FWHM and the incident angles as well as their corresponding footprint lengths. 2D µ-GIWAXS patterns obtained using (e) fixed-point and (f) fly scanning modes, respectively. (g) Radially integrated intensity plots along the MAPbI3 (110) diffraction ring, measured at a fixed point and by beam-direction-perpendicular fly scanning. (h) Statistical distribution of the curve differences calculated from (g). |
The unique advantages of µ-GIWAXS are most pronounced at very small incidence angles (e.g. 0.05°), as the footprint of the standard beam becomes overwhelmingly severe. However, the simple geometric model (Δ = d/sinα) becomes insufficient in this context, where refraction, Yoneda enhancement and the sharp transition in penetration depth complicate the footprint-broadening interpretation (Yoneda, 1963
; Steele et al., 2023
). The incidence angles employed in this study (0.3–9°) are all above αc (∼0.22°), ensuring the validity of the geometric model.
Equivalent footprint lengths can be achieved using different combinations of beam size and incidence angle, but with fundamentally different depth sensitivities. For instance, a footprint of ∼5.7 mm corresponds to the large beam (∼300 µm) at 3° or the microbeam (∼3 µm) at ∼0.03°—the latter below αc and thus confining the probe to the near-surface evanescent region. µ-GIWAXS therefore uniquely enables simultaneous minimization of footprint broadening and enhancement of surface sensitivity.
One thing that must be considered is the µ-GIWAXS data acquisition strategy. It is well known that the statistical reliability of XRD data from polycrystalline materials deteriorates with reduced beam size, as the interaction volume between X-rays and the material contains an insufficient number of crystallites to ensure statistically representative sampling. Therefore, special caution must be exercised when collecting µ-GIWAXS data. To address this, we propose fly scanning, which enlarges the X-ray–material interaction volume while preserving spectral resolution. This is achieved by continuously translating the sample stage back and forth perpendicular to the X-ray beam. Specifically, during exposure, the stage moves uniformly along the y-direction with a lateral scanning range of 1 mm (as illustrated in Fig. S4). The GIWAXS patterns obtained through fixed point measurement [Fig. 2
(e)] and fly scanning [Fig. 2
(f)] demonstrate almost identical peak widths for MAPbI3 (110) rings (FWHM of 0.109 nm−1 and 0.113 nm−1, respectively), confirming preserved spectral resolution (Fig. S5). Notably, diffraction rings in Fig. 2
(e) exhibit less continuous rings than those in Fig. 2
(f), evidencing improved data statistics. The improvement by fly scanning becomes particularly pronounced in their corresponding radially integrated intensity plots along the MAPbI3 (110) ring in Figs. 2
(e) and 2
(f). Through differential analysis of these plots (Fig. S6), we derived azimuthal intensity differential distribution curves [Fig. 2
(h)], which clearly reveal significantly reduced orientation dispersion achieved via fly scanning. Therefore, this strategy preserves spectral resolution and enhances data statistics by increasing the number of sampled crystallites. Furthermore, it minimizes X-ray irradiation time at individual positions, thereby mitigating radiation damage per unit area—a critical consideration in µ-GIWAXS measurements. The fly scanning inherently averages the diffraction signal over a larger lateral area, which may smooth out local structural variations. Therefore, fly scanning is better suited for obtaining statistically representative information (e.g. accurate orientation distribution), whereas fixed-point or fine-step scanning is preferable when spatial resolution is the primary concern.
3.2. Irradiation effect of microbeam
Prolonged X-ray irradiation is known to cause radiation damage in perovskites, leading to observable microstructural changes (Kirmani & Sellers, 2025
; Stuckelberger et al., 2020
). To achieve an X-ray microbeam, the beamline acceptance angle must be reduced, which inevitably decreases X-ray flux. Although the total flux is reduced, the smaller beam size results in significantly higher flux density. For instance, the flux density at BL03HB is estimated to be about three orders of magnitude higher than at BL14B1, potentially exacerbating radiation damage. While the protective role of nitrogen is recognized for conventional GIWAXS, it remains a question as to whether this protection holds under the extreme conditions of µ-GIWAXS (Francisco et al., 2024
; Svanström et al., 2021
). Therefore, evaluating the microbeam irradiation effects on perovskite films during µ-GIWAXS measurements becomes essential. Additionally, the presence of moisture and oxygen during experiments must be considered, as these factors may accelerate perovskite film decomposition (Zhu et al., 2023
; Zhang & Han, 2020
). To systematically investigate these effects, µ-GIWAXS measurements with an X-ray incidence angle of 0.5° were conducted continuously on the same spot of pure MAPbI3 and CsFAMAPbI3 films, both in air and nitrogen environments. Each µ-GIWAXS pattern was acquired with 30 s irradiation, with total measurement durations of 60 min in each case. Throughout the experiments, the beamline hutch was maintained at ∼20°C and ∼30% relative humidity.
Contour plots of the µ-GIWAXS 1D intensity profile as functions of X-ray irradiation time for MAPbI3 films measured in nitrogen [Fig. 3
(a)] and in air [Fig. 3
(b)] reveal their structural evolution during irradiation. In a nitrogen environment, the MAPbI3 (110) diffraction peak maintains mostly constant intensity with no detectable PbI2 formation. In contrast, in ambient air, progressive attenuation of the MAPbI3 (110) peak concurrent with PbI2 signal accumulation is observed. Notably, the time-dependent (110) peak position drift occurs exclusively in air. Quantitative analysis of these phenomena is present in Figs. 3
(c)–3(f). In Fig. 3
(c), the MAPbI3 (110) peak in air shifts from 9.980 to 9.965 nm−1 during 60 min irradiation, with no equivalent shift detected in nitrogen. Moreover, Fig. 3
(d) displays µ-GIWAXS 1D integrated plots for the MAPbI3 film in air at initial and final stages, with the corresponding µ-GIWAXS patterns shown in Fig. S7. In addition to the MAPbI3 (110) peak intensity drop and position shift, the MAPbI3 (211) peak located at q = 16.5 nm−1 severely decreases. The observed peak shift could align with the reported tetragonal-to-cubic phase transitions induced by thermal effects, proposing localized X-ray heating as a potential driver for this structural transformation. However, intermittent measurements of the irradiated MAPbI3 film at 20 min intervals (Fig. S8) reveal non-reversible (110) peak position changes, conclusively demonstrating that the shift originates from irradiation-induced decomposition in air rather than thermally activated reversible phase transitions. Moreover, the quantitative variation of diffraction peak intensities under irradiation reveals distinct degradation pathways for MAPbI3 in different environments. As shown in Fig. 3
(e), the (110) peak of MAPbI3 in air decreases nearly linearly with irradiation time, accompanied by a continuous rise in the PbI2 peak intensity, indicating progressive decomposition of MAPbI3 into PbI2. In contrast, the same perovskite phase in N2 exhibits negligible degradation with no detectable PbI2 formation. To verify the origin of the intensity variation in air, we further conducted a µ-GIWAXS measurement on another MAPbI3 film every 20 min in air. Fig. 3
(e) shows that the stable MAPbI3 (110) peak intensity throughout the measurement indicates that ambient conditions alone do not trigger rapid MAPbI3 degradation. Thus, the rapid decomposition observed in air arises from the synergistic effect of X-ray irradiation and environmental factors (e.g. moisture/oxygen), as neither factor alone induces significant degradation.
| Figure 3 (a, b) Contour plots of the µ-GIWAXS 1D intensity profile as functions of X-ray irradiation time for MAPbI3 films in (a) nitrogen and (b) air. (c) The MAPbI3 (110) diffraction peak position as functions of time derived from (a) and (b). (d) The µ-GIWAXS 1D intensity profile of the film measured in air at the beginning and at 60 min. (e) The evolution of MAPbI3 (110) and PbI2 peak intensity for the MAPbI3 film under X-ray irradiation in air and nitrogen environments. (f) The corresponding CsFAMAPbI3 (001) and PbI2 peak intensity variations for the CsFAMAPbI3 film under identical irradiation conditions. To enhance the visual clarity, the PbI2 content (q = 8.87–9.12 nm−1) has been magnified 30 times in (a), (b), (e) and (f). (g) The CsFAMAPbI3 (001) peak intensity as a function of irradiation time at beamlines BL03HB and BL17UM, with increased from the lowest (B1) to the highest (B4). |
For CsFAMAPbI3, known for its enhanced environmental stability, the contour plots of the µ-GIWAXS 1D intensity profile as functions of X-ray irradiation time in both N2 and air (Fig. S9) show no pronounced decline in the (001) peak intensity. However, Fig. 3
(f) reveals a gradual decomposition of the α-phase (001) peak in both environments, albeit with distinct PbI2 evolution trends: the PbI2 content remains constant in N2 but slightly decreases in air despite ongoing perovskite decomposition. This phenomenon can be explained by PbI2 from the decomposition of CsFAMAPbI3 decomposed further by X-rays into iodine and metallic Pb (Pb0) (Dang et al., 2017
). It should be noted that excess PbI2 and metallic Pb are known to exacerbate hysteresis in perovskite photovoltaic devices (Lin et al., 2021
; Zhao et al., 2022
). Therefore, one must be aware of the decomposition of that PbI2 in the investigation of irradiation effects of perovskites. In the MAPbI3 films with minimal amount of PbI2, X-ray irradiation in air predominantly drives MAPbI3-to-PbI2 decomposition, with PbI2 content increasing linearly over time [Fig. 3
(e)]. Prolonged irradiation may eventually activate PbI2 photodecomposition, offsetting its accumulation from perovskite degradation. To demonstrate the importance of PbI2 decomposition under X-ray irradiation, Fig. S10 tracks the (001) peak of an FAPbI3 film with 10% excess PbI2 under X-ray irradiation. Despite accelerated FAPbI3 decomposition in air, PbI2 content concurrently decreases (the degradation of PbI2 itself can outpace its generation) whereas, in N2, slower perovskite degradation coincides with stable PbI2 levels. These observations demonstrate that X-ray induced PbI2 decomposition is not negligible in both environments.
The irradiation-induced degradation of perovskite films exhibits a strong dependence on density, a critical factor requiring systematic analysis. Experiments were primarily conducted at the BL17UM beamline, which provides higher and a smaller beam spot compared with BL03HB. By employing an aluminium attenuator, different densities can be achieved at the sample spot. Fig. 3
(g) summarizes the evolution of the perovskite diffraction peak intensity from CsFAMAPbI3 films both in air and N2 under four distinct densities at BL03HB and BL17UM (B1 to B4, see Table S2). Notably, at the highest flux density (B4: 6.6 × 1020 photons s−1 cm−2), the peak intensity decreases by ∼50% within 10 s even in N2, indicating severe X-ray-induced decomposition. A clear trend emerges: higher flux densities accelerate degradation across both environments. To quantify irradiation-driven decomposition kinetics, all curves in Fig. 3
(g) were fitted with an exponential decay model, yielding the characteristic time constants (τ) listed in Table S2. From Fig. 3
(g) and Table S2, it can be concluded that a photon flux density of the order of ≤1015 photons s−1 cm−2 has a negligible impact on CsFAMAPbI3 films. However, with a density of over ∼1017 photons s−1 cm−2, the CsFAMAPbI3 films suffer noticeable decomposition after just 1 min of continuous irradiation. Since radiation damage is directly proportional to the absorbed dose (measured in Gray) and considering the significant variations in both X-ray energy and inherent to synchrotron experiments, we quantified this critical parameter through dose calculations. The absorbed dose in the thin films was computed using the formula =
, with the detailed calculation process provided in Note S1 in the supporting information (Berejnov et al., 2021
). The absorbed dose rate in Mode B1 was calculated to be 1.25 × 106 Gy s−1, yielding a total dose of approximately 108 Gy after 1 h of irradiation. This flux-dependent behavior underscores the necessity of optimizing irradiation parameters for perovskite characterization to minimize beam-induced artifacts. In radiation-sensitive applications, such as perovskite photovoltaics for space-based solar panels or X-ray detectors, the radiation tolerance of perovskite semiconductors is obviously a critical operational criterion, which requires further rigorous and comprehensive investigation. As for the µ-GIWAXS measurements, measuring the sample at a fixed point for a prolonged time should be avoided (Wei & Huang, 2019
; Tu et al., 2021
; Kirmani et al., 2022
). Based on the present study, an inert nitrogen atmosphere can effectively mitigate irradiation-induced degradation in perovskite samples when the absorbed dose rate remains at or below 106 Gy s−1 at an X-ray energy of 10 keV. For comparison, da Silva et al. reported observable damage above ∼104 Gy under 8.5 keV nano-beam condition (Francisco et al., 2024
). The comparatively higher tolerance in the present work is likely attributable to the harder X-ray energy (10 keV, below the Pb L-edges), which results in a lower photoelectric absorption cross-section. These dose thresholds are not universal material constants; they depend on beam energy, dose rate, incidence angle, film composition and thickness, and atmospheric environment. The values reported here are therefore specific to the present experimental conditions and serve as practical reference points for µ-GIWAXS experiments on halide perovskite thin films under comparable configurations.
Beyond photon flux density, perovskite film decomposition is also strongly influenced by the X-ray incidence angle during measurements. The relationship between decomposition and incidence angle is highly non-linear, governed by two competing factors: beam footprint (which determines the illuminated area) and X-ray penetration depth (where intensity attenuates to 1/e of the initial value) (Steele et al., 2023
; Francisco et al., 2024
). At smaller incidence angles (αi < αc), total external reflection confines the beam to a shallow surface layer (several nm). Despite intensity enhancement at the interface, the absorbed dose remains low due to minimal interaction volume. As αi increases toward and beyond αc, the beam footprint decreases proportionally to 1/sin(αi), concentrating the photon flux over a smaller area and increasing the photon density per unit area. Simultaneously, the X-ray penetration depth increases sharply, enabling deeper interaction and greater energy deposition throughout the film thickness. The combined effect of reduced footprint and increased penetration depth leads to maximum absorbed dose in this regime. When αi further increases, X-rays fully penetrate the film and reach the substrate. Although the reduced beam footprint leads to higher local photon density, the smaller irradiated volume results in a lower total absorbed X-ray dose within the perovskite film. Taking all these factors into account, the simplified process of the above different situations is illustrated in Table S3.
3.3. µ-GIWAXS mapping
Solution-processed perovskite films are prone to phase segregation and compositional inhomogeneity, which is particularly critical for large-area device fabrication (Gao et al., 2024
; Li et al., 2018
). These structural imperfections critically compromise material stability, severely degrading the performance and operational lifetime of perovskite-based devices. While conventional GIWAXS is extensively employed for microstructure characterization at macroscopic scales, its utility in resolving localized crystalline heterogeneity is inherently limited due to its large beam spot sizes, which is further elongated by the footprint effect along the incident direction (Epp, 2016
). By contrast, µ-GIWAXS is able to achieve spatially resolved structural analysis through planar raster scanning, achieving micrometre-scale spatial resolution (tens of micrometres) by sequentially collecting diffraction patterns from distinct sample regions.
To validate the spatial resolution of µ-GIWAXS mapping, a 13 mm × 15 mm planar scan was performed on a CsFAMAPbI3 film patterned with `SSRF' letters [Fig. 4
(a)]. The experimental setup, including schematic and photographic illustrations of the scanning configuration, is reported in Fig. 4
(b) and Fig. S11, respectively. Sample positioning was controlled using two stepper motors along two in-plane translation axes as illustrated in Fig. 4
(b): the X-axis, oriented perpendicular to the incident beam (lateral direction), and the Y-axis, oriented approximately parallel to the incident beam (longitudinal direction). The scanning was performed row by row in a left-to-right pattern along the X-direction. After completing a full row of 80 points, the sample stage was incremented stepwise in the Y-direction (from bottom to top) before moving to the start of the subsequent row, tracing a zigzag trajectory. It should be noted that the motors for both axes move strictly within the sample plane to prevent changes in the sample-to-detector distance and the resulting peak shifts in the diffraction patterns. Throughout the scanning process, the incident angle was consistently maintained at 3°. This angle balances spatial resolution and measurement reliability: at smaller grazing angles, the elongated footprint would not only smooth out the local microstructural variations targeted in this study but also cause overlap between adjacent scanning positions, resulting in repeated X-ray exposure and potential irradiation-induced degradation. The entire scan consisted of 80 rows with 50 points per row, and each point was exposed for 5 s. From the acquired µ-GIWAXS dataset, the intensity distribution of the CsFAMAPbI3 (001) diffraction peak was reconstructed into the spatial mapping shown in Fig. 4
(c). The `SSRF' pattern is clearly resolved in the mapping image, corresponding to regions where mechanical scratching during patterning severely degraded perovskite crystallinity. The perfect correspondence between the optical micrograph [Fig. 4
(a)] and the µ-GIWAXS mapping [Fig. 4
(c)] conclusively demonstrates the technique's capability for in-plane microstructure analysis with micrometre-scale spatial resolution.
| Figure 4 (a) A CsFAMAPbI3 film sample, inscribed with the letters `SSRF', as a sample for micro-structural analysis. (b) Schematic of the setup for the µ-mapping experiment. (c) The obtained µ-GIWAXS mapping image of the CsFAMAPbI3 (001) diffraction peak intensity from the CsFAMAPbI3 film. |
To further demonstrate the versatility of µ-GIWAXS mapping, we analyzed a 4 mm × 5 mm area using a 20 × 20 point grid, corresponding to step sizes of 200 µm and 250 µm, respectively. The scanning direction, X-ray incident angle, exposure time per point, and other experimental parameters were consistent with those described in the previous section. From the µ-GIWAXS datasets, spatially resolved mappings were reconstructed to visualize key microstructural parameters including CsFAMAPbI3 (001) peak intensity, PbI2 peak intensity and CsFAMAPbI3 (001) orientation as shown in Figs. 5
(a)–5(c).
| Figure 5 (a, b, c) Maps over a 4 mm × 5 mm area of a CsFAMAPbI3 film showing (a) (001) diffraction peak intensity, indicating the crystallinity of the α-phase perovskite; (b) (001) orientational order parameter, representing the degree of deviation from the out-of-plane orientation; and (c) PbI2 peak intensity. The statistical distributions derived from (a, b, c) is shown in (d). (e) Lattice strain maps over scanning areas of 4 mm × 5 mm of the CsFAMAPbI3 film. |
As shown in Fig. 5
(a), the CsFAMAPbI3 (001) peak intensity displays a highly heterogeneous spatial distribution, revealing spatially varying crystallization degree across the perovskite film surface. The peak intensity can differ by two times at different areas with interconnected high-intensity domains, which can be attributed to stochastic nucleation dynamics inherent to solution-processing. Fig. 5
(b) maps the spatial distribution of orientation order parameters within the perovskite thin films. Given that crystallographic alignment dictates critical optoelectronic performance metrics, such as carrier mobility and grain boundary density, the spatial mapping of crystal orientation provides valuable information for understanding local variations in device-relevant film properties (Li, Shen et al., 2023
; Li et al., 2024
). The orientation order parameter specifically quantifies angular deviation from the (001) preferential orientation (the calculation follows established crystallographic formalism) (Song et al., 2020
). The enhanced diffraction intensity on the film's right side correlates with a superior (001) orientation, which promotes the vertical transport of charge carriers. Conversely, the green regions in the map indicate a stronger (111) orientation, which is known to mitigate ion migration and enhance environmental stability. Comparison of Figs. 5
(a) and 5
(b) reveals that these two parameters—crystallinity and orientation—exhibit spatial inhomogeneity with distinct distributions, indicating that high crystallinity is not exclusively tied to a specific orientation. Notably, a left-to-right intensity gradient is observed in Fig. 5
(c), differing from the stochastic heterogeneity of other parameters. This pattern may be related to the scanning trajectory, where the elongated beam footprint could accumulate irradiation damage from the initial exposure area, causing irradiation-induced decomposition of the perovskite into PbI2. However, a detailed analysis suggests a fabrication-based origin is more plausible. Firstly, beam-footprint-driven damage would primarily produce a Y-direction gradient, contrary to the observed X-direction variation. Secondly, the decomposition extent is considered limited due to the short, non-continuous irradiation duration (≤100 s) and the material's self-healing properties (Wang, Song et al., 2024
). The poor correlation between the PbI2 and perovskite crystallinity distributions from Figs. 5
(a) and 5(c) also supports this conclusion. The observed pattern is more likely to originate from an inhomogeneous distribution of the cationic solution during the second fabrication step of perovskite film, leading to spatially incomplete conversion of PbI2. The statistical data for the regions scanned in Figs. 5(a)–5(c) are summarized in Fig. 5
(d), with detailed parameters (e.g. mean and standard deviation) provided in Table S4. This µ-GIWAXS mapping analysis serves as a powerful quantitative tool for microstructure characterization, demonstrating statistically comparable distributions.
Solution-processed mixed-cation perovskite films typically exhibit an inhomogeneous lattice strain distribution (Xiong et al., 2024
). By analyzing the µ-GIWAXS diffraction dataset acquired from the same scan as Figs. 5
(a)–5(c) (see Note S2 for calculation details), we mapped the corresponding strain distribution, as shown in Fig. 5
(e). Each point in the map corresponds to Δq, defined as Δq = q∥ − q⊥, where q∥ and q⊥ represent the in-plane and out-of-plane scattering vectors, respectively. Due to the inherent `missing wedge' in GIWAXS, the exact in-plane and out-of-plane directions are inaccessible (Qin et al., 2021
). Therefore, q∥ and q⊥ were extracted from the nearest accessible azimuthal angles, yielding an approximate measure of lattice strain. A negative Δq (Δq < 0) signifies tensile strain, mainly resulting from the perovskite's higher coefficient relative to the substrate. During cooling, the greater contraction of the perovskite imposes in-plane tension. Conversely, a positive Δq (Δq > 0) indicates compressive strain. Although less common, compressive strain—often linked to improved long-term stability—can be introduced via specific processing routes. The strain map reveals a non-uniform spatial distribution of tensile and compressive strains within the film, challenging the conventional assumption that only tensile strain is present. This heterogeneity may be attributed to multiple sources, including A-site cation size mismatch, local lattice distortions at residual PbI2 inclusions, compositional gradients inherent to the two-step deposition, and mismatch between the perovskite film and substrate (Xue et al., 2020
; Wang, He et al., 2024
; Chen et al., 2019
). The strain values here reflect the combined contribution of these factors; disentangling individual contributions would require more complementary spatially resolved compositional analysis (e.g. nano-XRF and TOF-SIMS). The co-existence of tensile and compressive strains can lead to local lattice distortions, potentially compromising the structural integrity (Zhu et al., 2019
). Statistical distributions of Δq corresponding to Fig. 5
(e) are provided in Fig. S12. It shows that most of the lattice points exhibit Δq < 0, and the integrated value of Δq is also negative, consistent with the notion that perovskite films prepared on glass substrates are dominated by tensile strain.
Finally, the influence of the scanning protocol was also evaluated. Fig. S13 presents perovskite crystallinity maps obtained from a 1 mm × 1 mm scanning area (with a 50 µm step size) and a 200 µm×200 µm area (with a 10 µm step size). The overlap between scanned points resulted in artificial continuous features, which may compromise measurement reliability. To mitigate this issue and minimize repeated X-ray irradiation, careful optimization of the X-ray incident angle and the use of a sufficiently large step size are recommended.
4. Conclusions
This study establishes synchrotron-based µ-GIWAXS as a transformative advancement over conventional GIWAXS by mitigating footprint effects, thereby achieving superior spectral resolution and measurement accuracy. By systematic tuning of the grazing incidence angle, we elucidated the relationship between footprint length and diffraction peak broadening, demonstrating consistency with geometrical broadening models. The implementation of beam-direction-perpendicular fly scanning further enhanced data quality while markedly reducing localized X-ray exposure.
Through systematic irradiation experiments, we evaluated the distinct degradation pathways for MAPbI3 and Cs0.05FA0.85MA0.1PbI3 films arising from the synergistic effect of X-ray irradiation and environmental moisture/oxygen. The concurrent photodecomposition of PbI2 was also identified. Under the present experimental conditions (10 keV, 0.5° incidence), a nitrogen atmosphere effectively suppresses radiation damage at absorbed dose rates up to 106 Gy s−1.
Furthermore, µ-GIWAXS enabled high-resolution spatial mapping of perovskite films, revealing key characteristics such as crystallinity variations, preferential orientations, PbI2 and lattice strain distributions. These capabilities directly address structural heterogeneity challenges in solution-processed films, particularly for large-area devices. Our findings underscore µ-GIWAXS as an indispensable tool for advancing perovskite film characterization, supporting both performance optimization and scalable manufacturing.
5. Related literature
The following references, not cited in the main body of the paper, have been cited in the supporting information: Rolston et al. (2018
); Wang et al. (2019
).
Supporting information
Figs. S1 to S13, Tables S1 to S4 and Notes S1 and S2. DOI: https://doi.org/10.1107/S1600577526004303/vl5054sup1.pdf
Footnotes
‡These authors contributed equally to this work.
Acknowledgements
The authors would like to acknowledge beamlines BL03HB, BL14B1 and BL17UM at Shanghai Synchrotron Radiation Facility (SSRF) for providing beam time.
Conflict of interest
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
Data availability
The data that support the findings of this study are available from the corresponding author upon reasonable request.
Funding information
The following funding is acknowledged: National Natural Science Foundation of China (grant No. 12075303 to Xingyu Gao; grant No. 11675252 to Xingyu Gao; grant No. U1632265 to Xingyu Gao); Chinese Academy of Sciences, Shanghai Branch (grant No. JCYJ-SHFY-2022-002 to Xingyu Gao); Shanghai Science and Technology Development Foundation (grant No. 24ZR1474100 to Zhenhuang Su).
References
Aydin, E., Allen, T. G., De Bastiani, M., Razzaq, A., Xu, L., Ugur, E., Liu, J. & De Wolf, S. (2024). Science 383, eadh3849. CrossRef PubMed Google Scholar
Berejnov, V., Rubinstein, B., Melo, L. G. A. & Hitchcock, A. P. (2021). J. Synchrotron Rad. 28, 834–848. Web of Science CrossRef CAS IUCr Journals Google Scholar
Birkholz, M. (2006). Thin film analysis by X-ray scattering. John Wiley & Sons. Google Scholar
Chan, P. F., Luo, Y., Cheng, H., Wang, X., Zhao, N., Wu, J., Su, C.-J., Lin, J.-M., Zhang, H. & Lu, X. (2025). Small 21, e06340. CrossRef PubMed Google Scholar
Chen, B., Rudd, P. N., Yang, S., Yuan, Y. & Huang, J. (2019). Chem. Soc. Rev. 48, 3842–3867. CrossRef CAS PubMed Google Scholar
Dang, Z., Shamsi, J., Palazon, F., Imran, M., Akkerman, Q. A., Park, S., Bertoni, G., Prato, M., Brescia, R. & Manna, L. (2017). ACS Nano 11, 2124–2132. CrossRef CAS PubMed Google Scholar
Duan, Y., Chen, Y., Wu, Y., Liu, Z., Liu, S. & Peng, Q. (2024). Adv. Funct. Mater. 34, 2315604. CrossRef Google Scholar
Epp, J. (2016). Materials characterization using nondestructive evaluation (nde) methods edited by G. Hübschen, I. Altpeter, R. Tschuncky & H.-G. Herrmann, pp. 81–124: Woodhead Publishing. Google Scholar
Francisco, M. C., d, S., Rodrigo, S., Maria, G. D. G., Verônica, C. T., Ana, F. N. & Hélio, C. N. T. (2024). Energy Mater. 4, 400058. Google Scholar
Gao, H., Xiao, K., Lin, R., Zhao, S., Wang, W., Dayneko, S., Duan, C., Ji, C., Sun, H., Bui, A. D., Liu, C., Wen, J., Kong, W., Luo, H., Zheng, X., Liu, Z., Nguyen, H., Xie, J., Li, L., Saidaminov, M. I. & Tan, H. (2024). Science 383, 855–859. CrossRef CAS PubMed Google Scholar
Goldsmith, C. C. & Walker, G. A. (1983). Adv. X-ray Anal. 27, 229–238. Google Scholar
Guo, Z., Yuan, M., Chen, G., Liu, F., Lu, R. & Yin, W.-J. (2024). Adv. Sci. 11, 2305799. CrossRef Google Scholar
He, R., Wang, W., Yi, Z., Lang, F., Chen, C., Luo, J., Zhu, J., Thiesbrummel, J., Shah, S., Wei, K., Luo, Y., Wang, C., Lai, H., Huang, H., Zhou, J., Zou, B., Yin, X., Ren, S., Hao, X., Wu, L., Zhang, J., Zhang, J., Stolterfoht, M., Fu, F., Tang, W. & Zhao, D. (2023). Nature 618, 80–86. CrossRef CAS PubMed Google Scholar
Hexemer, A. & Müller-Buschbaum, P. (2015). IUCrJ 2, 106–125. Web of Science CrossRef CAS PubMed IUCr Journals Google Scholar
Holder, C. F. & Schaak, R. E. (2019). ACS Nano 13, 7359–7365. CrossRef CAS PubMed Google Scholar
Huang, Z., Bai, Y., Huang, X., Li, J., Wu, Y., Chen, Y., Li, K., Niu, X., Li, N., Liu, G., Zhang, Y., Zai, H., Chen, Q., Lei, T., Wang, L. & Zhou, H. (2023). Nature 623, 531–537. CrossRef CAS PubMed Google Scholar
Jiang, Z. (2015). J. Appl. Cryst. 48, 917–926. Web of Science CrossRef CAS IUCr Journals Google Scholar
Jiang, Z., Ren, F., Zhou, Q., Wang, B., Sun, Z., Li, M., Zhao, Z., Liu, Z. & Chen, W. (2025). Adv. Energy Mater. 15, 2500598. CrossRef Google Scholar
Kirmani, A. R., Durant, B. K., Grandidier, J., Haegel, N. M., Kelzenberg, M. D., Lao, Y. M., McGehee, M. D., McMillon-Brown, L., Ostrowski, D. P., Peshek, T. J., Rout, B., Sellers, I. R., Steger, M., Walker, D., Wilt, D. M., VanSant, K. T. & Luther, J. M. (2022). Joule 6, 1015–1031. CrossRef CAS Google Scholar
Kirmani, A. R. & Sellers, I. R. (2025). Joule 9, 101852. CrossRef Google Scholar
Kong, H., Zhou, H., Xu, Q., Liu, K., Zhang, K., Wang, X., Wang, W., Wang, Z., Wang, S., Wang, Y., Tang, L., Yu, F., He, J. & Wang, Q. (2025). J. Synchrotron Rad. 32, 234–245. Web of Science CrossRef CAS IUCr Journals Google Scholar
Li, B., Shen, T. & Yun, S. (2023). Mater. Horiz. 10, 13–40. CrossRef PubMed Google Scholar
Li, D., Sun, X., Zhang, Y., Guan, Z., Yue, Y., Wang, Q., Zhao, L., Liu, F., Wei, J. & Li, H. (2024). Adv. Sci. 11, 2401184. CrossRef Google Scholar
Li, G., Su, Z., Canil, L., Hughes, D., Aldamasy, M. H., Dagar, J., Trofimov, S., Wang, L., Zuo, W., Jerónimo-Rendon, J. J., Byranvand, M. M., Wang, C., Zhu, R., Zhang, Z., Yang, F., Nasti, G., Naydenov, B., Tsoi, W. C., Li, Z., Gao, X., Wang, Z., Jia, Y., Unger, E., Saliba, M., Li, M. & Abate, A. (2023). Science 379, 399–403. CrossRef CAS PubMed Google Scholar
Li, N., Pratap, S., Körstgens, V., Vema, S., Song, L., Liang, S., Davydok, A., Krywka, C. & Müller-Buschbaum, P. (2022). Nat. Commun. 13, 6701. CrossRef PubMed Google Scholar
Li, P., Liang, C., Bao, B., Li, Y., Hu, X., Wang, Y., Zhang, Y., Li, F., Shao, G. & Song, Y. (2018). Nano Energy 46, 203–211. CrossRef CAS Google Scholar
Liang, Z., Zhang, Y., Xu, H., Chen, W., Liu, B., Zhang, J., Zhang, H., Wang, Z., Kang, D.-H., Zeng, J., Gao, X., Wang, Q., Hu, H., Zhou, H., Cai, X., Tian, X., Reiss, P., Xu, B., Kirchartz, T., Xiao, Z., Dai, S., Park, N.-G., Ye, J. & Pan, X. (2023). Nature 624, 557–563. CrossRef CAS PubMed Google Scholar
Lilliu, S., Dane, T. G., Alsari, M., Griffin, J., Barrows, A. T., Dahlem, M. S., Friend, R. H., Lidzey, D. G. & Macdonald, J. E. (2016). Adv. Funct. Mater. 26, 8221–8230. CrossRef CAS Google Scholar
Lin, W.-C., Lo, W.-C., Li, J.-X., Wang, Y.-K., Tang, J.-F. & Fong, Z.-Y. (2021). NPJ Mater. Degradation 5, 13. Google Scholar
Nozaki, H., Matsuo, H., Hayashi, Y., Kadoura, H. & Fujita, S. (2020). Mater. Lett. 277, 128383. CrossRef Google Scholar
Park, J., Kim, J., Yun, H.-S., Paik, M. J., Noh, E., Mun, H. J., Kim, M. G., Shin, T. J. & Seok, S. I. (2023). Nature 616, 724–730. CrossRef CAS PubMed Google Scholar
Pflüger, M., Soltwisch, V., Probst, J., Scholze, F. & Krumrey, M. (2017). IUCrJ 4, 431–438. Web of Science CrossRef PubMed IUCr Journals Google Scholar
Qin, M., Chan, P. F. & Lu, X. (2021). Adv. Mater. 33, 2105290. CrossRef Google Scholar
Rolston, N., Bush, K. A., Printz, A. D., Gold–Parker, A., Ding, Y., Toney, M. F., McGehee, M. D. & Dauskardt, R. H. (2018). Adv. Energy Mater. 8, 1802139. CrossRef Google Scholar
Sáez-Martínez, P., Gárate-Lagos, J., Camargo, S., Torres-Roquer, F., Queralt, I. & Salazar-Kuri, U. (2021). Forensic Sci. Int. 328, 111054. PubMed Google Scholar
Sakhatskyi, K., Turedi, B., Matt, G. J., Wu, E., Sakhatska, A., Bartosh, V., Lintangpradipto, M. N., Naphade, R., Shorubalko, I., Mohammed, O. F., Yakunin, S., Bakr, O. M. & Kovalenko, M. V. (2023). Nat. Photon. 17, 510–517. CrossRef CAS Google Scholar
Smilgies, D.-M. (2009). J. Appl. Cryst. 42, 1030–1034. Web of Science CrossRef CAS IUCr Journals Google Scholar
Song, J., Zhou, G., Chen, W., Zhang, Q., Ali, J., Hu, Q., Wang, J., Wang, C., Feng, W., Djurišić, A. B., Zhu, H., Zhang, Y., Russell, T. & Liu, F. (2020). Adv. Mater. 32, 2002784. CrossRef Google Scholar
Song, Y.-H., Li, B., Wang, Z.-J., Tai, X.-L., Ding, G.-J., Li, Z.-D., Xu, H., Hao, J.-M., Song, K.-H., Feng, L.-Z., Hu, Y.-L., Yin, Y.-C., Zhu, B.-S., Zhang, G., Ju, H., Zheng, G., Hu, W., Lin, Y., Fan, F. & Yao, H.-B. (2025). Nature 641, 352–357. CrossRef CAS PubMed Google Scholar
Steele, J. A., Solano, E., Hardy, D., Dayton, D., Ladd, D., White, K., Chen, P., Hou, J., Huang, H., Saha, R. A., Wang, L., Gao, F., Hofkens, J., Roeffaers, M. B. J., Chernyshov, D. & Toney, M. F. (2023). Adv. Energy Mater. 13, 2300760. Web of Science CrossRef Google Scholar
Stuckelberger, M. E., Nietzold, T., West, B. M., Luo, Y., Li, X., Werner, J., Niesen, B., Ballif, C., Rose, V., Fenning, D. P. & Bertoni, M. I. (2020). J. Phys. Chem. C 124, 17949–17956. CrossRef CAS Google Scholar
Svanström, S., García Fernández, A., Sloboda, T., Jacobsson, T. J., Rensmo, H. & Cappel, U. B. (2021). Phys. Chem. Chem. Phys. 23, 12479–12489. PubMed Google Scholar
Tu, Y., Wu, J., Xu, G., Yang, X., Cai, R., Gong, Q., Zhu, R. & Huang, W. (2021). Adv. Mater. 33, 2006545. CrossRef Google Scholar
Wang, C., He, B., Hui, W., Su, Z., Chen, L., Du, W., Zhang, L., Zhang, J., Fu, M., Wang, H., Sun, B., Wen, W., Wang, S., Lou, L., Hou, C., Zheng, G., Song, L., Chen, Y. & Gao, X. (2024). Adv. Funct. Mater. 34, 2403690. CrossRef Google Scholar
Wang, H., Zhu, C., Liu, L., Ma, S., Liu, P., Wu, J., Shi, C., Du, Q., Hao, Y., Xiang, S., Chen, H., Chen, P., Bai, Y., Zhou, H., Li, Y. & Chen, Q. (2019). Adv. Mater., 31, 1904408. CrossRef Google Scholar
Wang, L., Song, Y., Wang, J., Bi, W., Ding, L., Liu, H., Yang, X., Wang, Y., Yuan, S., Dong, Q., Yang, D. & Fang, Y. (2024). Appl. Mater. Interfaces 16, 44202–44209. CrossRef CAS Google Scholar
Wang, Z.-J., Wang, S.-S., Su, Z.-H., Yu, L., Wang, Y.-Z., Sun, B., Wen, W. & Gao, X.-Y. (2024). Nucl. Sci. Tech. 35, 119. CrossRef Google Scholar
Wei, H. & Huang, J. (2019). Nat. Commun. 10, 1066. CrossRef PubMed Google Scholar
Wei, Q., Zheng, D., Liu, L., Liu, J., Du, M., Peng, L., Wang, K. & Liu, S. (2024). Adv. Mater. 36, 2406295. CrossRef Google Scholar
Williamson, G. K. & Hall, W. H. (1953). Acta Metall. 1, 22–31. CrossRef CAS Web of Science Google Scholar
Xing-Yu, G., Tie-Ying, Y., Wen, W., Guang-Zhi, Y., Xiao-Long, L., Mei, G., Yue-Liang, G., Li, L., Yi, L., He, L., Xing-Min, Z., Bin, Z., Ting-Kun, L., Ying-Guo, Y., Zhong, L. & Xing-Tai, Z. (2015). Nucl. Sci. Tech. 26, 020101. Google Scholar
Xiong, Q., Huang, X., Wang, C., Zhou, Q., Gang, Y., Li, T., Hu, C., Zhang, N., Wang, X., Wu, J., Su, Z., Gao, X., Li, X., Zheng, N. & Gao, P. (2024). Joule 8, 817–834. CrossRef CAS Google Scholar
Xue, D.-J., Hou, Y., Liu, S.-C., Wei, M., Chen, B., Huang, Z., Li, Z., Sun, B., Proppe, A. H., Dong, Y., Saidaminov, M. I., Kelley, S. O., Hu, J.-S. & Sargent, E. H. (2020). Nat. Commun. 11, 1514. CrossRef PubMed Google Scholar
Yalcinkaya, Y., Hermes, I. M., Seewald, T., Amann–Winkel, K., Veith, L., Schmidt–Mende, L. & Weber, S. A. L. (2022). Adv. Energy Mater. 12, 2202442. CrossRef Google Scholar
Yin, R., Wu, Y., Huang, Z., Vasenko, A. S., Xu, S. & Zhou, H. (2025). Adv. Funct. Mater. 35, 2419184. CrossRef Google Scholar
Yoneda, Y. (1963). Phys. Rev. 131, 2010–2013. CrossRef Web of Science Google Scholar
Zhang, L., Fu, M., Jiang, X., Zhang, Z., Wang, C., Su, Z., He, B., Tang, L., Zheng, G., Gao, X. & He, J. (2025). Adv. Sci. 12, e12117. CrossRef Google Scholar
Zhang, S. & Han, G. (2020). Prog. Energy 2, 022002. CrossRef Google Scholar
Zhang, S., Ye, F., Wang, X., Chen, R., Zhang, H., Zhan, L., Jiang, X., Li, Y., Ji, X., Liu, S., Yu, M., Yu, F., Zhang, Y., Wu, R., Liu, Z., Ning, Z., Neher, D., Han, L., Lin, Y., Tian, H., Chen, W., Stolterfoht, M., Zhang, L., Zhu, W.-H. & Wu, Y. (2023). Science 380, 404–409. CrossRef CAS PubMed Google Scholar
Zhao, Y., Ma, F., Qu, Z., Yu, S., Shen, T., Deng, H.-X., Chu, X., Peng, X., Yuan, Y., Zhang, X. & You, J. (2022). Science 377, 531–534. CrossRef CAS PubMed Google Scholar
Zhou, Y., Herz, L. M., Jen, A. K. Y. & Saliba, M. (2022). Nat. Energy 7, 794–807. CrossRef CAS Google Scholar
Zhu, C., Niu, X., Fu, Y., Li, N., Hu, C., Chen, Y., He, X., Na, G., Liu, P., Zai, H., Ge, Y., Lu, Y., Ke, X., Bai, Y., Yang, S., Chen, P., Li, Y., Sui, M., Zhang, L., Zhou, H. & Chen, Q. (2019). Nat. Commun. 10, 815. CrossRef PubMed Google Scholar
Zhu, H., Teale, S., Lintangpradipto, M. N., Mahesh, S., Chen, B., McGehee, M. D., Sargent, E. H. & Bakr, O. M. (2023). Nat. Rev. Mater. 8, 569–586. CrossRef Google Scholar
This is an open-access article distributed under the terms of the Creative Commons Attribution (CC-BY) Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original authors and source are cited.

journal menu
access



