Wannier function

Wannier functions of triple- and single-bonded nitrogen dimers in palladium nitride.

The Wannier functions are a complete set of orthogonal functions used in solid-state physics. They were introduced by Gregory Wannier in 1937.[1][2] Wannier functions are the localized molecular orbitals of crystalline systems.

The Wannier functions for different lattice sites in a crystal are orthogonal, allowing a convenient basis for the expansion of electron states in certain regimes. Wannier functions have found widespread use, for example, in the analysis of binding forces acting on electrons.

Definition

Example of a localized Wannier function of titanium in barium titanate (BaTiO3)

Although, like localized molecular orbitals, Wannier functions can be chosen in many different ways,[3] the original,[1] simplest, and most common definition in solid-state physics is as follows. Choose a single band in a perfect crystal, and denote its Bloch states by

${\displaystyle \psi _{\mathbf {k} }(\mathbf {r} )=e^{i\mathbf {k} \cdot \mathbf {r} }u_{\mathbf {k} }(\mathbf {r} )}$

where u_k(r) has the same periodicity as the crystal. Then the Wannier functions are defined by

${\displaystyle \phi _{\mathbf {R} }(\mathbf {r} )={\frac {1}{\sqrt {N}}}\sum _{\mathbf {k} }e^{-i\mathbf {k} \cdot \mathbf {R} }\psi _{\mathbf {k} }(\mathbf {r} )}$,

where

• R is any lattice vector (i.e., there is one Wannier function for each Bravais lattice vector);
• N is the number of primitive cells in the crystal;
• The sum on k includes all the values of k in the Brillouin zone (or any other primitive cell of the reciprocal lattice) that are consistent with periodic boundary conditions on the crystal. This includes N different values of k, spread out uniformly through the Brillouin zone. Since N is usually very large, the sum can be written as an integral according to the replacement rule:

${\displaystyle \sum _{\mathbf {k} }\longrightarrow {\frac {N}{\Omega }}\int _{\text{BZ}}d^{3}\mathbf {k} }$

where "BZ" denotes the Brillouin zone, which has volume Ω.

Properties

On the basis of this definition, the following properties can be proven to hold:[4]

• For any lattice vector R' ,

${\displaystyle \phi _{\mathbf {R} }(\mathbf {r} )=\phi _{\mathbf {R} +\mathbf {R} '}(\mathbf {r} +\mathbf {R} ')}$

In other words, a Wannier function only depends on the quantity (r − R). As a result, these functions are often written in the alternative notation

${\displaystyle \phi (\mathbf {r} -\mathbf {R} ):=\phi _{\mathbf {R} }(\mathbf {r} )}$

• The Bloch functions can be written in terms of Wannier functions as follows:

${\displaystyle \psi _{\mathbf {k} }(\mathbf {r} )={\frac {1}{\sqrt {N}}}\sum _{\mathbf {R} }e^{i\mathbf {k} \cdot \mathbf {R} }\phi _{\mathbf {R} }(\mathbf {r} )}$,

where the sum is over each lattice vector R in the crystal.

• The set of wavefunctions ${\displaystyle \phi _{\mathbf {R} }}$ is an orthonormal basis for the band in question.

