We study how to efficiently solve the eigenvalue problems in computing band structure of three-dimensional dispersive metallic photonic crystals with face-centered cubic lattices based on the lossless Drude model. The discretized Maxwell equations result in large-scale standard eigenvalue problems whose spectrum contains many zero and cluster eigenvalues, both prevent existed eigenvalue solver from being efficient. To tackle this computational difficulties, we propose a hybrid Jacobi–Davidson method (hHybrid) that integrates harmonic Rayleigh–Ritz extraction, a new and hybrid way to compute the correction vectors, and a FFT-based preconditioner. Intensive numerical experiments show that the hHybrid outperforms existed eigenvalue solvers in terms of timing and convergence behaviors.
- Clustered eigenvalues
- Hybrid Jacobi–Davidson method
- Three-dimensional dispersive metallic photonic crystals
- Zero eigenvalues
ASJC Scopus subject areas
- Hardware and Architecture
- Physics and Astronomy(all)