---
title: Patch Euler Number in Kagome Systems
url: https://www.emergentmind.com/topics/patch-euler-number
type: topic
---

# Patch Euler Number in Kagome Systems

The patch Euler number is a local topological invariant defined for a pair of bands in a three-band, real-symmetric Hamiltonian with time-reversal symmetry \(T^2=+1\) when the corresponding global Euler number is ill-defined on the full Brillouin zone. In the kagome-lattice setting, it is attached to a subregion \(R\subset BZ\) on which the principal bands \(2\)–\(3\) are isolated from band \(1\), and it measures the net \(SO(2)\) vortex charge of a real two-band frame over that patch. Its central use is diagnostic: if the patch Euler number is non-zero, the enclosed principal Dirac points cannot annihilate under continuous deformations that keep \(R\) free of additional degeneracies with the third band [2507.19238].

## 1. Topological setting in the three-band kagome model

In the “base II” representation, the kagome tight-binding Bloch Hamiltonian is a \(3\times 3\) real-symmetric matrix,
$$
H(k)=-
\begin{bmatrix}
E_A & t_{AB}\cos(k\cdot\delta_1) & t_{AC}\cos(k\cdot\delta_2)\\
t_{AB}\cos(k\cdot\delta_1) & E_B & t_{BC}\cos(k\cdot\delta_3)\\
t_{AC}\cos(k\cdot\delta_2) & t_{BC}\cos(k\cdot\delta_3) & E_C
\end{bmatrix},
$$
with \(k=(k_x,k_y)\), nearest-neighbor vectors \(\delta_1,\delta_2,\delta_3\), and tunable parameters \(E_{A,B,C}\) and \(t_{AB,AC,BC}\). Because \(H(k)\) is real, time-reversal symmetry takes the form \(H(k)=H(-k)\), and the eigenvectors can be chosen real and orthonormal at each \(k\) [2507.19238].

The bands are labeled by increasing energy,
$$
E_1(k)\le E_2(k)\le E_3(k).
$$
The model generically features Dirac points between bands \(2\)–\(3\), termed principal Dirac points, and also between bands \(1\)–\(2\), termed adjacent Dirac points. The motion of the adjacent Dirac points is not ancillary: it determines whether the direct sum of bands \(2\)–\(3\) defines a globally meaningful rank-2 real bundle and, consequently, whether a pair of principal Dirac points can be annihilated.

The ordinary Euler number of two bands is globally defined only when those bands remain isolated from the third band over the entire Brillouin zone. In the kagome model, parameter regimes exist where bands \(2\)–\(3\) touch band \(1\) at adjacent Dirac points somewhere in \(BZ\). In that situation the global two-band bundle over the full \(2\)-torus does not exist, so the global Euler number is ill-defined. The patch Euler number is the local replacement for that missing invariant.

## 2. Real two-band geometry on a subregion

To define the invariant, one fixes a subregion \(R\subset BZ\) such that the principal bands \(2\)–\(3\) are isolated from band \(1\) for all \(k\in R\). Equivalently, there is no degeneracy between \(\{2,3\}\) and band \(1\) anywhere in \(R\). The direct sum of the eigenspaces of bands \(2\) and \(3\) then defines a rank-2 real subbundle \(E\to R\) [2507.19238].

On this patch one chooses a smooth, real, orthonormal frame \(\{u_1(k),u_2(k)\}\) spanning \(E_k\). In that frame, the real Berry connection is
$$
A_\mu^{ab}(k)=\langle u_a(k)\mid \partial_\mu u_b(k)\rangle,
\qquad \mu\in\{k_x,k_y\},\ a,b\in\{1,2\},
$$
and the \(SO(2)\) Berry, or Euler, connection is its off-diagonal component,
$$
\omega_\mu(k)=A_\mu^{12}(k)=-A_\mu^{21}(k).
$$
In differential-form notation, \(\omega=\omega_x\,dk_x+\omega_y\,dk_y\) is the only independent component of the skew-symmetric connection matrix.

The associated \(SO(2)\) curvature is
$$
f(k)=\partial_{k_x}\omega_{k_y}(k)-\partial_{k_y}\omega_{k_x}(k),
$$
or, equivalently, \(f=d\omega\). Under a local \(O(2)\) frame rotation \(u\to X(k)u\) with \(X(k)\) a local rotation by \(\theta(k)\), the connection transforms as \(\omega\to \omega+d\theta\), while the curvature remains gauge invariant, \(f\to f\). This local \(SO(2)\) geometry is the geometric input from which the patch Euler number is extracted.

The restriction to a patch is not merely technical. On \(R\), one can choose a smooth, real, orthonormal frame except possibly along gauge branch cuts, or Dirac strings, that connect Dirac points. The curvature remains well defined and gauge invariant on \(R\), and if \(\omega\) is smooth on \(\partial R\), Stokes’s theorem converts the area integral into a boundary integral.

## 3. Definition and computation

