Heat capacity of a 2D electron gas and a 2D phonon gas

Solution Preview

---

Electrons in metals form a degenerate Fermi gas, because the density of electrons is far larger than the number of quantum states per unit volume with energy below k T. Then, since each electron must occupy a different quantum state, most electrons will have energies far larger than k T. The thermal energy of k T which is available for thermal fluctuations is then only available for the electrons in the top energy levels around k T of the maximum energy level occupied. Only a small fraction of all the electrons are in these quantum states, so most electrons do not participate in thermal processes.

The density of states of a 2d Fermi gas can be computed as follows. The total number of quantum states in an area A and area Ap in momentum space is given by:

2 A Ap/h^2

Here the factor 2 comes from the two possible spin states the spin 1/2 particle can have. The area in momentum space between the magnitude of the momentum of P and P + dP is given by 2 pi P dP (i.e. the circumference times the width). The number of quantum states with momenta between P and P + dP is thus given by:

4 pi A/h^2 P dP

Since 2 P dP can be written as d(P^2), we can write this as:

2 pi A/h^2 d(P^2 )

Since P^2/(2m) is the energy E, we can rewrite this in terms of the range in energy dE:

4 pi m A/h^2 dE

The density of states is therefore:

4 pi m A/h^2 = m A/(pi hbar^2)

So, we see that the density of states does not depend in the energy in the 2d case. The total number of states below some arbitrary energy E is thus obtained by simply multiplying the density of states by E:

N(E) = m A/(pi hbar^2) E (1)

If we have M electrons at 0 K then these M electrons will occupy the M lowest energy levels. So, this means that the Fermi energy Ef will be such that N(Ef) = M:

M = m A/(pi hbar^2) Ef --------->

M/A = m/(pi hbar^2) Ef

Since the density of states does not depend on the energy, the electrons are uniformly distributed as a function of energy between zero and Ef, therefore the average energy of the electrons is Ef/2.

To compute the heat capacity of the 2d electron gas, we must consider the Fermi-Dirac distribution at finite temperature. We can write the total energy E and the number of particles N in terms of the temperature and the chemical potential as:

E = integral over e of e g(e) de/[exp(beta (e - mu)) + 1]

N = integral over e of g(e) de/[exp(beta (e - mu)) + 1]

Here

g(e) = m A/(pi hbar^2)

is the density of states obtained earlier in the problem. By "e" we denote the single particle electron energy (we use the lower case to avoid confusion with the total energy E). The integrals over e are from zero to infinity.

To compute the heat capacity near T = 0, we need to know how E behaves near zero. The above integral expression for E, gives E in terms of T and mu, but mu will depend on T also. So, we need to consider the expression for N in terms of mu and T also and solve for mu in terms of N and T.

So, let's first evaluate ...