Difference between revisions of "Reflectivity"

From GISAXS
Jump to: navigation, search
(References)
 
(22 intermediate revisions by the same user not shown)
Line 3: Line 3:
 
'''[[X-ray]] Reflectivity''' ('''XRR''' or '''XR''') and '''[[neutron]] reflectivity''' ('''NR''') are techniques which measure the intensity of reflected radiation as a function of angle (where, by definition for '''specular reflectivity''', the incident and exit angles are equal; <math>\scriptstyle \theta_i = \theta_f</math>). A plot of reflectivity (''R'') versus angle yields the '''reflectivity curve'''. For XR and NR, the data is typically plotted as a function of the [[momentum transfer]] parallel to the film normal:
 
'''[[X-ray]] Reflectivity''' ('''XRR''' or '''XR''') and '''[[neutron]] reflectivity''' ('''NR''') are techniques which measure the intensity of reflected radiation as a function of angle (where, by definition for '''specular reflectivity''', the incident and exit angles are equal; <math>\scriptstyle \theta_i = \theta_f</math>). A plot of reflectivity (''R'') versus angle yields the '''reflectivity curve'''. For XR and NR, the data is typically plotted as a function of the [[momentum transfer]] parallel to the film normal:
 
:<math>
 
:<math>
q_z = \frac{4\pi}{\lambda} \sin \theta
+
q_z = \frac{4\pi}{\lambda} \sin \theta_i
 
</math>
 
</math>
  
 
[[Image:Reflectivity example-real density.png|thumb|center|400px|Reflectivity calculation for a thin film.]]
 
[[Image:Reflectivity example-real density.png|thumb|center|400px|Reflectivity calculation for a thin film.]]
 
==Off-Specular Reflectivity==
 
TBD
 
  
 
==Mathematical form==
 
==Mathematical form==
 +
===Ideal interface===
 
In its simplest form, the Fresnel reflectivity can be given by:
 
In its simplest form, the Fresnel reflectivity can be given by:
 
:<math>
 
:<math>
R_F = \frac{k_z-\tilde{k}_z}{k_z+\tilde{k}_z}
+
r_F = \frac{k_z-\tilde{k}_z}{k_z+\tilde{k}_z}
 
</math>
 
</math>
Where:
+
Where <math>\scriptstyle k_z</math> is the vertical (''z'') component of the momentum of the incident wave, and <math>\scriptstyle \tilde{k}_z</math> is the vertical component of the transmitted wave:
 
:<math>
 
:<math>
 
\tilde{k}_z = - \sqrt{ n^2 k^2_0 - |k_{\parallel}|^2 }
 
\tilde{k}_z = - \sqrt{ n^2 k^2_0 - |k_{\parallel}|^2 }
 
</math>
 
</math>
And ''n'' is the complex [[refractive index]] of the substrate. The idealized uncorrelated roughness can be characterized by a mean standard deviation of the height <math>\scriptstyle \sigma = \sqrt{ \left \langle h^2 \right \rangle }</math>:
+
and ''n'' is the complex [[refractive index]] of the substrate. The reflected intensity (which is what is measured experimentally) is <math>R_F = |r_f|^2</math>. Close to the critical angle (<math>\scriptstyle \theta_i \approx \theta_c</math>), the beam is strongly refracted and:
 +
:<math>
 +
\tilde{k}_z = k \sqrt{ \theta_i - \theta_c }
 +
</math>
 +
Thus, when the incident angle is less than the critical angle (<math>\scriptstyle \theta_i < \theta_c</math>), the transmitted wave-vector is purely imaginary:
 +
\tilde{k}_z = i k \sqrt{ \theta_c - \theta_i }
 +
This corresponds to the wave being totally externally reflected (evanescent probing only). Above the critical angle, the reflected intensity falls off rapidly: <math> R_F \sim \theta_i^{-4}</math> (c.f. [[Fresnel plot]]).
 +
 
 +
===Rough interface===
 +
A flat interface with a small amount of roughness can be modelled using idealized uncorrelated roughness, which is characterized by a mean standard deviation of the height <math>\scriptstyle \sigma = \sqrt{ \left \langle h^2 \right \rangle }</math>:
 
