# Free electron model

(Redirected from Electron gas)
In three dimensions, the density of states of a gas of fermions is proportional to the square root of the kinetic energy of the particles.

In solid-state physics, the free electron model is a simple model for the behaviour of valence electrons in a crystal structure of a metallic solid. It was developed principally by Arnold Sommerfeld who combined the classical Drude model with quantum mechanical Fermi–Dirac statistics and hence it is also known as the Drude–Sommerfeld model.

The free electron empty lattice approximation forms the basis of the band structure model known as nearly free electron model. Given its simplicity, it is surprisingly successful in explaining many experimental phenomena, especially

## Ideas and assumptions

As in the Drude model, valence electrons are assumed to be completely detached from their ions (forming an electron gas). As in an ideal gas, electron-electron interactions are completely neglected. The electrostatic fields in metals are weak because of the screening effect.

The crystal lattice is not explicitly taken into account. A quantum-mechanical justification is given by Bloch's Theorem: an unbound electron moves in a periodic potential as a free electron in vacuum, except for the electron mass m becoming an effective mass m* which may deviate considerably from m (one can even use negative effective mass to describe conduction by electron holes). Effective masses can be derived from band structure computations. While the static lattice does not hinder the motion of the electrons, electrons can be scattered by impurities and by phonons; these two interactions determine electrical and thermal conductivity (superconductivity requires a more refined theory than the free electron model).

According to the Pauli exclusion principle, each phase space element (Δk)3(Δx)3 can be occupied only by two electrons (one per spin quantum number). This restriction of available electron states is taken into account by Fermi–Dirac statistics (see also Fermi gas). Main predictions of the free-electron model are derived by the Sommerfeld expansion of the Fermi–Dirac occupancy for energies around the Fermi level.

## Energy and wave function of a free electron

Plane wave traveling in the x-direction

For a free particle the potential is ${\displaystyle V({\mathbf {r}})=0}$. The Schrödinger equation for a free particle — in our case an electron — is[1][2][3]

${\displaystyle -{\frac {\hbar ^{2}}{2m}}\nabla ^{2}\Psi ({\mathbf {r}},t)=i\hbar {\frac {\partial }{\partial t}}\Psi ({\mathbf {r}},t)}$

The Schrödinger equation can be solved by separation of variables (see below) to yield a plane wave solution

${\displaystyle \Psi ({\mathbf {r}},t)={\frac {1}{\sqrt {\Omega _{r}}}}e^{i{\mathbf {k}}\cdot {\mathbf {r}}-i\omega t},}$

where ${\displaystyle {\mathbf {k}}}$ is the wave vector, ${\displaystyle \omega }$ is the frequency, and ${\displaystyle \Omega _{r}}$ is the volume of space where the electron can be found. The electron has kinetic energy

${\displaystyle E=\hbar \omega ={\frac {\hbar ^{2}k^{2}}{2m}}}$

which relates ${\displaystyle \omega }$ and ${\displaystyle {\mathbf {k}}}$. For solid state and condensed matter physics, the spatial form

${\displaystyle \psi _{\mathbf {k}}({\mathbf {r}})={\frac {1}{\sqrt {\Omega _{r}}}}e^{i{\mathbf {k}}\cdot {\mathbf {r}}}}$

of the plane wave is of major interest. It is the basis of electronic band structure models that are widely used in solid-state physics for model Hamiltonians like the nearly free electron model, the tight binding model and various models that use a muffin-tin approximation. The eigenfunctions of these Hamiltonians are Bloch waves which are modulated plane waves.

## Dielectric function of the electron gas

On a scale much larger than the inter atomic distance a solid can be viewed as an aggregate of a negatively charged plasma of the free electron gas and a positively charged background of atomic cores. The background is the rather stiff and massive background of atomic nuclei and core electrons which we will consider to be infinitely massive and fixed in space. The negatively charged plasma is formed by the valence electrons of the free electron model that are uniformly distributed over the interior of the solid. If an oscillating electric field is applied to the solid, the negatively charged plasma tends to move a distance x apart from the positively charged background. As a result, the sample is polarized and there will be an excess charge at the opposite surfaces of the sample. The surface charge density is

${\displaystyle \rho _{s}=-nex}$

where n is the number density of electrons. This produces a restoring electric field in the sample

${\displaystyle E={\frac {nex}{\epsilon _{0}}}}$

The dielectric constant of the sample at a fixed frequency ${\displaystyle \omega }$ is expressed as

${\displaystyle \epsilon ={\frac {D}{\epsilon _{0}E}}=1+{\frac {P}{\epsilon _{0}E}}}$

where ${\displaystyle D}$ is the electric displacement and ${\displaystyle P}$ is the polarization density.

The electric field and polarization densities are

${\displaystyle E(t)=E_{0}e^{-i\omega t},\quad P(t)=P_{0}e^{-i\omega t}}$

