---
title: Hybrid Wannier Charge Centers
url: https://www.emergentmind.com/topics/hybrid-wannier-charge-center-wcc
type: topic
---

# Hybrid Wannier Charge Centers

Hybrid Wannier Charge Centers (WCCs) are the real-space centers of hybrid Wannier functions, namely states that are Wannier-localized along one crystallographic direction while remaining Bloch-like in the others. In contemporary band theory, they are equivalently the eigenphases of a non-Abelian Wilson loop, or Berry phases along a closed reciprocal-space cycle, expressed as positions modulo a lattice vector. Their evolution over a projected Brillouin zone provides a bulk representation of electric polarization, Chern and \(\mathbb Z_2\) topology, crystalline topological indices, and, in several extensions, axion response, higher-order topology, chemical bonding, and optical geometry [1312.6940][1610.08983][1912.11887].

## 1. Formal definition and geometric meaning

Choosing \(z\) as the Wannierized direction, hybrid Wannier functions may be written as
\[
|W_{n l_z}(k_x,k_y)\rangle = \frac{1}{2\pi}\int dk_z\, e^{i\mathbf{k}\cdot(\mathbf{r}-l_z c\,\hat z)} |u_{n\mathbf{k}}\rangle ,
\]
or, equivalently in the notation used for 3D insulators,
\[
|h_{ln\mathbf{k}}\rangle=\frac{1}{2\pi} \int^{\pi/c}_{-\pi/c} dk_z \, e^{-ik_z l c} \ket{\psi_{n\mathbf{k}k_z}} .
\]
The corresponding WCC is the position expectation value
\[
\bar z_n(k_x,k_y)=\langle W_{n0}|\hat z|W_{n0}\rangle ,
\]
or, equivalently,
\[
z_{ln}(\mathbf{k})=\bra{h_{ln\mathbf{k}}z\ket{h_{ln\mathbf{k}} ,
\qquad
z_{ln}(\mathbf{k})=z_{0n}(\mathbf{k})+lc .
\]
In this representation, the WCCs are the eigenvalues of the projected position operator \(P z P\), and their dependence on the transverse momentum \((k_x,k_y)\) defines “Wannier bands” or “WCC sheets” over the projected 2D Brillouin zone [1312.6940][1912.11887][2101.10798].

For an isolated band, the WCC is the Berry phase along the Wannierized direction written as a real-space coordinate:
\[
\bar z(k_\perp)=\frac{c}{2\pi}\phi(k_\perp),
\qquad
\phi=-\mathrm{Im}\,\ln \prod \langle u_k|u_{k+b}\rangle .
\]
In the multiband case, the same object appears as a Wilson-loop eigenphase. If the Wilson-loop eigenvalues are \(e^{i\varphi_n(\mathbf{k})}\), then
\[
z_n(\mathbf{k}) = \frac{c}{2\pi}\,\varphi_n(\mathbf{k}) \quad \mathrm{mod}\ c .
\]
This equivalence underlies the standard identification of hybrid WCCs with Wilson-loop spectra [1312.6940][1912.11887].

The physical interpretation is twofold. First, the sum of WCCs gives a polarization-like quantity, so WCC flow is a momentum-resolved form of Berry-phase polarization. Second, because WCCs are defined modulo a lattice vector, only their global evolution, winding, and connectivity are topologically meaningful. This is why WCC plots are typically read as spectral-flow diagrams rather than as absolute positions [1610.08983][2303.16306].

## 2. Wilson loops, overlap matrices, and numerical construction

In practical calculations, hybrid WCCs are extracted from overlap matrices between neighboring \(k\)-points on a closed loop. For a chosen band subspace,
\[
S_{mn}(\mathbf{k}_i,\mathbf{k}_j) = \langle u^{(m)}_{\mathbf{k}_i} \mid u^{(n)}_{\mathbf{k}_j} \rangle ,
\]
and the cumulative loop matrix is
\[
M = \prod_{i=1}^{J-1} S(\mathbf{k}_i,\mathbf{k}_{i+1}) .
\]
Its eigenvalues \(\lambda\) encode the Berry phases, and the hybrid WCCs are
\[
w_{\{\alpha\}} = \frac{\arg(\lambda)}{2\pi} \mod 1 .
\]
This is the Wilson-loop/HWCC construction implemented directly from first-principles Bloch states in the all-electron WIEN2k workflow, where no maximally localized Wannier-function construction is needed [2303.16306].

A closely related formulation used in Z2Pack proceeds by parallel transport. One computes overlap matrices
\[
M_{mn}^{(k_i,k_{i+1})} = \langle u_{m,k_i}|u_{n,k_{i+1}}\rangle ,
\]
performs a singular value decomposition \(M=V\Sigma W^\dagger\), forms the unitary transport step \(WV^\dagger\), and multiplies these along the loop. If the resulting non-Abelian phase matrix has eigenvalues \(\lambda_i\), the normalized WCCs are
\[
\bar x_i=-\frac{\arg(\lambda_i)}{2\pi},
\]
up to orientation convention [1610.08983].

These constructions operate on the occupied subspace as a whole. The gauge freedom is handled implicitly by working with overlap matrices and Wilson loops rather than by differentiating phases of individual bands. This makes the method numerically stable and naturally suited to composite occupied manifolds, including all-electron, tight-binding, and \(k\cdot p\) calculations [1610.08983][2303.16306].

Several practical caveats are standard. The phases are defined modulo \(2\pi\), so WCCs are defined modulo a lattice vector or modulo \(1\) in fractional coordinates. The physically meaningful quantity is therefore the continuous evolution of the WCC spectrum as the transverse momentum varies. The WIEN2k implementation also emphasizes that there is no exact one-to-one mapping between original Bloch-band indices and WCC indices, because the WCCs arise from diagonalizing the Wilson-loop matrix for the entire occupied subspace [2303.16306].

## 3. Topological diagnostics from WCC flow

The principal use of hybrid WCCs is topological diagnosis through spectral flow. For Chern insulators, the Chern number is the net change of the summed WCCs over a pumping cycle:
\[
C= \frac{1}{a_x} \left( \sum_n \bar x_n(k_y=2\pi)-\sum_n \bar x_n(k_y=0) \right),
\]
or equivalently
\[
C=\frac{1}{ea}\left(P_e^h(2\pi)-P_e^h(0)\right).
\]
This is the statement that the Chern number is a quantized charge pump encoded in WCC winding [1610.08983].

For time-reversal-invariant insulators, the relevant invariant is the parity of WCC crossings over half the Brillouin zone. In practice one follows the largest gap in the WCC spectrum and counts whether a Wannier band is crossed an odd or even number of times. The WIEN2k implementation summarizes the criterion as follows: if “one Wannier band is crossed when following the largest gap in the HWCC spectrum,” the phase is nontrivial. This is the Wilson-loop version of partner switching or adiabatic Thouless charge pumping [2303.16306][1610.08983].

In three dimensions, WCC sheets play the same diagnostic role as surface-state dispersions, but are obtained from bulk calculations alone. The “Wannier Center Sheets in Topological Insulators” framework shows that trivial, weak, strong, and crystalline topological insulators can be distinguished by the connectivity of WCC sheets over projected 2D Brillouin zones. In first-principles examples, Sb\(_2\)Se\(_3\) exhibits topologically trivial WCC sheets, KHgSb exhibits weak-topological behavior, and Bi\(_2\)Se\(_3\) exhibits strong-topological behavior [1312.6940].

The same logic extends to semimetals when the chosen manifold remains gapped. Z2Pack uses WCC winding on a small sphere to determine the chirality of a Weyl node, and can detect Weyl or Dirac structure from changes of Chern or \(\mathbb Z_2\) invariants between planes or curved surfaces in the Brillouin zone [1610.08983]. In layered organic conductors, WCC sheets localized along \(y\) and plotted over \((k_x,k_z)\) acquire especially clear geometric forms: flat but discontinuously cut sheets for nodal-line semimetals, oppositely winding double spirals around Dirac points, same-sense spirals around Weyl points, and partner switching with \((\nu_0;\nu_x,\nu_y,\nu_z)=(0;0,0,1)\) for the weak topological-insulator phase [2508.07164].

## 4. Crystalline, axion, and higher-order extensions

Mirror symmetry admits a symmetry-resolved WCC analysis in which the occupied space is split into mirror eigenspaces on a mirror-invariant plane. In the hybrid-Wannier representation, the mirror Chern number can be determined from the winding numbers of touching points between Wannier bands on mirror-pinned planes and from the Chern numbers of flat Wannier bands pinned to those planes. For type-1 mirrors, the mirror Chern numbers on the \(G\) and \(X\) planes are
\[
\mu_G=\frac{1}{2}\left(p_A\,\overline C_A+W_A\right)+ \frac{1}{2}\left(p_B\,\overline C_B+W_B\right),
\]
\[
\mu_X=\frac{1}{2}\left(p_A\,\overline C_A+W_A\right)- \frac{1}{2}\left(p_B\,\overline C_B+W_B\right),
\]
while for a type-2 mirror
\[
\mu = p_A\,\overline C_A + W_A .
\]
This formalism was illustrated for SnTe in both monolayer and bulk forms [2101.10798].

Axion or Chern-Simons topology can likewise be reformulated in WCC language. In the hybrid-Wannier representation, the axion coupling \(\theta\) admits a decomposition
\[
\theta=\theta_{z\Omega}+\theta_{\Delta xy},
\]
with
\[
\theta_{z\Omega}=-\frac{1}{c}\int d^2k\sum_n z_n \, \Omega_{xy,0n,0n} .
\]
For axion-quantizing symmetries, the WCC sheet geometry determines when \(\theta=0\) or \(\pi\), and connected or disconnected Wannier spectra can encode the same axion-\(\mathbb Z_2\) index depending on whether the protecting symmetry is \(z\)-preserving or \(z\)-reversing [1912.11887]. In the dynamical setting of Chern-Simons axion pumping, a nonzero second Chern number requires sheet-touching events at which \(2\pi\) quanta of Berry curvature are transferred from one WCC sheet to another, yielding \(\Delta\theta=2\pi\) over the cycle [1411.1753].

Higher-order topology introduces an additional layer. In the generalized separable BBH model, the hybrid WCCs are exactly the Wilson-loop eigenphases \(e^{i2\pi \nu_\mu^j(k_{\bar\mu})}\), and the nested-Wilson-loop structure yields exact Wannier-sector polarizations
\[
p_{\mu}^{j}=\frac{\mathcal{N}_\mu}{2}\,\text{mod}\,1,
\]
and quadrupole moment
\[
q_{xy}=\frac{\mathcal{N}_x\mathcal{N}_y}{2}\,\text{mod}\,1 .
\]
This establishes a \(\mathbb{Z}\times\mathbb{Z}\) Wannier topology controlled by the winding numbers of the constituent one-dimensional chiral chains [2208.05058].

At the same time, ordinary hybrid WCC flow is not universally sufficient for higher-order phases. In the strained Kane–Mele model, WCC evolution clearly diagnoses the quantum spin Hall phase for \(t_1/t<2\), but becomes topologically trivial in the second-order topological-insulator phase for \(t_1/t>2\). The paper explicitly concludes that ordinary WCC evolution is “incapable of capturing any essence of the second order topological phase,” and instead uses spin-resolved bulk polarization to identify the obstructed atomic-insulator character of the SOTI phase [2311.16011].

## 5. Generalizations beyond standard WCC flow

A notable three-dimensional extension is the “joint Wannier center” construction for charge partitioning in matter. Starting from the same non-Abelian Berry-phase matrices that define hybrid WCCs along one direction,
\[
X^\alpha(k_\beta)= \log \prod_{k_\alpha} \underline{\underline{O}^\alpha(k_\alpha,k_\beta)},
\]
the method considers the three matrices \(X^1(k)\), \(X^2(k)\), and \(X^3(k)\) associated with three linearly independent reciprocal directions and performs maximal joint diagonalization. The resulting expectation values define real-space points \(\vec d_i(k)\), called joint Wannier centers. The paper explicitly states that these are not standard hybrid WCCs, but a three-directional extension built from hybrid-Wannier data and used for chemical-bonding analysis [2407.17202].

Closely related Berry-center ideas also enter nonlinear optics. In quantized formal-polarization crystals, the relevant objects are band-resolved Berry-phase polarizations,
\[
\mathbf P_n=-e\int[d\mathbf k]\, \mathbf A_{nn},
\]
interpreted as Wannier centers of the initial and final states. The Brillouin-zone average of the optical shift vector obeys
\[
-e\int[d\mathbf{k}]\,\mathbf{R}_{mn} = \mathbf{Q}W_{mn}+ \mathbf{P}_{m} - \mathbf{P}_{n}.
\]
This identifies the average shift vector with a difference of Berry-phase centers modulo a polarization quantum and winding term. In this setting the paper does not explicitly use the term “hybrid WCC,” but it uses the closely related language of Wannier centers, fractional Wannier-center sectors, and Berry-phase polarization differences between optically connected bands [2606.16108].

A one-dimensional experimental counterpart is Wannier-center spectroscopy in artificial lattices on Cs/InAs(111)A. There the band-resolved local density of states is integrated over energy,
\[
p_x = \frac{ \int_{E_n^\text{min}}^{E_n^\text{max}} \rho(x, E) dE}{\int_\text{All bands} \rho(x, E) dE},
\]
to reconstruct the real-space Wannier-center location. In this purely one-dimensional setting there is no transverse momentum and hence no nontrivial hybrid-WCC flow, but the method directly measures the 1D analogue of the WCC position and uses it to distinguish trivial from boundary-obstructed topological phases [2407.14465].

## 6. Implementations, representative workflows, and limitations

The most widely used general-purpose software framework in this area is Z2Pack, which implements WCC tracking for \(k\cdot p\) models, tight-binding models, and ab initio calculations. Its workflow is based on choosing an appropriate closed manifold in the Brillouin zone, computing Wilson loops along a family of closed paths, and extracting invariants from WCC winding, parity of crossings, or symmetry-resolved subspace flow [1610.08983]. A complementary all-electron route is provided by the WIEN2k implementation built around `wcc.py` in BerryPI, which constructs `case.klist`, computes overlaps via Wien2wannier, stores them in `case.mmn`, and consolidates the resulting WCC spectrum into `wcc.csv` [2303.16306].

The method is robust, but its scope is not unlimited. WCCs are naturally defined modulo a lattice vector, so branch choices matter and only continuous evolution is physically meaningful [2303.16306]. Dense meshes are often required when curves approach closely or when largest-gap tracking becomes ambiguous; the WIEN2k implementation recommends refining both the loop mesh and the evolution-direction mesh if the flow appears jagged or unclear [2303.16306]. Z2Pack likewise emphasizes convergence control, especially near small direct gaps or near topological phase transitions [1610.08983].

Several conceptual limitations are also standard. Individual WCC branches are properties of the chosen occupied subspace and should not be interpreted as trajectories of individual Bloch bands, because the Wilson-loop eigenvectors generally mix several bands [2303.16306]. Standard hybrid WCCs are exact only for localization along one direction; three-directional constructions such as joint Wannier centers sacrifice exact one-directional localization in exchange for joint localization criteria [2407.17202]. Ordinary WCC flow can cleanly diagnose first-order \(\mathbb Z_2\) topology, but may fail for higher-order phases, as in the strained Kane–Mele example [2311.16011].

Taken together, these developments place hybrid WCCs at the center of a broad geometric framework. In their canonical form they are Wilson-loop eigenphases interpreted as real-space charge centers; in application they act as bulk diagnostics of topology, symmetry, and adiabatic transport; and in more recent extensions they serve as building blocks for chemical partitioning, experimental charge-center spectroscopy, and Berry-geometric descriptions of optical response [1312.6940][1610.08983][2407.17202].

Source: https://www.emergentmind.com/topics/hybrid-wannier-charge-center-wcc