Conduction Electrons in Metal

The conduction electrons in a metal are non-localized (i.e., they are not tied to any particular atoms). In conventional metals, each atom contributes a fixed number of such electrons (corresponding to its valency). To a first approximation, it is possible to neglect the mutual interaction of the conduction electrons, because this interaction is largely shielded out by the stationary ions. The conduction electrons can, therefore, be treated as an ideal gas. However, the concentration of such electrons in a metal far exceeds the concentration of particles in a conventional gas. It is, therefore, not surprising that conduction electrons cannot normally be analyzed using classical statistics: in fact, they are subject to Fermi-Dirac statistics (because electrons are fermions).

Recall, from Section 8.5, that the mean number of particles occupying state $s$ (energy $\epsilon_s$ ) is given by

$$\bar{n}_s = \frac{1}{{\rm e}^{ \beta (\epsilon_s-\mu)} + 1},$$ (8.134)

according to the Fermi-Dirac distribution. Here,

$$\mu \equiv -k T \alpha$$ (8.135)

is termed the Fermi energy of the system. This energy is determined by the condition that

$$\sum_r \bar{n}_r = \sum_r\frac{1}{{\rm e}^{ \beta (\epsilon_r-\mu)} + 1} =N,$$ (8.136)

where $N$ is the total number of particles contained in the volume $V$ . It is clear, from the previous equation, that the Fermi energy, $\mu$ , is generally a function of the temperature, $T$ .

Let us investigate the behavior of the so-called Fermi function,

$$F(\epsilon) = \frac{1}{{\rm e}^{ \beta (\epsilon-\mu)} + 1},$$ (8.137)

as $\epsilon$ varies. Here, the energy is measured from its lowest possible value $\epsilon=0$ . If the Fermi energy, $\mu$ , is such that $\beta \mu\ll -1$ then $\beta (\epsilon-\mu)\gg 1$ , and $F$ reduces to the Maxwell-Boltzmann distribution. However, for the case of conduction electrons in a metal, we are interested in the opposite limit, where

$$\beta \mu \equiv \frac{\mu}{k T} \gg 1.$$ (8.138)

In this limit, if $0<\epsilon\ll \mu$ then $\beta (\epsilon-\mu)\ll -1$ , so that $F(\epsilon)=1$ . On the other hand, if $\epsilon\gg \mu$ then $\beta (\epsilon-\mu)\gg 1$ , so that $F(\epsilon)= \exp[-\beta (\epsilon-\mu)]$ falls off exponentially with increasing $\epsilon$ , just like a classical Maxwell-Boltzmann distribution. Note that $F=1/2$ when $\epsilon=\mu$ . The transition region in which $F$ goes from a value close to unity to a value close to zero corresponds to an energy interval of order $k T$ , centered on $\epsilon=\mu$ . In fact, $F=3/4$ when $\epsilon=\mu-(\ln 3) k T$ , and $F=1/4$ when $\epsilon=\mu+(\ln 3) k T$ . The behavior of the Fermi function is illustrated in Figure 8.4.

Figure 8.4: The Fermi function.

In the limit as $T\rightarrow 0$ , the transition region becomes infinitesimally narrow. In this case, $F=1$ for $\epsilon<\mu$ , and $F=0$ for $\epsilon>\mu$ , as illustrated in Figure 8.4. This is an obvious result, because when $T=0$ the conduction electrons attain their lowest energy, or ground-state, configuration. Because the Pauli exclusion principle requires that there be no more than one electron per single-particle quantum state, the lowest energy configuration is obtained by piling electrons into the lowest available unoccupied states, until all of the electrons are used up. Thus, the last electron added to the pile has a quite considerable energy, $\epsilon=\mu$ , because all of the lower energy states are already occupied. Clearly, the exclusion principle implies that a Fermi-Dirac gas possesses a large mean energy, even at absolute zero.

Let us calculate the Fermi energy, $\mu=\mu_F$ , of a Fermi-Dirac gas at $T=0$ . The energy of each particle is related to its momentum ${\bf p} =\hbar {\bf k}$ via

$$\epsilon = \frac{p^{ 2}}{2 m} =\frac{\hbar^{ 2} k^{ 2}}{2 m},$$ (8.139)

where ${\bf k}$ is the de Broglie wavevector. Here, $m$ is the electron mass. At $T=0$ , all quantum states whose energy is less than the Fermi energy, $\mu_F$ , are filled. The Fermi energy corresponds to a so-called Fermi momentum, $p_F=\hbar k_F$ , which is such that

$$\mu_F = \frac{p_F^{ 2}}{2 m} = \frac{\hbar^{ 2} k_F^{ 2}}{2 m}.$$ (8.140)

Thus, at $T=0$ , all quantum states with $k<k_F$ are filled, and all those with $k>k_F$ are empty.

Now, we know, by analogy with Equation (7.183), that there are $(2\pi)^{-3} V$ allowable translational states per unit volume of ${\bf k}$ -space. The volume of the sphere of radius $k_F$ in ${\bf k}$ -space is $(4/3) \pi k_F^{ 3}$ . It follows that the Fermi sphere of radius $k_F$ contains $(4/3) \pi k_F^{ 3} (2\pi)^{-3} V$ translational states. The number of quantum states inside the sphere is twice this, because electrons possess two possible spin states for every possible translational state. Because the total number of occupied states (i.e., the total number of quantum states inside the Fermi sphere) must equal the total number of particles in the gas, it follows that

