Structure-Preserving Methods for Computing Complex Band Structures of Three Dimensional Photonic Crystals

Tsung-Ming Huang Department of Mathematics, National Taiwan Normal University Tiexiang Li School of Mathematics, Southeast University Jia-Wei Lin Department of Applied Mathematics, National Chiao Tung University Wen-Wei Lin Department of Applied Mathematics, National Chiao Tung University Guohua Liu School of Mathematics, Southeast University Tian Heng Department of Applied Mathematics, National Chiao Tung University

Numerical Analysis and Scientific Computing Numerical Linear Algebra mathscidoc:1908.25001

This work is devoted to the numerical computation of complex band structure $\k=\k(\omega)\in\mathbb C^3$ for positive $\omega$ of three dimensional isotropic dispersive or non-dispersive photonic crystals from the perspective of structured quadratic eigenvalue problems (QEPs). Our basic strategy is to fix two degrees of freedom in $\k\in\mathbb C^3$ and to view the remaining one as the eigenvalue of a quadratic operator pencil derived from Maxwell's equations. Then Yee's scheme is employed to discretize $\nabla\times$ and $\k\times$ operators in this quadratic operator pencil. Distinct from the others' works which either ignore or directly exploit the Hamiltonian structure of the spectrum of the resulting QEP, we reformulate this QEP into an equivalent $\top$-palindromic QEP to facilitate the use of superior structure-preserving algorithms. Ultimately we rely on the structured Arnoldi algorithm, namely the G$\top$SHIRA algorithm, to compute eigenvalues of a $\top$-skew-Hamiltonian pair which are near or in $[-2,2]$, a much narrower region than the whole positive real axis in the origin problem. Moreover, to accelerate the inner iterations of the G$\top$SHIRA algorithm, we propose the preconditioning technique, making most of the eigenmatrix, which can essentially be seen as the Kronecker product of three discrete Fourier transformation matrices, of the commutative discretized $\partial_x,\partial_y,\partial_z$ operators. The advantage of our method is discussed in detail and corroborated by several numerical results.
dispersive permittivity, complex band structure, gyroscopic quadratic eigenvalue problem, $\top$-palindromic quadratic eigenvalue problem, G$\top$SHIRA, FFT
[ Download ] [ 2019-08-19 11:07:07 uploaded by tianheng ] [ 764 downloads ] [ 0 comments ]
@inproceedings{tsung-mingstructure-preserving,
  title={Structure-Preserving Methods for Computing Complex Band Structures of Three Dimensional Photonic Crystals},
  author={Tsung-Ming Huang, Tiexiang Li, Jia-Wei Lin, Wen-Wei Lin, Guohua Liu, and Tian Heng},
  url={http://archive.ymsc.tsinghua.edu.cn/pacm_paperurl/20190819110708049652412},
}
Tsung-Ming Huang, Tiexiang Li, Jia-Wei Lin, Wen-Wei Lin, Guohua Liu, and Tian Heng. Structure-Preserving Methods for Computing Complex Band Structures of Three Dimensional Photonic Crystals. http://archive.ymsc.tsinghua.edu.cn/pacm_paperurl/20190819110708049652412.
Please log in for comment!
 
 
Contact us: office-iccm@tsinghua.edu.cn | Copyright Reserved