OSA's Digital Library

Optics Express

Optics Express

  • Editor: C. Martijn de Sterke
  • Vol. 15, Iss. 16 — Aug. 6, 2007
  • pp: 10294–10302
« Show journal navigation

Lasing dynamics of a silicon photonic crystal microcavity

Shouyuan Shi and Dennis W. Prather  »View Author Affiliations

Optics Express, Vol. 15, Issue 16, pp. 10294-10302 (2007)

View Full Text Article

Acrobat PDF (247 KB)

Browse Journals / Lookup Meetings

Browse by Journal and Year


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools



In this paper we propose a novel silicon microcavity design based on the dispersion engineered photonic crystals (PhCs). With the unique self-collimation property of PhCs, we optimize the passive cavity by tuning the design parameters, such as coupling gap size and array size, to achieve higher Q factor and drop efficiency. Highest cavity mode below the band edge is of particular interest. The strong mode confinement in the low index active material offers an opportunity to realize a lasing mechanism. To investigate the lasing dynamics we introduce the rate equations of atomic system into the electromagnetic polarization to fully describe the nonlinearity of active medium. With these auxiliary differential equations we solve the time evolutions of the electromagnetic waves and atomic populations by using the FDTD method.

© 2007 Optical Society of America

1. Introduction

Silicon, the leading material in microelectronics during the last past decades, also promises to be the key material in the future. Silicon based microchip has been proven as an extremely powerful technology platform for a variety of applications of microelectronics, and now becomes the basis of complex microprocessor, larger memory circuits and other digital and analog electronics. The rapid progress of microprocessors may soon or later be severely limited by the transmission bandwidth capability of electronic connections. As a result, there is definitely a strong demand for fast communication links. To eliminate the bottleneck of electronic circuit and establish fast data communication between circuit boards, between chips on a board, or even with a single chip, the research of silicon photonics has attracted more and more attention due to its intrinsic high data transmission rate.

To this end, we will design a novel optical cavity based on the dispersion engineered photonic crystals (PhCs). To observe the lasing dynamics, we will incorporate the rate equations of a four level atomic system into the device to simulate the gain and absorption of the active material. By solving the Maxwell’s equations with these auxiliary differential equations with the Finite-difference Time-domain (FDTD) method, we will track the time evolutions of the electromagnetic waves and atomic populations.

2. Study of dispersion based photonic crystal cavity

The light propagation in a PhC is most appropriately interpreted through a dispersion diagram, which characterizes the relationship between the frequency of the wave, ω, and its associated wavevector, k. Dispersion surfaces provide the spatial variation of the spectral properties of a certain band within the photonic crystal structure. Electromagnetic wave propagates along the direction normal to the dispersion surface as shown in Fig. 1, which stems from the relation of the group velocity v g=∇k ω(k). The equi-frequency contour (EFC), which can be obtained from the dispersion surface at certain frequency, can lead to beam divergence or convergence as shown in the figure. The ability to shape the EFCs, and thereby engineer the dispersion properties of the PhC, opens up a new paradigm for the design of optical devices [5

5. D. W. Prather, S. Shi, D. Pustai, C. Chen, S. Venkataraman, A. Sharkawy, G. Schneider, and J. Murakowski, “Routing Optical Waves Without Waveguides,” Opt. Lett. 29, 50–52 (2004) [CrossRef] [PubMed]


7. K. K. Tsia and A. W. Poon, “Dispersion-guided resonances in two dimensional photonic-crystal embedded microcavities,” Opt. Express 12, 5711–5722 (2004) [CrossRef] [PubMed]

]. For the applications of self-collimation, we desire a flat EFC, in which case the wave is only allowed to propagate along those directions normal to the sides of the straight curvatures. As such, it is possible to vary the incident wavevector over a wide range of angles and yet maintain a narrow range of propagating angles within the PhC. Based on the dispersion waveguiding property in photonic crystals, a novel class of photonic-crystal-embedded microcavity (PCEM) coupled with waveguide [7

7. K. K. Tsia and A. W. Poon, “Dispersion-guided resonances in two dimensional photonic-crystal embedded microcavities,” Opt. Express 12, 5711–5722 (2004) [CrossRef] [PubMed]

] has been theoretically investigated by using FDTD algorithm, in which both optical resonant mode and quality factor are particularly considered.

