Maximally localized Wannier functions: Theory and applications
Maximally localized Wannier functions: Theory and applications
Maximally localized Wannier functions: Theory and applications
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
41<br />
FIG. 31 (Color online) Surface density-of-states (SDOS) of a<br />
semi-infinite crystal of Sb 2 Te 3 terminated with a [111] surface.<br />
Lighter (warmer) colors represent a higher SDOS. The surface<br />
states can be seen around Γ as (red) lines crossing at E = 0.<br />
From Zhang et al. (2009a).<br />
From these, the density of states (DOS) projected onto<br />
a given atomic plane P can be obtained(Grosso <strong>and</strong> Parravicini,<br />
2000) as<br />
N P l<br />
(k ∥ , ϵ) = − 1 π Im ∑ n∈P<br />
G nn<br />
ll (k ∥ , ϵ + iη), (105)<br />
where the sum over n is restricted to the orbitals ascribed<br />
to the chosen plane <strong>and</strong> η is a positive infinitesimal.<br />
The projection of the DOS onto the outermost atomic<br />
plane is shown in Fig. 31 as a function of energy ϵ <strong>and</strong><br />
momentum k ∥ for the (111) surface of Sb 2 Te 3 . The same<br />
method has been used to find the dispersion of the surface<br />
b<strong>and</strong>s in the TI alloy Bi 1−x Sb x (Zhang et al., 2009b) <strong>and</strong><br />
in ternary compounds with a honeycomb lattice (Zhang<br />
et al., 2011b).<br />
B. B<strong>and</strong> derivatives<br />
The first <strong>and</strong> second derivatives of the energy eigenvalues<br />
with respect to k (b<strong>and</strong> velocities <strong>and</strong> inverse effective<br />
masses) appear in a variety of contexts, such as the calculation<br />
of transport coefficients (Ashcroft <strong>and</strong> Mermin,<br />
1976; Grosso <strong>and</strong> Parravicini, 2000). There is therefore<br />
considerable interest in developing simple <strong>and</strong> accurate<br />
procedures for extracting these parameters from a firstprinciples<br />
b<strong>and</strong> structure calculation.<br />
A direct numerical differentiation of the eigenenergies<br />
calculated on a grid is cumbersome <strong>and</strong> becomes unreliable<br />
near b<strong>and</strong> crossings. It is also very expensive if<br />
a Brillouin zone integration is to be carried out, as in<br />
transport calculations. A number of efficient interpolation<br />
schemes, such as the method implemented in the<br />
BoltzTraP package (Madsen <strong>and</strong> Singh, 2006), have<br />
been developed for this purpose, but they are still prone<br />
to numerical instabilities near b<strong>and</strong> degeneracies (Uehara<br />
<strong>and</strong> Tse, 2000). Such instabilities can be avoided by using<br />
a tight-binding parametrization to fit the first-principles<br />
b<strong>and</strong> structure (Mazin et al., 2000; Schulz et al., 1992).<br />
As shown by Graf <strong>and</strong> Vogl (1992) <strong>and</strong> Boykin (1995),<br />
both the first <strong>and</strong> the second derivatives are easily computed<br />
within tight-binding methods, even in the presence<br />
of b<strong>and</strong> degeneracies, <strong>and</strong> the same can be done in a firstprinciples<br />
context using WFs.<br />
Let us illustrate the procedure by calculating the b<strong>and</strong><br />
gradient away from points of degeneracy; the treatment<br />
of degeneracies <strong>and</strong> second derivatives is given in Yates<br />
et al. (2007). The first step is to take analytically the<br />
derivative ∂ α = ∂/∂k α of Eq. (98),<br />
H W k,α ≡ ∂ α H W k<br />
= ∑ R<br />
e ik·R iR α ⟨0|H|R⟩. (106)<br />
The actual b<strong>and</strong> gradients are given by the diagonal elements<br />
of the rotated matrix,<br />
]<br />
∂ α ϵ nk =<br />
[U † k HW k,αU k (107)<br />
where U k is the same unitary matrix as in Eq. (99).<br />
It is instructive to view the columns of U k as orthonormal<br />
state vectors in the J-dimensional “tightbinding<br />
space” defined by the WFs. According to<br />
Eq. (99) the n-th column vector, which we shall denote<br />
by ||ϕ nk ⟩⟩, satisfies the eigenvalue equation Hk W||ϕ nk⟩⟩ =<br />
ϵ nk ||ϕ nk ⟩⟩. Armed with this insight, we now recognize<br />
in Eq. (107) the Hellmann-Feynman result ∂ α ϵ nk =<br />
⟨⟨ϕ nk ||∂ α Hk W||ϕ nk⟩⟩.<br />
1. Application to transport coefficients<br />
Within the semiclassical theory of transport, the electrical<br />
<strong>and</strong> thermal conductivities of metals <strong>and</strong> doped<br />
semiconductors can be calculated from a knowledge of<br />
the b<strong>and</strong> derivatives <strong>and</strong> relaxation times τ nk on the<br />
Fermi surface. An example is the low-field Hall conductivity<br />
σ xy of non-magnetic cubic metals, which in the constant<br />
relaxation-time approximation is independent of τ<br />
<strong>and</strong> takes the form of a Fermi-surface integral containing<br />
the first <strong>and</strong> second b<strong>and</strong> derivatives (Hurd, 1972).<br />
Previous first-principles calculations of σ xy using various<br />
interpolation schemes encountered difficulties for materials<br />
such as Pd, where b<strong>and</strong> crossings occur at the<br />
Fermi level (Uehara <strong>and</strong> Tse, 2000). A <strong>Wannier</strong>-based<br />
calculation free from such numerical instabilities was carried<br />
out by Yates et al. (2007), who obtained carefullyconverged<br />
values for σ xy in Pd <strong>and</strong> other cubic metals.<br />
A more general formalism to calculate the electrical<br />
conductivity tensor in the presence of a uniform magnetic<br />
field involves integrating the equations of motion of<br />
a wavepacket under the field to find its trajectory on the<br />
Fermi surface (Ashcroft <strong>and</strong> Mermin, 1976). A numerical<br />
implementation of this approach starting from the<br />
nn