Huang T.-MLiao WLin W.-WWEICHUNG WANG2021-07-212021-07-2120213770427https://www.scopus.com/inward/record.uri?eid=2-s2.0-85105690843&doi=10.1016%2fj.cam.2021.113581&partnerID=40&md5=9b95d6ab2a950ce051d1b7311f5c694bhttps://scholars.lib.ntu.edu.tw/handle/123456789/572154Numerical 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. ? 2021 Elsevier B.V.Eigenvalues and eigenfunctions; Linear systems; Matrix algebra; Photonic crystals; Contour integral based eigensolv; Contour integrals; Discrete single-curl operator; Dispersive photonic crystal; Eigen-value; Eigensolvers; Fast matrix–vector multiplication; Metallic photonic crystals; Property; The maxwell equation; Maxwell equationsjournal article10.1016/j.cam.2021.1135812-s2.0-85105690843