Consider a silicon photonic crystal slab perforated by a square lattice with air holes back-filled by the gain medium, i.e. Er-doped glasses, as shown in the inset of Fig.2 (a). The hole has radius of 0.3a, where a is the lattice constant. The silicon and glass have refractive indices of 3.5 and 1.5, respectively. The dispersion surface of first band diagram is plotted in Fig. 1(a). By carefully selecting the frequency, one can obtain a flat curvature within certain angular range at specified frequency, i.e. 0.18c/a, as depicted in Fig.1 (b), where the blue circle is the dispersion contour in free space with frequency of 0.18c/a. Such flatness of the curvature offers self-collimation along ΓM direction within a wide incident angular range.

Fig. 1. (a) The dispersion surface of the first band for silicon PhC with circular holes filled with Er-doped glass. (b) Flat EFCs perpendicular to ΓM direction at frequency of 0.18c/a.

Fig. 2. (a) Schematic of silicon laser based on dispersion engineered photonic crystals where the active material is introduced by backfilling the air holes of the PhC. (b) The highest cavity mode below the band edge.

We simulate the proposed device by using Finite-difference time-domain (FDTD) method. 6×6 array of glass filled cylinders is initially simulated; the highest resonant mode below the band edge of first dispersion band, is shown in Fig. 2(b). As we can observe that most optical mode are well confined within the low index materials. This unique property may benefit the lasing mechanism and lower the threshold optical pumping. The cavity property is largely dependent on the design parameters, such as gap size between the waveguide and resonator, and array size. Extensive study has been done to investigate the passive cavity performance. The gap between the cavity and waveguide is critical design parameter to achieve optimal coupling. First, we fixed the number of array and tune the gap size to study the Q factor and drop efficiency from the waveguide to the resonator. The gap is continually changed from the 150nm to 300nm. As we can see from Fig. 3 that the Q factor increases as the gap size increases due to reduction of the coupling loss, while the drop efficiency expectedly decreases. An optimal gap size of 160nm is chosen in the following design. In this case we measured Q factor of 520 and drop efficiency of 32.2%. In addition, we fix the gap size and increase the number of array size from 6 to 11 and find as the number of array increases the resonant frequency slightly shift to high end, and the Q-factor increase accordingly, which offer us another degree of freedom to tune the Q factor. In the case of array size of 10, the Q factor is 1900 and the drop efficiency remains as high as up to 56%.

Fig. 3. Q factor and drop efficiency as a function of gap size between the waveguide and resonator.

With the appropriate design of the microcavity based on the dispersion engineering of photonic crystals, we will further consider to backfill the air holes in photonic crystals with gain medium to achieve active operation.

3. Rate equations of a four-level atomic system

To simulate a laser dynamics in an optical cavity associated with a gain or active material, a conventional rate equation model [9

9. S. Chang and A. Taflove, “Finite-difference time-domain model of lasing action in a four-level two-electron atomic system,” Opt. Express 12, 3827–3833 (2004) [CrossRef] [PubMed]

] can be employed for the simulations of the nonlinearity and dispersion properties, i.e. a uniform glass host containing a dispersion of highly doped erbium ions. Coupled with Maxwell’s equations, we are able to model the time evolution of the atomic energy level populations as well as the optical signal propagation, amplification and absorption in the devices. In this paper we proposed a simplified realistic four-level atomic system as sketched in Fig. 3, where the principal transitions induced by the presence of the pump and signal beams are included. The absorption of pump radiation at 980nm wavelength promotes the electrons from level 0 to level 3. After decaying to level 2, these electrons provide the signal gain at 1500nm via the transition to level 1. Here we describe the time domain population dynamics by using a four-level rate equation formulism:


Fig. 4. Populations in the simplified four-level atomic system.

