Berry-Phase Polarization and Charge Pumping

The dipole moment of a chosen crystal cell is not, by itself, the bulk electric polarization of an infinite periodic solid. Charges can be reassigned across a cell boundary without changing the crystal, and occupied Bloch states are extended rather than attached to unique cells. The physically appropriate bulk object is therefore an equivalence class of polarization values, together with changes tracked along a specified insulating path.

Berry-phase polarization joins the ionic charge positions to the holonomy of the complete occupied electronic subspace. It is defined modulo a polarization quantum. When an insulating Hamiltonian executes a closed, charge-conserving adiabatic cycle, a continuous polarization branch can wind by an integer number of quanta. That winding is a Chern number and measures the charge transported across the sample.

This page owns that complete material workflow. The abstract connection and curvature remain with the Berry-geometry pages; localized-basis existence and obstruction remain with Wannier theory; static Brillouin-zone Chern response remains with Chern Numbers in Band Theory.

Required background. You should be able to identify a fixed-rank occupied Bloch subspace and use the convention $\mathcal A=+i\langle u|\partial u\rangle$. Bloch’s Theorem supplies the Brillouin-zone and reciprocal-sewing conventions, while Berry Connection supplies the gauge transformation law.

Helpful background. Conventions for Quantum Matter fixes charge, cell, and origin language. Zak Phase Preview provides one-dimensional holonomy intuition, and Wannier Functions explains localized frames and their obstructions. Before using the pump sections, review Adiabatic Approximation as a Method. Projectors is a useful compact review of the occupied-subspace algebra used below.

Declare a Crystalline Polarization Problem

Fix the convention ledger before evaluating a Berry phase. At minimum, record:

1. the primitive lattice vectors $\mathbf a_i$, reciprocal vectors $\mathbf b_j$, and the orientation for which $\mathbf a_i\cdot\mathbf b_j=2\pi\delta_{ij}$ and $\Omega_c>0$;
2. the real-space origin, basis positions, cell assignment, and boundary normal used for any surface statement;
3. $e>0$, electron charge $q_e=-e$, and ionic or pseudopotential-core charges $+Z_\kappa e$ at positions $\boldsymbol\tau_\kappa$;
4. whether spin is explicit in each occupied spinor; no factor of two is implicit on this page;
5. the state, filling, occupied rank, and projector $P(\mathbf k)$, including the minimum occupied-empty direct gap and the global or mobility gap that makes the intended state insulating;
6. the Bloch phase, cell normalization, reciprocal sewing, gauge, mesh, and Wilson-loop orientation;
7. the target: a polarization class, a branch-resolved $\Delta\mathbf P$, a response derivative, interface charge, current integral, or pumped charge;
8. the structural or Hamiltonian path, its direction, its endpoint relation, and the branch reference used to unwrap polarization;
9. the symmetry, if any, claimed to quantize a component or relate two branches; and
10. finite-size, temperature, drive-rate, leakage, screening, disorder, interaction, and numerical uncertainties relevant to the final claim.

The band projector must have fixed rank and remain directly separated from its complement wherever the occupied-subspace formula is used. Internal crossings among occupied bands are allowed. A positive direct separation alone does not make the physical state insulating: valence and conduction extrema may still overlap indirectly, leaving electron and hole pockets at the declared filling.

Electronic Polarization from an Occupied Subspace

Write the Bloch states as

$$|\psi_{n\mathbf k}\rangle = e^{i\mathbf k\cdot\mathbf r} |u_{n\mathbf k}\rangle,$$

with cell-periodic spinors normalized in the declared cell convention. Collect an orthonormal frame for the $N_{\mathrm{occ}}$ occupied states into $U(\mathbf k)$. The occupied projector is

$$P(\mathbf k) = U(\mathbf k)U^\dagger(\mathbf k).$$

It is unchanged by $U\mapsto UV$ for any $V(\mathbf k)\in U(N_{\mathrm{occ}})$, so it remains well defined when individual occupied eigenvectors cross or exchange labels. The frame’s non-Abelian Berry connection is

$$[\mathcal A_i(\mathbf k)]_{mn} = i\langle u_{m\mathbf k}|\partial_{k_i}u_{n\mathbf k}\rangle_c.$$

Changing the occupied frame by $U(\mathbf k)\mapsto U(\mathbf k)V(\mathbf k)$, with $V\in U(N_{\mathrm{occ}})$, gives

$$\mathcal A_i \mapsto V^\dagger\mathcal A_iV + iV^\dagger\partial_{k_i}V.$$

Individual diagonal connections and energy-sorted band labels are therefore not reliable at occupied-band crossings. The total occupied projector and the trace holonomy are the correct objects.

