COUPLED-MODE CALCULATION WITH THE R-MATRIX PROPAGATOR FOR THE DISPERSION OF SURFACE-WAVES ON A TRUNCATED PHOTONIC CRYSTAL

We describe a modal expansion method with an R-matrix propagation algo rithm that can be applied to calculate the dispersion of surface and g uided waves supported by a finite thickness photonic crystal. Easy to implement, the R matrix has the advantage of inherent numerical stabil ity and allows photonic crystals of complex structure, which are many waves in thickness to be analyzed. This method is also computationally simpler that the supercell method when looking for surface modes. It can also be applied to many other dispersion problems including multil ayer gratings and waveguides with internal structure. We apply this me thod to calculate dispersion of surface waves supported by a slablike photonic crystal structure, which consists of several rows of an infin itely periodic array of square dielectric cylinders. In the bulk, the periodic array has a photonic band gap and for the finite thickness st ructure, we calculate dispersion of surface waves that propagate along the boundaries at frequencies within the band gap. We vary the thickn ess of the cylinders in the outer layer and calculate the change in di spersion. When the overall thickness of the photonic slab is decreased , coupling and splitting between surface modes at each interface is ob served.