On gaps in the spectra of quasiperiodic Schrödinger operators with discontinuous monotone potentials

Abstract.

We show that, for one-dimensional discrete Schrödinger operators, stability of Anderson localization under a class of rank one perturbations implies absence of intervals in spectra. The argument is based on well-known result of Gordon and Simon, combined with a way to consider perturbations whose ranges are not necessarily cyclic. The main application of the results is showing that a class of quasiperiodic operators with sawtooth-like potentials, for which such a version of stable localization is known, has Cantor spectra. We also obtain several results on gap filling under rank one perturbations for some general (not necessarily monotone) classes of quasiperiodic operators with discontinuous potentials.

1. Introduction and statements of the results

1.1. Motivation and some of the main results

Let $f\colon\mathbb{R}\to[-\infty,+\infty)$ be a $1$-periodic function, and $H(x)$ be a family of discrete quasiperiodic Schrödinger operators on $\ell^{2}(\mathbb{Z})$:

(1.1) $$(H(x)\psi)(n)=\psi(n+1)+\psi(n-1)+f(x+n\alpha)\psi(n),\quad x\in[0,1).$$

Let $\gamma>0$. We will say that $f$ is $\gamma$-monotone if

(1.2) $$f(y)-f(x)\geq\gamma(y-x),\quad 0\leq x<y<1.$$

Two typical examples of $\gamma$-monotone potentials are $f(x)=\{x\}=x-\lfloor x\rfloor$ (the sawtooth potential) and $f(x)=\lambda\cot(\pi x)$, $\lambda>0$ (the Maryland model). We will say that $f$ is sawtooth-type if $f$ is $\gamma$-monotone and $f|_{[0,1)}$ extends to a (bounded) continuous map from $[0,1]$ to $\mathbb{R}$. Similarly, $f$ is Maryland-type if it is $\gamma$-monotone and $f|_{[0,1)}$ extends to a homeomorphism between $[0,1]$ and $[-\infty,+\infty]$.

In a series of works [13, 16, 15, 14], it is shown that a large class of one-dimensional quasiperiodic operators with $\gamma$-monotone potentials satisfies Anderson localization. As an example, suppose that $\alpha$ is Diophantine:

$$\operatorname{{dist}}(n\alpha,\mathbb{Z})\geq C|n|^{-\tau},\quad\forall n\in\mathbb{Z}\setminus\{0\}.$$

Then the results of [13, 16] show that Anderson localization holds for sawtooth-type and Maryland-type $f$ for almost all $x\in[0,1)$ (see Proposition 2.6 and Remark 2.7 below for more complete statements).

Denote by $\Sigma_{f}$ the almost sure spectrum of the operator (1.1). Shortly after the results of [13] were announced, a natural question was raised: suppose that $f$ is sawtooth-type. Can $\Sigma_{f}$ contain an interval, or should it be a Cantor set? In the referee report on [13], it was suggested that, for sufficiently large $\gamma$ and perhaps under some smoothness assumptions on $f$, one should be able to show that $\Sigma_{f}$ is an interval.

One of our main results is Theorem 2.9, which provides some answers to this question in the regime of localization for $\gamma$-monotone potentials. It is convenient to postpone the complete formulation of this result until Section 2. However, we state the following as a separate theorem (in fact, contained in Corollary 2.2 of the main Theorem 2.9).

Theorem 1.1.

Let $f$ be sawtooth-type function and suppose that $\alpha$ is Diophantine. Then $\sigma(H(x))$ is a Cantor set $($in particular, it contains no intervals$)$.

We note that the question of Cantor spectrum for quasiperiodic operators has been an active research area. Some recent results in this direction include [1, 2, 4, 8, 10]. See also the review [7] and references therein.

It is interesting to compare these results with those for Maryland-type potentials, which sometimes do not rely on localization and/or $\gamma$-monotonicity. We will call $f$ weak Maryland-type if $f|_{[0,1)}$ extends to a continuous map from $[0,1]$ to $[-\infty,+\infty]$ and $f(0)=-\infty$, $f(1-0)=+\infty$. We will also call $f$ weak sawtooth-type if $f|_{[0,1)}$ extends to a (bounded) continuous map from $[0,1]$ to $\mathbb{R}$, with

$$-\infty<f(0)<f(1-0)<+\infty.$$

It was shown in [16] that, for weak Maryland-type $f$ and $\alpha\in\mathbb{R}\setminus\mathbb{Q}$, one has $\sigma(H(x))=\mathbb{R}$. One can argue by identifying $\mathbb{R}\cup\{\infty\}=\overline{\mathbb{R}}$ with a circle, somewhat informally, that the range of $f$, in the case of weak sawtooth-type potentials, has a “gap” $\mathbb{R}\setminus f([0,1))$, which is not present in the Maryland-type case. In other words, one “gap” in the values of $f$ produces multiple gaps in the spectra. In the second part of our paper, we partially confirm this intuition by demonstrating the following gap filling phenomenon. Let $f$ be a weak sawtooth-type function. For $t\in\mathbb{R}\setminus[f(0),f(1-0)]$, consider the following family:

(1.3) $$f_{t}(x):=\begin{cases}f(x),&x\in(0,1);\\ t,&x=0.\end{cases}$$

Let $H_{t}(x)$ be $H(x)$ with $f$ replaced by $f_{t}$:

(1.4) $$(H_{t}(x)\psi)(n)=\psi(n+1)+\psi(n-1)+f_{t}(x+n\alpha)\psi(n),\quad x\in[0,1),\quad t\in\mathbb{R}\setminus f((0,1)).$$

Theorem 1.2.

Let $f$ be a weak sawtooth-type function and $\alpha\in\mathbb{R}\setminus\mathbb{Q}$. Then

(1.5) $$\overline{\bigcup\limits_{\mathbb{R}\setminus[f(0),f(1-0)]}\sigma(H_{t}(0))}=\mathbb{R}.$$

The closure in the left hand side of (1.5) would not be necessary if one included $t=\infty$ into consideration by allowing infinite potentials. Again, we postpone the more general setting until Section 4 (see, in particular, Theorems 4.2 and 4.4), where this inclusion is performed systematically. In some situations, one can even prove that every point of $\mathbb{R}$ is an actual eigenvalue of some operator from our family, see Corollary 4.3. Note that, if one considers the union over all $t\in\mathbb{R}$ in (1.5), the result immediately follows from spectral averaging (for example, [17, Theorem 11.8]), which is why we assume $f(0)<f(1-0)$.

The general version of Theorem 1.2 also holds in higher dimensions for operators on $\mathbb{Z}^{d}$. A similar question about Theorem 1.1 appears to be more difficult. While the general argument behind Theorem 1.1 is abstract, it relies on two properties of the operator: Anderson localization that is stable under a family of rank one perturbations (that is, replacing $f$ by $f_{t}$) and the fact that a standard basis vector is cyclic on a sufficiently large subspace. In Theorem 1.1, both of these properties rely on the one-dimensional structure of the operators, the first one through [14], and the second one through Lemma 2.1 below. In an upcoming work, we intend to address the first question in a perturbative setting on $\mathbb{Z}^{d}$.

1.2. Stability of localization and absence of intervals

We will now state the main abstract result that will be used to prove Theorem 1.1. It relies on the following classical result that, as far as the authors are aware, has been discovered independently in [11, 12, 9], see also [5].

Proposition 1.3.

Let $A$ be a self-adjoint operator on a Hilbert space with simple spectrum and a cyclic vector $\varphi$. Let $P$ be the orthogonal projection onto the span of $\{\varphi\}$. Then, for every $t$ in some dense $G_{\delta}$-subset of $\mathbb{R}$, all eigenvalues of the operator $A+tP$ are isolated.

For Anderson-type Hamiltonians whose point spectra are unions of intervals, this result implies that Anderson localization is unstable under rank one perturbations: indeed, for the set of values of $t$ under consideration, $H$ will not have any eigenvalues on any such interval, and therefore the spectral measure must contain a continuous component.

Our goal is a somewhat unusual application of this result in a reverse direction: if an operator has a “very stable” kind of localization, then its spectrum cannot contain an interval. More precisely, let $H$ be a discrete Schrödinger operator on $\ell^{2}(\mathbb{Z})$:

(1.6) $$H=\Delta+V,\quad(H\psi)(n)=\psi(n+1)+\psi(n-1)+V(n)\psi(n)$$

with real-valued potential $V$. Let also $P$ be a projection operator on the standard basis vector $e_{0}$: that is,

$$(P\psi)(n)=\psi(0)e_{0}.$$

The following is our main abstract result.

Theorem 1.4.

Define $H$ and $P$ as above. Let $J\subset\mathbb{R}$ be an open interval such that, for some $t_{0}>0$, the spectrum of $H+tP$ is purely point in $J$ for all $t\in[0,t_{0})$. Then $\sigma(H)\cap J$ does not contain an interval.

Remark 1.5.

Since $P$ is a rank one operator, the conclusion of the theorem also holds for $H+tP$ for all $t$.

Clearly, the result immediately follows from Proposition 1.3 if we assume that $e_{0}$ is a cyclic vector for $H$. The essential part of the proof is making sure that this assumption is satisfied on a sufficiently large subspace for one-dimensional Schrödinger operators with point spectra (we also do not a priori assume that the spectrum of $H$ is simple, although the spectrum on $J$ will be such since it is purely point).

1.3. Structure of the paper

In Section 2, we prove Theorem 1.4. Afterwards, we introduce the complete setting for quasiperiodic operators with $\gamma$-monotone potentials and prove Theorem 2.9, which implies Theorem 1.1 from the Introduction. In Section 3, we discuss Schrödinger operators with infinite coupling, which are used in the proof of Theorem 1.2. In Section 4, we state the main results on gap filling (Theorems 4.2 and 4.4), which imply Theorem 1.2.

The proofs of the gap filling results are based on rational approximations, for which it can be established by a topological argument (Proposition 4.14). The main difficulty is continuity, since the usual arguments such as [3] rely on compactness. For the operator family (1.6), two obstructions to compactness are related to the presence of discontinuities and the fact that the potential may be unbounded. In Subsections 4.3 and 4.4 we resolve both of these issues by constructing an extended operator family, using the results of Section 3 to deal with unbounded potentials and auxiliary Cantor-type set in the parameter space in order to deal with discontinuities.

1.4. Acknowledgements

I. K. was supported by the NSF grants DMS–1846114, DMS–2052519, and the 2022 Sloan Research Fellowship. L. P. was supported by the EPSRC grant EP/V051636/1. R. S. was supported by the NSF grant DMS–2306327.

