## Low-threshold lasing action in photonic crystal slabs enabled by Fano resonances |

Optics Express, Vol. 19, Issue 2, pp. 1539-1562 (2011)

http://dx.doi.org/10.1364/OE.19.001539

Acrobat PDF (2405 KB)

### Abstract

We present a theoretical analysis of lasing action in photonic crystal surface-emitting lasers (PCSELs). The semiclassical laser equations for such structures are simulated with three different theoretical techniques: exact finite-difference time-domain calculations, an steady-state ab-initio laser theory and a semi-analytical coupled-mode formalism. Our simulations show that, for an exemplary four-level gain model, the excitation of dark Fano resonances featuring arbitrarily large quality factors can lead to a significant reduction of the lasing threshold of PCSELs with respect to conventional vertical-cavity surface-emitting lasers. Our calculations also suggest that at the onset of lasing action, most of the laser power generated by finite-size PCSELs is emitted in the photonic crystal plane rather than the vertical direction. In addition to their fundamental interest, these findings may affect further engineering of active devices based on photonic crystal slabs.

© 2011 OSA

## 1. Introduction

1. T. H. Maiman, “Stimulated optical radiation in ruby,” Nature **187**, 493–494 (1960). [CrossRef]

2. E. Yablonovitch, “Inhibited spontaneous emission in solid-state physics and electronics,” Phys. Rev. Lett. **58**, 2059–2062 (1987). [CrossRef] [PubMed]

5. H. Yokoyama, “Physics and device applications of optical microcavities,” Science **256**, 66–70 (1992). [CrossRef] [PubMed]

6. N. M. Lawandy, R. M. Balachandran, A. S. L. Gomes, and E. Sauvain, “Laser action in strongly scattering media,” Nature **368**, 436–438 (1994). [CrossRef]

32. M. P. Nezhad, A. Simic, O. Bondarenko, B. Slutsky, A. Mizrahi, L. Feng, V. Lomakin, and Y. Fainman, “Room-temperature subwavelength metallo-dielectric lasers,” Nat. Photonics **4**, 395–399 (2010). [CrossRef]

37. S. Fan and J. D. Joannopoulos, “Analysis of guided resonances in photonic crystal slabs,” Phys. Rev. B **65**, 235112 (2002). [CrossRef]

8. O. Painter, R. K. Lee, A. Scherer, A. Yariv, J. D. O’Brien, P. D. Dapkus, and I. Kim, “Two-dimensional photonic band-gap defect mode laser,” Science **284**, 1819–1821 (1999). [CrossRef] [PubMed]

10. S. Noda, M. Yokoyama, M. Imada, A. Chutinan, and M. Mochizuki, “Polarization mode control of two-dimensional photonic crystal laser by unit cell structure design,” Science **293**, 1123–1125 (2001). [CrossRef] [PubMed]

12. M. Imada, A. Chutinan, S. Noda, and M. Mochizuki, “Multidirectionally distributed feedback photonic crystal lasers,” Phys. Rev. B **65**, 195306 (2002). [CrossRef]

16. H.- G. Park, S.- H. Kim, S.- H. Kwon, Y. -G. Ju, J.- K. Yang, J.- H. Baek, S. -B. Kim, and Y.- H. Lee, “Electrically driven single-cell photonic crystal laser,” Science **305**, 1444–1447 (2004). [CrossRef] [PubMed]

19. P. Bermel, E. Lidorikis, Y. Fink, and J. D. Joannopoulos, “Active materials embedded in photonic crystals and coupled to electromagnetic radiation,” Phys. Rev. B **73**, 165125 (2006). [CrossRef]

22. B. Bakir, C. Seassal, X. Letartre, P. Viktorovitch, M. Zussy, L. Cioccio, and J. Fedeli, “Surface-emitting micro-laser combining two-dimensional photonic crystal membrane and vertical Bragg mirror,” Appl. Phys. Lett. **88**, 081113 (2006). [CrossRef]

23. H. Altug, D. Englund, and J. Vuckovic, “Ultra-fast photonic-crystal nanolasers,” Nat. Phys. **2**, 485–488 (2006). [CrossRef]

26. H. Matsubara, S. Yoshimoto, H. Saito, Y. Jianglin, Y. Tanaka, and S. Noda, “GaN photonic crystal surface-emitting laser at blue-violet wavelengths,” Science **319**, 445–447 (2008). [CrossRef]

31. Y. Kurosaka, S. Iwahashi, Y. Liang, K. Sakai, E. Miyai, W. Kunishi, D. Ohnishi, and S. Noda, “On-chip beam-steering photonic-crystal lasers,” Nat. Photonics **4**, 447–450 (2010). [CrossRef]

*Q*/

*V*

_{mode}microphotonic cavities (

*Q*and

*V*

_{mode}being the corresponding quality factor and the modal volume, respectively) that can be introduced in PhCs simply by inducing local variations in the geometry or the dielectric constant of an otherwise perfectly periodic structure [33, 35

35. S. Noda, M. Fujita, and T. Asano, “Spontaneous-emission control by photonic crystals and nanocavities,” Nat. Photonics **1**, 449–458 (2007). [CrossRef]

10. S. Noda, M. Yokoyama, M. Imada, A. Chutinan, and M. Mochizuki, “Polarization mode control of two-dimensional photonic crystal laser by unit cell structure design,” Science **293**, 1123–1125 (2001). [CrossRef] [PubMed]

26. H. Matsubara, S. Yoshimoto, H. Saito, Y. Jianglin, Y. Tanaka, and S. Noda, “GaN photonic crystal surface-emitting laser at blue-violet wavelengths,” Science **319**, 445–447 (2008). [CrossRef]

36. A. A. Erchak, D. J. Ripin, S. Fan, P. Rakich, J. D. Joannopoulos, E. P. Ippen, G. S. Petrich, and L. A. Kolodziejski, “Enhanced coupling to vertical radiation using a two-dimensional photonic crystal in a semiconductor light-emitting diode,” Appl. Phys. Lett. **78**, 563–565 (2001). [CrossRef]

37. S. Fan and J. D. Joannopoulos, “Analysis of guided resonances in photonic crystal slabs,” Phys. Rev. B **65**, 235112 (2002). [CrossRef]

38. T. Ochiai and K. Sakoda, “Dispersion relation and optical transmittance of a hexagonal photonic crystal slab,” Phys. Rev. B **63**, 125107 (2001). [CrossRef]

*Q*-factor diverges despite lying above the light line). As we show below, in actual finite-size PhC slabs, these

*dark*Fano resonances retain some of the properties of their infinite periodic counterparts and, thus, in the limit of large size PhC slabs, can display arbitrarily large

*Q*values. Since these dark modes typically have photon lifetimes much longer than those of other modes, we expect them to dominate the lasing properties of PCSELs.

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

19. P. Bermel, E. Lidorikis, Y. Fink, and J. D. Joannopoulos, “Active materials embedded in photonic crystals and coupled to electromagnetic radiation,” Phys. Rev. B **73**, 165125 (2006). [CrossRef]

41. A. F. Oskooi, D. Roundy, M. Ibanescu, P. Bermel, J. D. Joannopoulos, and S. G. Johnson, “MEEP: A flexible free-software package for electromagnetic simulations by the FDTD method,” Comput. Phys. Commun. **81**, 687–702 (2010). [CrossRef]

21. H. E. Türeci, A. D. Stone, and B. Collier, “Self-consistent multimode lasing theory for complex or random lasing media,” Phys. Rev. A **74**, 043822 (2006). [CrossRef]

25. H. E. Türeci, L. Ge, S. Rotter, and A. D. Stone, “Strong interactions in multimode random lasers,” Science **320**, 643–646 (2008). [CrossRef] [PubMed]

10. S. Noda, M. Yokoyama, M. Imada, A. Chutinan, and M. Mochizuki, “Polarization mode control of two-dimensional photonic crystal laser by unit cell structure design,” Science **293**, 1123–1125 (2001). [CrossRef] [PubMed]

12. M. Imada, A. Chutinan, S. Noda, and M. Mochizuki, “Multidirectionally distributed feedback photonic crystal lasers,” Phys. Rev. B **65**, 195306 (2002). [CrossRef]

22. B. Bakir, C. Seassal, X. Letartre, P. Viktorovitch, M. Zussy, L. Cioccio, and J. Fedeli, “Surface-emitting micro-laser combining two-dimensional photonic crystal membrane and vertical Bragg mirror,” Appl. Phys. Lett. **88**, 081113 (2006). [CrossRef]

26. H. Matsubara, S. Yoshimoto, H. Saito, Y. Jianglin, Y. Tanaka, and S. Noda, “GaN photonic crystal surface-emitting laser at blue-violet wavelengths,” Science **319**, 445–447 (2008). [CrossRef]

