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 _mathrm{H}) is given as the winding number of an eigenvector of a (2 times 2) transfer matrix, as a function of the quasi-momentum (kin (0,2pi )) . This method is computationally efficient (of order (mathcal {O}(n^4)) in the resolution of the desired image). It also shows that for the honeycomb lattice the solution for (sigma _mathrm{H}) for flux (p/q) in the (r) -th gap conforms with the Diophantine equation (r=sigma _mathrm{H}cdot p+ scdot q) , which determines (sigma _mathrm{H}mod q) . A window such as (sigma _mathrm{H}in (-q/2,q/2)) , or possibly shifted, provides a natural further condition for (sigma _mathrm{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 _mathrm{H}in (-q,q)) .