High performance 4.7 THz GaAs quantum cascade lasers based on four quantum wells

GaAs/AlGaAs quantum cascade lasers based on four quantum well structures operating at 4.7 THz are reported. A large current density dynamic range is observed, leading to a maximum operation temperature of 150 K for the double metal waveguide device and a high peak output power more than 200 mW for the single surface plasmon waveguide device. A continuous wave, single mode, third order distributed feedback laser with a low electrical power dissipation and a narrow far-field beam pattern, which is required for a local oscillator in astronomy heterodyne spectrometers, is also demonstrated.

Quantum cascade lasers (QCLs) are unipolar semiconductor lasers using carrier recycling in the active region [1]. Making use of intersubband transitions in quantum wells (QWs) enables us to design a large optical gain coefficient covering a wide spectral range. In addition, unique waveguide structures make it possible to lase in the terahertz (THz) spectral region [2][3][4], leading to novel spectroscopic applications such as a chemical remote sensing and imaging [5,6].
As demonstrated [7,8], THz QCL is especially a good candidate to be used as a local oscillator in THz heterodyne spectrometer, that is a powerful tool for radio astronomy that studies space and planetary objects by detecting various spectral lines. There are several important spectral lines in THz region for astronomical observation: for example, the CII fine structure line at 1.9 THz, and the OH fine structure line at 1.8, 2.5, and 3.5 THz [9]. In particular, the fine structure line of the neutral atomic oxygen at 4.745 THz is of great interest since it is an important cooling line in the interstellar medium that relates to star bursts [9][10][11]. Compared to other existing THz coherent sources such as a gas laser, THz QCL combines high output power, high spectral purity, spectral tunability, low electrical power dissipation and compactness, which are highly required for usage in airborne platform laboratory. In the operation frequency range of the state-of-art THz QCLs, 4.745 THz is close to the highest operation frequency of these devices (5.4 THz) [12]. Unlike around 3 THz where QCL performance peaks and a number of studies have been reported [13], there are few reports in the higher frequency region [10,12,14,15]. Here we exploit THz QCLs operating around 4.6-4.8 THz for application to local oscillators in the heterodyne spectrometer.
It has been generally accepted that the laser performance of high frequency THz QCLs is limited by a nonradiative energy relaxation via a thermally activated longitudinal optical (LO) phonon scattering whereas a parasitic current is the source for limited performance of the low frequency THz QCLs, resulting in a sweet spot operation frequency around 3 THz [13,16]. A thermally activated LO phonon scattering time (τ) for electrons in the laser upper state is expressed as exponentially decreases at elevated temperature because thermally activated electrons tend to emit an LO phonon because of the smaller energy (= E LO −E ul ) difference between E LO (≈36 meV for GaAs) and E ul (≈20 meV). Hence threshold current density ( J th ) rapidly increases, resulting in a lower characteristics temperature (T 0 ) of J th . When J th exceeds the maximum current density ( J max ), laser action ceases by decreasing the injection efficiency due to formation of electric field domain boundaries. As a consequence, the maximum operation temperature (T max ) would be lower because the current density dynamic range defined as ( J max − J th )/J th significantly decreases with temperature. For improvement upon it, we take a simple approach to increase the current density dynamic range of THz QCLs: a four QW structure based on GaAs/Al 0.15 Ga 0.85 As material, where one QW is bridged between the active QWs and the electron extraction QW, is employed because it exhibits an inherently smaller parasitic current ( J para ) [18]. J para is attributed to the wavefunction coupling between the ground state (denoted by 'g' in figure 1(a)) in the electron extraction QW and the excited state ('e' in figure 1(a)) of the adjacent electron extraction QW in the next period. Because it occurs before pumping the laser upper state, J para provides a minimum of J th . Since the bridged QW spatially separates the active QWs and the electron extraction QW, J para would be expected to decrease due to their smaller wavefunction coupling. For the four QWs structure, the injection barrier thickness is an important parameter to increase the current density dynamic range since it determines J max and J para . Although a thinner injection barrier increases J max , J para also becomes large enough to limit J th because of the short lifetime (≈0.25 ps) of the excited state in the electron extraction QW. Figure 1(c) shows the computed J max , J para and current density dynamic range as a function of coupling energy when the ground state in the electron extraction QW is aligned with the laser upper state as depicted in figure 1(b). The Kazarinov-Suris model [19,20] was here used with an assumption that J th is limited by J para . We also used the computed upper state lifetime (= 5 ps) neglecting a radiative current and assumed that the in-plane electron momentum scattering time is 1 ps [21]. The normalized oscillator strength of the laser transition ¢ = ( ) f 0.67 was selected since a high selective carrier injection into the laser upper state would be possible for a large E ul . Although employing a thinner injection barrier increases J max , it leads to a smaller current density dynamic range since J para also increases rapidly with the coupling energy. In this work the injection barrier thickness (5.1 nm) with a coupling energy (ħΩ) of 0.85 meV is selected keeping J max large for a higher output power together with a large current density dynamic range. The designed emission energy is 19.2 meV. Starting from the injection barrier, the layer thickness in nanometer is 5.  Computed parasitic current density ( J para ), maximum current density ( J max ) and current density dynamic range (( J max −J para )/J para ) as a function of the coupling energy between the injection layer and the active layer. We used the computed upper state lifetime assuming that it is not dependent on the injection barrier thickness.
The laser structure was grown on SI-GaAs (001) substrate by solid-source MBE. The laser structure started with a 200 nm thick Al 0.5 Ga 0.5 As etching stop layer, followed by a 400 nm thick n-GaAs bottom contact layer (n=3.0×10 18 cm −3 ). Then 200 periods of GaAs/AlGaAs four QW active region (total thickness: 11.7 μm) was grown. At the end, a 40 nm thick n-GaAs top contact layer (n=3.8×10 18 cm −3 ) and a 20 nm thick n-InGaAs cap layer (n=1.0×10 19 cm −3 ) was grown on top of the active region. X-ray diffraction spectra show that the measured period thickness is in good agreement to the design. The MBE grown wafers were then processed into ridge laser structures with a double metal waveguide or a single surface plasmon waveguide. For a double metal waveguide, in order to decrease the waveguide losses, the bottom n-GaAs contact layer was etched away after removing the GaAs substrate. The processed wafers were cleaved into laser devices with ridge widths of 160 μm and cavity lengths of 1.0-3.0 mm. The optical power measurements were carried out by the Thomas Keating Ltd THz power meter. The power meter was placed in front of the cryostat window as close as possible without any collection optics because of the large surface of the detector element. Also we did not correct for the transmission of the polymethylpentene cryostat window. Figure 2(a) shows temperature dependent light (L)-current (I)-voltage (V ) characteristics of the Cu-Cu double metal waveguide device with a cavity length of 2.0 mm. Current started flowing by subband alignment between the ground state in the electron extraction QW and the lower states of the laser transition. Lasing started with J th =200 A cm −2 right after pumping the parasitic transport channel (V=7-9 V). J max reached 920 A cm −2 where the negative differential resistance voltage appeared at V=14.2 V, which is in close agreement with the designed one (V=13.8 V). The device provided a large current density dynamic range of 3.6, which is significantly larger than previously reported at high frequency THz QCLs (≈1 for 4.4 THz QCL [22] and ≈2 for 4.7 THz and 5.4 THz QCLs [12,15]). The maximum output peak power (P) was 11 mW with the slope efficiency of 7 mW A −1 at 10 K. T max was 150 K, which, to our knowledge, is the highest operation temperature for higher frequency THz QCLs [23]. The device equipped with the Au-Au double metal waveguide with a ridge width of 160 μm and a cavity length of 2.25 mm long exhibited similar L-I-V characteristics with lower T max (=142 K) (not shown). Figure 2(b) depicts J th as a function of temperature. Fitting data from 90 to 140 K via the empirical equation ( J th =J 0 exp (T/T 0 )) expressing T 0 of J th results in T 0 =60 K, which is lower than the value of T 0 (80 K) of the 3 THz QCLs [18] based on the same four QW scheme. As derived from equation (1), this is due to the larger thermally activated laser upper state lifetime originated from the smaller E LO −E ul . To evaluate T 0 , we computed temperature dependence of J th by the density matrix (DM) model taking into account of the kinetic energy balance for electrons [24]. In this model, electron-impurity scattering, interface roughness scattering, and alloy disorder scattering are included as an elastic scattering mechanism. The electron temperature (T e ) is defined as an averaged electron temperature over all subbands in one period of the active region. Here the total loss (= α w +α m , where α m is mirror loss and α w is waveguide loss) is assumed to be 30 cm −1 [25] which is independent on temperature. Also the electron extraction is supposed to be done not by resonant tunneling but by resonant LO phonon scattering. T 0 =69 K obtained by simulated J th from 90 to 140 K is in close agreement with the experimental value (= 60 K). In addition, the simulated J th shows a weak temperature dependence in the temperature range from 10 to 50 K comparable to the experimental J th ; the computed T e is almost constant (≈50 K) until the lattice temperature reaches 50 K, thus suppressing to exponentially decrease the upper state lifetime as shown in equation (1) if T e increases with the lattice temperature. Thus those facts imply that our DM model predicts right T e and reproduces thermally activated carrier scattering processes in the active region. However, it underestimates the absolute value of J th by a factor of two in this temperature range most probably due to underestimation of the elastic scattering rate. This leads to the larger simulated current density dynamic range (= 5.1) compared to the experimental one (= 3.6).
In addition to the observed T max , a remarkable feature of the large current density dynamic range is the capability of high P when a large out-coupling waveguide such as a single surface plasmon waveguide is employed. Figure 3(a) shows L-I-V curves of the device equipped with a single surface plasmon waveguide. The device had 150 μm wide ridge and 1.0 mm long cavity. Current pulses of 75 ns width with a repetition frequency of 10 kHz were used. One of the facets was high-reflection coated by Al 2 O 3 (200 nm)/Au (100 nm) films. I-V characteristics exhibited the same current density dynamic range (≈3.6) as that of the double metal waveguide. J th was 300 A cm −2 at 10 K. As expected, the device delivered P=220 mW at 10 K, which is among the highest peak power so far obtained from the same device size [22,26]. Moreover, as depicted in figure 3(b), another device having 200 μm wide and 6.5 mm long, with uncoated facets emitted P≈600 mW from a single facet, corresponding to P≈1.2 W for both facets. The slope efficiency of the 1.0 mm long device was 303 mW A −1 at 10 K. Taking into account of the ratio α m /(α w +α m ), the slope efficiency seems to be larger than expected one of the double metal waveguide device, probably due to the lower collected power by the non-ideal far-field pattern of the double metal waveguide. Emission spectrum of the 1.0 mm long device covers the broad spectral range from 4.6 to 5.0 THz including the target frequency (4.74 THz) (not shown). T max of the 1.0 mm long device was 105 K.
Finally we fabricated third order distributed feedback (DFB) laser structures [27] for continuous wave (CW) and single mode operation which is required for a local oscillator in a heterodyne spectrometer. To decrease electrical power dissipation, the active region thickness was thinned down to 9 μm by wet etching after removing the GaAs substrate. As depicted in figure 4(a), the DFB cavity is designed as a lateral corrugated grating with a square tooth shape and fabricated by standard photolithography and inductively coupled reactive ion etching. The width of the DFB cavity is 20 μm and the total length limited by the number of the grating period (=19) is ≈800 μm. Figure 4(b) shows emission spectra measured at 10 K. Due to the low transverse mode selectivity, we observed a few mode emissions from the as-processed devices. To suppress the unwanted modes, a thin Ni film (10 nm), highly absorptive in the THz region, was deposited on the device surface by a self-aligned method. The thin Ni film cleaned emission spectra since it increases waveguide loss on the absorbing end part of the DFB cavity. Although the procedure yielded a 10% increase in J th and a reduced P as depicted in figure 4(c), we achieved a single mode operation with side mode suppression above 20 dB. This technique is thus important to increase a yield of the single mode laser operating at the target frequency. Figure 4 displays (d) far-field laser beam pattern and (e) CW emission spectra of one of the QCLs in the DFBs array. All these measurements were done at 40 K. A single hot spot was observed in the angle area of 30°×15°. Judging from the far field pattern and the estimated effective index of the DFB cavity (= 3.06), the observed laser mode is attributed to the higher band edge third order DFB mode. The laser emission frequency was 4.745 THz with an output power of 0.5 mW, hitting the target frequency of the neutral atomic oxygen fine structure line. The temperature range of 40-50 K with an electrical dissipation of 1 W is accessible by a stirling-type cryostat which is useful in space-based observatories. We achieved the highest CW output power of 0.9 mW at 45 K with an electrical power dissipation of 1.2 W at 4.788 THz with the DFB QCLs being cooled in the stirling-type cryostat. The maximum CW operation temperature was 75 K.
The long-term stability of the laser emission frequency is a crucial point for an astronomical heterodyne spectrometer. The laser emission frequency of third order DFB THz QCLs is known to be sensitive to adsorption of molecules on the DFB cavity side wall [28] because of a strong leakage of the optical field by the narrow DFB cavity width; residual gas molecules in the pumped cryostat adsorb and condense on the cavity side wall, leading to a change of the effective index of the waveguide. To avoid it, a small metal cage (shown in figure 5(a)) was used here to trap gas molecules on it. Figures 5(b) and (c) depicts long-term time trace of detector signal from a Michelson interferometer with a fixed mirror position. The pressure inside the cryostat was 9.0×10 −8 mbar. Without the cage, fringes on the detector signal were visible, revealing presence of a long term frequency drift of the laser frequency. The drift frequency between peaks corresponds to 500 MHz. On the other hand, those fringes were found to disappear with the cage and the long-term laser frequency stability was dramatically improved. This demonstrates that protecting the cavity side wall is important for long-term laser frequency stability of third order DFB THz QCLs. Another issue is a lower beam coupling to a heterodyne mixer device in a heterodyne spectrometer. A waveguide integrated mixer is usually designed to couple to a specific beam pattern such as a Gaussian beam. As shown in figure 4(d), the far-field beam pattern of the third order DFB THz QCL is in a crescent shape, which contains a small portion of the Gaussian beam component, resulting in a lower coupling to the mixer waveguide. A large number of periods of the third order DFB with the perfect phase matching condition [29] should improve it. Other waveguides based on a photonic crystal [30], a patch antenna array [31] and a waveguide integrated block [32,33] generating a better beam shape could be thus suited, in particular, toward a high power THz local oscillator.
In conclusion, we have studied the characteristics of high frequency GaAs/AlGaAs THz QCL based on the four QW design. The observed J th and J max at 10 K were 200 A cm −2 and 920 A cm −2 , respectively. The resulting current density dynamic range (= 3.6) was about three times larger than the typical 3 THz QCL, leading to T max =150 K and P=220 mW at 10 K, which is among the best laser performance at the high frequency THz region. It shows that a large current density dynamic range improves the laser performance of THz QCL. We also demonstrated a 4.74 THz single mode, CW, third order DFB QCLs with a narrow far-field beam pattern (beam divergence of 30°×15°) and a low electrical dissipation (0.8 W) for use as a local oscillator in a To suppress unwanted cavity modes, a thin Ni film was directly deposited on the device surface by a self-aligned method. (b) Emission spectra and (c) light output-current-voltage characteristics before and after the Ni deposition. We found that the thin Ni film cleans emission spectra by increasing the waveguide loss on the absorbing end part. A single mode, continuous wave emission spectra (e) at 4.745-4.743 THz and (d) its far-field beam pattern measured at 40 K under the quasi-continuous wave operation. heterodyne spectrometer. The CW operation temperature around 50 K as well as the electrical power dissipation of 1 W is accessible by a stirling-type cryostat which is useful in space-based observatories.