2.5 Structure factor

Definition of structure factor

The structure factor describes the amplitude and phase of a wave diffracted from a crystal lattice. It connects the positions of atoms within a unit cell to the intensity of the diffracted beam you'd actually measure in an experiment. In other words, it's the bridge between "where atoms sit" and "what the diffraction pattern looks like."

The structure factor provides a quantitative measure of how atoms in a crystal collectively scatter incident radiation, whether that's X-rays, neutrons, or electrons.

Mathematical representation

The structure factor is a complex number (it has both amplitude and phase) and is written as $F_{hkl}$, where $h$, $k$, and $l$ are the Miller indices of the diffracting planes.

The core formula is:

$F_{hkl} = \sum_{j=1}^{N} f_j \exp[2\pi i(hx_j + ky_j + lz_j)]$

• $f_j$ is the atomic form factor of the $j$-th atom (its individual scattering strength)
• $x_j, y_j, z_j$ are the fractional coordinates of the $j$-th atom in the unit cell
• $N$ is the total number of atoms in the unit cell

Each atom contributes a scattered wave with its own amplitude ($f_j$) and phase (set by its position). The structure factor sums all of these contributions. If the waves add constructively, $|F_{hkl}|$ is large and you get a strong diffraction peak. If they cancel, $F_{hkl} = 0$ and that reflection is absent.

Physical interpretation

• The magnitude $|F_{hkl}|$ determines how intense the diffracted beam is.
• The phase of $F_{hkl}$ encodes information about where atoms sit relative to each other inside the unit cell.
• Together, amplitude and phase let you reconstruct the full electron density of the crystal, though getting the phase experimentally is notoriously difficult (more on that below).

Role in diffraction

The structure factor is what connects atomic-scale geometry to the diffraction patterns you collect in the lab. It determines both the intensity and the presence or absence of every diffraction peak.

Relationship to atomic positions

The structure factor is directly sensitive to where atoms are in the unit cell. Even small shifts in atomic coordinates change the phase terms $hx_j + ky_j + lz_j$, which changes $F_{hkl}$ and therefore the measured intensity. This sensitivity is exactly what makes diffraction so powerful for determining crystal structures.

Influence on diffracted intensity

The measured intensity of a diffracted beam is proportional to the squared modulus of the structure factor:

$I_{hkl} \propto |F_{hkl}|^2$

Larger $|F_{hkl}|$ means a stronger peak. When $F_{hkl} = 0$ exactly, you get a systematic absence: zero intensity at that reflection due to complete destructive interference among the scattered waves.

Derivation for crystals

The structure factor formula follows naturally from treating the crystal as a periodic array of identical unit cells and applying Fourier analysis to the scattering process.

Periodic arrangement of atoms

Crystals have translational symmetry: the same unit cell repeats in three dimensions. Because of this periodicity, you only need to know the contents of one unit cell to describe the entire crystal. The lattice periodicity is what restricts diffraction to discrete spots (the Laue condition), and the unit cell contents determine the intensity at each spot via the structure factor.

Fourier transform of atomic positions

The diffraction pattern is the Fourier transform of the electron density distribution in the crystal. In integral form, the structure factor is:

$F_{hkl} = \int_V \rho(\mathbf{r}) \exp[2\pi i(\mathbf{G}_{hkl} \cdot \mathbf{r})] \, d\mathbf{r}$

• $\rho(\mathbf{r})$ is the electron density at position $\mathbf{r}$
• $\mathbf{G}_{hkl}$ is the reciprocal lattice vector corresponding to $(h, k, l)$
• $V$ is the unit cell volume

When you approximate the electron density as a sum of individual atomic contributions centered at each atomic site, this integral reduces to the discrete sum formula given earlier. That's the connection: the sum over atoms is just the discretized version of the Fourier transform of $\rho(\mathbf{r})$.

Calculation methods

There are two main approaches to computing structure factors, and which one you'd use depends on the complexity of the structure.

Direct summation

1. List every atom in the unit cell with its fractional coordinates $(x_j, y_j, z_j)$ and form factor $f_j$.
2. For a given reflection $(hkl)$, compute the phase $\phi_j = 2\pi(hx_j + ky_j + lz_j)$ for each atom.
3. Sum up $f_j \exp(i\phi_j)$ over all atoms.

This is straightforward and works well for small unit cells. For structures with hundreds or thousands of atoms per unit cell, it becomes computationally expensive.

Fourier transform approach

1. Construct the electron density $\rho(\mathbf{r})$ on a grid from the known atomic positions and form factors.
2. Apply a Fourier transform to $\rho(\mathbf{r})$ to obtain $F_{hkl}$ for all reflections simultaneously.
3. Use fast Fourier transform (FFT) algorithms to speed up the computation.

This approach scales much better for large, complex structures and is the standard in modern crystallographic software.

Dependence on atomic form factors

The atomic form factor $f_j$ quantifies how strongly a single atom scatters radiation. It's a critical ingredient in the structure factor because it weights each atom's contribution.

Definition of atomic form factor

The form factor is the Fourier transform of the electron density of an isolated atom. It describes scattering amplitude as a function of scattering angle (or equivalently, $\sin\theta / \lambda$). Each element has a characteristic form factor curve.

