Localisation and Delocalisation for a Simple Quantum Wave Guide with Randomness

In this paper, we consider Schrödinger operators on M×Zd2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$M\times {\mathbb {Z}}^{d_{2}}$$\end{document}, with M={M1,…,M2}d1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$M=\{M_{1},\ldots ,M_{2}\}^{d_{1}}$$\end{document} (‘quantum wave guides’) with a ‘Γ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Gamma $$\end{document}-trimmed’ random potential, namely a potential which vanishes outside a subset Γ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Gamma $$\end{document} which is periodic with respect to a sub-lattice. We prove that (under appropriate assumptions) for strong disorder these operators have pure point spectrum outside the set Σ0=σ(H0,Γc)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Sigma _{0}=\sigma (H_{0,\Gamma ^{c}})$$\end{document} where H0,Γc\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H_{0,\Gamma ^{c}} $$\end{document} is the free (discrete) Laplacian on the complement Γc\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Gamma ^{c} $$\end{document} of Γ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Gamma $$\end{document}. We also prove that the operators have some absolutely continuous spectrum in an energy region E⊂Σ0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathcal {E}}\subset \Sigma _{0}$$\end{document}. Consequently, there is a mobility edge for such models. We also consider the case -M1=M2=∞\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$-M_{1}=M_{2}=\infty $$\end{document}, i.e. Γ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Gamma $$\end{document}-trimmed operators on Zd=Zd1×Zd2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbb {Z}}^{d}={\mathbb {Z}}^{d_{1}}\times {\mathbb {Z}}^{d_{2}}$$\end{document}. Again, we prove localisation outside Σ0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Sigma _{0} $$\end{document} by showing exponential decay of the Green function GE+iη(x,y)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$G_{E+i\eta }(x,y) $$\end{document} uniformly in η>0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta >0 $$\end{document}. For all energies E∈E\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E\in {\mathcal {E}}$$\end{document} we prove that the Green’s function GE+iη\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$G_{E+i\eta } $$\end{document} is not (uniformly) in ℓ1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\ell ^{1}$$\end{document} as η\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta $$\end{document} approaches 0. This implies that neither the fractional moment method nor multi-scale analysis can be applied here.