The authors would also like to thank S. Jitomirskaya for valuable discussions.

2. Absence of intervals in spectra

We will start from proving Theorem 1.4. Afterwards, we will establish its consequences for Schödinger operators with $\gamma$-monotone potentials, the most notable being Theorem 2.9 and Corollary 2.10.

2.1. Proof of Theorem 1.4

In the following lemma, $\mathbf{1}_{J}(H)$ denotes the spectral projector of a self-adjoint operator $H$ associated to an interval $J\subset\mathbb{R}$, and $e_{0}$ is the standard basis vector of $\ell^{2}(\mathbb{Z})$ associated to the origin. Suppose that $\varphi$ is an eigenfunction of $H$ with $\varphi(0)=0$. Then, the functions

$$\varphi_{\pm}:=\varphi|_{\ell^{2}(\mathbb{Z}_{\pm})}$$

are eigenfunctions of the respective half-line restrictions of $H$. Since an eigenfunction of $H$ cannot vanish on a half-line, one has

(2.1) $$\max\{\|\varphi_{+}\|,\|\varphi_{-}\|\}\leq\sqrt{1-\delta^{2}}\|\varphi\|,\quad\delta=\delta(\varphi)=\frac{1}{\|\varphi\|}\min\{\|\varphi_{+}\|,\|\varphi_{-}\|\}>0.$$

In our notation, $\|\cdot\|$ will always denote the $\ell^{2}$-norm.

Lemma 2.1.

Let $H$ be a discrete one-dimensional Schrödinger operator (1.6). Suppose that $J\subset\mathbb{R}$ is an open interval such that $\sigma(H)\cap J$ is pure point and $\mathbf{1}_{J}(H)e_{0}=0$ $($in other words, every eigenfunction with eigenvalue in $J$ vanishes at the origin$)$. Let $\varphi$, $\psi$ be two distinct normalized eigenfunctions with eigenvalues in $J$:

$$H\varphi=\lambda\varphi,\,\,H\psi=\mu\psi,\,\,\lambda,\mu,\in J,\,\,\lambda\neq\mu,\,\,\|\varphi\|=\|\psi\|=1,\,\,\varphi(0)=\psi(0)=0.$$

Then

(2.2) $$m:=\min\left\{\max\{\|\varphi_{-}\|,\|\psi_{+}\|\},\max\{\|\varphi_{+}\|,\|\psi_{-}\|\}\right\}\leq\max\{\|\varphi_{+}\|,\|\varphi_{-}\|\}\leq\sqrt{1-\delta^{2}},$$

where $\delta=\delta(\varphi)>0$ does not depend on $\mu$. Moreover,

(2.3) $$\operatorname{{dist}}(\lambda,\partial J)\leq\frac{|\lambda-\mu|}{\sqrt{1-m^{2}}}.$$

Proof.

Note that the inequalities in (2.2) follow from (2.1). Without loss of generality, assume that $m=\max\{\|\varphi_{-}\|,\|\psi_{+}\|\}$. Let

$$f:=a\varphi_{-}+b\psi_{+},$$

where $a,b$ are chosen in such a way that $(\Delta f)(0)=0$. Note that such choice is always possible with $a,b\neq 0$. Though not important, one can also choose $a,b\in\mathbb{R}$, since all eigenfunctions can be chosen to be real-valued. Denote

$$Q:=\mathbf{1}_{J}(H).$$

Suppose, $\theta\in\mathrm{Ran}(Q)$ is another eigenfunction of $H$ (with an eigenvalue in $J$). Then all six functions $\varphi_{-},\varphi_{+},\psi_{-},\psi_{+},\theta_{-},\theta_{+}$ are mutually orthogonal, either due to disjoint supports, or because of being distinct eigenfunctions of the half-line operator. We have

$$Qf=a\langle\varphi_{-},\varphi\rangle\varphi+b\langle\psi_{+},\psi\rangle\psi=a\|\varphi_{-}\|^{2}\varphi+b\|\psi_{+}\|^{2}\psi,$$

$$\frac{\|Qf\|^{2}}{\|f\|^{2}}=\frac{a^{2}\|\varphi_{-}\|^{4}+b^{2}\|\psi_{+}\|^{4}}{a^{2}\|\varphi_{-}\|^{2}+b^{2}\|\psi_{+}\|^{2}}\leq m^{2}<1.$$

As a consequence,

$$\|(1-Q)f\|^{2}\geq(1-m^{2})\|f\|^{2}.$$

Due to the choice of $f$, we also have

$$(H-\lambda)f=b(\mu-\lambda)\psi_{+},\quad\text{therefore}\quad\|(H-\lambda)f\|^{2}\leq(\mu-\lambda)^{2}\|f\|^{2}.$$

On the other hand, from the spectral theorem one always has

$$\|(H-\lambda)f\|^{2}\geq\operatorname{{dist}}(\lambda,\partial J)^{2}\|(1-Q)f\|^{2}.$$

Combining the above, we finally obtain

$$\operatorname{{dist}}(\lambda,\partial J)^{2}(1-m^{2})\leq(\mu-\lambda)^{2},$$

which completes the proof. ∎

Corollary 2.2.

Under the assumptions of Lemma 2.1, an eigenvalue of $H$ contained in $J$ cannot be a limit point of other eigenvalues of $H$. As a consequence, $\sigma(H)\cap J$ cannot contain an interval.

Proof.

Let $\lambda\in J$ be an eigenvalue of $H$ and $\mu_{j}\to\lambda$. In view of (2.2), one has (2.3) with $m$ independent of $\mu_{j}$, which leads to a contradiction. ∎

Proof of Theorem 1.4. Let $\mathcal{C}\subset\ell^{2}(\mathbb{Z})$ be the cyclic subspace for $H$ corresponding to the vector $e_{0}$. Recall that in the case of unbounded operators it is defined as

$$\mathcal{C}=\overline{\mathrm{span}\{(H-z)^{-1}e_{0}\colon z\in\mathbb{C}\setminus\mathbb{R}\}}.$$

Both $\mathcal{C}$ and $\mathcal{C}^{\perp}$ are invariant subspaces for $(H-z)^{-1}$ for all $z\in\mathbb{C}\setminus\mathbb{R}$ and for all spectral projections of $H$. In particular, the decomposition $\ell^{2}(\mathbb{Z})=\mathcal{C}\oplus\mathcal{C}^{\perp}$ defines a decomposition of $H$ into a direct sum of two unbounded self-adjoint operators $H_{c}\oplus H_{0}$, with $\sigma(H)=\sigma(H_{c})\cup\sigma(H_{0})$. Note that this can be applied for any $t$, and both $\mathcal{C}$ and $H_{0}$ do not depend on $t$:

$$H+tP=(H_{c}+tP)\oplus H_{0}.$$

As a consequence, the spectra of both operators in $J$ are purely point for $t\in[0,\varepsilon_{0})$. From Proposition 1.3, it follows that $\sigma(H_{c})$ cannot contain an interval. Let $I$ be a connected component of $J\setminus\sigma(H_{c})$. Then, all eigenvectors of $H$ with eigenvalues in $I$ must come from eigenvectors of $H_{0}$, which are orthogonal to $e_{0}$ and therefore vanish at the origin. From Corollary 2.2, $\sigma(H)$ cannot contain an interval in $I$, which completes the proof.  ∎

While not required for our main results, we would like to note that the assumption of $\sigma(H)\cap J$ being pure point is, in fact, not necessary in Lemma 2.1. An appropriate modification is described in the following lemma.

Lemma 2.3.

Let $H$ be a discrete one-dimensional Schrödinger operator (1.6), $\lambda<\lambda_{1}<\lambda_{2}$. Let $\varphi$ be an eigenfunction with eigenvalue $\lambda$:

$$H\varphi=\lambda\varphi,\quad\|\varphi\|=1,\quad\varphi(0)=0.$$

Let $Q:=\mathbf{1}_{[\lambda_{1},\lambda_{2}]}(H)$ be the spectral projection of $H$ associated to the interval $[\lambda_{1},\lambda_{2}]$. Suppose that $Qe_{0}=0$, and that

$$\theta\in\mathrm{Ran}\,Q,\quad\|\theta\|=1.$$

Then $\langle\theta_{+},\varphi_{+}\rangle=0$.

Proof.

First, we will establish the following claim: under the assumptions of the lemma, there exists $\theta^{\prime}\in\mathrm{Ran}\,Q$, $\|\theta^{\prime}\|=1$, satisfying

(2.4) $$|\langle\varphi_{+},\theta_{+}\rangle|\leq\frac{\lambda_{2}-\lambda_{1}}{\lambda_{2}-\lambda}|\langle\varphi_{+},\theta^{\prime}_{+}\rangle|.$$

It is easy to see that the conclusion of the lemma follows from iterating the claim. In order to establish the claim, note that

$$\lambda\langle{\varphi_{+},\theta_{+}}\rangle=\langle H\varphi_{+},\theta_{+}\rangle=\langle\varphi_{+},H\theta_{+}\rangle=\langle\varphi_{+},(H\theta)_{+}\rangle.$$

On the other hand,

$$(H\theta)_{+}=(\lambda_{2}\theta)_{+}+\tilde{\theta}_{+},\quad\text{where}\quad\tilde{\theta}=Q\tilde{\theta},\quad\|\tilde{\theta}\|\leq(\lambda_{2}-\lambda_{1})\|\theta\|=(\lambda_{2}-\lambda_{1}).$$

Let $\theta^{\prime}:=\frac{\tilde{\theta}_{+}}{\|\tilde{\theta}_{+}\|}$. Then

$$|(\lambda-\lambda_{2})\langle\varphi_{+},\theta_{+}\rangle|=|\langle\varphi_{+},\tilde{\theta}_{+}\rangle|\leq(\lambda_{2}-\lambda_{1})|\langle\varphi_{+},\theta^{\prime}_{+}\rangle|,$$

which implies (2.4) and completes the proof.  ∎

Corollary 2.4.

Let $H$ be a discrete one-dimensional Schrödinger operator (1.6). Suppose that $J\subset\mathbb{R}$ is an open interval such that $\mathbf{1}_{J}(H)e_{0}=0$. Let $\varphi$ be an eigenfunction of $H$ with the eigenvalue $\lambda\in J$, and let $\theta\in\mathrm{Ran}\,\mathbf{1}_{J\setminus\{\lambda\}}(H)$. Then $\langle\theta_{+},\varphi_{+}\rangle=0$.

Proof.

Follows from applying Lemma 2.3 to intervals of the form $[\lambda+1/n,+\infty)\cap J$ and $(-\infty,\lambda-1/n]\cap J$. ∎