In one dimension, include the reciprocal endpoint sewing and define

$$\gamma_{\mathrm{occ}} = \oint_{\mathrm{BZ}} dk\, \operatorname{Tr}\mathcal A_k \pmod{2\pi}.$$

For a forward mesh $k_0,k_1,\ldots,k_N=k_0+G$, form overlap matrices

$$[M_j]_{mn} = \langle u_m(k_j)|u_n(k_{j+1})\rangle_c.$$

The final state must be the representative sewn back to the initial fiber. A stable numerical Wilson loop uses the unitary polar factors

$$\widetilde M_j = M_j(M_j^\dagger M_j)^{-1/2}, \qquad \mathcal W = \widetilde M_0\widetilde M_1\cdots\widetilde M_{N-1}.$$

With the connection and forward-link convention used here,

$$\gamma_{\mathrm{occ}} = -\operatorname{Arg}\det\mathcal W \pmod{2\pi}.$$

The minus sign is not optional bookkeeping: a forward overlap is $1-i\mathcal A_k\,dk+O(dk^2)$. Reversing the loop or defining links in the opposite direction changes the corresponding sign convention.

Under a sewn occupied-frame change, $\mathcal W$ transforms by conjugation at the base point. Its determinant and unordered eigenvalue set are invariant. Across a parameter path, Wilson eigenphases may cross and permute; follow the subspace or the full center set continuously rather than sorting principal values independently at each step.

Add Ions and Identify the Polarization Quantum

For $e>0$, the one-dimensional electronic polarization per transverse cell is

$$P_{\mathrm{el}} = -\frac{e}{2\pi}\gamma_{\mathrm{occ}} \pmod e.$$

It is modulo $e$, not modulo $e/a$. In $d$ dimensions, for a globally smooth, reciprocally sewn occupied frame in a fixed cell,

$$\mathbf P_{\mathrm{el}} = -\frac{e}{(2\pi)^d} \int_{\mathrm{BZ}} d^d k\, \operatorname{Tr}\boldsymbol{\mathcal A}(\mathbf k).$$

Each explicitly included occupied spinor contributes once. A calculation that uses spinless spatial orbitals to represent two exactly degenerate occupied spin states must state and insert the factor of two.

The ionic contribution is

$$\mathbf P_{\mathrm{ion}} = \frac{e}{\Omega_c} \sum_\kappa Z_\kappa\boldsymbol\tau_\kappa,$$

so the total is

$$\mathbf P = \mathbf P_{\mathrm{ion}} + \mathbf P_{\mathrm{el}}.$$

Moving one electron center by a lattice vector does not change the infinite crystal, but it changes the representative by

$$\mathbf P_q(\mathbf R) = \frac{e\mathbf R}{\Omega_c}.$$

Consequently,

$$\mathbf P \sim \mathbf P + \frac{e\mathbf R}{\Omega_c}.$$

For particles of another conserved charge, replace $e$ by the charge quantum appropriate to the transported species and state the sign. Ionic and electronic pieces separately depend on origin and cell assignment. For a neutral cell, their origin shifts cancel in the total. Reassigning whole charges between neighboring cells still changes the displayed representative by a polarization quantum.

The full-dimensional expression also assumes that the relevant occupied bundle admits the needed smooth sewing. A static Chern insulator has no global smooth periodic frame, so a naive bulk polarization vector built from that formula is not licensed. Hybrid Wannier centers or slice-resolved quantities may remain meaningful, but they require their own direction and boundary ledger.

Wilson Loops, Wannier Centers, and Boundary-Charge Limits

When localized Wannier functions exist, the electronic polarization has the equivalent center representation

$$\mathbf P_{\mathrm{el}} = -\frac{e}{\Omega_c} \sum_{n=1}^{N_{\mathrm{occ}}} \overline{\mathbf r}_n \pmod{\frac{e\mathbf R}{\Omega_c}}.$$

In one dimension,

$$\sum_n\overline x_n = \frac{a}{2\pi}\gamma_{\mathrm{occ}} \pmod a.$$

This is a representation of the occupied charge center, not a statement that each individual Wannier center is an observable. A $U(N_{\mathrm{occ}})$ frame rotation can exchange or redistribute centers while preserving the occupied projector and their total modulo a lattice vector. Wannier Functions owns localization, gauge choice, and obstruction theory. Wannierization Workflows owns numerical construction and full-zone validation of the localized subspace and any interpolated position or velocity operators used downstream.

For a two-dimensional band structure, a Wilson loop along $k_x$ at fixed $k_y$ produces hybrid centers $\overline x_a(k_y)$. Their unordered flow can diagnose a Chern obstruction or a symmetry-protected pattern. That flow is not an extra set of electron trajectories. It is a gauge-invariant spectrum of an occupied-subspace holonomy.

