In this paper we propose and analyze a hp-adaptive discontinuous finite element \ud method for computing the band structure of 2D periodic photonic crystals. The \ud problem can be reduced to the computation of the discrete spectrum of each member in a family of periodic Hermitian eigenvalue problems on the primitive cell, \ud parametrised by a two-dimensional parameter - the quasimomentum. We propose \ud a residual-based error estimator and show that it is reliable and efficient for all \ud eigenvalue problems in the family. In particular we prove that if the error estimator converges to zero then the distance of the computed eigenfunction from the \ud true eigenspace also converges to zero and the computed eigenvalue converges to a \ud true eigenvalue. The results hold for eigenvalues of any multiplicity. We illustrate \ud the benefits of the resulting hp-adaptive method in practice, both for fully periodic \ud crystals and also for the computation of eigenvalues in the band gaps of crystals \ud with defects
To submit an update or takedown request for this paper, please submit an Update/Correction/Removal Request.