Temporal behaviour of field in high quality factor photonic crystal microcavity structure

Temporal behaviour of incident pulse in high-quality (Q) factor photonic crystal microcavities are studied by two dimensional finite difference time domain calculations. For high-Q mode excitation, two periods of oscillation are observed in addition to the exponential decay corresponding to the cavity mode photon life time. Long and short period oscillations correspond to beats with low-Q mode (for short pulse widths) and to lattice periodicity for all pulse widths, respectively. For low-Q mode and off-resonant excitations, long period oscillations correspond to coupling to bandedge states. © 2005 Optical Society of America OCIS codes: (160.6000) Semiconductors, including MQW, (320.7120) Ultrafast phenomena, (230.1150) All-optical devices, (060.1810) Couplers, switches, and multiplexers References and links 1. J. Vuc̆kovic, and Y.Yamamoto, ”Photonic crystal microcavities for cavity quantum electrodynamics with a single quantum dot,” Appl. Phys. Letts. 82, 2374-2376 (2003). 2. K.J.Vahala, ”Optical microcavities,” Nature 424, 839-846, (2003). 3. D.G.Angelakis, E. Paspalakis, and P.L.Knight, ”Coherent phenomena in photonic crystals,” Phys. Rev. A 64, 013801, (2001). 4. K. Srinivasan, P. E. Barclay, O.Painter, J. Chen, A. Y. Cho, C. Gmachl, ”Experimental demonstration of a high quality factor photonic crystal microcavity,” Appl. Phys. Lett. 83, 1915-1917 (2003). 5. Y.Akahane, T.Asano, B-S.Song, and S.Noda, ”High-Q photonic nanocavity in a two-dimensional photonic crystal,” Nature 425, 944-947, (2003). 6. R.Colombelli etal, ”Quantum cascade surface-emitting photonic crystal laser,” Science 302, 1374-1377, (2003). 7. Y.Xu, R.K.lee, and A.Yariv, ”Propagation and second-harmonic generation of electromagnetic waves in a coupled-resonator optical waveguide,” J. Opt. Soc. Am. B 17, 387-400 (2000). 8. S.F.Migaleev, and Y.S.Kivshar, ”Nonlinear transmission and light localization in photonic-crystal waveguides,” J.Opt. Soc. Am B 19., 2241-2249 (2002). 9. V.M.Aplakov, and M.E.Raikh, ”Strongly localized mode at the intersection of the phase slips in a photonic crystal without band gap,” Phys. Rev. Lett. 90, 253901, (2003). 10. P.R.Villeneuve, S.Fan, S.Johnson, and J.D.Joannopoulos, ”Three-dimensional photon confinement in photonic crystals of low-dimensional periodicity,” IEE Proc. Optoelectron. 145, 384-390, (1998). 11. G.Subramania, S.Y.Lin, J.R.Wendt, and J.M.Rivera, ”Tuning the microcavity resonant wavelength in a twodimensional photonic crystal by modifying the cavity geometry,” Appl. Phys. Lett. 83, 4491-4493, (2003). 12. H-Y Ryu, M.Notomi, and Y-H. Lee, ”High-quality-factor and small-mode-volume hexapole modes in photoniccrystal-slab nanocavities,” Appl. Phys. Lett. 83, 4294-4296, (2003). (C) 2005 OSA 24 January 2005 / Vol. 13, No. 2 / OPTICS EXPRESS 460 #6073 $15.00 US Received 16 December 2004; revised 9 January 2005; accepted 10 January 2005 13. T. Asano, W. Kunishi, K. Kiyota, B-S. Song, A. Noda, ”Dynamic properties of artificial defects in twodimensional photonic crystal slab,” Tech. Digest, Intl. Symposium on Photonic and Electromagnetic Crystal Structures (PECS V), March 2004, Kyoto, Japan, p166. 14. J. Erland, and I. Balslev, ”Theory of quantum beat and polarization interference in four-wave mixing,” Phys. Rev. A 48, R1765-R1768 (1993).


