Photonic Crystal - Computing Photonic Band Structure

Computing Photonic Band Structure

The photonic band gap (PBG) is essentially the gap between the air-line and the dielectric-line in the dispersion relation of the PBG system. To design photonic crystal systems, it is essential to engineer the location and size of the bandgap; this is done by computational modeling using any of the following methods.

1. Plane wave expansion method
2. Finite element method.
3. Finite difference time domain method
4. Order-n spectral method
5. KKR method
6. Bloch wave – MoM method

Essentially these methods solve for the frequencies (normal models) of the photonic crystal for each value of the propagation direction given by the wave vector, or vice-versa. The various lines in the band structure, correspond to the different cases of n, the band index. For an introduction to photonic band structure, see Joannopoulos.

The plane wave expansion method, can be used to calculate the band structure using an eigen formulation of the Maxwell's equations, and thus solving for the eigen frequencies for each of the propagation directions, of the wave vectors. It directly solves for the dispersion diagram. Electric field strength values can also be calculated over the spatial domain of the problem using the eigen vectors of the same problem. For the picture shown to the right, corresponds to the band-structure of a 1D DBR with air-core interleaved with a dielectric material of relative permittivity 12.25, and a lattice period to air-core thickness ratio (d/a) of 0.8, is solved using 101 planewaves over the first irreducible Brillouin zone.

In order to speed up the calculation of the frequency band structure, the Reduced Bloch Mode Expansion (RBME) method can be used. The RBME method applies "on top" of any of the primary expansion methods mentioned above. For large unit cell models, the RBME method can reduce the time for computing the band structure by up to two orders of magnitude.