Design of An Electrically Pumped Sigesn/Gesn/Sigesn Double-Heterostructure Midinfrared Laser
Design of An Electrically Pumped Sigesn/Gesn/Sigesn Double-Heterostructure Midinfrared Laser
midinfrared laser
G. Sun, R. A. Soref, and H. H. Cheng
[This article is copyrighted as indicated in the article. Reuse of AIP content is subject to the terms at: http://scitation.aip.org/termsconditions. Downloaded to ] IP:
137.132.123.69 On: Fri, 13 Mar 2015 04:36:52
JOURNAL OF APPLIED PHYSICS 108, 033107 共2010兲
I. INTRODUCTION pler because the exact layer thickness of the DHS active
region is not as critical as it is in a multiple QW device and
Silicon-germanium-tin technology is an excellent new should therefore be easier to implement. Our simulation in-
group-IV technology for the monolithic integration of active cludes the effects of nonradiative Auger recombination be-
and passive waveguided components on silicon or silicon- tween electrons and holes because that physical mechanism
on-insulator 共SOI兲.1,2 Within group-IV photonics, the is the key determinant of the diode’s temperature of opera-
group-IV electrically injected laser diode is a key component tion. Our result indicates that this DHS laser would operate
that is missing thus far, and that laser is the focus of this below room temperature. Nevertheless, with its simplicity
article. GeSn is probably the best-known Sn-containing alloy for device growth, it could very well lead to the first electri-
and for this material the infrared region is most natural for cally injected group-IV laser that can be monolithically inte-
radiative band-to-band emission because the direct band gap grated on Si or SOI. In addition, once the material growth
of this crystalline alloy is narrower than that of elemental Ge, becomes mature, such a GeSn/SiGeSn lattice-matched struc-
the semiconductor in which optically pumped Ge-on-Si las- ture can be extended to room-temperature QW structures that
ing at 1.6 m has already been observed.3 This article pre- require lower carrier densities for adequate optical gain be-
sents design-and-simulation results on an electrically in- cause of the reduced density of states and the simultaneous
jected P-type semiconductor/Intrinsic semiconductor/N-type Auger suppression.
semiconductor 共PIN兲 double heterostructure 共DHS兲 diode in
which a layer of GeSn is the active region and a SiGeSn
ternary layer forms each of the two cladding regions. The II. BAND STRUCTURE OF GESN/GESISN DHS
compositions of GeSn and SiGeSn are chosen in such a way Energy-band theory,8 plus a few FTIR absorption
that not only they form type-I band alignment but also are experiments,9 indicate that the band gap of unstrained crys-
lattice matched. The laser would be situated on a relaxed talline GeSn changes from indirect to direct as the percent of
buffer layer of GeSn-upon-silicon or SOI 共Ref. 4兲 whose ␣-Sn is increased. Since band offsets between ternary Sn-
lattice parameter is the same as that of the GeSn/SiGeSn containing alloys and Si or Ge are not known experimentally,
laser, hence the entire laser structure is unstrained. The re- we follow the assumptions made in Ref. 1 and calculate the
sults reported here compliment the previous simulations of a conduction band minima for the lattice-matched heterostruc-
1.55 m GeSn-quantum-well 共QW兲 laser,5 a strained GeSn/ ture consisting of Ge1−zSnz and a ternary Ge1−x−ySixSny based
GeSiSn QW laser,6 and a terahertz Ge/SiGeSn quantum cas- on Jaros’ band offset theory,10 which is in good agreement
cade laser.7 with experiment for many heterojunction systems. For ex-
The band-to-band laser diode simulated in this paper has ample, this theory predicts an average valence band offset,
an emission wavelength in the 1.7 to 1.8 m midwave in- ⌬E,a = 0.48 eV for a Ge/Si heterointerface 共higher energy
frared. Compared to the prior art5–7 this device is much sim- on the Ge side兲, close to the accepted value of ⌬E,a
= 0.5 eV. The basic ingredients of our band alignment cal-
a兲
Electronic mail: hhcheng@ntu.edu.tw. culation are the average 共between heavy, light, and split-off
TABLE I. Band parameters at various valleys used in the band alignment calculation 共Refs. 9–11兲.
hole bands兲 valence band offset between the two materials EX共Ge1−xSix兲 = 0.931 + 0.018x + 0.206x2 共5兲
and the compositional dependence of the band structure of
the ternary alloy. For the Ge/ ␣-Sn interface Jaros’ theory 共in electron volt兲 for Ge1−xSix alloys.14 On the other hand, the
predicts ⌬E,a = 0.69 eV 共higher energy on the ␣-Sn side兲. empirical pseudopotential calculations of Chelikovsky and
Thus, relative to the average valence band of Ge, the average Cohen place this minimum at 0.90 eV in ␣-Sn, virtually the
valence band position for Ge1−x−ySixSny is simply a linear same as its value in pure Ge.15 We thus assume, as in Ref. 1,
interpolation that the position of this minimum in ternary Ge1−x−ySixSny
alloys is independent of the Sn concentration y, and thus is
⌬E,a共Ge1−x−ySixSny兲 = − 0.48x + 0.69y. 共1兲 also given by Eq. 共5兲. Obviously, the calculation of band
structures outlined above is an approximation that is subject
Similarly, with these spin-orbit splitting values ⌬so共Ge兲
to experimental corrections as more measurements become
= 0.295 eV, ⌬so共Si兲 = 0.043 eV, ⌬so共Sn兲 = 0.800 eV,11 the
available. This implies that the compositions of Ge1−zSnz and
spin-orbit splitting for Ge1−x−ySixSny is
Ge1−x−ySixSny are necessarily adjusted in the DHS to arrive
⌬so共Ge1−x−ySixSny兲 = 0.295共1 − x − y兲 + 0.043x + 0.800y. at the band structure that is being proposed here. But it
should be pointed out that the laser behavior depends only on
共2兲
the band structure, and the results obtained in this design
The top of the valence band for Ge1−x−ySixSny can then be should be valid albeit at slightly different binary and ternary
determined as compositions.
The ␣-Sn composition dependence of the conduction
E共Ge1−x−ySixSny兲 = E,a共Ge1−x−ySixSny兲 band gaps for Ge1−zSnz at the three valleys L, ⌫, and X is first
⌬so共Ge1−x−ySixSny兲 calculated using Eqs. 共4兲 and 共5兲 to establish the crossing
+ . 共3兲 point where the ⌫-point band gap drops below that of the
3
L-point. Figure 1共a兲 shows that for ␣-Sn composition greater
The minima of the conduction band at points L and ⌫ can than z ⱖ 6%, Ge1−zSnz becomes direct band gap which can
then be calculated by evaluating the compositional depen- serve as the active layer in the DHS. Fixing at this compo-
dence of the band gaps of the ternary alloy as sition Ge0.94Sn0.06, we then looked for a lattice matched
Ge1−x−ySixSny that can be used as cladding layers for
E共Ge1−x−ySixSny兲 = EGe共1 − x − y兲 + ESix + ESny Ge0.94Sn0.06. The requirement for Ge1−x−ySixSny being the
− bGeSi共1 − x − y兲x − bGeSn共1 − x cladding is that it must form type-I band alignment with
Ge0.94Sn0.06. Such a simultaneous requirement can be satis-
− y兲y − bSiSnxy, 共4兲 fied by the additional degree of freedom in the Ge1−x−ySixSny
where both x and y can be tuned. Using the Vagard’s law for
where EGe, ESi, and ESn are the band gap of Ge, Si, and ␣-Sn,
the lattice constant of Ge1−x−ySixSny
respectively, at those points, and the bowing parameters
bGeSi, bGeSn, bSiSn have been discussed in Refs. 12 and 13. a共Ge1−x−ySixSny兲 = aGe共1 − x − y兲 + aSix + aSny, 共6兲
These values at L and ⌫ points have been given in Table I.
Finally, for the indirect conduction band minimum near where the lattice constants aGe = 5.64613 Å, aSi
the X-point, Weber and Alonso find = 5.43095 Å, and aSn = 6.48920 Å for Ge, Si, and ␣-Sn,11
[This article is copyrighted as indicated in the article. Reuse of AIP content is subject to the terms at: http://scitation.aip.org/termsconditions. Downloaded to ] IP:
137.132.123.69 On: Fri, 13 Mar 2015 04:36:52
033107-3 Sun, Soref, and Cheng J. Appl. Phys. 108, 033107 共2010兲
type-I confinement with Ge0.94Sn0.06, i.e., the band alignment respectively, where kB is the Boltzmann constant, and the
provides carrier confinement for both elections and holes electron and hole energies in the conduction and valence
when Ge0.94Sn0.06 is sandwiched by the Ge1−x−ySixSny band
cladding layers. In the material growth front, there are
currently challenges in incorporating high Sn compositions mr
Ee = 共E − Eg兲,
into the alloy, we have therefore chosen to design lasers me
based on relatively small Sn compositions. In particular,
we propose a laser made of the lattice matched mr
Ehh = 共E − Eg兲, 共9兲
Ge0.94Sn0.06 / Ge0.75Si0.15Sn0.1 DHS 共Fig. 2兲. It can be esti- mhh
mated from Fig. 1 that the band gap of the active region
are computed through the reduced effective mass
Ge0.94Sn0.06 is Eg,⌫共Ge0.94Sn0.06兲 = 0.619 eV, while the con-
duction band offset at ⌫-point ⌬Ec,⌫ = 149 meV, and that of memhh
valence band ⌬E = 54 meV. The laser device shall be grown mr = , 共10兲
me + mhh
on a relaxed Ge0.94Sn0.06 buffer on Si substrates to ensure
that the entire structure is strain free. in which me and mhh are the electron and heavy hole effec-
tive mass, respectively. Obviously, the spontaneous emission
rate rsp,e-lh共E兲dE for the radiative e-lh process can be ob-
III. CARRIER LIFETIME tained similarly. The total spontaneous emission rate per unit
volume can then be evaluated by integrating over the photon
Since the active region of Ge0.94Sn0.06 in this DHS has
energy E above the band gap Eg as
small amount of Sn, we shall use ⌫-point Ge parameters in
the following calculations by considering that the only effect
of incorporating the 6% Sn is to lower the ⌫-point of the Rrad = 冕 Eg
⬁
关rsp,e-hh共E兲 + rsp,e-lh共E兲兴dE. 共11兲
band gap to yield a direct band gap. The laser device under
consideration has a PIN DHS where the active Ge0.94Sn0.06 The radiative lifetime can then be obtained by rad = n / Rrad
region is undoped. The band-to-band lasing transitions occur where n is the carrier density. The result for a carrier density
as stimulated emissions triggered by recombination of of n = 1018 / cm3, is shown in Fig. 3 for a range of tempera-
electron-hole pairs that are injected into this region. The ac- ture. It can be seen that the radiative lifetime for a fixed
tive region in the DHS is always wide enough for the carrier carrier density is rather insensitive to temperature change,
dynamics to be treated as bulk-like. The analysis proceeds as showing a slight increase with the temperature.
[This article is copyrighted as indicated in the article. Reuse of AIP content is subject to the terms at: http://scitation.aip.org/termsconditions. Downloaded to ] IP:
137.132.123.69 On: Fri, 13 Mar 2015 04:36:52
033107-4 Sun, Soref, and Cheng J. Appl. Phys. 108, 033107 共2010兲
me mhh + me
= , ␣c = , 共14兲
mhh mhh + 2me
the threshold energy required for the second electron to par-
ticipate in the Auger process
mhh + 2me
ET = Eg , 共15兲
mhh + me
兩M a兩2 is the matrix element of coulomb interaction between
two electrons evaluated at ET,16 and introducing kT
= 冑共2me / ប2兲ET
冋 册
active region of the Ge0.94Sn0.06 / Ge0.75Si0.15Sn0.1 DHS.
ប2
⫻exp − k共k + 2kT兲 dk. 共16兲
2mekBT
The other important recombination mechanism is the
nonradiative Auger process where the recombination of an Other Auger processes involving different particles can be
electron-hole pair takes place by transferring energy and mo- evaluated following a similar procedure. Among them, the
mentum to a third particle which could be either an electron above illustrated e-hh Auger process involving a second
or a hole. We shall use the example of the Auger process electron dominates, and the result for the Auger recombina-
involving an electron as the third particle to illustrate the tion lifetime obtained by aug = n / Raug is shown in Fig. 3 for
calculation procedure whose recombination rate per unit vol- n = p = 1018 / cm3. Clearly, the Auger lifetime decreases rap-
ume can be approximated as Raug = Cn2 p with an Auger co- idly with the increase in temperature as suggested by Eq.
efficient C, and for nondegenerate semiconductor, the elec- 共13兲. At higher temperature, it surpasses radiative rate to be-
tron and hole concentrations come the dominant recombination mechanism in the active
region.
n = Nc exp − 冉 Ec − E fc
k BT
冊
, IV. OPTICAL GAIN AND THRESHOLD CURRENT
For a given injected carrier concentration n = p, we can
冉 冊
calculate the optical gain at a photon energy E due to e-hh
E f − E recombination as16
p = N exp − , 共12兲
冉 冊
k BT
e2ប2兩M b兩2 2mr 3/2
␥e-hh = 共E − Eg兲1/2关f c共Ee兲
where Nc and N are the effective density of states for con- 0m20cn̄E ប2
duction and valence band, respectively. According to Ref. 16, + f v共Ehh兲 − 1兴. 共17兲
the Auger rate per unit volume due to the recombination of
an e-hh pair involving a second electron as the third particle A similar expression can be obtained for the e-lh process
can be calculated as ␥e-lh. Thus, the total optical gain is then ␥ = ␥e-hh + ␥e-lh. The
profiles of optical gain are shown in Fig. 4. At a fixed tem-
冉 冊
perature, the optical gain increases with the carrier concen-
me␣2c 兩M a兩2 n2 p ET − Eg
Raug = I0 exp − , 共13兲 tration as the population inversion increases and the peak
4 ប 共1 + 2兲 Nc N
4 3 3/2 2
k BT shifts to higher photon energy as higher energy states in con-
duction and valence bands are occupied 关Fig. 4共a兲兴. For the
where under the assumption of parabolic bands same carrier concentration, greater optical gain is always ob-
[This article is copyrighted as indicated in the article. Reuse of AIP content is subject to the terms at: http://scitation.aip.org/termsconditions. Downloaded to ] IP:
137.132.123.69 On: Fri, 13 Mar 2015 04:36:52
033107-5 Sun, Soref, and Cheng J. Appl. Phys. 108, 033107 共2010兲
FIG. 5. 共Color online兲 Maximum optical gain as a function of the nominal FIG. 6. 共Color online兲 Nominal current density at threshold as a function of
current density for the temperature range from 100 to 300 K. temperature for different losses.
tained at lower temperatures as shown in Fig. 4共b兲. require a lower density of carriers to achieve optical gain
The carriers that are responsible for the optical gain are because of their reduced density of states and their simulta-
excited by external pumping. A preferred pumping method is neous suppression of the Auger process. The recent result by
current injection capable of sustaining the carriers character- expitaxial technique exhibiting very fine control of GeSiSn
ized by their lifetime. It is convenient to use the nominal alloy layers has indeed opened the pathway to developing
current density Jn to describe the pumping current in the GeSn QW lasers. There have already been proposals for la-
DHS, which is defined as current per unit volume given in sers based on strained GeSn/GeSiSn QW structures.5,6 By
the unit of ampere per centimeter square per micron, and reducing the active layer thickness to ⬃100 Å in our DHS,
when multiplied by the active region thickness, yields a we plan to investigate lattice-matched strain free GeSn/
pumping current density. At lasing threshold where the GeSnSi multiple QW lasers. These encouraging efforts will
stimulated emission is just about to take place, the nominal potentially lead to room temperature electrically injected la-
current density is related to carrier recombination rates when ser diodes that can be monolithically integrated on Si or SOI.
radiative and Auger processes are dominant The laser simulated here emits at a photon energy that is
Jn = e关Rrad共n兲 + Raug共n兲兴. 共18兲 higher than the ⌫ band gap energy of the active bulk layer.
Thus this DHS laser is really capable of operating anywhere
Since both the optical gain and the recombination rates de- within the 1.7 to 3.0 m mid IR band because the emission
pend on carrier density n, we can now relate the optical gain wavelength desired is chosen simply by adjusting the Sn
with the nominal current density. Figure 5 shows the maxi- content of the active layer to a value in the 6%–12% range as
mum optical gain obtained from the gain profile in Fig. 4 as seen in Fig. 1共a兲. Then the longer-wave laser is designed by:
a function of the nominal current density for a temperature 共1兲 taking the GeSn buffer lattice to match the new active
range from 100 to 300 K. GeSn lattice, 共2兲 choosing the ternary cladding composition
In order to achieve lasing, the optical gain must be high to match the buffer lattice and to provide type-I offsets to the
enough to compensate for the loss. There are various loss active region, and 共3兲 calculating as before the carrier life-
mechanisms such as free-carrier absorption, waveguide loss, time, gain and threshold properties versus temperature. It is
and mirror loss that contribute to the required infrared gain. clear that the longer wave DHS would, once again, require
We have examined the nominal current density necessary to cooling in the 100 to 200 K regime.
generate sufficient optical gain to overcome a range of losses Looking to the future when a multiple-QW structure is
共20– 100 cm−1兲 under different temperatures as shown in substituted for the present three layer structure, we envision
Fig. 6. that room-temperature CW operation will be feasible for the
At T = 300 K, for an active region of 0.1 m, the cur- new GeSn QW laser diode at shorter-than-DHS wavelengths
rent density that is required to compensate for loss of of operation—an emission that could be shifted towards the
20 cm−1 exceeds 10 KA/ cm2, and needs to be much 1.55 m telecom range if desired because of the blueshift in
higher if more loss exists. While it is always desirable to subband-to-subband energy induced by the quantum confine-
have lasers operate at or above room temperature, such a ment.
large pumping current makes it rather unlikely with the
Ge0.94Sn0.06 / Ge0.75Si0.15Sn0.1 DHS. Looking at the estimates
V. CONCLUSION
in Fig. 6, if we pick a reasonable “ceiling” for current den-
sity, we conclude that the GeSn DHS laser requires cryo- We propose a simple group-IV laser made of
genic cooling in the 100–200 K range. Ge0.94Sn0.06 / Ge0.75Si0.15Sn0.1 DHS where the direct
We have found that the strong temperature dependence Ge0.94Sn0.06 active region is confined by two Ge0.75Si0.15Sn0.1
of threshold current in bulk heterostructure long wave lasers cladding layers. The compositions of both active and clad-
is due to the high density of carriers and the Auger recom- ding layers are determined to yield lattice matching and
bination process. One possible solution to the room- type-I band alignment between them. The entire DHS would
temperature problem is to work with GeSn QW lasers that be grown on a relaxed Ge0.94Sn0.06 buffer on either Si or SOI,
[This article is copyrighted as indicated in the article. Reuse of AIP content is subject to the terms at: http://scitation.aip.org/termsconditions. Downloaded to ] IP:
137.132.123.69 On: Fri, 13 Mar 2015 04:36:52
033107-6 Sun, Soref, and Cheng J. Appl. Phys. 108, 033107 共2010兲
2
hence the laser structure is strain free. Both radiative and R. A. Soref, J. Kouvetakis and J. Menendez, Mater. Res. Soc. Symp. Proc.,
958, L01 共2007兲.
Auger processes are included in the simulation of carrier 3
J. Liu, X. Sun, R. Camacho-Aguilera, L. C. Kimerling, and J. Michel, Opt.
recombination rate. A relationship between the optical gain Lett. 35, 679 共2010兲.
4
and the nominal current density that injects carriers into the J. Tolle, R. Roucka, V. D’Costa, J. Menendez, A. Chizmeshya, and J.
active region is established for a wide range of operating Kouvetakis, Mater. Res. Soc. Symp. Proc. 891, 579 共2006兲.
5
G.-E. Chang, S.-W. Chang, and S. L. Shuang, Opt. Express 17, 11246
temperature. Assuming different losses, we have determined 共2009兲.
the threshold current density as a function of operating tem- 6
Y.-H. Zhu, Q. Xu, W.-J. Fan, and J.-W. Wang, J. Appl. Phys. 107, 073108
perature. The result indicates that the current density that is 共2010兲.
7
needed to compensate for the loss at room-temperature op- G. Sun, H. H. Cheng, J. Menendez, J. B. Khurgin, and R. A. Soref, Appl.
Phys. Lett. 90, 251105 共2007兲.
eration is too high, and cooling is necessary. However, the 8
R. A. Soref and C. H. Perry, J. Appl. Phys. 69, 539 共1991兲.
proposed DHS laser simplifies the device growth which 9
H. P. L. de Guevara, A. G. Rodríguez, H. Navarro-Contreras, and M. A.
could otherwise be much more complex for lasers based on Vidal, Appl. Phys. Lett. 91, 161909 共2007兲.
10
quantum structures. The implementation of this laser will M. Jaros, Phys. Rev. B 37, 7112 共1988兲.
11
S. Adachi, Properties of Group-IV, III-V, and II-VI Semiconductors 共Wiley,
lead to the first electrically injected group-IV laser that is England, 2005兲.
integrated on Si or SOI. 12
V. R. D’Costa, C. S. Cook, A. G. Birdwell, C. L. Littler, M. Canonico, S.
Zollner, J. Kouvetakis, and J. Menendez, Phys. Rev. B 73, 125207 共2006兲.
13
ACKNOWLEDGEMENTS V. R. D’Costa, C. S. Cook, J. Menéndez, J. Tolle, J. Kouvetakis, and S.
Zollner, Solid State Commun. 138, 309 共2006兲.
14
This work was supported in part by the Air Force Office J. Weber and M. I. Alonso, Phys. Rev. B 40, 5683 共1989兲.
15
M. L. Cohen and T. K. Bergstresser, Phys. Rev. 141, 789 共1966兲.
of Scientific Research, Dr. Gernot Pomrenke, Program Man- 16
G. P. Agrawal and N. K. Dutta, Long-Wavelength Semiconductor Lasers
ager. 共Van Nostrand Reinhold, New York, 1986兲.
17
Y.-Y. Fang, J. Xie, J. Tolle, R. Roucka, V. R. D’Costa, A. V. G.
1
V. R. D’Costa, Y.-Y. Fang, J. Tolle, J. Kouvetakis, and J. Menéndez, Thin Chizmeshya, J. Menendez, and J. Kouvetakis, J. Am. Chem. Soc. 130,
Solid Films 518, 2531 共2010兲. 16095 共2008兲.
[This article is copyrighted as indicated in the article. Reuse of AIP content is subject to the terms at: http://scitation.aip.org/termsconditions. Downloaded to ] IP:
137.132.123.69 On: Fri, 13 Mar 2015 04:36:52