Generalized Aubry-André formula and continuity of the intersection spectrum of the Almost Mathieu operator

Abstract: We consider the spectrum of the Almost Mathieu operator (AMO) and show that the moments of the restriction of the Lebesgue measure to the intersection spectrum $\text{Leb}|{Σ{α,λ}}$ are polynomials in coupling $λ$ with coefficients that are trigonometric polynomials in frequency $α$. The statement can be considered as a generalization of the Aubry-André formula for the measure of the spectrum of AMO. As a corollary, we obtain that the restriction of the Lebesgue measure to the intersection spectrum that we denote by $μ^{-}_{α, λ}$ depends continuously on the parameters (frequency $α$ and coupling $λ$) in weak-* topology. Moreover, we prove that the dependence is not just continuous but analytic in $λ$ and $C^{\infty}$ in $α$ in a sense that an integral of an analytic test function $\varphi(x)$ with respect to $μ^{-}_{α, λ}$ has the same kind of dependence. In particular, this implies that the Lebesgue measure of the part of the spectrum $Σ_{α,λ}$ that lies between two gaps depends analytically on the coupling constant $λ$ and $C^{\infty}$ on the frequency $α$ in an open domain (away from the critical coupling $λ=1$) where these gaps do not bifurcate.

• The paper introduces a generalized Aubry-André formula, extending spectral measure analysis to all even moments with explicit polynomial representations.
• It demonstrates that the AMO's intersection spectrum measure depends continuously on parameters (α, λ) in the weak-* topology, ensuring analytic regularity except at critical points.
• Using trace formulas and symmetry properties, the study reconciles prior spectral ambiguities and offers tools for analyzing quantum transport and localization.

Generalized Aubry-André Formula and Continuity of the Intersection Spectrum of the Almost Mathieu Operator

Introduction and Historical Context

The Almost Mathieu Operator (AMO), defined on $\ell^2(\mathbb{Z})$ by

$$(H_{\alpha, \lambda, \theta} \psi)_n = \psi_{n+1} + \psi_{n-1} + 2\lambda \cos(2\pi(n\alpha + \theta))\psi_n,$$

is foundational in the modern spectral theory of ergodic Schrödinger operators. Its spectral phenomena, including the Hofstadter butterfly—an intricate fractal structure revealed numerically by Hofstadter—have played a central role in the mathematics and physics of quasicrystals and electron dynamics in magnetic fields. Spectral measures and the nature of the spectrum (absolutely continuous, singular continuous, pure point) for various regimes of coupling $\lambda$ and irrational frequencies $\alpha$ have motivated conjectures with broad implications, such as the Aubry-André conjecture and the Ten Martini Problem.

Figure 1: Hofstadter butterfly, $\lambda=1$, exhibiting the fractal structure of the AMO spectrum as a function of frequency $\alpha$.

Intersection Spectrum and Parameter Continuity

While prior work established the continuity of the union spectrum $$\Sigma_{\alpha,\lambda}^+ = \bigcup_\theta \Sigma_{\alpha,\lambda,\theta}$$ in $\alpha$ and $\lambda$ in the Hausdorff metric, the behavior of the intersection spectrum

$$\Sigma^-_{\alpha, \lambda} = \bigcap_{\theta} \Sigma_{\alpha,\lambda, \theta}$$

was less transparent. This paper shows that continuity for the intersection spectrum should be understood through the family of measures $$(H_{\alpha, \lambda, \theta} \psi)_n = \psi_{n+1} + \psi_{n-1} + 2\lambda \cos(2\pi(n\alpha + \theta))\psi_n,$$0, i.e., Lebesgue measure restricted to $$(H_{\alpha, \lambda, \theta} \psi)_n = \psi_{n+1} + \psi_{n-1} + 2\lambda \cos(2\pi(n\alpha + \theta))\psi_n,$$1. It is proven that $$(H_{\alpha, \lambda, \theta} \psi)_n = \psi_{n+1} + \psi_{n-1} + 2\lambda \cos(2\pi(n\alpha + \theta))\psi_n,$$2 depends continuously on $$(H_{\alpha, \lambda, \theta} \psi)_n = \psi_{n+1} + \psi_{n-1} + 2\lambda \cos(2\pi(n\alpha + \theta))\psi_n,$$3 in the weak-* topology, thus reconciling prior confusion about spectral continuity at rational and irrational frequencies.

Figure 2: Hofstadter butterfly for $$(H_{\alpha, \lambda, \theta} \psi)_n = \psi_{n+1} + \psi_{n-1} + 2\lambda \cos(2\pi(n\alpha + \theta))\psi_n,$$4, magnified near $$(H_{\alpha, \lambda, \theta} \psi)_n = \psi_{n+1} + \psi_{n-1} + 2\lambda \cos(2\pi(n\alpha + \theta))\psi_n,$$5 with the intersection spectrum $$(H_{\alpha, \lambda, \theta} \psi)_n = \psi_{n+1} + \psi_{n-1} + 2\lambda \cos(2\pi(n\alpha + \theta))\psi_n,$$6 indicated. Weak convergence of measures is visible as $$(H_{\alpha, \lambda, \theta} \psi)_n = \psi_{n+1} + \psi_{n-1} + 2\lambda \cos(2\pi(n\alpha + \theta))\psi_n,$$7.

