An efficient contour integral based eigensolver for 3D dispersive photonic crystal

Tsung Ming Huang, Weichien Liao, Wen Wei Lin, Weichung Wang*

*Corresponding author for this work

Research output: Contribution to journalArticlepeer-review

1 Citation (Scopus)


Numerical simulations play a significant role in studying the properties of dispersive metallic photonic crystals. The dispersive photonic crystals are modeled by the Maxwell equations, and the equations are then discretized by the widely-used Yee's scheme. After applying certain similarity transformations to the discretized system, the original simulation problem becomes a non-Hermitian eigenvalue problem with clustered eigenvalues. An efficient contour integral (CI) based eigensolver is developed to overcome the difficulties of applying existing methods to solve eigenvalues in designated regions. This efficient method combines the contour integral, the fast matrix–vector multiplication, and efficient linear system solver. The numerical results illustrate the efficiency of our algorithm.

Original languageEnglish
Article number113581
JournalJournal of Computational and Applied Mathematics
Publication statusPublished - 2021 Oct 15


  • Contour integral based eigensolver
  • Discrete single-curl operator
  • Dispersive photonic crystal
  • Fast matrix–vector multiplication
  • The Maxwell equations

ASJC Scopus subject areas

  • Computational Mathematics
  • Applied Mathematics


Dive into the research topics of 'An efficient contour integral based eigensolver for 3D dispersive photonic crystal'. Together they form a unique fingerprint.

Cite this