Corollary 2.5.

Let $H$ be a discrete one-dimensional Schrödinger operator (1.6). Suppose that $J\subset\mathbb{R}$ is an open interval such that $\mathbf{1}_{J}(H)e_{0}=0$. Let $\varphi$, $\psi$ be two distinct normalized eigenfunctions with eigenvalues in $J$:

$$H\varphi=\lambda\varphi,\,\,H\psi=\mu\psi,\,\,\lambda,\mu,\in J,\,\,\lambda\neq\mu,\,\,\|\varphi\|=\|\psi\|=1,\,\,\varphi(0)=\psi(0)=0.$$

Define $m$ as in Lemma 2.1. Then

(2.5) $$\operatorname{{dist}}(\lambda,\partial J)\leq\frac{|\lambda-\mu|}{\sqrt{1-m^{2}}}.$$

As a consequence, an eigenvalue of $H$ contained in $J$ cannot be a limit point of other eigenvalues of $H$.

2.2. Quasiperiodic operators with $\gamma$-monotone potentials

In this section, we will apply Theorem 1.4 to quasiperiodic operators with $\gamma$-monotone potentials. In order to state the most general known results on localization for such operators, define

$$\beta(\alpha):=\limsup_{k\to+\infty}\frac{\log q_{k+1}}{q_{k}},$$

where $\frac{p_{k}}{q_{k}}$ is the sequence of continued fraction approximants to $\alpha$. In addition to $\gamma$-monotonicity, we will also require

(2.6) $$\int_{0}^{1}\log(1+|f(x)|)\,dx<+\infty,$$

in order for the Lyapunov exponent of the operator family $H(x)$ to exist. The following results, in its final form, is established in [14].

Proposition 2.6.

Assume that $f$ is $\gamma$-monotone and $\eqref{eq_log_integrable}$ holds. Suppose that $L(E)>\beta(\alpha)$. Then, every polynomially bounded solution of the eigenvalue equation

$$(H(x)\psi)(n)=E\psi(n),\quad n\in\mathbb{Z},$$

decays exponentially. As a consequence, the set

$$\{E\in\mathbb{R}\colon L(E)>\beta(\alpha)\}$$

can only support purely point spectrum of the operator family $\eqref{eq_h_def_qp}$.

Remark 2.7.

Without much detail, we will discuss several situations in which the above result is applicable. The function $f$ is always assumed to be $\gamma$-monotone.

1. $L(\cdot)$ is continuous in $E$ and one always has $L(E)\geq\log(\gamma/2)$. As a consequence, the results guarantee complete localization for large $\gamma$.

2. If $f$ is unbounded, then $L(E)>0$ for Lebesgue almost every $E$ [18]. The same conclusion holds if $f$ is bounded and only has finitely many discontinuities on each trajectory of the irrational rotation [6].

3. The integrated density of states is Lipschitz continuous. As a consequence, if the conclusion of the previous part holds, then for almost every $x$ the zero set of $L(E)$ does not contribute to spectral measure of $H(x)$. If $\beta(\alpha)=0$ (which is weaker than the Diophantine condition), then the previous claim implies localization for almost every $x$ for a class of $\gamma$-monotone potentials which includes sawtooth-type and Maryland-type.

4. Instead of $[0,1)$, one can also consider $\gamma$-monotone functions on $(0,1]$ (for example, by considering $-f(-x)$ instead of $f(x)$). As a consequence, each function from the family $f_{t}$ defined in (1.3) is $\gamma$-monotone in this sense.

5. In the case $f(0)=-\infty$, the operator $H(n\alpha)$ will have infinite coupling. In other words, the infinite value $V(n)=-\infty$ of the potential will enforce the Dirichlet condition $\psi(n)=0$, and the operator will split into two half-line operators. The conclusion of Proposition 2.6 still holds in this case, but the operator may not necessarily have simple spectrum. We discuss the case of infinite coupling in detail in Section 3.

In order to state the second main result, we will need to be somewhat careful about the dependence of $\sigma(H(x))$ on $x$ as a set. Let $f$ be $\gamma$-monotone. We will call $x_{0}\in[0,1)$ generic if, for every $y_{0}\in(x_{0}+\alpha\mathbb{Z})\,\,\mathrm{mod}\,\,1$, we have one of the following:

1. $y_{0}\in(0,1)$ and $f$ is continuous at $y_{0}$.

2. $y_{0}=0$ and $f(y_{0})=f(y_{0}+0)=-\infty$, $f(1-0)=+\infty$.

It is easy to see that $f$ is Maryland-type if and only if every point $x_{0}\in[0,1)$ is generic and that set of non-generic points of a $\gamma$-monotone function $f$ is at most countable.

Lemma 2.8.

Let $f$ be $\gamma$-monotone. Then, the following holds.

1. If $x,y\in[0,1)$ and $x$ is generic, then $\sigma(H(x))\subset\sigma(H(y))$.

2. As a consequence, $\sigma(H(x))=:\Sigma_{f}$ does not depend on $x$ for generic $x$.

3. For every $x\in[0,1)$ we have $\sigma(H(x+0))=\sigma(H(x-0))=\Sigma_{f}$.

4. For every $x\in[0,1)$, we have $\sigma_{\mathrm{ess}}(H(x))=\Sigma_{f}$.

Most of the claims of the lemma are known, especially if one avoids infinite coupling. We postpone the proof of the Lemma 2.8 until the end of Section 3, where these issues can be dealt with systematically. Note that the operator $H(x-0)$ (with the potential $V(n)=f(x+n\alpha-0)$) may not belong to the original family $\{H(x)\}_{x\in[0,1)}$. However, a natural way to extend the family is discussed in Section 4.

We are now ready to state the first main result of the paper. As mentioned earlier, the “if” part has already been established in [16] and is only included for completeness.

Theorem 2.9.

Assume that $f$ is $\gamma$-monotone and $\eqref{eq_log_integrable}$ holds. For a fixed $x\in[0,1)$, define $H(x)$ as in (1.1). Let $J\subset\mathbb{R}$ be an open energy interval such that $L(E)>\beta(\alpha)$ for all $E\in J$. Then, $\sigma(H(x))\cap J$ contains an interval if and only if $f$ is Maryland-type.

Proof.

Suppose that $f$ is not Maryland-type, and let $x_{0}$ be a non-generic point of $[0,1)$ with respect to $f$. Due to Lemma 2.8, it is sufficient to show that $\sigma(H(x_{0}))\cap J$ does not contain an interval. There exists

$$y_{0}=(x_{0}+m\alpha)\,\,\mathrm{mod}\,\,1\in[0,1)$$

such that one of the following holds

Case 1: $f(y_{0})<f(y_{0}+0)$. Let $I:=[f(y_{0}),f(y_{0}+0)]$ and consider

$$f_{t}(x):=\begin{cases}f(x),&x\in[0,1)\setminus\{y_{0}\};\\ t,&x=y_{0}.\end{cases}$$

for $t\in I$. Define $H_{t}(x)$ to be $H(x)$ with $f$ replaced by $f_{t}$. For every $t\in I$, $f_{t}$ is a $\gamma$-monotone function, and the corresponding operator family $H_{t}(\cdot)$ has the same Lyapunov exponent as $H(\cdot)$. Therefore, Proposition 2.6 is applicable, and $\sigma(H_{t}(x_{0}))\cap J$ are pure point. By Theorem 2.9, $\sigma(H_{t}(x_{0}))\cap J$ cannot contain intervals. For $t=f(y_{0})$, this implies that $\sigma(H(x_{0}))$ cannot contain an interval.

Case 2: $y_{0}\in(0,1)$, $f(y_{0}-0)<f(y_{0})$. One can repeat the construction from Case 1 with $I=[f(y_{0}-0),f(y_{0})]$. Alternatively, one can consider $g(x)=-f(-x)$ as will be later done in Case 4.

Case 3: $y_{0}=0$, $f(y_{0})>-\infty$. Similarly, one can repeat Case 1 with $I=[-\infty,f(y_{0})]$.

Case 4: $y_{0}=0$, $f(y_{0})=-\infty$, $f(1-0)<+\infty$. This case can be reduced to the previous case in two steps. First, replace $f(x)$ by $g(x):=-f(1-x)$. Denote the original operator by $H^{f}$, and the new operator by $H^{g}$. It is easy to see that $H^{f}(x)$ is unitarily equivalent to $-H^{g}(1-x)$, and $g$ is $\gamma$-monotone on $(0,1]$. Define the new function

$$u(x):=\begin{cases}g(x),&x\in(0,1);\\ -f(1-0),&x=0.\end{cases}$$

It is easy to see that $u(x)$ is now $\gamma$-monotone on $[0,1)$ and is within the class of operators considered in Case 2 (with $J$ replaced by $-J$). As a consequence, $\sigma(H^{u}(0))\cap(-J)$ cannot contain an interval. Since $H^{u}(0)$ is a finite rank perturbation of $H^{g}(1)$ and the latter is unitarily equivalent to $-H^{f}(0)$, we have that $\sigma(H^{f}(0))\cap J$ also cannot contain an interval. ∎

The following corollary implies Theorem 1.1, in view of Remark 2.7.

Corollary 2.10.

Let $f$ be $\gamma$-monotone, satisfying $\eqref{eq_log_integrable}$ and not of Maryland-type. Suppose that $L(E)>\beta(\alpha)$ for Lebesgue almost every $E\in\mathbb{R}$. Then $\sigma(H(x))$ does not contain any intervals. As a consequence, $\sigma_{\mathrm{ess}}(H(x))$ is a Cantor set, and $\sigma(H(x))$ is a Cantor set for generic $x$.

Proof.

As mentioned before ([13, 16, 14]), $L(E)$ is continuous. Therefore, one can apply Theorem 2.9 on each open interval of the set $\{E\colon L(E)>\beta(\alpha)\}$. Since each isolated point of $\sigma(H(x))$ is an eigenvalue of finite multiplicity, we have that $\sigma_{\mathrm{ess}}(H(x))$ is a closed set that does not contain any intervals or isolated points, and therefore is a Cantor set. ∎

Remark 2.11.

As mentioned above, for non-generic $x$ the spectrum $\sigma(H(x))$ may contain isolated eigenvalues in the gaps of the Cantor set $\Sigma_{f}=\sigma_{\mathrm{ess}}(H(x))$. This would typically happen if, for some $x$, one has $f(x-0)\neq f(x)\neq f(x+0)$ (where one assumes $+\infty=-\infty$). Theorem 4.2 provides a large class of such examples.

