X-Ray and Neutron Scattering 2

Fundamentals of Elastic Neutron-Nucleus Interaction

• Elastic neutron scattering from an atomic nucleus involves an incoming plane wave interacting with a fixed nuclear target, generating an outgoing spherical wave:   $\psi_i(x) = \psi_0 e^{i k_0 x}$   $\psi_f(\mathbf{r}) = \psi_0 \frac{f(\lambda, 2\theta)}{r} e^{i k r}$   

  

• Quantum mechanical description of the interaction is governed by the time-independent Schrödinger equation:   $\nabla^2 \psi(\mathbf{r}) + \frac{2m}{\hbar^2} \left( E - V(\mathbf{r}) \right) \psi(\mathbf{r}) = 0$

• The wave solution combining the incident plane wave and outgoing scattered spherical wave takes the form:   $\psi(\mathbf{r}) = \psi_0 e^{i k_0 x} + \psi_0 \frac{f(\lambda, 2\theta)}{r} e^{i k r}$

• Low-energy asymptotic expansion of the scattering amplitude $f(\lambda, 2\theta)$ yields:   $f(\lambda, 2\theta) = -b + i k b^2 + O(k^2)$   where $b$ is the Fermi scattering length, describing the inherent ability of a specific nucleus to scatter a neutron.

• The effective interaction potential $V(\mathbf{r})$ is represented by the Fermi pseudo-potential:   $V(\mathbf{r}) = \frac{2\pi\hbar^2}{m} b \,\delta(\mathbf{r})$

• Physical characteristics of the Fermi scattering length $b$:

  • Spatial magnitude: Typical value is on the order of $10^{-15}\,\text{m}$ ($1\,\text{fm}$).

  • Isotropic angular distribution: Deflection is uniform in all spatial directions due to nuclear dimensions being orders of magnitude smaller than thermal neutron wavelengths.

  • Energy independence: For thermal neutrons, $b$ is independent of the neutron wavelength.

  • Point-source behavior: Strong and highly localized interaction behaves as scattering from a spatial point source.

  • Complex scattering length representation: $b = b' + i b''$, where the imaginary component $b''$ accounts for absorption processes via $(n,\gamma)$ nuclear resonances (notably pronounced in samarium Sm, gadolinium Gd, boron B, and cadmium Cd).

  • Phase shift sign conventions: A positive sign ($+$, shift $0^\circ$) indicates in-phase scattering, whereas a negative sign ($-$, shift $180^\circ$) indicates out-of-phase scattering.

Nuclear Spin Effects and Incoherent Neutron Scattering

• The total scattering length $b$ depends on both the nuclear composition and the coupling between the nuclear spin $I$ and the neutron spin $s$ ($s = 1/2$):   $b = b_0 + \frac{1}{2} b_N \mathbf{I} \cdot \mathbf{s}$   where $b_0$ represents the nuclear constituent contribution, and $b_N$ represents the spin-dependent contribution.

• For any nucleus with a non-zero spin ($I \neq 0$), the interaction yields two distinct scattering lengths corresponding to the parallel ($b^+$) and antiparallel ($b^-$) coupled spin states:   $b^+ = b_0 + \frac{1}{2} b_N I$   $b^- = b_0 - \frac{1}{2} b_N (I + 1)$

• Incoherent neutron scattering arises fundamentally from the random spatial distribution of these nuclear spin orientation states across an ensemble of identical nuclei.

• Detailed Case Study: Hydrogen ($^1\text{H}$):

  • Hydrogen nuclear spin: $I = 1/2$; neutron spin: $s = \pm 1/2$. Total spin states $S = I + s$

  • Triplet state ($S = 1$): $3$ degenerate states, probability $w^+ = 3/4$, scattering length $b^+ = 1.085 \times 10^{-14}\,\text{m}$.

  • Singlet state ($S = 0$): $1$ state, probability $w^- = 1/4$, scattering length $b^- = -4.750 \times 10^{-14}\,\text{m}$.

  • Mean (coherent) scattering length $b_{\text{coh}} = \langle b \rangle$:     $b_{\text{coh}} = \langle b \rangle = \frac{3}{4} b^+ + \frac{1}{4} b^- = -0.374 \times 10^{-14}\,\text{m} = -3.74\,\text{fm}$

  • Mean square value $\langle b^2 \rangle$:     $\langle b^2 \rangle = \frac{3}{4} (b^+)^2 + \frac{1}{4} (b^-)^2$

  • Root-mean-square variation (incoherent scattering length $b_{\text{inc}}$):     $b_{\text{inc}} = \sqrt{\langle b^2 \rangle - \langle b \rangle^2} = 2.527 \times 10^{-14}\,\text{m} = 25.27\,\text{fm}$

• Neutron Scattering Cross Sections:

  • Definition: Cross section $\sigma$ quantifies the ratio of total scattered rate to incident flux:     $\sigma = \frac{\text{number of scattered neutrons/time}}{\text{incident neutron flux}} = \left[ \frac{\text{time}^{-1}}{\text{time}^{-1} \text{area}^{-1}} \right] = \text{area}$

  • Standard Unit: $1\,\text{barn} = 10^{-28}\,\text{m}^2 \approx \text{nuclear radius}^2$

  • Total scattering cross section expansion into coherent and incoherent parts:     $\sigma = 4\pi \langle b^2 \rangle = 4\pi \langle b \rangle^2 + 4\pi b_{\text{inc}}^2 = \sigma_{\text{coh}} + \sigma_{\text{inc}}$

  • Differential cross section integral:     $\sigma = \int d\varphi \int |b|^2 \sin(2\theta) \, d(2\theta) = 4\pi |b|^2$

• Cross Section Values across Selected Nuclides:   

  

  • Coherent cross sections ($\sigma_{\text{coh}}$) provide structured Bragg diffraction data, whereas incoherent cross sections ($\sigma_{\text{inc}}$) generate a isotropic, featureless background signal.

  • Specific cross section metrics:

  • $^1\text{H}$: $\sigma_{\text{coh}} = 1.8\,\text{barn}$, $\sigma_{\text{inc}} = 80.2\,\text{barn}$ (dominated by incoherent spin scattering).

  • $^2\text{H}$ (Deuterium): $\sigma_{\text{coh}} = 5.6\,\text{barn}$, $\sigma_{\text{inc}} = 2.0\,\text{barn}$ (ideal for isotopic substitution in structural studies).

  • Carbon ($\text{C}$): $\sigma_{\text{coh}} = 5.6\,\text{barn}$, $\sigma_{\text{inc}} = 0.0\,\text{barn}$.

  • Oxygen ($\text{O}$): $\sigma_{\text{coh}} = 4.2\,\text{barn}$, $\sigma_{\text{inc}} = 0.0\,\text{barn}$.

  • Aluminum ($\text{Al}$): $\sigma_{\text{coh}} = 1.5\,\text{barn}$, $\sigma_{\text{inc}} = 0.0\,\text{barn}$ (used routinely for sample container windows due to low scattering).

  • Vanadium ($\text{V}$): $\sigma_{\text{coh}} = 0.02\,\text{barn}$, $\sigma_{\text{inc}} = 5.0\,\text{barn}$ (nearly pure incoherent scatterer, utilized for sample containers and intensity normalization).

  • Iron ($\text{Fe}$): $\sigma_{\text{coh}} = 11.5\,\text{barn}$, $\sigma_{\text{inc}} = 0.4\,\text{barn}$.

  • Cobalt ($\text{Co}$): $\sigma_{\text{coh}} = 1.0\,\text{barn}$, $\sigma_{\text{inc}} = 5.2\,\text{barn}$.

  • Copper ($\text{Cu}$): $\sigma_{\text{coh}} = 7.5\,\text{barn}$, $\sigma_{\text{inc}} = 0.5\,\text{barn}$.

  • Argon-36 ($^{36}\text{Ar}$): $\sigma_{\text{coh}} = 24.9\,\text{barn}$, $\sigma_{\text{inc}} = 0.0\,\text{barn}$.

  • Cadmium ($\text{Cd}$): Complex scattering length $b = 4.87 - 0.7i\,\text{fm}$.

