Manual Indexing and Lattice Parameter Calculation of Diffraction Data

XRD errors include: (1) Systematic errors - focus errors, sample holder errors, beam divergence, slit errors, sample absorption, zero-point shift, (2) Sample-related errors - preferred orientation, strain, crystallite size, amorphous content, impurities. Parallel beam methods avoid some errors (e.g., sample height errors) that affect focused beam methods. Proper sample preparation minimizes sample-related errors. Lattice parameter determination: measure diffraction angles precisely, especially near 90° for high precision. Methods: external standard method (using NIST SRM 640c, 640d, 640e), internal standard method (using known standard in sample), Nelson-Riley extrapolation, least-squares fitting, or full-pattern fitting. Temperature affects lattice parameters; thermal expansion must be considered. High-angle peaks are most important for precision. Crystallite size determination from XRD peak broadening has limitations: peak broadening depends on crystallite size, strain, and instrumental factors. Peak broadening becomes significant below 1μm; below 0.1μm, peaks may disappear. Strain and residual stress also cause broadening. Scherrer equation (D = Kλ/(βcosθ)) is commonly used but requires standard materials for calibration. Methods: Scherrer equation, Williamson-Hall method, fundamental parameters method (FP method). Crystal structure analysis workflow: (1) Single-crystal XRD preferred for initial structure determination, (2) Powder XRD requires careful sample preparation, (3) Indexing (determining lattice parameters) is critical and difficult, (4) Determine lattice parameters using Le Bail or Pawley methods, (5) Generate initial structural model, (6) Rietveld refinement for structure determination. Long measurement times (overnight) may be needed for high-quality data. Rietveld refinement: compare calculated pattern (from structural model) with measured pattern, refine structural parameters to minimize difference. Requires initial structural model from single-crystal data or literature. After refinement, verify physical reasonableness (no overlapping atoms, reasonable bond lengths). Compare with complementary techniques (NMR, neutron diffraction) for validation. Software: Rietan-FP, TOPAS, GSAS.

Accurate lattice parameter determination requires correcting systematic errors from sample positioning. Vertical and horizontal displacements cause angular shifts following ΔD/D = tan(θ)(Δz) and ΔD/D = tan(θ)(Δx), respectively, where Δz and Δx represent displacement magnitudes. These errors differ from actual lattice parameter changes, requiring correction methods. The Nelson-Riley method plots lattice parameter versus [sin(2θ)/sin(θ)] and performs linear regression, eliminating displacement errors with the zero intercept yielding true lattice parameters. This method becomes increasingly accurate at higher 2θ angles where displacement effects diminish. For grazing incidence measurements, refraction correction Δ(2θ) = 2α + 2β adjusts for beam deflection at the sample surface, with corrections varying significantly near the critical angle. These precision techniques enable reliable structural characterization essential for materials science research.

The method of least squares finds the best-fit straight line through experimental data by minimizing the sum of squared errors. For data points (x_i, y_i) related by y = a + bx, the normal equations are Σy = aN + bΣx and Σxy = aΣx + bΣx². Applying this to XRD data with eight diffraction lines, solving yields the best-fit equation a_M = 3.160 - 0.023 × cos²(theta). Plotting this line and extrapolating to cos²(theta) = 0 gives the precise lattice parameter a_C = 3.160 Å, demonstrating how statistical methods eliminate systematic errors in crystallographic measurements.

This lecture explains how to determine the lattice parameter of cubic crystals from X-ray diffraction data using Bragg's law (λ = 2D sinθ) and the relationship D = a/√(h²+k²+l²), where a is the lattice parameter. The instructor demonstrates that peaks at higher Bragg angles provide more accurate lattice parameter values due to reduced error propagation (Δa/a = -cotθ × Δθ). To obtain the most precise value, the lecture introduces extrapolation methods: the Bradley-J function (plotting a versus cos²θ) and the Nelson-Riley function (plotting a versus cos²θ/(sinθ + cosθ)), which extrapolate to θ = 90° to eliminate systematic errors. The instructor also briefly covers indexing non-cubic crystals (tetragonal system) and other X-ray diffraction applications including phase identification, phase quantification, and crystallite size determination.

Refined lattice parameters are reported with their standard errors. For example, a = 4.917 ± 0.002 Å. The error values indicate the precision of the refinement. Parameters are typically reported to 3-4 decimal places depending on the refinement quality.

Add phases via Import > Phase > from CIF file, adding one phase at a time unless certain all phases are present. Select GSAS recommended origin shift for symmetry-centered origins. In Phase Data window, enable Phase Fraction refinement with default value 1. Set Max Cycles to 10 iterations per run. Run refinement via Data Tree > Calculate > Refine. Monitor weighted residual Rw value—ideally approaching zero, values below 10% indicate good fit. Accept refined results when prompted, or reject if refinement deteriorates. For proper plot visualization, click on a phase first in Data Tree, then PWDR, then Home. The magnifying glass enables zooming into specific peaks. Plot elements include: blue plus marks for experimental data, red lines for background, green lines for calculated patterns, teal lines for differences, and vertical colored marks for expected peak positions.

