Spectrum of OpN(cos(p+p₀) + λ·cos(z+ζ₀))

λ (coupling)
p₀ (quasi-momentum)
ζ₀ (Bloch phase)