Neutron Magnetic Scattering and Form Factors

• Neutrons possess an intrinsic spin magnetic moment $\boldsymbol{\nu}$ that interacts dipolar-wise with the magnetic field $\mathbf{B}$ created by unpaired valence electrons in magnetic atoms.

• Fundamental Magnetic Constants:

  • Neutron magnetic moment: $\mu = -1.913 \mu_N$

  • Nuclear magneton: $\mu_N = \frac{e\hbar}{2m_p} = 5.051 \times 10^{-27}\,\text{J\,T}^{-1}$

• Magnetic field generated by electrons at position $\mathbf{R}$ with spin $\boldsymbol{\nu}_e$ and velocity $\mathbf{v}_e$:   $\mathbf{B} = \frac{\mu_0}{4\pi} \left[ \nabla \times \left( \frac{\boldsymbol{\nu}_e \times \mathbf{R}}{R^3} \right) - \frac{e \mathbf{v}_e \times \mathbf{R}}{R^3} \right]$

• Magnetic interaction potential $V_m(\mathbf{r})$:   $V_m(\mathbf{r}) = \boldsymbol{\nu} \cdot \mathbf{B}$

• Neutron Magnetic Form Factor $f(q)$:

  • Defined as the Fourier transform of the spatial magnetization distribution of an isolated magnetic atom.   

    

  • Decay characteristics: The magnetic form factor drops off much more rapidly with increasing scattering vector $q$ than the X-ray atomic form factor.

  • Physical Origin: Unpaired magnetic electrons are situated strictly in outer valence orbitals, yielding a broader spatial distribution in real space, which transforms to a narrower function in reciprocal $q$-space.

• Total Structure Factor $F(\mathbf{q})$ for Neutron Scattering on Atomic Ensembles:   $F(\mathbf{q}) = \sum_j b_j e^{i \mathbf{q} \cdot \mathbf{r}_j} + \sum_j f_j(q) e^{i \mathbf{q} \cdot \mathbf{r}_j}$   where $b_j$ is the nuclear scattering length and $f_j(q)$ is the magnetic form factor of atom $j$.

Elastic X-ray Scattering by a Single Electron (Thomson Formula)

• Problem Decomposition Hierarchy:

  1. Elementary interaction: Scattering of electromagnetic wave by a single free electron (Thomson formula).

  2. Atomic summation: Summing scattering amplitudes of all bound electrons within an isolated atom (Atomic form factor).

  3. Crystal summation: Coherent summation over all atomic sites in a periodic 3D lattice (Structure factor and lattice interference).

• Kinematic Theory Assumptions:

  • Assumes a single scattering event per photon (first-order Born approximation of dynamical diffraction theory).

  • Ignores refractive index alterations, wave attenuation by diffraction, and multiple internal scattering inside the crystal.

  • Applicable to real, imperfect, mosaic, or sub-micrometer crystals ($< 1\,\mu\text{m}$).

• Classical Thomson Scattering Derivation:

  • Incident electric field drives a free electron into acceleration, acting as a classical oscillating dipole dipole antenna.   

    

  • Classical electron radius $r_e$:     $r_e = \frac{e^2}{4\pi\epsilon_0 m c^2} = 2.81794 \times 10^{-15}\,\text{m}$

  • Radiated electric field amplitude $E$ at distance $r$ and emission angle $\psi$ relative to acceleration axis:     $E = -E_0 r_e \frac{e^{i(\mathbf{k} \cdot \mathbf{r} - \omega t)}}{r} \sin\psi$

  • Scattered wave field amplitude along scattering angle $2\theta$:     $E_{\text{scat}} = -E_0 \frac{r_e}{r} \cos(2\theta) e^{i(\mathbf{k} \cdot \mathbf{r} - \omega t)}$     where $\mathbf{k} = \frac{2\pi}{\lambda} \mathbf{s}$ is the scattered wavevector, and $\mathbf{s}$ is the propagation unit vector.

• Polarization Factor $P$:

  • Intensity $I$ of scattered wave:     $I = |E_{\text{scat}}|^2 = I_0 \frac{r_e^2}{r^2} P$

  • Values of $P$ based on incident beam polarization state:

  • Electric field parallel to scattering plane: $P = \cos^2(2\theta)$

  • Electric field perpendicular to scattering plane: $P = 1$

  • Unpolarized incident radiation: $P = \frac{1}{2} \left( 1 + \cos^2(2\theta) \right)$

• Irradiance and Photon Flux Measurement:

  • Irradiance (or intensity) measured across detector surface area $\Delta A = r^2 \Delta \Omega$:     $\frac{\Delta W}{\Delta t} = I \Delta A = I_0 \frac{r_e^2}{r^2} P \Delta A$

  • Differential cross-section per solid angle $\Delta \Omega$:     $\frac{\Delta W}{\Delta t \,\Delta \Omega} = I_0 r_e^2 P$

Atomic Scattering Factor and Electron Density Distribution

• Geometry of Atomic Scattering and Phase Shift Derivation:   

  

  • Path length difference between two scattering centers separated by spatial vector $\mathbf{r}$:     $\text{path difference} = b - a = \mathbf{r} \cdot \mathbf{s} - \mathbf{r} \cdot \mathbf{s}_0 = \mathbf{r} \cdot (\mathbf{s} - \mathbf{s}_0)$

  • Phase difference $\varphi$:     $\varphi = \frac{2\pi}{\lambda} \mathbf{r} \cdot (\mathbf{s} - \mathbf{s}_0) = \mathbf{q} \cdot \mathbf{r}$

• Definition of Scattering Vector $\mathbf{q}$ (Momentum Transfer):$\mathbf{q} = \frac{2\pi}{\lambda} (\mathbf{s} - \mathbf{s}_0) = \mathbf{k} - \mathbf{k}_0$

  • Magnitude of scattering vector $q$:     $q = |\mathbf{q}| = \frac{4\pi}{\lambda} \sin\theta$

  • Physical momentum transferred to photon: $\hbar \mathbf{q}$.

