Newton-type method with nonequivalence deflation for nonlinear eigenvalue problems arising in photonic crystal modeling

Tsung Ming Huang, Wen Wei Lin, Volker Mehrmann

Research output: Contribution to journalArticle

4 Citations (Scopus)

Abstract

The numerical simulation of the band structure of three-dimensional dispersive metallic photonic crystals with face-centered cubic lattices leads to large-scale nonlinear eigenvalue problems, which are very challenging due to a high-dimensional subspace associated with the eigenvalue zero and the fact that the desired eigenvalues (with smallest real part) cluster and are close to the zero eigenvalues. For the solution of the nonlinear eigenvalue problem, a Newton-type iterative method is proposed and the nullspace-free method is applied to exclude the zero eigenvalues from the associated generalized eigenvalue problem. To find the successive eigenvalue/eigenvector pairs, we propose a new nonequivalence deflation method to transform converged eigenvalues to infinity, while all other eigenvalues remain unchanged. The deflated problem is then solved by the same Newtontype method, which is used as a hybrid method that combines with the Jacobi-Davidson and the nonlinear Arnoldi methods to compute the clustering eigenvalues. Numerical results illustrate that the proposed method is robust even for the case of computing many clustering eigenvalues in very large problems.

Original languageEnglish
Pages (from-to)B191-B218
JournalSIAM Journal on Scientific Computing
Volume38
Issue number2
DOIs
Publication statusPublished - 2016 Jan 1

    Fingerprint

Keywords

  • Dispersive metallic photonic crystals
  • Jacobi-Davidson method
  • Maxwell equation
  • Newton-type method
  • Nonequivalence deflation
  • Nonlinear Arnoldi method
  • Nonlinear eigenvalue problem
  • Shift-invert residual Arnoldi method

ASJC Scopus subject areas

  • Computational Mathematics
  • Applied Mathematics

Cite this