Band Structure
The band structure for the Hofstadter model is computed by constructing a set of simultaneous timeindependent Schrödinger equations for the basis sites, and then using the plane wave ansatz to formulate each equation as a difference relation. For simplicity, we start by discussing the band structure derivation for the conventional Hofstadter model, and then proceed to outline the general case.
To recap from the previous section, for Landau gauge in the ydirection \(\mathbf{A}=By \hat{\mathbf{e}}_x\), the general formula for the Peirels phase acquired by an electron hopping from site \(i=(X_i,Y_i)\) to site \(j=(X_j,Y_j)\) is given as
where \(\Delta X = X_j  X_i\), \(\Delta Y = Y_j  Y_i\), and we have defined the flux density with respect to the UC area, such that \(n_\phi=BA_\mathrm{UC}/\phi_0 = Ba^2/\phi_0\). In what follows, let us index the UC coordinates using \((m,n)\). The timeindependent Schrödinger equation for the conventional nearestneighbor (NN) squarelattice Hofstadter model is then given as
where \(\Psi_{m,n}\) is the wavefunction at unit cell \((m,n)\) and \(t\) is the hopping amplitude. Since we have taken Landau gauge in the ydirection, it is reasonable to assume plane wave behavior in the xdirection. Hence we can simplify the equation further by invoking the plane wave ansatz \(\Psi_{m,n} = e^{\mathrm{i}\mathbf{r}\cdot\mathbf{k}}\psi_{n}\), such that
which simplifies to
where
All distances are measured in units of the lattice constant \(a=1\). Given that there are \(q\) sites in our MUC, this is a \(q\)dimensional difference equation, known as the Harper equation, which is a special case of an almost Mathieu operator [Harper55]. This equation may be readily written in matrix form, such that \(\mathbf{H} \boldsymbol{\psi} = E \boldsymbol{\psi}\), where \(\mathbf{H}\) is a \(q\times q\) Hamiltonian matrix and \(\boldsymbol{\psi}=(\psi_0, \psi_1, \dots, \psi_{q1})^\intercal\) is a vector of length \(q\). For the case of \(n_\phi=1/4\), the Hamiltonian matrix may be written explicitly as
and hence there are 4 energy bands in the spectrum, as shown in the figures below.
Note
It is not possible to define a highsymmetry path for the generalized Hofstadter model as a continuous function of lattice anisotropy and obliqueness. Instead, we use reference paths to plot the 2D band structures.
Lattice 
UC Lattice Vectors 
UC Basis Vectors 
Reference Path 


\(\mathbf{a}_1 = a (1, 0)\), \(\mathbf{a}_2=a(0, 1)\) 
\(\mathbf{a}_{\text{b},1}=(0, 0)\) 
\(\Gamma\to X \to M \to Y \to \Gamma\) 