• Interference of Waves Scattered by Multiple Electrons:

  • For two electrons with phase difference $\varphi$:     $E_{\text{total}} = E_{\text{electron}} + E_{\text{electron}} e^{i \varphi}$     $I = E_{\text{total}} E_{\text{total}}^* = 2 E_{\text{electron}}^2 \left( 1 + \cos\varphi \right)$

  • At $\varphi = 0^\circ$: Complete constructive interference ($I = 4 E_{\text{electron}}^2$).

  • At $\varphi = 90^\circ$: Intermediate interference ($I = 2 E_{\text{electron}}^2$).

  • At $\varphi = 180^\circ$: Complete destructive interference ($I = 0$).

  • General continuous superposition over many electrons located at positions $\mathbf{r}_j$:     $E_{\text{total}} = E_{\text{electron}} \sum_j e^{i \mathbf{q} \cdot \mathbf{r}_j}$

• Atomic Scattering Factor (Atomic Form Factor) $f(\mathbf{q})$:

  • Defined as the 3D Fourier transform of the spatial atomic electron density distribution $\rho_{\text{atom}}(\mathbf{r})$:     $f(\mathbf{q}) = \int \rho_{\text{atom}}(\mathbf{r}) e^{i \mathbf{q} \cdot \mathbf{r}} \, d^3r$

  • Under isotropic spherical charge symmetry assumptions:     $f(q) = \int \rho_{\text{atom}}(r) e^{i q r} \, d^3r$

  • Limiting values:

  • Forward scattering limit ($q = 0$): $f(0) = Z$ (total atomic number / electron count).

  • High-angle limit ($q \to \infty$): $f(\infty) = 0$ (complete destructive interference among diffuse electron cloud components).

Diffraction from 3D Crystal Ensembles and Slit Interference Function

• Periodic Crystal Representation:

  • Crystal lattice vectors defined by repeating unit cells: $\mathbf{R}_n = n_1 \mathbf{a}_1 + n_2 \mathbf{a}_2 + n_3 \mathbf{a}_3$

  • Total position of atom $\alpha$ in unit cell $n$: $\mathbf{R}_{n\alpha} = \mathbf{r}_\alpha + \mathbf{R}_n$

  • Crystal domain size dimensions: $N_1 \mathbf{a}_1$, $N_2 \mathbf{a}_2$, $N_3 \mathbf{a}_3$, where $N_i$ represents the total number of unit cells along crystal axis $i$.

• Total Scattered Electric Field Amplitude:   $E_{\text{scat}} = E_{\text{electron}} \left( \sum_\alpha f_\alpha e^{i \mathbf{q} \cdot \mathbf{r}_\alpha} \right) \left( \sum_{n_1=0}^{N_1-1} e^{i \mathbf{q} \cdot n_1 \mathbf{a}_1} \right) \left( \sum_{n_2=0}^{N_2-1} e^{i \mathbf{q} \cdot n_2 \mathbf{a}_2} \right) \left( \sum_{n_3=0}^{N_3-1} e^{i \mathbf{q} \cdot n_3 \mathbf{a}_3} \right)$

• Slit Interference Function $L_N(\mathbf{q} \cdot \mathbf{a}_i)$:

  • Geometrical finite sum transformation:     $\left| \sum_{n=0}^{N-1} e^{i n \mathbf{q} \cdot \mathbf{a}} \right|^2 = \left| \frac{e^{i N \mathbf{q} \cdot \mathbf{a}} - 1}{e^{i \mathbf{q} \cdot \mathbf{a}} - 1} \right|^2 = \frac{\sin^2\left( \frac{N \mathbf{q} \cdot \mathbf{a}}{2} \right)}{\sin^2\left( \frac{\mathbf{q} \cdot \mathbf{a}}{2} \right)} = L_N(qa)$

  • Total Crystal Intensity Formula:     $I = I_{\text{electron}} |F(\mathbf{q})|^2 L_{N_1}(\mathbf{q} \cdot \mathbf{a}_1) L_{N_2}(\mathbf{q} \cdot \mathbf{a}_2) L_{N_3}(\mathbf{q} \cdot \mathbf{a}_3)$

• Fundamental Mathematical Properties of $L_N(qa)$:   

  

  • Principal Maxima Conditions: Occur when $\mathbf{q} \cdot \mathbf{a}_1 = 2\pi h$, $\mathbf{q} \cdot \mathbf{a}_2 = 2\pi k$, $\mathbf{q} \cdot \mathbf{a}_3 = 2\pi l$.

  • Maximum Peak Height: Scale quadratic with repeating unit count: $\max(L_N) = N^2$

  • Full Width at Half Maximum (FWHM) $\Gamma$: Scales inversely with unit cell count: $\Gamma \propto N^{-1}$

  • Integrated Area: Scales linearly with domain length: $\int L_N(qa) \, dq \propto N$

  • Infinite Crystal Limit ($N \to \infty$): $L_N(qa)$ converges rigorously to a periodic set of Dirac delta functions.

• Laue Fringes in Finite Thin Films:

  • Side maxima appearing between primary Bragg peaks are termed Laue fringes.

  • Observed experimentally in ultra-thin layered samples (such as a 5-layer cadmium stearate Langmuir-Blodgett film) where small finite layer count $N$ leaves visible inter-peak oscillations.

Structure Factor Theory, Centrosymmetry, and CsCl Example

• Form of Structure Factor $F_{hkl}$:   $F_{hkl} = \sum_\alpha f_\alpha e^{2\pi i (h x_\alpha + k y_\alpha + l z_\alpha)}$   where $(x_\alpha, y_\alpha, z_\alpha)$ are fractional atomic coordinates within the unit cell ($x_\alpha, y_\alpha, z_\alpha \in [0, 1]$).

  • Forward limit $F_{000}$ at $q = 0$: $F_{000} = \sum_\alpha Z_\alpha$ (total number of electrons per unit cell).

• Centrosymmetric Crystal Simplification:

  • In centrosymmetric space groups containing identical atomic pairs at $\mathbf{r}_\alpha$ and $-\mathbf{r}_\alpha$:     $F(\mathbf{q}) = \sum_\alpha f_\alpha \left( e^{i \mathbf{q} \cdot \mathbf{r}_\alpha} + e^{-i \mathbf{q} \cdot \mathbf{r}_\alpha} \right) = 2 \sum_\alpha f_\alpha \cos(\mathbf{q} \cdot \mathbf{r}_\alpha)$

  • Result: Structure factors $F_{hkl}$ in centrosymmetric crystals are strictly real numbers (phase angle $\phi$ is constrained to $0^\circ$ or $180^\circ$, meaning signs are purely positive or negative).

• Worked Example: Cesium Chloride (CsCl):

  • Crystal Structure: Primitive cubic lattice with two-atom basis:

  • Cesium ($\text{Cs}$) at fractional coordinate $(0, 0, 0)$

  • Chlorine ($\text{Cl}$) at fractional coordinate $(\frac{1}{2}, \frac{1}{2}, \frac{1}{2})$

  • Structure Factor Expansion:     $F_{hkl} = f_{\text{Cs}} e^{2\pi i (0)} + f_{\text{Cl}} e^{2\pi i (\frac{h}{2} + \frac{k}{2} + \frac{l}{2})} = f_{\text{Cs}} + f_{\text{Cl}} e^{i \pi (h + k + l)}$

  • Diffraction Selection Rules:

  • When $h + k + l$ is an even integer ($e.g., 110, 200$): $e^{i \pi (\text{even})} = +1$       $F_{hkl} = f_{\text{Cs}} + f_{\text{Cl}}$

  • When $h + k + l$ is an odd integer ($e.g., 100, 111$): $e^{i \pi (\text{odd})} = -1$       $F_{hkl} = f_{\text{Cs}} - f_{\text{Cl}}$