Polarization constrains boundary charge only after the boundary is declared. For an insulating interface with normal $\mathbf n$ directed from region 1 to region 2,

$$\sigma_{\mathrm b} = (\mathbf P_1-\mathbf P_2)\cdot\mathbf n$$

modulo the charge per surface primitive cell. An outer surface is the special case with vacuum polarization chosen as the reference. Changing termination can add an integer surface charge; reconstruction, metallic surface bands, adsorbates, defects, electrodes, and mobile screening charge can change the measured value. Bulk polarization therefore does not predict an absolute end charge without a termination and free-charge ledger.

Track Polarization Along a Gapped Path

Let an insulating Hamiltonian depend continuously on $\lambda\in[0,1]$. The fixed-rank occupied subspace must remain separated from empty states, and the physical filling must remain insulating, along the complete path. Define

$$\mathcal F_{k_i\lambda} = \partial_{k_i}\mathcal A_\lambda - \partial_\lambda\mathcal A_{k_i} - i[\mathcal A_{k_i},\mathcal A_\lambda].$$

At fixed lattice vectors, the electronic change is

$$\Delta P_{\mathrm{el},i} = \frac{e}{(2\pi)^d} \int_0^1 d\lambda \int_{\mathrm{BZ}}d^d k\, \operatorname{Tr}\mathcal F_{k_i\lambda}.$$

Add the ionic displacement on the same continuous branch. If the cell itself changes, use reduced polarization coordinates and distinguish proper from improper strain responses; the fixed-cell expression above is not sufficient.

Endpoint Berry phases alone do not select the branch of $\Delta\mathbf P$. Sample the path densely enough that neighboring representatives can be joined continuously, inspect occupied-subspace overlaps and the minimum gap, and add or subtract polarization quanta only to preserve that continuity. A nearest-branch rule is unsafe if one step can move by more than half a quantum; refine the path instead.

For a spatially uniform bulk at fixed geometry, charge conservation gives

$$\mathbf J(t) = \frac{d\mathbf P}{dt}, \qquad \Delta\mathbf P = \int_{t_i}^{t_f}dt\,\mathbf J(t).$$

This relation is why switching current can determine a polarization change. It does not turn a leaky terminal-current trace directly into bulk polarization: capacitive background, free-carrier conduction, electrode screening, domain fractions, and incomplete switching require a forward model.

A Closed Adiabatic Cycle Pumps Charge

Now let a translation-invariant one-dimensional electronic Hamiltonian depend cyclically on $\lambda\in[0,2\pi]$:

$$H(k,\lambda+2\pi) = H(k,\lambda).$$

Before applying the result, require conserved $U(1)$ charge, fixed filling and occupied rank, a gap for every $(k,\lambda)$, and adiabatic following. Hold the ionic coordinates fixed for the displayed electronic theorem, or track their current separately. Orient the parameter torus by $dk\wedge d\lambda$ and define

$$C_{k\lambda} = \frac{1}{2\pi} \int_0^{2\pi}d\lambda \int_{\mathrm{BZ}}dk\, \operatorname{Tr}\mathcal F_{k\lambda}.$$

This is an integer for the continuous occupied bundle over the torus. With the electron and orientation conventions fixed above, the electronic contribution to conventional positive charge transported toward $+x$ in one cycle is

$$Q_{\mathrm{el},+x} = eC_{k\lambda}.$$

The sign follows from $P_{\mathrm{el}}=-e\int\mathcal A_k/(2\pi)$ and $\int\mathcal F_{k\lambda}=-\partial_\lambda\int\mathcal A_k$ after the $k$-boundary term is sewn away. Defining $\mathcal F_{\lambda k}$, reversing the torus orientation, or asking for electron rather than conventional-charge transport reverses the displayed sign. For $C_{k\lambda}=+1$, the aggregate electronic center flow is by $-a$, while conventional charge $+e$ moves toward $+x$.

The Hamiltonian returns to itself, so the endpoint polarization classes agree. The continuously unwrapped branch need not return to the same representative:

$$P_{\mathrm{el}}(2\pi) = P_{\mathrm{el}}(0) + eC_{k\lambda}$$

in the one-dimensional convention. If ions move, the total change is

$$\Delta P_{\mathrm{total}} = eC_{k\lambda} + \Delta P_{\mathrm{ion}}.$$

A closed structural cycle can carry an additional ionic polarization quantum. The branch winding, rather than a supposedly absolute endpoint polarization, is the pump invariant. A generic static Zak phase is neither quantized nor a Chern number.

