X-Ray and Neutron scattering 4

Small-Angle X-Ray Scattering (SAXS)

Overview and Foundations

• Definition: Small-Angle X-ray Scattering (SAXS) is an analytical technique used to determine the structural characteristics of materials at length scales ranging from $1\text{ nm}$ to $100\text{ nm}$.

• Key Historical Figures:

  • Otto Kratky ($1902 - 1995$): Invented the Kratky collimation system (Kratky camera) for low-angle measurement.

  • Günther Porod ($1919 - 1984$): Formulated Porod's Law to analyze internal surface areas.

  • Andre Guinier ($1911 - 2000$): Developed the Guinier approximation to determine the radius of gyration.

  • Otto Glatter ($1945 - \text{present}$): Developed mathematical inversion methods for pair distance distribution function (PDDF) analysis.

  • Reference Text: Small Angle X-ray Scattering, O. Glatter and O. Kratky, Academic Press, London ($1982$).

• Primary Applications:

  • Determination of size, shape, and internal morphology of nano-sized particles in dilute or concentrated media.

  • Measurement of specific internal surface area in porous or multiphase systems.

  • Analysis of long-range periodicities ($2\pi/d$) in materials with large lattice parameters, such as lyotropic liquid crystals, thermotropic liquid crystals, and colloidal crystals.

Theoretical Kinematics and Form Factor

Non-Correlated Particles

• Intensity Equation: For a sample containing $N$ non-correlated particles within the total sample volume:

    $I(\vec{q}) = N \cdot |E_{tot}(\vec{q})|^2$

    where: * $E_{tot}(\vec{q})$ is the intra-particle form factor (scattering amplitude of all electrons contained within a single coherence volume). * $N$ is the total number of scattering particles.

• Key Assumptions:

  • At most one particle resides within the coherence volume of the X-ray beam.

  • Inter-particle interference is absent; total scattered intensity is the incoherent sum of individual particle contributions.

Summation over Scattering Centers (Kinematical Theory)

• The total scattering amplitude $E_{tot}(\vec{q})$ is the sum of scattered waves from all electron density fluctuations within the coherence volume $V$:

    $E_{tot}(\vec{q}) = \text{const} \int_V d\vec{r}\,\Delta\rho(\vec{r}) e^{-i\vec{q}\vec{r}}$

• Electron Density Decomposition:

    $\rho(\vec{r}) = \rho_0 + \Delta\rho(\vec{r})$

    where: * $\rho_0$ is the average electron density of the surrounding solvent or matrix. * $\Delta\rho(\vec{r})$ is the local electron density fluctuation relative to the matrix. * Uniform density background (matrix) contributes only to unscattered forward transmission at $\vec{q} = 0$; only spatial variations $\Delta\rho(\vec{r})$ give rise to observable SAXS intensity.

Pair Distance Distribution Function (PDDF)

• Single Particle Scattering Amplitude:

    $E_S(\vec{q}) = \int_{V_S} d\vec{r}\,\Delta\rho(\vec{r}) e^{-i\vec{q}\vec{r}}$

