Haldane model ($H(\mathbf{k})=d_0(\mathbf{k})I+\mathbf{d}(\mathbf{k})\cdot \mathbf{\sigma }$)

F. Duncan M. Haldane's 1988 tight-binding model adds complex second-neighbour hopping and opposite sublattice offsets to a honeycomb lattice. Directed phases break time-reversal symmetry while each primitive cell has zero net magnetic flux. Smaller loops can enclose magnetic flux; ordinary Bloch states remain available.^[ref]

Extend the nearest-neighbour model of graphene, ignoring physical spin and interactions. With bond length $a$, $t_1>0$, $t_2\ge 0$ and ordered sublattice basis $(A,B)$, chooseThe following Bloch Hamiltonian defines $\varphi$; reversing directed second-neighbour hopping sends $\varphi \mapsto -\varphi$:The Pauli matrices act on sublattice amplitudes, as in sublattice pseudospin in graphene. The identity matrix term $d_{0}I$ shifts both energy eigenvalues equally at each momentum; it changes neither the eigenvectors nor the direct separation $E_{+}-E_{-}$. A genuine band insulator additionally needs ${\min }_{\mathbf{k}}{E}_{+}>{\max }_{\mathbf{k}}{E}_{-}$. We use $t_2/t_1=0.08$ in the figures, within the original paper's sufficient nonoverlap restriction $|t_2/t_1|<1/3$.^[ref]

The only zeros of $f$, modulo the reciprocal lattice, are the two Dirac valleys in graphene,At these momenta ${\sum }_j\sin ({\mathbf{K}}_{\tau }\cdot {\mathbf{b}}_j)=\tau 3\sqrt{3}/2$. Thus the two Dirac mass coefficients, measured in energy units, are $m_{\tau }=M-\tau \mu$, with $\mu =3\sqrt{3}{t}_2\sin \varphi$. For $\mathbf{k}={\mathbf{K}}_{\tau }+\mathbf{q}$,The two Dirac cones have opposite orientation, encoded by $\tau$. A scalar sublattice offset $M$ gives equal-sign Dirac masses; imaginary second-neighbour hopping at $M=0$ gives opposite-sign Dirac masses. Both mechanisms open valley gapsAn observed band gap alone therefore cannot establish a Chern insulator. In particular, the gapless nearest-neighbour model of graphene and its sublattice-mass deformation should not be confused with the topologically nontrivial Haldane model.

A band gap closes when all three components of $\mathbf{d}$ vanish. Since $d_x=d_y=0$ only at the two valleys, the phase boundaries are exactly $M=\pm 3\sqrt{3}{t}_2\sin \varphi$. Usually one valley becomes massless; at $M=0$, $\sin \varphi =0$, both do. The occupied-band Chern number is undefined at every such touching, even if a numerical momentum grid misses its location.

For $\mathbf{d}\ne 0$, put $\hat{\mathbf{d}}=\mathbf{d}/|\mathbf{d}|$. The lower-band orthogonal projector is $P_{-}=(I-\hat{\mathbf{d}}\cdot \mathbf{\sigma })/2$. Applying the Pauli matrix multiplication law to $i\mathrm{Tr}(P_{-}[{\partial }_{k_x}{P}_{-},{\partial }_{k_y}{P}_{-}])$ givesThus the Chern number is the signed number of times the Brillouin zone covers a two-sphere under $\hat{\mathbf{d}}$. The lower energy eigenstate's actual Bloch vector is $-\hat{\mathbf{d}}$. The Berry curvature has the displayed positive sign for this lower band and the stated Berry connection; switching to the upper band changes its sign.

A useful local calculation exposes what each valley contributes. Write $p=\hslash v|\mathbf{q}|$. The linear Dirac Hamiltonian givesEach unbounded continuum cone gives a half-integer in this calculation. It is not a separately realizable half-integer Chern number of an isolated lattice energy band. The full honeycomb lattice supplies both valleys and the large-momentum completion. Starting with a trivial large-$|M|$ state, each single-valley mass crossing changes the lattice integer by the corresponding unit. There are no other band touchings. This establishes the exact lattice resultfor $t_2>0$, excluding the boundaries. This sign belongs to our explicit bond vectors, orientation and occupied-band convention.

For an original worked example, choose $t_1=1$, $t_2=0.08$, $\varphi =\pi /2$ and $M=0$. Then $\mu \simeq 0.415692$, $(m_{+},m_{-})=(-\mu ,+\mu )$, both valley gaps are $2\mu \simeq 0.831384$, and $C=-1$. Choosing instead $t_2=0$, $M=\mu$ gives exactly the same valley-gap magnitudes but $C=0$. The first state's two valley contributions reinforce; the second state's cancel.

An independent eigenvector-overlap calculation on periodic $31\times 31$, $61\times 61$ and $121\times 121$ momentum grids gives $C=-1$ for the first example and $C=0$ for the second. Reversing $\varphi$ in the first example gives $C=+1$. Direct integration of the analytic lattice Berry curvature agrees. For the overlap calculation one can remove the embedding phase from $f$, replacing it by $1+e^{i\mathbf{k}\cdot {\mathbf{a}}_1}+e^{i\mathbf{k}\cdot {\mathbf{a}}_2}$; this makes the cell-basis Hamiltonian periodic and prevents a false boundary-seam contribution. The local Berry curvature can depend on this basis convention; the integrated Chern number agrees. All mass-zero parameters are excluded before evaluating any integer.

Gregor Jotzu and collaborators realized an artificial Haldane model in 2014. They loaded spin-polarized potassium-40 fermions into a laser-made honeycomb lattice. Piezoelectric mirrors shook the optical lattice along two axes; their relative phase controlled time-reversal symmetry breaking. Optical phase control adjusted the sublattice offset.^[ref]

A magnetic field gradient drove Bloch oscillations. Band mapping measured interband transfer to locate valley-gap closures. Opposite-force drift measurements distinguished cancelling from reinforcing Berry curvature contributions; reversing shaking handedness reversed the drift. Interband transfer and harmonic confinement affected this readout, including drift beyond ideal transition lines. This neutral-atom Floquet Hamiltonian realization did not measure electrical sheet conductivity. Magnetic fields were used for preparation and forcing.^[ref]