The continuous Brillouin-zone integral describes the infinite periodic band system. On a finite ring, the charge integrated at one fixed boundary twist can vary continuously. Twist averaging recovers the Chern number, while equality at one fixed twist requires a thermodynamic-limit or twist-insensitivity argument.

Exact quantization is an adiabatic filled-state statement. Finite drive rate, thermal excitations, gap inhomogeneity, particle loss, coupling to reservoirs, and imperfect cycle closure produce platform-dependent corrections. Do not attach a universal error bar merely from the minimum gap; validate convergence with cycle time and inspect the actual transition channels.

For an open finite chain, topological pumping generally appears through edge spectral flow or exponentially small edge avoided crossings. Following one unique isolated open-system ground state around the full cycle is not the same theorem as the periodic bulk Chern number. State whether edge occupations are reset by reservoirs, followed diabatically across tiny edge gaps, or measured as accumulated boundary charge.

Symmetry Can Quantize a Polarization Class

A symmetry quantizes polarization only if it maps the declared state and charge distribution back to themselves while constraining the polarization class. In one dimension, inversion sends

$$P \mapsto -P \pmod e.$$

An inversion-symmetric gapped state therefore satisfies

$$P = 0 \quad\text{or}\quad \frac e2 \pmod e.$$

The claim concerns the total ionic plus electronic polarization, with an inversion center, origin, cell, and filling declared. The electronic SSH Zak phase can be $0$ or $\pi$ in a symmetry-compatible convention, but changing orbital embedding or unit cell changes its representative. The ionic reference and boundary termination are still required before inferring end charge.

Time reversal alone does not quantize ordinary electric polarization in one dimension because electric polarization is even under time reversal. Mirrors and rotations can constrain selected components, while nonsymmorphic or magnetic symmetries require their complete spatial and antiunitary actions. Symmetry of Bloch States owns those representation data.

A quantized Thouless pump does not need a point symmetry at each intermediate time. Its integer follows from charge conservation, cyclicity, the gap, and the Chern class over the complete cycle. Symmetry can constrain endpoints or parts of a cycle, but it is not a substitute for those pump conditions.

Interacting, Disordered, and Finite-System Extensions

Crystal momentum is unavailable in a disordered sample, and a generic interacting ground state is not a list of occupied one-particle bands. On a one-dimensional periodic system of length $L$, let $\hat X=\sum_j\hat x_j$ denote the coordinate sum only inside the large-gauge exponential

$$\hat U_X = \exp\!\left( \frac{2\pi i}{L}\hat X \right).$$

The bare $\hat X$ is not a single-valued position observable on the periodic ring; $\hat U_X$ is the well-defined periodic operator.

For a localized insulating state with $z=\langle\Psi|\hat U_X|\Psi\rangle\ne0$, an electronic polarization representative is

$$P_{\mathrm{el}} = -\frac{e}{2\pi} \operatorname{Im}\log z \pmod e.$$

Add ions and retain the same branch discipline. A nonzero value of $z$ in one small finite system does not by itself prove a thermodynamic insulating phase; study size scaling, the charge gap or localization criterion, and sensitivity to boundary twists.

For an interacting or disordered pump, impose a boundary twist $\theta\in[0,2\pi]$ and a cyclic parameter $\lambda$. If a unique many-body ground state remains separated by a gap over the complete $(\theta,\lambda)$ torus, set

$$\mathcal A_\mu^{\mathrm{MB}} = i\langle\Psi(\theta,\lambda)| \partial_\mu\Psi(\theta,\lambda)\rangle,$$

and

$$\mathcal F_{\theta\lambda}^{\mathrm{MB}} = \partial_\theta\mathcal A_\lambda^{\mathrm{MB}} - \partial_\lambda\mathcal A_\theta^{\mathrm{MB}}.$$

Orient the torus by $d\theta\wedge d\lambda$ and define

$$C_{\mathrm{MB}} = \frac{1}{2\pi} \int_0^{2\pi}d\lambda \int_0^{2\pi}d\theta\, \mathcal F_{\theta\lambda}^{\mathrm{MB}}.$$

Then $eC_{\mathrm{MB}}$ is the boundary-twist-averaged adiabatically pumped charge in the matching orientation convention. Pumped charge at one fixed twist approaches that value only with thermodynamic twist insensitivity; verify its size and twist dependence. Degenerate ground-state manifolds, fractional pumps, open-system steady states, and nonconserved charge require a different non-Abelian or dynamical statement; do not force them into the unique-ground-state formula. Boundary Conditions on Lattices owns the twist construction.