Based on classical electron oscillator model, the net macroscopic polarization P(t) induced with the presence of applied electric field E(t) for an isotropic medium can be described by the following equations,


where ΔN(t) is instantaneous population density difference between energy levels 1 and 2 of atomic transition, which is given by ΔN(t)=N1(t)-N2(t). κ=γrc, where γr=1/τ21 is the radiative decay rate for this from E2 and E1 transition and γc=e 2 ω 2 a/(m6πε0c3) is the classical rate. Δωa is the total energy decay which describes the transition linewidth. It can be easily derived from the above equation that the amplification line shape Lorentizen and homogeneous broadened, and can be considered as a quantum mechanically correct equation for the induced polarization density in a real atomic system. To excite the system, an optical pumping source is introduced into system, where the source is homogeneously distributed across the device. The amplification takes place when the external pumping mechanism produces population inversion ΔN(t)<0.

By doing so, then the effects of the nonlinear and active medium on the propagation of electromagnetic waves are taken into account through the polarization response of the medium. Coupled with the Maxwell’s equation in our optical devices, we are able to simulate the proposed photonic crystal laser cavity and more appropriate to understand the time evolution of such system.

4. Dynamics of electromagnetic wave in cavity system

Based on the theory we developed, in this section we will introduce Finite-difference Time-domain (FDTD) method to solve the coupled equations of atomic system and Maxwell’s equations for the application of lasing dynamics in the microcavity system.

As is well known, the FDTD method [10

10. A. Taflove, “Computational Electromagnetics: The Finite-Difference Time Domain Method,” (1995)

] has become one of the more widely used numerical techniques for solving electromagnetic boundary value problems. In short, the FDTD method directly discretize Maxwell’s equations in both time and space by using central difference technique, and iteratively updates the electrical and magnetic field components at all points in the computational space using the leapfrog fashion. The complicated material properties can be also incorporated in computation. To this end, the method is able to easily deal with the interactions between electromagnetic fields and objects having complex shape and/or inhomogeneous, dispersive and passive/gain media.

As shown in the schematic view of proposed novel silicon laser in Fig.2, the holes are back filled with an active material, i.e. Er doped Glass. The medium is uniformly pumped with an external source. For the Er3+ Ion, the typical lifetime is on the order of 10-3~10-2s. Numerically, with rigorous EM algorithm in time domain to simulate the lasing dynamics is far beyond the state of art computational capability, particularly in high dimension system. To phenomenally investigate the lasing dynamics in the microcavity with the backfilled gain medium, we scale the lifetime of electron population transition accordingly and yet maintain the reasonable timescale associated to different relaxation process. To this end we shorten the lifetime of atomic transition to reduce the computation time that needed to achieve the steady state. The lifetime are given by τ10=10-12 s, τ 21=10-10 s, and τ32=10-13s. The transition frequency associated with the energy levels E2 and E1 is chosen as 200THz and linewidth is taken to be 6THz. The pump rate into level E3 is chosen Wp = 2×108/s. The initial state of the simulation is included as follow: all the electrons are on the ground state, so there is no field in the cavity and no spontaneous emission. In this paper we choose N0=3×1024/m3. After the electrons are pumped, the system start to evolve both in population density and electromagnetic waves. To monitor the EM dynamics a detector is placed in the waveguide.

In the FDTD simulation, a total time of 66ps is simulated. The time step Δt is chosen as 3.3×10-17 s, so total 2 million FDTD iterations is used. The evolution of electromagnetic waves starts from an initial small noise over the photonic crystal cavity. A detector is placed in the waveguide to monitor the lasing dynamics, and the total field in the output straight waveguide is illustrated in Fig. 5(a). In the figure we can obviously lasing performance. For the clarity, we zoom the transient plot near 60ps as shown in Fig. 5(b), which indicates nearly single lasing model.

Fig. 5. (a) Lasing dynamics by monitoring output in the straight dielectric waveguide, (b) steady output of single optical mode as shown within the time window in figure (a).

