An efficient and unified method for band structure calculations of 2D anisotropic photonic-crystal fibers

In this article, band structure calculations of two dimensional (2D) anisotropic photonic-crystal fibers (PhCFs) are considered. In 2D PhCFs, Maxwell's equations for the transversal electric and magnetic mode become decoupled, but the difficulty, arising from the anisotropic permittivity epsilon\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\varvec{\varepsilon }}}$$\end{document} and/or permeability mu,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\varvec{\mu }}},$$\end{document} plaguing the frequency-domain finite difference method, especially the original Yee's scheme, is our top concern. To resolve this difficulty, we re-establish the connection between the lowest order finite element method with the quasi-periodic condition and Yee's scheme using 2D non-orthogonal mesh, whereby the decoupled Maxwell's equations in 2D anisotropic PhCFs are readily discretized into a generalized eigenvalue problem (GEP). Moreover, we spell out the nullspace of the resulting GEP, if it exists, and explicitly construct the Moore-Penrose pseudoinverse of the singular coefficient matrix, whose smallest positive eigenvalues can be solved by the inverse Lanczos method. Extensive band structures of 2D PhCFs are calculated and benchmarked against reliable results to demonstrate the accuracy and efficiency of our method.