Numerically, report the minimum gap, smallest overlap singular value, mesh and path refinement, reciprocal sewing residual, Wilson unitarization rule, branch-tracking algorithm, and current-integral cross-check. A rounded integer is not evidence of convergence if overlaps are nearly singular or curvature is unresolved.

Worked Audit: A Rice–Mele Pump

Consider a spinless two-orbital chain with one occupied band,

$$H(k,\lambda) = \mathbf d(k,\lambda)\cdot\boldsymbol\sigma,$$

where

$$\begin{aligned} d_x&=t_1(\lambda)+t_2(\lambda)\cos(ka),\\ d_y&=t_2(\lambda)\sin(ka),\\ d_z&=\Delta(\lambda), \end{aligned}$$

and

$$\begin{aligned} t_1(\lambda)&=t+\delta_0\cos\lambda,\\ t_2(\lambda)&=t-\delta_0\cos\lambda,\\ \Delta(\lambda)&=\Delta_0\sin\lambda. \end{aligned}$$

Take $t=1$, $\delta_0=0.2$, and $\Delta_0=0.4$ in hopping units. The complete ten-field audit is:

1. Crystal and orientation: a one-dimensional cell of length $a$, periodic bulk boundary conditions, $+x$ along increasing cell index, and the torus oriented by increasing $k$ followed by increasing $\lambda$.
2. Charge and state: $e>0$, one spinless electron in the lower band per cell, and a rigid neutralizing ionic background with no ionic winding.
3. Projector and gaps: the rank-one lower-band projector; the minimum band gap on the cycle is $0.8$ in the chosen units, attained near $k=\pi/a$.
4. Gauge and sewing: a periodic orbital gauge with the final $k$ link sewn to the initial fiber; forward overlaps are unitarized before multiplication.
5. Target: the unwrapped electronic polarization, the pump Chern number, and the electronic contribution to charge transported toward $+x$.
6. Path: $\lambda:0\to2\pi$ with the functions written above. The loop encircles the only nearby degeneracy at $k=\pi/a$, $\delta=0$, and $\Delta=0$, but never crosses it.
7. Branch: begin with any representative $P_0$ and follow the determinant Wilson phase continuously; the endpoint representative is $P_0-e$.
8. Symmetry: inversion is available at selected endpoints, but no point symmetry is assumed throughout the cycle or used to quantize the pump.
9. Boundary and probe: the invariant is computed in the periodic bulk; open-chain charge accumulation additionally needs an edge-occupation or reservoir protocol, while a current measurement integrates $J=dP/dt$.
10. Validation: projector-curvature quadrature gives $C=-0.9987$, $-0.9997$, and $-0.9999$ on $40^2$, $80^2$, and $160^2$ meshes. At the finest mesh the minimum overlap singular value exceeds $0.994$, the sewing residual is below $10^{-11}$, and the sampled minimum gap remains $0.8$. Finite drive rate, temperature, and system size remain separate errors.

The converged occupied-band result in the stated orientation is therefore

$$C_{k\lambda}=-1, \qquad Q_{\mathrm{el},+x}=-e.$$

The continuously followed electronic polarization falls by one quantum even though its principal-value endpoint equals its initial equivalence class. Reversing the cycle gives $C=+1$ and reverses the transported charge. Two explicitly occupied spin copies would pump $-2e$ for the displayed cycle only if both follow the same gapped path.

Now displace the loop so that it no longer encloses the degeneracy while keeping the gap open. It is contractible and gives $C=0$. By contrast, setting $\Delta=0$ while changing the sign of $\delta$ crosses the degeneracy. That path does not define a pump Chern number; a numerically rounded value obtained while missing the closing is a failed audit, not an approximate topological pump.

Worked Audit: Ferroelectric Switching and Measured Charge

Suppose a fixed-cell insulating path connects a centrosymmetric reference to a polar structure. Let the polarization quantum along the measurement direction be $P_q=0.80\ \mathrm{C\,m^{-2}}$, and choose the followed reference branch as $P_{\mathrm{ref}}=0$. At the polar endpoint a calculation reports

$$P_{\mathrm{ion}} = 0.46\ \mathrm{C\,m^{-2}}, \qquad P_{\mathrm{el}}^{\mathrm{principal}} = -0.73\ \mathrm{C\,m^{-2}}.$$

The accompanying ten-field audit records:

