Commun. Comput. Phys., 16 (2014), pp. 1355-1388.

An Efficient Calculation of Photonic Crystal Band Structures Using Taylor Expansions

Dirk Klindworth 1*, Kersten Schmidt 1

1 Research Center MATHEON, TU Berlin, Department of Mathematics, Strasse des 17. Juni 136, 10623 Berlin, Germany.

Received 24 May 2013; Accepted (in revised version) 26 June 2014
Available online 2 September 2014


In this paper we present an efficient algorithm for the calculation of photonic crystal band structures and band structures of photonic crystal waveguides. Our method relies on the fact that the dispersion curves of the band structure are smooth functions of the quasi-momentum in the one-dimensional Brillouin zone. We show the derivation and computation of the group velocity, the group velocity dispersion, and any higher derivative of the dispersion curves. These derivatives are then employed in a Taylor expansion of the dispersion curves. We control the error of the Taylor expansion with the help of a residual estimate and introduce an adaptive scheme for the selection of nodes in the one-dimensional Brillouin zone at which we solve the underlying eigenvalue problem and compute the derivatives of the dispersion curves. The proposed algorithm is not only advantageous as it decreases the computational effort to compute the band structure but also because it allows for the identification of crossings and anti-crossings of dispersion curves, respectively. This identification is not possible with the standard approach of solving the underlying eigenvalue problem at a discrete set of values of the quasi-momentum without taking the mode parity into account.

AMS subject classifications: 35Q61, 65N30, 65Z05, 78-04, 78M10
PACS: 42.70.Qs, 71.20.-b, 02.30.Jr, 02.70.Dh
Key words: Photonic crystals, band structure calculation, group velocity, Taylor expansion, adaptivity.

*Corresponding author.
Email: (D. Klindworth), (K. Schmidt)

The Global Science Journal