research papers
Probing local atomic structure in thin films by grazing-incidence total X-ray scattering with a convergent X-ray beam
aX-ray Research Laboratory, Rigaku Corporation, Matsubara-cho 3-9-12, Akishima, Tokyo, Japan
*Correspondence e-mail: [email protected]
Grazing-incidence total X-ray scattering (GI-TXS) is an important technique for evaluating the short-range order and medium-range order structures of thin films. A high-quality pair distribution function requires high-intensity total scattering data over a wide scattering vector magnitude (Q) range, which necessitates the use of high-energy X-rays. However, a lower incident angle leads to a wider footprint at the sample position, resulting in a trade-off between achieving high-intensity total scattering and using high-energy X-rays in conventional laboratory setups. In this study, we performed GI-TXS measurements on an indium tin oxide (ITO) thin film using a convergent incident X-ray beam to maximize the X-ray flux at the sample position. By incorporating the convergence angle distribution and thin film absorption parameters, the substrate contribution to the observed total scattering intensity was quantitatively estimated and corrected. As a result, the total scattering profile of the ITO thin film was successfully obtained up to Q = 21 Å−1 using a laboratory X-ray diffractometer.
Keywords: grazing-incidence total X-ray scattering; GI-TXS; convergent incident X-ray beams; thin film local structure; pair distribution functions.
1. Introduction
Thin films are key materials used in a wide range of fields, including electronic devices, optical devices, catalysts and energy storage. Parameters such as film thickness, crystallinity, strain etc. determine the physical properties of thin films, such as electrical conductivity and transmittance. These properties directly affect the performance of practical devices. The atomic-scale structure of thin films can be single crystalline or polycrystalline, and these states have been used for device manufacturing. There have been many studies evaluating the crystalline structure of thin films using Rietveld refinement applied to grazing-incidence X-ray diffraction (GI-XRD) data, beginning with the first report by Quaas et al. (1998
). Compared with conventional powder XRD (i.e. 2θ/θ measurements), GI-XRD enhances the thin film signal relative to the substrate signal by limiting the X-ray penetration depth near the critical angle. Several studies have reported that amorphous and partially crystalline thin films can be used in practical devices (Schropp et al., 2007
; Tsuruma et al., 2021
; Büschges et al., 2025
; Büschges et al., 2026
). In non-crystalline materials, device properties depend on short-range order (SRO) and medium-range order (MRO) structures rather than long-range order. Therefore, it is important to establish appropriate experimental methods to quantitatively evaluate SRO and MRO structures of thin films.
Pair distribution function (PDF) analysis enables the evaluation of the SRO and MRO of materials, and it has been widely applied to various research fields: batteries (Bréger et al., 2007
; Ohara et al., 2016
), ferroelectric materials (Petkov et al., 2006
; Yoneda & Kohara, 2009
) and oxide glasses (Onodera et al., 2020
; Hashimoto et al., 2022
). It is necessary to measure total scattering data up to a sufficiently high momentum transfer, typically on the order of Q > 20 Å−1, because the maximum measurable momentum transfer Qmax determines the real-space resolution of peaks in PDF profiles (Peterson et al., 2003
). To obtain total scattering intensity with a high Qmax for thin film materials, several studies have reported thin film total scattering measurements performed using high-energy synchrotron radiation (Dippel et al., 2019
; Dippel et al., 2020
; Hoffman et al., 2024
; Zakutayev et al., 2024
). In the case of laboratory instruments, grazing-incidence total X-ray scattering (GI-TXS) techniques are necessary to obtain the structure factor S(Q) over a wide scattering vector magnitude (Q) range (Matsubara et al., 1988
; Bylin et al., 2024
). Bylin et al. have demonstrated GI-TXS measurement using a collimated incident X-ray beam and an absorption correction scheme based on a quasi-parallel beam approximation. With a parallel beam, typical incident angles in grazing-incidence geometry for high-energy X-rays are significantly smaller than 1°, leading to beam footprints on the sample that are larger by orders of magnitude. For example, the footprint size is approximately 380 times larger than the incident beam height when θin = 0.15°, and the brilliance of the incident beam is very low. A focused incident beam can solve this problem and enhance the incident flux on the sample to enable high-quality data collection, but it introduces a large convergence angle of the beam, which significantly affects the incident angle at laboratory sources. As a result, the incident angle can no longer be represented by a single value, which complicates the absorption correction in grazing-incidence geometry. Therefore, we propose a method for calculating the absorption factor of a film–substrate system by explicitly incorporating the non-linear dependence on the incident angle distribution.
In this paper, we present a GI-TXS measurement procedure and a data processing approach using a convergent incident beam to maximize the incident X-ray flux at the sample position. Furthermore, we estimate the substrate contribution to the observed total scattering intensity using the thickness, density and refractive index of thin film materials to correct for X-ray penetration effects of the convergent beam. We demonstrate GI-TXS measurements and obtain the PDF for an indium tin oxide (ITO) thin film deposited on a glass substrate to validate the present method.
2. Experimental
2.1. Theory
To apply GI-TXS to thin film samples, it is essential to account for X-ray behaviour in matter. The X-ray refractive index of a material, which is slightly less than 1, can be calculated using
where δ and β are defined as follows:
is the classical electron radius (2.818 × 10−15 m), λ is the X-ray wavelength (Å),
N0 is Avogadro's number (mol−1), ρ is the mass density (g cm−3),
ci is the molar fraction of the ith atom,
Zi is the atomic number of the ith atom,
Mi is the atomic mass of the ith atom, and
and
are the anomalous dispersion terms of the atomic scattering factor of the ith atom. As shown in equation (1
), the X-ray refractive index is a complex number. Here, δ and β mainly represent the refraction and absorption effects of X-rays in matter, respectively.
For a parallel beam, the absorption factor for a thin film, including refraction effects in matter (Matsubara et al., 1988
), is given by
where t is the thickness of the thin film, θin is the incident angle to the thin film and Λ is the X-ray penetration depth (Omote, 2002
):
For a convergent beam focused at the sample position, the thin film is illuminated by a range of incident angles within (
is half of the convergence angle of the optics). Therefore, the absorption factor depends not only on refraction in matter but also on the convergence angle distribution, as follows:
where is the weight factor for the convergent beam,
is the convergence angle within
and
is the bin size. In this study,
and
are 0.14° and 0.001°, respectively.
is defined as
where is the beam profile intensity at
.
Finally, we obtain the total scattering intensity of the thin film from the following equation:
where is the total scattering intensity from both the thin film and substrate, and
is the total scattering intensity contributed only by the substrate. The normalized intensity
is obtained from
, where a is a scaling coefficient determined using the Krogh-Moe–Norman formula (Krogh-Moe, 1956
; Norman, 1957
):
This normalization ensures that the resulting structure factor satisfies the correct asymptotic behaviour at high Q and is consistent with the physical normalization condition. Here, is the Compton scattering intensity referenced from tabulated values, ρ is the number density of the sample, and
and
are defined as follows (Egami & Billinge, 2012
; Bylin et al., 2024
):
where ci and fi are the molar concentration and atomic scattering factor of the ith element, respectively.
Q is the magnitude of the scattering vector, defined as
where is the scattering angle and λ is the X-ray wavelength. The structure factor
is then obtained from
as follows:
The PDF [] is obtained by Fourier transforming
:
2.2. Total scattering and X-ray reflectivity measurements
GI-TXS measurement was performed using a SmartLab high-resolution X-ray diffractometer (Rigaku Corporation) equipped with a HyPix-3000HE detector (Rigaku Corporation), which is optimized for high-energy X-ray detection. Incident Ag Kα (λ = 0.5608 Å) X-rays were monochromated by an elliptic d-space graded multilayer mirror focused on the sample position (Rigaku Innovative Technologies). The incident X-ray distribution of the convergent beam was characterized with a Si(12 0 0) crystal at the sample position as shown in Fig. 1
. The FWHM is about 0.19°. An ITO thin film (GEOMATEC Co., Ltd), whose thickness is designed to be 90 nm on a glass substrate, was prepared as a sample. Total scattering data were collected at a grazing-incidence angle of 0.14° over a scattering angle range of 2.0° ≤ 2θ ≤ 140.0°, corresponding to Qmin = 0.39 Å−1 and Qmax = 21.06 Å−1. The angular step size was 0.05°, with a counting time of 15 s per step. The thin film total scattering profile is extracted by applying equations (4)–(7), and other intensity corrections are performed using established procedures (Thijsse, 1984
).
|
Figure 1
Convergent beam profile characterized with a Si(12 0 0) crystal (i.e. ω = 38.173°) at the sample position. |
To evaluate the film parameters of ITO, an X-ray reflectivity measurement was also performed on a SmartLab. A parallel Cu Kα1 (λ = 1.5405 Å) beam was obtained using a parabolic d-space graded multilayer mirror and a Ge(220) double-crystal monochromator. The ITO thin film density, thickness and roughness were determined by X-ray reflectivity (XRR) analysis (plugin in SmartLab Studio II software, Rigaku Corporation).
3. Results and discussion
3.1. GI-TXS measurement for SiO2 glass plate
To validate the GI-TXS measurement, an SiO2 glass plate is one of the simplest samples because it does not require an absorption factor to subtract the background profile (i.e. air scattering). Fig. 2
shows a comparison of S(Q) obtained from GI-TXS and Debye–Scherrer geometry using an SiO2 glass plate with a thickness of 1 mm and an SiO2 glass rod with a diameter of 0.5 mm as a sample. The S(Q) from GI-TXS is consistent with that from Debye–Scherrer geometry over a wide Q range. This result indicates that GI-TXS measurements can be applied to evaluate the local structure in thin films.
|
Figure 2
Comparison of SiO2 glass S(Q) obtained from GI-TXS (red) and Debye–Scherrer geometry (black). |
3.2. ITO thin film analysis
Fig. 3
shows the observed and calculated XRR profile after refinement. The corresponding refinement results indicate that the ITO thin film can be described as a single layer with a density of 6.71 g cm−3, a film thickness of 76.6 nm and a surface roughness of 2.0 nm, as summarized in Table 1
. These refined parameters were subsequently used to calculate the absorption factor for correcting the substrate contribution.
|
||||||||||||||||||||||
|
Figure 3
Comparison of X-ray reflectivity profiles for the ITO thin film (top) and corresponding deviation curve (bottom). The red solid line shows the observed intensity. Blue filled circles show the calculated profile obtained using the refined ITO thin film parameters. The reliability factor R = |
Fig. 4
compares the extracted total scattering intensity of the thin film obtained using two different absorption correction schemes. The convergent beam correction incorporates the experimentally determined incident X-ray beam distribution shown in Fig. 1
, whereas the conventional parallel beam correction assumes a single incident angle (θin = 0.14°), as commonly adopted in previous studies and shown here as a reference. The discrepancy between the two results suggests that the single-angle approximation is insufficient under the present experimental conditions, where the incident angle distribution is broad due to beam convergence. The observed intensity Iobs indicates a broad peak related to the glass substrate at around 2θ = 10°. The absorption correction taking into consideration the convergent beam effect successfully extracts the ITO thin film intensity Itf from Iobs. On the other hand, the absorption correction for the parallel beam fails to extract the true film signal, as Itf still contains a broad peak originating from the glass substrate. Fig. 5
shows the observed data after polarization correction, normalized to the atomic scattering factor. Itf is consistent with the atomic scattering factor over the entire Q range. These results indicate that an absorption correction is required to consider the contribution from the beam profile at each convergence angle.
|
Figure 4
Comparison of absorption correction results: (top) absorption correction for a convergent beam using equations (5) and (6); (bottom) absorption correction for a parallel beam using equation (3) |
|
Figure 5
Total scattering intensity of the ITO thin film normalized by the sum profile of the atomic scattering factor 〈 f 2〉 and the Compton scattering intensity Iinc. Red line: total scattering intensity of the ITO thin film; green line: sum profile of 〈 f 2〉 and Iinc; black solid line: 〈 f 2〉; black dashed line: Iinc. |
Fig. 6
shows a comparison of S(Q) for the ITO thin film and powder In2O3. In the high-Q region (Q > 8 Å−1), the frequency and magnitude of the oscillations of the ITO thin film, which are related to the local structure of ITO, are consistent with those of powder In2O3. In contrast, the intensity ratio of reflections in the low-Q region 2 < Q < 4 Å−1 is quite different from that of powder In2O3. This suggests that ITO forms an oriented thin film. Rietveld analysis was performed on GI-TXS profiles of the ITO thin film. The observed GI-TXS profile is well reproduced by the calculated profile obtained from Rietveld refinement, as shown in Fig. S1 in the supporting information. The preferred orientation was quantified using the March–Dollase function, with the [001] direction selected as the preferred orientation vector based on the ITO(004) reflection. The refined March–Dollase parameter of 0.3931 (12), which is significantly smaller than unity, is characteristic of the plate-like crystallites and indicates a preferred orientation along the c-axis direction.
|
Figure 6
Comparison of S(Q) for the ITO thin film (red) and powder In2O3 (black). The inset also shows a zoomed-in view over 8 < Q < 21 Å−1. For clarity, an offset has been applied to the S(Q) of powder In2O3. |
The obtained G(r) of the ITO thin film is different from that of powder In2O3, as shown in Fig. 7
. Since the intensity ratio of the reflections differs from that of the powder In2O3, we conclude that the ITO forms an oriented thin film on the glass substrate. To estimate the partial correlation contributions, the Faber–Ziman weighting factor (Faber & Ziman, 1965
) of ITO for X-ray G(r) is assumed below:
where M and O denote metal and oxygen, respectively. The Faber–Ziman weighting factors are determined by both atomic concentrations and the atomic numbers in the X-ray G(r). Since the Sn concentration is much lower than that of In in typical ITO, the contribution of Sn-related partial correlations to the total G(r) is relatively small. Therefore, as a practical approximation, In and Sn partial correlations are treated as a single metal correlation. The equation indicates that M–M and M–O partial correlations mainly contribute to the observed G(r). At r < 4 Å, the peak positions of the ITO thin film are close to those of powder In2O3. In contrast, the G(r) peak shapes of the ITO thin film show different features compared with those of powder In2O3 in the region of r > 6 Å. These tendencies are consistent with previous reports on the comparison of G(r) between textured thin film materials and the homogeneous state (Dippel et al., 2019
; Roelsgaard et al., 2019
; Hoffman et al., 2024
). Therefore, the observed differences in the peak intensities are interpreted as being consistent with the presence of preferred orientation in the ITO thin film. A more detailed study of the texture effect in G(r) is left for future work.
|
Figure 7
Comparison of G(r) between ITO thin film (red) and powder In2O3 (black). |
4. Conclusion
In this study, we propose a GI-TXS measurement protocol and data analysis process using a convergent beam at the sample position. The thin film S(Q) can be obtained by correcting reflection and absorption effects using the fundamental physical parameters of the thin film and the X-ray refractive indices. The proposed method enables reliable measurements of thin film S(Q) up to Qmax = 21 Å−1 using a laboratory diffractometer. The improvement arises from the non-linear absorption factor that accounts for the incident angle distribution of the convergent beam, thereby overcoming the limitations of the conventional single-angle approximation in the grazing-incidence geometry. It is confirmed that the obtained S(Q) of the ITO thin film in the high-Q region agrees well with that of the powder In2O3 sample in both the frequency and amplitude of the oscillations. This is reasonable because the short-range structure is expected to be similar in both the thin film and powder states, and the oscillations in this region reflect the short-range structure. However, it should be noted that, in the present out-of-plane total X-ray scattering measurements, the scattering profile contains contributions from different reciprocal-lattice vectors (G vector) for the incident and scattered X-rays at each 2θ position. Therefore, the resulting S(Q) represents a mixture of different Q-vector directions. In the future, we will attempt to measure the total scattering profile along the same Q vector via in-plane measurements. Despite this limitation, the present GI-TXS method can be widely applied to evaluate the SRO and MRO structures of amorphous and partially crystalline thin films.
Supporting information
Supporting information: Fig. S1 which shows the Rietveld refinement results. DOI: https://doi.org/10.1107/S1600576726007181/jur5014sup1.pdf
Acknowledgements
The authors acknowledge Dr Katsuhiko Inaba for valuable discussions about the basic theory of GI-XRD. The authors also acknowledge Dr Beverly Vincent, Rigaku, for his careful reading of the manuscript and suggestions.
Conflict of interest
The authors declare no conflicts of interest.
Data availability
The data are available upon request.
References
Bréger, J., Kang, K., Cabana, J., Ceder, G. & Grey, C. P. (2007). J. Mater. Chem. 17, 3167.
Google Scholar
Büschges, M. I., Dietz, C., Trouillet, V., Dippel, A., Igoa Saldaña, F. & Schneider, J. J. (2026). Adv. Elect Mater. 12, e00722.
Google Scholar
Büschges, M. I., Trouillet, V., Dippel, A. & Schneider, J. J. (2025). Adv. Mater. Int. 12, 2400758.
Google Scholar
Bylin, J., Kapaklis, V. & Pálsson, G. K. (2024). J. Appl. Cryst. 57, 1373–1383.
CrossRef
CAS
IUCr Journals
Google Scholar
Dippel, A.-C., Gutowski, O., Klemeyer, L., Boettger, U., Berg, F., Schneller, T., Hardtdegen, A., Aussen, S., Hoffmann-Eifert, S. & Zimmermann, M. V. (2020). Nanoscale 12, 13103–13112.
CrossRef
CAS
PubMed
Google Scholar
Dippel, A.-C., Roelsgaard, M., Boettger, U., Schneller, T., Gutowski, O. & Ruett, U. (2019). IUCrJ 6, 290–298.
CrossRef
CAS
PubMed
IUCr Journals
Google Scholar
Egami, T. & Billinge, S. J. L. (2012). Underneath the Bragg Peaks: Structural Analysis of Complex Materials, 2nd ed. Amsterdam: Elsevier.
Google Scholar
Faber, T. E. & Ziman, J. M. (1965). Philos. Mag. 11, 153–173.
CrossRef
CAS
Web of Science
Google Scholar
Hashimoto, H., Onodera, Y., Tahara, S., Kohara, S., Yazawa, K., Segawa, H., Murakami, M. & Ohara, K. (2022). Sci. Rep. 12, 516.
CrossRef
PubMed
Google Scholar
Hoffman, J. M., Thompson, N. B., Borkiewicz, O., He, X., Amsterdam, S., Xie, Z., -, Taggart, A., Mulfort, K. L., Martinson, A. B. F., Chen, L. X., Ruett, U. & Tiede, D. M. (2024). IUCrJ 11, 120–128.
CrossRef
CAS
PubMed
IUCr Journals
Google Scholar
Krogh-Moe, J. (1956). Acta Cryst. 9, 951–953.
CrossRef
CAS
IUCr Journals
Web of Science
Google Scholar
Matsubara, E., Waseda, Y. & Itozaki, H. (1988). Trans. JIM 29, 1–7.
CrossRef
CAS
Google Scholar
Norman, N. (1957). Acta Cryst. 10, 370–373.
CrossRef
CAS
IUCr Journals
Web of Science
Google Scholar
Ohara, K., Mitsui, A., Mori, M., Onodera, Y., Shiotani, S., Koyama, Y., Orikasa, Y., Murakami, M., Shimoda, K., Mori, K., Fukunaga, T., Arai, H., Uchimoto, Y. & Ogumi, Z. (2016). Sci. Rep. 6, 21302.
CrossRef
PubMed
Google Scholar
Omote, K. (2002). J. Surf. Anal. 9, 203–209.
CAS
Google Scholar
Onodera, Y., Kohara, S., Salmon, P. S., Hirata, A., Nishiyama, N., Kitani, S., Zeidler, A., Shiga, M., Masuno, A., Inoue, H., Tahara, S., Polidori, A., Fischer, H. E., Mori, T., Kojima, S., Kawaji, H., Kolesnikov, A. I., Stone, M. B., Tucker, M. G., McDonnell, M. T., Hannon, A. C., Hiraoka, Y., Obayashi, I., Nakamura, T., Akola, J., Fujii, Y., Ohara, K., Taniguchi, T. & Sakata, O. (2020). NPG Asia Mater. 12, 85.
Google Scholar
Peterson, P. F., Božin, E. S., Proffen, Th. & Billinge, S. J. L. (2003). J. Appl. Cryst. 36, 53–64.
Web of Science
CrossRef
CAS
IUCr Journals
Google Scholar
Petkov, V., Gateshki, M., Niederberger, M. & Ren, Y. (2006). Chem. Mater. 18, 814–821.
Web of Science
CrossRef
CAS
Google Scholar
Quaas, M., Eggs, C. & Wulff, H. (1998). Thin Solid Films 332, 277–281.
CrossRef
CAS
Google Scholar
Roelsgaard, M., Dippel, A.-C., Borup, K. A., Nielsen, I. G., Broge, N. L. N., Röh, J. T., Gutowski, O. & Iversen, B. B. (2019). IUCrJ 6, 299–304.
CrossRef
CAS
PubMed
IUCr Journals
Google Scholar
Schropp, R. E. I., Carius, R. & Beaucarne, G. (2007). MRS Bull. 32, 219–224.
CrossRef
CAS
Google Scholar
Thijsse, B. J. (1984). J. Appl. Cryst. 17, 61–76.
CrossRef
CAS
Web of Science
IUCr Journals
Google Scholar
Tsuruma, Y., Kawashima, E., Nagasaki, Y., Sekiya, T., Imamura, G. & Yoshikawa, G. (2021). Sci. Rep. 11, 9435.
CrossRef
PubMed
Google Scholar
Yoneda, Y. & Kohara, S. (2009). J. Korean Phys. Soc. 55, 741–745.
CrossRef
CAS
Google Scholar
Zakutayev, A., Jankousky, M., Wolf, L., Feng, Y., Rom, C. L., Bauers, S. R., Borkiewicz, O., LaVan, D. A., Smaha, R. W. & Stevanovic, V. (2024). Nat. Synth. 3, 1471–1480.
CrossRef
CAS
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.

menu
access