1. Crystal and orientation: a fixed three-dimensional primitive cell, positive cell volume, the origin at the reference inversion center, and $+z$ along the measured polar axis.
2. Charge and state: $e>0$, the listed core charges $Z_\kappa e$, twelve explicitly spinful occupied spinors, charge neutrality, and zero temperature.
3. Projector and gaps: a rank-twelve occupied projector at all 21 structural images; the minimum direct and indirect gaps are $1.20\ \mathrm{eV}$ and $0.92\ \mathrm{eV}$, respectively.
4. Gauge and sewing: reciprocal-sewn occupied-subspace Wilson strings along $k_z$, first on an $8\times8\times12$ mesh and then on a $12\times12\times18$ mesh, with no energy-sorted band tracking.
5. Target: $\Delta P_z$ from the centrosymmetric reference to the polar structure and the corresponding ideal electrode charge $A\Delta P_z$.
6. Path: a fixed-cell interpolation $\lambda:0\to1$ whose occupied rank and insulating gaps remain open at every image.
7. Branch: $P_{\mathrm{ref}}=0$ is followed continuously. The endpoint’s separately computed ionic and electronic representatives are the values above.
8. Symmetry: inversion quantizes the reference class; the selected branch is $0$ rather than $P_q/2$. The polar endpoint breaks inversion, so its polarization is not symmetry-quantized.
9. Boundary and probe: electrodes normal to $z$ integrate switching current; termination, screening, domains, and leakage belong to the forward model.
10. Validation: the refined mesh changes the final branch by less than $0.004\ \mathrm{C\,m^{-2}}$, the smallest overlap singular value is $0.985$, the sewing residual is below $10^{-10}$, and an independent current integration agrees within the same tolerance.

The principal endpoint total is $-0.27\ \mathrm{C\,m^{-2}}$. Intermediate gapped-path samples continue from the reference branch toward $+0.53\ \mathrm{C\,m^{-2}}$. Adding one quantum is therefore required:

$$P_{\mathrm{polar}}^{\mathrm{continuous}} = -0.27+0.80 = 0.53\ \mathrm{C\,m^{-2}}.$$

Thus $\Delta P=0.53\ \mathrm{C\,m^{-2}}$ on the stated path. That is a branch-resolved change from the chosen reference, not an absolute dipole assigned to one cell. If a proposed straight path closes the gap, choose and justify another insulating path or stop.

For electrode area $A$, ideal uniform switching transfers $Q_{\mathrm{electrode}}=A\Delta P$. A measured current trace must still be separated into polarization, ordinary leakage, linear capacitive background, domain-wall motion, incomplete switched volume, and electrode-screening contributions. Two-Dimensional Magnets and Ferroelectrics owns material-specific order and evidence; this page owns the branch and charge ledger used to interpret them.

Common Claim Failures

• Treating the dipole of one arbitrarily chosen periodic cell as an absolute bulk observable.
• Reporting only the electronic Berry phase while omitting ions, charge sign, spin multiplicity, origin, or orbital embedding.
• Summing separately labeled occupied-band phases through an internal crossing instead of using the full occupied projector or Wilson loop.
• Omitting reciprocal sewing from the last overlap link.
• Applying the insulating polarization formula to a metal, an indirectly overlapping filling, or a path whose gap closes.
• Choosing endpoint principal branches independently and calling their difference the transported charge.
• Calling a generic Zak phase topological or quantized without the symmetry and convention that enforce a discrete class.
• Inferring an absolute end charge from bulk polarization without termination, reconstruction, and screening data.
• Equating a finite-rate current close to one electron per cycle with a proven Chern pump before checking cycle closure, gap, filling, and rate convergence.
• Using the clean occupied-band torus in a disordered or interacting problem where a boundary-twist many-body construction is required.

Exercises

1. Translate a charge center by a lattice vector

In one dimension, an occupied electron center is shifted from $\bar x$ to $\bar x+a$ without changing the infinite crystal. Find the change in its Berry phase and electronic polarization.

Solution

The center relation gives

$$\Delta\gamma = \frac{2\pi}{a}a = 2\pi.$$

Therefore

$$\Delta P_{\mathrm{el}} = -\frac{e}{2\pi}\Delta\gamma = -e.$$

The two values represent the same bulk polarization class because the one-dimensional quantum is $e$.

2. Audit a large gauge transformation

For one occupied band, apply $|u_k\rangle\mapsto e^{-ikma}|u_k\rangle$, with integer $m$. Find the changes in $\mathcal A_k$, $\gamma$, and $P_{\mathrm{el}}$.

Solution

Here $\chi(k)=-kma$, so

$$\mathcal A_k' = \mathcal A_k-\partial_k\chi = \mathcal A_k+ma.$$

Across a Brillouin zone of length $2\pi/a$,

$$\gamma' = \gamma+2\pi m, \qquad P_{\mathrm{el}}' = P_{\mathrm{el}}-me.$$

This is a branch change, not a change of the periodic bulk state.

3. Check origin independence of a neutral cell