:<math>
 
:<math>
R_S = R_F e^{ -2 \sigma^2 k_z \tilde{k}_z }
+
r_S = r_F e^{ -2 \sigma^2 k_z \tilde{k}_z }
 
</math>
 
</math>
 
For a substrate with a single continuous layer of thickness ''h'' (e.g. a uniform thin film), the reflectivity becomes:
 
For a substrate with a single continuous layer of thickness ''h'' (e.g. a uniform thin film), the reflectivity becomes:
 
:<math>
 
:<math>
R_S = \frac{  r_{01} + r_{12} e^{ 2i \tilde{k}^1_z h }  }{  1 + r_{01} r_{12} e^{ 2i \tilde{k}^1_z h }  }
+
r_S = \frac{  r_{01} + r_{12} e^{ 2i \tilde{k}^1_z h }  }{  1 + r_{01} r_{12} e^{ 2i \tilde{k}^1_z h }  }
 
</math>
 
</math>
where <math>\tilde{k}^j_z=-\sqrt(n^2_j k^2_0 - |k_{\parallel}|^2 }</math> is the perpendicular component of the wave-vector (in medium ''j''). The reflectivity coefficients are:
+
where
 +
:<math>\tilde{k}^j_z = - \sqrt{ n^2_j k^2_0 - |k_{\parallel}|^2 }</math>
 +
is the perpendicular component of the wave-vector (in medium ''j''). The reflectivity coefficients are:
 
:<math>
 
:<math>
 
f_{ij} = \frac{ \tilde{k}^i_z - \tilde{k}^j_z  }{ \tilde{k}^i_z + \tilde{k}^j_z }
 
f_{ij} = \frac{ \tilde{k}^i_z - \tilde{k}^j_z  }{ \tilde{k}^i_z + \tilde{k}^j_z }
 
</math>
 
</math>
Where <math>\scriptstyle r_{01}</math> and <math>\scriptstyle r_{12}</math> for the [[vacuum]]-layer and layer-substrate interfaces, respectively. This is called the 'one-box model'.
+
Where <math>\scriptstyle r_{01}</math> and <math>\scriptstyle r_{12}</math> for the [[Material:Vacuum|vacuum]]-layer and layer-substrate interfaces, respectively. This is called the 'one-box model'.
 +
 
 +
===Multiple interfaces: Parratt formalism===
 +
The reflectivity from a stack of different layers can be computed using the so-called Parratt recursion formalism, which accounts for the effect of reflection from each internal interface, including their coherent interference to yield the final reflection intensity. This formalism can, in fact, be used to simulate or fit an arbitrary electron-density distribution (in the film-normal direction). Any continuous variation of density can be simulated by discretizing the system into a finite number of layers of finite thickness.
 +
 
 +
The Parratt formalism is recursive. The reflected and transmitted amplitude of layer ''j'' are given by <math>\scriptstyle R_j</math> and <math>\scriptstyle T_j</math>, such that:
 +
:<math>
 +
X_j = \frac{R_j}{T_j} = e^{ - 2i k_{z,j} z_j } \frac{  r_{j,j+1} + X_{j+1} e^{2 i k_{z,j} z_j}  }{  1 + r_{j,j+1} X_{j+1} e^{2 i k_{z,j} z_j}  }
 +
</math>
 +
where
 +
:<math>
 +
r_{j,j+1} = \frac{ k_{z,j} - k_{z,j+1} }{ k_{z,j} + k_{z,j+1} }
 +
</math>
 +
is the Fresnel coefficient for interface <math>\scriptstyle j</math>. The resultant system of equations can be solved by setting <math>T_1 = 1</math> (i.e. the incident wave is normalized to unity) and <math>R_{N+1} = 0</math> (i.e. there is no reflection from the bottom of the substrate).
 +
 
 +
