• Photonics Research
  • Vol. 9, Issue 7, 1234 (2021)
Bahareh Marzban1、2, Daniela Stange2、3, Denis Rainko2、3, Zoran Ikonic4, Dan Buca2、3, and Jeremy Witzens1、2、*
Author Affiliations
  • 1Institute of Integrated Photonics, RWTH Aachen University, 52074 Aachen, Germany
  • 2Jülich-Aachen Research Alliance (JARA), Fundamentals of Future Information Technologies, Germany
  • 3Peter Grünberg Institute, 52428 Jülich, Germany
  • 4University of Leeds, School of Electronic and Electrical Engineering, Woodhouse, Leeds LS2 9JT, UK
  • show less
    DOI: 10.1364/PRJ.416505 Cite this Article Set citation alerts
    Bahareh Marzban, Daniela Stange, Denis Rainko, Zoran Ikonic, Dan Buca, Jeremy Witzens. Modeling of a SiGeSn quantum well laser[J]. Photonics Research, 2021, 9(7): 1234 Copy Citation Text show less
    Modeled microdisk laser. (a) Schematic cross-section of the disk with labeled dimensions. (b) Strain across the layer stack, along a vertical axis near the periphery of the disk, where the optical field of whispering gallery ground modes is maximal. (c) Measured thresholds reported by Stange et al. [29] at a 1064 nm pump wavelength (Nd:YAG).
    Fig. 1. Modeled microdisk laser. (a) Schematic cross-section of the disk with labeled dimensions. (b) Strain across the layer stack, along a vertical axis near the periphery of the disk, where the optical field of whispering gallery ground modes is maximal. (c) Measured thresholds reported by Stange et al. [29] at a 1064 nm pump wavelength (Nd:YAG).
    Flat-band energy diagram at 20 K. (a) Conduction band energy levels, with quantization levels in the wells (the gray curve shows the non-quantized EΓ band edge energy, the blue curves the quantized EΓ1 and EΓ2 levels, and the red curve EL1 inside the wells and the non-quantized EL inside the barriers). (b) Valence band energy levels, with quantization levels in the wells (the gray curve shows the non-quantized EHH band edge, the blue curves the EHH1, EHH2, and EHH3 levels, and the red curve ELH1 inside the wells and the non-quantized ELH inside the barriers).
    Fig. 2. Flat-band energy diagram at 20 K. (a) Conduction band energy levels, with quantization levels in the wells (the gray curve shows the non-quantized EΓ band edge energy, the blue curves the quantized EΓ1 and EΓ2 levels, and the red curve EL1 inside the wells and the non-quantized EL inside the barriers). (b) Valence band energy levels, with quantization levels in the wells (the gray curve shows the non-quantized EHH band edge, the blue curves the EHH1, EHH2, and EHH3 levels, and the red curve ELH1 inside the wells and the non-quantized ELH inside the barriers).
    Net gain as a function of pump power density. The experimentally determined thresholds are indicated by the vertical lines. (a) Lifetimes equal to 217 ps and 80 ps at 50 K and 100 K, respectively. Positive net gains equal to 267 cm−1 and 342 cm−1 are observed at the experimental thresholds at 50 K and 100 K, respectively. (b) Directness ΔEL−Γ reduced by 30 meV and lifetimes uniformly set to 80 ps. The net gain at threshold is close to zero for both temperatures and rolls over at 100 K. Pump power densities indicated on the x axes can be converted into carrier concentrations averaged across the thin film stack by taking the assumed lifetimes and stack thickness (690 nm) into account. QW dependent carrier concentrations at threshold are reported in Fig. 6.
    Fig. 3. Net gain as a function of pump power density. The experimentally determined thresholds are indicated by the vertical lines. (a) Lifetimes equal to 217 ps and 80 ps at 50 K and 100 K, respectively. Positive net gains equal to 267cm1 and 342cm1 are observed at the experimental thresholds at 50 K and 100 K, respectively. (b) Directness ΔELΓ reduced by 30 meV and lifetimes uniformly set to 80 ps. The net gain at threshold is close to zero for both temperatures and rolls over at 100 K. Pump power densities indicated on the x axes can be converted into carrier concentrations averaged across the thin film stack by taking the assumed lifetimes and stack thickness (690 nm) into account. QW dependent carrier concentrations at threshold are reported in Fig. 6.
    (a) Field profiles for the TE3, TM0, and TM3 modes, with Er, Ez, and Ep the E-field components along the radial, z, and azimuthal directions. While TM3 is a quasi-TM mode, it also has strong in-plane field components Er and Ep away from the disk’s symmetry plane. This explains why it is the quasi-TM mode with the highest gain at low temperatures (50 K), at which the TE material net gain is the highest. At higher temperatures (100 K), the TM material net gain is highest, and TM0 becomes the dominant lasing mode. (b) TE and TM confinement factors ΓTE and ΓTM as a function of QW number for each of these three modes. The pronounced off-symmetry-plane ΓTE for TM3 is readily apparent.
    Fig. 4. (a) Field profiles for the TE3, TM0, and TM3 modes, with Er, Ez, and Ep the E-field components along the radial, z, and azimuthal directions. While TM3 is a quasi-TM mode, it also has strong in-plane field components Er and Ep away from the disk’s symmetry plane. This explains why it is the quasi-TM mode with the highest gain at low temperatures (50 K), at which the TE material net gain is the highest. At higher temperatures (100 K), the TM material net gain is highest, and TM0 becomes the dominant lasing mode. (b) TE and TM confinement factors ΓTE and ΓTM as a function of QW number for each of these three modes. The pronounced off-symmetry-plane ΓTE for TM3 is readily apparent.
    Carrier concentrations versus vertical position for (a) 20 kW/cm2 pump intensity at 50 K and (d) 150 kW/cm2 pump intensity at 100 K along the axis indicated in Fig. 1. Barrier layers and buffer are shown as gray regions. Conduction and valence band energies at (b), (c) 50 K and (e), (f) 100 K, with the same color coding as in Fig. 2. In addition, quasi-Fermi levels are shown in green. To reduce the number of curves, the red curves in (b) and (e) show the EL band edge in the barriers, but the EL1 subband edge in the wells. Similarly, the purple curves in (c) and (f) show ELH in the barriers and ELH1 in the wells.
    Fig. 5. Carrier concentrations versus vertical position for (a) 20  kW/cm2 pump intensity at 50 K and (d) 150  kW/cm2 pump intensity at 100 K along the axis indicated in Fig. 1. Barrier layers and buffer are shown as gray regions. Conduction and valence band energies at (b), (c) 50 K and (e), (f) 100 K, with the same color coding as in Fig. 2. In addition, quasi-Fermi levels are shown in green. To reduce the number of curves, the red curves in (b) and (e) show the EL band edge in the barriers, but the EL1 subband edge in the wells. Similarly, the purple curves in (c) and (f) show ELH in the barriers and ELH1 in the wells.
    (a), (b) Electron and (d), (e) hole concentrations for each of the wells at (a), (d) 50 K and at the 20 kW/cm2 threshold pump intensity and (b), (e) 100 K and 150 kW/cm2. Electrons are classified as (total) Γ- and L-valley electrons, holes as (total) HHs and LHs. Confined Γ-valley electrons, HHs, and LHs are also separately shown. (c), (f) Carrier densities integrated over the entire stack and broken down in Γ- and L-valley electrons inside the wells, total electron density in the barriers, and total electron density in the buffer. Similarly, holes are broken down in HHs and LHs inside the wells, total holes in the barriers, and total holes in the buffer.
    Fig. 6. (a), (b) Electron and (d), (e) hole concentrations for each of the wells at (a), (d) 50 K and at the 20  kW/cm2 threshold pump intensity and (b), (e) 100 K and 150  kW/cm2. Electrons are classified as (total) Γ- and L-valley electrons, holes as (total) HHs and LHs. Confined Γ-valley electrons, HHs, and LHs are also separately shown. (c), (f) Carrier densities integrated over the entire stack and broken down in Γ- and L-valley electrons inside the wells, total electron density in the barriers, and total electron density in the buffer. Similarly, holes are broken down in HHs and LHs inside the wells, total holes in the barriers, and total holes in the buffer.
    Gain at 50 K and 100 K, respectively, at 20 kW/cm2 and 150 kW/cm2 threshold pump intensities. (a) Material gain for each well, for both polarizations. (b) Modal gain for each well. The TE3/TM3 modes have the highest gain at 50 K, and the TE3/TM0 modes at 100 K. (c) Total modal gain of the TE3 and TM0 modes, respectively, at 50 K and 100 K. (d) Breakdown of the modal gain for the TM0 mode at 100 K across quantum wells.
    Fig. 7. Gain at 50 K and 100 K, respectively, at 20  kW/cm2 and 150  kW/cm2 threshold pump intensities. (a) Material gain for each well, for both polarizations. (b) Modal gain for each well. The TE3/TM3 modes have the highest gain at 50 K, and the TE3/TM0 modes at 100 K. (c) Total modal gain of the TE3 and TM0 modes, respectively, at 50 K and 100 K. (d) Breakdown of the modal gain for the TM0 mode at 100 K across quantum wells.
    Calculated absorption losses at 50 K and 100 K and at respective lasing thresholds 20 kW/cm2 and 150 kW/cm2. (a) FCA as predicted by the modified Drude model and (b) IVBA, calculated for both TE and TM polarizations.
    Fig. 8. Calculated absorption losses at 50 K and 100 K and at respective lasing thresholds 20  kW/cm2 and 150  kW/cm2. (a) FCA as predicted by the modified Drude model and (b) IVBA, calculated for both TE and TM polarizations.
     Ge0.867Sn0.133(gain medium)Si0.048Ge0.822Sn0.130(barrier layers)Ge0.9Sn0.1(buffer layer)
    Γ-valley mobility at 5×1017cm3 total carrier concentration3500cm2/(Vs)1850cm2/(Vs)3190cm2/(Vs)
    L-valley mobility at 5×1017cm3 total carrier concentration775  cm2/(Vs)750cm2/(Vs)775cm2/(Vs)
    Average hole mobility at 5×1017cm3 total carrier concentration500cm2/(Vs)500cm2/(Vs)500cm2/(Vs)
    Carrier lifetime (final model)80 ps80 ps80 ps
    Background p-type doping2×1017cm32×1017cm32×1017cm3
    Absorption coefficient at 1064 nm5×104cm16×104cm14×104cm1
    Refractive index at 1064 nm4.54.54.5
    Refractive index at 2.5 μm4.24.24.1
    Luttinger parameterγ1=19.29γ1=18.96γ1=16.21
    γ2=7.36γ2=7.21γ2=5.84
    γ3=8.80γ3=8.64γ3=7.26
    Top surface recombination velocity10,000 cm/s
    Bottom surface recombination velocity
    Deformation potentials, bandgaps, and bandgap bowingTable 1 in Ref. [37]
    Table 1. Summary of Assumed Material Properties in the 50–100 K Range
    ParameterValue
    Gain broadening25 meV
    Bandgap narrowing at n=2×1017cm324 meV
    Bandgap narrowing at n=2×1018cm337 meV
    Correction to Drude model absorption, electrons2×
    Correction to Drude model absorption, holes2×
    Roughness correlation length (Lc)500 nm
    Roughness amplitude (σ)55 nm
    Table 2. Summary of Gain Modeling Assumptions
    Bahareh Marzban, Daniela Stange, Denis Rainko, Zoran Ikonic, Dan Buca, Jeremy Witzens. Modeling of a SiGeSn quantum well laser[J]. Photonics Research, 2021, 9(7): 1234
    Download Citation