A one-dimensional cell contains charge $+e$ at $\tau$ and one occupied electron center at $\bar x$. Show that a common origin shift leaves the total polarization unchanged, while translating only the electron center by $a$ changes its representative by a quantum.

Solution

In the one-dimensional charge-per-boundary convention,

$$P = \frac{e\tau}{a} - \frac{e\bar x}{a}.$$

Under $\tau\mapsto\tau-x_0$ and $\bar x\mapsto\bar x-x_0$, the two origin shifts cancel because the cell is neutral. Under $\bar x\mapsto\bar x+a$ alone, $P\mapsto P-e$, which is the same equivalence class.

4. Unwrap a sewn multiband Wilson loop

A spinless one-dimensional pump has a rank-two occupied subspace. At five successive values of $\lambda$, a code multiplies the unitary polar factors of the forward overlap matrices in the explicit order $\mathcal W=\widetilde M_0\widetilde M_1\cdots\widetilde M_{N-1}$. The final link includes reciprocal sewing. It reports the principal phases

$$\operatorname{Arg}\det\mathcal W = 0, \quad -0.60\pi, \quad +0.80\pi, \quad +0.40\pi, \quad 0.$$

Use $\gamma=-\operatorname{Arg}\det\mathcal W$ and a nearest-neighbor lift on a sufficiently fine $\lambda$ mesh. Assume every resolved lifted increment has magnitude below $\pi$ and the winding is unchanged when the $\lambda$ mesh is doubled. Find $P_{\mathrm{el}}/e$, the pump Chern number, and the transported conventional charge for fixed ions. An independent two-dimensional $(k,\lambda)$ mesh audit gives

$$\begin{array}{c|ccc} N_k\times N_\lambda & 40^2 & 80^2 & 160^2 \\ \hline C_{k\lambda} & -0.992 & -0.998 & -0.9995 \\ \min\sigma(M_j) & 0.910 & 0.976 & 0.994 \\ \text{sewing residual} & 8\times10^{-4} & 5\times10^{-6} & 3\times10^{-8} \end{array}$$

Assess convergence, and explain what would fail if the final reciprocal-sewing link were omitted.

Solution

The principal values of $\gamma$ are

$$0, \quad +0.60\pi, \quad -0.80\pi, \quad -0.40\pi, \quad 0.$$

The stated nearest-neighbor and refinement tests select the lift

$$\gamma = 0, \quad 0.60\pi, \quad 1.20\pi, \quad 1.60\pi, \quad 2\pi.$$

Therefore

$$\frac{P_{\mathrm{el}}}{e} = -\frac{\gamma}{2\pi} = 0, \quad -0.30, \quad -0.60, \quad -0.80, \quad -1.00.$$

The cycle has $\Delta P_{\mathrm{el}}=-e$, so in the page’s $dk\wedge d\lambda$ orientation $C_{k\lambda}=-1$ and $Q_{\mathrm{el},+x}=-e$. Because the ions are fixed, this is also the total polarization change and transported conventional charge.

The Chern estimate approaches $-1$, the smallest raw-overlap singular value approaches one, and the sewing residual falls rapidly. The independently converged Chern audit selects the same winding as the fine-$\lambda$ branch lift. Together with an open gap, sub-$\pi$ resolved phase increments, and stability under both $k$- and $\lambda$-mesh refinement, those trends support convergence; rounding one coarse-mesh Chern estimate to an integer would not. Under an occupied-frame $U(2)$ change, a correctly sewn Wilson loop transforms by conjugation at the base point, so its determinant phase and eigenvalue set are invariant. Without the final sewing link, the product connects inequivalent endpoint frames, the gauge matrices do not close by conjugation, and its phase is not a valid Brillouin-zone holonomy.

5. Reverse a pump

A gapped spinless cycle has $C_{k\lambda}=+1$ in the convention of this page. What are the conventional charge and aggregate electron-center displacements? What changes when the cycle orientation is reversed?

Solution

The forward cycle transports $Q_{\mathrm{el},+x}=+e$. Conventional positive charge moving toward $+x$ is equivalent to the occupied electron-center set flowing by one lattice vector toward $-x$. Reversing the cycle changes the torus orientation, so $C=-1$, $Q_{\mathrm{el},+x}=-e$, and the center flow reverses.

6. Separate inversion quantization from an SSH edge claim

A student finds an electronic Zak phase $\pi$ for an inversion-symmetric SSH chain and concludes that the absolute polarization is $e/2$ and every termination carries charge $e/2$. Repair the claim.

Solution

Inversion constrains the total polarization class to $0$ or $e/2$ modulo $e$ after the inversion center, origin, cell, ionic charges, filling, and embedding are declared. The electronic Zak representative changes with unit cell and embedding. Boundary charge is fixed only modulo an integer charge and also depends on termination, edge occupation, reconstruction, and screening. The bulk result licenses a class and an interface-difference statement, not one absolute charge for every cut.

