Nearly free electron model

In solid-state physics, the nearly free electron model (or NFE model and quasi-free electron model) is a quantum mechanical model of physical properties of electrons that can move almost freely through the crystal lattice of a solid. The model is closely related to the more conceptual empty lattice approximation. The model enables understanding and calculation of the electronic band structures, especially of metals.

This model is an immediate improvement of the free electron model, in which the metal was considered as a non-interacting electron gas and the ions were neglected completely.

Mathematical formulation

edit
Dispersion relation for the 2D nearly free electron model as a function of the underlying crystalline structure.

The nearly free electron model is a modification of the free-electron gas model which includes a weak periodic perturbation meant to model the interaction between the conduction electrons and the ions in a crystalline solid. This model, like the free-electron model, does not take into account electron–electron interactions; that is, the independent electron approximation is still in effect.

As shown by Bloch's theorem, introducing a periodic potential into the Schrödinger equation results in a wave function of the form

 

where the function   has the same periodicity as the lattice:

 

(where   is a lattice translation vector.)

Because it is a nearly free electron approximation we can assume that

  where   denotes the volume of states of fixed radius   (as described in Gibbs paradox).[clarification needed]

A solution of this form can be plugged into the Schrödinger equation, resulting in the central equation:

 

where   is the total energy, and the kinetic energy   is characterized by

 

which, after dividing by  , reduces to

 

if we assume that   is almost constant and  

The reciprocal parameters   and   are the Fourier coefficients of the wave function   and the screened potential energy  , respectively:

   

The vectors   are the reciprocal lattice vectors, and the discrete values of   are determined by the boundary conditions of the lattice under consideration.

Before doing the perturbation analysis, let us first consider the base case to which the perturbation is applied. Here, the base case is  , and therefore all the Fourier coefficients of the potential are also zero. In this case the central equation reduces to the form

 

This identity means that for each  , one of the two following cases must hold:

  1.  ,
  2.  

If   is a non-degenerate energy level, then the second case occurs for only one value of  , while for the remaining  , the Fourier expansion coefficient   is zero. In this case, the standard free electron gas result is retrieved:

 

If   is a degenerate energy level, there will be a set of lattice vectors   with  . Then there will be   independent plane wave solutions of which any linear combination is also a solution:

 

Now let   be nonzero and small. Non-degenerate and degenerate perturbation theory, respectively, can be applied in these two cases to solve for the Fourier coefficients   of the wavefunction (correct to first order in  ) and the energy eigenvalue   (correct to second order in  ). An important result of this derivation is that there is no first-order shift in the energy   in the case of no degeneracy, while there is in the case of degeneracy (and near-degeneracy), implying that the latter case is more important in this analysis. Particularly, at the Brillouin zone boundary (or, equivalently, at any point on a Bragg plane), one finds a twofold energy degeneracy that results in a shift in energy given by:[clarification needed]

 .

This energy gap between Brillouin zones is known as the band gap, with a magnitude of  .

Results

edit

Introducing this weak perturbation has significant effects on the solution to the Schrödinger equation, most significantly resulting in a band gap between wave vectors in different Brillouin zones.

Justifications

edit

In this model, the assumption is made that the interaction between the conduction electrons and the ion cores can be modeled through the use of a "weak" perturbing potential. This may seem like a severe approximation, for the Coulomb attraction between these two particles of opposite charge can be quite significant at short distances. It can be partially justified, however, by noting two important properties of the quantum mechanical system:

  1. The force between the ions and the electrons is greatest at very small distances. However, the conduction electrons are not "allowed" to get this close to the ion cores due to the Pauli exclusion principle: the orbitals closest to the ion core are already occupied by the core electrons. Therefore, the conduction electrons never get close enough to the ion cores to feel their full force.
  2. Furthermore, the core electrons shield the ion charge magnitude "seen" by the conduction electrons. The result is an effective nuclear charge experienced by the conduction electrons which is significantly reduced from the actual nuclear charge.

See also

edit

References

edit
  • Ashcroft, Neil W.; Mermin, N. David (1976). Solid State Physics. Orlando: Harcourt. ISBN 0-03-083993-9.
  • Kittel, Charles (1996). Introduction to Solid State Physics (7th ed.). New York: Wiley. ISBN 0-471-11181-3.
  • Elliott, Stephen (1998). The Physics and Chemistry of Solids. New York: Wiley. ISBN 0-471-98194-X.