In addition, snapshots of 2D magnetic fields at three instant time of 10, 30 and 50ps are plotted in Fig. 6, where the lasing dynamics of electromagnetic field is gradually established within the cavity. At last, the high intensity distribution of EM field can be observed in the cavity area, particularly in the low index material region and steady results are outputted from both ends of straight waveguide.

Fig. 6. Snapshots of 2D field distribution at three instant time of 10, 30 and 50ps.

Figure 7 shows the calculated time evolution of the electron population inversion between energy level 2 and 1 at the position of one of central holes in the cavity. The plotted population inversion is normalized to the initial population density of the ground state N0. During the simulation of the population evolution, the ground state population density almost remains unchanged compared to higher-level populations. In the beginning the population inversion linearly increases, so it leads to significant amplification of electromagnetic wave in the cavity. In the meanwhile the desired cavity mode is established gradually. At last, a convergent population inversion, -6.9×10-4, can be observed. The population dynamics are well consistent with the lasing output from waveguide. In addition to this, we also calculate the optical output intensity in the waveguide by varying the pumping rate from 0.05×108/s to 8×108/s to study the lasing threshold. To obtain a stable lasing behavior, particularly near the threshold point, due to lower available gain a longer simulation time is required to actively excite the optical mode inside the microcavity. At the pumping rate of 0.05×108/s, 24 million FDTD iterations are used. Figure 8 shows the optical output intensity as the function of the pumping rate. The lasing threshold can be predicted from the output optical intensity as the pumping rate reduces. As shown in the inset of Fig. 8, and we have measured the lasing threshold to be 0.07×108/s for the system.

Fig. 7. Normalized population inversion at pumping rate of 2×108/s.
Fig. 8. Output intensity as the function of the pumping rate.

5. Conclusion

In this paper we proposed a novel silicon laser design based on the dispersion engineering of photonic crystals. By carefully designing the photonic crystals, we took advantages of the self-collimation property to suppress the optical cavity modes. The relative high concentration of optical mode in lower index gain material will offer benefit to improve the lasing performance. Detailed studies have been performed to investigate the passive cavity design in terms of the coupling gap and the array size of air holes to maximize Q factor while maintains the high drop efficiency. To simulate the lasing dynamics in such optical cavity, we incorporated the rate equations of a four-level atomic system into Maxwell’s with the assistance of the induced electric polarization contributed from the response of a collection of atoms. Finally, we numerically implemented time evolutions of electromagnetic waves and atomic populations by using FDTD method to the application of our engineered dispersion based photonic crystal laser.

References and links


B. Jalali, M. Paniccia, and G. Reed, “Silicon photonics,” IEEE Microwave Magazine 7, 58–68 (2006) [CrossRef]


L. Pavesi, “Will silicon bethe photonic material of the third millenium,” J. Phys. Condens. Matter 15, R1169–R1196 (2003) [CrossRef]


H. Rong, R. Jones, A. Liu, O. Cohen, D. Hak, A. Fang, and M. Paniccia, “A continues-wave Raman silicon laser,” Nature 433, 725–728 (2005) [CrossRef] [PubMed]


O. Boyraz and B. Jalali, “Demonstartion of a silicon Raman laser,” Opt. Express 12, 5269–5273 (2004) [CrossRef] [PubMed]


D. W. Prather, S. Shi, D. Pustai, C. Chen, S. Venkataraman, A. Sharkawy, G. Schneider, and J. Murakowski, “Routing Optical Waves Without Waveguides,” Opt. Lett. 29, 50–52 (2004) [CrossRef] [PubMed]


M. R. Newton, K. A. Morey, Y. H. Zhang, R. J. Snow, M. Diwekar, J. Shi, and H. S. White, “Anisotropic diffusion in face-centered cubic opals,” Nano. Lett. 4, 875–880 (2004) [CrossRef]


K. K. Tsia and A. W. Poon, “Dispersion-guided resonances in two dimensional photonic-crystal embedded microcavities,” Opt. Express 12, 5711–5722 (2004) [CrossRef] [PubMed]


D.W. Prather, A. Sharkawy, and S. Shi, “Handbook of Nanoscience, Engineering and Technology,” 1,b (2002)


