X-ray reflectivity

Last updated November 15, 2022

X-ray reflectivity (sometimes known as X-ray specular reflectivity, X-ray reflectometry, or XRR) is a surface-sensitive analytical technique used in chemistry, physics, and materials science to characterize surfaces, thin films and multilayers.^[1]^[2]^[3]^[4] It is a form of reflectometry based on the use of X-rays and is related to the techniques of neutron reflectometry and ellipsometry.

The basic principle of X-ray reflectivity is to reflect a beam of X-rays from a flat surface and to then measure the intensity of X-rays reflected in the specular direction (reflected angle equal to incident angle). If the interface is not perfectly sharp and smooth then the reflected intensity will deviate from that predicted by the law of Fresnel reflectivity. The deviations can then be analyzed to obtain the density profile of the interface normal to the surface.

History

The technique appears to have first been applied to X-rays by Lyman G. Parratt in 1954.^[5] Parratt's initial work explored the surface of copper-coated glass, but since that time the technique has been extended to a wide range of both solid and liquid interfaces.

Approximation

When an interface is not perfectly sharp, but has an average electron density profile given by $\rho _{e}(z)$ , then the X-ray reflectivity can be approximated by:^[2]^: 83

R(Q)/R_{F}(Q)=\left|{\frac {1}{\rho _{\infty }}}{\int \limits _{-\infty }^{\infty }{e^{iQz}\left({\frac {d\rho _{e}}{dz}}\right)dz}}\right|^{2}

Here $R(Q)$ is the reflectivity, $Q=4\pi \sin(\theta )/\lambda$ , $\lambda$ is the X-ray wavelength (typically copper's K-alpha peak at 0.154056 nm), $\rho _{\infty }$ is the density deep within the material and $\theta$ is the angle of incidence. Below the critical angle $Q<Q_{c}$ (derived from Snell's law), 100% of incident radiation is reflected through total external reflection, $R=1$ . For $Q\gg Q_{c}$ , $R\sim Q^{-4}$ . Typically one can then use this formula to compare parameterized models of the average density profile in the z-direction with the measured X-ray reflectivity and then vary the parameters until the theoretical profile matches the measurement.

Oscillations

For films with multiple layers, X-ray reflectivity may show oscillations with Q (angle/wavelength), analogous to the Fabry-Pérot effect, here called Kiessig fringes.^[6] The period of these oscillations can be used to infer layer thicknesses, interlayer roughnesses, electron densities and their contrasts, and complex refractive indices (which depend on atomic number and atomic form factor), for example using the Abeles matrix formalism or the recursive Parratt-formalism as follows:

X_{j}={\frac {R_{j}}{T_{j}}}={\frac {r_{j,j+1}+X_{j+1}e^{2ik_{j+1,z}d_{j}}}{1+r_{j,j+1}X_{j+1}e^{2ik_{j+1,z}d_{j}}}}e^{-2ik_{j,z}d_{j}}

where X_j is the ratio of reflected and transmitted amplitudes between layers j and j+1, d_j is the thickness of layer j, and r_j,j+1 is the Fresnel coefficient for layers j and j+1

r_{j,j+1}={\frac {k_{j,z}-k_{j+1,z}}{k_{j,z}+k_{j+1,z}}}

where k_j,z is the z component of the wavenumber. For specular reflection where the incident and reflected angles are equal, Q used previously is two times k_z because $Q=k_{\text{incident}}+k_{\text{reflected}}$ . With conditions R_N+1 = 0 and T₁ = 1 for an N-interface system (i.e. nothing coming back from inside the semi-infinite substrate and unit amplitude incident wave), all X_j can be calculated successively. Roughness can also be accounted for by adding the factor

r_{j,j+1,{\text{rough}}}=r_{j,j+1,{\text{ideal}}}e^{-2k_{j,z}k_{j+1,z}\sigma _{j}^{2}}

where $\sigma$ is a standard deviation (aka roughness).

Thin film thickness and critical angle can also be approximated with a linear fit of squared incident angle of the peaks $\theta ^{2}$ in rad² vs unitless squared peak number $N^{2}$ as follows:

\theta ^{2}=\left({\frac {\lambda }{2d}}\right)^{2}N^{2}+\theta _{c}^{2}

.

Curve fitting

X-ray reflectivity measurements are analyzed by fitting to the measured data a simulated curve calculated using the recursive Parratt's formalism combined with the rough interface formula. The fitting parameters are typically layer thicknesses, densities (from which the index of refraction $n$ and eventually the wavevector z component $k_{j,z}$ is calculated) and interfacial roughnesses. Measurements are typically normalized so that the maximum reflectivity is 1, but normalization factor can be included in fitting, as well. Additional fitting parameters may be background radiation level and limited sample size due to which beam footprint at low angles may exceed the sample size, thus reducing reflectivity.

Several fitting algorithms have been attempted for X-ray reflectivity, some of which find a local optimum instead of the global optimum. The Levenberg-Marquardt method finds a local optimum. Due to the curve having many interference fringes, it finds incorrect layer thicknesses unless the initial guess is extraordinarily good. The derivative-free simplex method also finds a local optimum. In order to find global optimum, global optimization algorithms such as simulated annealing are required. Unfortunately, simulated annealing may be hard to parallelize on modern multicore computers. Given enough time, simulated annealing can be shown to find the global optimum with a probability approaching 1,^[7] but such convergence proof does not mean the required time is reasonably low. In 1998,^[8] it was found that genetic algorithms are robust and fast fitting methods for X-ray reflectivity. Thus, genetic algorithms have been adopted by the software of practically all X-ray diffractometer manufacturers and also by open source fitting software.

Fitting a curve requires a function usually called fitness function, cost function, fitting error function or figure of merit (FOM). It measures the difference between measured curve and simulated curve, and therefore, lower values are better. When fitting, the measurement and the best simulation are typically represented in logarithmic space.

From mathematical standpoint, the $\chi ^{2}$ fitting error function takes into account the effects of Poisson-distributed photon counting noise in a mathematically correct way:

F=\sum _{i}{\frac {(x_{simul,i}-x_{meas,i})^{2}}{x_{meas,i}}}

.

However, this $\chi ^{2}$ function may give too much weight to the high-intensity regions. If high-intensity regions are important (such as when finding mass density from critical angle), this may not be a problem, but the fit may not visually agree with the measurement at low-intensity high-angle ranges.

Another popular fitting error function is the 2-norm in logarithmic space function. It is defined in the following way:

F={\sqrt {\sum _{i}(\log x_{simul,i}-\log x_{meas,i})^{2}}}

.

Needless to say, in the equation data points with zero measured photon counts need to be removed. This 2-norm in logarithmic space can be generalized to p-norm in logarithmic space. The drawback of this 2-norm in logarithmic space is that it may give too much weight to regions where relative photon counting noise is high.

Open source software

Diffractometer manufacturers typically provide commercial software to be used for X-ray reflectivity measurements. However, several open source software packages are also available: GenX^[9]^[10] is a commonly used open source X-ray reflectivity curve fitting software. It is implemented in the Python programming language and runs therefore on both Windows and Linux. Motofit^[11]^[12] runs in the IGOR Pro environment, and thus cannot be used in open-source operating systems such as Linux. Micronova XRR^[13] runs under Java and is therefore available on any operating system on which Java is available. Reflex^[14]^[15] is a standalone software dedicated to the simulation and analysis of X-rays and neutron reflectivity from multilayers. REFLEX is a user-friendly freeware program working under Windows, Mac and Linux platforms.

Related Research Articles

In 3D computer graphics, ray tracing is a technique for modeling light transport for use in a wide variety of rendering algorithms for generating digital images.

X-ray crystallography is the experimental science determining the atomic and molecular structure of a crystal, in which the crystalline structure causes a beam of incident X-rays to diffract into many specific directions. By measuring the angles and intensities of these diffracted beams, a crystallographer can produce a three-dimensional picture of the density of electrons within the crystal. From this electron density, the mean positions of the atoms in the crystal can be determined, as well as their chemical bonds, their crystallographic disorder, and various other information.

<span class="mw-page-title-main">Reflectance</span> Capacity of an object to reflect light

The reflectance of the surface of a material is its effectiveness in reflecting radiant energy. It is the fraction of incident electromagnetic power that is reflected at the boundary. Reflectance is a component of the response of the electronic structure of the material to the electromagnetic field of light, and is in general a function of the frequency, or wavelength, of the light, its polarization, and the angle of incidence. The dependence of reflectance on the wavelength is called a reflectance spectrum or spectral reflectance curve.

A logarithmic spiral, equiangular spiral, or growth spiral is a self-similar spiral curve that often appears in nature. The first to describe a logarithmic spiral was Albrecht Dürer (1525) who called it an "eternal line". More than a century later, the curve was discussed by Descartes (1638), and later extensively investigated by Jacob Bernoulli, who called it Spira mirabilis, "the marvelous spiral".

The Phong reflection model is an empirical model of the local illumination of points on a surface designed by the computer graphics researcher Bui Tuong Phong. In 3D computer graphics, it is sometimes referred to as "Phong shading", particularly if the model is used with the interpolation method of the same name and in the context of pixel shaders or other places where a lighting calculation can be referred to as “shading”.

In physics and chemistry, Bragg's law, Wulff–Bragg's condition or Laue–Bragg interference, a special case of Laue diffraction, gives the angles for coherent scattering of waves from a crystal lattice. It encompasses the superposition of wave fronts scattered by lattice planes, leading to a strict relation between wavelength and scattering angle, or else to the wavevector transfer with respect to the crystal lattice. Such law had initially been formulated for X-rays upon crystals. However, It applies to all sorts of quantum beams, including neutron and electron waves at atomic distances, as well as visible light at artificial periodic microscale lattices.

Specular reflection, or regular reflection, is the mirror-like reflection of waves, such as light, from a surface.

The bidirectional reflectance distribution function is a function of four real variables that defines how light is reflected at an opaque surface. It is employed in the optics of real-world light, in computer graphics algorithms, and in computer vision algorithms. The function takes an incoming light direction, $, and outgoing direction,, and returns the ratio of reflected radiance exiting along to the irradiance incident on the surface from direction . Each direction is itself parameterized by azimuth angle and zenith angle, therefore the BRDF as a whole is a function of 4 variables. The BRDF has units sr -1, with steradians (sr) being a unit of solid angle.$

<span class="mw-page-title-main">Specular highlight</span>

A specular highlight is the bright spot of light that appears on shiny objects when illuminated. Specular highlights are important in 3D computer graphics, as they provide a strong visual cue for the shape of an object and its location with respect to light sources in the scene.

Rietveld refinement is a technique described by Hugo Rietveld for use in the characterisation of crystalline materials. The neutron and X-ray diffraction of powder samples results in a pattern characterised by reflections at certain positions. The height, width and position of these reflections can be used to determine many aspects of the material's structure.

<span class="mw-page-title-main">Bragg plane</span>

In physics, a Bragg plane is a plane in reciprocal space which bisects a reciprocal lattice vector, $, at right angles. The Bragg plane is defined as part of the Von Laue condition for diffraction peaks in x-ray diffraction crystallography.$

In condensed matter physics and crystallography, the static structure factor is a mathematical description of how a material scatters incident radiation. The structure factor is a critical tool in the interpretation of scattering patterns obtained in X-ray, electron and neutron diffraction experiments.

$<span class="mw-page-title-main">Fiber diffraction</span> Subarea of scattering, an area in which molecular structure is determined from scattering data$

Fiber diffraction is a subarea of scattering, an area in which molecular structure is determined from scattering data. In fiber diffraction the scattering pattern does not change, as the sample is rotated about a unique axis. Such uniaxial symmetry is frequent with filaments or fibers consisting of biological or man-made macromolecules. In crystallography fiber symmetry is an aggravation regarding the determination of crystal structure, because reflexions are smeared and may overlap in the fiber diffraction pattern. Materials science considers fiber symmetry a simplification, because almost the complete obtainable structure information is in a single two-dimensional (2D) diffraction pattern exposed on photographic film or on a 2D detector. 2 instead of 3 co-ordinate directions suffice to describe fiber diffraction.

Neutron reflectometry is a neutron diffraction technique for measuring the structure of thin films, similar to the often complementary techniques of X-ray reflectivity and ellipsometry. The technique provides valuable information over a wide variety of scientific and technological applications including chemical aggregation, polymer and surfactant adsorption, structure of thin film magnetic systems, biological membranes, etc.

In geometry, various formalisms exist to express a rotation in three dimensions as a mathematical transformation. In physics, this concept is applied to classical mechanics where rotational kinematics is the science of quantitative description of a purely rotational motion. The orientation of an object at a given instant is described with the same tools, as it is defined as an imaginary rotation from a reference placement in space, rather than an actually observed rotation from a previous placement in space.

<span class="mw-page-title-main">Plane mirror</span> Mirror with a flat reflecting surface