$$2 \frac{V}{(2\pi)^3}\left(\frac{4}{3} \pi k_F^{ 3}\right) = N.$$ (8.141)

The previous expression can be rearranged to give

$$k_F = \left(3\pi^{ 2} \frac{N}{V}\right)^{1/3}.$$ (8.142)

Hence,

$$\lambda_F \equiv \frac{2\pi}{k_F} =\frac{2\pi}{(3\pi^{ 2})^{1/3}} \left(\frac{V}{N}\right)^{1/3},$$ (8.143)

which implies that the de Broglie wavelength, $\lambda_F$ , corresponding to the Fermi energy, is of order the mean separation between particles $(V/N)^{1/3}$ . All quantum states with de Broglie wavelengths $\lambda\equiv 2\pi/k > \lambda_F$ are occupied at $T=0$ , whereas all those with $\lambda<\lambda_F$ are empty.

According to Equations (8.140) and (8.142), the Fermi energy at $T=0$ takes the form

$$\mu_F = \frac{\hbar^{ 2}}{2 m}\left(3\pi^{ 2} \frac{N}{V}\right)^{2/3}.$$ (8.144)

It is easily demonstrated that $\mu_F\gg k T$ for conventional metals at room temperature. (See Exercise 17.)

The majority of the conduction electrons in a metal occupy a band of completely filled states with energies far below the Fermi energy. In many cases, such electrons have very little effect on the macroscopic properties of the metal. Consider, for example, the contribution of the conduction electrons to the specific heat of the metal. The heat capacity at constant volume, $C_V$ , of these electrons can be calculated from a knowledge of their mean energy, $\overline{E}$ , as a function of $T$ : that is,

$$C_V = \left(\frac{\partial \overline{E}}{\partial T}\right)_V.$$ (8.145)

If the electrons obeyed classical Maxwell-Boltzmann statistics, so that $F\propto \exp(-\beta \epsilon)$ for all electrons, then the equipartition theorem would give

$$\overline{E} = \frac{3}{2} N k T,$$ (8.146)

$$C_V =\frac{3}{2} N k.$$ (8.147)

However, the actual situation, in which $F$ has the form shown in Figure 8.4, is very different. A small change in $T$ does not affect the mean energies of the majority of the electrons, with $\epsilon\ll \mu$ , because these electrons lie in states that are completely filled, and remain so when the temperature is changed. It follows that these electrons contribute nothing whatsoever to the heat capacity. On the other hand, the relatively small number of electrons, $N_{\rm eff}$ , in the energy range of order $k T$ , centered on the Fermi energy, in which $F$ is significantly different from 0 and 1, do contribute to the specific heat. In the tail end of this region, $F\propto \exp(-\beta \epsilon)$ , so the distribution reverts to a Maxwell-Boltzmann distribution. Hence, from Equation (8.147), we expect each electron in this region to contribute roughly an amount $(3/2) k$ to the heat capacity. Hence, the heat capacity can be written

$$C_V\simeq \frac{3}{2} N_{\rm eff} k.$$ (8.148)

However, because only a fraction $k T/\mu$ of the total conduction electrons lie in the tail region of the Fermi-Dirac distribution, we expect

$$N_{\rm eff} \simeq \frac{k T}{\mu} N.$$ (8.149)

It follows that

$$C_V \simeq \frac{3}{2} N k \frac{k T}{\mu}.$$ (8.150)

Because $k T\ll \mu$ in conventional metals, the molar specific heat of the conduction electrons is clearly very much less than the classical value $(3/2) R$ . This accounts for the fact that the molar specific heat capacities of metals at room temperature are about the same as those of insulators. Before the advent of quantum mechanics, the classical theory predicted incorrectly that the presence of conduction electrons should raise the heat capacities of monovalent metals by 50 percent [i.e., $(3/2) R$ ] compared to those of insulators.

Note that the specific heat (8.150) is not temperature independent. In fact, using the superscript $e$ to denote the electronic specific heat, the molar specific heat can be written

$$c_V^{(e)} = \gamma T,$$ (8.151)

where $\gamma$ is a (positive) constant of proportionality. At room temperature $c_V^{(e)}$ is completely masked by the much larger specific heat, $c_V^{(L)}$ , due to lattice vibrations. However, at very low temperatures $c_V^{(L)}=A T^{ 3}$ , where $A$ is a (positive) constant of proportionality. (See Section 7.14.) Clearly, at low temperatures, $c_V^{(L)}=A T^{ 3}$ approaches zero far more rapidly that the electronic specific heat, as $T$ is reduced. Hence, it should be possible to measure the electronic contribution to the molar specific heat at low temperatures.

The total molar specific heat of a metal at low temperatures takes the form

$$c_V = c_V^{(e)} + c_V^{(L)} = \gamma T + A T^{ 3}.$$ (8.152)

Hence,

$$\frac{c_V}{T} = \gamma + A T^{ 2}.$$ (8.153)

It follows that a plot of $c_V/T$ versus $T^{ 2}$ should yield a straight-line whose intercept on the vertical axis gives the coefficient $\gamma$ . Figure 8.5 shows such a plot. The fact that a good straight-line is obtained verifies that the temperature dependence of the heat capacity predicted by Equation (8.152) is indeed correct.

Figure: The low-temperature heat capacity of potassium, plotted as $c_V/T$ versus $T^{ 2}$ . The straight-line shows the fit $c_V/T = 2.08+2.57 T^{ 2}$ . From C. Kittel, and H. Kroemer, Thermal Physics (W.H. Freeman & co., New York NY, 1980).