The Colored Hofstadter Butterfly for the Honeycomb Lattice
  We rely on a recent method for determining edge spectra and we use it to compute the Chern numbers for Hofstadter models on the honeycomb lattice having rational magnetic flux per unit cell. Based on the bulk-edge correspondence, the Chern number $\sigma_H$ is given as the winding number of an eigenvector of a $2 \times 2$ transfer matrix, as a function of the quasi-momentum $k \in (0,2 \pi)$. This method is computationally efficient (of order $O(n^4)$ in the resolution of the desired image). It also shows that for the honeycomb lattice the solution for $\sigma_H $ for flux $p/q$ in the $r$-th gap conforms with the Diophantine equation $r=\sigma_H\cdot p+ s\cdot q$, which determines $\sigma_H \mod q$. A window such as $\sigma_H \in(-q/2,q/2)$, or possibly shifted, provides a natural further condition for $\sigma_H$, which however turns out not to be met. Based on extensive numerical calculations, we conjecture that the solution conforms with the relaxed condition $\sigma_H\in(-q,q)$. 
Andrea Agazzi, Jean-Pierre Eckmann and Gian Michele Graf
math-ph cond-mat.mes-hall math.MP math.SP