For a smooth gauge on \(R\) and its boundary, the patch Euler number is
$$
e_{\mathrm{patch}}(R)=\frac{1}{2\pi}\int_R f(\mathbf{k})\,d^2k.
$$
When the gauge is smooth on \(\partial R\), Stokes’s theorem gives the equivalent boundary formula
$$
e_{\mathrm{patch}}(R)=\frac{1}{2\pi}\oint_{\partial R}\omega(\mathbf{k})\cdot d\mathbf{k}.
$$
If unavoidable gauge discontinuities occur along \(\partial R\), for example because a Dirac string crosses the boundary, the gauge-covariant expression is
$$
e_{\mathrm{patch}}(R)=\frac{1}{2\pi}\left[\int_R f(\mathbf{k})\,d^2k-\oint_{\partial R}\omega(\mathbf{k})\cdot d\mathbf{k}\right].
$$
In the formulation used in the paper, this subtraction produces a gauge-invariant integer or half-integer result. When a globally smooth \(\omega\) exists on \(R\cup \partial R\), the boundary term vanishes and the simple area integral is sufficient [2507.19238].

A practical numerical workflow is explicitly given:

1. Diagonalize \(H(k)\) on a mesh in \(R\) to obtain real, normalized states for the principal bands.
2. Enforce a continuous gauge on \(R\) as much as possible and record Dirac strings if they are unavoidable.
3. Compute \(A_\mu^{ab}(k)=\langle u_a\mid \partial_\mu u_b\rangle\) numerically, for example by finite differences.
4. Form \(\omega_\mu=A_\mu^{12}\) and \(f=\partial_{k_x}\omega_{k_y}-\partial_{k_y}\omega_{k_x}\).
5. Integrate \(f\) over \(R\), and if the gauge is not smooth along \(\partial R\), subtract the boundary contribution.

The sign of \(e_{\mathrm{patch}}\) is gauge dependent, but its vanishing or non-vanishing in a fixed setup is the physically relevant datum for annihilation. The formalism also admits half-integer values in certain conventions when adjacent Dirac strings pierce the patch and boundary corrections are essential.

## 4. Dirac points, vortex charges, and annihilation obstruction

Locally, the \(SO(2)\) frame can be parameterized by an angle \(\theta(k)\), so that \(\omega=d\theta\) and \(f=d\omega\). Around an isolated principal Dirac point, \(\theta\) winds by \(\pm 2\pi\) along a small loop \(C\), yielding
$$
\frac{1}{2\pi}\oint_C \omega\cdot dk
=
\frac{1}{2\pi}\oint_C d\theta
=
w\in\mathbb{Z}.
$$
Each principal Dirac point therefore contributes an integer vortex charge \(q_i=\pm 1\). The sign convention is explicit: counterclockwise increase of \(\theta\) gives \(+1\), while clockwise gives \(-1\) [2507.19238].

When \(R\) encloses only principal Dirac points and remains free of third-band degeneracies, the patch Euler number reduces to the charge sum
$$
e_{\mathrm{patch}}(R)=\sum_{i\in R} q_i.
$$
Adjacent Dirac strings that cross \(R\) or \(\partial R\) can reverse one local contribution by flipping the gauge; the boundary correction in the definition of \(e_{\mathrm{patch}}\) precisely accounts for this. In that sense the invariant measures the net \(SO(2)\) vortex charge of the two-band frame on the chosen patch.

The obstruction criterion is direct. If
$$
e_{\mathrm{patch}}(R)\neq 0,
$$
then the enclosed principal Dirac points cannot be pairwise annihilated by any continuous parameter deformation that keeps \(R\) free of additional degeneracies with the third band. If two principal Dirac points were to annihilate and open a gap between bands \(2\) and \(3\) over \(R\), the rank-2 bundle over \(R\) would have to deform into a trivial direct sum of two rank-1 trivial bundles, a triviality guaranteed by \(T^2=+1\). A non-zero patch Euler number signals a twisting that cannot be unwound without crossing a degeneracy with band \(1\).

## 5. Quaternionic charges, non-abelian braiding, and parity

The paper gives a second description in homotopy-theoretic terms. The space of real orthonormal frames of the three-band system, modulo sign flips of individual eigenvectors, is homotopy equivalent to \(SU(2)/Q_8\), where
$$
Q_8=\{\pm 1,\pm i,\pm j,\pm k\}.
$$
Loops in \(k\)-space that encircle Dirac points are therefore assigned quaternionic charges. The identifications used are [2507.19238]:

- **Principal Dirac point**: \(\pm i\)
- **Adjacent Dirac point**: \(\pm k\)
- **One principal and one adjacent in the loop**: \(\pm j\)
- **Two same-type non-annihilable points in the loop**: \(-1\)
- **Trivial loop**: \(1\)

In this language, non-abelian braiding in \(k\)-space acts by conjugation,
$$
q_i \to s q_i s^{-1},
$$
where \(s\) is the charge of the braid path around the other Dirac point. Braiding a principal Dirac point around an adjacent one can therefore change the sign or type of the principal charge. This explains why a configuration that was previously non-annihilable may become annihilable after adjacent strings cross the patch.