${\displaystyle {\begin{aligned}\int _{\text{crystal}}\phi _{\mathbf {R} }(\mathbf {r} )^{*}\phi _{\mathbf {R'} }(\mathbf {r} )d^{3}\mathbf {r} &={\frac {1}{N}}\sum _{\mathbf {k,k'} }\int _{\text{crystal}}e^{i\mathbf {k} \cdot \mathbf {R} }\psi _{\mathbf {k} }(\mathbf {r} )^{*}e^{-i\mathbf {k'} \cdot \mathbf {R'} }\psi _{\mathbf {k'} }(\mathbf {r} )d^{3}\mathbf {r} \\&={\frac {1}{N}}\sum _{\mathbf {k,k'} }e^{i\mathbf {k} \cdot \mathbf {R} }e^{-i\mathbf {k'} \cdot \mathbf {R'} }\delta _{\mathbf {k,k'} }\\&={\frac {1}{N}}\sum _{\mathbf {k} }e^{i\mathbf {k} \cdot \mathbf {(R-R')} }\\&=\delta _{\mathbf {R,R'} }\end{aligned}}}$

Wannier functions have been extended to nearly periodic potentials as well.[5]

Localization

The Bloch states ψ_k(r) are defined as the eigenfunctions of a particular Hamiltonian, and are therefore defined only up to an overall phase. By applying a phase transformation e^iθ(k) to the functions ψ_k(r), for any (real) function θ(k), one arrives at an equally valid choice. While the change has no consequences for the properties of the Bloch states, the corresponding Wannier functions are significantly changed by this transformation.

One therefore uses the freedom to choose the phases of the Bloch states in order to give the most convenient set of Wannier functions. In practice, this is usually the maximally-localized set, in which the Wannier function ϕ_R is localized around the point R and rapidly goes to zero away from R. For the one-dimensional case, it has been proved by Kohn[6] that there is always a unique choice that gives these properties (subject to certain symmetries). This consequently applies to any separable potential in higher dimensions; the general conditions are not established, and are the subject of ongoing research.[7]

A Pipek-Mezey style localization scheme has also been recently proposed for obtaining Wannier functions.[8] Contrary to the maximally localized Wannier functions (which are an application of the Foster-Boys scheme to crystalline systems), the Pipek-Mezey Wannier functions do not mix σ and π orbitals.

Rigorous results

The existence of exponentially localized Wannier functions in insulators was proved mathematically in 2006.[7]

Modern theory of polarization

Wannier functions have recently found application in describing the polarization in crystals, for example, ferroelectrics. The modern theory of polarization is pioneered by Raffaele Resta and David Vanderbilt. See for example, Berghold,[9] and Nakhmanson,[10] and a power-point introduction by Vanderbilt.[11] The polarization per unit cell in a solid can be defined as the dipole moment of the Wannier charge density:

${\displaystyle \mathbf {p_{c}} =-e\sum _{n}\int \ d^{3}r\,\,\mathbf {r} |W_{n}(\mathbf {r} )|^{2}\ ,}$

where the summation is over the occupied bands, and W_n is the Wannier function localized in the cell for band n. The change in polarization during a continuous physical process is the time derivative of the polarization and also can be formulated in terms of the Berry phase of the occupied Bloch states.[4][12]

Wannier interpolation

Wannier functions are often used to interpolate bandstructures calculated ab initio on a coarse grid of k-points to any arbitrary k-point. This is particularly useful for evaluation of Brillouin-zone integrals on dense grids and searching of Weyl points, and also taking derivatives in the k-space. This approach is similar in spirit to the tight binding approximation, but in contrast allows for an exact description of bands in a certain energy range. Wannier interpolation schemes have been derived for spectral properties,[13] anomalous Hall conductivity,[14] orbital magnetization,[15] thermoelectric and electronic transport properties,[16] gyrotropic effects,[17] shift current,[18] spin Hall conductivity [19] [20] and other effects.

See also

• Orbital magnetization

References

1. Wannier Gregory H (1937). "The Structure of Electronic Excitation Levels in Insulating Crystals". Physical Review. 52 (3): 191–197. Bibcode:1937PhRv...52..191W. doi:10.1103/PhysRev.52.191.
2. Wannier, Gregory H. (1 September 1962). "Dynamics of Band Electrons in Electric and Magnetic Fields". Reviews of Modern Physics. 34 (4). American Physical Society (APS): 645–655. Bibcode:1962RvMP...34..645W. doi:10.1103/revmodphys.34.645. ISSN 0034-6861.
3. Marzari et al.: An Introduction to Maximally-Localized Wannier Functions
4. A Bohm, A Mostafazadeh, H Koizumi, Q Niu and J Zqanziger (2003). The Geometric Phase in Quantum Systems. Springer. pp. §12.5, p. 292 ff. doi:10.1007/978-3-662-10333-3. ISBN 978-3-540-00031-0.
5. MP Geller and W Kohn Theory of generalized Wannier functions for nearly periodic potentials Physical Review B 48, 1993
6. W. Kohn (1959). "Analytic Properties of Bloch Waves and Wannier Functions". Physical Review. 115 (4): 809–821. Bibcode:1959PhRv..115..809K. doi:10.1103/PhysRev.115.809.
7. Brouder, Christian; Panati, Gianluca; Calandra, Matteo; Mourougane, Christophe; Marzari, Nicola (25 January 2007). "Exponential Localization of Wannier Functions in Insulators". Physical Review Letters. 98 (4). American Physical Society (APS): 046402. arXiv:cond-mat/0606726. Bibcode:2007PhRvL..98d6402B. doi:10.1103/physrevlett.98.046402. ISSN 0031-9007. PMID 17358792. S2CID 32812449.
8. Jónsson Elvar Ö., Lehtola Susi, Puska Martti, Jónsson Hannes (2017). "Theory and Applications of Generalized Pipek–Mezey Wannier Functions". Journal of Chemical Theory and Computation. 13 (2): 460–474. arXiv:1608.06396. doi:10.1021/acs.jctc.6b00809. PMID 28099002. S2CID 206612913.
9. Berghold, Gerd; Mundy, Christopher J.; Romero, Aldo H.; Hutter, Jürg; Parrinello, Michele (15 April 2000). "General and efficient algorithms for obtaining maximally localized Wannier functions". Physical Review B. 61 (15). American Physical Society (APS): 10040–10048. Bibcode:2000PhRvB..6110040B. doi:10.1103/physrevb.61.10040. ISSN 0163-1829.
10. Nakhmanson, S. M.; Calzolari, A.; Meunier, V.; Bernholc, J.; Buongiorno Nardelli, M. (10 June 2003). "Spontaneous polarization and piezoelectricity in boron nitride nanotubes". Physical Review B. 67 (23): 235406. arXiv:cond-mat/0305329v1. Bibcode:2003PhRvB..67w5406N. doi:10.1103/physrevb.67.235406. ISSN 0163-1829. S2CID 119345964.
11. D Vanderbilt Berry phases and Curvatures in Electronic Structure Theory.
12. C. Pisani (1994). Quantum-mechanical Ab-initio Calculation of the Properties of Crystalline Materials (Proceedings of the IV School of Computational Chemistry of the Italian Chemical Society ed.). Springer. p. 282. ISBN 978-3-540-61645-0.
13. Yates, Jonathan R.; Wang, Xinjie; Vanderbilt, David; Souza, Ivo (2007-05-21). "Spectral and Fermi surface properties from Wannier interpolation". Physical Review B. 75 (19). American Physical Society (APS): 195121. arXiv:cond-mat/0702554. Bibcode:2007PhRvB..75s5121Y. doi:10.1103/physrevb.75.195121. ISSN 1098-0121. S2CID 31224663.
14. Wang, Xinjie; Yates, Jonathan R.; Souza, Ivo; Vanderbilt, David (2006-11-21). "Ab initiocalculation of the anomalous Hall conductivity by Wannier interpolation". Physical Review B. 74 (19): 195118. arXiv:cond-mat/0608257. Bibcode:2006PhRvB..74s5118W. doi:10.1103/physrevb.74.195118. ISSN 1098-0121. S2CID 30427871.
15. Lopez, M. G.; Vanderbilt, David; Thonhauser, T.; Souza, Ivo (2012-01-31). "Wannier-based calculation of the orbital magnetization in crystals". Physical Review B. 85 (1): 014435. arXiv:1112.1938. Bibcode:2012PhRvB..85a4435L. doi:10.1103/physrevb.85.014435. ISSN 1098-0121. S2CID 44056938.
16. Pizzi, Giovanni; Volja, Dmitri; Kozinsky, Boris; Fornari, Marco; Marzari, Nicola (2014-01-01). "BoltzWann: A code for the evaluation of thermoelectric and electronic transport properties with a maximally-localized Wannier functions basis". Computer Physics Communications. 185 (1): 422–429. arXiv:1305.1587. Bibcode:2014CoPhC.185..422P. doi:10.1016/j.cpc.2013.09.015. ISSN 0010-4655. S2CID 6140858. Retrieved 2020-07-13.
17. Tsirkin, Stepan S.; Puente, Pablo Aguado; Souza, Ivo (2018-01-29). "Gyrotropic effects in trigonal tellurium studied from first principles". Physical Review B. 97 (3): 035158. arXiv:1710.03204. Bibcode:2018PhRvB..97c5158T. doi:10.1103/physrevb.97.035158. ISSN 2469-9950. S2CID 55517213.
18. Ibañez-Azpiroz, Julen; Tsirkin, Stepan S.; Souza, Ivo (2018-06-26). "Ab initio calculation of the shift photocurrent by Wannier interpolation". Physical Review B. 97 (24): 245143. arXiv:1804.04030. Bibcode:2018PhRvB..97x5143I. doi:10.1103/physrevb.97.245143. ISSN 2469-9950. S2CID 67751414.
19. Qiao, Junfeng; Zhou, Jiaqi; Yuan, Zhe; Zhao, Weisheng (2018-12-03). "Calculation of intrinsic spin Hall conductivity by Wannier interpolation". Physical Review B. 98 (21): 214402. arXiv:1810.07637. Bibcode:2018PhRvB..98u4402Q. doi:10.1103/physrevb.98.214402. ISSN 2469-9950. S2CID 119223848.
20. Ryoo, Ji Hoon; Park, Cheol-Hwan; Souza, Ivo (2019-06-07). "Computation of intrinsic spin Hall conductivities from first principles using maximally localized Wannier functions". Physical Review B. 99 (23): 235113. arXiv:1906.07139. Bibcode:2019PhRvB..99w5113R. doi:10.1103/physrevb.99.235113. ISSN 2469-9950. S2CID 189928182.

Further reading

• Karin M Rabe; Jean-Marc Triscone; Charles H Ahn (2007). Physics of Ferroelectrics: a Modern Perspective. Springer. p. 2. ISBN 978-3-540-34590-9.

External links

• Wannier Gregory H (1937). "The Structure of Electronic Excitation Levels in Insulating Crystals". Physical Review. 52 (3): 191–197. Bibcode:1937PhRv...52..191W. doi:10.1103/PhysRev.52.191.
• Wannier90 computer code that calculates maximally localized Wannier functions
• Wannier Transport code that calculates maximally localized Wannier functions fit for Quantum Transport applications
• WannierTools: An open-source software package for novel topological materials
• WannierBerri - a python code for Wannier interpolation and tight-binding calculations

See also

• Bloch's theorem
• Hannay angle
• Geometric phase