In solidstate physics, the tightbinding model (or TB model) is an approach to the calculation of electronic band structure using an approximate set of wave functions based upon superposition of wave functions for isolated atoms located at each atomic site. The method is closely related to the LCAO method used in chemistry. Tightbinding models are applied to a wide variety of solids. The model gives good qualitative results in many cases and can be combined with other models that give better results where the tightbinding model fails. Though the tightbinding model is a oneelectron model, the model also provides a basis for more advanced calculations like the calculation of surface states and application to various kinds of manybody problem and quasiparticle calculations.
Contents

Introduction 1

Historical background 2

Mathematical formulation 3

Translational symmetry and normalization 3.1

The tight binding Hamiltonian 3.2

The tight binding matrix elements 3.3

Evaluation of the matrix elements 4

Connection to Wannier functions 5

Second quantization 6

Example: onedimensional sband 7

Table of interatomic matrix elements 8

See also 9

References 10

Further reading 11

External links 12
Introduction
The name "tight binding" of this electronic band structure model suggests that this quantum mechanical model describes the properties of tightly bound electrons in solids. The electrons in this model should be tightly bound to the atom to which they belong and they should have limited interaction with states and potentials on surrounding atoms of the solid. As a result the wave function of the electron will be rather similar to the atomic orbital of the free atom to which it belongs. The energy of the electron will also be rather close to the ionization energy of the electron in the free atom or ion because the interaction with potentials and states on neighboring atoms is limited.
Though the mathematical formulation^{[1]} of the oneparticle tightbinding Hamiltonian may look complicated at first glance, the model is not complicated at all and can be understood intuitively quite easily. There are only three kinds of elements that play a significant role in the theory. Two of those three kinds of elements should be close to zero and can often be neglected. The most important elements in the model are the interatomic matrix elements, which would simply be called the bond energies by a chemist.
In general there are a number of atomic energy levels and atomic orbitals involved in the model. This can lead to complicated band structures because the orbitals belong to different pointgroup representations. The reciprocal lattice and the Brillouin zone often belong to a different space group than the crystal of the solid. Highsymmetry points in the Brillouin zone belong to different pointgroup representations. When simple systems like the lattices of elements or simple compounds are studied it is often not very difficult to calculate eigenstates in highsymmetry points analytically. So the tightbinding model can provide nice examples for those who want to learn more about group theory.
The tightbinding model has a long history and has been applied in many ways and with many different purposes and different outcomes. The model doesn't stand on its own. Parts of the model can be filled in or extended by other kinds of calculations and models like the molecular electronics, for example, tightbindinglike models are applied in which the role of the atoms in the original concept is replaced by the molecular orbitals of conjugated systems and where the interatomic matrix elements are replaced by inter or intramolecular hopping and tunneling parameters. These conductors nearly all have very anisotropic properties and sometimes are almost perfectly onedimensional.
Historical background
By 1928, the idea of a molecular orbital had been advanced by

Crystalfield Theory, Tightbinding Method, and JahnTeller Effect in E. Pavarini, E. Koch, F. Anders, and M. Jarrell (eds.): Correlated Electrons: From Models to Materials, Jülich 2012, ISBN 9783893367962
External links
Further reading

N. W. Ashcroft and N. D. Mermin, Solid State Physics (Thomson Learning, Toronto, 1976).

Stephen Blundell Magnetism in Condensed Matter(Oxford, 2001).

S.Maekawa et al. Physics of Transition Metal Oxides (SpingerVerlag Berlin Heidelberg, 2004).

John Singleton Band Theory and Electronic Properties of Solids (Oxford, 2001).

^ ^{a} ^{b} ^{c}

^ ^{a} ^{b}

^ As an alternative to neglecting overlap, one may choose as a basis instead of atomic orbitals a set of orbitals based upon atomic orbitals but arranged to be orthogonal to orbitals on other atomic sites, the socalled Löwdin orbitals. See

^ Orfried Madelung, Introduction to SolidState Theory (SpringerVerlag, Berlin Heidelberg, 1978).

^

