We derive a light-intensity-dependent dielectric constant for gain medium based on
the conventional rate equation model. A scattering-matrix method in conjunction with
an efficient iteration procedure is proposed to simulate photonic crystal lasers (PCLs).
The light output vs pumping (L-I) curve, lasing mode profile, and chirping effect of
lasing wavelength can be calculated. We check our method in a 1D DBR laser and the
L-I curve agrees well with results by the rate equation model. Our method can be
extended to 3D systems. More complex 2D and 3D PCLs will be simulated in the
The planewave based transfer matrix method has been developed with rational function interpolation to efficiently simulate photonic crystal devices. Cavities embedded in three-dimensional layer-by-layer photonic crystal are systematically studied as an example to show the power of transfer matrix method with the relation between resonant frequencies and the cavity size obtained.