3. Operators with infinite potentials

Theorem 2.9 partially relies on Lemma 2.8. The proof of this lemma is fairly standard, but requires certain care when considering strong limits as the values of the potential approach infinity. A similar issue appears when one considers gap filling such as in Theorem 1.2: an operator with infinite coupling belongs to a family of rank one perturbations that are used to fill the gaps. In the proof of the gap filling result, we will need to use rational approximations and obtain several statements about continuity of spectra, based on the same ideas as in [3] and later developments, some of which were implemented in [16]. In the case of quasiperiodic operators with discontinuous/unbounded potentials, the usual issue is lack of compactness of the set of operators under consideration, either due to multiple ways of approaching a discontinuity in the limit, or due to the value of the potential approaching infinity.

In order to take full advantage of the above arguments, we will extend the class of operators under consideration, making the parameter space compact in the natural topology. This will include allowing the potential to take infinite values. For the convenience of the reader, we discuss the relevant theory (which we do not believe to be completely novel) in this section.

In order to avoid overusing of the word “generalized”, the following conventions will be used without further clarification in this and later sections.

• The potential of a Schrödinger operator can take values in $\overline{\mathbb{R}}=\mathbb{R}\cup\{\infty\}$, which is identified with a circle. The positive orientation of the circle is inherited from the positive direction on $\mathbb{R}$. There will be no distinction between $V(n)=+\infty$ and $V(n)=-\infty$.

• The main topology for these operators will be the $\overline{\mathbb{R}}$-strong topology. Following Proposition 3.1 and Corollary 3.2, it is equivalent to pointwise convergence of potentials in $\overline{\mathbb{R}}$.

• The (generalized) spectrum of a Schrodinger operator $H$ is denoted by $\bar{\sigma}(H)$ and is a closed subset of $\overline{\mathbb{R}}$. We will always have $\sigma(H)=\bar{\sigma}(H)\cap\mathbb{R}$.

3.1. Potentials with infinite values

Let $\overline{\mathbb{R}}=\mathbb{R}\cup\{\infty\}$ be the extended real line and $V\colon\mathbb{Z}^{d}\to\overline{\mathbb{R}}$. Denote

$$I(V):=\{\mathbf{n}\in\mathbb{Z}^{d}\colon V(\mathbf{n})=\infty\}.$$

The discrete Schrödinger operator with potential $V$ is defined by the expression

(3.1) $$(H\psi)(\mathbf{n})=\mathbf{1}_{\mathbb{Z}^{d}\setminus I(V)}(\mathbf{n})(\Delta\psi)(\mathbf{n})+V(\mathbf{n})\psi(\mathbf{n}),\quad\mathbf{n}\in\mathbb{Z}^{d}.$$

on the domain

$$\operatorname{{Dom}}(H)=\left\{\psi\in\ell^{2}(\mathbb{Z})\colon\sum_{\mathbf{n}\in\mathbb{Z}}|V(\mathbf{n})\psi(\mathbf{n})|^{2}<+\infty\right\}.$$

Here we assume $0\cdot\infty=0$. The closure of $\operatorname{{Dom}}(H)$ in $\ell^{2}(\mathbb{Z}^{d})$ is $\ell^{2}(\mathbb{Z}^{d}\setminus I(V))$. It is convenient to consider $\ell^{2}(I(V))$ as the eigenspace of $H$ corresponding to the eigenvalue $\infty$, which makes the spectral measure of $H$ a measure on $\overline{\mathbb{R}}$ rather than $\mathbb{R}$. We denote

$$\bar{\sigma}(H):=\begin{cases}\overline{\sigma(H)},&I(V)=\varnothing;\\ \sigma(H)\cup\{\infty\},&I(V)\neq\varnothing,\end{cases}$$

where $\overline{\sigma(H)}$ means the closure in $\overline{\mathbb{R}}$. Equivalently, one can consider the (generalized) Cayley transform

$$U_{H}:=i{(H-i\mathbf{1}_{\mathbb{Z}^{d}\setminus I(V)})(H+i\mathbf{1}_{\mathbb{Z}^{d}\setminus I(V)})^{-1}}\oplus i\mathbf{1}_{I(V)}$$

and define the spectral measure of $H$ on $\overline{\mathbb{R}}$ as the pullback of the spectral measure of $U$ under the map $t\mapsto i\frac{t-i}{t+i}$. As a consequence, the inclusion of $\infty$ into the spectral theory of $H$ can be done by the means of the usual spectral theory of $U_{H}$.

Denote by $\mathcal{H}$ the set of all operators $H$ of the form (3.1), and let

$$\mathcal{U}_{H}:=\{U_{H}\colon H\in\mathcal{H}\}$$

be the set of their Cayley transforms. Our goal is to relate pointwise convergence of the potentials with convergence of the spectra of the corresponding Schrödinger operators. Since the spectra can be completely described in terms of $U_{H}$, this can be done by the means of strong operator topology on the set $\mathcal{U}_{H}$. The key property is the following proposition.

Proposition 3.1.

Let $H_{j}=\Delta+V_{j}$, $H=\Delta+V$. Suppose that $V_{j}(\mathbf{n})$ converges to $V(\mathbf{n})$ in $\overline{\mathbb{R}}$ for every $\mathbf{n}\in\mathbb{Z}$. Then $U_{H_{j}}$ converges to $U_{H}$ strongly.

Proof.

Instead of Cayley transforms, one can equivalently work with (generalized) resolvents. If $H\in\mathcal{H}$, we will use the notation $(H+i\mathbf{1})^{-1}$ to denote the resolvent of $H$, which is originally a bounded operator on $\ell^{2}(\mathbb{Z}^{d}\setminus I(V))$, extended by zero into $\ell^{2}(I(V))$:

$$(H+i\mathbf{1})^{-1}:=(H+i\mathbf{1}_{\mathbb{Z}^{d}\setminus I(V)})^{-1}\oplus\mathbf{0}_{I(V)}.$$

With this convention, we have

$$U_{H}=i\mathbf{1}+2(H+i\mathbf{1})^{-1}.$$

As a consequence, the strong convergence of $U_{H_{j}}$ is equivalent to the strong convergence of the (generalized) resolvents defined above. It is easy to check that the projected resolvent identity

(3.2) $$(H_{j}+i\mathbf{1})^{-1}\mathbf{1}_{\mathbb{Z}^{d}\setminus I(V)}-\mathbf{1}_{\mathbb{Z}^{d}\setminus I(V_{j})}(H+i\mathbf{1})^{-1}=(H_{j}+i\mathbf{1})^{-1}(V-V_{j})(H+i\mathbf{1})^{-1}$$

holds with the convention $0\cdot\infty=0$, which eliminates all entries of $V(\mathbf{n})-V_{j}(\mathbf{n})$ with $\mathbf{n}\in I(V_{j})\cup I(V)$.

Suppose that $V_{j}\to V$ pointwise on $\mathbb{Z}^{d}$. Our goal is to show strong convergence of $U_{H_{j}}$ or, equivalently, $(H_{j}+i\mathbf{1})^{-1}$. Since the resolvents are uniformly bounded, it is sufficient to check the convergence on a dense set. Suppose first that $e_{\mathbf{n}}$ is a standard basis vector with $\mathbf{n}\in I(V)\setminus I(V_{j})$. Since $V_{j}(\mathbf{n})\to\infty$, we can assume that, eventually, $V_{j}(\mathbf{n})\neq 0$, and write

$$(H_{j}+i\mathbf{1})^{-1}e_{\mathbf{n}}=\frac{1}{V_{j}(\mathbf{n})}(H_{j}+i\mathbf{1})^{-1}\left((H_{j}+i\mathbf{1})e_{\mathbf{n}}-ie_{\mathbf{n}}-\Delta e_{\mathbf{n}}\right)\to 0\,\,\text{as}\,\,j\to+\infty.$$

For $\mathbf{n}\in I(V_{j})$, one already has $(H_{j}+i\mathbf{1})^{-1}e_{\mathbf{n}}=0$. By linearity, strong convergence is now established on a dense subset of $\ell^{2}(I(V))$.

For $\mathbf{n}\in\mathbb{Z}^{d}\setminus I(V)$, let $\psi_{\mathbf{n}}:=(H+iI)e_{\mathbf{n}}$. Since linear combinations of $e_{\mathbf{n}}$ are dense in $\operatorname{{Dom}}H$ both in $\ell^{2}$ and $H$-norms, we have that the linear combinations of $\psi_{\mathbf{n}}$ are dense in $\ell^{2}(\mathbb{Z}^{d}\setminus I(V))$. By applying (3.2) to $\psi_{\mathbf{n}}$ and using $(V-V_{j})e_{\mathbf{n}}\to 0$ and that the resolvents are bounded, we arrive to convergence of $(H_{j}+i\mathbf{1})^{-1}$ on a dense subset of $\ell^{2}(\mathbb{Z}^{d}\setminus I(V))$, which completes the proof. ∎

The standard diagonal procedure argument immediately applies the converse statement.

Corollary 3.2.

The set $\mathcal{H}$ is compact in the topology of pointwise $\overline{\mathbb{R}}$ convergence of potentials from Proposition 3.1. As a consequence, the set $\mathcal{U}$ is also compact in the strong operator topology, and the converse of Proposition 3.1 also holds: strong topology on $\mathcal{U}$ induces pointwise $\overline{\mathbb{R}}$ convergence for potentials on $\mathcal{H}$.

We will refer to this convergence of (generalized) Schrödinger operators on $\mathcal{H}$ as the $\overline{\mathbb{R}}$-strong convergence. Finally, the following is the standard result of semicontinuity of the spectra.

Corollary 3.3.

Suppose that $H_{j}\to H$ $\overline{\mathbb{R}}$-strongly, and $E\in\bar{\sigma}(H)$. Then there exists a sequence $E_{j}\in\bar{\sigma}(H_{j})$ such that $E_{j}\to E$ in $\overline{\mathbb{R}}$.

As in the usual setting, the proof can be performed by contradiction, using the resolvents of $U_{H}$ and $U_{H_{j}}$ near the point $i\frac{E-i}{E+i}$. Corollary 3.3, as is well known, states that new spectrum cannot appear in strong limit. A converse statement is not necessarily true, but the following weaker claim is useful, see also [3].

Proposition 3.4.

Suppose that $H_{j}\to H$ $\overline{\mathbb{R}}$-strongly. Assume also that a sequence $\psi_{j}\colon\mathbb{Z}^{d}\to[-1,1]$ converges pointwise to $\psi\colon\mathbb{Z}^{d}\to[-1,1]$ and satisfies a sequence of eigenvalue equations

$$H_{j}\psi_{j}=E_{j}\psi_{j},\quad E_{j}\in\overline{\mathbb{R}},\quad E_{j}\to E\in\overline{\mathbb{R}},\quad\psi_{j}(0)\geq 1/2.$$

Then $E\in\bar{\sigma}(H)$.

Proof.

Suppose first that $E_{j}\to\infty$. Then $\psi_{j}(\mathbf{n})\to 0$ for all $\mathbf{n}\in\mathbb{Z}^{d}\setminus I(V)$. Since $\psi_{j}(0)\geq 1/2$, we have that $0\in I(V)$, and therefore $\infty\in\bar{\sigma}(H)$.

It remains to consider the case $E\in\mathbb{R}$. In this case, $\psi_{j}(\mathbf{n})\to 0$ for all $\mathbf{n}\in I(V)$, and one can pass to the limit in the eigenvalue equation, thus obtaining $H\psi=E\psi$. By the standard Weyl sequence argument, we have that $\psi\in\bar{\sigma}(H)$. ∎

We will now discuss some results related to compactness.

Lemma 3.5.

Let $S\subset\mathbb{Z}^{d}$, and suppose that $|V(\mathbf{n})|\geq M$ for all $\mathbf{n}\in S$. Let

$$V_{\infty}:=\begin{cases}\infty,&\mathbf{n}\in S\\ V(\mathbf{n}),&\mathbf{n}\in\mathbb{Z}^{d}\setminus S,\end{cases}$$

and define $H_{\infty}$ using (3.1) with $V$ replaced by $V_{\infty}$. Then

$$\|U_{H}-U_{H_{\infty}}\|\leq\frac{2+2^{d}}{M}.$$

Proof.

Recall the projected resolvent identity (3.2):

$$R:=(H_{\infty}+i\mathbf{1})^{-1}\mathbf{1}_{\mathbb{Z}^{d}\setminus I(V)}-\mathbf{1}_{\mathbb{Z}^{d}\setminus I(V_{\infty})}(H+i\mathbf{1})^{-1}=(H_{\infty}+i\mathbf{1})^{-1}(V-V_{\infty})(H+i\mathbf{1})^{-1}=0,$$

since $V=V_{\infty}$ on $\mathbb{Z}^{d}\setminus(I(V)\cup I(V_{\infty}))=\mathbb{Z}^{d}\setminus I(V_{\infty})$. Moreover, both resolvents vanish on $\ell^{2}(I(V))$. Therefore,

$$\|U_{H_{\infty}}-U_{H}\|=\|R-\mathbf{1}_{I(V_{\infty})\setminus I(V)}(H+i\mathbf{1})^{-1}\mathbf{1}_{I(V_{\infty})\setminus I(V)}\|=\|\mathbf{1}_{I(V_{\infty})\setminus I(V)}(H+i\mathbf{1})^{-1}\mathbf{1}_{I(V_{\infty})\setminus I(V)}\|.$$

Suppose that $\operatorname{{supp}}\psi\subset I(V_{\infty})\setminus I(V)$, and $\varphi(\mathbf{n}):=V^{-1}(\mathbf{n})\psi(\mathbf{n})$. Then

$$\|(H+i\mathbf{1})^{-1}\psi\|=\|(H+i\mathbf{1})^{-1}\left((H+i\mathbf{1})\varphi-i\varphi-\Delta\varphi\right)\|\leq(2+2^{d})\|\varphi\|\leq\frac{2+2^{d}}{M}\|\psi\|,$$

since $|V(\mathbf{n})|\geq M$ on $\operatorname{{supp}}\psi$ and $\|(H+i\mathbf{1})^{-1}\|\leq 1$. ∎

Lemma 3.6.

Let $V_{1},V_{2}\colon\mathbb{Z}^{d}\to\overline{\mathbb{R}}$ be such that $V_{1}(\mathbf{n})=V_{2}(\mathbf{n})$ for all but finitely many $\mathbf{n}\in\mathbb{Z}^{d}$. Define $H_{1}$, $H_{2}$ as in (3.1) with $V_{1}$, $V_{2}$, respectively. Then $U_{H_{1}}-U_{H_{2}}$ has finite rank.

Proof.

In the resolvent identity

$$(H_{1}+i\mathbf{1})^{-1}\mathbf{1}_{\mathbb{Z}^{d}\setminus I(V_{2})}-\mathbf{1}_{\mathbb{Z}^{d}\setminus I(V_{1})}(H_{2}+i\mathbf{1})^{-1}=(H_{1}+i\mathbf{1})^{-1}(V_{2}-V_{1})(H_{2}+i\mathbf{1})^{-1}$$

the right hand side is of finite rank. Since $I(V_{1})$ and $I(V_{2})$ can differ by at most finitely many points, the left hand side is also a finite rank perturbation of $U_{H_{1}}-U_{H_{2}}$. ∎

Clearly, if $V_{1}$ and $V_{2}$ are two finite potentials on $\mathbb{Z}^{d}$ such that $|V_{1}(\mathbf{n})-V_{2}(\mathbf{n})|\to 0$ as $|\mathbf{n}|\to+\infty$, then $H_{1}-H_{2}=(\Delta+V_{1})-(\Delta+V_{2})$ is compact. In the infinite potential case, in addition to allowing the values of the potentials approach one another, one can also allow both of them approach the infinity.

In order to make it precise, let $V_{1},V_{2}\colon\mathbb{Z}^{d}\to\overline{\mathbb{R}}$. We will say that $V_{1}\sim V_{2}$ if for every $\varepsilon>0$ and $M>0$ there exists $N=N(\varepsilon,M)$ such that $|V_{1}(\mathbf{n})-V_{2}(\mathbf{n})|<\varepsilon$ or $|V_{1}(\mathbf{n})|>M$ and $|V_{2}(\mathbf{n})|>M$ for all $\mathbf{n}$ with $|\mathbf{n}|>N$.

Lemma 3.7.

Let $V_{1},V_{2}\colon\mathbb{Z}^{d}\to\overline{\mathbb{R}}$. Define $H_{1}$, $H_{2}$ as in (3.1) with $V_{1}$, $V_{2}$, respectively. Suppose that $V_{1}\sim V_{2}$. Then $U_{H_{1}}-U_{H_{2}}$ is compact.

Proof.

Fix $\varepsilon,M>0$ and find $N$ using the fact that $V_{1}\sim V_{2}$. Let

$$V_{j}^{\prime}(\mathbf{n}):=\begin{cases}0,&|\mathbf{n}|\leq N\\ \infty,&|\mathbf{n}|>N\,\,\text{and}\,\,|V_{j}(\mathbf{n})|>M\\ V_{j}(\mathbf{n}),&|\mathbf{n}|>N\,\,\text{and}\,|\,V_{j}(\mathbf{n})|\leq M,\end{cases}$$

and define $H_{j}^{\prime}$, $j=1,2$, accordingly. From Lemmas 3.5, 3.6, and several applications of the triangle inequality, we arrive to

$$U_{H_{1}}-U_{H_{2}}=(U_{H_{1}^{\prime}}-U_{H_{2}^{\prime}})+A_{\varepsilon,M}+B_{\varepsilon,M},$$

where $A_{\varepsilon,M}$ has finite rank, and, say $\|B_{\varepsilon,M}\|\leq\frac{2^{d+2}}{M}$. On the other hand, since $I(V_{1}^{\prime})=I(V_{2}^{\prime})$, the resolvent identity implies that $\|U_{H_{1}^{\prime}}-U_{H_{2}^{\prime}}\|\leq\varepsilon$. We thus have obtained a sequence of finite rank approximations converging to $U_{H_{1}}-U_{H_{2}}$ in the operator norm, which completes the proof. ∎

Corollary 3.8.

Let $f_{0}\colon[0,1)\to[-\infty,+\infty)$ be $\gamma$-monotone, and define

$$f_{1}(x):=f_{0}(x+0),\quad x\in[0,1),$$

(3.3) $$(H_{f_{j}}(x)\psi)(\mathbf{n})=(\Delta\psi)(\mathbf{n})+f_{j}(x+\mathbf{n}\cdot\alpha)\psi(\mathbf{n}),\quad x\in[0,1),\quad\mathbf{n}\in\mathbb{Z}^{d}$$

with the conventions adopted above for infinite values of $f_{j}$. Then the difference

$$U_{H_{f_{0}}(x)}-U_{H_{f_{1}}(x)}$$

is compact, and therefore $\sigma_{\mathrm{ess}}(H_{f_{0}}(x))=\sigma_{\mathrm{ess}}(H_{f_{1}}(x))$.

Proof.

From $\gamma$-monotonicity, it is easy to see that the change in the potential between $f_{0}$ and $f_{1}$ is within the assumptions of Lemma 3.7. ∎

Proof of Lemma 2.8. Note that all of the claims hold for operators on $\mathbb{Z}^{d}$ considered in Corollary 3.8. The first claim follows from the fact that, for generic $x$ and any $y$, there is a sequence $\mathbf{n}_{j}$ such that $H(y+\mathbf{n}_{j}\cdot\alpha)\to H(x)$ $\overline{\mathbb{R}}$-strongly. The second claim follows from the first one. This argument can also be applied to every $x$ if $f$ is either left or right continuous, which implies the third claim. Finally, Corollary 3.8 implies the fourth claim.

4. Structure of the hull, rational approximations, and gap filling

4.1. General properties of gap filling

In the natural topology, $\overline{\mathbb{R}}$ is homeomorphic to a circle. We will assume that positive orientation is inherited from that of $\mathbb{R}$. For $t_{1},t_{2}\in\overline{\mathbb{R}}$, a positively-oriented arc that starts at $t_{1}$ and ends at $t_{2}$, will be denoted by $(t_{1},t_{2})$ and called an interval. One can also consider closed intervals $[t_{1},t_{2}]$. For example, one has have

$$(1,-2]=\overline{\mathbb{R}}\setminus(-2,1]=(1,+\infty)\cup\{\infty\}\cup(-\infty,-2];$$

$$\quad(0,0)=\overline{\mathbb{R}}\setminus\{0\},\quad[\infty,1)=\{\infty\}\cup(-\infty,1).$$

