Simulation and optimization of 2.6–2.8 μm GaSb-based VCSELs

We present the simulation results of threshold operation of mid-infrared GaSb-based vertical-cavity surface-emitting lasers (VCSELs) obtained with the use of comprehensive fully self-consistent optical-electrical-thermal-recombination numerical model. The results show that by a proper design of VCSEL structure and composition of the active region it is theoretically possible to achieve room-temperature (RT) threshold operation for wavelength of 2.8 μm which is about 0.2 μm longer than those reported so far in the literature for III-V VCSELs with type-I quantum wells. Calculated values of the RT threshold current were equal to 2.5 and 4.0 mA for tunnel-junction diameters of 2 and 4 μm, respectively.


Introduction
Mid-infrared vertical-cavity surface-emitting lasers (VCSELs) have a great potential as low cost, ultra-low threshold, small beam divergence light sources in distant air monitoring. Although there are several important gases, like C 2 O, N 2 O, H 2 S, CH 4 in the midinfrared spectral range, the longest wavelength reported to date for electrically pumped room-temperature (RT) continues-wave (CW) VCSELs is 2.62 lm (Arafin 2010). In this work we present the results of numerical analysis carried out for devices designed to emit This article is part of the Topical Collection on Numerical Simulation of Optoelectronic Devices 2016.
Guest edited by Yuh-Renn Wu, Weida Hu, Slawomir Sujecki, Silvano Donati, Matthias Auf der Maur and Mohamed Swillam. the light from the 2.6-2.8 lm wavelength range. Due to the fact that the single-fundamental-mode operation is particularly suitable for absorption spectroscopy, we concentrate on suppression of the higher-order transverse modes in the modelled devices.

The laser structure
The structure under consideration was the GaSb-based quantum-well (QW) GaInAsSb/ GaSb VCSEL (Fig. 1) similar to the device proposed in ).Its intentionally undoped active region is assumed to be composed of 8-nm Ga 1Àx In x As 0:15 Sb 0:85 QWs (indium content x is equal to 0.43 and 0.49 for 1.6% and 2.0% compressively strained QWs, respectively) separated by 10-nm GaSb barriers. The active region is embedded in Al 0:15 Ga 0:85 As 0:01 Sb 0:99 waveguide and sandwiched by p-and n-type GaSb spacers. Above p-type spacer, the tunnel junction (TJ) composed of p þþ -GaSb and n þþ À InAs 0:91 Sb 0:09 is located. To minimize the absorption loss, the TJ is situated at the standing-wave node. Upper spacer is manufactured from n-GaSb. The 3-k cavity with several n þ -GaSb current spreading layers situated in cavity nodes is terminated on both sides by distributed-Braggreflectors (DBRs): the 4-pair a À Si/SiO 2 top DBR and the 24-pair AlAs 0:08 Sb 0:92 /GaSb bottom n-type DBR. Bottom DBR diameter is assumed to be equal to 100 lm, whereas the upper DBR diameter is larger by 10 lm than TJ diameter. The top contact is produced in a form of a ring of 5 lm width. It is separated from the top spacer with the 200-nm thick highly-doped n þ À InAs 0:91 Sb 0:09 contact layer. The whole bottom 100-lm diameter surface of the n-type GaSb substrate is covered by the bottom contact.

The model
Our three-dimensional optical-electrical-thermal-gain self-consistent VCSEL thresholdmodel has been adapted to simulate the CW operation of the GaSb-based VCSELs. The detailed description of the numerical algorithm applied in our model can be found in  (Sarzała 2012), whereas in (Piskorski 2015) the method to calculate the shifts in the conduction and valence bands due to strain effects is presented. Values of all necessary material parameters can be obtained with the use of formulas given in (Piskorski 2016).

