Skip to content

Polarization and Thouless Pumping

Bulk polarization in a crystal is not an absolute dipole moment: it is defined modulo the charge transported by shifting one electron through a lattice period. Its change along a gapped path is physical, and the charge transported in a closed adiabatic cycle is an integer Chern number. The gap, cycle orientation, occupied subspace, and unit cell are part of the statement.

Required background. Berry geometry supplies the connection and curvature conventions.

Helpful background. Background responses and invertible phases supplies the effective-response interpretation.

Consider a one-dimensional insulator with lattice period aa, cell-periodic occupied states unk|u_{nk}\rangle, and electron charge e-e. With kk oriented from 00 to 2π/a2\pi/a and Ak=inoccunkkunk\mathcal A_k=i\sum_{n\in\mathrm{occ}}\langle u_{nk}|\partial_k u_{nk}\rangle, the electronic polarization per cell is

Pel=e2πBZAkdk(mode).P_{\rm el}=-\frac{e}{2\pi}\int_{\rm BZ}\mathcal A_k\,dk \pmod e.

A periodic frame transformation with determinant winding ww changes the integral by 2πw-2\pi w, hence PelPel+ewP_{\rm el}\mapsto P_{\rm el}+ew. This is the polarization quantum, not an error. Ionic charges and their chosen cell positions must be added before comparing with a measured total polarization Resta 1994, §§ II–III.

Inversion symmetry quantizes PP to 00 or e/2e/2 modulo ee: inversion sends PPP\mapsto-P, so 2P=02P=0 modulo ee. Moving the spatial origin or redefining the unit cell can shift the electronic and ionic pieces; the total symmetry-quantized class remains the invariant after the same convention is used on both sides.

Let H(k,t)H(k,t) be periodic in tt with period TT and gapped for every (k,t)(k,t). Differentiating the Berry-phase formula and using Ωkt=kAttAk\Omega_{kt}=\partial_k\mathcal A_t-\partial_t\mathcal A_k gives

Q+x=ΔP=e2π0TdtBZdkΩkt=eCkt,CktZ.Q_{+x}=\Delta P =\frac{e}{2\pi}\int_0^Tdt\int_{\rm BZ}dk\,\Omega_{kt} =eC_{kt}, \qquad C_{kt}\in\mathbb Z.

This sign corresponds to the stated (k,t)(k,t) orientation and positive current toward +x+x. Reversing the cycle reverses CktC_{kt} and the transported charge. Quantization follows because the occupied states form a bundle over the (k,t)(k,t) torus Thouless 1983.

For a finite open chain, boundary charge changes as the pump proceeds and edge levels may cross the chemical potential. Those crossings are the boundary realization of bulk transport; the instantaneous boundary charge is termination dependent. In a disordered or interacting ring, a twisted boundary phase replaces kk, and a unique many-body gap throughout the (θ,t)(\theta,t) torus is the relevant hypothesis.

Write the cycle as H(s)H(s) with s=t/Tcyc[0,1]s=t/T_{\mathrm{cyc}}\in[0,1]. For instantaneous many-body eigenstates m(s)|m(s)\rangle, a sufficient finite-system adiabatic condition is

maxm0,sm(s)sH(s)0(s)Tcyc[Em(s)E0(s)]21,\max_{m\ne0,\,s} \frac{\hbar\,|\langle m(s)|\partial_sH(s)|0(s)\rangle|} {T_{\mathrm{cyc}}[E_m(s)-E_0(s)]^2}\ll1,

with a nonzero minimum gap Δmin\Delta_{\min} throughout the cycle. Thus TcycT_{\mathrm{cyc}} must be long compared with the transition-matrix-element scale, not merely with /Δmin\hbar/\Delta_{\min}, and short compared with heating or decoherence times. A bulk gap closing permits CC to change. Partial filling, Landau–Zener excitation, noncyclic driving, or coupling to reservoirs can make transported charge noninteger. Quantization of a change does not make a particular branch of absolute polarization unique.

An inversion-symmetric one-dimensional insulator has P=e/2P=e/2 modulo ee. Show that this is compatible with inversion.

Solution

Inversion requires P=PP=-P modulo ee, or 2P=ne2P=ne. For P=e/2P=e/2, e/2=e/2e-e/2=e/2-e, so the two values differ by one polarization quantum and represent the same bulk class.

  • Raffaele Resta, “Macroscopic Polarization in Crystalline Dielectrics: The Geometric Phase Approach,” Reviews of Modern Physics 66 (1994) 899–915, doi:10.1103/RevModPhys.66.899.
  • David J. Thouless, “Quantization of Particle Transport,” Physical Review B 27 (1983) 6083–6087, doi:10.1103/PhysRevB.27.6083.