S. Chang and A. Taflove, “Finite-difference time-domain model of lasing action in a four-level two-electron atomic system,” Opt. Express 12, 3827–3833 (2004) [CrossRef] [PubMed]


A. Taflove, “Computational Electromagnetics: The Finite-Difference Time Domain Method,” (1995)


A. Nagra and R. A. York, “FDTD analysis of wave propagation in nonlinear absorbing and gain media,” IEEE Trans. Antennas Propag. 46, 334–340 (1998) [CrossRef]


X. Jiang and C. M. Soukoulis, “Time dependent theory for random lasers,” Phys. Rev. Lett. 85, 70–73 (2000) [CrossRef] [PubMed]


P. Sebbah and C. Vanneste, “Random laser in the localized regime,” Phys. Rev. B 66, 144202 (2002) [CrossRef]

OCIS Codes
(000.4430) General : Numerical approximation and analysis
(140.3380) Lasers and laser optics : Laser materials
(140.4780) Lasers and laser optics : Optical resonators
(190.4360) Nonlinear optics : Nonlinear optics, devices

ToC Category:
Photonic Crystals

Original Manuscript: May 30, 2007
Revised Manuscript: July 23, 2007
Manuscript Accepted: July 23, 2007
Published: July 31, 2007

Shouyuan Shi and Dennis W. Prather, "Lasing dynamics of a silicon photonic crystal microcavity," Opt. Express 15, 10294-10302 (2007)

Sort:  Year  |  Journal  |  Reset  


  1. B. Jalali, M. Paniccia, and G. Reed, "Silicon photonics," IEEE Microw. Mag. 7, 58-68 (2006). [CrossRef]
  2. L. Pavesi, "Will silicon bethe photonic material of the third millenium," J. Phys.: Condens. Matter 15, R1169-R1196 (2003). [CrossRef]
  3. H. Rong, R. Jones, A. Liu, O. Cohen, D. Hak, A. Fang, and M. Paniccia, "A continues-wave Raman silicon laser," Nature 433, 725-728 (2005). [CrossRef] [PubMed]
  4. O. Boyraz and B. Jalali, "Demonstartion of a silicon Raman laser," Opt. Express 12, 5269-5273 (2004). [CrossRef] [PubMed]
  5. D. W. Prather, S. Shi, D. Pustai, C. Chen, S. Venkataraman, A. Sharkawy, G. Schneider, and J. Murakowski, "Routing optical waves without waveguides," Opt. Lett. 29, 50-52 (2004). [CrossRef] [PubMed]
  6. M. R. Newton, K. A. Morey, Y. H. Zhang, R. J. Snow, M. Diwekar, J. Shi, and H. S. White, "Anisotropic diffusion in face-centered cubic opals," Nano. Lett. 4, 875-880 (2004). [CrossRef]
  7. K. K. Tsia and A. W. Poon, "Dispersion-guided resonances in two dimensional photonic-crystal embedded microcavities," Opt. Express 12, 5711-5722 (2004). [CrossRef] [PubMed]
  8. D. W. Prather, A. Sharkawy, and S. Shi, Handbook of Nanoscience, Engineering and Technology (CRC Press, 2002) 1, b.
  9. S. Chang and A. Taflove, "Finite-difference time-domain model of lasing action in a four-level two-electron atomic system," Opt. Express 12, 3827-3833 (2004). [CrossRef] [PubMed]
  10. A. Taflove, Computational Electromagnetics: The Finite-Difference Time Domain Method (Artech House, Boston 1995).
  11. A. Nagra and R. A. York, "FDTD analysis of wave propagation in nonlinear absorbing and gain media," IEEE Trans. Antennas Propag. 46, 334-340 (1998). [CrossRef]
  12. X. Jiang and C. M. Soukoulis, "Time dependent theory for random lasers," Phys. Rev. Lett. 85, 70-73 (2000). [CrossRef] [PubMed]
  13. P. Sebbah and C. Vanneste, "Random laser in the localized regime," Phys. Rev. B 66,144202 (2002). [CrossRef]

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited