World Library  
Flag as Inappropriate
Email this Article

Kronig-Penney model

Article Id: WHEBN0002413252
Reproduction Date:

Title: Kronig-Penney model  
Author: World Heritage Encyclopedia
Language: English
Subject: Tight binding, Nearly free electron model, Ralph Kronig
Collection:
Publisher: World Heritage Encyclopedia
Publication
Date:
 

Kronig-Penney model

In quantum mechanics, the particle in a one-dimensional lattice is a problem that occurs in the model of a periodic crystal lattice. The potential is caused by ions in the periodic structure of the crystal creating an electromagnetic field so electrons are subject to a regular potential inside the lattice. This is an extension of the free electron model that assumes zero potential inside the lattice.

Problem definition

When talking about solid materials, the discussion is mainly around crystals - periodic lattices. Here we will discuss a 1D lattice of positive ions. Assuming the spacing between two ions is a, the potential in the lattice will look something like this:

The mathematical representation of the potential is a periodic function with a period a. According to Bloch's theorem, the wavefunction solution of the Schrödinger equation when the potential is periodic, can be written as:

\psi (x) = e^{ikx} u(x). \,\!

Where u(x) is a periodic function which satisfies:

u(x+a)=u(x) \,\!

When nearing the edges of the lattice, there are problems with the boundary condition. Therefore, we can represent the ion lattice as a ring following the Born-von Karman boundary conditions. If L is the length of the lattice so that L >> a, then the number of ions in the lattice is so big, that when considering one ion, its surrounding is almost linear, and the wavefunction of the electron is unchanged. So now, instead of two boundary conditions we get one circular boundary condition:

\psi (0)=\psi (L). \,\!

If N is the number of Ions in the lattice, then we have the relation: aN = L. Replacing in the boundary condition and applying Bloch's theorem will result in a quantization for k:

\psi (0) = e^{ik \cdot 0} u(0) = e^{ikL} u(L) = \psi (L) \,\!
u(0) = e^{ikL} u(L)=e^{ikL} u(N a) \rightarrow e^{ikL} = 1 \,\!
\Rightarrow kL = 2\pi n \rightarrow k = {2\pi \over L} n \qquad \left( n=0, \pm 1, \pm 2, ..., \pm {N \over 2} \right). \,\!

Kronig–Penney model

The Kronig–Penney model (named after Ralph Kronig and William Penney) is a simple, idealized quantum-mechanical system that consists of an infinite periodic array of rectangular potential barriers.

The potential function is approximated by a rectangular potential:


Using Bloch's theorem, we only need to find a solution for a single period, make sure it is continuous and smooth, and to make sure the function u(x) is also continuous and smooth.

Considering a single period of the potential:
We have two regions here. We will solve for each independently:


\mathrm{For} \quad 0 < x <(a-b) :
{-\hbar^2 \over 2m} \psi_{xx} = E \psi \,\!
\Rightarrow \psi = A e^{i \alpha x} + A' e^{-i \alpha x} \quad \left( \alpha^2 = {2mE \over \hbar^2} \right) \,\!


\mathrm{For} \quad -b  :
{-\hbar^2 \over 2m} \psi_{xx} = (E+V_0)\psi \,\!
\Rightarrow \psi = B e^{i \beta x} + B' e^{-i \beta x} \quad \left( \beta^2 = {2m(E+V_0) \over \hbar^2} \right). \,\!

To find u(x) in each region, we need to manipulate the electron's wavefunction:

\psi(0
\Rightarrow u(0

And in the same manner:

u(-b

To complete the solution we need to make sure the probability function is continuous and smooth, i.e.:

\psi(0^{-})=\psi(0^{+}) \qquad \psi'(0^{-})=\psi'(0^{+}). \,\!

And that u(x) and u'(x) are periodic

u(-b)=u(a-b) \qquad u'(-b)=u'(a-b). \,\!

These conditions yield the following matrix:

\begin{pmatrix} 1 & 1 & -1 & -1 \\ \alpha & -\alpha & -\beta & \beta \\ e^{i(\alpha-k)(a-b)} & e^{-i(\alpha+k)(a-b)} & -e^{-i(\beta-k)b} & -e^{i(\beta+k)b} \\ (\alpha-k)e^{i(\alpha-k)(a-b)} & -(\alpha+k)e^{-i(\alpha+k)(a-b)} & -(\beta-k)e^{-i(\beta-k)b} & (\beta+k)e^{i(\beta+k)b} \end{pmatrix} \begin{pmatrix} A \\ A' \\ B \\ B' \end{pmatrix} = \begin{pmatrix} 0 \\ 0 \\ 0 \\ 0 \end{pmatrix}. \,\!

For us not to have the trivial solution, the determinant of the matrix must be 0. This leads us to the following expression:

\cos(k a) = \cos(\beta b) \cos[\alpha(a-b)]-{\alpha^2+\beta^2 \over 2\alpha \beta} \sin(\beta b) \sin[\alpha(a-b)]. \,\!

To further simplify the expression, we perform the following approximations:

b \rightarrow 0 \ ; \ V_0 \rightarrow \infty \ ; \ V_0 b = \mathrm{constant} \,\!
\Rightarrow \beta^2 b = \mathrm{constant} \ ; \ \alpha^2 b \rightarrow 0 \,\!
\Rightarrow \beta b \rightarrow 0 \ ; \ \sin(\beta b) \rightarrow \beta b \ ; \ \cos(\beta b) \rightarrow 1. \,\!

The expression will now be:

\cos(k a) = \cos(\alpha a)-P{\sin(\alpha a) \over \alpha a} \qquad \left( P={m V_0 ba\over \hbar^2} \right). \,\!

Kronig–Penney model: Alternative Solution

An alternative treatment to a similar problem is given. Here we have a delta periodic potential: V(x) = A\cdot\sum_{n=-\infty}^{\infty}\delta(x-n\cdot a). A is some constant, and a is the lattice constant (the spacing between each site). Since this potential is periodic, we could expand it as a Fourier series: V(x) = \sum_K \tilde{V}(K)\cdot e^{i\cdot K \cdot x}, where \tilde{V}(K) = \frac{1}{a}\int_{-a/2}^{a/2}dx\,V(x)\,e^{-i\cdot K\cdot x} = \frac{1}{a}\int_{-a/2}^{a/2}dx\sum_{n=-\infty}^{\infty}A\cdot \delta(x-na)\,e^{-i\,K\,x} = \frac{A}{a}.

The wave-function, using Bloch's theorem, is equal to \psi_k(x) = e^{i k x} u_k(x) where u_k(x) is a function that is periodic in the lattice, which means that we can expand it as a Fourier series as well: u_k(x)=\sum_{K} \tilde{u}_k(K)e^{i K x}. Thus the wave function is: \psi_k(x)=\sum_{K}\tilde{u}_k(K)\,e^{i(k+K)x}.

Putting this into the Schroedinger equation, we get:

\left[\frac{\hbar^2(k+K)^2}{2m}-E_k\right]\cdot\tilde{u}_k(K)+\sum_{K'}\tilde{V}(K-K')\,\tilde{u}_k(K')=0

or rather:

\left[\frac{\hbar^2(k+K)^2}{2m}-E_k\right]\cdot\tilde{u}_k(K)+\frac{A}{a}\sum_{K'}\tilde{u}_k(K')=0

Now we define a new function:

f(k):=\sum_{K'}\tilde{u}_k(K')

Plug this into the Schroedinger equation:

\left[\frac{\hbar^2(k+K)^2}{2m}-E_k\right]\cdot\tilde{u}_k(K)+\frac{A}{a}f(k)=0

Solving this for \tilde{u}_k(K) we get:

\tilde{u}_k(K)=\frac{\frac{2m}{\hbar^2}\frac{A}{a}f(k)}{\frac{2mE_k}{\hbar^2}-(k+K)^2}=\frac{\frac{2m}{\hbar^2}\frac{A}{a}}{\frac{2mE_k}{\hbar^2}-(k+K)^2}\,f(k)

We sum this last equation over all values of K to arrive at:

\sum_{K}\tilde{u}_k(K)=\sum_{K}\frac{\frac{2m}{\hbar^2}\frac{A}{a}}{\frac{2mE_k}{\hbar^2}-(k+K)^2}\,f(k)

Or:

f(k)=\sum_{K}\frac{\frac{2m}{\hbar^2}\frac{A}{a}}{\frac{2mE_k}{\hbar^2}-(k+K)^2}\,f(k)

Conveniently, f(k) cancel outs and we get:

1=\sum_{K}\frac{\frac{2m}{\hbar^2}\frac{A}{a}}{\frac{2mE_k}{\hbar^2}-(k+K)^2}

Or:

\frac{\hbar^2}{2m}\frac{a}{A}=\sum_{K}\frac{1}{\frac{2mE_k}{\hbar^2}-(k+K)^2}

To save ourselves some unnecessary notational effort we define a new variable:

\alpha^2:=\frac{2mE_k}{\hbar^2}

and finally our expression is:

\frac{\hbar^2}{2m}\frac{a}{A}=\sum_{K}\frac{1}{\alpha^2-(k+K)^2}

Now, K is a reciprocal lattice vector, which means that a sum over K is actually a sum over integer multiples of \frac{2\pi}{a}:

\frac{\hbar^2}{2m}\frac{a}{A}=\sum_{n=-\infty}^{\infty}\frac{1}{\alpha^2-(k+\frac{2\pi n}{a})^2}

We can juggle this expression a little bit to make it more suggestive:

\frac{\hbar^2}{2m}\frac{a}{A}=\sum_{n=-\infty}^{\infty}\frac{1}{\alpha^2-(k+\frac{2\pi n}{a})^2}=-\frac{1}{2\alpha}\sum_{n=-\infty}^{\infty}\left[\frac{1}{(k+\frac{2\pi n}{a})-\alpha}-\frac{1}{(k+\frac{2\pi n}{a})+\alpha}\right]=-\frac{a}{4\alpha}\sum_{n=-\infty}^{\infty}\left[\frac{1}{\pi n + \frac{k a}{2}-\frac{\alpha a}{2}}-\frac{1}{\pi n +\frac{k a}{2}+\frac{\alpha a}{2}}\right]=
=-\frac{a}{4\alpha}\left[\left(\sum_{n=-\infty}^{\infty}\frac{1}{\pi n + \frac{k a}{2}-\frac{\alpha a}{2}}\right)-\left(\sum_{n=-\infty}^{\infty}\frac{1}{\pi n +\frac{k a}{2}+\frac{\alpha a}{2}}\right)\right]

If we use a nice identity of a sum of the cotangent function (Equation 18) which says:

cot(x)=\sum_{n=-\infty}^{\infty}\frac{1}{n\pi+x}

and plug it into our expression we get to:

\frac{\hbar^2}{2m}\frac{a}{A}=-\frac{a}{4\alpha}\left[\cot\left(\frac{k a}{2}-\frac{\alpha a}{2}\right)-\cot\left(\frac{k a}{2}+\frac{\alpha a}{2}\right)\right]

We use the sum of \cot and then, the product of \sin (which is part of the formula for the sum of \cot) to arrive at:

\cos(k a)=\cos(\alpha a)+\frac{m A}{\hbar^2 \alpha a}\sin(\alpha a)

This equation shows the relation between the energy (through \alpha) and the wave-vector, k, and as you can see, since the left hand side of the equation can only range from -1 to 1 then there are some limits on the values that \alpha(and thus, the energy) can take, that is, at some ranges of values of the energy, there is no solution according to these equation, and thus, the system will not have those energies: energy gaps. These are the so called band-gaps, which can be shown to exist in any shape of periodic potential (not just delta or square barriers).

See also

External links

  • 1-D periodic potential applet
  • "The Wolfram Demonstrations Project.
This article was sourced from Creative Commons Attribution-ShareAlike License; additional terms may apply. World Heritage Encyclopedia content is assembled from numerous content providers, Open Access Publishing, and in compliance with The Fair Access to Science and Technology Research Act (FASTR), Wikimedia Foundation, Inc., Public Library of Science, The Encyclopedia of Life, Open Book Publishers (OBP), PubMed, U.S. National Library of Medicine, National Center for Biotechnology Information, U.S. National Library of Medicine, National Institutes of Health (NIH), U.S. Department of Health & Human Services, and USA.gov, which sources content from all federal, state, local, tribal, and territorial government publication portals (.gov, .mil, .edu). Funding for USA.gov and content contributors is made possible from the U.S. Congress, E-Government Act of 2002.
 
Crowd sourced content that is contributed to World Heritage Encyclopedia is peer reviewed and edited by our editorial staff to ensure quality scholarly research articles.
 
By using this site, you agree to the Terms of Use and Privacy Policy. World Heritage Encyclopedia™ is a registered trademark of the World Public Library Association, a non-profit organization.
 


Copyright © World Library Foundation. All rights reserved. eBooks from Project Gutenberg are sponsored by the World Library Foundation,
a 501c(4) Member's Support Non-Profit Organization, and is NOT affiliated with any governmental agency or department.