Adaptive curvilinear coordinates in a plane-wave solution of Maxwell’s equations in photonic crystals

A method is described to compute the modes propagating at a given frequency in dielectric systems that are periodic in two dimensions and uniform in the third dimension, using a plane-wave basis expressed in a system of generalized curvilinear coordinates. The coordinates are adapted to the structure under consideration by increasing the effective plane-wave cutoff in the vicinity of the interfaces between dielectrics, where the electromagnetic fields vary most rapidly. The favorable efficiency and convergence properties of the method are shown by comparison with the conventional plane-wave formulation of Maxwell's equations. Although the method is developed to study propagation in photonic crystal fibers, it is also applicable more generally to plane-wave modal solutions of structured dielectrics.

[1]  M. Skorobogatiy.  Modeling the impact of imperfections in high-index-contrast photonic waveguides. , 2004, Physical review. E, Statistical, nonlinear, and soft matter physics.

[2]  F Benabid,et al.  Ultrahigh efficiency laser wavelength conversion in a gas-filled hollow core photonic crystal fiber by pure stimulated rotational Raman scattering in molecular hydrogen. , 2004, Physical review letters.

[3]  P. Russell,et al.  Scaling laws and vector effects in bandgap-guiding fibres. , 2004, Optics express.

[4]  P. Roberts,et al.  Robust photonic band gaps for hollow core guidance in PCF made from high index glass. , 2003, Optics express.

[5]  P. Russell,et al.  Photonic Crystal Fibers , 2003, Science.

[6]  Steven G. Johnson,et al.  Geometric variations in high index-contrast waveguides, coupled mode theory in curvilinear coordinates. , 2002, Optics express.

[7]  Fetah Benabid,et al.  Particle levitation and guidance in hollow-core photonic crystal fiber. , 2002, Optics express.

[8]  D. Sorensen.  Numerical methods for large eigenvalue problems , 2002, Acta Numerica.

[9]  Steven G. Johnson,et al.  Block-iterative frequency-domain methods for Maxwell's equations in a planewave basis. , 2001, Optics express.

[10]  Knight,et al.  Single-Mode Photonic Band Gap Guidance of Light in Air. , 1999, Science.

[11]  H. Dr.  GENERALIZED-GRADIENT FUNCTIONALS IN ADAPTIVE CURVILINEAR COORDINATES , 1996 .

[12]  Gérard Granet,et al.  Differential covariant formalism for solving Maxwell's equations in curvilinear coordinates: oblique scattering from lossy periodic surfaces , 1995 .

[13]  Nicorovici,et al.  Lattice sums for off-axis electromagnetic scattering by gratings. , 1994, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[14]  Weizhang Huang,et al.  Moving mesh partial differential equations (MMPDES) based on the equidistribution principle , 1994 .

[15]  Ross C. McPhedran,et al.  An electromagnetic theory of dielectric waveguides with multiple embedded cylinders , 1994 .

[16]  Gygi.  Electronic-structure calculations in adaptive coordinates. , 1993, Physical review. B, Condensed matter.

[17]  J. Joannopoulos,et al.  Accurate theoretical analysis of photonic band-gap materials. , 1993, Physical review. B, Condensed matter.

[18]  T. Arias,et al.  Iterative minimization techniques for ab initio total energy calculations: molecular dynamics and co , 1992 .

[19]  F. Gygi,et al.  Adaptive Riemannian metric for plane-wave electronic-structure calculations , 1992 .

[20]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[21]  R. Holland.  Finite-Difference Solution of Maxwell's Equations in Generalized Nonorthogonal Coordinates , 1983, IEEE Transactions on Nuclear Science.

[22]  Robert D. Russell,et al.  Variational Mesh Adaptation Methods for Axisymmetrical Problems , 2003, SIAM J. Numer. Anal..