If $\mathcal{A}$ is an arc and $f\colon\mathcal{A}\to\overline{\mathbb{R}}$ is continuous, we will say that $f$ is monotone if it is locally monotone in coordinate charts preserving orientation (essentially, it corresponds to a usual definition of a monotone function on $\mathbb{R}$, with a convention that passing through $\infty$ happens in the correct direction). If both $\mathcal{A}$ and $\mathcal{B}$ are arcs not equal to $\overline{\mathbb{R}}$, then $f\colon\mathcal{A}\to\mathcal{B}$ (not necessarily continuous) will be called monotone if it preserves ordering on the arcs induced by the orientation. Note that, for a fixed $f$, whether or not it is monotone may depend on the choice of the arcs $\mathcal{A}$ and $\mathcal{B}$, even if all these choices are of the same orientation.

Let $A$ be a self-adjoint operator and $P$ be a rank one projection whose range is cyclic for $A$. Suppose, for simplicity, that $\sigma(A)=\sigma_{\mathrm{ess}}(A)$, and let

$$A_{t}:=A+tP.$$

In applications, $A$ will be a Schrödinger operator, and $P$ a projection onto one of the standard basis vectors. Therefore, we will be able to also consider $t=\infty$, by defining

$$A_{\infty}:=\left.(1-P)A(1-P)\right|_{\operatorname{{ran}}(1-P)},$$

in accordance with the conventions in Section 3. Let

$$G_{j}=(g_{j}^{-},g_{j}^{+})\subset\overline{\mathbb{R}}$$

be a spectral gap of $A$, including, possibly, one infinite gap considered as a subset of $\overline{\mathbb{R}}$ (in other words, $(g_{j}^{-},g_{j}^{+})$ are the connected components of the complement of the closure of $\sigma(A)$ in $\overline{\mathbb{R}}$). The general theory of rank one perturbations (see, for example, Chapter 11 of [17]) implies that $\sigma(A_{t})$ with $t\in\overline{\mathbb{R}}\setminus\{0\}$ contains $\sigma(A)$ as well as at most one simple eigenvalue in each $(g_{j}^{-},g_{j}^{+})$. This eigenvalue is described by a continuous monotone function $\lambda_{j}\colon(t_{j}^{-},t_{j}^{+})\to\overline{\mathbb{R}}$, where $(t_{j}^{-},t_{j}^{+})\subset\overline{\mathbb{R}}\setminus\{0\}$ is a (maximal open) interval on which $(g_{j}^{-},g_{j}^{+})$ contains an eigenvalue: that is,

$$\lambda_{j}(t_{j}^{\pm}\mp 0)=g_{j}^{\pm}.$$

The functions $\lambda_{j}$ are strictly monotone in $t$ (in the sense of maps between arcs in $\overline{\mathbb{R}}$). We also have no eigenvalues in $G_{j}$ for $t\in\overline{\mathbb{R}}\setminus(t_{j}^{-},t_{j}^{+})$. For $t=t_{j}^{\pm}$, the operator $A_{t}$ may or may not have an eigenvalue at the corresponding endpoint $g_{j}^{\pm}$. In the above notation, we also mean that $A_{\infty}$ has a simple eigenvalue $\lambda=\infty$ with the eigenspace $\operatorname{{ran}}P$, following the conventions of Section 3.

To summarize the above, as $t$ runs over the interval $(0,0)=(0,+\infty)\cup\{\infty\}\cup(-\infty,0)$, each gap will eventually have one eigenvalue appear at one of its endpoints, move continuously and monotonically along the gap, and then disappear at the other endpoint. Eigenvalues in different gaps may appear and disappear at different times, but no eigenvalue can re-appear in the same gap. The starting and ending point of the range for $t$ being $t=0$ corresponds to the fact that we originally assume that $A$ has no isolated eigenvalues.

4.2. Statements of the results

We will consider Schrödinger operators associated to $1$-periodic maps $f\colon\mathbb{R}\to[-\infty,+\infty)$:

(4.1) $$(H(x)\psi)(\mathbf{n})=(\Delta\psi)(\mathbf{n})+f(x+\mathbf{n}\cdot\alpha)\psi(\mathbf{n}),\quad x\in[0,1),\quad\mathbf{n}\in\mathbb{Z}^{d},$$

where the operator is considered in the sense of (3.1) for infinite values of $f$. Both $[0,1)$ and $[-\infty,+\infty)$ can be identified with circles, making $f|_{[0,1)}$ a circle map. We will consider a class of $f$ such that the corresponding circle maps have one discontinuity (which, without loss of generality, can be placed at the origin).

Definition 4.1.

We say that a $1$-periodic function $f\colon\mathbb{R}\to[-\infty,+\infty)$ has a simple discontinuity if $f|_{[0,1)}$ extends to a continuous map from from $[0,1]$ to $[-\infty,+\infty]$, with $f(0)\neq f(1-0)$ and at least one value among $f(0)$ and $f(1-0)$ being finite.

Every sawtooth-type potential defined in the introduction has a simple discontinuity. The Maryland-type ones do not, since the corresponding circle maps will be continuous. Let $\mathcal{A}$ be a closed arc in $\overline{\mathbb{R}}$ (which is associated to $[-\infty,+\infty)$) whose endpoints are $f(0)$ and $f(1-0)$. Note that there are two choices of such arc (more on this later). Our goal will be to consider a family of operators, parametrized by $t\in\mathcal{A}$, where the value $f(0)$ will be replaced by $t$. More precisely, for $t\in\mathcal{A}$, define

$$f_{t}(x):=\begin{cases}f(x),&x\in(0,1);\\ t,&x=0.\end{cases}$$

Our goal is to state a gap filling result of the kind described in Theorem 1.2. In the case of a sawtooth-type potential, say, $f(x)=\{x\}$, there is an obvious choice between two arcs (the points of the arc should complement the range of $f$). In general, since $f$ is not assumed to be monotone, we will need to specify it more carefully. For an arc $\mathcal{A}$ with endpoints $f(0)$ and $f(1-0)$, let

(4.2) $$\varphi_{\mathcal{A}}\colon[-1,0]\to\mathcal{A},\quad\varphi_{\mathcal{A}}(-1)=f(1-0),\quad\varphi_{\mathcal{A}}(0)=f(0)$$

be a continuous parametrization of $\mathcal{A}$. Consider a new map $\tilde{f}\colon[-1,1)\to[-\infty,+\infty)$ defined by

(4.3) $$\tilde{f}(y):=\begin{cases}f(y),&x\in[0,1);\\ \varphi_{\mathcal{A}}(y),&y\in[-1,0).\end{cases}$$

Clearly, if one now identifies $[-1,1)$ and $[-\infty,+\infty)$ with circles, $\tilde{f}$ will become a continuous circle map. We are now ready to state the gap filling result.

Theorem 4.2.

Suppose that $f$ has a simple discontinuity. Define $H_{t}(x)$ by $\eqref{eq_h_def_qp_zd}$ with $f$ replaced by $f_{t}$, $t\in\mathcal{A}$. Suppose that the associated circle map $\tilde{f}$, defined above, is not homotopic to a constant. Then

(4.4) $$\bigcup_{t\in\mathcal{A}}\bar{\sigma}(H_{t}(0))=\overline{\mathbb{R}}.$$

Note that, for every $f$ with a simple discontinuity, at least one choice of the arc $\mathcal{A}$ will produce $\tilde{f}$ not homotopic to a constant. As described in the previous section, (4.4) holds for any family of rank one perturbations with $\mathcal{A}$ replaced by $\overline{\mathbb{R}}$ ([17, Theorem 11.8]).

The following can also be obtained as a corollary of Proposition 2.6, assuming that, in addition to the above, the operator is in the complete spectral localization regime.

Corollary 4.3.

Suppose that $f$ and $\mathcal{A}$ satisfy the assumptions of Theorem 4.2 and Proposition 2.6. Suppose also that $L(E)>\beta(\alpha)$ for all $E\in\mathbb{R}$. Then, for every $E\in\overline{\mathbb{R}}$, there exists $x\in[0,1)$ and $t\in\mathcal{A}$ such that $E$ is an eigenvalue of $H_{t}(x)$.

The proof of Theorem 4.2 is based on the fact that a version of gap filling holds for rational approximations of the operator. In the process, we use a version of continuity of the spectra under rational approximations, extending the known results of [3]. Instead of considering separate cases as in [16], we use a systematic construction of a compactification of the hull that involves operators with infinite potentials. The fact that gap filling holds for rational frequencies uses a “topological” argument: if $u\colon\mathbb{T}\to U(n)$ is a continuous family of unitary matrices that is not homotopic to a constant, then the union of their spectra must be equal to $\mathbb{T}$. As a consequence, one also has the following “topological” version of the result of [16] for Maryland-type potentials, which also works in every dimension.

Theorem 4.4.

Suppose that $f\colon\mathbb{R}\to[-\infty,+\infty)$ is $1$-periodic, and the circle map associated to $f|_{[0,1)}$ is continuous and not homotopic to a constant. Define $H(x)$ as in (4.1). Then $\sigma(H(x))=\mathbb{R}$.

4.3. Structure of the hull

The proof of Theorem 4.2 will be based on rational approximations, which requires some kind of continuity of spectra. The usual arguments of [3] require compactness of the parameter space, which is not present in the case of potentials with simple discontinuity, for two reasons: unboundedness of the potentials and the presence of discontinuity. The constructions of Section 3 allow to consider unbounded potentials, since $\mathcal{U}_{H}$ is compact. However, the presence of discontinuity still requires to extend the parameter space: for example, $H(\mathbf{n}\cdot\alpha-0)$ does not necessarily belong to the original operator family. Finally, the role of the function $\tilde{f}$ in Theorem 4.2 suggests that the extended operator family (that is, the one that includes $H_{t}(0)$ and all it’s translations) can also play a role in the proof.

In this section, we will construct a family of operator that will extend the original family $\{H(x)\colon x\in[0,1)\}$ and serve the purposes described above. It will be convenient to construct an auxiliary Cantor set whose gaps are labeled by $\{\mathbf{n}\cdot\alpha\}$, $\mathbf{n}\in\mathbb{Z}^{d}$. Since we will be dealing only with aspects of topology, the exact construction of such a set would not matter (but some maps will be needed to be fixed in a way continuous with respect to $\alpha$).

Let $h\colon[0,1]\to[0,2]$ be defined as follows:

$$h(y):=y+3^{-d}\sum_{\mathbf{n}\colon\{\mathbf{n}\cdot\alpha\}\leq y}2^{-|\mathbf{n}|_{1}}.$$

