“…where 𝑁 𝑘 , 𝑁 𝐸 , 𝐸 lim are specified for each simulation. In all figures, the dark areas correspond to the essential spectrum; see Figures 1,2,8,9,12,15,16,17. Outside the dark areas, in particular for (𝑘, 𝐸) ∈ 𝒦 (see (7.1)), we seek edge state curves by studying where (𝑘, 𝐸) ↦ log |Δ(𝑘, 𝐸)| takes on very large negative values. In Figure 12, we consider 4 different edges, three of zigzag-type and one of armchair-type, with numerical parameters and 𝑁 𝑘 = 𝑁 𝐸 = 1000 and specified values of 𝐸 lim .…”