31. Y. Kurosaka, S. Iwahashi, Y. Liang, K. Sakai, E. Miyai, W. Kunishi, D. Ohnishi, and S. Noda, “On-chip beam-steering photonic-crystal lasers,” Nat. Photonics **4**, 447–450 (2010). [CrossRef]

## 2. Methods

### 2.1. General framework

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

19. P. Bermel, E. Lidorikis, Y. Fink, and J. D. Joannopoulos, “Active materials embedded in photonic crystals and coupled to electromagnetic radiation,” Phys. Rev. B **73**, 165125 (2006). [CrossRef]

34. A. S. 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]

*N*(

_{i}*i*= 0, 1, 2,3). Maxwell’s equations for isotropic media are given by

*∂*

**B**(

**r**,

*t*)/

*∂t*= −∇ ×

**E**(

**r**,

*t*) and

*∂*

**D**(

**r**,

*t*)/

*∂t*= ∇ ×

**H**(

**r**,

*t*), where

**B**(

**r**,

*t*) =

*μμ*

_{o}**H**(

**r**,

*t*),

**D**(

**r**,

*t*) =

*ɛɛ*

_{o}**E**(

**r**,

*t*) +

**P**(

**r**,

*t*) and

**P**(

**r**,

*t*) is the dispersive electric polarization density that corresponds to the transitions between two atomic levels,

*N*

_{1}and

*N*

_{2}. The vector

**P**introduces gain in Maxwell’s equation and its time evolution can be shown [47] to follow that of a homogeneously broadened Lorentzian oscillator driven by the coupling between the population inversion and the external electric field. Thus,

**P**obeys the equation of motion where Γ

*stands for the linewidth of the atomic transitions at*

_{m}*ω*= (

_{m}*E*

_{2}–

*E*

_{1})/

*h̄*, and accounts for both the nonradiative energy decay rate, as well as dephasing processes that arise from incoherently driven polarizations.

*E*

_{1}and

*E*

_{2}correspond to the energies of

*N*

_{1}and

*N*

_{2}, respectively.

*σ*is the coupling strength of

_{m}**P**to the external electric field and Δ

*N*(

**r**,

*t*) =

*N*

_{2}(

**r**,

*t*) –

*N*

_{1}(

**r**,

*t*) is the population inversion driving

**P**. Positive inversion is attained when Δ

*N*(

**r**,

*t*) > 0, in which case the medium is amplifying; when Δ

*N*(

**r**,

*t*) < 0, the medium is absorbing. In order to model realistic gain media, only conditions favorable to the former are considered.

*N*(

**r**,

*t*) > 0, the plus sign corresponds to radiation while the minus sign represents excitation.

*R*

_{p}is the external pumping rate that transfers electrons from the ground state to the third excited level, and is proportional to the incident pump power.

*τ*is the nonradiative decay lifetimes from level

_{ij}*i*to

*j*(

*i*>

*j*) so that the energy associated with the decay term

*R*

_{p}.

*N*

_{0}to

*N*

_{3}at a rate

*R*

_{p}. These electrons then decay nonradiately into

*N*

_{2}after a short lifetime

*τ*

_{32}. By enforcing

*τ*

_{21}≫ (

*τ*

_{32},

*τ*

_{10}), a metastable state is formed at

*N*

_{2}favoring a positive population inversion between

*N*

_{2}and

*N*

_{1}(i.e. Δ

*N*> 0), which are separated by energy

*h̄ω*. In this regime, a net decay of electrons to

_{m}*N*

_{1}occurs through stimulated emission and nonradiative relaxation. Lastly, electrons decay nonradiatively and quickly to

*N*

_{0}. Lasing occurs for pumping rates beyond a given threshold

*N*levels may stand for clusters of closely spaced but distinct levels, where the relaxation processes among them are much faster than that with all other levels. More generally, although we will focus on the particular case of four-level gain media, we expect the lasing properties of the PCSELs to be influenced primarily by the EM properties of the passive dielectric structure, rather than the microscopic details of the gain mechanism. Hence, our results should be broadly applicable to any active device describable by semiclassical laser theory. The specific effects of optimizing lasing action will depend on the gain medium. Thus, as we show below, for a four-level gain medium, it leads to an arbitrary reduction of the lasing threshold; whereas in a semiconductor laser, it may lead to lasing thresholds close to the limit set by the transparency condition.

_{i}### 2.2. Finite-difference time-domain simulations of active media

48. K. S. Yee, “Numerical solution of initial boundary value problems involving Maxwell’s equations in isotropic media,” IEEE Trans. Antennas Propag. **14**, 302–307 (1966). [CrossRef]

*x*= Δ

*y*= Δ

*z*, where Δ

*x*, Δ

*y*and Δ

*z*are the space increments in the

*x*,

*y*and

*z*directions) should also be made fine enough so that the PCSEL structures in consideration are well-represented. For numerical stability, Von Neumann analysis places an upper bound on the size of the time step, Δ

*t*≤

*S*Δ

*x*/

*c*, where

*c*is the speed of light. The Courant factor

*S*is typically chosen to be

**E**and

**H**fields at alternate time steps. For simplicity, we show explicitly the set of discretized equations implemented for a one dimensional (1D) setup assuming a non-magnetic and isotropic medium, and denote any functions of space and time as

*F*(

^{n}*i*) =

*F*(

*i*Δ

*z*,

*n*Δ

*t*).

**H**is first updated as Next, we update the polarization density

**P**at

*n*+ 1 from the two previous instances of

**P**, and the previous

*N*and

_{i}**E**according to Eq. (1). Note that components of

**P**reside at the same locations as those of

**E**. We can then use these updated

**H**and

**P**values to retrieve

**E**at

*n*+ 1: Lastly,

*N*at

_{i}*n*+ 1 requires

*N*at

_{i}*n*and both the previous and updated

**E**and

**P**values at

*n*and

*n*+ 1. Since the population densities of the four levels are interdependent in Eq. (2) to (5), they must be solved simultaneously by setting up the following system of equations: where

**Ã**and

**B̃**are tensors that couple the population in the four atomic levels and updated

*N*values at

_{i}*n*+ 1 in Eq. (11) can be computed by inverting

**Ã**. Note that the atomic population density

*N*is a scalar and depends generally on all three components of

_{i}**E**and

**P**. The above cycle is repeated at each time step until steady state is reached, allowing the full temporal development of the laser mode to be tracked.

**P**and

*N*, all physical quantities of the active material are tracked at all points in the computational domain and at all times. The numerical results computed in this

_{i}*ab initio*way are exact apart from discretization of space-time, and allow for nonlinear interactions between the media and the fields. In our simulations, the electric, magnetic and polarization fields are initialized to zero except for background noise while the total electron density is initialized to the ground state level. The computational domain is truncated with Bloch periodic boundary conditions or perfectly matched layers (PML), which are artificial absorbing material designed so that the computational grid’s boundaries are reflectionless in the limit Δ