Introduction
Photonic crystal (PC) microcavities that can exhibit high quality factor (Q) and small modal volumes (V) are considered ideal for both fundamental cavity Quantum Electrodynamics (cQED) and various device applications [1,2,3,4]. With recent demonstration of high-Q cavities [4,5,6], it may be interesting to study the possible effects of these structures on the incident field. In order to characterize the microcavities, transmission measurements to estimate the coupling efficiency and the spectrum are studied. The cavity Q is determined from the spectral width (∆λ ) of the transmission spectrum (Q≈ λ /∆λ for the cavity mode resonance at λ ). While this is valid for a low-Q cavity, in general, the cavity Q may be obtained, by estimating the cavity photon life time (τ) from Q = ω 0 τ where ω 0 is the resonant cavity mode frequency. However, estimation of τ by the use of a short pulse may not be trivial due to the possible excitation of multiple modes (due to the finite spectral width of the short pulse). Thus, the validity of estimation of Q-factor from the cavity mode photon life time using a short pulse is suspected. In addition, to study the feasibility of applying these cavities for quantum information processing and photonic switching applications, it would be interesting to study the effect of cavity on the time dependence of the incident field. Thus, we numerically studied the pulse width dependence of the pulse evolution. Though this may be of interest, to our knowledge, the time dependence of the incident pulse profile is not studied till date. Considerable work has been reported on the nonlinear propagation in coupled-resonator optical waveguide (CROW) structures [7,8]. However, the pulse profile is not reported in these structures also.
Several high-Q PC microcavity designs based on fractional edge dislocations or phase slips, graded PC lattices, super defects, and hexapole modes have been reported till date [1,9,4,10,11,12]. Though super defects, with modified set of holes as the defect, may lead to high-Q but at the expense of large modal volume (V) [10,11]. Recent demonstration shows a defect with 3 missing holes and suppressed edge effects leading to a large Q of about 45,000 and a reasonable mode volume [5]. However, these may not be suitable to introduce a nonlinear medium with refractive index different from the substrate material or Quantum Dots (QDs) or atoms for cQED studies. Similarly, the high-Q and small V hexapole modes suffer from field cancellation (a node for the electric field) at the center of the defect [12]. To overcome these aspects, it would be essential to design PC structures that demonstrate high-Q and small V and which can employ an acceptor defect with variable refractive index.
Defects in photonic crystals are generally referred to as acceptor (reduced air hole regions) and donor type (increased air hole region) in a regular PC lattice. While acceptor defects are found to confine the field in the dielectric regions, the mode volume is increased [1]. Figure 1 shows a schematic of the 2D PC structure, being studied here, that is formed by a hexagonal lattice of air holes defined in a dielectric medium. In this structure, we have increased the radius of the central hole but filled it with a dielectric material to form the central defect where the incident field is concentrated. Thus, we effectively have an acceptor type defect (increased dielectric region or reduced air region) so that the field is confined in the dielectric region. In addition, the size and shape of the neighbouring holes are modified to reduce the in-plane radiation losses to increase the Q-factor. We studied the dependence of Q on the radius and the refractive index of the dielectric defect material (central hole) to further optimize the structure for high-Q. The optimum defect radius is about 0.2µm (1.2×r, r is the radius of regular air hole) and the defect refractive index is about 3.9 (Sb-based semiconductor materials for GaAs PCs and PbSe or PbS for PCs fabricated on Si-slab). Details of these studies will be presented elsewhere. Briefly, from the calculated mode profiles, the Q value peaks when the defect size is comparable to that of the mode (about 1.2×r). That is, the mode is nicely confined in the defect region. For a reasonable 2D PC of size (15a) 2 , a being the lattice periodicity, an inplane Q-factor of about 1.7×10 5 and a small mode volume, in terms of the cubic wavelength inside the medium, of about 0.42×(λ /n) 3 were observed. These designs were extended to 3D structures and a Q-factor of about 5,000 was achieved by modifying the neighbouring holes to reduce the vertical radiation losses. In this paper, we present two-dimensional finite difference time domain (FDTD) calculations on this structure to study the time dependence of the incident power to study the temporal behaviour of the pulse profile.

