Theoretical Foundation: The Scherrer Equation in Powder XRD
In an idealized, infinite three-dimensional crystal lattice, constructive interference between coherently scattered X-ray photons satisfies Bragg's law with mathematical perfection. In such hypothetical systems, scattering at any angle differing infinitesimally from the exact Bragg condition undergoes complete destructive interference due to the infinite sum of identical phase-shifted lattice planes.
In real nanomaterials, thin films, metal-organic frameworks (MOFs), and catalyst nanoparticles, the finite physical boundary of the crystal domain limits the total number of parallel reflecting planes. As this boundary decreases, cancellation of scattered waves adjacent to the Bragg angle becomes incomplete, yielding finite diffraction line breadth.
- D: Mean volume-weighted crystallite dimension along the crystallographic normal to reflecting lattice planes (h, k, l).
- K: Dimensionless shape factor reflecting domain morphology (0.94 for spheres; 1.00 for regular cubes).
- λ: Characteristic wavelength of the incident primary monochromatic X-ray beam.
- βsample: Corrected structural line broadening at half-maximum intensity, strictly expressed in dimensional radians.
- θ: Exact Bragg diffraction angle in radians or degrees (corresponding to half the recorded detector angle 2θ).
Mathematical Derivation from Finite Lattice Sums
Consider a crystal slab consisting of m parallel lattice planes separated by interplanar spacing d. The total thickness normal to the planes is D = m · d. Let the exact Bragg condition be satisfied at θB where path difference 2d · sin(θB) = λ.
As the incident angle shifts away from θB to an angle θ1 = θB + Δθ, the phase difference between scattering from the top surface plane and the bottom plane of the stack reaches 2π radians, causing the total scattered intensity to plunge to zero:
2m · d · sin(θ₂) = (m - 1) λ
Subtracting these two boundary conditions yields:
Because the full width at half-maximum β is approximately equal to half the total peak envelope (β ≈ Δθ = θ1 - θ2), introducing the geometric factor K establishes Paul Scherrer's classic relation: D = (K · λ) / (β · cos(θ)).
Crucial Instrumental Broadening Deconvolution (βinst)
A pervasive laboratory reporting error is substituting raw experimental full width at half maximum directly into the Scherrer formula. Diffractometers contribute physical instrumental broadening (βinst) stemming from finite X-ray tube target focal spot size, Soller slit axial divergence, beam footprint displacement, and detector channel resolution.
To isolate true specimen broadening (βsample), an unbroadened, macrocrystalline, defect-free reference standard (such as NIST SRM 660c LaB6 or annealed silicon wafer) must be analyzed across identical optic configurations:
Worked Laboratory Walkthrough: Zinc Oxide (101) Reflection
Sample: Hydrothermally synthesized ZnO wurtzite nanoparticles recorded on a Cu-Ka diffractometer.
- Diffractometer Parameters: Reflection centroid 2θ = 36.25°, recorded raw FWHM = 0.280°. Instrumental resolution at this angular coordinate: βinst = 0.080°. Radiation: Cu-Ka1 (λ = 1.54056 Å = 0.154056 nm). Shape factor: K = 0.94.
- Instrumental Deconvolution (Gaussian Model):
βdeg = √((0.280°)² - (0.080°)²) = √(0.0784 - 0.0064) = √0.0720 = 0.2683° - Unit Conversion to Natural Radians:
βrad = 0.2683° × (π / 180°) = 0.0046835 radians - Bragg Angle Determination:
θ = 36.25° / 2 = 18.125° → cos(18.125°) = 0.95036 - Scherrer Equation Computation:
D = (0.94 × 0.154056 nm) / (0.0046835 rad × 0.95036) = 32.53 nm (325.3 Å) - Dislocation Density (δ):
δ = 1 / D² = 1 / (32.53 × 10^-9 m)² = 9.45 × 10^14 lines/m²
Williamson-Hall Analysis: Decoupling Size from Microstrain
In many synthesized materials (such as ball-milled metals, core-shell architectures, or doped metal-organic frameworks), observed peak broadening is a combined convolution of **size broadening** and **lattice microstrain broadening**:
When plotting βcos(θ) on the y-axis against 4sin(θ) on the x-axis for multiple reflections:
- The slope yields intrinsic microstrain (ε).
- The y-intercept yields true strain-free crystallite size (D).
Frequently Asked Questions (PXRD Crystallographic Sizing)
What is the physical difference between crystallite size, grain size, and particle size?
Crystallite size denotes coherent, defect-free single crystalline sub-domains. A single grain may consist of multiple crystallites separated by small-angle boundaries, whereas a physical particle represents an agglomerate of multiple crystalline grains.
Why does the Scherrer equation become unreliable above 100 nm?
As crystallites exceed 100 nm, physical diffraction peak broadening converges toward the diffractometer's internal instrumental resolution. When peak breadth is dominated by optics rather than crystal dimensions, deconvolution subtraction creates massive relative error.
Why must FWHM be converted from degrees to radians?
The Scherrer equation is derived from Fourier wave mechanics where angular broadening β is a natural radian measure (arc length / radius). Entering raw degrees 2θ will underestimate crystallite size by a factor of 57.3 (π / 180).
How does lattice microstrain distort Scherrer crystallite results?
Lattice strain broadens peaks proportional to tan(θ). Applying the basic Scherrer formula to strained materials artificially treats strain-induced broadening as size truncation, underestimating crystallite size. Williamson-Hall plotting is recommended.
How should I determine the instrumental broadening βinst value?
Acquire a diffraction scan of a certified line standard (such as NIST SRM 660c LaB6 or high-purity strain-free silicon powder) under identical optical conditions. Fit the peak closest to your sample reflection to determine βinst.
What shape factor (K) should be selected for MOF and zeolitic crystals?
For cubic or octahedral MOFs (e.g., HKUST-1, UiO-66), K values between 0.89 and 0.94 are standard. If crystals exhibit pronounced anisotropic growth, shape factors along individual reflection planes must be adjusted.