Spectral Analysis of the Schrödinger Operator for the Incommensurate System

Many novel and unique physical phenomena of incommensurate systems can be illustrated and predicted using the spectra of the associated Schrödinger operators. However, the absence of periodicity in these systems poses significant challenges for obtaining the spectral information. In this paper, by embedding the system into higher dimensions together with introducing a regularization technique, we prove that the spectrum of the Schrödinger operator for the incommensurate system can be approximated by the spectra of a family of regularized Schrödinger operators, which are elliptic, retain periodicity, and enjoy favorable analytic and spectral properties. We also show the existence of Bloch-type solutions to the Schrödinger equation for the incommensurate system, which can be well approximated by the Bloch solutions to the equations associated with the regularized operators. Our analysis provides a theoretical support for understanding and computing the incommensurate systems.

💡 Research Summary

The paper addresses the longstanding challenge of analyzing the spectral properties of Schrödinger operators that describe incommensurate (quasi‑periodic) layered materials, such as twisted bilayer graphene and other Moiré systems. In these systems two periodic potentials V₁(r) and V₂(r) are defined on lattices R₁ and R₂ that are not mutually commensurate, so the combined potential lacks any global translational symmetry. Consequently, the usual Bloch theorem does not apply and standard techniques for periodic operators cannot be used to obtain the spectrum or eigenfunctions.

To overcome this difficulty the authors embed the original d‑dimensional problem into a higher‑dimensional space ℝ^{2d} by treating the two layers as independent coordinates (r, r′). In this enlarged space they define an extended Schrödinger operator

  \tilde H = -½ ∑{i=1}^d (∂{r_i}+∂_{r′_i})² + V₁(r) + V₂(r′).

Because the potential now depends separately on r and r′, \tilde H is periodic with respect to the product lattice \tilde R = R₁ × R₂. However, the principal symbol of \tilde H vanishes on the subspace ξ+ξ′ = 0, making the operator degenerate elliptic. As a result, \tilde H defined on the natural Sobolev domain H²(ℝ^{2d}) is not closed and its spectrum would be the whole complex plane.

The first major contribution is a rigorous functional‑analytic treatment that restores self‑adjointness. By examining the symmetric part \tilde D = ∑ (∂{r_i}+∂{r′i})² and using Fourier analysis, the authors show that \tilde D is essentially self‑adjoint. Since the potentials are bounded, the Kato–Rellich theorem then guarantees that \tilde H is essentially self‑adjoint on H²(ℝ^{2d}). Its unique self‑adjoint extension is taken as the closure, equipped with the graph norm ‖u‖G = (‖u‖²{L²}+‖\tilde H u‖²{L²})^{½}.

Next, the Bloch–Floquet transform U_{\tilde R} is applied. This unitary map decomposes L²(ℝ^{2d}) into a direct integral over the Brillouin zone \tilde Γ* of fiber spaces L²( \tilde Γ). In each fiber the operator becomes a family of “k‑dependent” operators

  \tilde H(k) = -½ (∇r+∇{r′}+i(k+k′))² + ½|k+k′|² + V₁(r)+V₂(r′),

acting on periodic functions on the torus T^{2d}. The decomposition shows that the spectrum of \tilde H is the union over k of the spectra of these fiber operators.

Because \tilde H(k) still suffers from degeneracy, the authors introduce a regularization parameter δ>0 and define

  \tilde H_δ = \tilde H - (δ²/2) ∑{i=1}^d (∂{r_i} - ∂_{r′_i})².

For any positive δ the operator \tilde H_δ is uniformly elliptic, self‑adjoint, and has a purely discrete spectrum that varies continuously with δ. Using Fredholm theory and perturbation arguments, they prove that as δ → 0⁺ the spectra σ( \tilde H_δ ) converge (in the sense of Hausdorff distance) to σ( \tilde H ), i.e., the regularized family approximates the original degenerate operator.

Having established spectral approximation, the paper turns to eigenfunctions. It proves the existence of Bloch‑type generalized eigenfunctions for the original incommensurate operator H = -½Δ + V₁ + V₂, and shows that for each eigenvalue λ there exists a sequence of Bloch solutions ũ_{δ,λ} of the regularized problem \tilde H_δ ũ = λ ũ that converges (in L² and stronger Sobolev norms) to a solution of H u = λ u as δ → 0⁺. This provides a rigorous justification for using Bloch‑type ansätze in numerical simulations of incommensurate systems.

The paper’s structure is as follows: Section 2 introduces the incommensurate model, the higher‑dimensional embedding, and the Bloch‑Floquet machinery. Section 3 proves essential self‑adjointness of \tilde H and carries out the direct‑integral decomposition. Section 4 presents the regularization, establishes uniform ellipticity, and proves spectral convergence. Section 5 links the spectra of \tilde H and the original operator H, and demonstrates the existence and approximation of Bloch‑type solutions. Section 6 discusses implications, noting that the framework eliminates the need for large supercell constructions, provides a solid mathematical foundation for plane‑wave methods, and can be extended to arbitrary dimensions and multiple layers.

In summary, the authors deliver a comprehensive analytical framework that (i) restores periodicity by embedding into a product space, (ii) resolves the degeneracy through a controlled regularization, (iii) proves that the regularized spectra converge to the true spectrum of the incommensurate Schrödinger operator, and (iv) guarantees the existence of Bloch‑type solutions that can be approximated by those of the regularized, fully elliptic operators. This work bridges a critical gap between rigorous spectral theory and practical computational approaches for a broad class of quasi‑periodic materials, offering a mathematically sound alternative to supercell approximations and paving the way for accurate, efficient simulations of phenomena such as flat‑band physics, correlated insulating states, and unconventional superconductivity in twisted multilayer structures.