A plane mirror is a mirror with a flat (planar) reflective surface. For light rays striking a plane mirror, the angle of reflection equals the angle of incidence. The angle of the incidence is the angle between the incident ray and the surface normal. Therefore, the angle of reflection is the angle between the reflected ray and the normal and a collimated beam of light does not spread out after reflection from a plane mirror, except for diffraction effects.

Time-domain thermoreflectance is a method by which the thermal properties of a material can be measured, most importantly thermal conductivity. This method can be applied most notably to thin film materials, which have properties that vary greatly when compared to the same materials in bulk. The idea behind this technique is that once a material is heated up, the change in the reflectance of the surface can be utilized to derive the thermal properties. The reflectivity is measured with respect to time, and the data received can be matched to a model with coefficients that correspond to thermal properties.

<span class="mw-page-title-main">Transfer-matrix method (optics)</span>

The transfer-matrix method is a method used in optics and acoustics to analyze the propagation of electromagnetic or acoustic waves through a stratified medium. This is for example relevant for the design of anti-reflective coatings and dielectric mirrors.

<span class="mw-page-title-main">Grazing-incidence small-angle scattering</span>

Grazing-incidence small-angle scattering (GISAS) is a scattering technique used to study nanostructured surfaces and thin films. The scattered probe is either photons or neutrons. GISAS combines the accessible length scales of small-angle scattering and the surface sensitivity of grazing incidence diffraction (GID).

In a Fourier transform (FT), the Fourier transformed function $is obtained from by:$

References

↑ Holý, V.; Kuběna, J.; Ohlídal, I.; Lischka, K.; Plotz, W. (1993-06-15). "X-ray reflection from rough layered systems". Physical Review B. American Physical Society (APS). 47 (23): 15896–15903. Bibcode:1993PhRvB..4715896H. doi:10.1103/physrevb.47.15896. ISSN 0163-1829. PMID 10005989.
1 2 J. Als-Nielsen, D. McMorrow, Elements of Modern X-Ray Physics, Wiley, New York, (2001).
↑ J. Daillant, A. Gibaud, X-Ray and Neutron Reflectivity: Principles and Applications. Springer, (1999).
↑ M. Tolan, X-Ray Scattering from Soft-Matter Thin Films, Springer, (1999).
↑ Parratt, L. G. (1954-07-15). "Surface Studies of Solids by Total Reflection of X-Rays". Physical Review. American Physical Society (APS). 95 (2): 359–369. Bibcode:1954PhRv...95..359P. doi:10.1103/physrev.95.359. ISSN 0031-899X.
↑ Kiessig, Heinz (1931). "Untersuchungen zur Totalreflexion von Röntgenstrahlen". Annalen der Physik (in German). Wiley. 402 (6): 715–768. Bibcode:1931AnP...402..715K. doi:10.1002/andp.19314020607. ISSN 0003-3804.
↑ Granville, V.; Krivanek, M.; Rasson, J.-P. (1994). "Simulated annealing: a proof of convergence". IEEE Transactions on Pattern Analysis and Machine Intelligence. Institute of Electrical and Electronics Engineers (IEEE). 16 (6): 652–656. doi:10.1109/34.295910. ISSN 0162-8828.
↑ Dane, A.D.; Veldhuis, A.; Boer, D.K.G.de; Leenaers, A.J.G.; Buydens, L.M.C. (1998). "Application of genetic algorithms for characterization of thin layered materials by glancing incidence X-ray reflectometry". Physica B: Condensed Matter. Elsevier BV. 253 (3–4): 254–268. Bibcode:1998PhyB..253..254D. doi:10.1016/s0921-4526(98)00398-6. ISSN 0921-4526.
↑ Bjorck, Matts. "GenX - Home". genx.sourceforge.net.
↑ Björck, Matts; Andersson, Gabriella (2007-11-10). "GenX: an extensible X-ray reflectivity refinement program utilizing differential evolution". Journal of Applied Crystallography. International Union of Crystallography (IUCr). 40 (6): 1174–1178. doi:10.1107/s0021889807045086. ISSN 0021-8898.
↑ "Main Page - Motofit". motofit.sourceforge.net.
↑ Nelson, Andrew (2006-03-12). "Co-refinement of multiple-contrast neutron/X-ray reflectivity data using MOTOFIT". Journal of Applied Crystallography. International Union of Crystallography (IUCr). 39 (2): 273–276. doi:10.1107/s0021889806005073. ISSN 0021-8898.
↑ "jmtilli/micronovaxrr". GitHub. 2017-07-25.
↑ Vignaud, Guillaume; Gibaud, Alain (2019-02-01). "REFLEX: a program for the analysis of specular X-ray and neutron reflectivity data". Journal of Applied Crystallography. International Union of Crystallography (IUCr). 52 (1): 201–213. doi:10.1107/s1600576718018186. ISSN 1600-5767. S2CID 104467965.
↑ "Main Page - Reflex". reflex.irdl.fr/Reflex/reflex.html.