and the polarization density with n electron density is

${\displaystyle P=-nex}$

The force F of the oscillating electric field causes the electrons with charge e and mass m to accelerate with an acceleration a

${\displaystyle F=ma=m{\frac {d^{2}x}{dt^{2}}}=-eE}$

which, after substitution of E, P and x, yields a harmonic oscillator equation.

After a little algebra the relation between polarization density and electric field can be expressed as

${\displaystyle P=-{\frac {ne^{2}}{m\omega ^{2}}}E}$

The frequency dependent dielectric function of the solid is

${\displaystyle \epsilon (\omega )=1-{\frac {ne^{2}}{\epsilon _{0}m\omega ^{2}}}}$

At a resonance frequency ${\displaystyle \omega _{p}}$, called the plasma frequency, the dielectric function changes sign from negative to positive and real part of the dielectric function drops to zero.

${\displaystyle \omega _{p}={\sqrt {\frac {ne^{2}}{\epsilon _{0}m}}}}$

The plasma frequency represents a plasma oscillation resonance or plasmon. The plasma frequency can be employed as a direct measure of the square root of the density of valence electrons in a solid. Observed values are in reasonable agreement with this theoretical prediction for a large number of materials.[4] Below the plasma frequency, the dielectric function is negative and the field cannot penetrate the sample. Light with angular frequency below the plasma frequency will be totally reflected. Above the plasma frequency the light waves can penetrate the sample.

## Solution of the Schrödinger equation

### The Schrödinger equation

For a free particle the potential is ${\displaystyle V({\mathbf {r}})=0}$, so the Schrödinger equation for the free electron is[1][2][3]

${\displaystyle -{\frac {\hbar ^{2}}{2m}}\nabla ^{2}\Psi ({\mathbf {r}},t)=i\hbar {\frac {\partial }{\partial t}}\Psi ({\mathbf {r}},t)}$

This is a type of wave equation that has numerous kinds of solutions. One way of solving the equation is splitting it in a time-dependent oscillator equation and a space-dependent wave equation like

${\displaystyle i\hbar {\frac {\partial }{\partial t}}\Psi ({\mathbf {r}},t)=E\Psi ({\mathbf {r}},t)}$

and

${\displaystyle -{\frac {\hbar ^{2}}{2m}}\nabla ^{2}\Psi ({\mathbf {r}},t)=E\Psi ({\mathbf {r}},t)}$

and substituting a product of solutions like

${\displaystyle \Psi ({\mathbf {r}},t)=\psi ({\mathbf {r}})e^{\alpha t}}$

The Schrödinger equation can be split in a time dependent part and a time independent part.Which is derived.

### Solution of the time dependent equation

The peculiar time dependent part of the Schrödinger equation is, unlike the Klein–Gordon equation for pions and most of the other well known wave equations, a first order in time differential equation with a 90° out of phase driving mechanism, while most oscillator equations are second order in time differential equations with 180° out of phase driving mechanisms.

The equation that has to be solved is

${\displaystyle i\hbar {\frac {\partial }{\partial t}}e^{\alpha t}=Ee^{\alpha t}}$.

The complex (imaginary) exponent is proportional to the energy

${\displaystyle \alpha =-{\frac {iE}{\hbar }}}$

The imaginary exponent can be transformed to an angular frequency

${\displaystyle E=i\hbar \alpha =\hbar \omega }$

The wave function now has a stationary and an oscillating part

${\displaystyle \Psi ({\mathbf {r}},t)=\psi ({\mathbf {r}})e^{-i\omega t}}$

The stationary part is of major importance to the physical properties of the electronic structure of matter.

### Solution of the time independent equation

The wave function of free electrons is in general described as the solution of the time independent Schrödinger equation for free electrons

${\displaystyle -{\frac {\hbar ^{2}}{2m}}\nabla ^{2}\psi ({\mathbf {r}})=E\psi ({\mathbf {r}})}$
${\displaystyle \nabla ^{2}={\frac {\partial ^{2}}{{\partial x}^{2}}}+{\frac {\partial ^{2}}{{\partial y}^{2}}}+{\frac {\partial ^{2}}{{\partial z}^{2}}}}$

The wave function can be factorized for the three Cartesian directions

${\displaystyle \psi ({\mathbf {r}})=\phi _{x}(x)\phi _{y}(y)\phi _{z}(z)}$

Now the time independent Schrödinger equation can be split in three independent parts for the three different Cartesian directions

${\displaystyle -{\frac {\hbar ^{2}}{2m}}{\frac {\partial ^{2}}{{\partial x}^{2}}}\phi _{x}(x)=E_{x}\phi _{x}(x)}$

As a solution an exponential function is substituted in the time independent Schrödinger equation

${\displaystyle \phi _{x}(x)=N_{x}e^{\kappa x}}$

The solution of