It will be convenient to consider the range $[0,2]$ of $h$ as a subset of the circle that is identified with, say, $[0,3)$. Denote that circle by $\mathcal{S}$, so that $h\colon[0,1]\to\mathcal{S}$. Let also $\mathcal{C}_{\alpha}$ be the closure of the range of $h$ in $\mathcal{S}$. We will fix the orientation on $\mathcal{S}$ so that $h\colon[0,1]\to[0,2]$ is monotone (here, $[0,1]$ and $[0,2]$ are considered as arcs of $\mathcal{S}$).

The set $\mathcal{C}_{\alpha}$ is a Cantor set if and only if at least one component of $\alpha$ is irrational. In all cases, the gaps or $\mathcal{C}_{\alpha}$ (that is, connected components of $\mathcal{S}\setminus\mathcal{C}_{\alpha}$) are in one-to-one correspondence with the numbers $\{\mathbf{n}\cdot\alpha\}$. Let $\mathcal{G}_{\mathbf{n}}$ be the gap associated to $\mathbf{n}\in\mathbb{Z}$. Note that, if the components of $\alpha$ are rationally independent, different values of $\mathbf{n}$ may correspond to the same value of $\{\mathbf{n}\cdot\alpha\}$ and therefore to the same gap.

Lemma 4.5.

Suppose that $f$ has a simple discontinuity. Define $H(x)$ as in $\eqref{eq_h_def_qp_zd}$. Then the closure of the set

(4.5) $$\{H(x)\colon x\in[0,1)\}$$

in the $\overline{\mathbb{R}}$-strong topology is homeomorphic to $\mathcal{C}_{\alpha}$. The homeomorphism can be defined on the dense subset

(4.6) $$\{H(x)\colon x\in[0,1)\setminus(\mathbb{Z}^{d}\cdot\alpha)\}$$

by the map $H(x)\mapsto h(x)$, and is extended to the closure by continuity.

Proof.

The closure of (4.5) can be described by

$$\{H(x\pm 0)\colon x\in[0,1)\},$$

where one has $H(x)=H(x+0)$ for all $x$ and $H(x)=H(x-0)$ for $x\in[0,1)\setminus\mathbb{Z}^{d}\cdot\alpha$. For $s,s_{j}\in\{+0,-0\}$, it is easy to see that $H(x_{j}s_{j})\to H(xs)$ if and only if $x_{j}\to xs$ with $x_{j}\neq x$ for $j$ large enough, or $x_{j}s_{j}=xs$ for $j$ large enough. The latter is equivalent to the topology on $\mathcal{C}_{\alpha}$. ∎

The $\mathbb{Z}^{d}$-action $x\mapsto x+\mathbf{n}\cdot\alpha$ can be extended to $\mathcal{C}_{\alpha}$ by continuity. For $\theta\in\mathcal{C}_{\alpha}$, we will denote this action in the same way: $\theta\mapsto\theta+\mathbf{n}\cdot\alpha$. Note that, even though $\mathcal{C}_{\alpha}$ is a subset of a circle, the action is not a circle rotation (and, in general, cannot be continuously conjugated to one). For $\theta\in\mathcal{C}_{\alpha}$, we will define the corresponding element of the closure of (4.6) by $H(\theta)$.

We are deliberately using the same notation for $H(x)$ and $H(\theta)$, meaning the following: the set $[0,1)$ of values of $x$ is naturally identified with the set $\{H(x+0)\colon x\in[0,1)\}$ which is identified with a subset of $\mathcal{C}_{\alpha}$ of values of $\theta$.

Proposition 4.6.

Suppose that $\alpha$ has at least one irrational component. There exists $\Sigma\subset\overline{\mathbb{R}}$ such that $\bar{\sigma}(H(\theta))=\Sigma$ for all $\theta\in\mathcal{C}_{\alpha}$.

Proof.

Follows in a standard way from Proposition 3.1, since the trajectories of the action are dense in $\mathcal{C}_{\alpha}$ (one can apply the usual proof to the family $\{U_{H(\theta)}\colon\theta\in\mathcal{C}_{\alpha}\}$ in order to avoid dealing with infinite potentials). ∎

4.4. Rank one perturbations and the extended operator family

As in the previous section, suppose that $f\colon[0,1)\to\overline{\mathbb{R}}$ has a simple discontinuity. Assume that $\mathcal{A}$ is an arc with endpoints $f(0),f(1-0)$. As in the statement of Theorem 4.2, let

$$f_{t}(x):=\begin{cases}f(x),&x\in(0,1);\\ t,&x=0,\end{cases}$$

for $t\in\mathcal{A}$. Let also, as before, $H_{t}(x)$ be the operator (4.1) with $f$ replaced by $f_{t}$.

Recall that

$$\varphi_{\mathcal{A}}\colon[-1,0]\to\mathcal{A},\quad\varphi_{\mathcal{A}}(-1)=f(1-0),\quad\varphi_{\mathcal{A}}(0)=f(0),$$

defined in (4.2), is a continuous parametrization of $\mathcal{A}$. Note that the change of the value $f(0)$ by $t$ only affects the operators $H(x)$ with $x=\mathbf{n}\cdot\alpha\pm 0$ and interpolates between $H(\mathbf{n}\cdot\alpha-0)$ and $H(\mathbf{n}\cdot\alpha+0)=H(\mathbf{n}\cdot\alpha)$. These two operators correspond to the endpoints of the gap $\mathcal{G}_{\mathbf{n}}$ of $\mathcal{C}_{\alpha}$. We will use the points of $\mathcal{G}_{\mathbf{n}}$ to parametrize this interpolation. More precisely, for each $\mathbf{n}$, let

$$\varphi_{\mathbf{n}}\colon[-1,0]\to\overline{\mathcal{G}_{\mathbf{n}}}$$

be a monotone parametrization of the gap $\mathcal{G}_{\mathbf{n}}$ the appropriately rescaled arc length (where the orientation and arc length on $\mathcal{S}$ are determined from its identification with $[0,3)$ as done in the definition). For $\theta\in\mathcal{G}_{\mathbf{n}}$, define

$$H(\theta):=H_{\varphi_{\mathcal{A}}(\varphi_{\mathbf{n}}^{-1}(\theta))}(\varphi_{\mathbf{n}}(1)).$$

More precisely, it is defined through replacing $f(0)$ by $\varphi_{\mathcal{A}}(\varphi_{\mathbf{n}}^{-1}(\theta))$ in the expression for the corresponding matrix element of the operator:

$$H(\theta)=\left(H(\mathbf{n}\cdot\alpha-0)-f(1-0)\langle e_{-\mathbf{n}},\cdot\rangle e_{-\mathbf{n}}\right)+\varphi_{\mathcal{A}}(\varphi_{\mathbf{n}}^{-1}(\theta))\langle e_{-\mathbf{n}},\cdot\rangle e_{-\mathbf{n}},\quad\theta\in\overline{\mathcal{G}_{\mathbf{n}}}.$$

The result of this construction is a family of operators $\{H(\theta)\colon\theta\in\mathcal{S}\}$, with the map $\theta\mapsto H(\theta)$ being a homeomorphism between $\mathcal{S}$ and the above set of operators in the $\overline{\mathbb{R}}$-strong topology. The action of $\mathbb{Z}^{d}$ naturally extends to $\mathcal{S}$ from the conjugation action on $H(\theta)$.

Remark 4.7.

The family $\{H(\theta)\colon\theta\in\mathcal{S}\}$ is an extension of the family $\{H(\theta)\colon\theta\in\mathcal{C}_{\alpha}\}$ which, in turn, is the closure of the family $\{H(x)\colon x\in[0,1)\}$. As mentioned above, the latter inclusion allows to uniquely associate some value $\theta\in\mathcal{C}_{\alpha}\subset\mathcal{S}$ to each value $x\in[0,1)$.

For $t\in\mathcal{A}$, denote by

$$\mathcal{S}_{\alpha,t}:=\{\varphi_{\mathbf{n}}(\varphi_{\mathcal{A}}^{-1}(t))\colon\mathbf{n}\in\mathbb{Z}^{d}\}$$

the set of parameters corresponding to operators with the same choice of $t$. It is easy to see that $\mathcal{S}_{\alpha,t}$ are invariant with respect to the translation, and

$$\cup_{t\in\mathcal{A}}\mathcal{S}_{\alpha,t}=\cup_{\mathbf{n}\in\mathbb{Z}^{d}}\overline{\mathcal{G}_{\mathbf{n}}},\quad\mathcal{S}\setminus\left(\cup_{t\in\mathrm{int}\,\mathcal{J}}\mathcal{S}_{\alpha,t}\right)=\mathcal{C}_{\alpha}.$$

In other words, $S_{\alpha,t}$ contains one point in each gap $\mathcal{G}_{\mathbf{n}}$. If $\alpha$ has at least one irrational component, these gaps are dense, and we have

$$\overline{\mathcal{S}_{\alpha,t}}=\mathcal{S}_{\alpha,t}\cup\mathcal{C}_{\alpha}.$$

Remark 4.8.

It is easy to see that, in the case of a simple discontinuity, all operators $H(\theta)$ are different for different $\theta\in\mathcal{S}$. If one repeats the steps of the proof for the “trivial” case of $f$ such as the associated circle map is continuous (as in Theorem 4.4), the homeomorphism claim in Lemma 4.5 will fail, since the operators at the endpoints of the gaps will actually be equal. The set (4.5) will be closed, and there will be no need to add new points.

4.5. Continuity of spectra

We will now discuss continuity in $\alpha$. Whenever the dependence on $\alpha$ is important, we will use the notation $H(\alpha,\theta)$.

Lemma 4.9.

Suppose that $\theta_{j}\to\theta\in\mathcal{S}$, and $\alpha_{j}\to\alpha$, where the components of $\alpha$ are rationally independent. Then $H(\alpha_{j},\theta_{j})\to H(\alpha,\theta)$ $\overline{\mathbb{R}}$-strongly. If $\theta_{j}\in\mathcal{S}_{\alpha_{j},t_{j}}$, then either $\theta\in\mathcal{C}_{\alpha}$ or $\theta\in\mathcal{S}_{\alpha,t}$ and $t_{j}\to t$. If $\theta\in\mathcal{C}_{\alpha}$, then one can replace $\theta_{j}$ by (any) closest to it point of $\mathcal{C}_{\alpha_{j}}$, thus making $\theta_{j}\in\mathcal{C}_{\alpha_{j}}$.

Proof.

Follows from the fact that, for each $\mathbf{n}\in\mathbb{Z}^{d}$, the endpoints of each gap $G_{\mathbf{n}}$ are continuous in $\alpha$ at all $\alpha$ with rationally independent components. ∎

Remark 4.10.

