4.5 The ”Tight-Binding”-Method: LCAO

Even the perfect crystal can not be calculated exactly since

Therefore we will discuss two approximations which describe two limits of a solid.
Expanding into plane waves
The Eigenvectors of free particles are plane waves. Describing a solid not taking into account the periodic lattice potential, we reach the model of a free electron gas:

ψnk(r) = eikrc nk. (4.27)

Since all constants cn obviously display the periodicity of the lattice, the Bloch-Theorem (4.2) is fulfilled.
This approximation can be improved significantly by taking into account a perturbation of the periodic lattice. We substitute Eq. (4.27) by

ψnk(r) = eikrc nk(r). (4.28)

The functions cnk(r) display the lattice periodicity and are constant on a defined set of intervals; they therefore obey the Bloch condition. This so-called Kronig-Penney model can explain several aspects, e.g.

This approach describes only these solids adequately which consist of almost free electrons. These electrons show a great mobility and will not be scattered effectively by the lattice atoms. We precisely described the properties of metal electrons (S-electrons) which exist e.g. in alkaline end earth-alkaline metals. The more plane waves are used for the description of the solid (the more accurate we approximate the lattice potential) the better we will describe the electronic state of these solids. The expansion into plane waves is quite simple on a computer. Taking into account more expansion coefficients (increasing the computational time) one gets quite good results even for systems for which the above requirements do not met optimally.
Expanding into atomic wave functions
Often plane waves are not at all a good approximation for the electronic states of a solid:



Ionic crystals strong influence of the lattice potential


Insulators no free electrons


Semiconductorsstrong covalent bonds, atomic like orbitals


Here we find a strong dominance of the atomic character of the electronic wave functions which defines the properties of the solid. Considering the bad conductivity of these solids is a strong hint that plane waves will not be an effective tool to approximate these solids:

bad electronic conductivity
no free electrons
no good description by plane waves

Therefore we will use an approximation which implements the lattice periodicity as a ”small” perturbation for the atomic orbitals.

CONCEPT:

  1. The lattice periodic potential is described by the function
    U(r) = U(r + R) (4.29)
  2. We expand at every lattice point R into atomic wave
    Hatψn = Enψn (4.30)

    These functions form for each point a CONS (Complete orthonormal system). To evaluate an arbitrary function one has to expand it in every point into the CONS. For a solid and using the above approximations we only have to expand the state function at the positions of atoms.

  3. The periodicity of the crystal is implemented by the Bloch theorem:
    ψ(r + R) = eikRψ(r) (4.31)

    All calculations can therefore be restricted to the first Brillouin zone.

The Hamiltonian for the first Brillouin zone (centered around the position of the atom R) can be written as

H = Hat + ΔU(r) (4.32)

The kinetic energy is completely included in Hat. ΔU(r) describes the additional influence of the neighboring atoms for the position R.
The atomic wave functions have a range which follows the above assumptions (as displayed in Fig. 4.1). Thus they do not overlap strongly with the orbitals of their neighboring atoms.
The complete wave function is summed up using the atomic wave functions of Eq. (4.30):

ψnk(r) = ReikRΦ n(r - R) (4.33)

PIC
Fig. 4.1 Along a series of atom positions the range of the atomic wave functions Φ, the periodic lattice potential U(r) and the perturbation of the atomic wave unction ΔU(r) is plotted.

For this function the Bloch theorem holds:

ψnk(r + R) = ReikR Φn(r + R - R) = eikR Reik(R-R)Φ n(r - (R- R)) = eikR R̄eikR̄Φ n(r -R̄) = eikRΦ nk(r) (4.34)

HINT: Eq. (4.34) is the very general rule for the Bloch condition using localized states.

To include the perturbation ΔU(r) of the atomic energy levels we use in Eq. (4.33) the linear combination of atomic orbitals

Φ(r) = nbnψn(r) (4.35)

Thus we replace the purely atomic states with mixed atomic states. The calculation of these coefficients bn is the main task in the following.


With frame Back Forward as PDF

© J. Carstensen (Quantum Mech.)