*z*→ 0 [49

49. S. D. Glauber, “An anisotropic perfectly matched layer absorbing medium for the truncation of FDTD lattices,” IEEE Trans. Antennas Propag. **44**, 1630–1639 (1996). [CrossRef]

*a*, where we checked that the relative differences in frequency and

*Q*values between the 20 pixels per

*a*case and the 40 pixels per

*a*case is less than 5%.

### 2.3. Coupled-mode theory formalism applied to lasing media

50. M. Soljac̆ić, M. Ibanescu, S. G. Johnson, Y. Fink, and J. D. Joannopoulos, “Optimal bistable switching in non-linear photonic crystals,” Phys. Rev. E **66**, 055601 (2002). [CrossRef]

56. J. Bravo-Abad, A. W. Rodriguez, J. D. Joannopoulos, P. T. Rakich, S. G. Johnson, and M. Soljac̆ić, “Efficient low-power terahertz generation via on-chip triply-resonant nonlinear frequency mixing,” Appl. Phys. Lett. **96**, 101110 (2010). [CrossRef]

*P*

_{in}, can be seen as the result of the rate of work performed by the current induced in the system by the active media,

**J**(

**r**,

*t*), against the electric field of the cavity,

**E**(

**r**,

*t*). Noticing that

**J**(

**r**,

*t*) actually comes from the temporal variation of the polarization density, i.e.

**J**(

**r**,

*t*) =

*∂*

**P**(

**r**,

*t*)/

*∂t*[see definition of

**P**in Eq. (1)],

*P*

_{in}can be written as [57]

**E**(

**r**,

*t*) =

**E**

_{0}(

**r**)

*a*(

*t*)exp(−

*iω*) and

_{c}t**P**(

**r**,

*t*) =

**E**

_{0}(

**r**)

*P*(

*t*)exp(−

*iω*), respectively.

_{m}t**E**

_{0}(

**r**) is the normalized cavity mode profile (∫

*d*

^{3}

*rɛ*

_{0}

*n*

^{2}(

**r**)|

**E**

_{0}(

**r**)|

^{2}= 1) and

*a*(

*t*) is the corresponding slowly-varying wave amplitude, normalized so that |

*a*(

*t*)|

^{2}is the energy stored in the resonant mode [43].

*ω*and

_{c}*ω*stand for the resonant frequency of the cavity and the considered atomic transition [see Eq. (1)], respectively.

_{m}*ω*=

_{c}*ω*, using first order perturbation theory in Maxwell’s equations, one can show from energy conservation arguments [43] that the temporal evolution of the electric field amplitude

_{m}*a*(

*t*) is governed by the equation where

*τ*and

_{ex}*τ*are, respectively, the decay rates due to

_{IO}*external*losses (mainly absorption and radiation losses) and due to the decay into the waveguide. The confinement factor

*A*denotes the active part of the structure, accounts for the fact that only the active region drives the temporal evolution of the cavity mode amplitude [as seen in Eq. (13)].

*P*(

*t*). Within the first-order perturbation theory approach we are describing here, from Eq. (1) one can obtain a simple first-order differential equation for the temporal evolution of

*P*(

*t*) by making three main assumptions:

*(i)*we assume that the linewidth of the atomic transition is much smaller than its frequency, i.e. Γ

*≪*

_{m}*ω*;

_{m}*(ii)*we apply the rotating-wave approximation (RWA), i.e. we consider just the terms that oscillate as exp(±

*iω*); and,

_{m}t*(iii)*we apply the slowly-varying envelope approximation (SVEA) [47], i.e. we assume

*d*

^{2}

*P*(

*t*)/

*d*

^{2}

*t*≪

*ω*(

_{m}dP*t*)/

*dt*. For the structures that we will consider below, these approximations give good agreement between CMT and FDTD.

*P*(

*t*) is given by Here 〈Δ

*N*(

*t*)〉 is defined as the population inversion Δ

*N*(

**r**,

*t*) averaged over the mode profile in the gain region of the system Notice also that in Eq. (14), the definitions of

*σ*and

_{m}*ω*are the same as those give in Eq. (1).

_{m}**E**

_{0}(

**r**)

**|**

^{2}over the active region

*A*, following a derivation similar to that given above and after some straightforward algebra, one can obtain the following set of equations governing the time evolution of the average population densities 〈

*N*(

_{i}*t*)〉 (using the same definition of average as in Eq. (15), and with

*i*=0,...,3) where the parameter

*ξ*

_{2}is given by

**E**

_{0}(

**r**), and their corresponding decay rates, one can compute all the relevant physical quantities characterizing lasing action in such structure just by solving the linear system of first-order differential equations given by Eqs. (13), (14), (15)–(19). In particular, once such system of equations have been solved, the total emitted power can be easily computed from

*P*(

_{e}*t*) = 2|

*a*(

*t*)|

^{2}/

*τ*. Although the generalization of the approach described here to the case where more than a single mode is lasing is straightforward, for simplicity, we only consider the single-mode case in this paper.

_{IO}### 2.4. Steady-state ab-initio laser theory

21. H. E. Türeci, A. D. Stone, and B. Collier, “Self-consistent multimode lasing theory for complex or random lasing media,” Phys. Rev. A **74**, 043822 (2006). [CrossRef]

25. H. E. Türeci, L. Ge, S. Rotter, and A. D. Stone, “Strong interactions in multimode random lasers,” Science **320**, 643–646 (2008). [CrossRef] [PubMed]

*dN*/

_{i}*dt*= 0) [21

21. H. E. Türeci, A. D. Stone, and B. Collier, “Self-consistent multimode lasing theory for complex or random lasing media,” Phys. Rev. A **74**, 043822 (2006). [CrossRef]

45. L. Ge, R. Tandy, A. D. Stone, and H. E. Türeci, “Quantitative Verification of Ab Initio Self-Consistent Laser Theory,” Opt. Express **16**, 16895 (2008). [CrossRef] [PubMed]

*ω*is to be determined self-consistently [21

_{L}**74**, 043822 (2006). [CrossRef]

25. H. E. Türeci, L. Ge, S. Rotter, and A. D. Stone, “Strong interactions in multimode random lasers,” Science **320**, 643–646 (2008). [CrossRef] [PubMed]

*ω*differs from the resonance frequency of the passive cavity. In the present work, the gain center will be tuned exactly to the cavity frequency, so we will find that

_{m}*ω*=

_{L}*ω*, but the SALT can also handle the more general detuned case.

_{m}*ɛ*, and an active part that describes the effect of the inverted gain medium. The nonlinearity arises from the spatial hole-burning term

*h*(

*r⃗*); in the single-mode regime this describes the saturation effect of the mode pattern on the gain medium, reducing the inversion and thus changing its own spatial profile. In the multimode regime, this term would also describe mode coupling, which increases the thresholds of the higher-order modes.

**r**) reduces at infinity to a superposition of purely outgoing waves with frequency

*ω*. This non-Hermitian boundary condition is required for a rigorously correct description of an open, steady-state system containing an amplifying medium [21

_{L}**74**, 043822 (2006). [CrossRef]

*ω*, the CF states {

*u*(

_{n}*r⃗*;

*ω*)|

*n*= 1, 2,···} are the discrete solutions to with outgoing boundary conditions. The

*η*’s are complex eigenvalues that can be roughly interpreted as various values of the complex dielectric constant which can produce a resonance pole (purely outgoing solution) at the given real

_{n}*ω*[42]. Equation (30) can be solved by variants of existing numerical techniques; in the present paper, we use the finite-element method (FEM). It can be shown that the CF states are self-orthogonal: Comparing Eq. (30) to Eq. (28), we see that the threshold lasing mode corresponds to a CF state with

*ω*=

*ω*. Hence, the lasing threshold is found by sweeping over a range of frequencies near the gain center

_{L}*ω*, and finding CF eigenvalues obeying Of these, the solution with smallest

_{m}*δ*

_{0}yields the threshold mode, and the matching

*ω*is the lasing frequency

*ω*at threshold. From

_{L}*R*

_{p}-dependence of

*ω*and the mode profile: where

_{L}*I*(

*R*

_{p}) is the mode intensity and

*u*is the CF state corresponding to the mode at threshold. SPA-SALT still includes the non-linear saturation effect on the lasing mode, but by eliminating the need to recompute

_{L}*ω*and the lasing mode profile for each

_{L}*R*

_{p}, the above-threshold SALT calculations is greatly speeded up. It gives excellent results for high-Q laser cavities such as PCSELs, although for low-Q cavities (e.g. random lasers) it is much less accurate [42]. Inserting Eq. (33) into Eq. (28), and using Eq. (31), now yields the simple approximate expression We can obtain the total power output by inserting the above equations back into the physical electric field Eq. (23), and computing the flux of the Poynting vector at infinity. A brief calculation gives, in the physical limit

*L*is the out-of-plane height, and we have assumed that the gain medium is distributed uniformly over the region

_{z}*A*. Aside from small differences in the integrands, this expression agrees with the CMT result derived in Eq. (22). In particular, it states that the mode intensities are approximately linear in

*R*

_{p}for

## 3. Results and discussion

### 3.1. Passive properties—bandstructures

*x*and

*z*directions. In this system, a one-wavelength thick cavity with

*n*= 3.55 (e.g. as in InGaAsP) is enclosed by 25 and 30 bilayers of quarter-wave distributed Bragg reflectors (DBRs) on the top and bottom sides of the structure, respectively. The dielectric mirrors consist of alternate layers of dielectric with

*n*= 3.17 and

*n*= 3.51 (e.g. as in the InP-based lattice matched InP and InGaAlAs, which offers a relatively larger refractive index contrast of Δ

*n*= 0.34 at 1.55

*μm*wavelength; allowing broadband, high reflectivity and low penetration depth DBRs to be attained with fewer layers). Pink shaded regions in Fig. 1(a) represent the continuum of bands guided in the DBRs, while the red line represents the air light line that separates the modes that are propagating in air from those that are evanescent in air. Only transverse magnetic (TM) modes with electric field oriented along the

*z*direction are considered.

*n*= 3.55 center layer (not shown) and that of the multilayer cladding (bottom edge of the lower continuum region). Thus, this mode is guided within the cavity layer via total internal refraction, just as in regular dielectric waveguide slab, with no means of coupling to air. It is the portion of the second mode which lies above the air light line [plotted as a blue line in Fig. 1(a)] that is useful for laser operation. In fact, it is most often desirable to operate at the frequency that corresponds to

*k*= 0 (the so-called Γ point) so that the power is vertically emitted through the surface in the longitudinal (

_{x}*y*) direction. This mode resides in the lowest photonic bandgap of the periodic claddings and, therefore, is trapped within the cavity layer by the high reflectivities (> 96%) of the DBRs. From our calculations, we find that

*Q*, which measures the loss of the VCSEL in the

*y*direction, is 7500 at

*k*= 0 and may generally be increased further by adding more bilayers of the claddings. Thus, VCSEL structures similar to the one described here, resemble a conventional laser cavity in which the eigenmodes are formed in the longitudinal direction due to feedback from the dielectric mirrors and in which the number of the modes increases with the cavity thickness. Notice that the group velocity (

_{x}*v*=

_{g}*dω*/

*dk*) is near zero for small values of

_{x}*k*, which maximizes the wave-matter interaction inside the cavity and, at the same time, enhances the lateral modal confinement.

_{x}*air-bridge*type PhC slabs with 1D corrugation and punctured 2D square lattice of air cylinders respectively. These PhC slabs can support Fano resonances. As mentioned in the introduction, these guided resonances appear in the system when periodic air perturbations, introduced in an otherwise uniform slab, enable the coupling between the guided modes supported by the slab and the external radiation continuum, with the strength of this coupling measured by

*Q*of the slab structures. One major difference between these PhC slabs and VCSEL-like structures is that in the former light confinement occurs in the in-plane periodic directions due to Bragg diffractions, and in the out-of-plane direction due to index guiding. It is this presence of index guiding in the third dimension that limits the photon lifetime at frequencies above the air light line, leading to far-field radiation. Since discrete translational symmetries exists due to in-plane periodicity, the projected band diagrams are plotted with respected to the lateral wave vectors along the irreducible Brillouin zone. We shall briefly examine the geometries of the two slab structures separately, before drawing the similarities between them when operated as band-edge mode lasers.

*a*-thick

*n*= 3.17 (e.g. as in InP) slab with a set of 1D periodic grooves along the

*x*-direction. These grooves are 0.15

*a*deep and 0.1

*a*wide, and extend uniformly in the

*z*direction. Only modes with electric field oriented along

*z*are considered. On the other hand, the PhC slab shown in inset of Fig. 1(c) consist of a 0.3

*a*-thick

*n*= 3.17 (e.g. as in InP) slab punctured with a 2D square lattice of circular air cylinders in the lateral directions, with both depth and radius being equal to 0.25

*a*. In this case, only transverse-electric-like (TE-like) modes, with the electric fields primarily horizontal near the center of the slab, are excited. As in the case of the VCSEL, the modes above the light line at Γ are the most desirable for lasing, since they allow the power to be coupled vertically out of the slab surface. Moreover, in this structure, the zero in-plane group velocity facilitates formation of standing waves, as in any conventional cavity, leading to lateral feedback of the eigenmodes. In fact, in the finite size devices that we will be considering next, Δ

*k*

_{||}≠ 0 so that the dispersion curves near Γ may be well approximated by the second order Taylor expansion, in which case,

*v*becomes directly proportional to the curvature of the bands. Hence, flat dispersion curves having high density of photonic states and low

_{g}*v*are favorable for enhancing light-matter interaction, which is essential for lasing to take place. Note that a VCSEL, on the other hand, has the same direction of periodicity, feedback, and power emission.

_{g}12. M. Imada, A. Chutinan, S. Noda, and M. Mochizuki, “Multidirectionally distributed feedback photonic crystal lasers,” Phys. Rev. B **65**, 195306 (2002). [CrossRef]

*K*= 2

*π*/

*a*is the Bragg grating vector, and we only consider

*q*= 1 to ensure vertical outcoupling. All other higher lying frequency modes result in additional out-of-plane emission directions at oblique angles from the slab surfaces. For the corrugated slab, the phase matching conditions in the reciprocal space also implies that the waves traveling in the +

*x*direction are coupled to those in the −

*x*direction within each unit cell, forming an in-plane feedback mechanism, similar to a 1D cavity. These lateral standing waves are in turn coupled into

*y*because the Bragg condition is also satisfied along the slab normal, enabling perpendicular surface emission. For the slab shown in Fig. 1(c), phase matching at Γ again couples waves in the four equivalent Γ –

*X*directions of a unit cell to the waves emitting in

*z*. Here, the main feedback mechanism is provided separately by waves traveling in the ±

*x*and ±

*y*directions. Further coupling of waves between these orthogonal directions is facilitated by higher order waves traveling in the Γ –

*M*directions (see inset of Fig. 1(c) for the definitions of directions in the reciprocal space of a square lattice). Due to the ease of fabrication resulting from the connected nature of the defect-free lattice, as well as other advantages mentioned at the beginning of this section, PhC slab structures hold great potential as laser devices. The key is its ability to excite a single lateral and longitudinal mode over a large 2D lasing area, as a result of multidimensional distributed feedback mechanism described above. Intuitively, we may treat each unit cell as an individual cavity in-sync with its neighbors, to produce coherent laser oscillations, and desired properties of the lasing mode may be affected simply by tuning the design of each lattice cell. This approach has been experimentally realized to control the polarization of the lasing mode [10

**293**, 1123–1125 (2001). [CrossRef] [PubMed]

*Q*, low threshold laser: the existence of band-edge modes with infinite photon lifetime, i.e. with no means of coupling out of the slab. This phenomenon occurs for the lower band edge in Fig. 1(b) and for singly degenerate modes in 2D periodic PhC slabs, corresponding to the two lowest band-edge modes at Γ in Fig. 1(c). The absence of radiative components at these points in the band diagram is a result of in-phase superpositions of the forward and backward traveling waves, with in-plane electric field vectors adding destructively. This same feature can be explained using the symmetry mismatch existing between the guided modes in the PhC slab and the diffracted radiation field in air [38

38. T. Ochiai and K. Sakoda, “Dispersion relation and optical transmittance of a hexagonal photonic crystal slab,” Phys. Rev. B **63**, 125107 (2001). [CrossRef]

*Q*above the air light line can only be achieved in PhC slabs, this property being absent in VCSELs, or conventional microcavity structures that use high reflectivity mirrors for mode trapping.

*Q*of the two bands above the light line at small values of

*k*, plotted against frequency, in the vicinity of the bandgap for the PCSEL structure shown in Fig. 1(b), with grooves 0.1

_{x}*a*wide and 0.15

*a*deep. We see from the figure that the two band-edge modes differ drastically. The

*Q*of the lower frequency mode diverges rapidly as

*k*→ 0, while that of the next-ordered band remains finite. This is clearly illustrated by the electric field profiles in the unit cell, depicted in the two leftmost panels of Fig. 2(c) for the lower (left) and upper (right) band-edge modes. The unbounded

_{x}*Q*mode, whose radiative electric field component is anti-symmetric about the groove, interferes destructively with itself in the far-field, resulting in no net outcoupling to air. For

*k*away from Γ, this symmetry mismatch is lost, and

_{x}*Q*decreases rapidly but remains large. On the other hand, the second mode is symmetric and vertical emission out of the slab is possible. Note that despite this leakage, most of the electric field is confined within the slab, forming a standing wave pattern due to the lateral feedback mechanism described previously, a signature of Fano resonances. Apart from mode symmetries, the resonances in the slabs are also influenced by the size of the grooves, which may be regarded as periodic dielectric perturbations to an otherwise uniform slab. Results for 1D periodic grooves with depth 0.1

*a*and 0.05

*a*are also shown in Fig. 2(a). Consistent with predictions from the perturbation theory [33], the bandgap decreases with the grooves size while

*Q*increases, approaching the slab waveguide limit of infinity when no grooves are present.

### 3.2. Passive properties—finite structures

*Q*factor, as a function of frequency, for finite slabs with lateral sizes ranging from 20 to 320 periods. These results were obtained from FDTD calculations, with the boundary of the computational domain padded with absorbing boundary conditions (PMLs) to mimic the behavior of a slab in free space. A couple of key observations are in order:

*(i)*The lower band-edge mode of the finite PhC slab no longer possesses an unbounded

*Q*, owing to the fact that an additional loss channel is opened up: energy can now leak from the sides of the slab. This can be observed in the top right panel of Fig. 2(c) for a 20

*a*long PhC slab. These lateral losses dominate in the lower band-edge mode. The bottom right panel of Fig. 2(c) shows the symmetric mode, where both vertical and lateral power emission appears equally dominant. It is thus no surprise that the net

*Q*of the lower band-edge mode remains higher than that of the symmetric mode [see Fig. 2(b)].

*(ii)*The

*Q*of the lasing structure increases with the number of periods, so the lasing threshold correspondingly decreases. We shall quantify the losses in Fig. 3, as functions of the number of periods.

*(iii)*The resonant frequencies of the upper band-edge mode are different in the finite and infinite slabs, due to the presence of lateral losses in the former. In the finite system, increasing

*ω*leads to a corresponding increase in

*k*and

_{x}*v*, and hence a decrease in

_{g}*Q*. In the infinite system, there are no lateral losses, so

*Q*increases with frequency near the band-edge. For the lower band-edge, mode symmetry considerations ensure that

*Q*remains a maximum for both infinite and finite slabs.

*Q*values of the corrugated slab in order to understand how the lateral size of the device,

*L*, affects the outcoupling of Fano resonances. Figure 3 compiles the total

_{x}*Q*(

*Q*

_{tot}) of the two band-edge modes presented in Fig. 2(b) for PCSEL structures having 0.05

*a*deep grooves, with

*L*ranging from 20

_{x}*a*to 320

*a*. In order to operate the device at typical optical communication wavelength (∼ 1.55

*μ*m), we set

*a*= 675 nm here and in subsequent results. Since a larger PhC slab provides longer confinement time,

*Q*increases with

*L*for both symmetric and anti-symmetric modes. The anti-symmetric mode has higher

_{x}*Q*, due to its reduced vertical emission, as already observed in Fig. 2(c). For

*L*> 100

_{x}*μ*m, the total

*Q*of both modes tends towards that of their infinite counterpart [see Fig. 2(b)]:

*Q*

_{sym}saturates at 1964, whereas

*Q*

_{anti-sym}is unbounded. Therefore, the anti-symmetric mode holds great potential for low-threshold laser operation. Using approximate analytic relationships that govern

*Q*’s dependence on

*L*(unique for each mode), curves fitted to the calculated data are also plotted in Fig. 3. We shall specify these relationships in the next paragraph.

_{x}### 3.3. Lasing—infinite periodic structures

*τ*

_{10}=

*τ*

_{32}= 5 ×10

^{−14}s and

*τ*

_{21}= 5 × 10

^{−12}s (so a metastable state can form at

*N*

_{2}). For the coupling constant, we take

*σ*= 1 ×10

_{m}^{−4}C

^{2}/kg (this value was obtained assuming that the Purcell effect is negligible); for the total electron density, we take

*N*

_{tot}= 5×10

^{23}m

^{−3}. These values are realistic, and similar to those used in Ref. [9

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

*ω*, in order to select the mode that we wish to lase; in the FDTD calculations, this frequency is set to the frequency of the corresponding passive mode. The gain linewidth Γ

_{m}*is taken to be 0.002 (2*

_{m}*πc*/

*a*), which is sufficiently narrow to avoid exciting neighboring modes.

*infinite*slab. The computational domain is similar to the one shown in Fig. 1(b), with periodic boundary conditions along the left and right boundaries and PML absorbers along the top and bottom boundaries. Figure 4 shows the resulting plot of output power versus

*R*

_{p}. Three different structures, with groove depths of 0.05

*a*, 0.1

*a*, and 0.15

*a*, are simulated; the groove width is kept at 0.1

*a*, and slab thickness at 0.3

*a*. The filled circles in this plot are the results of FDTD calculations (time-stepping until steady-state laser operation was achieved); the solid lines are the CMT predictions, with parameters fitted from separate FDTD calculations of the passive structure’s Fano resonance frequency, electric field mode profile, and

*Q*.

*ω*is situated at the resonance frequency of the symmetric mode, as determined by the passive-structure FDTD calculations presented earlier. (Since this is the infinite structure, the anti-symmetric mode does not allow power emission.) For groove depths [0.05

_{m}*a*,0.1

*a*, 0.15

*a*], we have

*ω*= [0.435, 0.444, 0.456] ×(2

_{m}*πc*/

*a*) and

*Q*values [1964, 451, 230]; see Fig. 2(a). As expected, the laser threshold is inversely proportional to

*Q*; physically speaking, higher input pump rates are needed to overcome larger losses. Moreover, the three structures exhibit very similar rates of growth of output power,

*dP*/

*dR*

_{p}. As we shall see, this is not true for finite structures.

*R*

_{p}above the lasing threshold, and the FDTD results are very close to linear. The match remains excellent for

*R*

_{p}as much as an order of magnitude above the lasing threshold. This shows that the CMT model that we have developed greatly complements the FDTD approach. The CMT is particularly useful for

*R*

_{p}near threshold, where FDTD computations are very time-consuming due to the temporal turn-on delay before lasing action begins. For larger

*R*

_{p}, the results begin to deviate; the influence of the gain media on the fields can no longer be taken to be linear, so second order corrections to CMT are required and the lasing modes are no longer accurately described by the modes of the linear (passive) cavity.

### 3.4. Lasing—finite structures

*R*

_{p}versus total power output obtained by FDTD and CMT, for finite slabs. (In the FDTD calculations, the PML absorber is now placed along all four boundaries of the computational domain.) Three different slab widths are used:

*L*= 20

_{x}*a*, 40

*a*, and 80

*a*. We fix the groove depth at 0.05

*a*, with all other parameters kept the same as in Fig. 4.

*a*, which is relatively low; in particular, the groove depth of 0.05

*a*corresponds to one pixel. The low resolution is due to computational limitations, especially near threshold, where very many time steps are required to bring the laser to its steady state. The CMT parameters were fitted from passive FDTD calculations performed with the same spatial resolution, and the fact that the CMT results agree well with FDTD indicates that the two methods face consistent discretization errors, as expected.

*a*within the dielectric structure. Convergence was tested by halving the mesh size, whereupon ∼1% deviation in the computed CF eigenvalues was observed.

*η*(

_{n}*ω*)} of the finite PCSEL structure, for several different values of the total slab length

*L*. Recall from Section 2.4 that the CF eigenvalues are the discrete complex contributions to the dielectric function from the gain medium, required to produce a resonance at a given real frequency

_{x}*ω*. For each

*ω*, we find that one particular CF state has

*η*(

_{n}*ω*) lying significantly closer to the real axis (i.e. requiring less amplification) than all the others. The spatial structure of this CF state, shown in the inset of Fig. 6 for

*L*= 20

_{x}*a*, corresponds closely to the anti-symmetric Fano mode of the passive structure (see Fig. 2). The complex CF eigenvalues depend on

*ω*, and increasing

*ω*causes them to drift leftwards in the complex

*η*plane (reflecting increasing mode confinement). The CF eigenvalue spectra in Fig. 6 are plotted at the optimal threshold lasing frequency

*ω*of each structure. In the SALT calculations, we did not choose

_{L}*ω*using the FDTD calculations of the passive structure, as we did for the FDTD and CMT calculations; this would be inappropriate, as the SALT calculation uses a different and finer grid (this discretization mismatch is further exacerbated by the fact that Γ

_{m}*is chosen to be small). Instead, the thresholds were found using the self-consistent procedure described in Section 2.4, and the optimal choice of gain center*

_{m}*ω*is the one that minimizes

_{m}*η*(

_{n}*ω*) purely imaginary, so that

_{L}*ω*≃

_{m}*ω*.

_{L}*L*on the optimal lasing frequency

_{x}*ω*, the lasing threshold

_{L}*dP*/

*dR*

_{p}. The increase of

*ω*with

_{L}*L*agrees with the FDTD results shown in Fig. 2, while the numerical values of

_{x}*ω*differ from the FDTD predictions by ∼ 4%, a very acceptable deviation considering the difference in resolution between the two calculations. Likewise,

_{L}*L*. The power slope, given by Eq. (36), was found to increase approximately linearly with

_{x}*L*. The power output calculated by the SALT is shown by the dashed curves in Fig. 5. Apart from the aforementioned 9% difference in

_{x}*R*

_{p}.

*ω*and Γ

_{m}*. For experimental systems, in which the gain parameters are not fully controllable, such calculations are useful for estimating the range of single-mode operation.*

_{m}### 3.5. 3D slab lasers

*a*, while the radius and depth of the air cylinders are 0.15

*a*. This design operates based on the same principles as the simpler 1D periodic grooves design, so that the physical concepts explored previously may be equally applied in this case. In Fig. 8, the magnetic and electric field profiles of TE-like excitations are provided for the four modes at Γ, two being non-degenerate (two lowest frequency modes) and the remaining pair is degenerate. Analogous to the anti-symmetric modes that exist in corrugated slabs, the non-degenerate modes of such infinite periodic slabs have infinite photon lifetime, which again may be attributed to mode symmetry mismatch with the radiative continuum [37

37. S. Fan and J. D. Joannopoulos, “Analysis of guided resonances in photonic crystal slabs,” Phys. Rev. B **65**, 235112 (2002). [CrossRef]

*E*components are of opposite directions relative to the air cylinder. The same holds for

_{x}*E*. Hence, no coupling to air is observed. On the other hand, vertical radiation of the electric field occurs for the degenerate modes in Fig. 8(c) indicating low finite

_{y}*Q*values. We apply the CMT approach to calculate the power output from a unit cell for the degenerate mode at three radii of the air cylinders shown in Fig. 8(d). As in the corrugated slab,

*Q*increases for smaller air cylinders leading to lower threshold pump rate. Practical considerations with regards to size of the structural periodic perturbations include the ease of fabrication as dimensions scale down, and also the need for close proximity to the gain layers for enhancement of the band-edge effects.

*a*-thick PhC slabs with air cylinders 0.15

*a*deep and diameter 0.3

*a*. Three sizes of the PhC region are studied: 15

*a*× 15

*a*, 25

*a*× 25

*a*and 35

*a*× 35

*a*. In order to excite PhC states and to model realistic conditions in similar PhC lasers operating with an optically or electrically pumped central area, the 2D PhC region has to be surrounded by un-pumped regions. This is achieved in our simulations by truncating the finite size PhC in air and extending the uniform dielectric slab into the PML, which surround the whole computational domain. The lasing mode considered is that of the first singly-degenerate mode shown in Fig. 8(a) and may be compared to the field profiles presented in Fig. 9 for the slab with 15

*a*× 15

*a*PhC region. The slope [as defined in Eq. (22)] and threshold pump rate improves for the larger PhC, consistent with what we would have expected, while the PCSEL remains single-mode. The primary losses for the small-sized PhC considered here is through the lateral leakage into the absorbing boundaries. Further note that magnitudes of the output power, and hence slope, are significantly higher than those found for the corrugated slabs in Fig. 5. This can be understood as a consequence of an additional dimension present in the current calculations. We confirm that the singly degenerate modes for these finite slabs have high

*Q*s relative to the degenerate modes and hence, should be the ones within the first set of frequencies at Γ that, in practice, are selected for lasing. Such mode selection (singly-degenerate) should be possible owing to the distinct frequencies and field profiles of the four modes at Γ, and its lower pumping rate requirements. We also note here that upon successful coupling to the desired laser modes,

*Q*may further be increased by enhancing the confinement geometrically:

*(i)*in the lateral directions by employing PhC heterostructures, and

*(ii)*in the out-of-plane direction via the addition of DBRs on top or below the PhC slab [58

58. T. Lu, C. Kao, H. Kuo, G. Huang, and S. Wang, “CW lasing of current injection blue GaN-based vertical cavity surface emitting laser,” Appl. Phys. Lett. **92**, 141102 (2008). [CrossRef]

## 4. Conclusions

*Q*factors displayed by these dark states lead to a significant reduction of the corresponding lasing threshold with respect to conventional VCSELs. In addition, our calculations suggest that, for input pump rates close to the threshold, PCSEL structures emit most of their lasing power in the plane of periodicity rather than in the vertical direction. However, notice that this lasing power can be directed into the out-of-plane direction simply by perturbing the symmetry of dark Fano resonances [10

**293**, 1123–1125 (2001). [CrossRef] [PubMed]

**85**, 70–73 (2000). [CrossRef] [PubMed]

## Acknowledgments

## References and links

1. | T. H. Maiman, “Stimulated optical radiation in ruby,” Nature |

2. | E. Yablonovitch, “Inhibited spontaneous emission in solid-state physics and electronics,” Phys. Rev. Lett. |

3. | S. Haroche and D. Kleppner, “Cavity quantum electrodynamics,” Phys. Today |

4. | Y. Yamamoto (ed.), |

5. | H. Yokoyama, “Physics and device applications of optical microcavities,” Science |

6. | N. M. Lawandy, R. M. Balachandran, A. S. L. Gomes, and E. Sauvain, “Laser action in strongly scattering media,” Nature |

7. | D. S. Wiersma and A. Lagendijk, “Light diffusion with gain and random lasers,” Phys. Rev. E |

8. | O. Painter, R. K. Lee, A. Scherer, A. Yariv, J. D. O’Brien, P. D. Dapkus, and I. Kim, “Two-dimensional photonic band-gap defect mode laser,” Science |

9. | X. Jiang and C. M. Soukoulis, “Time dependent theory for random lasers,” Phys. Rev. Lett. |

10. | S. Noda, M. Yokoyama, M. Imada, A. Chutinan, and M. Mochizuki, “Polarization mode control of two-dimensional photonic crystal laser by unit cell structure design,” Science |

11. | H. Cao, Y. Ling, J. Y. Xu, C. Q. Cao, and P. Kumar, “Photon statistics of random lasers with resonant feedback,” Phys. Rev. Lett. |

12. | M. Imada, A. Chutinan, S. Noda, and M. Mochizuki, “Multidirectionally distributed feedback photonic crystal lasers,” Phys. Rev. B |

13. | J. C. Johnson, H.- J. Choi, K. P. Knutsen, R. D. Schaller, P. Yang, and R. J. Saykally, “Single gallium nitride nanowire lasers,” Nat. Mater. |

14. | D. J. Bergman and M. L. Stockman, “Surface Plasmon Amplification by Stimulated Emission of Radiation: Quantum Generation of Coherent Surface Plasmons in Nanosystems,” Phys. Rev. Lett. |

15. | X. Duan, Y. Huang, R. Agarwal, and C. M. Lieber, “Single-nanowire electrically driven lasers,” Nature |

16. | H.- G. Park, S.- H. Kim, S.- H. Kwon, Y. -G. Ju, J.- K. Yang, J.- H. Baek, S. -B. Kim, and Y.- H. Lee, “Electrically driven single-cell photonic crystal laser,” Science |

17. | T. Baba, D. Sano, K. Nozaki, K. Inoshita, Y. Kuroki, and F. Koyama, “Observation of fast spontaneous emission decay in GaInAsP photonic crystal point defect nanocavity at room temperature,” Appl. Phys. Lett. |

18. | H. Y. Ryu, M. Notomi, E. Kuramochi, and T. Segawa, “Large spontaneous emission factor (> 0.1) in the photonic crystal monopole-mode laser,” Appl. Phys. Lett. |

19. | P. Bermel, E. Lidorikis, Y. Fink, and J. D. Joannopoulos, “Active materials embedded in photonic crystals and coupled to electromagnetic radiation,” Phys. Rev. B |

20. | C. Gmachl, F. Capasso, E.E. Narimanov, J.U. Nöckel, A. D. Stone, J. Faist, D. Sivco, and A. Cho, “High power directional emission from lasers with chaotic resonators,” Science |

21. | H. E. Türeci, A. D. Stone, and B. Collier, “Self-consistent multimode lasing theory for complex or random lasing media,” Phys. Rev. A |

22. | B. Bakir, C. Seassal, X. Letartre, P. Viktorovitch, M. Zussy, L. Cioccio, and J. Fedeli, “Surface-emitting micro-laser combining two-dimensional photonic crystal membrane and vertical Bragg mirror,” Appl. Phys. Lett. |

23. | H. Altug, D. Englund, and J. Vuckovic, “Ultra-fast photonic-crystal nanolasers,” Nat. Phys. |

24. | M. T. Hill, Y. -S. Oei, B. Smalbrugge, Y. Zhu, T. de Vries, P. J. van Veldhoven, F. W. M. van Otten, T. J. Eijkemans, J. P. Turkiewicz, H. de Waardt, E. J. Geluk, S. -H. Kwon, Y.- H. Lee, R. Notzel, and M. K. Smit, “Small-divergence semiconductor lasers by plasmonic collimation,” Nat. Photonics |

25. | H. E. Türeci, L. Ge, S. Rotter, and A. D. Stone, “Strong interactions in multimode random lasers,” Science |

26. | H. Matsubara, S. Yoshimoto, H. Saito, Y. Jianglin, Y. Tanaka, and S. Noda, “GaN photonic crystal surface-emitting laser at blue-violet wavelengths,” Science |

27. | S. Gottardo, R. Sapienza, P. D. García, A. Blanco, D. S. Wiersma, and C. López, “Resonance-driven random lasing,” Nat. Photonics |

28. | N. I. Zheludev, S. L. Prosvirnin, N. Papasimakis, and V. A. Fedotov, “Lasing spaser,” Nat. Photonics |

29. | M. A. Noginov, G. Zhu, A. M. Belgrave, R. Bakker, V. M. Shalaev, E. E. Narimanov, S. Stout, E. Herz, T. Suteewong, and U. Wiesner, “Demonstration of a spaser-based nanolaser,” Nature |

30. | R. F. Oulton, V. J. Sorger, T. Zentgraf, R.-M. Ma, C. Gladden, L. Dai, G. Bartal, and X. Zhang, “Plasmon lasers at deep subwavelength scale,” Nature |

31. | Y. Kurosaka, S. Iwahashi, Y. Liang, K. Sakai, E. Miyai, W. Kunishi, D. Ohnishi, and S. Noda, “On-chip beam-steering photonic-crystal lasers,” Nat. Photonics |

32. | M. P. Nezhad, A. Simic, O. Bondarenko, B. Slutsky, A. Mizrahi, L. Feng, V. Lomakin, and Y. Fainman, “Room-temperature subwavelength metallo-dielectric lasers,” Nat. Photonics |

33. | J. D. Joannopoulos, S. G. Johnson, R. D. Meade, and J. N. Winn, |

34. | A. S. Nagra and R. A. York, “FDTD analysis of wave propagation in nonlinear absorbing and gain media,” IEEE Trans. Antennas Propag. |

35. | S. Noda, M. Fujita, and T. Asano, “Spontaneous-emission control by photonic crystals and nanocavities,” Nat. Photonics |

36. | A. A. Erchak, D. J. Ripin, S. Fan, P. Rakich, J. D. Joannopoulos, E. P. Ippen, G. S. Petrich, and L. A. Kolodziejski, “Enhanced coupling to vertical radiation using a two-dimensional photonic crystal in a semiconductor light-emitting diode,” Appl. Phys. Lett. |

37. | S. Fan and J. D. Joannopoulos, “Analysis of guided resonances in photonic crystal slabs,” Phys. Rev. B |

38. | T. Ochiai and K. Sakoda, “Dispersion relation and optical transmittance of a hexagonal photonic crystal slab,” Phys. Rev. B |

39. | A. Yariv and P. Yeh, |

40. | A. Taflove and S. C. Hagness, |

41. | A. F. Oskooi, D. Roundy, M. Ibanescu, P. Bermel, J. D. Joannopoulos, and S. G. Johnson, “MEEP: A flexible free-software package for electromagnetic simulations by the FDTD method,” Comput. Phys. Commun. |

42. | Li Ge, D. Y. Chong, and A. D. Stone, “Steady-state Ab Initio Laser Theory: Generalizations and Analytic Results,” arxiv:1008.0628. |

43. | H. A. Haus, |

44. | L. Ge, “Steady-state Ab Initio Laser Theory and its Applications in Random and Complex Media,” Yale PhD thesis, 2010; also A. Cerjan, Y.-D. Chong, L. Ge and A. D. Stone, unpublished. |

45. | L. Ge, R. Tandy, A. D. Stone, and H. E. Türeci, “Quantitative Verification of Ab Initio Self-Consistent Laser Theory,” Opt. Express |

46. | H. E. Li and K. Iga, |

47. | A. E. Siegman, |

48. | K. S. Yee, “Numerical solution of initial boundary value problems involving Maxwell’s equations in isotropic media,” IEEE Trans. Antennas Propag. |

49. | S. D. Glauber, “An anisotropic perfectly matched layer absorbing medium for the truncation of FDTD lattices,” IEEE Trans. Antennas Propag. |

50. | M. Soljac̆ić, M. Ibanescu, S. G. Johnson, Y. Fink, and J. D. Joannopoulos, “Optimal bistable switching in non-linear photonic crystals,” Phys. Rev. E |

51. | W. Suh, Z. Wang, and S. Fan, “Temporal coupled-mode theory and the presence of non-orthogonal modes in lossless multimode cavities,” IEEE J. Quantum Electron. |

52. | J. Bravo-Abad, A. Rodriguez, P. Bermel, S. G. Johnson, J. D. Joannopoulos, and M. Soljac̆ić, “Enhanced non-linear optics in photonic-crystal microcavities,” Opt. Express |

53. | R. E. Hamam, M. Ibanescu, E. J. Reed, P. Bermel, S. G. Johnson, E. Ippen, J. D. Joannopoulos, and M. Soljac̆ić, “Purcell effect in nonlinear photonic structures: A coupled mode theory analysis,” Opt. Express |

54. | J. Bravo-Abad, E. P. Ippen, and M. Soljac̆ić, “Ultrafast photodetection in an all-silicon chip enabled by two-photon absorption,” Appl. Phys. Lett. |

55. | H. Hashemi, A. W. Rodriguez, J. D. Joannopoulos, M. Soljac̆ić, and S. G. Johnson, “Nonlinear harmonic generation and devices in doubly-resonant Kerr cavities,” Phys. Rev. A |

56. | J. Bravo-Abad, A. W. Rodriguez, J. D. Joannopoulos, P. T. Rakich, S. G. Johnson, and M. Soljac̆ić, “Efficient low-power terahertz generation via on-chip triply-resonant nonlinear frequency mixing,” Appl. Phys. Lett. |

57. | J. A. Kong, |

58. | T. Lu, C. Kao, H. Kuo, G. Huang, and S. Wang, “CW lasing of current injection blue GaN-based vertical cavity surface emitting laser,” Appl. Phys. Lett. |

**OCIS Codes**

(250.7270) Optoelectronics : Vertical emitting lasers

(270.3430) Quantum optics : Laser theory

(230.5298) Optical devices : Photonic crystals

**ToC Category:**

Quantum Optics

**History**

Original Manuscript: October 20, 2010

Revised Manuscript: December 8, 2010

Manuscript Accepted: December 10, 2010

Published: January 13, 2011

**Citation**

Song-Liang Chua, Yidong Chong, A. Douglas Stone, Marin Soljačić, and Jorge Bravo-Abad, "Low-threshold lasing action in photonic crystal slabs enabled by Fano resonances," Opt. Express **19**, 1539-1562 (2011)

http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-19-2-1539

Sort: Year | Journal | Reset

### References

- T. H. Maiman, “Stimulated optical radiation in ruby,” Nature 187, 493–494 (1960). [CrossRef]
- E. Yablonovitch, “Inhibited spontaneous emission in solid-state physics and electronics,” Phys. Rev. Lett. 58, 2059–2062 (1987). [CrossRef] [PubMed]
- S. Haroche, D. Kleppner, “Cavity quantum electrodynamics,” Phys. Today 42, 24–30 (1989). [CrossRef]
- Y. Yamamoto (ed.), Coherence, Amplification, and Quantum effects in Semiconductor Lasers (John Wiley & Sons, 1991), and references therein.
- H. Yokoyama, “Physics and device applications of optical microcavities,” Science 256, 66–70 (1992). [CrossRef] [PubMed]
- N. M. Lawandy, R. M. Balachandran, A. S. L. Gomes, E. Sauvain, “Laser action in strongly scattering media,” Nature 368, 436–438 (1994). [CrossRef]
- D. S. Wiersma, A. Lagendijk, “Light diffusion with gain and random lasers,” Phys. Rev. E 54, 4256–4265 (1996). [CrossRef]
- O. Painter, R. K. Lee, A. Scherer, A. Yariv, J. D. O’Brien, P. D. Dapkus, I. Kim, “Two-dimensional photonic band-gap defect mode laser,” Science 284, 1819–1821 (1999). [CrossRef] [PubMed]
- X. Jiang, C. M. Soukoulis, “Time dependent theory for random lasers,” Phys. Rev. Lett. 85, 70–73 (2000). [CrossRef] [PubMed]
- S. Noda, M. Yokoyama, M. Imada, A. Chutinan, M. Mochizuki, “Polarization mode control of two-dimensional photonic crystal laser by unit cell structure design,” Science 293, 1123–1125 (2001). [CrossRef] [PubMed]
- H. Cao, Y. Ling, J. Y. Xu, C. Q. Cao, P. Kumar, “Photon statistics of random lasers with resonant feedback,” Phys. Rev. Lett. 86, 4524–4527 (2001). [CrossRef] [PubMed]
- M. Imada, A. Chutinan, S. Noda, M. Mochizuki, “Multidirectionally distributed feedback photonic crystal lasers,” Phys. Rev. B 65, 195306 (2002). [CrossRef]
- J. C. Johnson, H.- J. Choi, K. P. Knutsen, R. D. Schaller, P. Yang, R. J. Saykally, “Single gallium nitride nanowire lasers,” Nat. Mater. 1, 106–110 (2002). [CrossRef]
- D. J. Bergman, M. L. Stockman, “Surface Plasmon Amplification by Stimulated Emission of Radiation: Quantum Generation of Coherent Surface Plasmons in Nanosystems,” Phys. Rev. Lett. 90, 027402 (2003). [CrossRef] [PubMed]
- X. Duan, Y. Huang, R. Agarwal, C. M. Lieber, “Single-nanowire electrically driven lasers,” Nature 421, 241–245 (2003). [CrossRef] [PubMed]
- H.- G. Park, S.- H. Kim, S.- H. Kwon, Y. -G. Ju, J.- K. Yang, J.- H. Baek, S. -B. Kim, Y.- H. Lee, “Electrically driven single-cell photonic crystal laser,” Science 305, 1444–1447 (2004). [CrossRef] [PubMed]
- T. Baba, D. Sano, K. Nozaki, K. Inoshita, Y. Kuroki, F. Koyama, “Observation of fast spontaneous emission decay in GaInAsP photonic crystal point defect nanocavity at room temperature,” Appl. Phys. Lett. 85, 3889–3891 (2004). [CrossRef]
- H. Y. Ryu, M. Notomi, E. Kuramochi, T. Segawa, “Large spontaneous emission factor (> 0.1) in the photonic crystal monopole-mode laser,” Appl. Phys. Lett. 84, 1067–1069 (2004). [CrossRef]
- P. Bermel, E. Lidorikis, Y. Fink, J. D. Joannopoulos, “Active materials embedded in photonic crystals and coupled to electromagnetic radiation,” Phys. Rev. B 73, 165125 (2006). [CrossRef]
- C. Gmachl, F. Capasso, E.E. Narimanov, J.U. Nöckel, A. D. Stone, J. Faist, D. Sivco, A. Cho, “High power directional emission from lasers with chaotic resonators,” Science 280, 1556–64 (1998). [CrossRef] [PubMed]
- H. E. Türeci, A. D. Stone, B. Collier, “Self-consistent multimode lasing theory for complex or random lasing media,” Phys. Rev. A 74, 043822 (2006). [CrossRef]
- B. Bakir, C. Seassal, X. Letartre, P. Viktorovitch, M. Zussy, L. Cioccio, J. Fedeli, “Surface-emitting micro-laser combining two-dimensional photonic crystal membrane and vertical Bragg mirror,” Appl. Phys. Lett. 88, 081113 (2006). [CrossRef]
- H. Altug, D. Englund, J. Vuckovic, “Ultra-fast photonic-crystal nanolasers,” Nat. Phys. 2, 485–488 (2006). [CrossRef]
- M. T. Hill, Y. -S. Oei, B. Smalbrugge, Y. Zhu, T. de Vries, P. J. van Veldhoven, F. W. M. van Otten, T. J. Eijkemans, J. P. Turkiewicz, H. de Waardt, E. J. Geluk, S. -H. Kwon, Y.- H. Lee, R. Notzel, M. K. Smit, “Small-divergence semiconductor lasers by plasmonic collimation,” Nat. Photonics 1, 589–594 (2007). [CrossRef]
- H. E. Türeci, L. Ge, S. Rotter, A. D. Stone, “Strong interactions in multimode random lasers,” Science 320, 643–646 (2008). [CrossRef] [PubMed]
- H. Matsubara, S. Yoshimoto, H. Saito, Y. Jianglin, Y. Tanaka, S. Noda, “GaN photonic crystal surface-emitting laser at blue-violet wavelengths,” Science 319, 445–447 (2008). [CrossRef]
- S. Gottardo, R. Sapienza, P. D. García, A. Blanco, D. S. Wiersma, C. López, “Resonance-driven random lasing,” Nat. Photonics 2, 429–432 (2008). [CrossRef]
- N. I. Zheludev, S. L. Prosvirnin, N. Papasimakis, V. A. Fedotov, “Lasing spaser,” Nat. Photonics 2, 351–354 (2008). [CrossRef]
- M. A. Noginov, G. Zhu, A. M. Belgrave, R. Bakker, V. M. Shalaev, E. E. Narimanov, S. Stout, E. Herz, T. Suteewong, U. Wiesner, “Demonstration of a spaser-based nanolaser,” Nature 460, 1110–1112 (2009). [CrossRef] [PubMed]
- R. F. Oulton, V. J. Sorger, T. Zentgraf, R.-M. Ma, C. Gladden, L. Dai, G. Bartal, X. Zhang, “Plasmon lasers at deep subwavelength scale,” Nature 461, 629 (2009). [CrossRef] [PubMed]
- Y. Kurosaka, S. Iwahashi, Y. Liang, K. Sakai, E. Miyai, W. Kunishi, D. Ohnishi, S. Noda, “On-chip beam-steering photonic-crystal lasers,” Nat. Photonics 4, 447–450 (2010). [CrossRef]
- M. P. Nezhad, A. Simic, O. Bondarenko, B. Slutsky, A. Mizrahi, L. Feng, V. Lomakin, Y. Fainman, “Room-temperature subwavelength metallo-dielectric lasers,” Nat. Photonics 4, 395–399 (2010). [CrossRef]
- J. D. Joannopoulos, S. G. Johnson, R. D. Meade, J. N. Winn, Photonic Crystals: Molding the Flow of Light (Princeton Univ. Press, 2008).
- A. S. Nagra, R. A. York, “FDTD analysis of wave propagation in nonlinear absorbing and gain media,” IEEE Trans. Antennas Propag. 46, 334–340 (1998). [CrossRef]
- S. Noda, M. Fujita, T. Asano, “Spontaneous-emission control by photonic crystals and nanocavities,” Nat. Photonics 1, 449–458 (2007). [CrossRef]
- A. A. Erchak, D. J. Ripin, S. Fan, P. Rakich, J. D. Joannopoulos, E. P. Ippen, G. S. Petrich, L. A. Kolodziejski, “Enhanced coupling to vertical radiation using a two-dimensional photonic crystal in a semiconductor light-emitting diode,” Appl. Phys. Lett. 78, 563–565 (2001). [CrossRef]
- S. Fan, J. D. Joannopoulos, “Analysis of guided resonances in photonic crystal slabs,” Phys. Rev. B 65, 235112 (2002). [CrossRef]
- T. Ochiai, K. Sakoda, “Dispersion relation and optical transmittance of a hexagonal photonic crystal slab,” Phys. Rev. B 63, 125107 (2001). [CrossRef]
- A. Yariv, P. Yeh, Photonics: Optical Electronics in Modern Communications (Oxford University Press, New York, NY, 2007).
- A. Taflove, S. C. Hagness, Computational Electrodynamics: The Finite-Difference Time-Domain Method (Artech House, 2005).
- A. F. Oskooi, D. Roundy, M. Ibanescu, P. Bermel, J. D. Joannopoulos, S. G. Johnson, “MEEP: A flexible free-software package for electromagnetic simulations by the FDTD method,” Comput. Phys. Commun. 81, 687–702 (2010). [CrossRef]
- Li Ge, D. Y. Chong, A. D. Stone, “Steady-state Ab Initio Laser Theory: Generalizations and Analytic Results,” arxiv:1008.0628.
- H. A. Haus, Waves and Fields in Optoelectronics (Prentice-Hall, 1984).
- L. Ge, “Steady-state Ab Initio Laser Theory and its Applications in Random and Complex Media,” Yale PhD thesis, 2010; also A. Cerjan, Y.-D. Chong, L. Ge and A. D. Stone, unpublished.
- L. Ge, R. Tandy, A. D. Stone, H. E. Türeci, “Quantitative Verification of Ab Initio Self-Consistent Laser Theory,” Opt. Express 16, 16895 (2008). [CrossRef] [PubMed]
- H. E. Li, K. Iga, Vertical-cavity surface-emitting laser devices (Springer, 2003).
- A. E. Siegman, Lasers (Univ. Science Books, 1986).
- K. S. Yee, “Numerical solution of initial boundary value problems involving Maxwell’s equations in isotropic media,” IEEE Trans. Antennas Propag. 14, 302–307 (1966). [CrossRef]
- S. D. Glauber, “An anisotropic perfectly matched layer absorbing medium for the truncation of FDTD lattices,” IEEE Trans. Antennas Propag. 44, 1630–1639 (1996). [CrossRef]
- M. Soljac̆ić, M. Ibanescu, S. G. Johnson, Y. Fink, J. D. Joannopoulos, “Optimal bistable switching in non-linear photonic crystals,” Phys. Rev. E 66, 055601 (2002). [CrossRef]
- W. Suh, Z. Wang, S. Fan, “Temporal coupled-mode theory and the presence of non-orthogonal modes in lossless multimode cavities,” IEEE J. Quantum Electron. 40, 1511–1518 (2004). [CrossRef]
- J. Bravo-Abad, A. Rodriguez, P. Bermel, S. G. Johnson, J. D. Joannopoulos, M. Soljac̆ić, “Enhanced non-linear optics in photonic-crystal microcavities,” Opt. Express 15, 16161–16176 (2007). [CrossRef] [PubMed]
- R. E. Hamam, M. Ibanescu, E. J. Reed, P. Bermel, S. G. Johnson, E. Ippen, J. D. Joannopoulos, M. Soljac̆ić, “Purcell effect in nonlinear photonic structures: A coupled mode theory analysis,” Opt. Express 16, 12523–12537 (2008). [CrossRef] [PubMed]
- J. Bravo-Abad, E. P. Ippen, M. Soljac̆ić, “Ultrafast photodetection in an all-silicon chip enabled by two-photon absorption,” Appl. Phys. Lett. 94, 241103 (2009). [CrossRef]
- H. Hashemi, A. W. Rodriguez, J. D. Joannopoulos, M. Soljac̆ić, S. G. Johnson, “Nonlinear harmonic generation and devices in doubly-resonant Kerr cavities,” Phys. Rev. A 79, 013812 (2009). [CrossRef]
- J. Bravo-Abad, A. W. Rodriguez, J. D. Joannopoulos, P. T. Rakich, S. G. Johnson, M. Soljac̆ić, “Efficient low-power terahertz generation via on-chip triply-resonant nonlinear frequency mixing,” Appl. Phys. Lett. 96, 101110 (2010). [CrossRef]
- J. A. Kong, Electromagnetic wave theory (EMW Publishing, 2005).
- T. Lu, C. Kao, H. Kuo, G. Huang, S. Wang, “CW lasing of current injection blue GaN-based vertical cavity surface emitting laser,” Appl. Phys. Lett. 92, 141102 (2008). [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.