A finite element based fast eigensolver for three dimensional anisotropic photonic crystals

Abstract The standard Yee's scheme for the Maxwell eigenvalue problem places the discrete electric field variable at the midpoints of the edges of the grid cells. It performs well when the permittivity is a scalar field. However, when the permittivity is a Hermitian full tensor field it would generate un-physical complex eigenvalues or frequencies. In this paper, we propose a finite element method which can be interpreted as a modified Yee's scheme to overcome this difficulty. This interpretation enables us to create a fast FFT eigensolver that can compute very effectively the band structure of the anisotropic photonic crystal with SC and FCC lattices. Furthermore, we overcome the usual large null space associated with the Maxwell eigenvalue problem by deriving a null-space free discrete eigenvalue problem which involves a crucial Hermitian positive definite linear system to be solved in each of the iteration steps. It is demonstrated that the CG method without preconditioning converges in 37 iterations even when the dimension of a matrix is as large as 5 , 184 , 000 .

[1]  Markus Richter.  Optimization of Photonic Band Structures , 2010 .

[2]  Qing Huo Liu,et al.  Mixed Finite-Element Method for Resonant Cavity Problem With Complex Geometric Topology and Anisotropic Lossless Media , 2016, IEEE Transactions on Magnetics.

[3]  Ekmel Ozbay,et al.  2D Anisotropic Photonic Crystals of Hollow Semiconductor Nanorod with Liquid Crystals , 2013 .

[4]  Ben-Yuan Gu,et al.  Creation of partial band gaps in anisotropic photonic-band-gap structures , 1998 .

[5]  Daniele Boffi,et al.  Finite element approximation of eigenvalue problems , 2010, Acta Numerica.

[6]  K. Yee.  Numerical solution of initial boundary value problems involving maxwell's equations in isotropic media , 1966 .

[7]  Wen-Wei Lin,et al.  Singular Value Decompositions for Single-Curl Operators in Three-Dimensional Maxwell's Equations for Complex Media , 2014, SIAM J. Matrix Anal. Appl..

[8]  So-Hsiang Chou,et al.  Eigendecompositions and fast eigensolvers for Maxwell equations , 2016 .

[9]  Peter Monk,et al.  Finite Element Methods for Maxwell's Equations , 2003 .

[10]  Kunimasa Saitoh,et al.  Full-vectorial imaginary-distance beam propagation method based on a finite element scheme: application to photonic crystal fibers , 2002 .

[11]  Steven G. Johnson,et al.  Symmetry-protected topological photonic crystal in three dimensions , 2016, Nature Physics.

[12]  Q. Liu,et al.  The Mixed Finite-Element Method With Mass Lumping for Computing Optical Waveguide Modes , 2016, IEEE Journal of Selected Topics in Quantum Electronics.

[13]  Daniele Boffi,et al.  Modified edge finite elements for photonic crystals , 2006, Numerische Mathematik.

[14]  Zabel,et al.  Photonic band structures of optically anisotropic periodic arrays. , 1993, Physical review. B, Condensed matter.

[15]  Gang Bao,et al.  An adaptive edge element method with perfectly matched absorbing layers for wave scattering by biperiodic structures , 2010, Math. Comput..

[16]  D. Dobson.  An Efficient Method for Band Structure Calculations in 2D Photonic Crystals , 1999 .

[17]  Wen-Wei Lin,et al.  Electromagnetic field behavior of 3D Maxwell's equations for chiral media , 2019, J. Comput. Phys..

[18]  H. Whitney.  Geometric Integration Theory , 1957 .

[19]  Vivette Girault,et al.  Finite Element Methods for Navier-Stokes Equations - Theory and Algorithms , 1986, Springer Series in Computational Mathematics.

[20]  VOLKER MEHRMANN,et al.  A Newton-Type Method with Nonequivalence Deflation for Nonlinear Eigenvalue Problems Arising in Photonic Crystal Modeling , 2016, SIAM J. Sci. Comput..

[21]  J. Nédélec.  Mixed finite elements in ℝ3 , 1980 .

[22]  Liang Fu,et al.  Weyl points and line nodes in gyroid photonic crystals , 2012, Nature Photonics.

[23]  M. Shashkov,et al.  Mimetic Discretizations for Maxwell's Equations , 1999 .

[24]  Anders Bjarklev,et al.  Fabrication of Photonic Crystal Fibres , 2003 .

[25]  Dominique Pagnoux,et al.  Complete Analysis of the Characteristics of Propagation into Photonic Crystal Fibers, by the Finite Element Method , 2000 .

[26]  Roy A. Nicolaides,et al.  Convergence analysis of a covolume scheme for Maxwell's equations in three dimensions , 1998, Math. Comput..

[27]  Saulius Juodkazis,et al.  Microcavities in polymeric photonic crystals , 2001 .

[28]  Guo-Zhen Yang,et al.  LARGE ABSOLUTE BAND GAP IN 2D ANISOTROPIC PHOTONIC CRYSTALS , 1998 .

[29]  Weichung Wang,et al.  Eigendecomposition of the Discrete Double-Curl Operator with Application to Fast Eigensolver for Three-Dimensional Photonic Crystals , 2013, SIAM J. Matrix Anal. Appl..

[30]  Kurt Busch,et al.  Liquid-Crystal Photonic-Band-Gap Materials: The Tunable Electromagnetic Vacuum , 1999 .

[31]  Willy Dörfler,et al.  Photonic Crystals: Mathematical Analysis and Numerical Approximation , 2011 .

[32]  Q. Liu,et al.  Mixed Spectral-Element Method for 3-D Maxwell's Eigenvalue Problem , 2015 .

[33]  Wen-Wei Lin,et al.  Matrix representation of the double-curl operator for simulating three dimensional photonic crystals , 2013, Math. Comput. Model..

[34]  Robert Rogowski,et al.  Photonic band gap analysis using finite-difference frequency-domain method. , 2004, Optics express.

[35]  Wen-Wei Lin,et al.  Eigenvalue solvers for three dimensional photonic crystals with face-centered cubic lattice , 2014, J. Comput. Appl. Math..

[36]  J. Nédélec.  A new family of mixed finite elements in ℝ3 , 1986 .

[37]  Q. Liu,et al.  Spectral element method for band structures of two-dimensional anisotropic photonic crystals. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

[38]  R. Rumpf,et al.  Finite-difference frequency-domain algorithm for modeling electromagnetic scattering from general anisotropic objects , 2014 .

[39]  Qing Huo Liu,et al.  Accurate determination of band structures of two-dimensional dispersive anisotropic photonic crystals by the spectral element method. , 2009, Journal of the Optical Society of America. A, Optics, image science, and vision.

[40]  Partha P. Banerjee,et al.  Computational Methods for Electromagnetic and Optical Systems , 2000 .