# 31.4: The Total Scattering Intensity is Related to the Periodic Structure of the Electron Density in the Crystal

- Page ID
- 14576

## Systematic Absences

Diffraction data from a face-centered cubic crystal will be missing reflections (spots) from all *hkl* planes in which *h + k*, *h + l*, and *k + l* are odd numbers. The reason for these absent reflections is that the arrangement of the lattice points leads to total destructive interference of the diffracted X-rays. Although the lack of data may seem like a hindrance to determining a crystal structure, the systematic absence of reflections from specific families of planes actually reveals the presence of symmetry features in the crystal lattice that help identify the lattice type. Crystals with centering features (body-centered cels, face-centered cells, etc.), glide planes, and screw axes will produce diffraction patterns with systematic absences.

## The Structure Factor

Systematic reflection absences, as well as variations in reflection intensity, are described mathematically by the structure factor, a mathematical function that incorporates all of the variables that affect how X-rays are scattered by the electrons of the atoms in the crystal lattice. (We do not concern ourselves with the effect of the nucleus on scattering because the mass of the nucleus is so large relative to the mass of an electron.)

### The Atomic/Ionic Scattering Factor

The construction of a structure factor begins with the fundamental assumption that the total scattering intensity is related to the electron density of the atoms in the crystal. But the scattering also depends on the scattering angle and the wavelength of the radiation. To account for all of these parameters, an **atomic scattering factor**, \(f\), has been developed:

\[f = 4\pi \int_0^{\infty} \rho(r)\dfrac{\sin \, kr}{kr}r^2 \, dr \nonumber\]

In this equation, \(r\) is the radius of an idealized spherical atom, \(\rho(r)\) is the electron density of this idealized atom. Because the X-ray wavelength used in diffraction experiments is of similar size as the size of the atoms it scatters from, there will be interference among the waves scattered by the electrons on a single atom. The more electrons in an atom and the more diffuse the electron cloud around the atom, the greater the destructive interference among the scattered waves. To account for this intra-atomic interference, the \(\dfrac{\sin \, kr}{kr}\) term is included in the scattering factor. Here, \(k\) is equal to \(\dfrac{4\pi \sin \theta}{\lambda}\), where \(\theta\) is the scattering angle and \(\lambda\) is the wavelength of the X-ray. A plot of \(f\) versus \(\dfrac{\sin \theta}{\lambda}\) is shown in figure \(\PageIndex{1}\). Notice that at \(\theta\) = 0, \(f\) is equal to the number of electrons in the particle. Also notice that, although the isoelectronic ions Ca^{2+} and Cl^{-} have an identical \(f\) value at \(\theta\) = 0, their \(f\) values differ at larger angles because the ions have different electron densities. Finally, notice that H atoms do not scatter a great deal, thus making them difficult to "locate" among the scattering from all of the larger atoms.

### Scattering by Neighboring Atoms and Neighboring Planes

Along with the scattering factor for individual atoms, we must also take into consideration the interferences caused by neighboring atoms in the lattice. As our example, we will analyze a crystal that contains atoms of two elements, here designated **P** and **Q**. We will assume that the respective atoms have scattering factors \(f\)_{P} and \(f\)_{Q}. We will also assume that the distance between successive P atoms is the same as the distance between successive Q atoms. This distance is \(\dfrac{a}{h}\), where \(a\) is the length of the unit cell along the **a** axis. Finally, we will assume that the distance between **P** and **Q** atoms is \(x\). In this arrangement, the **P** atoms lie in neighboring *hkl *planes, and the **Q** atoms lie in neighboring *hkl* planes that are interleaved between the *hkl* planes of the **P** atoms. Figure \(\PageIndex{2}\) shows the arrangement of the atoms lying along the **a** axis.

In this arrangement, the difference in path length of X-rays scattered by successive P atoms, \(\Delta\)_{PP}, (and the difference in path length of X-rays scattered by successive Q atoms, (\(\Delta\)_{QQ}) is described by the van Laue equation

\(\Delta\)_{PP }= \(\Delta\)QQ = \(\dfrac{a}{h}(cos\alpha - cos\alpha_0) = n\lambda\)

The difference in path length for X-rays diffracted by neighboring **P** and **Q** atoms is thus

\(\Delta\)PQ = \(x(\cos\alpha - \cos\alpha_0) = \dfrac{nh\lambda x}{a}\)

Because of this difference in path length, there is a resulting a phase difference between the scattered beams from neighboring **P** and **Q** atoms due to the difference in the time required to travel the different distances:

\(\phi = 2\pi \dfrac{\Delta_{PQ}}{\lambda} = 2\pi \dfrac{h\lambda x}{\lambda a} = \dfrac{2\pi h x}{a} \)

The light scattered from neighboring **P **and **Q** atoms will have an amplitude dependent on the scattering factors of each atom and the angular frequency of the X-ray radiation, \(\omega\)

\[A = f_P \, cos \omega t + f_Q \, cos(\omega t + \phi) \nonumber\]

Replacing the cosine functions with exponential functions will simplify the following mathematical calculation, thus

\[A = f_P \, e^{i\omega t} + f_Q \, e^{i(\omega t + \phi)} \label{1}\]

Because intensity is proportional to amplitude squared,

\( I \propto |A|^2 = (f_P \, e^{i\omega t} + f_Q \, e^{i(\omega t + \phi)})(f_P \, e^{-i\omega t} + f_Q \, e^{-i(\omega t + \phi)} \)

which can be simplified to

\(f_P^2 + f_Pf_Qe^{i\phi} + f_Pf_Qe^{-i\phi} + f_Q^2 \)

Switching back to cosine functions, we get

\(f_P^2 + f_Q^2 + 2 f_Pf_Qcos\phi \)

The \(f_P^2\) term represents the constructive interference of the x-rays scattered from the parallel planes through the **P** atoms. The \(f_Q^2\) term represents the constructive interference of the x-rays scattered from the parallel planes through the **Q** atoms. The \(2 f_Pf_Qcos\phi \) term represents the interference from the scattering from the two sets of interleaved **P** and **Q** planes. From this equation, we can see that the intensity of the scattered light is independent of the angular frequency of the X-ray radiation because there is no term involving \(\omega\). Thus, in equation \(\ref{1}\) the frequency-dependent \(e^{i\omega t}\) must have a value of 1, and the equation can be rewritten and redefined as the structure factor along the **a** axis of the crystal, \(F(h)\)

\[F(h) = f_P + f_Q \, e^{i\phi)} = f_P + f_Q \, e^{\dfrac{2\pi h x}{a}} \label{2}\]

For a unit cell with sides of length \(a, \, b,\) and \(c\) that contains atoms of type \(j\) positioned at points \(x_j, y_j \) and \(z_j\), equation \(\ref{2}\) becomes

\[F(hkl) = \sum_j f_{j} e^{2\pi (hx_{j}/a + ky_{j}/b + lz_{j}/c)} \label{3}\]

Here \(f_j\) is the scattering factor for type \(j\) atoms, and *hkl* are the Miller indicies of the diffracting planes. By convention, \(x_j, y_j, \) and \(z_j\) are expressed in units of \(a, b,\) and \(c\), so that \(x_j^{'} =\dfrac{x_j}{a}\), \(y_j^{'} =\dfrac{y_j}{b}\), and \(z_j^{'} =\dfrac{z_j}{c}\), to give

\[F(hkl) = \sum_j f_j e^{2\pi (hx_{j}^{'} + ky_{j}^{'} + lz_{j}^{'})} \label{4}\]

### Combining the Influences

\(F(hkl)\) is the structure factor of a crystal, and the intensity of reflections from this crystal are proportional to \(|F(hkl)|^2\). The strcuture factor can thus be used to explain the systematic absences described in the introduction to this section. Whenever \(F(hkl)\) = 0 for a set of Miller planes, the reflection will be absent.

Example \(\PageIndex{1}\)

Derive an expression for the structure factor of a face-centered cubic unit cell of identical atoms. Determine which* hkl* planes will not produce reflections.

**Solution**

The coordinates of the lattice points in a face-centered cubic unit cell are (0,0,0), (1,0,0), (0,1,0), (0,0,1), (1,0,1), (1,1,0), (1,1,1), (0,1,1), (0,1/2,1/2), (1/2,1/2,0), (1/2,0,1/2), (1,1/2,1/2), (1/2,1,1/2), (1/2,1/2,1). The length of the sides of the cell is \(a\). Each corner lattice point is shared by 8 cells, so only 1/8 of each of these atoms is considered to be contributing to the scattering. Each face lattice point is shared by 2 cells, so only 1/2 of each of these atoms is considered to be contributing to the scattering. We will use equation \(\ref{4}\) to solve the problem, remembering that \(a = b = c\) for a cubic cell:

\[F(hkl) = \dfrac{1}{8} f(e^{2\pi i(0+0+0)} + e^{2\pi i(h+0+0)} + e^{2\pi i(0+k+0)} + e^{2\pi i(0+0+l)} + e^{2\pi i(h+0+l)} + e^{2\pi i(h+k+0)} + e^{2\pi i(h+k+l)} + e^{2\pi i(0+k+l)}) + \dfrac{1}{2}f(e^{2\pi i(0+k/2+l/2)} + e^{2\pi i(h/2+k/2+0)}+ e^{2\pi i(h/2+0+l/2)} + e^{2\pi i(h+k/2+l/2)} + e^{2\pi i(h/2+k+l/2)} + e^{2\pi i(h/2+k/2+l)}) \label{5}\]

Because \(e^{2\pi i} \) = 1 and \(e^{\pi i}\) = -1, equation \(\ref{5}\) becomes

\[F(hkl) = \dfrac{1}{8}f(1^0 + 1^h + 1^k + 1^l + 1^{h+l} + 1^{h+k} + 1^{h+k+l} + 1^{k+l}) + \dfrac{1}{2}f((-1)^{k+l} + (-1)^{h+k} + (-1)^{h+l} + (-1)^{2h+k+l} + (-1)^{h+2k+l} + (-1)^{h+k+2l}) \nonumber\]

Because \(1^n = 1\) for all \(n\),

\[F(hkl) = \dfrac{1}{8}f(8) + \dfrac{1}{2}f((-1)^{k+l} + (-1)^{h+k} + (-1)^{h+l} + (-1)^{2h+k+l} + (-1)^{h+2k+l} + (-1)^{h+k+2l}) \nonumber\]

Upon testing a few *hkl* planes, you will see that if \(h,k,l\) are a mixture of odd and even numbers, then \(F(hkl) = 0\), but if \(h,k,l\) are all odd or all even, the value of \(F(hkl) = 4f\). Three examples:

All odd \(F(111) = f (1 + \dfrac{1}{2}((-1)^2 + (-1)^2 + (-1)^2 + (-1)^4 + (-1)^4 + (-1)^4 = f(1 + 3) = 4f \)

All even \(F(222) = f (1 + \dfrac{1}{2}((-1)^4 + (-1)^4 + (-1)^4 + (-1)^8 + (-1)^8 + (-1)^8 = f(1 + 3) = 4f \)

Mix \(F(110) = f (1 + \dfrac{1}{2}((-1)^1 + (-1)^2 + (-1)^1 + (-1)^3 + (-1)^3 + (-1)^2 = f(1 + \dfrac{1}{2}((-1) + 1 + (-1) + (-1) + (-1) + 1)) = 1+ (-1) = 0 \)

For crystals comprised of atoms of more than one element, there are often reflections that are weaker than other reflections. These lower intensities are sometimes created when the scattering factors of different elements are not identical, leading to incomplete destructive interference, a smaller \(F(hkl)\), and thus a reflection of lower intensity. The lower intensity reflections are also a direct result of the structure factor. For example, for an ionically-bonded compound with the formula MX, which crystallizes in a face-centered cubic lattice, the structure factor can be shown to be

\[f(hkl) = f_+[1 + (-1)^{h+k} + (-1)^{h+l} + (-1)^{k+l}] + f_-[(-1)^{h+k+l} + (-1)^{h} + (-1)^{k} + (-1)^{l}] \nonumber\]

Inspection of this equation reveals that \(F(hkl) = 4(f_+ + f_-)\) if \(h, k, l\) are all even, and \(F(hkl) = 4(f_+ - f_-)\)if \(h, k, l\) are all odd. Thus, reflections from the all-even \(hkl\) planes are more intense than reflections from the all-odd \(hkl\) planes.

## See also

The structure factor. P. Coppens. *International Tables for Crystallography* (2006). Vol. B, ch. 1.2, pp. 10-24

## Contributors

- Tom Neils - Grand Rapids Community College