${\displaystyle {\frac {\partial ^{2}}{{\partial x}^{2}}}\phi _{x}(x)=\kappa ^{2}N_{x}e^{\kappa x}=-{\frac {2m}{\hbar ^{2}}}E_{x}N_{x}e^{\kappa x}}$

gives the exponent

${\displaystyle \kappa =ik_{x}=i{\sqrt {\frac {2mE_{x}}{\hbar ^{2}}}}}$

which yields the wave equation

${\displaystyle \psi ({\mathbf {r}})=N_{x}N_{y}N_{z}e^{i(k_{x}x+k_{y}y+k_{z}z)}=N_{\mathbf {k}}e^{i{\mathbf {k}}\cdot {\mathbf {r}}}}$

and the energy

${\displaystyle E={\frac {\hbar ^{2}}{2m}}(k_{x}^{2}+k_{y}^{2}+k_{z}^{2})}$

With the normalization

${\displaystyle 1=\int _{\Omega _{r}}|\psi _{\mathbf {k}}({\mathbf {r}})|^{2}d{\mathbf {r}}=\int _{\Omega _{r}}\psi _{\mathbf {k}}^{*}({\mathbf {r}})\psi _{{\mathbf {k}}'}({\mathbf {r}})d{\mathbf {r}}=\int _{\Omega _{r}}(N_{\mathbf {k}}^{*}e^{-i{\mathbf {k}}\cdot {\mathbf {r}}})(N_{\mathbf {k}}e^{i{\mathbf {k}}\cdot {\mathbf {r}}})d{\mathbf {r}}=N_{\mathbf {k}}^{2}|\Omega _{r}|_{vol}}$

and the wave vector magnitude

${\displaystyle k={\sqrt {k_{x}^{2}+k_{y}^{2}+k_{z}^{2}}}}$

we arrive at the plane wave solution with a wave function

${\displaystyle \psi _{\mathbf {k}}({\mathbf {r}})={\frac {1}{\sqrt {\Omega _{r}}}}e^{i{\mathbf {k}}\cdot {\mathbf {r}}}}$

for free electrons with a wave vector ${\displaystyle {\mathbf {k}}}$ and a kinetic energy

${\displaystyle E={\frac {\hbar ^{2}k^{2}}{2m}}}$

in which ${\displaystyle \Omega _{r}}$ is the volume of space occupied by the electron.

### The traveling plane wave solution

The product of the time independent stationary wave solution and time dependent oscillator solution

${\displaystyle \Psi ({\mathbf {r}},t)=\psi ({\mathbf {r}})e^{-i\omega t}}$

gives the traveling plane wave solution

${\displaystyle \Psi ({\mathbf {r}},t)={\frac {1}{\sqrt {\Omega _{r}}}}e^{i{\mathbf {k}}\cdot {\mathbf {r}}-i\omega t}}$

which is the final solution for the free electron wave function.

## Fermi energy

According to the Pauli principle, the electrons in the ground state occupy all the lowest-energy states, up to some Fermi energy ${\displaystyle E_{F}}$. Since the energy is given by

${\displaystyle E({\mathbf {k}})={\frac {\hbar ^{2}k^{2}}{2m}}}$,

this corresponds to occupying all the states with wave vectors ${\displaystyle |{\mathbf {k}}|, where ${\displaystyle k_{F}}$ is so-called Fermi wave vector, given by

${\displaystyle k_{F}=(3\pi ^{2}N_{e}/V)^{1/3}}$,

where ${\displaystyle N_{e}}$ is the total number of electrons in the system, and V is the total volume. The Fermi energy is then

 ${\displaystyle E_{F}={\frac {\hbar ^{2}}{2m}}\left({\frac {3\pi ^{2}N_{e}}{V}}\right)^{2/3}\,}$

In a nearly-free-electron model of a ${\displaystyle Z}$-valent metal, one can replace ${\displaystyle N_{e}}$ with ${\displaystyle NZ}$, where ${\displaystyle N}$ is the total number of metal ions.

## Density of states

The density of states (DOS) corresponds to electrons with a spherically-symmetric parabolic dispersion

${\displaystyle E({\mathbf {k}})={\frac {\hbar ^{2}k^{2}}{2m}}}$,

with two electrons (one of each spin) per each "quantum" of the phase space, ${\displaystyle \Delta V\Delta {\mathbf {k}}=(2\pi )^{3}}$ . In 3D, this corresponds to

${\displaystyle N(E)={\frac {V}{2\pi ^{2}}}\left({\frac {2m}{\hbar ^{2}}}\right)^{3/2}{\sqrt {E}}}$,

where ${\displaystyle V}$ is the total volume.

Combining these expressions for the Fermi energy and the DOS, one can show that the following relationship holds at the Fermi level:

${\displaystyle N(E_{F})={\frac {3ZN}{2}}{\frac {1}{E_{F}}},}$

where Z is the charge of each of the N metal ions in the crystal.