^
References
See also
Not all interatomic matrix elements are listed explicitly. Matrix elements that are not listed in this table can be constructed by permutation of indices and cosine directions of other matrix elements in the table.

E_{s,s} = V_{ss\sigma}

E_{s,x} = l V_{sp\sigma}

E_{x,x} = l^2 V_{pp\sigma} + (1  l^2) V_{pp\pi}

E_{x,y} = l m V_{pp\sigma}  l m V_{pp\pi}

E_{x,z} = l n V_{pp\sigma}  l n V_{pp\pi}

E_{s,xy} = \sqrt{3} l m V_{sd\sigma}

E_{s,x^2y^2} = \frac{\sqrt{3}}{2} (l^2  m^2) V_{sd\sigma}

E_{s,3z^2r^2} = [n^2  (l^2 + m^2) / 2] V_{sd\sigma}

E_{x,xy} = \sqrt{3} l^2 m V_{pd\sigma} + m (1  2 l^2) V_{pd\pi}

E_{x,yz} = \sqrt{3} l m n V_{pd\sigma}  2 l m n V_{pd\pi}

E_{x,zx} = \sqrt{3} l^2 n V_{pd\sigma} + n (1  2 l^2) V_{pd\pi}

E_{x,x^2y^2} = \frac{\sqrt{3}}{2} l (l^2  m^2) V_{pd\sigma} + l (1  l^2 + m^2) V_{pd\pi}

E_{y,x^2y^2} = \frac{\sqrt{3}}{2} m(l^2  m^2) V_{pd\sigma}  m (1 + l^2  m ^2) V_{pd\pi}

E_{z,x^2y^2} = \frac{\sqrt{3}}{2} n(l^2  m^2) V_{pd\sigma}  n(l^2  m^2) V_{pd\pi}

E_{x,3z^2r^2} = l[n^2  (l^2 + m^2)/2]V_{pd\sigma}  \sqrt{3} l n^2 V_{pd\pi}

E_{y,3z^2r^2} = m [n^2  (l^2 + m^2) / 2] V_{pd\sigma}  \sqrt{3} m n^2 V_{pd\pi}

E_{z,3z^2r^2} = n [n^2  (l^2 + m^2) / 2] V_{pd\sigma} + \sqrt{3} n (l^2 + m^2) V_{pd\pi}

E_{xy,xy} = 3 l^2 m^2 V_{dd\sigma} + (l^2 + m^2  4 l^2 m^2) V_{dd\pi} + (n^2 + l^2 m^2) V_{dd\delta}

E_{xy,yz} = 3 l m^2 nV_{dd\sigma} + l n (1  4 m^2) V_{dd\pi} + l n (m^2  1) V_{dd\delta}

E_{xy,zx} = 3 l^2 m n V_{dd\sigma} + m n (1  4 l^2) V_{dd\pi} + m n (l^2  1) V_{dd\delta}

E_{xy,x^2y^2} = \frac{3}{2} l m (l^2  m^2) V_{dd\sigma} + 2 l m (m^2  l^2) V_{dd\pi} + [l m (l^2  m^2) / 2] V_{dd\delta}

E_{yz,x^2y^2} = \frac{3}{2} m n (l^2  m^2) V_{dd\sigma}  m n [1 + 2(l^2  m^2)] V_{dd\pi} + m n [1 + (l^2  m^2) / 2] V_{dd\delta}

E_{zx,x^2y^2} = \frac{3}{2} n l (l^2  m^2) V_{dd\sigma} + n l [1  2(l^2  m^2)] V_{dd\pi}  n l [1  (l^2  m^2) / 2] V_{dd\delta}

E_{xy,3z^2r^2} = \sqrt{3} \left[ l m (n^2  (l^2 + m^2) / 2) ]V_{dd\sigma}  2 l m n^2 V_{dd\pi} + [l m (1 + n^2) / 2] V_{dd\delta} \right]

E_{yz,3z^2r^2} = \sqrt{3} \left[ m n (n^2  (l^2 + m^2) / 2) V_{dd\sigma} + m n (l^2 + m^2  n^2) V_{dd\pi} [ m n (l^2 + m^2) / 2 ]V_{dd\delta} \right]