Results
Temporal studies in photonic crystals performed till date are to study the group velocity and more recently the cavity Q [13]. In the two-dimensional PC structure studied here, the high-Q features are achieved by modifying the central hole (defect) and the six nearest neighbour holes (Fig. 1). Figure 2 shows the band structure of the 2D structure of air holes for transverse magnetic mode excitation . The band structure calculations show two defect modes at 1.53 and 1.414 µm wavelengths in the PBG region. This is verified to be equivalent to a 3D structure of air holes in a Si-slab of thickness 0.2µm. Figure 2 also shows the transmission spectra with the photonic bandgap (PBG) in the region of 1.3 to 1.8 µm and a defect mode peak at about 1.53µm. The second peak at 1.414 µm, however, is not seen in the spectrum due either to a weak peak or due to the position of the detector. It would further be interesting to note that while the cavity mode photon life time may be estimated unambiguously from the temporal data (inset of Figure 3(a)), to calculate the spectral width accurately, one requires a longer calculation. For temporal studies, we varied the gaussian pulse width in the range of 16fsec to 3psec to study the effect of pulse width on the cavity mode photon lifetime and the field (due to the possible excitation of multiple modes). By using 2D FDTD calculations, we calculated the transmitted power and the power at the center of the PC structure (defect region). The grid size employed in these calculations is 0.01x0.01 µm. Broad pulse calculations indicate a high-Q at 1.53µm wavelength. Figure 3(a) compares the time evolution of the transmitted pulse intensity for two pulse widths resonant at 1.53 µm. While the 100fsec pulse excitation shows a non-oscillatory, single exponential decay in the long time limit (inset of Fig. 3(a)) with a decay constant of about 133psec, the 17fsec pulse excitation shows periodic oscillations even in the long time limit.  Considering that these oscillations correspond to beats between two states, the period of oscillation corresponds to 33.3fsec (cτ ≈ 10µm) or to an energy separation of about 63meV. The solid line in Fig. 3(b) is a fit based on an expression for beats for a two-level system given by, exp(−tγ 1 )+α 2 exp(−tγ 2 )+2αCos(t∆E)exp(−0.5t(γ 1 + γ 2 )) where γ 1 and γ 2 are the decay rates of 1.53 µm and 1.414 µm modes, ∆E is the energy separation between the two modes, and α is the coupling strength. The fit parameters are 133 psec, 100 fsec, 65 meV, and 0.04 for γ −1 1 , γ −1 2 , ∆E, and α, respectively. As shown in Fig. 2 between the two modes at 1.53 and 1.414 µm matches quite well with that estimated from the fit and the observed period of oscillation in Fig. 3(b) indicating that the oscillations observed for short pulse excitation could be due to beats between these two modes. We further verified this by studying the pulse width dependence. These studies show that the oscillations are observed only for pulse widths shorter than 90fsec. This is consistent with the fact that for pulse widths shorter than 90fsec, the spectral width is broad enough to simultaneously excite both the modes and we observe the oscillations in the time evolution similar to the oscillations observed in multi-level electronic, excitonic or atomic systems [14]. We also performed similar studies on the short wavelength mode. Figure 4 shows the calculated time evolution of transmitted field intensity for 0.5psec broad pulse that is resonant at the 1.414 µm mode in a PC structure of size 12a 2 . It may be interesting to note that, though the pulse width is broad (thus the spectral width is narrow enough to not excite the 1.53µm mode), we see oscillations in the time evolution. The time decay shows two exponential decays in addition to these oscillations. The dashed lines shown in Fig. 4 are exponential fits. The short time decay constant is about 0.68 psec. Eventhough most of the field decayed before the long time limit exponential decay begins, the long time decay constant shows about 22.3psec decay constant similar to that observed for 1.53 µm mode. When we increase the size of the PC structure from (12a) 2 to (15a) 2 , the short time exponential decay constant at 1.414 µm falls from 0.68 psec to 0.45 psec (within the pulse width). The long period time constant of 130 psec is still similar to that for the 1.53µm mode. Though we do not understand the origin of this long time decay of 1.414 µm, it may be very interesting to study the possibility of coupling to the long wavelength mode.
Inset of Fig. 4 shows the time evolution for a short pulse (width is 17fsec) excitation resonant at 1.414µm. As may be expected, we observe oscillations but a non-exponential long time limit decay. A closer inspection shows two time periods. The longer oscillation period matches with the beat frequency corresponding to the simultaneous excitation of both modes as observed for the short pulse excitation at 1.53 µm. However, the decay is very fast unlike for the 1.53 µm excitation case. It may be interesting to note that the 1.414 µm mode is separated by similar energy from both the high energy PBG edge at 1.3 µm (0.954 eV) and the 1. mode. Consider that there are two decay paths for the 1.414 µm mode with time scales T 1 and T 2 corresponding to coupling to 1.53µm and PBG edge, respectively. The decay time (T) of the 1.414 µm mode would be given by (C 1 /T 1 )+(C 2 /T 2 ) where C 1 and C 2 are the coupling strengths. From the above results we see that the T 1 is about 130psec and the T is about 100fsec (comparable to the pulse width). So, for coupling to the 1.53 µm mode to be non-negligible and the decay to be a slower one, C 1 has to be, atleast, a 100 times larger than C 2 (coupling to the band edge states). Thus, from the fast decay time constant one may infer that the radiative coupling to the PBG edge may be more effective than coupling to the high-Q mode. In other words, the observed long period oscillations for the low-Q mode excitation could correspond to the coupling to the PBG edge.  Fig. 5 shows the corresponding close up of short period oscillations for short pulse excitation resonant at 1.414 µm with the oscillation period (0.45±0.05 µm) similar to that observed at 1.53 µm. This indicates that the short period oscillations observed could be related to the lattice periodicity (a=0.49µm).
We also studied the time evolution of field at two non-resonant wavelengths above and below the 1.53 µm wavelength. Fig. 6 shows the time dependence for 1.5 and 1.57 µm excitation. The pulse width of 167 fsec is chosen such that the two defect modes are not within the pulse spectral width. However, a small component (at 1.53 µm) in the tail of the Gaussian excitation pulse excites the 1.53 µm mode. Correspondingly, a very weak peak at 1.53 µm is visible in the fourier transform of the temporal data shown in Fig. 6. For the purpose of studying the coupling behaviour, from the temporal data we evaluate the beat frequency and the corresponding energy separation. This information is not available from the spectrum (fourier transform of the temporal data). As expected, the time evolution shows a fast decay at both the non-resonant wavelengths. In addition, the beat period of about 0.6psec for 167 pulse excitation (solid line in Fig. 6) corresponds to an energy separation of about 93 meV for 1.5µm (0.827 eV) excitation indicating beats with the low energy edge of the PBG at ≈1.8 µm (0.689 eV). Though a short pulse width (67 fsec) at these off-resonant wavelengths is considered such that the high-Q mode  Fig. 6) and short (67fsec) pulse excitations at 1.57 µm (0.79eV), long period oscillations correspond to beats with the low energy PBG edge. Clearly, the off-resonant and the resonant low-Q mode excitations, for all pulse widths, show a fast decay (similar to the pulse width). Only the resonant excitation of the high-Q mode at 1.53 µm shows a long exponential decay in the long time limit. This addresses, thus, the question about the use of short pulse to estimate the τ of a high-Q mode.

Conclusion
In conclusion, we have studied, by 2D FDTD calculations, the time dependence of the incident pulse in a high-Q photonic crystal cavity. These studies give an insight into the coupling behaviour of defect modes in high-Q PC cavities. Though the pulse spectral width covers both the modes, only the pulse resonant at the high-Q mode shows the long cavity mode decay in addition to the oscillations corresponding to beats with the low-Q defect mode. While long period oscillations in case of simultaneous excitation of both the defect modes correspond to beats, the short period oscillation for all pulse widths in the range 16fsec to 3 psec correspond to the lattice periodicity. Low-Q mode and off-resonant excitations couple to the band edge states and decay within the short pulse width. From the beat frequency (long period oscillations) and the fast decay, it may be inferred that the low-Q modes in this cavity couple to the band edge states. It may be interesting to study the effect of pulse widths in structures in which the two defect modes are close to each other and far away from the bandedge. These results also indicate that in spite of the apprehensions in using the more general definition based on the cavity mode life time to estimate Q, in addition to the oscillations corresponding to beats, only the resonant excitation of the high-Q mode shows the long time exponential decay corresponding to the cavity mode life time.