Rietveld refinement is a whole-pattern fitting technique developed by Professor Hugo Rietveld in the 1960s that fits crystallographic parameters (unit cell dimensions, atomic coordinates, crystallite size, microstrain, site occupancies, and thermal displacement parameters) to experimental X-ray diffraction data by minimizing the difference between calculated and measured patterns using least squares optimization; the method requires prior phase identification and structural models (typically from CIF files) as starting points, and relies on accurate mathematical description of peak profiles either through empirical functions or the fundamental parameters approach based on instrument configuration.

This extended segment covers the complete workflow from phase acceptance to quantitative analysis. Users learn to manually select preferred candidates from the candidate list and reject unwanted automatic suggestions using the Escape key. The software provides semi-quantitative estimates of phase amounts before initiating Rietveld refinement via the FP button, which launches FullProf with default settings. During refinement, users monitor progress, R-factor values, and comments. Post-refinement, the software displays refined profiles and difference plots where minimal residuals confirm successful refinement. Final quantitative results replace initial estimates in the match list.

Rietveld refinement is a multi-parameter nonlinear least squares curve-fitting technique for analyzing powder X-ray diffraction patterns, where every point in the diffraction profile is modeled as the sum of contributions from all reflections, and the method simultaneously refines structural parameters (atom positions, thermal factors, site occupancies) along with peak shape parameters (instrumental broadening, crystallite size, microstrain) to minimize the weighted difference between observed and calculated intensities.

This tutorial introduces the Rietveld refinement method using GSAS software for materials engineering, demonstrating how to set up and begin the refinement process by loading sample data (such as LaB6 reference patterns) and configuring initial parameters for phase identification and structural analysis.

The Williamson Hall plot method determines crystallite size (D) and lattice strain (ε) from XRD data by plotting βcosθ against cosθ, where β is the full width at half maximum (FWHM) converted to radians, θ is the Bragg angle, K is a constant (typically 0.9), and λ is the X-ray wavelength (e.g., 0.15406 nm for Cu-Kα radiation). The crystallite size is calculated from the y-intercept using D = Kλ/intercept, while the lattice strain is obtained from the slope of the linear fit.

Beyond basic characterization, X-ray diffraction reveals critical microstructural information. Peak broadening arises from crystallite size effects (Scherrer equation) and microstrain (β_strain = 4εsinθ/cosθ). The Williamson-Hall plot separates these contributions, yielding crystallite size from the y-intercept and microstrain from the slope. Degree of crystallinity measures long-range order in semi-crystalline materials through deconvolution of crystalline and amorphous peaks. Pole figures visualize preferred orientation by mapping crystallographic plane distributions, essential for understanding texture in rolled metals and epitaxial films.

The Williamson-Hall method determines crystallite size and microstrain from XRD data through a systematic workflow. First, extract peak positions (2θ) and FWHM values from the XRD pattern. Convert 2θ to radians and calculate θ, cosθ, and sinθ for each peak. Compute β cosθ and 4 sinθ values. Plot β cosθ versus 4 sinθ and perform linear regression to obtain slope and intercept. Calculate crystallite size using D = Kλ/intercept and microstrain using ε = slope/4. This method separates size broadening from strain broadening, providing quantitative structural parameters from diffraction data.

XRD peak broadening reveals crystallite size (Scherrer equation: D = Kλ/(β cosθ)) and microstrain (Williamson-Hall method). Size broadening is isotropic; strain broadening is anisotropic. Detection limits: ~30 nm minimum, ~1000 nm maximum. Texture analysis reveals preferred orientation affecting mechanical properties. Residual stress analysis uses peak shifts from lattice parameter changes. These parameters correlate structure with material properties.

The video explains the Williamson-Hall (W-H) plot method for calculating crystallite size and microstrain from X-ray diffraction (XRD) data. It asserts that this method provides more accurate crystallite size measurements compared to the Scherrer formula. The W-H plot method accounts for both crystallite size broadening and lattice strain effects in XRD peak broadening, allowing simultaneous determination of these two parameters. The author demonstrates the authentic procedure for applying the W-H plot, emphasizing clarity in each step. The method involves plotting the product of the diffraction angle cosine and the full width at half maximum (FWHM) of XRD peaks against the sine of half the diffraction angle. The slope of the resulting linear fit gives the microstrain, while the intercept relates to the crystallite size. This approach corrects for the limitations of the Scherrer equation, which assumes only size-induced broadening and neglects strain. The video is part of a series on XRD data analysis, indicating its context within broader materials characterization workflows. The author encourages viewers to subscribe for more content and provides links to related videos on XRD peak indexing, d-spacing calculation, Rietveld refinement, and other nanomaterial analysis techniques.