Introduction
Quantum waveguides are quantum mechanical structures which are confined in certain spaces dimensions, but unconfined in others. The last decades showed a growing interest in these systems in the mathematical literature. The book [8] by Exner and Kovařík gives an overview on the state of art (as of 2015) as well as an extensive list of references on waveguides.
A quantum waveguide with the simplest geometry is given by a particle in a (k-dimensional) strip in Z k+m or R k+m . Other examples are tubes or wires which are bended or twisted (see, for example, the discussion in Krejčiřík [21]). Of particular interest are waveguides with randomness either in the geometry of the system or in the potential energy.
In this paper, we consider waveguides with a simple geometry, namely on a strip in Z d , for example, on with a random potential V ω . The potential we consider is 'sporadic' or 'Γtrimmed', in the sense that V ω (x) = 0 for lattice points x ∈ Γ. Here, Γ is an L-periodic subset of the strip for a sub lattice L. In example (1), we may choose, for instance, For x ∈ Γ, the potentials are independent and identically distributed. Random operators H ω with such potentials are called 'Γ-trimmed'. Spectral theory for trimmed Anderson models (i.e. on Z d ) was done in the PhD-thesis of Obermeit [23] and the papers of Rojas-Molina [24], Elgart-Klein [5], Elgart-Sodin [6] and Kirsch-Krishna [16]. The present paper was inspired by [6].
We show that models as in (1), (2) have a mobility edge (or rather mobility edges). The measure theoretical nature of the spectrum depends on the energy region. More precisely, denote by H 0,Γ c the free (discrete) Laplacian on the set Γ c . Outside of the spectrum Σ 0 := σ(H 0,Γ c ) the operator H ω has dense point spectrum for high enough disorder. On the other hand, we prove that H ω has some absolutely continuous spectrum inside σ(H 0,Γ c ) regardless of the strength (or even existence) of the randomness.
The absolutely continuous spectrum comes from the existence of canonical extended states. More precisely, in an energy region inside Σ 0 we find periodic solutions of the free Schrödinger equation which vanish on the set Γ. Hence, these functions solve the Schrödinger equation with a random (in fact, with an arbitrary) potential on Γ as well.
Kotani and Simon consider random operators on a strip [20] with d 2 = 1. They give abstract conditions for absolutely continuous spectrum in terms of Lyapunov exponents. The examples for absolutely continuous spectrum they give are deterministic potentials. Their method does not extend to d 2 > 1.
To prove pure point spectrum we employ the multiscale analysis (see, for example, Dreifus-Klein [27] or Disertori et.al. [4] and references given there). The critical ingredient in our case is a Wegner Estimate. To prove this estimate, we use in an essential way that we work outside the spectrum σ(H 0,Γ c ). In fact, the estimate blows up when we approach σ(H 0,Γ c ).
From the multiscale bounds not only pure point spectrum follows but also dynamical localisation (see  or [26]).
We also consider the case −M 1 = M 2 = ∞, in other words a Γ-trimmed potential on Z d = Z d1 × Z d2 . Again we use multiscale analysis outside Σ 0 for high enough disorder, which implies uniform exponential decay of the Green function and existence of pure point spectrum.
Inside a subset E of Σ 0 we prove that the Green function G E+iε (x, y) is not only not uniformly exponentially decaying, but that for x / ∈ Γ even Under appropriate conditions, we even have E = Σ 0 . Consequently, we have a 'phase transition' at Σ 0 = σ(H 0,Γ c ) which manifests itself in the behaviour of the Green function.
For the Anderson model (with full randomness), it is expected that in higher dimension there is a mobility edge, namely a transition from pure point spectrum to absolutely continuous spectrum depending on the energy range and the strength of the disorder. All that is rigorously known (on Z d ) is the existence of pure point spectrum (see, for example, Aizenman-Warzel [2] or Kirsch [15]). However, on the Bethe tree Klein [17] proved the existence of absolutely continuous spectrum (see also , Aizenman-Warzel [1], Froese et al. [9]).
There are random Schrödinger operators with decaying randomness for which a mobility edge is known to exist (see Krishna [22], Kirsch et.al. [10,13], and Jaksic-Last [11]. These models are not ergodic. However, the models we consider here are either ergodic in Z d1 -direction (for the strip) or even ergodic with respect to a d-dimensional sub-lattice.
The systems we consider are about the simplest wave guides possible. We expect that the localisation results we have in this paper can be extended to more complicated wave guide systems using essentially the same technique. We also expect the delocalisation results to hold more generally, assuming there are big enough regions without (random) potential. However, such results would presumably require refined methods.

Setup
We consider quantum systems (wave guides) on By e ν we denote the standard basis of Z d . The lattice L and the subset L M1M2 are defined by and with M 1 , M 2 ∈ Z, M 1 < M 2 . In this paper, we always assume that d 1 , d 2 > 0. Then, we define 'cubes' X M1M2 by 'periodising' C 0 : and Informally, we consider X ∞ as Thus, X M1M2 is a strip ('waveguide') of restricted width in d 1 directions and unconfined in d 2 directions, X ∞ is (for the moment) just a complicated expression for Z d .
Sometimes we omit the indices M 1 , M 2 and ∞ if they are clear from the context or if they are irrelevant. For simplicity of arguments and, in fact, without loss of generality we assume that M 2 − M 1 is even.
The discrete Laplacian H 0 on X ∞ = Z d is given by: When we restrict H 0 to subsets of Z d we have to impose boundary conditions. In the following, we will most of the time work with either 'simple' boundary conditions or with 'periodic' boundary conditions.
is called the Laplacian on Λ with simple boundary conditions. Here Definition 2.2. Suppose the box Λ ⊂ Z d is given by for some d ≤ d, then we call the operator H Λ 0 defined by where the Laplacian on Λ with periodic boundary conditions.
For the operator H 0 on X M1M2 , we impose periodic boundary conditions. Now, we define the set Γ of 'active sites', i.e. the sites where the potential V ω may be nonzero. The active sites inside C 0 are denoted by Γ 0 , with Vol. 23 (2022) Localisation and Delocalisation 4135 and Γ := Γ 0 + L M1M2 (17) where we include again the case −M 1 = M 2 = ∞.
For some results, we can deal with the following less restrictive model In this article, we investigate spectral properties of operators H on X of the form: where the potential V is supported by Γ, i.e. V (n) = 0 for n / ∈ Γ. Most of the time we suppose that V is a random potential with independent, identically distributed random variables V ω (γ), γ ∈ Γ, but some of our results are independent of such an assumption.

Results
We prove localisation under fairly weak assumptions on Γ.
Theorem 3.1. Suppose that ∅ = Γ = X . Assume that the random variables V ω (n), n ∈ Γ are independent with a common distribution P 0 which has a bounded density ρ (with respect to Lebesgue measure) with compact support. If The multiscale analysis gives a form of localisation which is more than merely pure point spectrum, namely dynamical localisation. Dynamical localisation appears in various forms (for a detailed discussion see, for example, [18]), we use it in the following form which we take from [26]. We say that the random operator H satisfies dynamical localisation in the energy interval I if for any p > 0 for any compactly supported ϕ. That dynamical localisation follows from multiscale analysis is proved in [3], see also [26].
The following result shows that Theorem 3.1 is not an empty statement. and Theorem 3.1 and Corollary 3.4 are proved in Sect. 5. The proof of Proposition 3.5 is given in Sect. 4. These results reprove and extend previous results in [5,6,16,23,24]. Now, we turn to a class of examples for which we can prove the existence of absolutely continuous spectrum. Definition 3.6. We call Γ as in (16) and (17) In particular, if Γ is a strict single layer set, then Vol. 23 (2022) Localisation and Delocalisation 4137 Theorem 3.9 applies in particular to Γ-trimmed random potential as in Theorem 3.1. Thus, for such a random potential there is an energy region with pure point spectrum and a region with absolutely continuous spectrum. Consequently, there exists a mobility edge.
The proof of Theorem 3.9 is contained in Sect. 6. Unfortunately, the proof of Theorem 3.9 does not work for the case −M 1 = M 2 = ∞, i.e. for X = Z d . However, for this case we can at least show, that both the fractional moment method and the multiscale analysis cannot work. In fact, the spectral values in E belong to 'extended states' in an informal sense.
Let Γ ⊂ X ∞ be a one layer set and V ω be a random potential on Γ satisfying the assumptions of Theorem 3.1. Denote by G Vω E+iζ (x, y) the Green function (i.e. the kernel of the resolvent ( Then, we show Theorem 3.10. We assume X = Z d and Γ is a single layer set.
for all x / ∈ Γ and all ω.
Part 2 of Theorem 3.10 is actually a deterministic result, and it holds for any potential vanishing outside Γ.
We prove this theorem in Sect. 7.

The Random Operator and Its Spectrum
In this section, we consider operators with random potential. To emphasise this, we write Hypothesis 4.1. We suppose that the potentials V ω (γ), γ ∈ Γ are i.i.d. with a common distribution P 0 . We assume that the support S of P 0 is compact.
The assumption that S is compact can be relaxed considerably but we do not bother to do so.
We denote the corresponding probability space by (Ω, F, P) which can and will be taken to be S Γ , γ∈Γ B(S), γ∈Γ P 0 .
We write the lattice L as L = L 1 × L 2 with L 1 ⊂ Z d1 and L 2 ⊂ Z d2 and points x in X as x = (x 1 , x 2 ) with x i ∈ Z di . We define 'shift' operators T j , j ∈ L := L 2 on (Ω, F, P) by It is easy to see that the shift T j is measure preserving, i.e. P(T j −1 A) = P(A) for every A ∈ F.
The following result tells us that the family {T j } j∈L is ergodic: This result can be found in [7], for example. Define for j ∈ L the shift operator for (x 1 , x 2 ) ∈ X . The operators U j are unitary on 2 (X ); moreover, the operators H ω are ergodic in the sense  We now investigate the spectrum (as a set).

Definition 4.5.
A function W : X ∈ R is called an admissible potential (with respect to P 0 ) if We denote the set of admissible potentials by A.

Remark 4.6.
Taking there is a one-to-one correspondence τ between Ω and the set A of admissible potentials, namely τ (ω)(n) = γ∈Γ ω γ δ γ n . We may therefore identify Ω and A.
Set S n = supp ϕ n which is a finite set. By the Borel-Cantelli lemma, there is a vector j n ∈ L such that With ψ n (x) = ϕ n (x + j n ), we therefore get thus E ∈ σ(H ω ). 2. Since the set A has probability one (in the sense of Remark 4.6) This together with 1. proves the theorem.
Set V a = aχ Γ then V a is an admissible potential if a ∈ supp P 0 . We define H a = H 0 + V a and set E min (a) = inf σ(H a ) and E max (a) = sup σ(H a ). Suppose now, that W is an admissible potential then by monotonicity We will have a closer look at the operators H a . Let us denote by H 0,Γ and H 0,Γ c the operator H 0 restricted to 2 (Γ) and 2 (Γ c ), respectively, with simple boundary conditions.  Proof. Take ϕ ∈ 2 (Γ c ) with ϕ 2 (Γ c ) = 1 and define ϕ ∈ 2 (X ) by Then, From (47), the equalities (43) and (44) follow by the Min-Max-Principle.
To prove the strict inequality (45), we show that sup H a is strictly increasing with a.
The operator H a is periodic; therefore, E a = sup σ (H a ) is given by the sup h a where h a is the operator H a restricted to the periodic cell C 0 with periodic boundary conditions. The corresponding eigenfunction ψ a is strictly positive.
By the Hellmann-Feynman theorem This proves 45. The inequality in Remark 4.10 is proved in a similar way. Now, we consider the case that Γ = G × Z d2 . In this case, the operator H a separates in the following way.
Definition 4.11. The space X splits in a part Z ⊂ Z d1 and Z d2 , namely We denote the Laplacian on 2 (Z) (possibly with periodic boundary conditions) by H (1) 0 and the Laplacian on 2 (Z d2 ) by H (2) 0 . We also set H Then Consequently This proves the following Corollary: The above consideration also allows us to prove Proposition 3.8.

Localisation
In this section, we prove pure point spectrum with exponentially decaying eigenfunction for energies E / ∈ σ(H 0,Γ C ) for high enough disorder. For the (L-periodic) set Γ of 'active' sites, we merely assume that Γ = ∅. We may also assume that Γ = C 0 , since otherwise Γ is the whole space, a case which is known, of course.
In the following, we consider the case X = X ∞ , at the end of this section we comment on the case X M1M2 .
We will use multiscale analysis (see [15] and references given there). During the proof we need to consider boxes Λ which are unions of shifted C 0 . Recall that the unit cell C 0 is defined by: One of the crucial ingredients of (most versions of) multiscale analysis is the Wegner-Estimate. To prove this, we need the following result. By H Λ ω , we denote the operator H ω = H 0 + V ω restricted to 2 (Λ) with periodic boundary conditions. and with a constant C which is independent of Λ and E.
Since we consider periodic boundary conditions on Λ, we have as any eigenfunction on Λ can be periodically extended to an eigenfunction on X .
Now we turn to the Wegner estimate. By N (A, E) we denote the number of eigenvalues of the operator A up to E.
where |Λ| denote the volume (number of points) of Λ.
Proof. The proof is a combination of the proofs from [15] and [12]. We sketch the main ideas. We use the abbreviation H = H ω (Λ). Let g be a monotone C ∞ -function with 0 ≤ g(t) ≤ 1, g(t) = 0 for t ≤ −2ε and g(t) = 1 for t ≥ 2ε.
We obtain Let E n denote the eigenvalues of H ω (Λ) labelled in increasing order. These eigenvalues depend on the values v j : Thus, we may consider where ψ n is a normalised eigenfunction of H with eigenvalue E n . Above we used the Hellmann-Feynman theorem and, in the final step, Proposition 5.2.
Vol. 23 (2022) Localisation and Delocalisation 4143 Summarising we proved Suppose supp P 0 ⊂ [a, b] and denote by H(v j = c) the operator H with V j replaced by the value c, then We used that changing the potential at one site j is a rank one perturbation and 0 ≤ g(λ) ≤ 1.
Performing the integrals over the v k , k = j gives the desired result.
Once we have the Wegner estimate, the multiscale analysis follows the usual path. We need an initial length scale estimate and the induction step over growing length scales.
The initial length scale estimate follows directly from the Wegner estimate Theorem 5.3. In fact, as long as we are away from the spectrum of H 0,Γ c , we can make the right hand side of (65) as small as we like by taking ρ ∞ small. This corresponds to high disorder. Using a Combes-Thomas estimate, this allows us to prove the initial scale estimate. For details see Section 11.1 in [15].
The induction step follows the lines in [15] sections 9 and 10. The only difference being that we deal with periodic boundary conditions while [15] uses simple boundary conditions. If X = X M1M2 , we start the induction with a cube of the form:

Absolutely Continuous Spectrum
In this section we consider a special case of the above operators. We start with the following observation: were defined in Definition 4.11.
We discuss a class of examples for which Proposition 6.1 applies. We look at X or at the strip X M1M2 . In the latter case, we impose periodic boundary conditions and take M 2 − M 1 is even. This way we avoid to discuss various cases separately. Proof. Again by applying addition theorems and the fact that Φ L vanishes on G 0 , hence on G, we see that (83) holds. Moreover, since M 2 − M 1 is even Φ L satisfies periodic boundary conditions. Now, we are ready to prove Theorem 3.9.
Proof. (Theorem 3.9) Take E ∈ E, then E = e L + η for some η ∈ [−2d 2 , 2d 2 ]. Denote by E L the (one-dimensional) subspace of 2 (Z) generated by the eigenfunction Φ L . The (closed) subspace h L = E L ⊗ 2 (Z d2 ) of 2 (X ) is invariant under the operator H 0 + W and restricted to h L the operators H 0 + W and H 0 agree. Consequently, H 0 + W on h L is unitarily equivalent to H