Reciprocal Lattice Formalism in 2D and 3D

• Mathematical Purpose of Reciprocal Space:

  • Used to expand lattice-periodic real space functions (such as crystal electron density $\rho(\mathbf{r})$) as 3D Fourier series:     $\rho(\mathbf{r}) = \rho(\mathbf{r} + \mathbf{R}_{uvw}) = \sum_{hkl} \varrho_{hkl} e^{i \mathbf{G}_{hkl} \cdot \mathbf{r}}$     where $\mathbf{R}_{uvw} = u\mathbf{a} + v\mathbf{b} + w\mathbf{c}$ ($u, v, w \in \mathbb{Z}$).

  • Reciprocal lattice vector definition: $\mathbf{G}_{hkl} = h \mathbf{a}^* + k \mathbf{b}^* + l \mathbf{c}^*$ ($h, k, l \in \mathbb{Z}$).

• Two-Dimensional Reciprocal Lattice Real-Space Relations:   

  

  • Real unit cell parameters: $a, c$, unit cell area $A$, real angle $\beta$.

  • Reciprocal parameters: $\mathbf{a}^* \perp \mathbf{c}$, $a^* = \frac{2\pi}{A \sin\beta}$, $\mathbf{c}^* \perp \mathbf{a}$, $c^* = \frac{2\pi}{A \sin\beta}$, reciprocal angle \beta^* = 180^\circ - \beta$.\n - 2D reciprocal vector \mathbf{g}{hl} = h\mathbf{a}^* + l\mathbf{c}^ satisfies:\n    \mathbf{g}{hl} \perp (hl) \quad \text{and} \quad |\mathbf{g}{hl}| = \frac{2\pi}{d{hl}}\n\n- Three-Dimensional Reciprocal Basis Vector Definitions:\n  \mathbf{a}^ = 2\pi \frac{\mathbf{b} \times \mathbf{c}}{(\mathbf{a} \times \mathbf{b}) \cdot \mathbf{c}}\n  \mathbf{b}^* = 2\pi \frac{\mathbf{c} \times \mathbf{a}}{(\mathbf{a} \times \mathbf{b}) \cdot \mathbf{c}}\n  \mathbf{c}^* = 2\pi \frac{\mathbf{a} \times \mathbf{b}}{(\mathbf{a} \times \mathbf{b}) \cdot \mathbf{c}}\n  where V_{EZ} = (\mathbf{a} \times \mathbf{b}) \cdot \mathbf{c} is the real space unit cell volume.\n\n- Orthogonality and Scalar Product Metric Rules:\n  \mathbf{a} \cdot \mathbf{a}^* = \mathbf{b} \cdot \mathbf{b}^* = \mathbf{c} \cdot \mathbf{c}^* = 2\pi\n  \mathbf{a} \cdot \mathbf{b}^* = \mathbf{b} \cdot \mathbf{a}^* = \mathbf{c} \cdot \mathbf{a}^* = \mathbf{a} \cdot \mathbf{c}^* = \mathbf{b} \cdot \mathbf{c}^* = \mathbf{c} \cdot \mathbf{b}^* = 0\n\n- Fundamental Geometric Correlations:\n 1. Vector Direction: \mathbf{G}{hkl}$is strictly collinear with the netplane normal$\mathbf{n}{hkl}$($\mathbf{G}{hkl} \perp (hkl)).\n 2. Vector Length: Magnitude equals inverse interplanar spacing d{hkl}:\n     |\mathbf{G}{hkl}| = \frac{2\pi}{d{hkl}}\n\n\n# Equivalence of Laue, Bragg, and Vector Diffraction Conditions\n\n- 1. Laue Condition:\n - Statement: Constructive interference occurs if and only if the scattering vector \mathbf{q}$exactly equals a reciprocal lattice vector$\mathbf{G}{hkl} in both magnitude and direction:\n    \mathbf{q} = \mathbf{G}{hkl}\n\n- 2. Ewald Sphere Construction:\n  \n - Sphere radius k = |\mathbf{k}0| = |\mathbf{k}| = \frac{2\pi}{\lambda}.\n - Center placed at point -\mathbf{k}_0$relative to origin$(0,0,0) of reciprocal space.\n - A diffraction peak is formed whenever a reciprocal lattice point intersects the surface of the Ewald sphere.\n - Modifications to bring reflections into diffraction condition: Rotating crystal (orienting \mathbf{G}{hkl}$), or using polychromatic radiation (varying sphere radius between$k_{\text{min}}$and$k_{\text{max}}).\n\n- 3. Real-Space Bragg Condition Derivation from Laue Condition:\n - Equating vector magnitudes from \mathbf{q} = \mathbf{G}{hkl}:\n    |\mathbf{q}| = |\mathbf{G}{hkl}| \implies \frac{4\pi}{\lambda} \sin\theta = \frac{2\pi}{d_{hkl}}\n    \lambda = 2 d_{hkl} \sin\theta\n - Bragg Law with integer diffraction order n$(for coprime Miller indices$hkl):\n    n \lambda = 2 d_{hkl} \sin\theta\n - Directional Specular Constraint: \mathbf{q} \parallel \mathbf{n}{hkl}$, forcing incident and scattered beams to form specular reflection angles$\theta$with respect to the$(hkl) lattice plane series.\n\n- 4. Vector Wavevector Form (Laue Equation):\n - Expressing scattered vector as \mathbf{k} = \mathbf{k}_0 + \mathbf{G}{hkl}$under elastic condition ($|\mathbf{k}|^2 = |\mathbf{k}0|^2):\n    |\mathbf{k}_0 + \mathbf{G}{hkl}|^2 = |\mathbf{k}0|^2 + 2 \mathbf{k}_0 \cdot \mathbf{G}{hkl} + |\mathbf{G}{hkl}|^2 = |\mathbf{k}_0|^2\n    2 \mathbf{k}_0 \cdot \mathbf{G}{hkl} + |\mathbf{G}{hkl}|^2 = 0\n\n\n# Finite Crystal Size Effects and Peak Intensity Quantities\n\n- Reciprocal Lattice Convolution for Finite Crystals:\n - Real space bounded crystal density is modeled as an infinite lattice modulated by a finite crystal shape envelope function \Omega(\mathbf{r})$($\Omega = 1$inside crystal volume$V_C$,$\Omega = 0 outside).\n - According to Fourier transform convolution theorem:\n    \mathcal{F}{\text{real crystal}} = \mathcal{F}{\text{crystal lattice}} \otimes \mathcal{F}{\text{crystal shape}}\n - Reciprocal lattice points expand from zero-dimensional delta points into 3D diffuse intensity distributions governed by the Geometric Shape Factor G(\mathbf{q}):\n    G(\mathbf{q}) = \frac{1}{V{EZ}} \int \Omega(\mathbf{r}) e^{i \mathbf{q} \cdot \mathbf{r}} \, d^3r\n - Total amplitude distribution:\n    E_{\text{scat}}(\mathbf{q}) \propto E_{\text{electron}} F(\mathbf{q}) G(\mathbf{q} - \mathbf{G}{hkl})\n\n- Integrated Bragg Peak Intensity Derivation:\n - Integrating local intensity over 3D reciprocal volume element (\Delta q)^3$surrounding reciprocal point$\mathbf{G}{hkl}:\n    I(\mathbf{G}{hkl}) = C \int |F(\mathbf{q})|^2 |G(\mathbf{q})|^2 \, d^3q \approx C |F(\mathbf{G}{hkl})|^2 \int |G(\mathbf{q})|^2 \, d^3q\n - Evaluating shape factor volume integral invariant \int |G(\mathbf{q})|^2 \, d^3q = \frac{8\pi^3 V_C}{V_{EZ}^2} yields:\n    I(\mathbf{G}{hkl}) = C |F(\mathbf{G}{hkl})|^2 \frac{8\pi^3 V_C}{V_{EZ}^2}\n\n- Complete Experimental Bragg Peak Intensity Formula I_{hkl}:\n  I_{hkl} = I_0 r_e^2 \frac{\lambda^3 V_C}{V_{EZ}^2} \cdot T \cdot P \cdot L \cdot H_{hkl} \cdot f_T \cdot |F_{hkl}|^2\n - Metric definitions:\n - I_0: Primary beam incident intensity\n - r_e$: Classical electron radius ($2.81794 \times 10^{-15}\,\text{m})\n - V_C: Coherently illuminated sample volume\n - V_{EZ}: Unit cell volume\n - T: Transmission coefficient accounting for beam absorption attenuation\n - P: Polarization factor\n - L: Lorentz factor\n - H_{hkl}: Reflection multiplicity factor (for powder diffraction integration)\n - f_T: Debye-Waller factor (intensity reduction due to thermal atomic vibrations)\n - |F_{hkl}|^2: Square of structure factor amplitude\n\n- Lorentz Factor L(\theta) Derivation for Rotating Crystal Setup:\n - Quantifies relative angular sweep velocity of a reciprocal lattice point traversing the Ewald sphere.\n - Linear velocity v$of reciprocal point rotating at angular frequency$\omega$at distance$q = \frac{4\pi}{\lambda} \sin\theta:\n    v = \omega \left( \frac{\lambda}{2\pi} q \right) = 2 \omega \sin\theta\n - Normal velocity component v_n perpendicular to Ewald sphere surface:\n    v_n = v \cos\theta = 2 \omega \sin\theta \cos\theta = \omega \sin(2\theta)\n - Dimensionless Lorentz factor L(\theta):\n    L(\theta) = \frac{\omega}{v_n} = \frac{1}{2 \sin\theta \cos\theta} = \frac{1}{\sin(2\theta)}\n\n\n# The Crystallographic Phase Problem and Solution Strategies\n\n- The Phase Problem Definition:\n  \n - Real-space electron density is derived via inverse Fourier transformation of structure factors:\n    \rho(x,y,z) = \frac{1}{V_{EZ}} \sum_{hkl} F_{hkl} e^{-2\pi i (h x + k y + l z)}\n - Experimental diffraction measurements detect only peak intensities I_{hkl} \propto |F_{hkl}|^2$, capturing structure factor magnitudes $|F_{hkl}|$ while losing the phase angles $\phi_{hkl}$.