7. Infer an interface charge

Two one-dimensional insulators have representatives $P_1=0.10e$ and $P_2=0.60e$. The interface normal points from region 1 to region 2. Find the bound interface charge class and state what remains undetermined.

Solution

With the orientation used here,

$$Q_{\mathrm b} = P_1-P_2 = -0.50e \pmod e.$$

Equivalently the class can be represented by $+0.50e$. The integer part, microscopic distribution, and measured net charge require termination, interface reconstruction, defect, free-carrier, and screening information.

8. Choose the correct extension or stop

Audit four proposals: (a) compute band polarization at a metallic filling; (b) claim an exact pump from one fast cycle that crosses a gap closing; (c) evaluate a disordered interacting pump using a unique gapped ground state over boundary twist and cycle parameter; (d) follow one unique gapped open-chain ground state through a nontrivial cycle without discussing edge occupations.

Solution

(a) Stop: the occupied rank changes at the Fermi surface, so the insulating band-polarization formula is not defined for that ground state.

(b) Stop: the closing destroys the pump torus gap, and fast evolution adds nonadiabatic transitions. Near-integer transferred charge would be an empirical finite-protocol result, not the Chern theorem.

(c) Proceed with the many-body twist construction, provided the gap, charge conservation, orientation, size convergence, and numerical curvature are all validated. The relevant invariant is $C_{\mathrm{MB}}$, not a sum of band Chern numbers. It fixes the twist-averaged pumped charge; a fixed-twist result also requires demonstrated thermodynamic twist insensitivity.

(d) Repair the protocol: a nontrivial open pump generally has boundary spectral flow or tiny edge avoided crossings. Specify reservoirs, edge-state occupation, or diabatic passage through edge splittings, and compare with the periodic bulk invariant.

References

• M. V. Berry, “Quantal Phase Factors Accompanying Adiabatic Changes,” Proceedings of the Royal Society of London A 392, 45–57 (1984), doi:10.1098/rspa.1984.0023.
• R. D. King-Smith and D. Vanderbilt, “Theory of Polarization of Crystalline Solids,” Physical Review B 47, 1651–1654 (1993), doi:10.1103/PhysRevB.47.1651.
• Q. Niu and D. J. Thouless, “Quantised Adiabatic Charge Transport in the Presence of Substrate Disorder and Many-Body Interaction,” Journal of Physics A 17, 2453–2462 (1984), doi:10.1088/0305-4470/17/12/016.
• G. Ortiz and R. M. Martin, “Macroscopic Polarization as a Geometric Quantum Phase: Many-Body Formulation,” Physical Review B 49, 14202–14210 (1994), doi:10.1103/PhysRevB.49.14202.
• R. Resta, “Macroscopic Polarization in Crystalline Dielectrics: The Geometric Phase Approach,” Reviews of Modern Physics 66, 899–915 (1994), doi:10.1103/RevModPhys.66.899.
• R. Resta, “Quantum-Mechanical Position Operator in Extended Systems,” Physical Review Letters 80, 1800–1803 (1998), doi:10.1103/PhysRevLett.80.1800.
• M. J. Rice and E. J. Mele, “Elementary Excitations of a Linearly Conjugated Diatomic Polymer,” Physical Review Letters 49, 1455–1458 (1982), doi:10.1103/PhysRevLett.49.1455.
• D. J. Thouless, “Quantization of Particle Transport,” Physical Review B 27, 6083–6087 (1983), doi:10.1103/PhysRevB.27.6083.
• D. Vanderbilt, Berry Phases in Electronic Structure Theory, Cambridge University Press (2018), doi:10.1017/9781316662205.
• D. Vanderbilt and R. D. King-Smith, “Electric Polarization as a Bulk Quantity and Its Relation to Surface Charge,” Physical Review B 48, 4442–4455 (1993), doi:10.1103/PhysRevB.48.4442.
• H. Watanabe and M. Oshikawa, “Inequivalent Berry Phases for the Bulk Polarization,” Physical Review X 8, 021065 (2018), doi:10.1103/PhysRevX.8.021065.
• D. Xiao, M.-C. Chang, and Q. Niu, “Berry Phase Effects on Electronic Properties,” Reviews of Modern Physics 82, 1959–2007 (2010), doi:10.1103/RevModPhys.82.1959.
• J. Zak, “Berry’s Phase for Energy Bands in Solids,” Physical Review Letters 62, 2747–2750 (1989), doi:10.1103/PhysRevLett.62.2747.