Calculating Defect Modes in 2D Photonic Crystals Using the Plane Wave Expansion Method

The plane wave expansion method serves as a classical numerical approach for calculating photonic band structures, particularly effective for solving electromagnetic wave propagation in periodic dielectric structures. For defect mode calculations in 2D photonic crystals, this method employs Fourier expansions of electromagnetic fields and dielectric constants in reciprocal lattice space, transforming Maxwell's equations into eigenvalue problems to solve for defect states within specific frequency ranges. In computational implementation, the process begins with constructing a complete periodic photonic crystal structure, determining lattice type and medium distribution. The introduction of defects (point defects or line defects) then breaks the periodicity to form localized states. Through plane wave expansion, both dielectric constants and electromagnetic fields are expressed as Fourier series in reciprocal space, converting the problem into matrix eigenvalue solving. Key computational steps involve implementing Fourier transform algorithms and building the eigenvalue matrix using reciprocal lattice vectors. Critical aspects in defect mode calculation include: - Selecting appropriate plane wave truncation numbers to ensure computational accuracy - Correctly handling dielectric function modifications after defect introduction - Solving eigenvalue equations to obtain defect mode frequencies and field distributions - Analyzing quality factors and localization characteristics of defect modes This method effectively reveals how defects influence photonic band structures, providing theoretical foundations for designing optical devices like photonic crystal waveguides and resonators. Computational results are typically presented as band diagrams and field distribution plots, visually demonstrating the position and characteristics of defect-induced localized states within band gaps. The implementation typically involves constructing a Hamiltonian matrix using Fourier coefficients of the dielectric function and solving it using numerical eigenvalue solvers like MATLAB's eigs function or LAPACK routines.