In the case when $\alpha$ has rationally dependent components, the conclusion of Lemma 4.9 may not be true, and the set of such limit points is more complicated. In this case, multiple gaps with different labels are combined into one gap. If $\mathbf{n}\cdot\alpha\in\mathbb{Z}$, then every gap whose label is a multiple of $\mathbf{n}$ becomes a part of the gap with label $\mathbf{0}$. If $\alpha_{j}$ have rationally independent components, one can choose $\theta_{j}$ between the gaps with labels $\mathbf{0}$ and $\mathbf{n}$. Due to compactness, the operators $H(\theta_{j})$ must have a limit point. It is easy to see that it does not correspond to any point of $\mathcal{S}$. For example, if $d=1$ and $\alpha_{j}\to\alpha=p/q$, the limit point of $H(\alpha_{j},\theta_{j})$ may not be a $q$-periodic operator, since the choices of the $\pm$ will not necessarily be $q$-periodic.

We will now summarize some properties of invariance of spectra.

1. For $\alpha$ with at least one irrational component, the spectrum of $H(\theta)$ is constant in $\theta\in\mathcal{C}_{\alpha}$. We will denote it by $\Sigma(\alpha)$.

2. For all $\alpha$, the spectrum $H(\theta)$ is constant in $\theta\in\mathcal{S}_{\alpha,t}$ (since the corresponding operators are obtained from one another by translation and therefore are unitarily equivalent). We will use the notation $\Sigma(\alpha,t):=\sigma(H(\theta))$ for $\theta\in\mathcal{S}_{\alpha,t}$. Since $H(\theta)$ is a rank one perturbation of $H(\theta^{\prime})$ with some $\theta^{\prime}\in\mathcal{C}_{\alpha}$ and the spectrum of the latter has no isolated points, we have that $\Sigma(\alpha,t)$ contains $\Sigma(\alpha)$ and at most one eigenvalue in each gap of $\Sigma(\alpha)$. The notation $\Sigma(\alpha,t)$ is used in order to reflect the fact that, in general, these eigenvalues will not be constant in $t$.

3. For $\alpha$ with rationally dependent components, the spectrum may not be constant in $\theta$. We will denote it by $\Sigma(\alpha,\theta)$.

With all the preparations, the proof of the following continuity result can now be performed in the same language as in [3].

Theorem 4.11.

Suppose that the components of $\alpha$ are rationally independent, and let $\alpha_{j}\to\alpha$ be a sequence of vectors with all rational components. Then the following holds:

1. $E\in\Sigma(\alpha)$ if and only if there exists $E_{j}\in\Sigma(\alpha_{j},\theta_{j})$ with $\theta_{j}\in\mathcal{C}_{\alpha_{j}}$, $E_{j}\to E$.

2. $E\in\Sigma(\alpha,t)$ if and only if there exists $E_{j}\in\Sigma(\alpha_{j},t_{j})$ with $t_{j}\to t$ and $E_{j}\to E$.

3. As a consequence, the union spectra $\cup_{\theta\in\mathcal{C}_{\alpha}}\bar{\sigma}(H(\alpha,\theta))$ and $\cup_{\theta\in\mathcal{S}}\bar{\sigma}(H(\alpha,\theta))$ are continuous in $\alpha$ at all rationally independent points, with respect to the Hausdorff distance topology.

Proof.

As in [3], the “only if” part in (1) and (2) is standard. In both cases, Lemma 4.9 allows to find a sequence of operators with rational $\alpha_{j}$ converging to $H(\alpha,\theta)$ $\overline{\mathbb{R}}$-strongly, and the implication follows from Corollary 3.3.

For the “if” part in (1), note that, since $f$ cannot be identically equal to $\infty$, for large enough $j$ one will be able to find a sequence of bounded solutions of the eigenvalue equation

(4.7) $$H(\alpha_{j},\theta_{j}^{\prime})\psi=E\psi,\quad\|\psi\|_{\ell^{\infty}(\mathbb{Z}^{d})}\leq 1,\quad\psi(0)\geq 1/2,$$

where $\theta_{j}^{\prime}\in\mathcal{C}_{\alpha_{j}}$ is a translation of $\theta_{j}$. By the diagonal process, one can assume that $\theta_{j}^{\prime}\to\theta\in\mathcal{C}_{\alpha}$ and $\psi_{j}\to\psi$ pointwise, with the same properties. From Proposition 3.4, this implies $E\in\bar{\sigma}(H(\theta))$.

A similar argument works for (2), where one can find a converging sequence $\theta_{j}^{\prime}\in\mathcal{S}_{\alpha_{j},t_{j}}$ and a sequence $\psi_{j}$ satisfying (4.7). As in Lemma 4.9, the limit of $\theta_{j}^{\prime}$ will be either from $\mathcal{S}_{\alpha,t}$ or from $\mathcal{C}_{\alpha}$, which in both cases implies $E\in\Sigma(\alpha,t)$.

The proof of (3) goes along the same lines as (1), (2), as well as in [3]. If $E\in\cup_{\theta\in\mathcal{C}_{\alpha}}\bar{\sigma}(H(\alpha,\theta))$, one can find a sequence of rational approximants to $\alpha$ and conclude, by a similar diagonal process, that there exists $\theta\in\mathcal{C}_{\alpha}$ and a bounded non-trivial solution of the eigenvalue equation $H(\theta)\psi=E\psi$. Then, continuity follows from another application of the same diagonal process. The same holds for $\mathcal{C}_{\alpha}$ replaced by $\mathcal{S}$. ∎

Remark 4.12.

For $\alpha$ with rationally dependent components, the argument in (3) above appears to fail by the same reason as in Remark 4.10.

4.6. The case of rational $\alpha$

For rational $\alpha$ (that is, $\alpha$ with all components being rational), the operator $H(\theta)$ becomes a periodic operator. Suppose that, additionally, $\alpha=(p_{1}/q_{1},\ldots,p_{d}/q_{d})$, with $q_{1},\ldots,q_{d}$ pairwise coprime. Let $q=q_{1}\ldots q_{d}$. There exists $\mathbf{n}_{q}\in\mathbb{Z}^{d}$ such that $\mathbf{n}_{q}\cdot\alpha=1/q\,\,\mathrm{mod}\,\,\mathbb{Z}$. The set $\mathcal{C}_{\alpha}$ becomes a union of $q$ arcs, with the set of gap labels being $0,1/q,2/q,\ldots,(q-1)/q$. The operators $H(\theta)$ and $H(\theta+1/q)$ are unitarily equivalent (here, $H(\theta+1/q)$ is defined using the translation by $\mathbf{n}_{q}$). In terms of variable $x$, we also have $H(x)$ unitarily equivalent to $H(x+1/q)$, as well as $H(x-0)$ to $H(x-0-1/q)$. It is also convenient to represent $\mathcal{C}_{\alpha}$ in terms of variable $x$ as a disjoint union

$$[0,1/q]\sqcup[1/q,2/q]\sqcup\ldots\sqcup[(q-1)/q,1],$$

where the point $1/q$ on the interval $[0,1/q]$ denotes $H(1/q-0)$, and the one on $[1/q,2/q]$ corresponds to $H(1/q+0)$.

Let $h(\theta)$ be $H(\theta)$ restricted to $\{0,1,2,\ldots,q_{1}-1\}\times\ldots\times\{0,1,2,\ldots,q_{d}-1\}$ with periodic boundary conditions (in other words, the operator with quasimomentum zero). It is well known that $\bar{\sigma}(h(\theta))\subset\bar{\sigma}(H(\theta))$. The following is the gap filling result for rational $\alpha$.

Theorem 4.13.

Suppose that $\alpha=(p_{1}/q_{1},\ldots,p_{d}/q_{d})$, with $q_{1},\ldots,q_{d}$ pairwise coprime. Define $h(\theta)$ as above, and suppose that the map $\tilde{f}$, defined in (4.3), is of non-zero topological degree. Then

$$\bigcup_{\theta\in\mathcal{S}}\bar{\sigma}(h(\theta))=\overline{\mathbb{R}}.$$

Before proceeding to the proof, we will need to use a well known fact about families of unitary operators. Let $u\colon\mathbb{T}\to U(N)$ be a continuous map from the circle to the group of $N\times N$ unitary matrices. The map $t\mapsto\det(u(t))$ is a map from $\mathbb{T}$ to the unit circle in $\mathbb{C}$.

Proposition 4.14.

Let $u\colon\mathbb{T}\to U(N)$ be a continuous map from the circle to the group of $N\times N$ unitary matrices, and suppose that the map $t\mapsto\det(u(t))$ has a non-zero winding number. Then

(4.8) $$\bigcup_{t\in\mathbb{T}}\sigma(u(t))=\mathbb{T}.$$

Proof.

It is easy to see that, if $\eqref{eq_union_unitary}$ does not hold, then there is an open arc of $\mathbb{T}$ such that $\sigma(u(t))$ avoids that arc, and the complement of that arc can be continuously contracted into a point. As a consequence, the map $u$ is homotopic to a constant, and the corresponding determinant map will have a zero winding number. ∎

Proof of Theorem 4.13. Consider the map $\theta\mapsto U_{h(\theta)}$. In view of the above, it is sufficient to show that the corresponding determinant map has a non-zero winding number. Since the said winding number is a topological invariant, one can replace $\Delta$ by $s\Delta$, where $s$ continuously changes from $1$ to $0$ and reduce it to the case of a diagonal operator. In that case, each diagonal entry can be considered as a scalar unitary operator, with the winding number of the determinant equal to that of the map $\tilde{f}$. As a consequence, each diagonal entry will produce an equal non-zero contribution to the total winding number, which completes the proof.

Conclusion of the proofs of Theorems 4.2 and 4.4. Theorem 4.2 can now be obtained as a consequence of Theorem 4.13 and part 3) of Theorem 4.11. In view of Remark 4.8, the same arguments can also be applies in the simpler case of Theorem 4.4, with the convention $\mathcal{C}_{\alpha}=\mathcal{S}_{\alpha}$ associated with the original set of values of $x\in[0,1)$.  ∎

Proof of Corollary 4.3. From the proof of Theorem 4.11, we see that, for every $E\in\mathbb{R}$, there exists $\theta\in\mathcal{S}$ and a bounded solution to the eigenvalue equation $H(\theta)\psi=E\psi$. If $\theta\in\mathcal{C}_{\alpha}$, this solution will decay exponentially from Proposition 2.6 applied to the original operator $H$ or the one with $f$ replaced by $f(\cdot-0)$ (see Remark 2.7). If $\theta$ is in some gap, one can also apply Proposition 2.6 with $f$ replaced by $f_{t}$.  ∎