Quasiperiodic Elliptic Operators: Projection Method and Convergence Analysis
This work provides a rigorous numerical framework for quasiperiodic eigenvalue problems, addressing a key challenge in modeling photonic quasicrystals and other quasiperiodic systems.
The authors developed a high-accuracy projection method for computing eigenpairs of quasiperiodic elliptic operators, achieving spectral accuracy with rigorous convergence guarantees. Numerical experiments on 1D photonic quasicrystals and 2D/3D QEOs validated the method's accuracy and efficiency.
Quasiperiodic elliptic operators (QEOs) serve as fundamental models in both mathematics and physics, as exemplified by their role in the numerical modeling of one-dimensional photonic quasicrystals. However, distinct from periodic elliptic operators, approximating eigenpairs for QEOs poses significant challenges, particularly in capturing the full spectral structure (notably the continuous spectrum) and deriving convergence guarantees in the absence of compactness. In this paper, we develop a high-accuracy numerical method to compute eigenpairs of QEOs based on the projection method, which embeds quasiperiodic operators into a higher-dimensional periodic torus. To address the non-compactness issue, we construct a directional-derivative Hilbert space along irrational manifolds of a high-dimensional torus and characterize operators equivalent to QEOs within this space. By integrating spectral theory for non-compact operators into the Babuška-Osborn eigenproblem framework, we establish rigorous convergence analysis and prove that our method achieves spectral accuracy. Numerical experiments validate the accuracy and efficiency of the proposed method, including a one-dimensional photonic quasicrystal and two- and three-dimensional QEOs.