Compact 2D FDTD method for modeling photonic crystal fibers

Combining perfectly matched layer (PML) for the boundary treatment, we present an efficient compact 2-dimensional finite-difference time-domain (2D FDTD) method for modeling photonic crystal fibers. For photonic crystal fibers, if we assume that the propagation constant along the propagation direction is fixed, three-dimensional hybrid guided modes can be calculated by using only a two-dimension mesh. Because of using the real variable method, the computation time, i.e., it is of order N. Comparing with the plane wave expansion method, FDTD make the computation time and computer memory are significantly reduced. The numerical results for a triangular lattice photonic crystal fiber are in very good agreement with the results from the local basis function method. This method can easily be used for any complicated inclusions.