This page is based on this Wikipedia article
Text is available under the CC BY-SA 4.0 license; additional terms may apply.
Images, videos and audio are available under their respective licenses.

[1] Holý, V.; Kuběna, J.; Ohlídal, I.; Lischka, K.; Plotz, W. (1993-06-15). "X-ray reflection from rough layered systems". Physical Review B. American Physical Society (APS). 47 (23): 15896–15903. Bibcode:1993PhRvB..4715896H. doi:10.1103/physrevb.47.15896. ISSN 0163-1829. PMID 10005989.

[alsnielsen-2] 1 2 J. Als-Nielsen, D. McMorrow, Elements of Modern X-Ray Physics, Wiley, New York, (2001).

[3] J. Daillant, A. Gibaud, X-Ray and Neutron Reflectivity: Principles and Applications. Springer, (1999).

[4] M. Tolan, X-Ray Scattering from Soft-Matter Thin Films, Springer, (1999).

[5] Parratt, L. G. (1954-07-15). "Surface Studies of Solids by Total Reflection of X-Rays". Physical Review. American Physical Society (APS). 95 (2): 359–369. Bibcode:1954PhRv...95..359P. doi:10.1103/physrev.95.359. ISSN 0031-899X.

[6] Kiessig, Heinz (1931). "Untersuchungen zur Totalreflexion von Röntgenstrahlen". Annalen der Physik (in German). Wiley. 402 (6): 715–768. Bibcode:1931AnP...402..715K. doi:10.1002/andp.19314020607. ISSN 0003-3804.

[7] Granville, V.; Krivanek, M.; Rasson, J.-P. (1994). "Simulated annealing: a proof of convergence". IEEE Transactions on Pattern Analysis and Machine Intelligence. Institute of Electrical and Electronics Engineers (IEEE). 16 (6): 652–656. doi:10.1109/34.295910. ISSN 0162-8828.

[8] Dane, A.D.; Veldhuis, A.; Boer, D.K.G.de; Leenaers, A.J.G.; Buydens, L.M.C. (1998). "Application of genetic algorithms for characterization of thin layered materials by glancing incidence X-ray reflectometry". Physica B: Condensed Matter. Elsevier BV. 253 (3–4): 254–268. Bibcode:1998PhyB..253..254D. doi:10.1016/s0921-4526(98)00398-6. ISSN 0921-4526.

[9] Bjorck, Matts. "GenX - Home". genx.sourceforge.net.

[10] Björck, Matts; Andersson, Gabriella (2007-11-10). "GenX: an extensible X-ray reflectivity refinement program utilizing differential evolution". Journal of Applied Crystallography. International Union of Crystallography (IUCr). 40 (6): 1174–1178. doi:10.1107/s0021889807045086. ISSN 0021-8898.

[11] "Main Page - Motofit". motofit.sourceforge.net.

[12] Nelson, Andrew (2006-03-12). "Co-refinement of multiple-contrast neutron/X-ray reflectivity data using MOTOFIT". Journal of Applied Crystallography. International Union of Crystallography (IUCr). 39 (2): 273–276. doi:10.1107/s0021889806005073. ISSN 0021-8898.

[13] "jmtilli/micronovaxrr". GitHub. 2017-07-25.

[14] Vignaud, Guillaume; Gibaud, Alain (2019-02-01). "REFLEX: a program for the analysis of specular X-ray and neutron reflectivity data". Journal of Applied Crystallography. International Union of Crystallography (IUCr). 52 (1): 201–213. doi:10.1107/s1600576718018186. ISSN 1600-5767. S2CID 104467965.

[15] "Main Page - Reflex". reflex.irdl.fr/Reflex/reflex.html.

[1]

[2]

[3]

[4]

[5]

[6]

[7]

[8]

[9]

[10]

[11]

[12]

[13]

[14]

[15]