Relationship to electron density

• At zero scattering angle, $f_j$ equals the total number of electrons in the atom ($f_j(0) = Z$).
• As the scattering angle increases, $f_j$ decreases because the finite spatial extent of the electron cloud causes partial destructive interference.
• Core electrons, which are tightly localized, contribute to scattering at higher angles more than diffuse valence electrons do.

This angle dependence means that high-angle reflections are generally weaker, even before accounting for thermal effects.

Systematic absences

Systematic absences are reflections where $F_{hkl} = 0$ for every crystal of that structure type, not by accident but because of symmetry. They're one of the most useful features of diffraction patterns for identifying crystal symmetry.

Conditions for zero intensity

The structure factor vanishes when the symmetry of the unit cell forces exact cancellation among the scattered waves. The specific rules depend on the symmetry elements present:

• Body-centered cubic (BCC): $F_{hkl} = 0$ when $h + k + l$ is odd. This happens because the atom at $(0,0,0)$ and the one at $(1/2, 1/2, 1/2)$ scatter exactly out of phase for those reflections.
• Face-centered cubic (FCC): $F_{hkl} = 0$ unless $h, k, l$ are all even or all odd.
• Glide planes and screw axes produce their own characteristic absence conditions along specific rows or planes of reciprocal space.

Connection to crystal symmetry

Different space groups produce different patterns of systematic absences. By cataloging which reflections are missing from your diffraction data, you can narrow down the space group of the crystal. This is typically one of the first steps in structure determination.

Applications in crystallography

Structure determination

Structure factor amplitudes $|F_{hkl}|$ are extracted from measured intensities. To reconstruct the electron density, you also need the phases. Methods like Patterson synthesis (which uses $|F|^2$ and doesn't require phases) or direct methods (which estimate phases statistically) are used to get an initial structural model.

Refinement of atomic positions

Once you have an approximate structure, refinement adjusts atomic positions, occupancies, and displacement parameters to minimize the difference between calculated and observed structure factors. This is typically done by least-squares fitting. The quality of the fit is tracked by R-factors, which compare $|F_{\text{obs}}|$ and $|F_{\text{calc}}|$ across all measured reflections.

Temperature effects

Atoms in a crystal aren't stationary. They vibrate around their equilibrium positions, and this thermal motion smears out the electron density, reducing diffraction peak intensities.

Debye-Waller factor

The Debye-Waller factor (also called the temperature factor or atomic displacement parameter) accounts for this thermal smearing. It modifies each atom's form factor:

$f_j \rightarrow f_j \exp\!\left(-B_j \frac{\sin^2\theta}{\lambda^2}\right)$

Here $B_j$ is related to the mean-square displacement $\langle u_j^2 \rangle$ of atom $j$ from its equilibrium site: $B_j = 8\pi^2 \langle u_j^2 \rangle$.

Influence on peak intensity

• The exponential decay means high-angle reflections (large $\sin\theta/\lambda$) are suppressed more than low-angle ones.
• Higher temperatures increase $\langle u^2 \rangle$, making the falloff steeper.
• Correcting for this effect is essential for accurate structure determination. If you ignore it, your refined atomic positions and electron density maps will be distorted.

Experimental measurement

X-ray diffraction techniques

X-rays scatter off the electron density in the crystal. Single-crystal X-ray diffraction gives a full three-dimensional set of $|F_{hkl}|$ values and is the gold standard for structure determination. Powder X-ray diffraction compresses this into a one-dimensional pattern where peaks at the same $d$-spacing overlap, so extracting individual structure factors requires additional deconvolution methods (like Rietveld refinement).

Neutron diffraction techniques

Neutrons scatter off atomic nuclei rather than electrons, which gives them distinct advantages:

• They can distinguish between elements with similar electron counts (or even different isotopes of the same element).
• They're much better at locating light atoms like hydrogen, which are nearly invisible to X-rays.
• They can probe magnetic structure because neutrons have a magnetic moment that interacts with unpaired electrons.

Neutron and X-ray diffraction are complementary, and using both together gives a more complete picture of the crystal structure.

Interpretation of structure factor

Phase information

Diffraction experiments measure intensities, which are proportional to $|F_{hkl}|^2$. The phase of $F_{hkl}$ is lost in this process. This is the famous phase problem in crystallography: you need both amplitude and phase to reconstruct the electron density, but you can only directly measure the amplitude.

Several strategies exist to recover or estimate phases:

• Patterson methods: Work with $|F|^2$ directly to find interatomic vectors, useful for locating heavy atoms.
• Direct methods: Use statistical relationships among structure factors to estimate phases. This works well for small-to-medium structures.
• Molecular replacement: If a similar structure is already known, its phases can serve as a starting model.

Electron density maps

With both amplitudes and phases in hand, you calculate the electron density via inverse Fourier transform:

$\rho(\mathbf{r}) = \frac{1}{V} \sum_{hkl} F_{hkl} \exp[-2\pi i(\mathbf{G}_{hkl} \cdot \mathbf{r})]$

The resulting 3D map shows where electrons are concentrated, and atomic positions appear as peaks in this map. The structure is then refined iteratively: adjust the model, recalculate $F_{hkl}$, compare with experiment, and repeat until convergence.