The results
As a verification of correctness of the set of material data we compared the experimental and numerical results for the voltage-current characteristics (Figs. 2a and 3a) and ambient temperature dependences of the threshold current (Figs. 2b and 3b) for two VCSELs emitting in the mid-infrared wavelength range. The first structure with the emission wavelength of 2.4-lm was proposed in ), whereas the second one with the emission wavelength of 2.6-lm in ). The TJ diameters were equal to 6 lm and 5.5 lm, respectively. The main difference between these two VCSELs is the composition of the active region. In the 2.4-lm device barriers are made from AlGaAsSb, whereas in the 2.6-lm laser GaSb material is used in order to increase the valence band offset.
As can be seen in Figs. 2a and 3a, there is a very good agreement between experimental data for the voltage-current characteristics and theoretical results obtained with the use of our model. From that we can conclude that the electrical parameters (electrical conductivities, mobilities, free carrier concentrations) and the thermal conductivities have been correctly estimated. A slightly poorer agreement has been found for the ambient temperature dependences of the threshold current (Figs. 2b and 3b). However, the general trends are similar. The observed discrepancy can take its origin from the fact that in the real structure QWs might have slightly different material compositions and widths, whereas, in the simulation, they were identical to simplify the calculations. Moreover, we assumed that the carrier concentration is the same in every QW, which is not true in the real device. Both simplifications lead to broadening of the gain spectrum, and therefore have an influence on the ambient temperature dependences of the threshold current. Nevertheless, the gain and optical parameters seem to have reasonable values to perform the simulation of a device based on the antimonides.
In the above calculations five parameters were adjusted to fit the experimental data: (1) reverse saturation current density j s ¼ 1:0 A/m 2 and (2) shape factor b pn ¼ 40 V À1  (Vurgaftman 2001), its theoretical value is underestimated compared to the experimentally determined one (Chuang 2009). Blueshift, however, is the effect of annealing-induced interdiffusion on the electronic structure during the growth of the active region (Dier 2005). Energy shift caused by the above interdiffusion depends on the growth conditions, and therefore, it is specific for the given structure.
In Fig. 4 we present the values of the threshold current I th for various ambient temperatures T amb calculated for 2.6-lm VCSEL with the strain in the QW equal to ¼ À1:6%. As can be seen, with the increase of the T amb the threshold current increases for all considered TJ diameters. This behaviour can be explained with the use of Fig. 5 in which we present the gain spectra calculated for VCSELs designed to emit the light from the 2.6-2.8 lm wavelength range. For the 2.6-lm device the gain peak is shifted to a longer wavelength than the emission one. To achieve the threshold operation for higher T amb it is necessary to obtain the similar value of gain in the active region (in our calculations the threshold gain slightly increases which is a result of temperature influence on absorption coefficients taken into account). However, for higher temperatures, gain Calculated threshold current vs. ambient temperature for 2.6-lm VCSEL with 1.6% compressively strained QWs. Filled and empty symbols correspond to LP 01 and LP 11 mode operation, respectively calculated for a fixed carrier concentration decreases. Therefore, it is necessary to increase the carrier concentration in the active region which leads to the higher I th . From Fig. 4 it can be also seen that for wider apertures the temperature range, in which the LP 01 mode operation is achieved, decreases. This can be explained by the fact that for wider apertures the LP 11 mode losses calculated in respect to LP 01 mode become smaller. For the device with the same active region but designed to emit the wavelength of 2.7 lm by increasing the thickness of DBR layers and the optical length of the cavity, we obtained the stable LP 01 mode operation only for the TJ with the smallest aperture (Fig. 6a). As can be seen from Fig. 5, for 2.7-lm VCSEL and ¼ À1:6%, due to the shift of gain spectrum towards the longer wavelengths, the value of gain for the wavelength corresponding to LP 11 mode is higher than the one for LP 01 mode. This leads to the significant reduction of the LP 11 mode losses, and therefore we observe LP 11 mode operation for all considered apertures wider than 2 lm, for which the LP 11 mode losses are still enough high to provide LP 01 mode emission. We performed analogous calculations for 2.8-lm VCSEL with the strain in the QW equal to ¼ À1:6%. For this device LP 11 mode operation is more favourable than for the 2.7-lm VCSEL with the same active region. The  values of the threshold current which, we do not report here, follow the same trend where the threshold current increases and the temperature range, for which LP 01 mode operation is observed, decreases. As a result, for the 2.8-lm VCSEL, we obtained the stable fundamental mode operation only for device with TJ diameter of 2 lm One of methods to increase the LP 11 mode losses and therefore to obtain laser operation in the LP 01 mode is to shift the gain spectrum toward the longer wavelength. In VCSELs, for the situation when the wavelength corresponding to the gain peak is longer than the emission wavelength, the LP 01 mode operation is more favourable than the LP 11 one. In our calculations we shifted the gain spectrum for the 2.7-lm VCSEL (Fig. 5) by increasing the indium content in the GaInAsSb QWs from 0.43 to 0.49 which was followed by the increase in the value of compressive strain by 0.4%. As a result we obtained the stable LP 01 mode operation for TJ diameters up to 6 lm ( Fig. 6b) with the distinctly lower values of the threshold currents. Although the compressive strain of 1.6% in our initial structure is high, surface-emitting lasers with antimonide-based active regions with even higher strain in QWs were already reported. In (Rösener 2011) the solid-source molecular beam epitaxy was used to grow active region which consists of ten compressively strained GaInAsSb QWs for whose the compressive strain was determined to be of 1.8% with the use of highresolution x-ray diffraction. Higher strain (2.08%) for active region with ten GaInAsSb QWs was reported in (Holl 2016), where it is also mentioned that the thickness of each QW was equal to 10 nm. It is the same value which we used in calculations described in this paper. Although the higher strain leads to a reduction of the critical thickness at which relaxation of the QW layers occurs, it seems that active regions considered here can be realized with the use of present technology.
The structure with the 2% compressively-strained active-region was further modified in order to achieve emission around 2.8 lm. Similarly to the 2.7-lm devices, the thickness of DBR layers and n-GaSb layers in the cavity have been increased. As can be seen from Fig. 7 after the above modifications we obtained LP 01 mode operation above RT for TJ diameters up to 4 lm. Calculated values of the RT threshold current were equal to 2.5 and 4.0 mA for tunnel-junction diameters of 2 and 4 lm, respectively. The threshold current values and temperature-operation ranges for the corresponding TJ diameters are similar to those from Fig. 6a which is the result of the similar cavity-gain detuning for these lasers.

Conclusion
We performed the threshold analysis of an operation of GaSb-based VCSELs emitting in the 2.6-2.8 lm wavelength range. The simulation was performed for VCSELs with various TJ diamters and for ambient temperatures ranging from 250 K to the maximal operating temperature for a given device. The results show that by a proper design of VCSEL structure and composition of the active region it is theoretically possible to achieve the threshold operation for wavelengths about 0.2 lm longer than those reported so far in the literature. Moreover, obtaining the desired low-threshold stable single-fundamental-mode operation in these devices is also possible. Our results may be useful for designing the electrically-pumped continuous-wave III-V VCSELs able to achieve the emission wavelength of 3 lm.