A simple implementation of plane wave method is presented for modeling photonic crystals with arbitrary shaped ‘atoms’. The Fourier transform for a single ‘atom’ is first calculated either by analytical Fourier transform or numerical FFT, then the shift property is used to obtain the Fourier transform for any arbitrary supercell consisting of a finite number of ‘atoms’. To ensure accurate results, generally, two iterating processes including the plane wave iteration and grid resolution iteration must converge. Analysis shows that using analytical Fourier transform when available can improve accuracy and avoid the grid resolution iteration. It converges to the accurate results quickly using a small number of plane waves. Coordinate conversion is used to treat non-orthogonal unit cell with non-regular ‘atom’ and then is treated by standard numerical FFT.
Original Publication Citation
Guo, S., & Albin, S. (2003). Simple plane wave implementation for photonic crystal calculations. Optics Express, 11(2), 167-175. doi:10.1364/OE.11.000167
Guo, Shangping and Albin, Sacharia, "Simple Plane Wave Implementation for Photonic Crystal Calculations" (2003). Electrical & Computer Engineering Faculty Publications. 170.