• Phase Solution Methodologies:

  1. Test Structures (Trial and Error / Packing Algorithms):

  • Applicable when molecular geometry is known; unit cell content is packed algorithmically, enlarged (~130%), shrunk to cell dimensions, and optimized via DFT energy minimization and Molecular Dynamics to match calculated $|F_C|$ with experimental $|F_O|$.

  • Applied in polycyclic aromatic hydrocarbons such as ternaphthalene.

  1. Patterson Synthesis Function $P(X,Y,Z)$ (Heavy Atom Method):      

     

  • Formulated by A. L. Patterson (1934), performing a Fourier transform using purely measured raw intensities $|F_{hkl}|^2$ without phase requirements:        $P(X,Y,Z) = \frac{1}{V_{EZ}} \sum_{hkl} |F_{hkl}|^2 e^{-2\pi i \left( \frac{h}{a}X + \frac{k}{b}Y + \frac{l}{c}Z \right)}$

  • Real-space autocorrelation expression:        $P(X,Y,Z) = \frac{V_{EZ}}{a b c} \int_0^a \int_0^b \int_0^c \rho(x,y,z) \,\rho(x+X, y+Y, z+Z) \, dx \, dy \, dz$

  • Properties: Peaks in $P(X,Y,Z)$ correspond to interatomic vector distances $(X, Y, Z)$ between atoms, with peak heights proportional to the product of atomic numbers ($Z_A Z_B$).

  • Highly effective for metal-organic frameworks (MOFs) to extract heavy metal positions (e.g., Cu nodes) and linker alignment vectors.

  1. Direct Methods (Sayre Equation / Karle & Hauptman):

  • Developed by Jerome Karle and Herbert Hauptman (Nobel Prize in Chemistry 1985).

  • Uses statistical relationships between phase triplets among strong reflections:        $F_{hkl} \approx \sum_{h'k'l'} F_{h'k'l'} \cdot F_{(h-h', k-k', l-l')}$        $\phi_{hkl} \approx \phi_{h'k'l'} + \phi_{(h-h', k-k', l-l')}$

  • Highly automated for unit cells containing up to ~1000 non-hydrogen atoms, utilizing data oversampling (e.g., 186 spatial parameters constrained by thousands of measured reflections).

Anomalous Dispersion and Absolute Structure / Chirality Determination

• Physical Mechanism of Anomalous Dispersion:

  • Occurs when incident X-ray photon energy approaches an atomic absorption edge ($e.g.,$ Selenium Se K-edge at $E = 12.66\,\text{keV}$, $\lambda = 0.9795\,\text{\AA}$).

  • Resonant coupling modifies the atomic scattering factor with energy-dependent dispersion corrections $f'(\lambda)$ and $f''(\lambda)$:     $f(\mathbf{q}, \lambda) = f_0(\mathbf{q}) + f'(\lambda) + i f''(\lambda)$     where $f_0(\mathbf{q})$ is standard form factor, $f'(\lambda)$ is the real dispersion correction (photoelectric effect correction), and $f''(\lambda)$ is the imaginary absorption term.

• Breakdown of Friedel's Law in Non-Centrosymmetric Structures:

  • Standard Friedel Law states equal intensity for inverse reflections: $I_{hkl} = I_{\bar{h}\bar{k}\bar{l}}$.

  • In non-centrosymmetric structures containing an anomalous scatterer, Friedel symmetry breaks down:     $|F_{hkl}|^2 \neq |F_{\bar{h}\bar{k}\bar{l}}|^2$     $F_{hkl} = f_1 e^{i \mathbf{q} \cdot \mathbf{r}_1} + \left( f_2 + f'_2 + i f''_2 \right) e^{i \mathbf{q} \cdot \mathbf{r}_2}$     $F_{\bar{h}\bar{k}\bar{l}} = f_1 e^{-i \mathbf{q} \cdot \mathbf{r}_1} + \left( f_2 + f'_2 + i f''_2 \right) e^{-i \mathbf{q} \cdot \mathbf{r}_2}$

• Bijvoet Method for Chirality Determination:

  • Developed by J. M. Bijvoet (1951, 1954) to determine absolute molecular configuration (distinguishing $S$ and $R$ enantiomers).

  • Evaluates intensity differences between Friedel pairs: $\Delta F_{\text{exp}} = |F_{hkl}| - |F_{\bar{h}\bar{k}\bar{l}}|$

  • Bijvoet parameter $B$ compares calculated model differences against experimental data:     $B = \frac{\Delta F_{\text{calc}}}{\Delta F_{\text{exp}}}$

  • $B > 0$: Absolute configuration model is correct.

  • $B < 0$: Model configuration is inverted ($x_i, y_i, z_i \to -x_i, -y_i, -z_i$).

• Absolute Structure Determination using Laboratory Sources:

  • High-precision single crystal diffractometers utilizing $\text{Mo}\text{K}\alpha$ radiation ($E = 20.0\,\text{keV}$) can resolve absolute structures of light-element molecular crystals ($\text{O}, \text{C}, \text{H}$) without heavy anomalous scatterers by accurately measuring thousands of weak Friedel pair differences ($e.g.,$ space group $P2_12_12_1$, 5300 measured reflections, 2300 Friedel pairs, $R_1 = 3.4\%$, $wR_2 = 8.1\%$).

Experimental Single Crystal Diffraction Techniques (Laue and Rotating Crystal)

• Primary Single Crystal Measurement Strategies:

  1. Polychromatic / Fixed Crystal: Fix sample orientation $\mathbf{n}_{hkl}$, vary wavelength $\lambda$ (Laue Method).

  2. Monochromatic / Rotating Crystal: Fix wavelength $\lambda$, rotate crystal to alter orientation angles (Rotating Crystal / Modern Diffractometer).

• Historical Note: DNA structure determination by James Watson, Francis Crick, and Maurice Wilkins (Nobel Prize in Physiology or Medicine 1962) relied upon fiber X-ray diffraction data recorded by Rosalind Franklin.

• Laue Technique Setup and Characteristics:   

  

  • Historical Milestone: First experimental demonstration of X-ray diffraction by Max von Laue (1912, Nobel Prize in Physics 1914) using $\text{CuSO}_4$ (1st experiment) and $\text{ZnS}$ (5th experiment).

  • Experimental Parameters: Fixed single crystal exposed to continuous white radiation ($\lambda_{\text{min}} < \lambda < \lambda_{\text{max}}$).

  • $\lambda_{\min} \approx 0.31\,\text{\AA}$ ($40\,\text{keV}$ Bremsstrahlung limit).

  • $\lambda_{\max} \approx 4\,\text{\AA}$ (air absorption boundary).

  • Ewald Construction: Ewald spheres form a continuous nested volume between limiting radii $k_{\min} = \frac{2\pi}{\lambda_{\max}}$ and $k_{\max} = \frac{2\pi}{\lambda_{\min}}$. Multiple reflections (e.g., 330 and 120) satisfy diffraction simultaneously.

  • Advantages & Disadvantages: Captures vast reciprocal space volumes rapidly in one frame; reflection indexing is highly complex due to overlapping harmonic wavelengths.

  • Modern Applications: Crystal orientation alignment and Laue microdiffraction (mapping local grain orientations with microfocused beams).

• Rotating Crystal Method:   

  

  • Setup: Monochromatic X-ray beam hits a single crystal rotating about a fixed crystallographic axis surrounded by cylindrical film.

  • Ewald View: Reciprocal lattice points rotate through a static Ewald sphere surface, generating layer lines of reflections.

Modern Single Crystal Structure Workflow, Refinement, and Databases

• Step-by-Step Crystallographic Structure Solution Workflow:

  1. Crystal Selection & Mounting: Isolate sub-millimeter single crystal, mount on glass fiber or plastic loop, and center precisely in X-ray beam.

  2. Indexation & Unit Cell Parameter Determination: Measure 3D diffraction vectors $\mathbf{G}_{hkl}$ to define reciprocal basis vectors and cell dimensions ($a, b, c, \alpha, \beta, \gamma$).

  • Expected reflection counts: Inorganics ($10^3$), Organics ($10^4 - 10^5$), Biological macromolecules/proteins ($10^5 - 10^6$).

  1. Intensity Data Collection: Record complete 3D sphere of reflections using 2D area detectors.

  2. Phase Problem Solution: Extract initial rough atomic positions using automated black-box software ($e.g.,$ Olex2) executing Direct Methods or Patterson algorithms.

  3. Structure Model Refinement: Least-squares parameter optimization comparing observed ($F_O$) and calculated ($F_C$) structure factors.

• Atomic Refinement Parameters:

  • Isotropic refinement: 4 parameters per atom ($x, y, z$ coordinates + isotropic thermal displacement parameter $U$).

  • Anisotropic refinement: 9 parameters per atom ($x, y, z$ coordinates + 6 components of the anisotropic thermal tensor $U_{ij}$).

  • Residual Agreement Factor ($R$ -factor):     $R = \frac{\sum ||F_O| - |F_C||}{\sum |F_O|}$

  • $R < 5\%$: Excellent quality structure.

  • $R = 5\% - 10\%$: Acceptable structure solution.

• Exemplary Single Crystal Systems:

  • Biphenyl ($\text{C}_{12}\text{H}_{10}$):     Monoclinic, a = 8.12\,\text{\AA}, b = 5.64\,\text{\AA}, c = 9.47\,\text{\AA}, \beta = 95.4^\circ, Z = 2$, \rho_{\text{calc}} = 1.179\,\text{g\,cm}^{-3}$,$\rho_{\text{meas}} = 1.170\,\text{g\,cm}^{-3}, 164 electrons per unit cell, solved from 170 reflections.\n - Terthiophene derivative (\text{C}{28}\text{H}{40}\text{S}_3$, Formula Weight$472.78\,\text{g\,mol}^{-1}):\n - Low-Temperature Phase (100\,\text{K}$): Orthorhombic, space group$Pca2_1$,$a = 6.3465\,\text{nm}, b = 0.55248\,\text{nm}, c = 2.95109\,\text{nm}$,$V = 10.3475\,\text{nm}^3, Z = 16$, $\rho_{\text{calc}} = 1.214\,\text{g\,cm}^{-3}$, Herringbone packing motif.

  • Room-Temperature Phase ($296\,\text{K}$): Monoclinic, space group $P2_1/c$, $a = 0.55671\,\text{nm}, b = 6.3113\,\text{nm}, c = 0.78847\,\text{nm}, \beta = 107.604^\circ$, V = 2.6406\,\text{nm}^3, Z = 4$, \rho_{\text{calc}} = 1.189\,\text{g\,cm}^{-3}, Parallel stacking motif.\n\n- Standardized File Formats and Structural Databases:\n - CIF (Crystallographic Information File): World-standard format containing unit cell metrics, space group symmetry operators, fractional atomic coordinates, thermal displacement parameters, experimental details, R\n-factors, and literature references.\n - Primary Databases:\n - Cambridge Structural Database (CSD): Over 1 million organic and metal-organic structures (~50,000 new entries/year).\n - Inorganic Crystal Structure Database (ICSD): ~185,000 inorganic structures.\n - CRYSTMET: ~172,000 metals, alloys, and intermetallic phases.\n - Protein Data Bank (PDB): ~130,000 biological macromolecules (~10,000 new entries/year).\n - Crystallography Open Database (COD): ~376,000 open-access structures across all material classes.\n\n- PDB Diffraction Resolution Quality Standards:\n - > 4.0\,\text{\AA}$($< 20^\circ\,2\theta): Individual atomic coordinates meaningless; defines overall molecular envelope only.\n - 4.0 - 3.0\,\text{\AA}$($25^\circ - 30^\circ\,2\theta): Main polypeptide chain conformation reasonable; side chain placement problematic.\n - 2.5 - 2.0\,\text{\AA}$($35^\circ - 45^\circ\,2\theta$): Side chains resolved with few errors; bound solvent ($\text{H}2\text{O}) molecules resolvable.\n - 1.8 - 1.5\,\text{\AA}$($50^\circ - 60^\circ\,2\theta): Highly accurate structural models.\n - < 1.5\,\text{\AA}$($> 60^\circ\,2\theta): True atomic resolution.\n - AlphaFold AI Machine Learning: Predicts 3D protein folding structures directly from primary amino acid sequences trained on historical PDB structural data.\n\n\n# Dynamical Scattering Theory vs. Kinematical Scattering Theory\n\n- Conceptual Comparison Matrix:\n - Kinematical Theory: First-order approximation. Assumes single scattering per photon, zero optical refraction, and neglects wave attenuation by diffraction. Overestimates strong reflection intensities (extinction effect). Peak width \Delta q \to 0$, peak height$\propto N^2$, integrated peak intensity$I \propto |F{hkl}|^2$.

  • Dynamical Theory: Exact physical treatment. Solves Maxwell's equations inside a polarizable periodic medium, treating total internal wavefields as superpositions of incident and scattered spherical waves. Incorporates optical refraction and multiple internal scattering. Yields finite peak heights, peak width $\Delta q \propto |F_{hkl}|$, and integrated peak intensity $I \propto |F_{hkl}|$.

• Key Experimental Phenomena Explained Exclusively by Dynamical Theory:

  1. Optical Refraction Effects:      

     

  • X-rays possess a refractive index $n$ slightly less than unity ($n < 1$):        $n = 1 - \delta - i \beta$

  • Refraction alters the real space Bragg diffraction angle:        $\lambda = 2 d_{hkl} \sqrt{n^2 - \cos^2\theta}$

  • Shifts the peak position to slightly larger scattering vectors $Q$ compared to kinematical predictions.

  1. Borrmann Effect (Anomalous Transmission):

  • Occurs in large, highly perfect crystals under exact Bragg conditions.

  • Incident and diffracted waves interfere to form a standing wave field inside the crystal lattice.

  • For specific wave modes, standing wave nodes coincide precisely with atomic planes, suppressing photoelectric absorption and causing anomalously high X-ray transmission through thick crystals.

Darwin and Laue Formulations of Dynamical Scattering

• Darwin Formulation (C. G. Darwin, 1914):

  • Crystal Model: Stack of infinite parallel netplanes separated by spacing $d_{hkl}$. Uses Fresnel reflection and transmission equations at each plane.   

    

  • Darwin Equations for Infinite Perfect Crystal:

  • Reflection Angular Darwin Width $2\delta$:       $2\delta = \frac{\lambda^2 r_e \rho_{\text{atomic}}}{\pi \sin(2\theta)} |F_{hkl}|$

  • Peak Angular Shift $\Delta \Theta$:       $\Delta \Theta = \frac{\lambda^2 r_e \rho_{\text{atomic}}}{\pi \sin(2\theta)} F_0$

  • Integrated Dynamical Intensity $I_{\text{dyn}}$:       $I_{\text{dyn}} = I_0 \frac{8 \lambda^2 r_e \rho_{\text{atomic}}}{3\pi \sin(2\theta)} |F_{hkl}| \left( \frac{1 + \cos^2(2\theta)}{2} \right)$       where $\rho_{\text{atomic}}$ is atomic volume density, $F_{hkl}$ is structure factor, and $F_0 = F(q=0)$.

• Laue Formulation (Maxwell Equations in Periodic Media):

  • Represents crystal as a continuous polarizable medium with periodic electric polarizability $\chi(\mathbf{r})$ matching electron density periodicity:     $\chi(\mathbf{r}) = \frac{\lambda^2 r_e}{\pi} \rho_e(\mathbf{r}) = \chi(\mathbf{r} + \mathbf{R})$

  • Electromagnetic Wave Equation in Polarizable Medium:     $\nabla^2 \mathbf{E}(\mathbf{r}) + k^2 \mathbf{E}(\mathbf{r}) = \nabla \left( \nabla \cdot \mathbf{E}(\mathbf{r}) \right) - k^2 \chi(\mathbf{r}) \mathbf{E}(\mathbf{r})$

• Specific Physical Cases Demanding Dynamical Theory:

  • Primary extinction corrections in strong Bragg reflections.

  • Large, cm-sized highly perfect single crystals ($e.g.,$ semiconductor-grade silicon).

  • X-ray Topography (mapping internal crystal defects, dislocations, and strain fields).

  • Internal field standing wave phenomena ($e.g.,$ X-ray fluorescence and photoelectron generation).

  • Multiple Bragg beam simultaneous reflections.

  • Grazing Incidence X-ray Diffraction (GIXD) near the total external reflection critical angle.

Comparison of X-ray, Electron, and Neutron Interactions

• Fundamental Radiation Probe Comparison:

  • X-rays:

  • Interaction mechanism: Weak electromagnetic interaction with total electron density $\rho_e(\mathbf{r})$.

  • Penetration depth: Medium (micrometer $\mu\text{m}$ range).

  • Theoretical framework: Kinematical theory valid for small/imperfect crystals ($< 1\,\mu\text{m}$); dynamical theory required for large perfect crystals.

  • Electrons:

  • Interaction mechanism: Extremely strong Coulomb electrostatic interaction with both atomic electrons and positively charged atomic nuclei.

  • Penetration depth: Small (nanometer $\text{nm}$ range), rendering it highly surface-sensitive.

  • Theoretical framework: High scattering probability causes multi-scattering; dynamical theory is strictly mandatory ($e.g.,$ Low-Energy Electron Diffraction LEED).

  • Neutrons:

  • Interaction mechanism: Very weak short-range nuclear force interaction with atomic nuclei plus spin-dipolar magnetic interaction with unpaired valence electrons.

  • Penetration depth: Large (centimeter $\text{cm}$ range), enabling non-destructive bulk probing and extreme environment penetration.

  • Isotopic sensitivity: Each isotope exhibits unique scattering lengths $b$.

  • Theoretical framework: Kinematical scattering approach is almost universally sufficient.

• Neutron Interferometer Application (H. Rauch et al., 1974):   

  

  • Setup: Uses a monolithic, highly perfect silicon single crystal carved with three parallel identical crystal plates (splitter S, mirror M, analyzer A).

  • Mechanism: Coherently splits incoming neutron wave packets into two spatially separated trajectories (I and II). Passing an aluminum phase shifter sheet of thickness $\Delta D$ through one path modulates relative quantum phase $\Delta \phi$.

  • Result: Recombination at analyzer plate yields complementary sinusoidal intensity oscillations in the forward transmitted beam ($O$) and deviated diffracted beam ($H$) as a function of $\Delta D$, demonstrating macroscopic quantum coherence of neutrons.

Fundamentals of Elastic Neutron-Nucleus Interaction

• Elastic neutron scattering from an atomic nucleus involves an incoming plane wave interacting with a fixed nuclear target, generating an outgoing spherical wave:
  $\psi<em>i(x) = \psi</em>0 e^{i k<em>0 x}$ $\psi</em>f(\mathbf{r}) = \psi_0 \frac{f(\lambda, 2\theta)}{r} e^{i k r}$

Nuclear Spin Effects and Incoherent Neutron Scattering

• The total scattering length $b$ depends on both the nuclear composition and the coupling between the nuclear spin $I$ and the neutron spin $s$ ($s = \frac{1}{2}$):
  $b = b<em>0 + \frac{1}{2} b</em>N \mathbf{I} \cdot \mathbf{s}$
  where $b<em>0$ represents the nuclear constituent contribution, and $b</em>N$ represents the spin-dependent contribution.

• For any nucleus with a non-zero spin ($I \neq 0$), the interaction yields two distinct scattering lengths corresponding to the parallel ($b^+$) and antiparallel ($b^-$) coupled spin states:
  $b^+ = b<em>0 + \frac{1}{2} b</em>N I$
  $b^- = b<em>0 - \frac{1}{2} b</em>N (I + 1)$

• Incoherent neutron scattering arises fundamentally from the random spatial distribution of these nuclear spin orientation states across an ensemble of identical nuclei.

• Detailed Case Study: Hydrogen ($^{1}\text{H}$):

  • Hydrogen nuclear spin: $I = \frac{1}{2}$; neutron spin: $s = \pm \frac{1}{2}$. Total spin states $S = I + s$

  • Triplet state ($S = 1$): $3$ degenerate states, probability $w^+ = \frac{3}{4}$, scattering length $b^+ = 1.085 \times 10^{-14}\,\text{m}$.

  • Singlet state ($S = 0$): $1$ state, probability $w^- = \frac{1}{4}$, scattering length $b^- = -4.750 \times 10^{-14}\,\text{m}$.

  • Mean (coherent) scattering length $b<em>{\text{coh}} = \langle b \rangle$: $b</em>{\text{coh}} = \langle b \rangle = \frac{3}{4} b^+ + \frac{1}{4} b^- = -0.374 \times 10^{-14}\,\text{m} = -3.74\,\text{fm}$

  • Mean square value $\langle b^2 \rangle$:
    $\langle b^2 \rangle = \frac{3}{4} (b^+)^2 + \frac{1}{4} (b^-)^2$

  • Root-mean-square variation (incoherent scattering length $b<em>{\text{inc}}$): $b</em>{\text{inc}} = \sqrt{\langle b^2 \rangle - \langle b \rangle^2} = 2.527 \times 10^{-14}\,\text{m} = 25.27\,\text{fm}$

• Neutron Scattering Cross Sections:

  • Definition: Cross section $\sigma$ quantifies the ratio of total scattered rate to incident flux:
    $\sigma = \frac{\text{number of scattered neutrons/time}}{\text{incident neutron flux}} = \text{area}$

  • Standard Unit: $1\,\text{barn} = 10^{-28}\,\text{m}^2 \approx \text{nuclear radius}^2$

  • Total scattering cross section expansion into coherent and incoherent parts:
    $\sigma = 4\pi \langle b^2 \rangle = 4\pi \langle b \rangle^2 + 4\pi b<em>{\text{inc}}^2 = \sigma</em>{\text{coh}} + \sigma_{\text{inc}}$

  • Differential cross section integral:
    $\sigma = \int d\varphi \int |b|^2 \sin(2\theta) \, d(2\theta) = 4\pi |b|^2$

Neutron Magnetic Scattering and Form Factors

• Neutrons possess an intrinsic spin magnetic moment $\boldsymbol{\nu}$ that interacts dipolar-wise with the magnetic field $\mathbf{B}$ created by unpaired valence electrons in magnetic atoms.

• Fundamental Magnetic Constants:

  • Neutron magnetic moment: $\mu = -1.913 \mu_N$

  • Nuclear magneton: $\mu<em>N = \frac{e\hbar}{2m</em>p} = 5.051 \times 10^{-27}\,\text{J\,T}^{-1}$

  • Magnetic field generated by electrons at position $\mathbf{R}$ with spin $\boldsymbol{\nu}<em>e$ and velocity $\mathbf{v}</em>e$:
    $\mathbf{B} = \frac{\mu<em>0}{4\pi} \left[ \nabla \times \left( \frac{\boldsymbol{\nu}</em>e \times \mathbf{R}}{R^3} \right) - \frac{e \mathbf{v}_e \times \mathbf{R}}{R^3} \right]$

  • Magnetic interaction potential $V<em>m(\mathbf{r})$: $V</em>m(\mathbf{r}) = \boldsymbol{\nu} \cdot \mathbf{B}$

• Neutron Magnetic Form Factor $f(q)$:

  • Defined as the Fourier transform of the spatial magnetization distribution of an isolated magnetic atom.

Elastic X-ray Scattering by a Single Electron (Thomson Formula)

• Problem Decomposition Hierarchy:

  1. Elementary interaction: Scattering of electromagnetic wave by a single free electron (Thomson formula).

  2. Atomic summation: Summing scattering amplitudes of all bound electrons within an isolated atom (Atomic form factor).

  3. Crystal summation: Coherent summation over all atomic sites in a periodic 3D lattice (Structure factor and lattice interference).

• Classical Thomson Scattering Derivation:

  • Incident electric field drives a free electron into acceleration, acting as a classical oscillating dipole antenna.

  • Classical electron radius $r<em>e$: $r</em>e = \frac{e^2}{4\pi\epsilon_0 m c^2} = 2.81794 \times 10^{-15}\,\text{m}$

  • Radiated electric field amplitude $E$ at distance $r$ and emission angle $\psi$ relative to acceleration axis:
    $E = -E<em>0 r</em>e \frac{e^{i(\mathbf{k} \cdot \mathbf{r} - \omega t)}}{r} \sin(\psi)$

  • Scattered wave field amplitude along scattering angle $2\theta$:
    $E<em>{\text{scat}} = -E</em>0 \frac{r_e}{r} \cos(2\theta) e^{i(\mathbf{k} \cdot \mathbf{r} - \omega t)}$

• Polarization Factor $P$:

  • Intensity $I$ of scattered wave:
    $I = |E<em>{\text{scat}}|^2 = I</em>0 \frac{r_e^2}{r^2} P$

  • Values of $P$ based on incident beam polarization state:

  • Electric field parallel to scattering plane: $P = \cos^2(2\theta)$

  • Electric field perpendicular to scattering plane: $P = 1$

  • Unpolarized incident radiation: $P = \frac{1}{2} \left( 1 + \cos^2(2\theta) \right)$