• Spatial Autocorrelation Function $\gamma(\vec{r})$: Represents the spatial overlap of particle electron density with itself shifted by vector $\vec{r}$:

    $\gamma(\vec{r}) = \int d\vec{r}'\,\Delta\rho(\vec{r}') \Delta\rho(\vec{r}' - \vec{r})$

• Real-Space Intensity Relation:

    $I(\vec{q}) = \int d\vec{r}\,\gamma(\vec{r}) e^{-i\vec{q}\vec{r}}$

• Orientational Averaging (Isotropic Systems):

  • Setting $\vec{q} \cdot \vec{r} = q r \cos(\theta)$ and integrating over spherical coordinates ($d\vec{r} = r^2 \sin(\theta)\,dr\,d\theta\,d\varphi$):

        $I(q) = 4\pi \int_0^\infty dr\,p(r) \frac{\sin(qr)}{qr}$

*   **Pair Distance Distribution Function** $p(r)$: Defined as $p(r) = r^2 \gamma(r)$. The value $p(r)$ is directly proportional to the frequency of occurrence of internal interatomic/inter-electron distances within the interval $[r, r+dr]$.

• Inverse Transformation (The Inverse Scattering Problem): Calculates real-space structural parameters from measured scattering curves:

    $p(r) = \frac{1}{2\pi^2} \int_0^\infty dq\,I(q) q r \sin(qr)$

Form Factors for Specific Geometries

Homogeneous Sphere (Rayleigh, 1911)

• Scattering Intensity:

    $I(q) \propto \left[ 3 \frac{\sin(qR) - qR \cos(qR)}{(qR)^3} \right]^2$

    where $R$ is the sphere radius.

• Characteristic Minima: Sharp zeros occur in the scattering function at precise dimensionless values:

    $qR = 4.493, \quad 7.725, \quad 10.90$

Homogeneous Cylinder (Mittelbach & Porod, 1961)

• Scattering Intensity:

    $I(q) \propto \int_0^{\pi/2} d\alpha \left[ \frac{2 J_1\left(\frac{qd}{2} \sin(\alpha)\right)}{\frac{qd}{2} \sin(\alpha)} \cdot \frac{\sin\left(\frac{qL}{2} \cos(\alpha)\right)}{\frac{qL}{2} \cos(\alpha)} \right]^2 \sin(\alpha)$

    where: * $d$ is the cylinder diameter. * $L$ is the cylinder length. * $J_1$ is the first-order Bessel function of the first kind. * $\alpha$ is the angle between the cylinder axis and the scattering vector $\vec{q}$.

Monodisperse vs. Polydisperse Systems

• Monodisperse Systems: Uniform particle dimensions yield deep, sharply defined form factor minima.

• Polydisperse Systems: Size distributions cause superposition of out-of-phase oscillation patterns, resulting in progressive smearing and flattening of scattering minima.

Particle Parameters Derived from SAXS

Molecular Weight Determination

• Forward scattering intensity $I(0)$ is evaluated at $q = 0$:

    $I(0) = \left| \int_V d\vec{r}\,\Delta\rho(\vec{r}) \right|^2 = (n_e - n_0)^2$

    where: * $n_e$ is the total number of electrons in a single particle. * $n_0$ is the number of electrons in an equivalent volume of matrix/solvent.

• Knowing the stoichiometric elemental composition allows absolute conversion of $I(0)$ into particle molecular weight.

Radius of Gyration ($R_G$) & Guinier Approximation

• Definition: Radius of gyration ($R_G$, Streumassenradius) represents the root-mean-square distance of scattering centers from the electron density center of gravity:

    $R_G^2 = \frac{\int \vec{r}^2 \Delta\rho(\vec{r})\,d\vec{r}}{\int \Delta\rho(\vec{r})\,d\vec{r}}$

• Guinier Approximation Formula: Valid in the low-angle limit ($0 < q R_G < 1$):

    $I(q) = I(0) \exp\left( -\frac{q^2 R_G^2}{3} \right)$

• Guinier Plot: Plotting $\ln I(q)$ versus $q^2$ yields a straight line in the low-angle regime with:

  • $\text{Slope} = -\frac{R_G^2}{3}$

  • $y\text{-intercept} = \ln I(0)$

• Geometrical Relationships:

  • Sphere of radius $r$: $R_G = \sqrt{\frac{3}{5}} r$ (or $R_G = \frac{5}{3}r$ per structural approximation models).

  • Flat disc of thickness $t$: $R_G = t$

  • Needle/rod of diameter $c$: $R_G = c$

Porod's Law and Specific Surface Area

• Asymptotic Behavior: At high scattering angles ($q \rightarrow \infty$), scattering originates from abrupt boundaries between phases:

    $I(q) \propto q^{-4}$

• Porod Invariant $\widetilde{Q}$:

    $\widetilde{Q} = \int_0^\infty dq\,q^2 I(q)$

• Specific Internal Surface Area ($S/V$):

    $\frac{\pi}{\widetilde{Q}} \lim_{q \rightarrow \infty} \left[ I(q) q^4 \right] = \frac{S}{V}$

    Allows quantitative comparison of internal surface area per unit volume between distinct morphologies ($S_A/V > S_B/V$).

Correlated (Interacting) Particle Systems

• Inter-particle Correlation: In concentrated/saturated solutions or densely packed suspensions, inter-particle distance distributions distort the form factor profile.

• Factorization Model:

    $E(q) = S_p(q) \cdot E_S(q)$

    where: * $E_S(q)$ is the intra-particular form factor (particle geometry). * $S_p(q)$ is the inter-particular structure factor (spatial distribution/packing).

• Isotropic Hard Sphere Liquids: The radial pair distribution function $g(r)$ describes local structural correlation and liquid short-range order via Fourier transformation of $S_p(q)$.

Experimental Instrumentation for SAXS

Kratky Camera (Laboratory Source)

• Design Objective: Resolves weak scattering at tiny angles ($q \rightarrow 0$, low $2\theta$) without interference from the intense primary beam.

• Key Components:

  • Extreme collimation section defined by a precision slit $S$.

  • A block-collimator and beamstop system aligned to absorb the direct beam.

  • Evacuated flight tubes to eliminate background air scattering.

Synchrotron SAXS and Small-Angle Neutron Scattering (SANS)

• Synchrotron SAXS (e.g., Beamline ID02 at ESRF, Grenoble):

  • Sample-to-detector distance configurable from $0.8\text{ m}$ to $31\text{ m}$.

  • Achieves minimum scattering vector $q_{min} \approx 10^{-3}\text{ nm}^{-1}$, resolving maximum structures up to $d_{max} \approx 6\,\mu\text{m}$.

• Small-Angle Neutron Scattering (SANS) (e.g., ILL, Grenoble):

  • Detector flight tubes extending up to $40\text{ m}$.

Thin Films & Surfaces: Reciprocal Space Mapping

Classifications of Crystalline Order in Thin Films

• Amorphous: Lacks long-range translational or orientation order (e.g., $a\text{-Si}$, $a\text{-SiO}_2$).

• Random Polycrystalline: Randomly oriented micro-crystallites exhibiting full 3D powder averaging (e.g., sputtered metallic films).

• Textured: Preferred crystallographic out-of-plane orientation perpendicular to the substrate ($z$), but random azimuthal in-plane orientation ($xy$) (e.g., vapor-deposited organic semiconductors or metals).

• Mosaic: Preferred orientation with small angular distribution (tilt/twist) around ideal crystallographic axes (e.g., molecular films on rubbed polyimide).

• Perfect: Uniaxial single-crystalline epitaxy across the film (e.g., single-crystal GaN grown on GaAs).

Epitaxial Film Strain & Relaxation

• Lattice Matching Principle: Epitaxial growth requires matching substrate unit cell parameters:

    $\frac{|a - a_s|}{a} < 0.1\%$

    where $a$ is native film lattice parameter and $a_s$ is substrate lattice parameter.

• Strained Layer: Film matches substrate in-plane lattice parameter ($d_L^\parallel = a_s$), producing tetragonally distorted out-of-plane spacing ($d_L^\perp$).

• Relaxed Layer: Above critical thickness, misfit dislocations form, allowing the film to revert to its bulk native lattice constant a$.\n\n\n\n---\n\n## Reciprocal Space Broadening Mechanisms\n\n* **Finite Crystal Size**: Broadens Bragg nodes along both lateral (q_{xy}$) and vertical ($q_z) directions.\n* **Mosaicity**: Causes arc-like angular broadening perpendicular to the scattering vector \vec{q}.\n* **Lattice Parameter Variations** (\Delta d/d): Shifts and broadens reflection nodes along the radial scattering direction.\n\n---\n\n## Scanning Geometries in Reciprocal Space\n\n* **Coplanar Geometry**: Primary beam, surface normal, and diffracted beam remain within a single plane.\n* **Scan Modes**:\n * **Specular Scan** (\theta/2\Theta$or$2\theta/\omega$): Angles step synchronously ($\omega = \theta$). Probes electron density variations strictly normal to surface ($q_z).\n * **Rocking Curve** (\omega$-scan):$2\Theta$is fixed while sample angle$\omega$rotates. Measures mosaicity and lateral orientation distribution ($q_x).\n * **Detector Scan** (2\theta$-scan): Sample angle$\omega$is stationary while detector angle$2\Theta scans.\n\n\n\n---\n\n## Crystal Truncation Rods (CTRs)\n\n* **Physical Origin**: Abrupt termination of a periodic crystal lattice by a flat surface breaks 3D translational symmetry.\n* **Mathematical Formalism**: Fourier transform of a 3D semi-infinite crystal lattice multiplied by a step function \Theta(z).\n* **Reciprocal Space Signature**: Continuous lines of scattered intensity ( truncation rods ) extend along q_z perpendicular to the surface, passing through bulk Bragg points.\n* **Asymptotic Intensity Decay**: Intensity decays proportional to q_z^{-2} away from Bragg peaks.\n* **Structural Sensitivity**: CTR profiles are extremely sensitive to sub-angstrom surface features, including interface roughness, relaxation, and surface reconstruction (e.g., \text{Si}(100) reconstruction).\n\n# X-Ray Reflectivity (XRR)\n\n## Refractive Index of Materials for X-Rays\n\n* **Complex Index of Refraction**:\n\n    n = 1 - \delta + i\beta\n\n* **Refractive Index Decrement** \delta$: Describes phase shift velocity (typically$\delta \sim 10^{-5} - 10^{-6}):\n\n    \delta = \frac{\lambda^2}{2\pi} r_e \rho_e\n\n    where:\n * \lambda is X-ray wavelength.\n * r_e = \frac{e^2}{4\pi\varepsilon_0 m c^2} = 2.818 \times 10^{-15}\text{ m} is classical electron radius.\n * \rho_e is electron density of the solid.\n* **Absorption Index** \beta$: Describes linear photoelectric attenuation (typically$\beta \sim 10^{-7} - 10^{-8}):\n\n    \beta = \frac{\lambda}{4\pi} \mu_x\n\n    where \mu_x is linear absorption coefficient.\n\n---\n\n## Total External Reflection\n\n* Because n < 1$in condensed matter for hard X-rays, total external reflection occurs when X-rays impinge from vacuum/air ($n_1 = 1$) at grazing incident angles below a critical angle$\alpha_c$.

• Snell's Law at Grazing Incidence:

    $n_1 \cos(\alpha_i) = n_2 \cos(\alpha')$

• Critical Angle Derivation:

    $\cos(\alpha_c) = 1 - \delta$

    Applying Taylor expansion $\cos(\alpha_c) \approx 1 - \frac{\alpha_c^2}{2}$:

    $1 - \frac{\alpha_c^2}{2} = 1 - \delta \implies \alpha_c = \sqrt{2\delta}$

• Critical Scattering Vector $q_c$:

    $q_c = \frac{4\pi}{\lambda} \sin(\alpha_c)$

Material Parameters at Standard Characteristic Emission Wavelengths

Fresnel Reflectivity and Surface Roughness

Ideal Smooth Surface (Fresnel Formulas)

• For $s$-polarized light at grazing incidence ($\sin(\alpha_i) \approx \alpha_i$):

    $\frac{E_R^s}{E_0} = \frac{\sin(\alpha_i) - n_1 \sin(\alpha_1)}{\sin(\alpha_i) + n_1 \sin(\alpha_1)} \approx \frac{\alpha_i - \alpha_1}{\alpha_i + \alpha_1}$

    $\frac{E_T^s}{E_0} = \frac{2 \sin(\alpha_i)}{\sin(\alpha_i) + n_1 \sin(\alpha_1)} \approx \frac{2 \alpha_i}{\alpha_i + \alpha_1}$

• Reflectivity $R_F$ and Transmissivity $T_F$:

    $R_F = |r|^2 = \left| \frac{E_R}{E_0} \right|^2, \quad T_F = |t|^2 = \left| \frac{E_T}{E_0} \right|^2$

• Asymptotic Regimes:

  • For $q < q_c$: Total reflection plateau ($R_F \approx 1$).

  • For $q > 3 q_c$: Fresnel fall-off follows $R_F \propto q^{-4}$.

Real Rough Surfaces (Nevot-Croce Factor)

• Root-mean-square (RMS) surface roughness $\sigma_{RMS}$ dampens reflected amplitude via a Gaussian attenuation factor:

    $r_{rough} = r_{ideal} \exp\left( -\frac{1}{2} q^2 \sigma_{RMS}^2 \right)$

    $R_{real} = R_{ideal} \exp\left( -q^2 \sigma_{RMS}^2 \right) = C q^{-4} \exp\left( -q^2 \sigma_{RMS}^2 \right)$

Thin Film Interference: Kiessig Fringes

• Physical Mechanism: Constructive and destructive interference between waves reflected at the top surface and bottom substrate interface.

• Film Thickness Calculation: Spacing between consecutive fringe maxima ($\Delta q_z$) directly determines total film thickness $d$:

    $d \approx \frac{2\pi}{\Delta q_z}$

Reflectivity Data Analysis: Parratt Formalism

• Dynamical Recursion Scheme ($L.G. Parratt, 1954$): Exactly solves wave field propagation through multi-layered stratified media containing $N$ interfaces.

    $r_{total} = r_{01} + \frac{t_{01} t_{10} r_{12} p}{1 - r_{10} r_{12} p}$

    where: * $r_{ij}$ and $t_{ij}$ are Fresnel reflection and transmission coefficients at interface $ij$. * $p = e^{i \phi} = e^{i q_z d_1}$ is phase factor across layer thickness $d_1$

• Parameters Extracted via Model Fitting:

  1. Layer thickness $d$

  2. Surface RMS roughness $\sigma_{surf}$

  3. Interface RMS roughness $\sigma_{inter}$

  4. Layer mass density $\rho$ / electron density $\rho_e$

Exemplary Experimental Case Studies

• Thermally Oxidized Silicon Wafer ($\text{Si}$ + $\text{SiO}_2$ + surface water contamination):

  • $\text{SiO}_2\text{ layer}$: $d = 149.3 \pm 3.4\text{ nm}$, $\sigma = 0.33 \pm 0.02\text{ nm}$, $\rho_e = 0.674 \pm 0.054\text{ \AA}^{-3}$, \rho = 2.15 \pm 0.04\text{ g/cm}^3$.\n * \text{Top contamination layer (H}_2\text{O)}$:$d = 0.95 \pm 0.01\text{ nm}$,$\sigma = 0.44 \pm 0.03\text{ nm}$,$\rho_e = 0.377 \pm 0.041\text{ \AA}^{-3}$,$\rho = 1.12 \pm 0.1\text{ g/cm}^3$.

• Organic Polymer Thin Films (P3HT processed from different solvents):

  • Toluene-cast film: $\sigma_{surf} = 4.7\text{ nm}$, $\sigma_{inter} = 0.2\text{ nm}$.

  • Chloroform-cast film: $\sigma_{surf} = 0.9\text{ nm}$, $\sigma_{inter} = 1.2\text{ nm}$.

• Model-Independent Fitting: Uses discretized box profiles (e.g., 96 sub-boxes) to extract internal structural periodicities, such as the $d\text{-spacing} = 1.66\text{ nm}$ of P3HT polymer backbones.

Off-Specular (Diffuse) Reflectivity & Dynamic Effects

• Off-Specular Geometry ($\alpha_i \neq \alpha_f$): Generates non-zero in-plane scattering components ($q_x \neq 0$).

• Lateral Correlation Length $\xi_\parallel$: Characterizes in-plane surface height variations $h(x,y)$. Derived from diffuse peak width $\Delta q_x$:

    $\xi_\parallel = \frac{2\pi}{\Delta q_x}$

• Yoneda / Vineyard Peaks: Dynamic wavefield enhancements occurring in diffuse channels when either incident angle $\alpha_i = \alpha_c$ or exit scattering angle $\alpha_f = \alpha_c$ ($Y. Yoneda, 1963; G.H. Vineyard, 1982$).

Neutron Reflectivity (NR)

• Complementary Contrast: Uses nuclear scattering length density rather than electron density.

• Isotopic Substitution: Isotopic contrast matching (substituting Hydrogen with Deuterium, e.g., protonated vs. deuterated polymer layers) allows tracking interdiffusion at buried organic interfaces without altering chemical properties.

Grazing Incidence X-Ray Diffraction (GIXD)

Fundamental Principles

• Technique Description: Structural characterization tool optimized for thin films ($< 100\text{ nm}$ down to single monolayers).

• Evanescent Wavefield Propagation: Setting incident angle $\alpha_i \approx \alpha_c$ forms an evanescent wave at the interface that propagates parallel to the surface and decays exponentially in depth ($e^{-z/\Lambda}$).

• Penetration Depth $\Lambda$:

    $\Lambda = \frac{1}{k \text{Im}\left( \sqrt{\alpha_i^2 - 2\delta - 2i\beta} \right)}$

*   For $\alpha_i < \alpha_c$, penetration depth is restricted to $\Lambda \sim 5\text{ nm}$, eliminating background noise from the underlying bulk substrate.

• Absorption Factor $b_\mu$:

    $b_\mu = \frac{2 k \mu_x}{q_c^2}$

    Typical values at $8\text{ keV}$: Carbon ($0.001$), Silicon ($0.011$), Gold ($0.043$).

Wavevector Coordinates & Refraction Corrections

• Scattering Vector Coordinates:

    $q_{xy} = \frac{2\pi}{\lambda} \sqrt{\cos^2(\alpha_i) + \cos^2(\alpha_f) - 2\cos(\alpha_i)\cos(\alpha_f)\cos(\theta_f)}$

    $q_z \approx \frac{2\pi}{\lambda} \left( \sin(\alpha_i) + \sin(\alpha_f) \right)$

• Refraction Correction Formulas ($W.C. Marra, I. Robinson$):

    $k_{0z} = -\frac{2\pi}{\lambda} \sqrt{n^2 - \cos^2(\alpha_i)}$

    $k_z = \frac{2\pi}{\lambda} \sqrt{n^2 - \cos^2(\alpha_f)}$

    $q_z = k_z - k_{0z}$

Reciprocal Space Mapping and Structural Analysis

• Double Peaks: Arise from primary beam reflections occurring prior to or after Bragg scattering at the substrate interface.

• Polymorph & Phase Analysis: Identifies surface-induced crystal structures distinct from bulk phases (e.g., distinguishing pentacene surface-induced phase from bulk Campbell phase).

• Unit Cell Indexing: Combines out-of-plane specular peaks $(0,0,q_z)$ with in-plane GIXD peaks $(q_{xy}, q_z)$. For example, indexing 75 reflection nodes yields a triclinic cell:

  • $a = 5.067\text{ \AA}, \quad b = 8.064\text{ \AA}, \quad c = 8.882\text{ \AA}$

  • $\alpha = 91.64^\circ, \quad \beta = 93.34^\circ, \quad \gamma = 94.01^\circ$

• Structure Solution Correction Factors: Extracting structure factors from thin film intensity rods requires multiplication by geometry-specific corrections:

    $I_{corr} = I_{obs} \times (\text{rod interception factor}) \times (\text{Lorentz factor}) \times (\text{polarization factor}) \times (\text{area factor})$

• Monolayers: Probes 2D lattices in single molecular sheets (e.g., quinquethiophene monolayers on $\text{SiO}_x$).

Rotating GIXD & Alignment Error Analysis

• Sample Rotation ($\phi$-scan): Rotating the sample $360^ ext{o}$ around its surface normal maps full 3D reciprocal volumes and computes pole figures for epitaxial texture determination.

• Alignment Error Signatures ($V. Holzer, 2022$):

  • Sample $z$\text{-height error}: Shifts diffraction peaks perpendicular to Debye-Scherrer rings.

  • Incident angle $\alpha_i$ or tilt error: Shifts diffraction peaks along Debye-Scherrer rings.

Instrumentation Comparison

• Laboratory GIXD (e.g., Bruker D8 DISCOVER):

  • Multilayer mirror, long Soller slits, 6-axis Eulerian cradle ($\phi, \chi, x, y, z$), VANTEC 1D detector.

  • Measurement duration per map: Several days.

• Synchrotron GIXD (e.g., Beamline I07 at Diamond Light Source, UK):

  • Large 2D area detectors (PILATUS 2M).

  • Measurement duration per map: $1 - 30\text{ seconds}$.

  • Reference: O. Werzer, S. Kowarik, F. Gasser, Z. Jiang, J. Strzalka, C. Nicklin, R. Resel, Nat. Rev. Meth. Primers ($2024$).

Inelastic Scattering

Physical Characteristics of Phonons and Neutrons

Thermal Neutrons

• Mass: $m_n = 1.675 \times 10^{-27}\text{ kg}$

• Charge: $0$ (interacts directly with atomic nuclei via short-range nuclear forces).

• De Broglie Wavelength & Momentum: $p = m_n v = \hbar k$

  • For $\lambda = 1.65\text{ \AA}$: Wavevector $k = \frac{2\pi}{\lambda} = 3.80\text{ \AA}^{-1}$.

• Energy: $E = \frac{1}{2} m_n v^2 = \frac{\hbar^2 k^2}{2 m_n}$

  • For $\lambda = 1.65\text{ \AA}$: Energy $E = 25\text{ meV}$, which matches thermal energy at room temperature ($T \approx 300\text{ K}$).

Phonons

• Collective lattice vibrations propagating through crystalline solids.

• Momentum Range: $0 < k < \pi/a \sim 1\text{ \AA}^{-1}$ (Brillouin zone boundary).

• Energy Range: $0 - 100\text{ meV}$.

Energy Exchange Regimes

• Elastic Scattering: No energy exchange ($\Delta E = 0, |\vec{k}_f| = |\vec{k}_i|$).

• Quasi-Elastic Scattering: Continuous energy exchange centered near $\Delta E = 0$, caused by stochastic relaxations, diffusion, or unquantized motions.

• Inelastic Scattering: Quantized energy exchange involving creation or destruction of elementary excitations:

  • Inelastic Loss / Phonon Creation: $\hbar\omega < 0$ ($|\vec{k}_f| < |\vec{k}_i|$). Incoming particle transfers energy to the crystal lattice.

  • Inelastic Gain / Phonon Annihilation: $\hbar\omega > 0$ ($|\vec{k}_f| > |\vec{k}_i|$). Incoming particle absorbs energy from a thermal lattice vibration.

Conservation Laws and Reciprocal Space Kinematics

• Momentum Conservation:

    $\vec{q} = \vec{k}_f - \vec{k}_i = \vec{G}_{hkl} \pm \vec{k}_{phonon}$

    where $\vec{G}_{hkl}$ is a reciprocal lattice vector.

• Energy Conservation:

    $E_f = E_i \pm \hbar\omega_{phonon}$

• Translational Symmetry: Phonon modes can be probed in higher Brillouin zones where $\vec{q} = \vec{G}_{hkl} + \vec{k}_{phonon}$.

Phonon Dispersion Relations & Measurement Modes

• Dispersion Curve: Maps excitation frequency $\nu$ or energy $\hbar\omega$ against reduced wavevector $\vec{k}_{phonon}$ along principal symmetry directions (e.g., $[h00]$, $[\zeta\zeta 0]$, $[\zeta\zeta\zeta]$).

• Measurement Modes:

  • Constant-$q$ Scan: Momentum transfer $\vec{q}$ is fixed while energy transfer $\Delta E = E_f - E_i$ is swept.

  • Constant-Energy Scan: Energy transfer $\Delta E$ is held constant while momentum transfer $\vec{q}$ is stepped across the Brillouin zone.

• Branch Classifications: Acoustic modes (LA: Longitudinal Acoustic, TA: Transverse Acoustic) and Optical modes (LO: Longitudinal Optical, TO: Transverse Optical).

Instrumentation: Triple-Axis Spectrometer (TAS)

• Working Principle: Operates using three independent axes of rotation:

  1. Axis 1 (Monochromator): Selects incident wavevector $\vec{k}_i$ from the source beam.

  2. Axis 2 (Sample Stage): Orients the crystal to set the spatial direction of $\vec{k}_i$ relative to the lattice.

  3. Axis 3 (Analyzer): Analyzes magnitude and scattering angle of exit wavevector $\vec{k}_f$

• Allows selective access to any point in 4D $\left(\vec{q}, \hbar\omega\right)$ energy-momentum space.

Comparison: Inelastic Neutrons vs. Inelastic X-Rays

Epitaxial Growth Regimes

• Strained Layer (Pseudomorphic Growth): The film lattice distorts to match the substrate lattice constant ($a<em>s$), maintaining coherent interface registry. This regime occurs when lattice mismatch is small: $\frac{|a - a</em>s|}{a} < 0.1\%$

• Relaxed Layer: Above a critical film thickness, misfit dislocations form at the interface, relieving strain and allowing the film to revert to its native bulk lattice constant ($a$).

Reciprocal Space Broadening Mechanisms

• Finite Crystal Size: Broadens Bragg reflection nodes along both lateral ($q<em>{xy}$) and vertical ($q</em>z$) reciprocal space directions.

• Mosaicity: Causes arc-like angular broadening perpendicular to the scattering vector ($\mathbf{q}$).

• Lattice Parameter Variations ($\Delta d/d$): Shifts and broadens reflection nodes along the radial scattering direction.

Scanning Geometries in Reciprocal Space

• Coplanar Geometry: The primary beam, surface normal, and diffracted beam all reside within a single scattering plane.

• Scan Modes:

  • Specular Scan ($\theta / 2\theta$ or $\omega / 2\theta$): Incident and exit angles step synchronously ($\omega = \theta$). Probes structural and electron density variations strictly perpendicular to the surface ($q_z$).

  • Rocking Curve ($\omega$-scan): The detector angle $2\theta$ is held fixed while the sample angle $\omega$ rotates. Measures mosaicity and lateral orientation distribution ($q_x$).

  • Detector Scan ($2\theta$-scan): The sample angle $\omega$ remains stationary while the detector angle $2\theta$ scans.

Crystal Truncation Rods (CTRs)

• Physical Origin: The abrupt termination of a periodic crystal lattice at a flat surface breaks three-dimensional translational symmetry.

• Mathematical Formalism: Formulated as the Fourier transform of a three-dimensional semi-infinite crystal lattice multiplied by a step function $\Theta(z)$.

• Reciprocal Space Signature: Continuous lines of scattered intensity (truncation rods) extend along $q_z$ perpendicular to the surface, intersecting bulk Bragg points.

• Asymptotic Intensity Decay: Scattered intensity decays proportional to $q_z^{-2}$ away from Bragg peaks.

• Structural Sensitivity: CTR profiles are extremely sensitive to sub-angstrom surface features, including interface roughness, surface relaxation, and surface reconstruction (e.g., $\text{Si}(100)$ reconstruction).

X-Ray Reflectivity (XRR)

Refractive Index of Materials for X-Rays

• Complex Index of Refraction:
  $n = 1 - \delta + i\beta$

• Refractive Index Decrement ($\delta$): Describes phase velocity modification (typically $\delta \sim 10^{-5}$ to $10^{-6}$): $\delta = \frac{\lambda^2}{2\pi} r<em>e \rho</em>e$ where:

  • $\lambda$ is the X-ray wavelength.

  • $r<em>e = \frac{e^2}{4\pi \varepsilon</em>0 m c^2} = 2.818 \times 10^{-15}\,\text{m}$ is the classical electron radius.

  • $\rho_e$ is the electron density of the material.

• Absorption Index ($\beta$): Describes photoelectric absorption (typically $\beta \sim 10^{-7}$ to $10^{-8}$):
  $\beta = \frac{\lambda}{4\pi} \mu<em>x$ where $\mu</em>x$ is the linear absorption coefficient.

Total External Reflection

• Because $n < 1$ in condensed matter for hard X-rays, total external reflection occurs when X-rays impinge from vacuum or air ($n<em>1 = 1$) at grazing incidence angles below a critical angle $\alpha</em>c$.