Abstract

The dispersion relation of Bloch waves is derived from the properties of a single grating layer. A straightforward way to impose the Bloch condition leads to the calculation of the eigenvalues of the transfer matrix through the single grating layer. Unfortunately, the transfer-matrix algorithm is known to be unstable as a result of the growing evanescent waves. This problem appears again in the calculation of the eigenvalues, making unusable the transfer matrix in numerous practical problems. We propose two different algorithms to circumvent this problem. The first one takes advantage of scattering matrices, while the second one takes advantage of impedance matrices. Numerical evidence of the efficiency of the algorithms is given. Dispersion diagrams of simple cubic and woodpile photonic crystals are obtained by using, respectively, the scattering and impedance matrices.

© 2002 Optical Society of America

Stable and efficient Bloch-mode computational method for one-dimensional grating waveguides

Qing Cao, Philippe Lalanne, and Jean-Paul Hugonin
J. Opt. Soc. Am. A 19(2) 335-338 (2002)

Dispersion in photonic media and diffraction from gratings: a different modal expansion for the R-matrix propagation technique

J. Merle Elson and P. Tran
J. Opt. Soc. Am. A 12(8) 1765-1771 (1995)

Photonic bandgap calculations with Dirichlet-to-Neumann maps

Jianhua Yuan and Ya Yan Lu
J. Opt. Soc. Am. A 23(12) 3217-3222 (2006)

Anomalous refractive properties of photonic crystals

Boris Gralak, Stefan Enoch, and Gérard Tayeb
J. Opt. Soc. Am. A 17(6) 1012-1020 (2000)

Significance of Bloch impedance over wave impedance in photonic crystal waveguides

Chandra Prakash and Mrinal Sen
J. Opt. Soc. Am. B 38(6) 1997-2003 (2021)