Generalized Aubry-André Formula: Moments as Polynomials

A principal result of the paper is the discovery that moments of the measure $$(H_{\alpha, \lambda, \theta} \psi)_n = \psi_{n+1} + \psi_{n-1} + 2\lambda \cos(2\pi(n\alpha + \theta))\psi_n,$$8 are polynomials in $$(H_{\alpha, \lambda, \theta} \psi)_n = \psi_{n+1} + \psi_{n-1} + 2\lambda \cos(2\pi(n\alpha + \theta))\psi_n,$$9 whose coefficients are trigonometric polynomials in $\lambda$0. Explicitly, for even $\lambda$1,

$\lambda$2

can be written as

$\lambda$3

where $\lambda$4 has degree $\lambda$5 in $\lambda$6 and coefficients in $\lambda$7. Odd moments vanish by spectral symmetry. The leading coefficient and the polynomial structure are rigorously justified via trace formulas linked to the symmetry properties of the AMO at rational $\lambda$8.

This result subsumes the original Aubry-André formula for the Lebesgue measure of the spectrum, extending it to higher moments and thus providing more refined information about the spectral distribution.

Analytic and Smooth Dependence on Parameters

The study further establishes that the measure of any connected component of the spectrum (i.e., between two persistent gaps) is analytic in $\lambda$9 and $\alpha$0 in $\alpha$1, as long as the relevant gaps do not bifurcate within a parameter domain. The analyticity persists except at the non-analytic critical point $\alpha$2, but passage to normalized measures removes even this singularity.

The technical approach relies on a four-trace formula, reducing the analysis of measures and their moments to tractable computations with traces of finite- and infinite-dimensional operators possessing certain reflection symmetries. Exponential decay (Combes–Thomas-type) estimates for matrix elements ensure convergence and regularity.

Figure 3: Hofstadter butterfly for fixed $\alpha$3 near $\alpha$4, showing the role of spectral symmetry in the trace computations used for moments.

Figure 4: Magnified Hofstadter butterfly for $\alpha$5, highlighting a region between non-bifurcating gaps where analyticity in $\alpha$6 is established.

Explicit Computation of Moments and Spectral Identities

The paper provides explicit polynomials for low-order moments, reinforcing the generalized Aubry-André structure: $\alpha$7 with similar expressions for higher moments. Furthermore, duality relations between moments at $\alpha$8 and $\alpha$9 are established, illustrating a robust and deep algebraic structure in the AMO’s spectral measures.

Figure 5: Traces of transfer matrices $\lambda=1$0 at rational frequencies, underpinning the construction of moment polynomials and the four-trace formula.

Singularity and Non-analyticity in Frequency

While analyticity in $\lambda=1$1 is robust (and can be "restored" at $\lambda=1$2 by normalization), dependence on frequency $\lambda=1$3 is $\lambda=1$4 but not analytic for general test functions. The argument rests on explicit counterexamples motivated by invariant graph constructions over irrational rotations, specifically demonstrating non-analytic behavior in constructed quasi-periodic cocycles even when all formal derivatives exist.

Implications and Prospects for Future Research

The results provide both conceptual and computational advances:

• Parameter regularity: The established continuity and analyticity of spectral measures constitute a definitive answer to prior ambiguities regarding the limiting behavior of spectral quantities at rational/irrational $\lambda=1$5 and away from the critical coupling.
• Fine spectral invariants: By encoding higher moments of the intersection spectral measure as explicit polynomials, new invariants for spectral analysis and their duality relations become available. These may inform the study of quantum transport, localization, and the multifractal features of the spectrum.
• Extension to broader classes: The methodological framework—combining trace identities, symmetry, and cocycle dynamics—invites extension to more general quasi-periodic Schrödinger potentials and more complex spectral phenomena, as well as potential applications in mathematical physics beyond one-dimensional models.

Figure 6: Regions in $\lambda=1$6 parameter space for fixed rational approximation, clarifying the loci where spectral regularity results uniformly apply.

Conclusion

This work provides a rigorous and explicit framework for understanding the parameter dependence of spectral measures for the Almost Mathieu operator, generalizing the Aubry-André formula to all moments and proving continuity and analyticity properties previously unclear for intersection spectra. The approach—blending explicit combinatorial operator computations, symmetry exploitation, and dynamical systems reasoning—sets a new standard for achieving explicit spectral information and parameter regularity in quasi-periodic operator theory. The polynomial invariants and trace formulae provided offer powerful tools for further exploration of the AMO and related operators.