Introduction
Bragg’s law is the cornerstone of X‑ray diffraction (XRD), a technique that allows scientists to peer into the atomic arrangement of crystalline solids. Formulated in 1913 by William Henry Bragg and his son William Lawrence Bragg, the law provides a simple quantitative relationship between the wavelength of incident X‑rays, the spacing between lattice planes in a crystal, and the angle at which constructive interference occurs. By measuring the angles at which intense diffracted beams appear, researchers can calculate the interplanar distances (d‑spacings) and, ultimately, deduce the crystal structure, phase composition, and even strain or texture of a material. In this article we will unpack Bragg’s law from its basic definition to its practical implications, illustrate the concept with step‑by‑step reasoning, give real‑world examples, discuss the underlying physics, clarify common misunderstandings, and answer frequently asked questions Simple as that..
Detailed Explanation
At its heart, Bragg’s law describes the condition for constructive interference of X‑rays scattered by the periodic array of atoms in a crystal lattice. Now, when an X‑ray beam strikes a set of parallel lattice planes, each plane acts as a source of secondary spherical waves. If the path difference between waves reflected from successive planes equals an integer multiple of the X‑ray wavelength, the waves reinforce each other, producing a bright diffracted spot. If the path difference is not an integer multiple, the waves interfere destructively and the intensity drops to zero.
Mathematically, the law is expressed as
[ n\lambda = 2d\sin\theta ]
where
- (n) is the order of reflection (a positive integer, usually taken as 1 for the first‑order peak),
- (\lambda) is the wavelength of the incident X‑ray radiation,
- (d) is the perpendicular distance between adjacent lattice planes (the interplanar spacing), and
- (\theta) is the angle between the incident ray and the crystal plane (often called the Bragg angle; the detector is placed at (2\theta) relative to the incident beam).
The equation tells us that for a given wavelength and crystal spacing, only certain angles satisfy the condition for constructive interference. Still, conversely, if we know (\lambda) and measure (\theta), we can solve for (d). This inverse relationship is what makes XRD a powerful structural probe Easy to understand, harder to ignore..
Step‑by‑Step or Concept Breakdown
1. Visualizing the Geometry
Imagine a stack of equally spaced atomic planes (like the pages of a book). An incoming X‑ray wavefront hits the top plane at an angle (\theta). Part of the wave is reflected specularly (angle of incidence equals angle of reflection) from that plane, while another portion penetrates to the next plane, reflects there, and emerges later.
2. Calculating the Path Difference
The extra distance traveled by the wave that reflects from the second plane compared to the first is the sum of two segments: the distance down to the second plane and the distance back up. Using simple trigonometry, each segment equals (d\sin\theta). Hence the total path difference is
[ \Delta = 2d\sin\theta . ]
3. Applying the Interference Condition
Constructive interference occurs when (\Delta) equals an integer number of wavelengths:
[ \Delta = n\lambda \quad (n = 0,1,2,\dots). ]
Substituting (\Delta) gives Bragg’s law:
[ n\lambda = 2d\sin\theta . ]
4. Solving for Unknowns
- If (\lambda) and (\theta) are known (typical in a laboratory diffractometer where the X‑ray source is fixed), compute (d = \frac{n\lambda}{2\sin\theta}).
- If (d) and (\lambda) are known (e.g., from a reference crystal), predict the diffraction angles: (\theta = \arcsin!\left(\frac{n\lambda}{2d}\right)).
5. Interpreting the Diffraction Pattern
A powder diffractometer records intensity versus (2\theta). Each peak corresponds to a specific set of lattice planes ((hkl)) that satisfy Bragg’s law. The position ((2\theta)) yields (d_{hkl}); the intensity reflects the arrangement of atoms within those planes (structure factor). By indexing all observable peaks, one can reconstruct the unit cell parameters and, with additional analysis, the full crystal structure.
Real Examples
Example 1: Determining the Lattice Parameter of Silicon
Silicon crystallizes in the diamond cubic structure with a known lattice constant (a \approx 5.On the flip side, 43) Å. Using Cu Kα radiation ((\lambda = 1.5406) Å), the most intense peak appears from the (111) planes That's the part that actually makes a difference. Which is the point..
- Calculate the interplanar spacing for (111):
[ d_{111} = \frac{a}{\sqrt{h^{2}+k^{2}+l^{2}}} = \frac{5.On top of that, 43}{\sqrt{3}} \approx 3. 136\text{ Å}.
- Insert into Bragg’s law for first order ((n=1)):
[ \sin\theta = \frac{\lambda}{2d} = \frac{1.In real terms, 5406}{2\times 3. In real terms, 136} \approx 0. 2455.
- Hence (\theta \approx 14.2^{\circ}) and the detector angle (2\theta \approx 28.4^{\circ}).
Indeed, the Si (111) peak is observed near 28.4° 2θ in a standard XRD pattern, confirming the lattice constant.
Example 2: Identifying an Unknown Phase in a Mixture
A forensic laboratory receives a white powder suspected to be either calcium carbonate (calcite) or calcium sulfate (gypsum). Also, 4^{\circ}) and (36. Here's the thing — using Mo Kα radiation ((\lambda = 0. 7107) Å), they record a diffraction pattern. On the flip side, two strong peaks appear at (2\theta = 29. 0^{\circ}).
- For the first peak, assuming (n=1):
[ d = \frac{\lambda}{2\sin\theta} = \frac{0.7^{\circ})} \approx 2.In practice, 7107}{2\sin(14. 80\text{ Å}.
-
This spacing matches the (104) plane of calcite ((d_{104}=2.80) Å) Practical, not theoretical..
-
The second peak yields (d \approx 2.48) Å, corresponding to the (110) plane of calcite as well.
Since the observed spacings align with calcite and not with gypsum (whose strongest peaks appear near 2θ ≈ 23° and 31° for Mo Kα), the sample is identified as calcium carbonate That alone is useful..
Example 3: Measuring Thin‑Film Strain
A epitaxial GaN film grown on sapphire is expected to be tensile‑strained because its lattice constant (a = 3.189 Å) slightly exceeds that of sapphire’s basal plane (a = 4.758 Å, but the matching direction
The strain state of an epitaxial film can be extracted directly from the shift of its Bragg peaks relative to those of a relaxed bulk reference. For a wurtzite GaN layer on sapphire, the in‑plane lattice mismatch is accommodated by a slight elongation (tensile strain) of the GaN basal plane and a concomitant contraction of the out‑of‑plane c‑axis to maintain volume conservation.
Step 1 – Choose a suitable reflection
The (0002) planes of GaN are parallel to the substrate surface and give a strong, narrow peak in θ‑2θ scans. Their spacing is related to the c lattice parameter by
[ d_{0002}= \frac{c}{2}. ]
Step 2 – Determine the expected c for a fully relaxed GaN crystal
From bulk GaN, (c_{0}=5.185) Å, giving
[ d_{0002}^{\text{relaxed}} = \frac{5.185}{2}=2.5925;\text{Å}. ]
Step 3 – Relate in‑plane strain to out‑of‑plane strain
For a hexagonal crystal under biaxial strain in the ab plane, the elastic compliance relation yields
[ \frac{\Delta c}{c_{0}} = -2\frac{C_{13}}{C_{33}},\frac{\Delta a}{a_{0}}, ]
where (C_{13}=106) GPa and (C_{33}=398) GPa for GaN. Think about it: the ratio (-2C_{13}/C_{33}\approx -0. 53).
Step 4 – Calculate the in‑plane strain
The sapphire basal plane lattice constant is (a_{\text{sap}}=4.758) Å. The GaN lattice constant that matches the sapphire along the ([1\bar{1}00]) direction is (a_{\text{match}} = a_{\text{sap}}/\sqrt{3}=2.748) Å (because the GaN hexagonal cell is rotated 30° relative to sapphire). The actual GaN a in the film is forced to this value, giving
[ \frac{\Delta a}{a_{0}} = \frac{a_{\text{match}}-a_{0}}{a_{0}} = \frac{2.748-3.189}{3.189}= -0.138. ]
Thus the film experiences a compressive in‑plane strain of –13.8 % (the negative sign indicates compression relative to relaxed GaN) That alone is useful..
Step 5 – Predict the out‑of‑plane lattice parameter
[ \frac{\Delta c}{c_{0}} = -0.53 \times (-0.073, ] [ c = c_{0},(1+0.Day to day, 073)=5. 185\times1.138) = +0.073\approx5.566;\text{Å}.
Consequently
[ d_{0002}^{\text{strained}} = \frac{c}{2}=2.783;\text{Å}. ]
Step 6 – Convert the shifted spacing to a diffraction angle
Using Cu Kα radiation ((\lambda=1.5406) Å) and first‑order Bragg law:
[
Using Cu Kα radiation (λ = 1.5406 Å) and first‑order Bragg law:
[ \sin\theta = \frac{\lambda}{2d}\qquad\text{and}\qquad 2\theta = 2\arcsin!\left(\frac{\lambda}{2d}\right). ]
Relaxed GaN reference
(d_{0002}^{\text{relaxed}} = 2.5925;\text{Å})
[ \sin\theta_{\text{rel}} = \frac{1.Also, 5406}{2\times2. 5925}=0.Consider this: 2972;;\Rightarrow;; \theta_{\text{rel}} = \arcsin(0. 2972)=17.29^{\circ}, ] [ 2\theta_{\text{rel}} = 34.58^{\circ}.
Strained film
(d_{0002}^{\text{strained}} = 2.783;\text{Å})
[ \sin\theta_{\text{str}} = \frac{1.5406}{2\times2.2768)=16.That said, 2768;;\Rightarrow;; \theta_{\text{str}} = \arcsin(0. Because of that, 09^{\circ}, ] [ 2\theta_{\text{str}} = 32. On the flip side, 783}=0. 18^{\circ}.
The tensile out‑of‑plane expansion therefore shifts the (0002) reflection to lower angles by
[ \Delta(2\theta)=2\theta_{\text{rel}}-2\theta_{\text{str}}\approx 2.40^{\circ}. ]
In practice one records a high‑resolution θ‑2θ scan, fits the GaN (0002) peak (and, if present, the sapphire substrate peak for alignment) with a pseudo‑Voigt or Pearson VII function, and extracts the peak position to a precision of better than 0.01° 2θ. The measured shift is then converted back to a strain value using the inverse of the steps above:
[ \frac{\Delta c}{c_{0}} = \frac{2d_{\text{meas}}-c_{0}}{c_{0}},\qquad \frac{\Delta a}{a_{0}} = -\frac{C_{33}}{2C_{13}}\frac{\Delta c}{c_{0}}. ]
For the present example, inserting the measured (2\theta_{\text{str}}=32.Here's the thing — a compressive in‑plane strain of –13. 138), i.e. 073) and consequently (\Delta a/a_{0}=-0.Think about it: 783;\text{Å}), (\Delta c/c_{0}=+0. 18^{\circ}) yields (d_{\text{meas}}=2.8 % and a tensile out‑of‑plane strain of +7.3 %, exactly as predicted from the elastic model.
Sources of uncertainty
- Instrumental zero‑shift (checked using a Si standard).
- Sample tilt or surface roughness, which can broaden the peak and bias the fitted position; a rocking‑curve measurement helps to correct for tilt.
- Stress relaxation at the film surface or interface intermixing, which would cause the measured strain to deviate from the ideal elastic prediction; complementary techniques such as reciprocal‑space mapping or Raman spectroscopy can be used to cross‑validate the XRD result.
By carefully accounting for these factors, the shift of the (00
- reflection can be reliably correlated to the biaxial strain state of the GaN epilayer. The methodology outlined here—starting from the elastic stiffness constants, calculating the Poisson-mediated lattice response, converting the resulting $c$-axis expansion into an angular shift via Bragg’s law, and finally inverting the measured peak position to obtain quantitative strain values—provides a self-consistent framework for non-destructive stress analysis in wurtzite heterostructures.
When combined with reciprocal-space mapping to separate tilt, curvature, and composition gradients, and corroborated by Raman spectroscopy for local stress verification, high-resolution XRD remains the gold standard for determining the elastic strain tensor in nitride-based devices. Accurate strain quantification is, in turn, essential for predicting piezoelectric polarization fields, band-structure modifications, and ultimately the optical and electronic performance of quantum wells, high-electron-mobility transistors, and vertical power devices grown on lattice-mismatched substrates.
