Abstract
We designed, fabricated and characterised electrically injected quantum cascade lasers with photonic crystal reflectors emitting at terahertz frequencies (3.75 THz). These in-plane emitting structures display typical threshold current densities of 420 A/cm2 and output powers of up to 2 mW under pulsed excitation. The emission characteristics are shown to be robust, as with increasing current the emission remains singlemode with no drift in wavelength, this results from the narrow reflectivity band of the photonic crystal reflectors.
©2005 Optical Society of America
1. Introduction
Since their inception enormous progress has been made in both Quantum Cascade Lasers (QCLs) [1, 2] and Photonic Crystals (PhCs) [3, 4, 5, 6]. This is due in large part to the advances in semiconductor growth and processing technology. QCLs have steadily progressed towards low threshold, high power, room temperature operation in the Mid-Infrared (MIR) spectral range [2]. More recently, advances have been made for the far-infrared (30–300 μm, i.e. 1–10 THz) [7, 8, 9]. Likewise, PhC losses have been drastically reduced, for example at telecom wavelengths (1.55 μm) reported losses [10] are as low as 2.0 dB/cm [11] for high index contrast waveguide structures. Thus, PhCs are now becoming a viable solution for integrated optics applications.
Although extensive applications exist for Terahertz (THz) frequencies such as imaging chemical sensing and astronomy [12], progress is needed to create versatile and bright sources, focussing not only on generation but also on compact efficient light guiding solutions. The necessary material growth thicknesses of epitaxial layers for vertical refractive index guiding at THz frequencies make them impractical. PhCs enable the control of light propagation on its wavelength scale and as such can offer versatile, small all on-chip light guiding solutions for THz QCLs. There are two aspects which particularly suit a marriage between PhCs and QCLs: Firstly, as QCLs use only one carrier type for light generation, in our case electrons, no surface recombination occurs. Secondly, as QCLs are intersubband devices the polarization selection rules require that they emit in the Transverse Magnetic (TM) direction (i.e. electric field perpendicular to the growth direction). This results in emission that is predominantly in-plane, which is well adapted to planar PhCs. Moreover, at these long wavelengthis the impact of fabrication imperfections, which still limits PhC performances [10] is reduced.
Control of emission from QCLs through gratings was introduced at MIR wavelengths as early as 1997, when Faist et al. [13] used Distributed Feedback lasers (DFBs) to obtain single mode operation for gas sensing applications; DFBs used to obtain single mode operation normally result in linewidths limited by the experimental set-up resolution. Recently a DFB-QCL operating at 63 μm [14] was demonstrated by etching a grating into the upper cladding, periodically changing the penetration depth of the optical mode into the metallic plasmon guiding surface and thus its effective refractive index. As QCLs emit in TM they cannot easily be configured as a vertical cavity surface emitter. However, DFBs can be used on as a second order grating to diffract the emission vertically, such that single mode vertical emission can be obtained [15].
So far little work has been done to incorporate PhCs in QCL design. In the seminal work by Colombelli et al. [6], a MIR (8 μm) single plasmon guided QCL structure used PhCs to enable lasing through a feedback process, i.e. propagation in a slow light mode in the transmission band of the PhC structure, and vertical extraction by acting as diffraction grating. In this configration several microcavity lasers can be produced on a single chip to enable high output power. In contrast to this work where PhCs were used by tailoring their equi-frequency surfaces in their transmission bands, here PhCs are used at their ‘stopband’ energies, i.e. as reflectors, and are positioned on either side of the QCL active region. Using in-plane PhCs will allow their unique dispersion properties to be exploited in the future e.g. by making lensed mirrors.
2. Design and fabrication
To establish a suitable periodic perturbation of the dielectric is the major challenge in creating PhCs for THz-QCLs. The effect of a periodic perturbation of the dielectric on the dispersion curve of a homogeneous material are investigated by Plane Wave Expansion (PWE) calculations [16]. Figure 1 shows the dispersion curve for TM polarised light through a (a) homogeneous material with refractive index n = 3.9. Figure 1 (b) and (c) show the dispersion of PhCs consisting of a triangular lattice of pillars with n = 3.9. The computation of this effective index in the pillar region is a challenging task, we use a value of 3.9 that agrees well with the experimental data even though it is slightly larger than those used in standard DFBs lasers (3.6). The fill factor f is equal to 0.4 and is defined as the surface area of the PhC motif divided by the unit cell surface area. The dispersion curves in Fig. 1(b) and (c) are that of a PhC with a small (Δn = 0.3) and large (Δn = 2.6) refractive index contrast between the pillars and the surrounding medium, respectively. Comparing Fig. 1(b) to the homogeneous medium (Fig. 1(a)), we see a splitting of the degeneracy, resulting in small energy regions of high reflectivity ‘stopbands’ that appear at the high symmetry points and are strongly dependent on the crystallographic direction within the PhC. Increasing the refractive index contrast, as is the case in Fig. 1(c), broadens the stopbands until eventually a high reflectivity exists independent of the direction of light propagation i.e. ‘full bandgap’.
The preferred method to obtain a modulation of the refractive index is by etching the pillars or holes into the dielectric structure. PhC structures for telecommunication wavelengths use planar dielectric waveguides which provide good vertical confinement of the optical mode, reducing the necessary etch depths, and facilitating the in and out light coupling with standard techniques. At THz frequencies the necessary epitaxial growth thicknesses for dielectric waveguides are impractical, single plasmonic waveguiding at the metal semiconductor interface provides a standard solution to vertical light confinement at wavelengths greater than 15 μm [17]. Single plasmon waveguide are easily fabricated and enable high power structures with little out-of-plane losses. A pictorial representation of the PhC-QCL with a single plasmon waveguide is shown in Fig. 2(a) and (b). Figure 2(c) and (d) shows the electric field intensity versus depth into our QCL structure, calculated using the transfer matrix method. Growing a heavily doped buried layer (≈ 15 μm buried into the structure) helps to increase the confinement of the electric field in the active region, as well as serving as part of the bottom contact. The optical mode overlap with the gain region (top contact region profile in Fig. 2(c)) is Γ=0.35 [18]. A plasmon between the pillars (etch depth of 15 μm) is created by depositing a metal layer. A 5nm/50nm Ti/Au layer is deposited at this lower layer, which results in a surface plasmon, see also Fig. 3(a), this helps vertical confinement and mode matching (see Fig. 2).
Since ≈ 65 % of the light is propagating in the substrate, the difference between the mode profiles in and between the pillars is such that the modulation of the effective refractive index does not considerably effect the optical mode.
As seen from Fig. 1 even a small perturbation in the refractive index will give stopbands in certain directions. Our current technology enables 4 μm pillars with an etch depth of slightly more than 15 μm, i.e. stopping at the heavily doped layer.
An important point to note is that there is a mode mismatch at the interface between the pillar and the etched regions, that can generate out-of-plane losses, see Fig. 2(c) and (d). To maintain reasonable losses and a sufficiently large bandgap a compromise fill factor of 0.4 was used in our structures [20]. However, it should be noted that out-of-plane losses are expected to be less for TM polarisation than TE [5]. Hence, the fact that QCL emission is TM naturally reduces losses in these structures when compared to predominantly TE emitting structures.
To design the optimum PhC structure PWE and 2D finite difference time domain calculations of the band diagrams and transmission spectra were used. Pillars (as opposed to holes) consistently showed larger bandgaps for TM polarisation; it is well known that isolated areas of high index material exhibit larger bandgaps for TM polarised light in contrast to TE polarised light, which favours a network of connected high index material [21].
PWE calculations for triangular and square arrays were performed for different fill factors along all the main crystallographic directions for square and triangular arrays, see Fig. 1(d). The bandgap versus fill factor for triangular arrays give larger bandgaps for our structures; in general triangular PhC arrays display larger bandgaps due to their higher symmetry.
The GaAs/Al0.15Ga0.85As QCL structure was grown using molecular beam epitaxy on a semi-insulating GaAs substrate. It consists of 120 periods of the active regions, bounded by two heavily doped regions which serve as contacts for the structure, see Fig. 3(a). The active region was based on the ‘bound-to-continuum’ design combined with a longitudinal optical phonon extraction mechanism as described in reference [18].
The devices were made using standard fabrication methods: The ridge waveguide is obtained by wet etching through the whole active region down to the heavily doped. The bottom contact (area 220 μm × 1 mm) is composed of an alloyed (400°C 1 min) Ni/Ge/Au/Ni/Au (10/20/50/20/200 nm) layer. The top Schottky contact area (120 μm × 1 mm) consists of a Ti/Au (3/400 nm) layer. This contact, together with a 700 nm thick heavily doped (2 × 1018 cm-3) buried layer, forms a single plasmonic waveguide.
Due to the long wavelengths, the corresponding PhC lattice constants ‘a’ are large, ≈ 21 μm and a fill factor of 0.4 results in the smallest pillar radius ≈ 4 μm, which is well within the standard UV photolithography fabrication tolerances (≈ 0.5μm). However, due to the large number of growth periods in vertical direction required for the QCL and the long emission wavelength, large etch depths are needed for the PhC to impact on the optical mode. To obtain large etch depths a mask transfer into 800 nm of SiO2 was performed; a standard photolithographic procedure defined the pillar mask pattern in the photoresist. A CHF3 based Reactive Ion Etch (RIE) was used to transfer the photoresist mask into the SiO2. After removing the photoresist the pillars were etched by a Cl2 - based RIE. Approximately two hours are needed to obtain the etch depth of 15 μm, see Fig. 3(c).
A schematic of the final structure can be seen in Fig. 3(a) along with Scanning Electron Microscope (SEM) pictures of the top (Fig. 3(b)) and side (Fig. 3(c)) of the structures. Very good side wall verticality was obtained for the entire 15 μm depth.
3. Experimental results
The first test structure consists of a QCL bounded by 4 rows of a ΓK oriented triangular lattice PhC, i.e. light propagates in the ΓK direction of the PhC, with a period of 21 μm. Lasing action is observed under both continuous wave and pulsed (duty cycle of 0.008 %) excitation at 9 K. The current voltage and output light versus injected current curves for a structure under pulsed excitation are presented in Fig. 4(a). The output light peak power was measured by collecting the laser emission using a light pipe and sending it to a broadband, calibrated thermopile. The performance of the laser is good: a lasing threshold of Jth = 421 A/cm2 and a peak optical output power of 2 mW were measured.
The current light curve shows a linear increase of the output power up to a current density of 750 A/cm2. As shown in Fig. 4(a) lasing action was observed up to liquid nitrogen temperatures (77 K) in pulsed mode.
High resolution (0.03 meV) spectra were taken using a Fourier Transform Infra-Red (FTIR) spectrometer with a cryogenically cooled bolometer lock-in detection system. Figure 4(b) shows the spectrum taken at 9K at four different currents each marked with a symbol on Fig. 4(a). These spectra display singlemode emission at λ = 81.7 μm until roll over unlike their Fabry-Pérot counterparts that show blue shift of the order of the mode envelope due to a Stark shift of the lasing transition [18]. Hence, a stable singlemode emission is shown for the whole dynamic range of the PhC-QCL device.
Lasing action is observed at u = 0.26 in reduced energy units, that corresponds exactly to the degeneracy point at ΓK point, shown within the circle on Fig. 1(a) and (b); [22] i.e. the solution of the first order Bragg condition. There is a high reflectivity at this point as the degeneracy is broken by a small periodic perturbation, such as that introduced by our planar PhCs. The narrow stopband has a higher reflectivity than the cleaved or etched facets as is attested to by the lasing action that occurs from reflections at the PhC. Since the reflectivity is high within a narrow energy band along the ΓK direction singlemode emission with no wavelength drift is observed even when the gain curve shifts. In fact, as long as the width of the high reflectivity stop band is smaller than the free spectral range of the optical cavity the emission will remain singlemode.
Several tests were performed to ensure lasing resulted from reflections at the PhCs, and not from reflections of either the cleaved facet or the etched facet of the ridge waveguide, One confirmation is given by the fact that the sub-threshold Fabry-Pérot fringe spacings, seen as beating in the interferogram of Fig. 5(a) and as fringes spacing in Fig. 5(b), corresponds to the correct optical path distance of 1.1 mm (see Fig. 3(b)) when assuming an effective index for the guided mode of nvg = 3.9 [18]. Having a 1 mm long ridge waveguide (see Fig. 3) implies that the region of high reflectivity is inside the pillar structure. Moreover, the spacing of the sub-threshold Fabry-Pérot fringes of a sample cleaved just after the pillars on one side of the laser but 2.5 mm from the ridge waveguide on the other side of the laser were analysed. The fringes continue to correspond to an optical path length of 1.1 mm. This definitively precludes lasing due to reflections at the cleaved facets. To further rule out the possibility that lasing could be attributed to reflections at the etched facets of the ridge waveguide, a device was fabricated without pillars, i.e. a ridge waveguide bounded by a plateau etched to a depth of where the pillars were in the previous structure. In this case the sub-threshold Fabry-Pérot fringe spacings correspond to an optical path length of that of the cleaved facets (i.e. 1.4 mm), see Fig. 3(b). Thus, we conclude that the lasing action does not occur from reflections at either the cleaved facet, nor the etched facet, but at the pillars.
The second structure measured is the twin of the first structure, except that the PhC is aligned in the ΓM direction. Figure 6(a) shows the current versus voltage and the output light versus injected current curves. The structure exhibits a 24 % higher threshold, Jth = 546 A/cm2, than that measured in the ΓK device. The high optical saturation which occurs at 875 A/cm2 can be explained by the fact that the cavity is broadband and the peak gain can be exploited. The threshold currents and saturation densities of the structures are typical for Fabry-Pérot THz-QCLs [18]. When compared to the ΓK device, not only is the threshold higher, but more striking is the contrasting behaviour of the spectra, see Fig. 6(b). Here, the emission is clearly multimode, and shows a Stark blue shift of the lasing peak with increasing current. The lasing wavelength just above threshold is 82.2 μm and shifts to 79.7 μm just before roll over [23]. Hence, the small periodic dielectric perturbation introduced by our 2D-PhC has a negligible effect on the reflectivity seen by the lasing wavelength. These results concur with the fact that there is a higher threshold obtained from the ΓM than the ΓK orientated PhCs. The lasing behaviour in the ΓM orientated PhC structure suggests a flat low broadband reflectivity, thus unsurprisingly when the peak of the gain curves shift, the mode hops to the next lasing mode. The free spectral range between these lasing modes, corresponding to an optical path of 1.3 mm, suggest that lasing occurs due to reflection at the cleaved facets.
4. Conclusion
To conclude, we have designed and fabricated electrically injected all on chip THz PhC-QCLs. Two different symmetry directions of the PhC triangular lattice were tested, ΓM and ΓK. Dramatically different results were obtained: For light propagating in the ΓK direction singlemode emission with no wavelength drift with increasing current was observed whilst for light propagation in the ΓM direction multimode emission was blue shifted with increasing current. These results were explained by the high reflectivity points in the dispersion curve of the PhC. The advantage of using 2D-PhCs over 1D Bragg reflectors is that they enable a better transverse mode selection even when applied to a simple ridge structure, as demonstrated by the differenct emission spectrum when using ΓM and ΓK orientated PhCs. These structures also have the advantage that they can be easily placed in a coupled cavity configuration [24]. Furthermore they pave the way for exploitation of the unique in-plane dispersion properties which would allow lensing and self-collimation effects etc. Further work is being carried out to increase the PhC effect by means of a metal-metal confinement between the pillars.
We would like thank to the Swiss National Center of Competence in Research (NCCR) for providing funding through the Quantum Photonics project.
References and links
1 . J. Faist , F. Capasso , D. L. Sivco , C. Sirtori , A. L. Hutchinson , and A. Y. Cho , “ Quantum Cascade Laser ,” Science 264 , 553 – 556 ( 1994 ). [CrossRef] [PubMed]
2 . M. Beck , D. Hofstetter , T. Aellen , J. Faist , U. Oesterle , M. Ilegems , E. Gini , and H. Melchior “ Continuous Wave Operation of Mid-Infrared Semiconductor Laser at Room Temperature ,” Science 295 , 301 – 305 ( 2002 ). [CrossRef] [PubMed]
3 . S. John , “ Strong localization of photons in certain disordered dielectric superlattices ,” Phys. Rev. Lett. 58 , 2486 – 2489 ( 1987 ). [CrossRef] [PubMed]
4 . E. Yablonovitch , “ Inhibited Spontaneous Emission in Solid-state Physics and Electronics ,” Phys. Rev. Lett. 58 , 2059 – 2062 ( 1987 ). [CrossRef] [PubMed]
5 . K. Busch , S. Lölkes , R. B. Wehrspohn , and H. Föll Photonic Crystals: Advances in Design, Fabrication and Characterization ( Wiley-VCH, Weinheim, 2004 ).
6 . R. Colombelli , K. Srinivasan , M. Troccoli , O. Painter , C. F. Gmachl , D. M. Tennant , A. M. Sergent , D. L. Sivco , A. Y. Cho , and F. Capasso “ Quantum ascade Surface-Emitting Photonic Crystal Laser ,” Science 302 , 1374 – 1377 ( 2003 ). [CrossRef] [PubMed]
7 . M. Rochat , D. Hofstetter , M. Beck , and J. Faist “ Long-wavelength ( λ ≈ 16 μ m), room-temperature, single-frequency quantum-cascade lasers based on a bound-to-continuum transition ,” Appl. Phys. Lett. 79 , 4273 – 4271 ( 2001 ). [CrossRef]
8 . R. Köhler , A. Tredicucci , F. Beltram , H. E. Beere , E. H. Linfield , A. G. Davies , D. A. Richie , R. C. Lotti , and F. Rossi “ Terahertz semiconductor-heterostructure laser ,” Nature 417 , 156 – 159 ( 2002 ). [CrossRef] [PubMed]
9 . B. S. Williams , S. Kumar , Q. Hu , and J. L. Reno “ Operation of terahertz quantum-cascade lasers at 164 K in pulsed mode and at 117 K in continuous-wave mode ,” Opt. Express 13 , 3331 – 3339 ( 2005 ). http://www.opticsexpress.org/abstract.cfm?URI=OPEX-13-9-3331 [CrossRef] [PubMed]
10 . R. Ferrini , D. Leuenberger , M. Mulot , M. Qiu , J. Moosburger , M. Kamp , A. Forchel , S. Anand , and R. Houdre “ Optical Study of Two-Dimensional InP-Based Photonic Crystals by Internal Light Source Technique ,” IEEE J. Quantum Electron. 38 , 786 – 799 ( 2002 ). [CrossRef]
11 . E. Kuramochi , M. Notomi , S. Hughes , L. Ramunno , G. Kira , S. Mitsugi , A. Shinya , and T. Watanabe “ Scattering Loss of Photonic Crystal Waveguides and Nanocavities induced by Structural Disorder ,” CLEO Pacific Rim Tokyo Japan CTuE1-1 ( 2005 ).
12 . P. H. Siegel “ Terahertz Technology ” IEEE Trans. Microwave Theory Tech. 50 , 910 – 928 ( 2002 ). [CrossRef]
13 . J. Faist , C. Gmachl , F. Capasso , C. Sirtori , D. L. Sivco , J. N. Baillargeon , and A. Y. Cho “ Distributed feedback quantum cascade lasers ,” Appl. Phys. Lett. 70 , 2670 – 2672 ( 1997 ). [CrossRef]
14 . L. Mahler , R. Kohler , A. Tredicucci , F. Beltram , H. E. Beere , E. H. Linfield , D. A. Ritchie , and A. G. Davies “ Single-mode operation of terahertz quantum cascade lasers with distributed feedback resonators ,” Appl. Phys. Lett. 84 , 5446 – 5448 ( 2004 ). [CrossRef]
15 . D. Hofstetter , J. Faist , M. Beck , and U. Oesterle “ Surface-.emitting 10.1 μ m quantum-cascade distributed feedback lasers ,” Appl. Phys. Lett. 75 , 3769 – 3771 ( 1999 ). [CrossRef]
16 . M. Plihal and A. A. Maradudin “ Photonic band structure of two-dimensional systems: The triangular lattice ,” Phys. Rev. B 44 , 8565 – 8571 ( 1991 ). [CrossRef]
17 . C. Sirtori , C. Gmachl , F. Capasso , J. Faist , D. L. Sivco , A. L. Hutchinson , and A. Y. Cho “ Long-wavelength ( λ ≈ 8 – 11.5 μ m) semiconductor lasers with waveguides based on surface plasmons ,” Opt. Lett. 23 , 1366 – 1368 ( 1998 ). [CrossRef]
18 . G. Scalari , N. Hoyler , M. Giovannini , and J. Faist “ Terahertz bound-to-continuum quantum-cascade lasers based on optical-phonon scattering extraction ” Appl. Phys. Lett. 86 , 181101 – 3 ( 2005 ). [CrossRef]
19 . M. Meier , A. Mekis , A. Dodabalapur , A. Timko , R.E. Slusher , J. D. Joannopoulos , and O. Nalamasu “ Laser Action from two-dimensional distributed feedback in photonic crystals ,” Appl. Phys. Lett. 74 , 7 – 9 ( 1999 ). [CrossRef]
20 . R. Ferrini , R. Houdre , H. Benisty , M. Qui , and J. Moosburger “ Radiation losses in planar photonic crystals: two-dimensional representation of hole depth and shape by an imaginary dielectric constant ,” J. Opt. Soc. Am. B 20 , 469 – 478 ( 2003 ). [CrossRef]
21 . J. D. Joannopoulos , R. D. Meade , and J. N. Winn , Photonic Crystals: Moulding the Flow of Light ( Princeton University Press, New Jersey, 1995 ).
22 . Note that the point appears to be located at the ΓM point in Fig. 1(a) and (b) due to a peculiarity of band folding in the reduced brillouin zone scheme where ΓM fold onto ΓM, but ΓK, folds on to ΓK, ΓM ,ΓK
23 . The spacing between the modes can be seen clearly on the spectra taken at 0.8 A (blue line Fig. 5(b)). The optical path calculated from the sub-threshold Fabry-Pérot fringes using an n gr = 3.9 is 1.3 mm that suggests that the highest reflectivity is at the cleaved facets.
24 . S. Mahnkopf , R. Marz , M. Kamp , G.-H. Duan , F. Lelarge , and A. Forchel “ Tunable Photonic Crystal Coupled Cavity Lasers ,” IEEE J. Quantum Electron. 40 , 1306 – 1314 ( 2004 ). [CrossRef]