Berry phases and finite electric fields
Modern Theory of Polarization
Berry phase expression for the macroscopic polarization
Calculating the change in dipole moment per unit cell under PBC's, is a nontrivial task. In general one cannot define it as the first moment of the induced change in charge density δ(r), through
without introducing a dependency on the shape of Ω0, the chosen unit cell.[1]
Recently King-Smith and Vanderbilt[2], building on the work of Resta[3], showed that the electronic contribution to the difference in polarization ΔPe, due to a finite adiabatic change in the Hamiltonian of a system, can be identified as a geometric quantum phase or Berry phase of the valence wave functions. We will briefly summarize the essential results (for a review of geometric quantum phases in polarization theory see the papers of Resta[4][5]).
Central to the modern theory of polarization is the proposition of Resta[3] to write the electronic contribution to the change in polarization due to a finite adiabatic change in the Kohn-Sham Hamiltonian of the crystalline solid, as
with
where me and e are the electronic mass and charge, N is the number of unit cells in the crystal, Ω0 is the unit cell volume, M is the number of occupied bands, p is the momentum operator, and the functions ψ(λ)nk are the usual Bloch solutions to the crystalline Hamiltonian. Within Kohn-Sham density-functional theory, the potential V(λ) is to be interpreted as the Kohn-Sham potential V(λ)KS, where λ parameterizes some change in this potential, for instance due to the displacement of an atom in the unit cell.
King-Smith and Vanderbilt[2] have cast this expression in a form in which the conduction band states ψ(λ)mk no longer explicitly appear, and they show that the change in polarization along an arbitrary path, can be found from only a knowledge of the system at the end points
with
where f is the occupation number of the states in the valence bands, u(λ)nk is the cell-periodic part of the Bloch function ψ(λ)nk, and the sum n runs over all M occupied bands.
The physics behind the equation above becomes more transparent when this expression is written in terms of the Wannier functions of the occupied bands,
where Wn is the Wannier function corresponding to valence band n.
This shows the change in polarization of a solid, induced by an adiabatic change in the Hamiltonian, to be proportional to the displacement of the charge centers rn=⟨ W(λ)n|r| W(λ)n⟩, of the Wannier functions corresponding to the valence bands.
It is important to realize that the polarization in terms of Bloch or Wannier functions, and consequently the change in polarization, is only well-defined modulo feR/Ω0, where R is a lattice vector. This indeterminacy stems from the fact that the charge center of a Wannier function is only invariant modulo R, with respect to the choice of phase of the Bloch functions.
In practice one is usually interested in polarization changes |ΔPe| << |feR1/Ω0|, where R1 is the shortest nonzero lattice vector. An arbitrary term feR/Ω0 can therefore often be removed by simple inspection of the results. In cases where |ΔPe| is of the same order of magnitude as feR1/Ω0 any uncertainty can always be removed by dividing the total change in the Hamiltonian λ1→λ2 into a number of intervals.
Computational aspects
In general, the direct evaluation of Pe(λ) is useless, because there is no specific relationship between the phases of the eigenvectors u(λ)kn generated by a numerical diagonalization routine. This problem is circumvented by dividing the Brillouin zone integration in two parts, a two-dimensional integral and a line integral, and by transforming the above into three equations, which separately provide the components of Pe(λ) along the directions of three reciprocal lattice vectors G1, G2, and G3, which together span a unit cell of the reciprocal lattice. The component of Pe(λ) along for instance G1 can be found from
where the two-dimensional integral is taken over the area A, spanned by G2 and G3, and the line integral runs over a line segment parallel to G1. Interchanging the indices 1, 2 and 3 in the equation above yields the expressions for two other components of Pe(λ).
Thus the electronic part of the polarization Pe(λ) is given (modulo feR/Ω0) by the sum
where the lattice vectors Ri obey the relationship Ri·Gj=2πδij.
The integration over A in the above, is straightforward and can be performed by sampling a 2D Monkhorst-Pack mesh of k-points,[6] termed the perpendicular mesh or k⊥-mesh by King-Smith and Vanderbilt. However, to remove the influence of the random phase of the functions u(λ)nk, introduced by the diagonalization routine, King-Smith and Vanderbilt[2] propose to replace the line integral alias integration in the parallel or G|| direction by,
which is evaluated by calculating the cell-periodic parts of the wave functions at a string of J k-points, kj= k⊥+jG||/J (with j=0,..,J-1), and where for sufficiently large J one has that
Note: the determinant appearing in the equation above is the determinant of the M×M matrix formed by letting n and m run over all valence bands.
The crucial step, instrumental in removing the random phase, is that the functions u(λ)nkJ are not obtained from an independent diagonalization, but found through their relationship with the functions u(λ)nk0,
This way the product φJ(λ) becomes cyclic, and contains both u(λ)nkj as well as its complex conjugate for every k-point in the string, thus removing the random phase.
In practice Eq.\ (\ref{Polarization2}) is evaluated by way of the following summation over the $\mathbf{k}_{\perp}$-mesh,
\noindent where we used $\mathbf{R}_{i}\cdot\mathbf{G}_{i}=2\pi$. Since $\phi^{\left(\lambda\right)}_{J}\left(\mathbf{k}_{\perp}\right)$ as given by Eq.\ (\ref{CyclicForm}), is only well defined modulo $2\pi$ ($\text{Im}\{\ln z\}$, with $z \in \mathbb{C}$, is only defined up to an integer multiple of $2\pi$), evaluation of Eq.\ (\ref{Polarization3}) yields $({\bf P}_{e}^{\left(\lambda\right)})_i{}$ modulo $fe{\bf R}_{i}/\Omega_{0}N_{k_{\perp}}$ (instead of modulo $fe\mathbf{R}/\Omega_{0}$ as Eqs.\ (\ref{Polarization}) and (\ref{Wannier}) in Sec.\ \ref{BerryPhaseTechnique}). The change in polarization $\Delta \mathbf{P}_{e}$ however, can be obtained modulo $fe\mathbf{R}/\Omega_{0}$. To this end we compare $\phi^{\left(\lambda_{1}\right)}_{J}\left(\mathbf{k}_{\perp}\right)$ to $\phi^{\left(\lambda_{2}\right)}_{J}\left(\mathbf{k}_{\perp}\right)$, and remove the occurrence of an arbitrary multiple of $2\pi$, on the premise that $|\phi^{\left(\lambda_{2}\right)}_{J}\left(\mathbf{k}_{\perp}\right) - \phi^{\left(\lambda_{1}\right)}_{J}\left(\mathbf{k}_{\perp}\right)| \ll 2\pi$.
The formulation of the Berry phase technique as presented by Eqs.\ (\ref{Polarization})-(\ref{Polarization3}) can be applied in the context of any norm-conserving pseudopotential or all-electron scheme. Within the context of the ultrasoft pseudopotential formalism\ \cite{Vanderbilt90,Kresse94} however, Eqs.\ (\ref{Polarization})-(\ref{Polarization2}), (\ref{CyclicForm}), (\ref{CyclicFormLimit}) and\ (\ref{Polarization3}) must be modified \cite{Vanderbilt98}. This modified formulation, and information with respect to its implementation in {\sc VASP}, can be found below, in Appendix A.
Self-consistent response to finite electric fields
Related Tags and Sections
LBERRY, IGPAR, NPPSTR, LPEAD, IPEAD, LCALCPOL, LCALCEPS, EFIELD_PEAD, DIPOL
References
- ↑ P. Vogl, J. Phys. C: Solid State Phys. 11, 251 (1978).
- ↑ a b c R. D. King-Smith and D. Vanderbilt, Phys. Rev. B 47, 1651 (1993).
- ↑ a b R. Resta, Ferroelectrtics 136, 51 (1992).
- ↑ R. Resta, Rev. Mod. Phys. 66, 899 (1994).
- ↑ R. Resta, in Berry Phase in Electronic Wavefunctions, Troisième Cycle de la Physique en Suisse Romande, Ann%eacute;e Academique 1995-96, (1996).
- ↑ H. J. Monkhorst and J. D. Pack, Phys. Rev. B 13, 1588 (1976).
Cite error: <ref>
tag with name "Vanderbilt90" defined in <references>
is not used in prior text.
Cite error: <ref>
tag with name "Kresse94" defined in <references>
is not used in prior text.
Cite error: <ref>
tag with name "Vanderbilt98" defined in <references>
is not used in prior text.