\(\mathbf{a}_1 = a (1, 0)\), \(\mathbf{a}_2=a (1/2, \sqrt{3}/2)\) 
\(\mathbf{a}_{\text{b},1}=(0, 0)\) 
\(\Gamma \to K \to M \to K' \to \Gamma\) 

\(\mathbf{a}_1 = a (1, 0)\), \(\mathbf{a}_2=\alpha a (\cos \theta, \sin \theta)\) 
\(\mathbf{a}_{\text{b},1}=(0, 0)\) 
\(\begin{cases}\Gamma \to K \to M \to K' \to \Gamma, \;\;\; &\text{when }\theta\text{ is a multiple of }\pi/3\text{ but not }\pi/2 \\ \Gamma \to X \to M \to Y \to \Gamma, \;\;\; &\text{otherwise}\end{cases}\) 

\(\mathbf{a}_1 = a (1, 0)\), \(\mathbf{a}_2=\alpha a (\cos \theta, \sin \theta)\) 
\(\mathbf{a}_{\text{b},1}=(0, 0)\), \(\mathbf{a}_{\text{b},2}=(a_{1,x}/2, a_{2,y}/3)\) 
\(\Gamma \to K \to M \to K' \to \Gamma\) 

\(\mathbf{a}_1 = a (1, 0)\), \(\mathbf{a}_2=\alpha a (\cos \theta, \sin \theta)\) 
\(\mathbf{a}_{\text{b},1}=(0, 0)\), \(\mathbf{a}_{\text{b},2}=(a_{1,x}/2, 0)\), \(\mathbf{a}_{\text{b},3}=(a_{2,x}/2, a_{2,y}/2)\) 
\(\Gamma \to K \to M \to K' \to M' \to \Gamma\) 
The reference points are defined as: \(\Gamma = (0, 0)\), \(M = (1/2, 1/2)\), \(X = (1/2, 0)\), \(Y = (0, 1/2)\), \(K = (2/3, 1/3)\), \(K' = (1/3, 2/3)\), \(M' = (0, 1/2)\), in fractional units of the MUC reciprocal lattice vectors \((\mathbf{b}_1, \mathbf{b}_2)\). The difference between \(Y\) and \(M'\) is simply convention, depending on whether we are referring to rectangular or hexagonal Brillouin zones.
We emphasize that the paths defined above are only highsymmetry paths in special cases, where the corresponding symmetries are present, such as at zero magnetic field. In all other cases, these are simply reference paths through \(k\)space to facilitate easy comparison, and the full 3D band structure should also be examined.
In the general case, the procedure follows in a similar way. We start by writing down the timeindependent Schrödinger equation for each site in the basis, which we index using lowercase Greek letters \(\alpha,\beta\). This yields a set of \(N_\mathrm{b}\) simultaneous equations, which we can write in matrix form, such that \(\mathbf{H}\boldsymbol{\psi}=E\boldsymbol{\psi}\), where \(\boldsymbol{\psi}=(\psi^0, \psi^1, \dots, \psi^{N_\mathrm{b}1})^\intercal\) is a vector of length \(N_\mathrm{b}\), and the \(N_\mathrm{b}\times N_\mathrm{b}\) Hamiltonian matrix is given as
where \(H^{\alpha\beta}\) is the Hamiltonian for hoppings from sublattice \(\alpha\) to sublattice \(\beta\). Then, for each \(H^{\alpha\beta}\) we can write down a \(q\times q\) Harper matrix, as before. For Landau gauge in the ydirection \(\mathbf{A}=By \hat{\mathbf{e}}_x\), the general form of the Harper matrix is given as
where \(\Lambda_{l, n}\) is the diagonal function, and we have dropped the \(\alpha\beta\) superscripts for readability. The second matrix simply accounts for rolled over boundary terms. Since we are working in Landau gauge in the ydirection, the diagonal function may be written as
where \(\langle \dots \rangle^l_\kappa\) denotes the subset of \(\kappa\)th nearest neighbors with a net \(y\) unit cell displacement of \(l\), \(\theta_{ij}\) is the Peierls phase, \(\\\mathbf{k}\) is the momentum vector, and \(\\\mathbf{r}\) is the displacement vector. We emphasize that hoppings that are related by Hermitian conjugation, which are outside the scope of the unit cell, are not included in the diagonal function matrix. For example, for NN hopping on the triangular lattice, we include 3 of the 6 nearest neighbors in the diagonal function matrix, and the rest are captured when we add on the Hermitian conjugate to the Hamiltonian. However, for NN hopping on the kagome lattice, we include all 4 of 4 nearest neighbors for each basis site, since they are all within the same unit cell. Overall, we are left we left with a \(N_\mathrm{b}q\times N_\mathrm{b}q\) block Hamiltonian matrix, which yields \(N_\mathrm{b}q\) bands.
In HofstadterTools, we can analyze the resulting band structure by computing its key properties, which are listed in the tables below. These band properties may be selected by passing flags to the band_structure
program, and are grouped by computational expense. By default, HofstadterTools prints properties in the basic group only (for speed reasons). When computing topology and quantum geometry properties of bands, it is important to use manifestly gauge invariant expressions, so that we omit spurious Bloch phase factors and can compute the quantities quickly and accurately. To this end, we use the Fukui formula to compute the (first) Chern number [Fukui05] and the projector formalism to compute the quantum geometric tensor [Mera22].
Note
The Chern numbers of the bands may also be inferred by plotting the Wilson loops, which are the products of Berry phases around cycles of the Brillouin zone, as shown in the middle figure above. The magnitude of the Chern number corresponds to the number of windings of the Wilson loop and the sign of the Chern number corresponds to its direction.
In the projector formalism, the quantum geometric tensor is defined as
where \(\mathcal{P}_\mathbf{k} = \sum_n^{N_\mathrm{g}} \ket{u_n(\mathbf{k})} \bra{u_n(\mathbf{k})}\) is the band projector, \(\ket{u_n(\mathbf{k})}\) is the eigenvector of band \(n\) at momentum \(\mathbf{k}\), and \(N_\mathrm{g}\) is the number of touching bands in a band group. The real part of the quantum geometric tensor is given by the FubiniStudy metric \(g_{\mu\nu}(\mathbf{k})=\Re[\mathcal{R}_{\mu\nu}]\), which corresponds to the distance between eigenstates on the Bloch sphere, whereas the imaginary part of the quantum geometric tensor is given by the Berry curvature \(\mathcal{B}(\mathbf{k})=2 \Im [\mathcal{R}_{01}(\mathbf{k})]\). Crucially, since band geometry and topology are components of the same tensor, we can derive relations between them, namely
where we define \(\mathcal{D}\) as the determinant inequality saturation measure (DISM) and \(\mathcal{T}\) as the trace inequality saturation measure (TISM). It has been shown analytically that when the trace(determinant) inequality is saturated for a Chern band, the algebra of projected density operators is identical(isomorphic) to that in Landau levels [Roy14].
Using these band properties, we can perform several sanity checks on our computed band structures. In terms of band topology, we know that all of the Chern numbers in a Hofstadter spectrum must sum to zero. In terms of band geometry, we know that as we take the Landau level limit \(q\to\infty\), the TISM and DISM must monotonically approach zero from above.
Note
The band structures can also be checked by comparing against results in the literature. For example, the Chern numbers can be benchmarked against Fig.2.6 of [AidelsburgerPhD] and the values of the TISM can be benchmarked against Fig.3 of [Bauer22].
Alias 
Name 
Symbol/Formula 
Description 


Band Number 
\(n\) 
Bands are numbered in ascending order, with respect to energy, starting from zero. 

Band Group 
N/A 
Bands are considered grouped when they are touching within the band gap threshold. 

Isolated Band Flag 
N/A 
A band is isolated when it is not touching any other bands, i.e. it is a member of a band group of size one. 

Band Width 
\(W = \mathrm{max}(E_n)  \mathrm{min}(E_n)\) 
The band width for an isolated band (group) is the difference between the largest and smallest energies in that band (group). 

Band Gap 
\(\Delta = \mathrm{min}(E_{n+1})  \mathrm{max}(E_n)\) 
The band gap for an isolated band (group) is the difference between the smallest and largest energies of the subsequent and current bands. 

Gaptowidth Ratio 
\(\Delta/W\) 
The gaptowidth ratio for an isolated band (group) is defined as the ratio between the band gap and width. 
Alias 
Name 
Symbol/Formula 
Description 


Berry Curvature Fluctuations 
\(\hat{\sigma}_{\mathcal{B}} = \sigma_\mathcal{B} / \mu_\mathcal{B}\) 
For \(C=1\) bands, this definition is equivalent to Eq.(6) from [Jackson15]. 

(First) Chern Number 
\(C=\frac{1}{2\pi} \iint_\mathrm{BZ} \mathcal{B}(\mathbf{k}) \mathrm{d}^2 k\) 
The Chern number is an integer and is computed using the Fukui formula [Fukui05]. 
Alias 
Name 
Symbol/Formula 
Description 


FubiniStudy Metric Fluctuations 
\(\sigma_{g} = \sqrt{\frac{1}{2}\sum_{ij} \sigma^2_{g_{ij}}}\) 
This definition is equivalent to Eq.(8) from [Jackson15]. 

Mean of the Diagonal FubiniStudy Metric 
\(\mu_{g_{xx}}\) 
This quantity is equal to \(\mu_{g_{yy}}\) and is studied in [Bauer16]. 

Mean of the Offdiagonal FubiniStudy Metric 
\(\mu_{g_{xy}}\) 
This quantity is equal to \(\mu_{g_{yx}}\) and is studied in [Bauer16]. 

Standard Deviation of the Diagonal FubiniStudy Metric 
\(\sigma_{g_{xx}}\) 
This quantity is equal to \(\sigma_{g_{yy}}\) and is studied in [Bauer16]. 

Standard Deviation of the Offdiagonal FubiniStudy Metric 
\(\sigma_{g_{xy}}\) 
This quantity is equal to \(\sigma_{g_{yx}}\) and is studied in [Bauer16]. 

Brillouinzoneaveraged Trace Inequality Saturation Measure (TISM) 
\(\langle\mathcal{T}\rangle = \langle \mathrm{tr}\,g(\mathbf{k})  \mathcal{B}(\mathbf{k}) \rangle\) 
The TISM is nonnegative and defined in Eq.(10) of [Jackson15]. 

Brillouinzoneaveraged Determinant Inequality Saturation Measure (DISM) 
\(\langle\mathcal{D}\rangle = \langle \mathrm{det}\,g(\mathbf{k})  \mathcal{B}^2(\mathbf{k})/4 \rangle\) 
The DISM is nonnegative and defined in Eq.(9) of [Jackson15]. 