E_{zx,3z^2r^2} = \sqrt{3} \left[ l n (n^2  (l^2 + m^2) / 2) V_{dd\sigma} + l n (l^2 + m^2  n^2) V_{dd\pi}  [l n (l^2 + m^2) / 2] V_{dd\delta} \right]

E_{x^2y^2,x^2y^2} = \frac{3}{4} (l^2  m^2)^2 V_{dd\sigma} + [l^2 + m^2  (l^2  m^2)^2] V_{dd\pi} + [n^2 + (l^2  m^2)^2 / 4] V_{dd\delta}

E_{x^2y^2,3z^2r^2} = \sqrt{3} \left[ (l^2  m^2) [n^2  (l^2 + m^2) / 2] V_{dd\sigma} / 2 + n^2 (m^2  l^2) V_{dd\pi} + [(1 + n^2)(l^2  m^2) / 4 ]V_{dd\delta}\right]

E_{3z^2r^2,3z^2r^2} = [n^2  (l^2 + m^2) / 2]^2 V_{dd\sigma} + 3 n^2 (l^2 + m^2) V_{dd\pi} + \frac{3}{4} (l^2 + m^2)^2 V_{dd\delta}
where d is the distance between the atoms and l, m and n are the direction cosines to the neighboring atom.

\vec{\bold{r}}_{n,n'} = (r_x,r_y,r_z) = d (l,m,n)
The interatomic vector is expressed as
which, with a little patience and effort, can also be derived from the cubic harmonic orbitals straightforwardly. The table expresses the matrix elements as functions of LCAO twocentre bond integrals between two cubic harmonic orbitals, i and j, on adjacent atoms. The bond integrals are for example the V_{ss\sigma}, V_{pp\pi} and V_{dd\delta} for sigma, pi and delta bonds.

E_{i,j}(\vec{\bold{r}}_{n,n'}) = \langle n,iHn',j\rangle
In 1954 J.C. Slater and F.G. Koster published, mainly for the calculation of transition metal dbands, a table of interatomic matrix elements^{[1]}
Table of interatomic matrix elements
This example is readily extended to three dimensions, for example, to a bodycentered cubic or facecentered cubic lattice by introducing the nearest neighbor vector locations in place of simply n a.^{[6]} Likewise, the method can be extended to multiple bands using multiple different atomic orbitals at each site. The general formulation above shows how these extensions can be accomplished.

For k = 0 the energy is E = (E_0  2 \Delta)/ (1 + 2 S) and the state consists of a sum of all atomic orbitals. This state can be viewed as a chain of bonding orbitals.

For k = \pi / (2 a) the energy is E = E_0 and the state consists of a sum of atomic orbitals which are a factor e^{i \pi / 2} out of phase. This state can be viewed as a chain of nonbonding orbitals.

Finally for k = \pi / a the energy is E = (E_0 + 2 \Delta) / (1  2 S) and the state consists of an alternating sum of atomic orbitals. This state can be viewed as a chain of antibonding orbitals.

E(k)= \frac{E_02\Delta\,\cos(ka)}{1 + 2 S\,\cos(ka)}.
Thus the energy of this state k\rangle can be represented in the familiar form of the energy dispersion:

\frac{1}{N}\sum_n \langle n1Hn\rangle e^{+ika}=\Delta e^{ika}\frac{1}{N}\sum_n 1 = \Delta e^{ika} \ .

\frac{1}{N}\sum_n \langle n1n\rangle e^{+ika}= S e^{ika}\frac{1}{N}\sum_n 1 = S e^{ika} \ .
and

\frac{1}{N}\sum_n \langle nHn\rangle = E_0 \frac{1}{N}\sum_n 1 = E_0 \ ,
where, for example,

Hk\rangle=\frac{1}{\sqrt{N}}\sum_n e^{inka} H n\rangle

\langle k Hk\rangle =\frac{1}{N}\sum_{n,\ m} e^{i(nm)ka} \langle mHn\rangle =\frac{1}{N}\sum_n \langle nHn\rangle+\frac{1}{N}\sum_n \langle n1Hn\rangle e^{+ika}+\frac{1}{N}\sum_n\langle n+1Hn\rangle e^{ika} = E_0 2\Delta\,\cos(ka)\ ,
The energy E_{i} is the ionization energy corresponding to the chosen atomic orbital and U is the energy shift of the orbital as a result of the potential of neighboring atoms. The \langle n\pm 1Hn\rangle=\Delta elements, which are the Slater and Koster interatomic matrix elements, are the bond energies E_{i,j}. In this one dimensional sband model we only have \sigmabonds between the sorbitals with bond energy E_{s,s} = V_{ss\sigma}. The overlap between states on neighboring atoms is S. We can derive the energy of the state k\rangle using the above equation:

\langle nn\rangle= 1 \ ; \langle n \pm 1n\rangle= S \ .

\langle n\pm 1Hn\rangle=\Delta \

\langle nHn\rangle= E_0 = E_i  U \ .
where N = total number of sites and k is a real parameter with \frac{\pi}{a}\leqq k\leqq\frac{\pi}{a}. (This wave function is normalized to unity by the leading factor 1/√N provided overlap of atomic wave functions is ignored.) Assuming only nearest neighbor overlap, the only nonzero matrix elements of the Hamiltonian can be expressed as

k\rangle =\frac{1}{\sqrt{N}}\sum_{n=1}^N e^{inka} n\rangle
To find approximate eigenstates of the Hamiltonian, we can use a linear combination of the atomic orbitals
Here the tight binding model is illustrated with a sband model for a string of atoms with a single sorbital in a straight line with spacing a and σ bonds between atomic sites.
Example: onedimensional sband
This interaction Hamiltonian includes direct Coulomb interaction energy and exchange interaction energy between electrons. There are several novel physics induced from this electronelectron interaction energy, such as metalinsulator transitions (MIT), hightemperature superconductivity, and several quantum phase transitions.

\displaystyle H_{ee}=\frac{1}{2}\sum_{n,m,\sigma}\langle n_1 m_1, n_2 m_2\frac{e^2}{r_1r_2}n_3 m_3, n_4 m_4\rangle c^\dagger_{n_1 m_1 \sigma_1}c^\dagger_{n_2 m_2 \sigma_2}c_{n_4 m_4 \sigma_2} c_{n_3 m_3 \sigma_1}
In the strongly correlated electron system, it is necessary to consider the electronelectron interaction. This term can be written in
Here, hopping integral \displaystyle t corresponds to the transfer integral \displaystyle\gamma in tight binding model. Considering extreme cases of t\rightarrow 0, it is impossible for electron to hop into neighboring sites. This case is the isolated atomic system. If the hopping term is turned on (\displaystyle t>0) electrons can stay in both sites lowering their kinetic energy.

\displaystyle \langle i,j \rangle nearest neighbor index

\displaystyle t  hopping integral

\displaystyle\sigma  spin polarization

H = t \sum_{\langle i,j \rangle,\sigma}( c^{\dagger}_{i,\sigma} c^{}_{j,\sigma}+ h.c.),

c^\dagger_{i\sigma} , c_{j\sigma}  creation and annihilation operators
Using the atomic orbital as a basis state, the second quantization Hamiltonian operator in the tight binding framework can be written as:
Modern explanations of electronic structure like tJ model and Hubbard model are based on tight binding model.^{[5]} Tight binding can be understood by working under a second quantization formalism.
Second quantization
However it is not easy to calculate directly either Bloch functions or Wannier functions. An approximate approach is necessary in the calculation of electronic structures of solids. If we consider the extreme case of isolated atoms, the Wannier function would become an isolated atomic orbital. That limit suggests the choice of an atomic wave function as an approximate form for the Wannier function, the socalled tight binding approximation.
These real space wave functions {a_m\mathbf{(R_n,r)}} are called Wannier functions, and are fairly closely localized to the atomic site R_{n}. Of course, if we have exact Wannier functions, the exact Bloch functions can be derived using the inverse Fourier transform.

a_m\mathbf{(R_n,r)}=\frac{1}{\sqrt{N}}\sum_{\mathbf{k}}{e^{\mathbf{ik\cdot R_n}}\psi_m\mathbf{(k,r)}}=\frac{1}{\sqrt{N}}\sum_{\mathbf{k}}{e^{\mathbf{ik\cdot (rR_n)}}u_m\mathbf{(k,r)}}.
Using the Fourier transform analysis, a spatially localized wave function for the mth energy band can be constructed from multiple Bloch waves:
where R_{n} denotes an atomic site in a periodic crystal lattice, k is the wave vector of the Bloch wave, r is the electron position, m is the band index, and the sum is over all N atomic sites. The Bloch wave is an exact eigensolution for the wave function of an electron in a periodic crystal potential corresponding to an energy E_{m} (k), and is spread over the entire crystal volume.

\psi_m\mathbf{(k,r)}=\frac{1}{\sqrt{N}}\sum_{n}{a_m\mathbf{(R_n,r)}} e^{\mathbf{ik\cdot R_n}}\ ,
Bloch wave functions describe the electronic states in a periodic crystal lattice. Bloch functions can be represented as a Fourier series^{[4]}
Connection to Wannier functions
The tight binding model works particularly well in cases where the band width is small and the electrons are strongly localized, like in the case of dbands and fbands. The model also gives good results in the case of open crystal structures, like diamond or silicon, where the number of neighbors is small. The model can easily be combined with a nearly free electron model in a hybrid NFETB model.^{[2]}
The inter atomic overlap matrix elements \alpha_{m,l} should be rather small or neglectable. If they are large it is again an indication that the tight binding model is of limited value for some purposes. Large overlap is an indication for too short inter atomic distance for example. In metals and transition metals the broad sband or spband can be fitted better to an existing band structure calculation by the introduction of nextnearestneighbor matrix elements and overlap integrals but fits like that don't yield a very useful model for the electronic wave function of a metal. Broad bands in dense materials are better described by a nearly free electron model.
The inter atomic matrix elements \gamma_{m,l} can be calculated directly if the atomic wave functions and the potentials are known in detail. Most often this is not the case. There are numerous ways to get parameters for these matrix elements. Parameters can be obtained from chemical bond energy data. Energies and eigenstates on some high symmetry points in the Brillouin zone can be evaluated and values integrals in the matrix elements can be matched with band structure data from other sources.
As mentioned before the values of the \beta_mmatrix elements are not so large in comparison with the ionization energy because the potentials of neighboring atoms on the central atom are limited. If \beta_m is not relatively small it means that the potential of the neighboring atom on the central atom is not small either. In that case it is an indication that the tight binding model is not a very good model for the description of the band structure for some reason. The inter atomic distances can be too small or the charges on the atoms or ions in the lattice is wrong for example.
Evaluation of the matrix elements
denote the overlap integrals between the atomic orbitals m and l on adjacent atoms.

\alpha_{m,l}(\boldsymbol{R_n}) = \int \varphi_m^*(\boldsymbol{r}) \varphi_l(\boldsymbol{r  R_n}) \, d^3 r \ ,
The last terms
is the inter atomic matrix element between the atomic orbitals m and l on adjacent atoms. It is also called the bond energy or two center integral and it is the most important element in the tight binding model.

\gamma_{m,l}(\boldsymbol{R_n}) = \int \varphi_m^*(\boldsymbol{r}) \Delta U(\boldsymbol{r}) \varphi_l(\boldsymbol{r  R_n}) \, d^3 r \ ,
The next term
is the atomic energy shift due to the potential on neighboring atoms. This term is relatively small in most cases. If it is large it means that potentials on neighboring atoms have a large influence on the energy of the central atom.

\beta_m = \int \varphi_m^*(\boldsymbol{r})\Delta U(\boldsymbol{r}) \varphi_m(\boldsymbol{r}) \, d^3 r \ ,
The element
The tight binding matrix elements
where E_{m} is the energy of the mth atomic level, and \alpha_{m,l}, \beta_m and \gamma_{m,l} are the tight binding matrix elements.

\varepsilon_m(\boldsymbol{k}) = E_m  N\ b (0)^2 \left(\beta_m + \sum_{\boldsymbol{R_n}\neq 0}\sum_l \gamma_{m,l}(\boldsymbol{R_n}) e^{i \boldsymbol{k} \cdot \boldsymbol{R_n}}\right) \ ,


= E_m  \ \frac {\beta_m + \sum_{\boldsymbol{R_n}\neq 0}\sum_l e^{i \boldsymbol{k} \cdot \boldsymbol{R_n}} \gamma_{m,l}(\boldsymbol{R_n})}{\ \ 1 + \sum_{\boldsymbol{R_n \neq 0}}\sum_l e^{i \boldsymbol{k \cdot R_n}} \alpha_{m,l} (\boldsymbol{R_n})} \ ,
Here terms involving the atomic Hamiltonian at sites other than where it is centered are neglected. The energy then becomes



\approx E_m + b^*(0)\sum_{\boldsymbol{R_n}} e^{i \boldsymbol{k \cdot R_n}}\ \int d^3 r \ \varphi^* (\boldsymbol{rR_n})\Delta U (\boldsymbol{r}) \psi (\boldsymbol{r}) \ .

\varepsilon_m = \int d^3 r \ \psi^* (\boldsymbol{r})H(\boldsymbol{r}) \psi (\boldsymbol{r})


=\sum_{\boldsymbol{R_n}} b^* ( \boldsymbol{R_n})\ \int d^3 r \ \varphi^* (\boldsymbol{rR_n})H(\boldsymbol{r}) \psi (\boldsymbol{r}) \

=\sum_{\boldsymbol{R_{\ell}}} \ \sum_{\boldsymbol{R_n}} b^* ( \boldsymbol{R_n})\ \int d^3 r \ \varphi^* (\boldsymbol{rR_n})H_{\mathrm{at}}(\boldsymbol{rR_{\ell}}) \psi (\boldsymbol{r}) \ + \sum_{\boldsymbol{R_n}} b^*( \boldsymbol{R_n})\ \int d^3 r \ \varphi^* (\boldsymbol{rR_n})\Delta U (\boldsymbol{r}) \psi (\boldsymbol{r}) \ .
Using the tight binding form for the wave function, and assuming only the mth atomic energy level is important for the mth energy band, the Bloch energies \varepsilon_m are of the form
The tight binding Hamiltonian


\psi (\boldsymbol{r}) \approx \frac {1} {\sqrt{N}} \sum_{m,\boldsymbol{R_n}} e^{i \boldsymbol{k \cdot R_n}} \ \varphi_m (\boldsymbol{rR_n}) \ .
and

b_n (0) \approx \frac {1} {\sqrt{N}} \ ,
where α (R_{p} ) are the atomic overlap integrals, which frequently are neglected resulting in^{[3]}

b^*(0)b(0) = \frac {1} {N}\ \cdot \ \frac {1}{1 + \sum_{\boldsymbol{R_p \neq 0}} e^{i \boldsymbol{k \cdot R_p}} \alpha (\boldsymbol{R_p})} \ ,
so the normalization sets b(0) as



=N b^*(0)b(0)\sum_{\boldsymbol{R_p}} e^{i \boldsymbol{k \cdot R_p}}\ \int d^3 r \ \varphi^* (\boldsymbol{r}) \varphi (\boldsymbol{rR_p})\ ,



=N b^*(0)b(0)\sum_{\boldsymbol{R_p}} e^{i \boldsymbol{k \cdot R_p}}\ \int d^3 r \ \varphi^* (\boldsymbol{rR_p}) \varphi (\boldsymbol{r})\



= b^*(0)b(0)\sum_{\boldsymbol{R_n}} e^{i \boldsymbol{k \cdot R_n}}\sum_{\boldsymbol{R_{\ell}}} e^ {i \boldsymbol{k \cdot R_{\ell}}}\ \int d^3 r \ \varphi^* (\boldsymbol{rR_n}) \varphi (\boldsymbol{rR_{\ell}})



= \sum_{\boldsymbol{R_n}} b^* ( \boldsymbol{R_n})\sum_{\boldsymbol{R_{\ell}}} b ( \boldsymbol{R_{\ell}})\int d^3 r \ \varphi^* (\boldsymbol{rR_n}) \varphi (\boldsymbol{rR_{\ell}})

\int d^3 r \ \psi^* (\boldsymbol{r}) \psi (\boldsymbol{r}) = 1
Normalizing the wave function to unity:

b_m (\boldsymbol{R_p}) = e^{i\boldsymbol{k \cdot R_{p}}} b_m ( \boldsymbol{0}) \ .
or

b_m ( \boldsymbol{R_p+R_{\ell}}) = e^{i\boldsymbol{k \cdot R_{\ell}}}b_m ( \boldsymbol{R_p}) \ , (where in RHS we have replaced the dummy index \boldsymbol{R_n} with \boldsymbol{R_p} )
By substituting \boldsymbol{R_p}= \boldsymbol{R_n}  \boldsymbol{R_\ell}, we find

\sum_{m,\boldsymbol{R_n}} b_m ( \boldsymbol{R_n}) \ \varphi_m (\boldsymbol{rR_n+R_{\ell}})=e^{i\boldsymbol{k \cdot R_{\ell}}}\sum_{m,\boldsymbol{R_n}} b_m ( \boldsymbol{R_n}) \ \varphi_m (\boldsymbol{rR_n})\ .
where \bold{k} is the wave vector of the wave function. Consequently, the coefficients satisfy

\psi(\boldsymbol{r+R_{\ell}}) = e^{i\boldsymbol{k \cdot R_{\ell}}}\psi(\boldsymbol{r}) \ ,
The Bloch theorem states that the wave function in crystal can change under translation only by a phase factor:
Translational symmetry and normalization
where m refers to the mth atomic energy level and R_n locates an atomic site in the crystal lattice.

\psi(\boldsymbol{r}) = \sum_{m,\boldsymbol{R_n}} b_m ( \boldsymbol{R_n}) \ \varphi_m (\boldsymbol{rR_n}),
A solution \psi_r to the timeindependent single electron Schrödinger equation is then approximated as a linear combination of atomic orbitals \varphi_m( \boldsymbol{r R_n}):

H (\boldsymbol{r}) = \sum_{\boldsymbol{R_n}} H_{\mathrm{at}}(\boldsymbol{r  R_n}) +\Delta U (\boldsymbol{r}) \ .
We introduce the atomic orbitals \varphi_m( \boldsymbol{r} ), which are eigenfunctions of the Hamiltonian H_{at} of a single isolated atom. When the atom is placed in a crystal, this atomic wave function overlaps adjacent atomic sites, and so are not true eigenfunctions of the crystal Hamiltonian. The overlap is less when electrons are tightly bound, which is the source of the descriptor "tightbinding". Any corrections to the atomic potential \Delta U required to obtain the true Hamiltonian H of the system, are assumed small:
Mathematical formulation
The tightbinding model is typically used for calculations of electronic band structure and band gaps in the static regime. However, in combination with other methods such as the random phase approximation (RPA) model, the dynamic response of systems may also be studied.
Recently, in the research about strongly correlated material, the tight binding approach is basic approximation because highly localized electrons like 3d transition metal electrons sometimes display strongly correlated behaviors. In this case, the role of electronelectron interaction must be considered using the manybody physics description.
In this approach, interactions between different atomic sites are considered as perturbations. There exist several kinds of interactions we must consider. The crystal Hamiltonian is only approximately a sum of atomic Hamiltonians located at different sites and atomic wave functions overlap adjacent atomic sites in the crystal, and so are not accurate representations of the exact wave function. There are further explanations in the next section with some mathematical expressions.
between these points.
Brillouin zone but, rather, firstprinciples calculations are carried out only at highsymmetry points and the band structure is interpolated over the remainder of the Bloch's theorem. With the SK tightbinding method, electronic band structure calculations on a solid need not be carried out with full rigor as in the original SK tightbinding method sometimes referred to as the [1]
This article was sourced from Creative Commons AttributionShareAlike 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, EGovernment 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 nonprofit organization.