The interference between various reflections generically leads to oscillations in the reflectivity curve (so-called '''Kiessig fringes'''). A one-box model (uniform thin film on a substrate) will give rise to a set of fringes with a single periodicity, which is related to the film thickness ''h'' [[q value|via]] <math> \Delta q_z = 2 \pi / h</math>. Multiple layers of different thickness will give fringes of different periodicities, superimposed in the reflectivity curve. This can also give rise to beat-patterns in the curve. A multi-layer stack with a repeated thickness (e.g. alternating layers of two materials, always with the same thickness) will interfere constructively to give a strong peak at a ''q<sub>z</sub>'' corresponding to the layering distance. This is called a [[Bragg peak]], and can also be thought of as the [[Structure Factor]] of lamellar ordering.
 +
 
 +
There are many available [[Software#Reflectivity|software packages]] which implement the Parratt formalism. These packages can be used to fit experimental data, and thereby obtain a possible (thought not unambiguous) curve for the [[realspace]] average density ([[SLD]]) profile in the film-normal direction.
 +
 
 +
* L. G. Parratt [http://journals.aps.org/pr/abstract/10.1103/PhysRev.95.359 Surface Studies of Solids by Total Reflection of X-Rays] ''Phys. Rev.'' '''1954''', 95, 359. [http://dx.doi.org/10.1103/PhysRev.95.359 doi: 10.1103/PhysRev.95.359]
 +
 
 +
===Master-Formula===
 +
For sufficiently large incident angle (<math>\scriptstyle \theta_i \gg \theta_c</math>), a [[kinematical appoximation]] can be valid. In such a case, the [[electron-density distribution|electron-density profile]] <math>\scriptstyle \rho_e(z)</math> can be described using the so-called Master-Formula:
 +
:<math>
 +
R(q_z) = R_F \left|  \frac{1}{\rho_e(z \to \infty)}  \int \frac{\mathrm{d}\rho_e}{\mathrm{d}z} e^{i q_z z} \mathrm{d}z \right|^2
 +
</math>
 +
This formalism breaks down near the [[critical angle]], but can be useful for predicting the reflectivity at higher angles, since the closed-form equation is simple to evaluate. It should be noted that this is essentially a [[Fourier transform]] of the density variation. In fact, this form highlights that the reflectivity comes from the ''variation'' of the material's refractive index (only points with non-zero <math>\mathrm{d}\rho_e/\mathrm{d}z</math> contribute to the integral). This makes plain that the reflectivity comes from locations where the refractive-index changes abruptly (at interfaces), or gradually.
 +
* P. Fenter, F. Schreiber, V. Bulovic, and S. R. Forrest [http://www.sciencedirect.com/science/article/pii/S000926149700941X Thermally induced failure mechanism of organic light emitting device structures probed by X-ray specular reflectivity] ''Chem. Phys. Lett.'' '''1997''', 277, 521. [http://dx.doi.org/10.1016/S0009-2614(97)00941-X doi: 10.1016/S0009-2614(97)00941-X]
 +
* I.W. Hamley and J.S. Pedersen [http://scripts.iucr.org/cgi-bin/paper?pii=S0021889893006260 Analysis of neutron and X-ray reflectivity. I. Theory] ''J. Appl. Cryst.'' '''1994''', 27, 29-35. [http://dx.doi.org/10.1107/S0021889893006260 doi: 10.1107/S0021889893006260]
 +
* J. S. Pedersen and I. W. Hamley [http://scripts.iucr.org/cgi-bin/paper?gl0299 Analysis of neutron and X-ray reflectivity. II. Constrained least-squares methods] ''J. Appl. Cryst.'' '''1994''', 27, 36-49. [http://dx.doi.org/10.1107/S0021889893006272 doi: 10.1107/S0021889893006272]
 +
 
 +
==Off-Specular Reflectivity==
 +
By definition, specular reflectivity involves setting the incident and detection angles to be equal, <math>\scriptstyle \theta_i = \theta_f</math>. This inherently probes exactly along the ''q<sub>z</sub>'' axis in [[reciprocal-space]]. One can also measure the '''off-specular reflectivity'''; i.e., the intensity of the reflectivity signal when the incident and exit angles are not quite equal. This can be accomplished either by holding the incident angle fixed, and rocking the point detector (or, equivalently, using a 1D detector); or this can be accomplished by holding the detector fixed and rocking the sample or input angle. In all cases, one is measuring the intensity in reciprocal-space ''near'' (but not exactly along) the ''q<sub>z</sub>'' axis. Thus, off-specular reflectivity directly connects to other [[scattering]] geometries. E.g. the [[specular rod]] seen in [[GISAXS]] is off-specular reflectivity (in GI, the specular reflection is a single point on the reflectivity curve; data near this point is thus off-specular). Note that during a specular reflectivity scan, one may want to measure and subtract a [[background]] scan from the data. A typical way to obtain a background is to perform a corresponding scan with a small offset of the sample or detector. Thus, in specular reflectivity, the 'background' is actually an off-specular signal.
 +
 
 +
While the specular reflectivity (''q<sub>z</sub>'' axis) encodes information about the structure of the sample in the film-normal direction, the off-specular data encodes information about deviations from this idealized structure. I.e. it contains information about in-plane correlations/order. This may be as simple as probing the lengthscale of roughness, or even probing well-defined in-plane structures.
 +
 
 +
* R. M. Richardson, J. R. P. Webster and A. Zarbakhsh [http://scripts.iucr.org/cgi-bin/paper?gl0511 Study of Off-Specular Neutron Reflectivity Using a Model System] ''J. Appl. Cryst.'' '''1997''', 30, 943-947. [http://dx.doi.org/10.1107/S0021889897003440 doi: 10.1107/S0021889897003440]
 +
 
 +
==References==
 +
===Modeling===
 +
* L. G. Parratt [http://journals.aps.org/pr/abstract/10.1103/PhysRev.95.359 Surface Studies of Solids by Total Reflection of X-Rays] ''Phys. Rev.'' '''1954''', 95, 359. [http://dx.doi.org/10.1103/PhysRev.95.359 doi: 10.1103/PhysRev.95.359]
 +
* I.W. Hamley and J.S. Pedersen [http://scripts.iucr.org/cgi-bin/paper?pii=S0021889893006260 Analysis of neutron and X-ray reflectivity. I. Theory] ''J. Appl. Cryst.'' '''1994''', 27, 29-35. [http://dx.doi.org/10.1107/S0021889893006260 doi: 10.1107/S0021889893006260]
 +
* J. S. Pedersen and I. W. Hamley [http://scripts.iucr.org/cgi-bin/paper?gl0299 Analysis of neutron and X-ray reflectivity. II. Constrained least-squares methods] ''J. Appl. Cryst.'' '''1994''', 27, 36-49. [http://dx.doi.org/10.1107/S0021889893006272 doi: 10.1107/S0021889893006272]
 +
* Krassimir Stoev and Kenji Sakurai [http://rigaku.com/downloads/journal/Vol14.2.1997/stoev.pdf Recent theoretical models in grazing incidence x-ray reflectometry] ''The Rigaku Journal'' '''1997''', 14, 22-37.
 +
 
 +
===Data Fitting===
 +
* Nelson, A. [http://www.iucr.org/cgi-bin/paper?ce5001 Co-refinement of multiple contrast neutron / X-ray reflectivity data using MOTOFIT] ''Journal of Applied Crystallography'' '''2006''', 39, 273-276. [http://dx.doi.org/10.1107/S0021889806005073 doi: 10.1107/S0021889806005073]
 +
* Charles F. Laub and Tonya L. Kuhl [http://scitation.aip.org/content/aip/journal/jcp/125/24/10.1063/1.2403126 Fitting a free-form scattering length density profile to reflectivity data using temperature-proportional quenching] ''J. Chem. Phys.'' '''2006''', 125, 244702. [http://dx.doi.org/10.1063/1.2403126 doi: 10.1063/1.2403126]
 +
* A. van der Lee, F. Salah and B. Harzallah [http://scripts.iucr.org/cgi-bin/paper?S0021889807032207 A comparison of modern data analysis methods for X-ray and neutron specular reflectivity data] ''J. Appl. Cryst.'' '''2007''', 40, 820-833. [http://dx.doi.org/10.1107/S0021889807032207 doi: 10.1107/S0021889807032207]
 +
* S. M. Danauskas, D. Li, M. Meron, B. Lin and K. Y. C. Lee [http://scripts.iucr.org/cgi-bin/paper?S0021889808032445 Stochastic fitting of specular X-ray reflectivity data using ''StochFit''] ''J. Appl. Cryst.'' '''2008''', 41, 1187-1193. [http://dx.doi.org/10.1107/S0021889808032445 doi: 10.1107/S0021889808032445]
  
===Parratt formalism===
+
===Data Treatment===
TBD
+
*  Arijeet Das, Shreyashkar Dev Singh, R. J. Choudhari, S. K. Rai and Tapas Ganguli  [http://scripts.iucr.org/cgi-bin/paper?rg5148 Data-reduction procedure for correction of geometric factors in the analysis of specular X-ray reflectivity of small samples] ''J. Appl. Cryst.'' '''2018''' [https://doi.org/10.1107/S1600576718010579 doi: 10.1107/S1600576718010579]
  
 +
===Reviews===
 +
* T.P. Russell [http://www.sciencedirect.com/science/article/pii/S0920230705800027 X-ray and neutron reflectivity for the investigation of polymers] ''Materials Science Reports'' '''1990''', 5 (4), 171–271. [http://dx.doi.org/10.1016/S0920-2307(05)80002-7 doi: 10.1016/S0920-2307(05)80002-7]
 +
* H. Zabel [http://link.springer.com/article/10.1007%2FBF00324371 X-ray and neutron reflectivity analysis of thin films and superlattices] ''Applied Physics A'' '''1994''', 58 (3), 159-168. [http://dx.doi.org/10.1007%2FBF00324371 doi: 10.1007%2FBF00324371]
 +
* A. Gibaud, S. Hazra [http://www.researchgate.net/profile/Satyajit_Hazra/publication/232747386_X-ray_reflectivity_and_diffuse_scatterin/links/0912f50d1a7d2e520d000000.pdf X-ray reflectivity and diffuse scattering] ''Current Science'' '''2000''', [http://www.currentscience.ac.in/php/toc.php?vol=078&issue=12 78 (12)], 1467-1477.
  
 
==See Also==
 
==See Also==
 +
* [[Software#Reflectivity|Reflectivity Software]]
 
* [[Fresnel plot]]
 
* [[Fresnel plot]]
 
* [[Reflectivity_oscillations_below_critical_angle|Oscillations below the critical angle]]
 
* [[Reflectivity_oscillations_below_critical_angle|Oscillations below the critical angle]]
 
* [[DWBA]]: A formalism for modeling [[GISAXS]] data, including reflection effects.
 
* [[DWBA]]: A formalism for modeling [[GISAXS]] data, including reflection effects.
 +
* Frank Schreiber, Alexander Gerlach [http://neutrons.ornl.gov/mr/pdf/beamline_04A_xray_neutron_reflectivity_tutorial.pdf X-Ray and Neutron Reflectivity for the Investigation of Thin Films – A Short Tutorial]

Latest revision as of 12:28, 29 August 2018

Reflectivity refers to the measurement of the intensity of reflection off of a flat interface. The term both describes the physical phenomenon, as well as the experimental technique.

X-ray Reflectivity (XRR or XR) and neutron reflectivity (NR) are techniques which measure the intensity of reflected radiation as a function of angle (where, by definition for specular reflectivity, the incident and exit angles are equal; ). A plot of reflectivity (R) versus angle yields the reflectivity curve. For XR and NR, the data is typically plotted as a function of the momentum transfer parallel to the film normal:

Reflectivity calculation for a thin film.

Mathematical form

Ideal interface

In its simplest form, the Fresnel reflectivity can be given by:

Where is the vertical (z) component of the momentum of the incident wave, and is the vertical component of the transmitted wave:

and n is the complex refractive index of the substrate. The reflected intensity (which is what is measured experimentally) is . Close to the critical angle (), the beam is strongly refracted and:

Thus, when the incident angle is less than the critical angle (), the transmitted wave-vector is purely imaginary: \tilde{k}_z = i k \sqrt{ \theta_c - \theta_i } This corresponds to the wave being totally externally reflected (evanescent probing only). Above the critical angle, the reflected intensity falls off rapidly: (c.f. Fresnel plot).

Rough interface

A flat interface with a small amount of roughness can be modelled using idealized uncorrelated roughness, which is characterized by a mean standard deviation of the height :

For a substrate with a single continuous layer of thickness h (e.g. a uniform thin film), the reflectivity becomes:

where

is the perpendicular component of the wave-vector (in medium j). The reflectivity coefficients are:

Where and for the vacuum-layer and layer-substrate interfaces, respectively. This is called the 'one-box model'.

Multiple interfaces: Parratt formalism

The reflectivity from a stack of different layers can be computed using the so-called Parratt recursion formalism, which accounts for the effect of reflection from each internal interface, including their coherent interference to yield the final reflection intensity. This formalism can, in fact, be used to simulate or fit an arbitrary electron-density distribution (in the film-normal direction). Any continuous variation of density can be simulated by discretizing the system into a finite number of layers of finite thickness.

The Parratt formalism is recursive. The reflected and transmitted amplitude of layer j are given by and , such that:

where

is the Fresnel coefficient for interface . The resultant system of equations can be solved by setting (i.e. the incident wave is normalized to unity) and (i.e. there is no reflection from the bottom of the substrate).

The interference between various reflections generically leads to oscillations in the reflectivity curve (so-called Kiessig fringes). A one-box model (uniform thin film on a substrate) will give rise to a set of fringes with a single periodicity, which is related to the film thickness h via . Multiple layers of different thickness will give fringes of different periodicities, superimposed in the reflectivity curve. This can also give rise to beat-patterns in the curve. A multi-layer stack with a repeated thickness (e.g. alternating layers of two materials, always with the same thickness) will interfere constructively to give a strong peak at a qz corresponding to the layering distance. This is called a Bragg peak, and can also be thought of as the Structure Factor of lamellar ordering.

There are many available software packages which implement the Parratt formalism. These packages can be used to fit experimental data, and thereby obtain a possible (thought not unambiguous) curve for the realspace average density (SLD) profile in the film-normal direction.

Master-Formula

For sufficiently large incident angle (), a kinematical appoximation can be valid. In such a case, the electron-density profile can be described using the so-called Master-Formula:

This formalism breaks down near the critical angle, but can be useful for predicting the reflectivity at higher angles, since the closed-form equation is simple to evaluate. It should be noted that this is essentially a Fourier transform of the density variation. In fact, this form highlights that the reflectivity comes from the variation of the material's refractive index (only points with non-zero contribute to the integral). This makes plain that the reflectivity comes from locations where the refractive-index changes abruptly (at interfaces), or gradually.

Off-Specular Reflectivity

By definition, specular reflectivity involves setting the incident and detection angles to be equal, . This inherently probes exactly along the qz axis in reciprocal-space. One can also measure the off-specular reflectivity; i.e., the intensity of the reflectivity signal when the incident and exit angles are not quite equal. This can be accomplished either by holding the incident angle fixed, and rocking the point detector (or, equivalently, using a 1D detector); or this can be accomplished by holding the detector fixed and rocking the sample or input angle. In all cases, one is measuring the intensity in reciprocal-space near (but not exactly along) the qz axis. Thus, off-specular reflectivity directly connects to other scattering geometries. E.g. the specular rod seen in GISAXS is off-specular reflectivity (in GI, the specular reflection is a single point on the reflectivity curve; data near this point is thus off-specular). Note that during a specular reflectivity scan, one may want to measure and subtract a background scan from the data. A typical way to obtain a background is to perform a corresponding scan with a small offset of the sample or detector. Thus, in specular reflectivity, the 'background' is actually an off-specular signal.

While the specular reflectivity (qz axis) encodes information about the structure of the sample in the film-normal direction, the off-specular data encodes information about deviations from this idealized structure. I.e. it contains information about in-plane correlations/order. This may be as simple as probing the lengthscale of roughness, or even probing well-defined in-plane structures.

References

Modeling

Data Fitting

Data Treatment

Reviews

See Also