Kernel compensation method for Maxwell eigenproblem with mimetic finite difference discretization

We present a kernel compensation method for Maxwell eigenproblem for photonic crystals to avoid the infinite-dimensional kernels that cause many difficulties in the calculation of energy gaps. The quasi-periodic problem is first transformed into a periodic one on the cube by the Floquet-Bloch theory. Then the compensation operator is introduced in Maxwell's equation with the shifted curl operator. The discrete problem depends on the compatible discretization of the de Rham complex, which is implemented by the mimetic finite difference method in this paper. We prove that the compensation term exactly fills up the kernel of the original problem and avoids spurious eigenvalues. Also, we propose an efficient preconditioner and its FFT and multigrid solvers, which allow parallel computing. Numerical experiments for different three-dimensional lattices are performed to validate the accuracy and effectiveness of the method.