The annihilation criterion becomes an ordered-product condition: for a collection of Dirac points inside a loop \(L\), the ordered product of their charges along \(L\) must equal \(1\). If the product is \(-1\), corresponding to two same-type, same-orientation points, annihilation is obstructed. The gauge-sensitive flips of \(e_{\mathrm{patch}}\) induced by adjacent Dirac strings are thus the local, \(SO(2)\)-frame manifestation of a non-abelian quaternionic algebra.

The same section relates the Euler class to the second Stiefel–Whitney class. For an oriented rank-2 real bundle,
$$
w_2 \equiv e \pmod 2.
$$
Over a patch, odd \(e_{\mathrm{patch}}\) gives a non-zero \(w_2\), signaling a parity obstruction consistent with the \(-1\) product in the quaternionic picture. Even \(e_{\mathrm{patch}}\), by contrast, has \(w_2=0\), which allows trivialization and hence annihilation if no other obstruction is present.

## 6. Kagome tuning protocols and photonic realization

The kagome tight-binding model admits concrete parameter protocols in which the patch Euler number changes as adjacent Dirac points are created and their strings reconfigure. In the simplest nearest-neighbor limit, equal hoppings \(t>0\) and zero on-site energies give a real-symmetric, \(T\)-invariant Bloch Hamiltonian \(H_0(k)\). Deforming \(E_{A,B,C}\) and \(t_{AB,AC,BC}\) then moves both principal and adjacent Dirac points [2507.19238].

Two protocols are described explicitly:

| Protocol | Tuning | Outcome |
|---|---|---|
| Six-panel sequence (Fig. 6) | Start at \(E_A=1\), \(t_{AB}=1\), \(t_{AC}=-1\), \(E_B=0\), \(t_{BC}=-0.3\), \(E_C=0\); vary \(t_{BC}\to 0\to 0.2\to 0.7\), then \(E_A\to -0.3\to -1\) | Principal Dirac points first merge and bounce with \(e_{\mathrm{patch}}=-1\); after adjacent Dirac points appear and an adjacent Dirac string crosses \(R\), \(e_{\mathrm{patch}}=0\) and annihilation becomes allowed |
| Alternative photonic-friendly protocol | Tune only \(E_A\) from \(0.5\) to \(-1.5\) | Principal Dirac points merge and initially bounce, then annihilate at the \(BZ\) edge thanks to periodicity; \(e_{\mathrm{patch}}\) changes as adjacent strings reconfigure |

The figure sequence gives the clearest explicit demonstration. In Fig. 6(a\(\to\)c), principal Dirac points merge but bounce and no gap opens; the chosen rectangular patch \(R\) has \(e_{\mathrm{patch}}=-1\). In Fig. 6(d\(\to\)f), after tuning \(E_A\), adjacent Dirac points appear; an adjacent Dirac string crosses \(R\), one principal charge is reversed, \(e_{\mathrm{patch}}=0\), and the principal Dirac points annihilate.

The paper also argues that the same deformation can be implemented in realistic photonic systems. In photonic kagome lattices, specifically microcavity pillar arrays or EIT-induced lattices, \(E_A\) can be tuned in situ by optical pumping and effective hoppings can be controlled through geometry. Beyond tight-binding numerics, the authors solve a \(2\)D stationary Schrödinger equation with a structured potential \(U_0(x,y)+U_1(x,y)\). The resulting dispersions show a double Dirac point at \(\Gamma\) in the unperturbed case and separated Dirac points when \(E_A<0\), reproducing the annihilation sequence. In that setting, selecting a momentum-space patch \(R\) and computing \(e_{\mathrm{patch}}\) from the photonic bands provides a direct topological readout of whether annihilation is obstructed or allowed.

## 7. Terminological ambiguity and non-equivalent usages

The expression “patch Euler number” is specific in the kagome-band-topology setting, but the phrase is not standard in the vortex-patch literature on the \(2\)D incompressible Euler equations. One source states this explicitly and notes that the relevant object there is the Euler equation for vortex patches rather than the topological Euler characteristic or Euler number [1703.09674].

This ambiguity matters because several fluid-dynamical papers involve both “patch” and “Euler” while referring to entirely different quantities. In one interpretive usage, the phrase is attached to the patch problem for \(2\)D Euler as a quantitative measure of time analyticity, namely the radius \(R\) in bounds of the form \(C R^{-m} m!\) for derivatives of the Lagrangian flow [1907.13407]. Other works on Euler patches study winding numbers for particle trajectories in disk-like vortex patches [2008.05085], or global regularity of patch solutions for Loglog–Euler type active scalar equations [2510.18759]. These are mathematically separate from the momentum-space invariant \(e_{\mathrm{patch}}(R)\).

In the topological usage treated here, the patch Euler number is therefore not a fluid-dynamical observable attached to a moving characteristic-function domain. It is a local Euler invariant of a real rank-2 band bundle on a chosen region of the Brillouin zone, designed precisely for situations in which the full-zone Euler number cannot be defined because the principal two-band sector is not globally isolated from the third band.

Source: https://www.emergentmind.com/topics/patch-euler-number