# Modeling Quantum Well Lasers

Authors: Dan Alexandru Anghel; Andreea Rodica Sterian; Paul E. Sterian
Journal: Mathematical Problems in Engineering
 (2012)

---

## Abstract

In semiconductor laser modeling, a good mathematical model gives near-reality results. Three methods of modeling solutions from the rate equations are presented and analyzed. A method based on the rate equations modeled in Simulink to describe quantum well lasers was presented. For different signal types like step function, saw tooth and sinus used as input, a good response of the used equations is obtained. Circuit model resulting from one of the rate equations models is presented and simulated in SPICE. Results show a good modeling behavior. Numerical simulation in MathCad gives satisfactory results for the study of the transitory and dynamic operation at small level of the injection current. The obtained numerical results show the specific limits of each model, according to theoretical analysis. Based on these results, software can be built that integrates circuit simulation and other modeling methods for quantum well lasers to have a tool that model and analysis these devices from all points of view.


---

## Body

## 1. Introduction

Laser semiconductor diodes are key components in modern optical communications, storage, printing, medicine, and information processing. These types of devices were developed and evolved constantly in the direction of size reduction and integration. A quantum well (QW) laser improves the functioning characteristics of laser diodes in the direction of low threshold current and narrow emission band as well as emitted wavelength dependence on nanostructure dimension (quantum size effect) [1–8].A quantum well laser is a structure in which the active region of the device is so narrow that quantum confinement occurs, according to quantum mechanics. The wavelength of the light emitted by a quantum well laser is determined by the width of the active region rather than just the band gap of the material from which the device is realized. Consequently, much shorter wavelengths can be obtained from quantum well lasers than from conventional laser diodes using a particular semiconductor material. The efficiency of a quantum well laser is also greater than a conventional laser diode due to the stepwise form of its density of states function. In the optoelectronic integrated circuits (OEICs), the quantum well lasers will interact optically and electrically with other devices as integrated modulators, optical fibers and coherent amplifiers, detectors, and optical integrated guides [9–19]. A fully analytical method to describe these kind of complex systems cannot be realistic without questionable simplifying assumptions so that the implementations of refined mathematical methods and algorithms as well as circuit simulators become essentially in the applications development [20–24].


## 2. General Modeling Methods and Techniques

This paper aims at presenting and analyzing modeling techniques for quantum well lasers starting from the rate equations. The issue of modeling of quantum well semiconductor lasers has been addressed by several authors [1–3, 7, 13]. In operation of a laser, effects are involved that can impact a variety of performance parameters of the device. In some cases is crucial to obtain a clear picture of how a laser works in various conditions or parameters. A good laser model which allows calculation and testing of certain parameters of the device is an important result. The main models of laser with quantum wells are based on the description of the semiconductor lasers using the rate equations formalism.While the first models were based on one pair of equations to describe the density of photons and carriers in the active region, recent approaches include additional rate equations, to take into account, and carriers transport between active region and adjacent layers of the structure.Note that in most cases the rate equations lead to multiple solutions, although only one solution is correct. Javro and Kang [25] showed that incorrect solutions or without physical sense can be eliminated or avoided through a change of variables in the rate equations.However, the transformations used are available under certain conditions and for some cases give unrealistic solutions. These shortcomings are caused mainly due to the linear character of the gain-saturation coefficient. More general expression of the gain-saturation coefficient, proposed by Channin can be used to obtain models for operation regimes having a unique solution. Agrawal suggests another expression for this coefficient, which is also suitable. As it is shown, any of these two forms of the gain-saturation coefficient can be used to obtain models with a solution unique of the operation regime.In this paper, based on known rate equations models we try, firstly, the circuit simulation for two different modeling systems: Simulink and SPICE. The Simulink modeling technique is very simple. It is based on standard rate equations with the sets of parameters given directly in Simulink. The second model is based on the standard rate equations that use a nonlinear gain-saturation term proposed by Channin. In this second model we also describe the circuit level in SPICE with simulation results.A numerical experiment in MathCad which gives satisfactory results for the study of the transitory and dynamic operation at small level of the injection signal is shortly presented, also.


## 3. The Model with Linear Gain Saturation, in Simulink

One of the prevailing laser diode models is based on a set of rate equations. The rigorous derivation of these equations originates from Maxwell equations with a quantum mechanical approach for the induced polarization. However, the rate equation could also be derived by considering physical phenomena [1, 3, 5]:
(3.1)dNdt=IqVact-g0(N-N0)(1-εS)S-Nτp+Neτn,(3.2)dSdt=Γg0(N-N0)(1-εS)S+ΓβNτn-Sτp,(3.3)SPf=Γτpλ0Vactηhc=ϑ.
Equation (3.1) relates the rate of change in carrier concentration N to the injection current I, the carrier recombination rate and the stimulated emission rate. Equation (3.2) relates the rate of change in photon density S to the photons loss, the rate of coupled recombination into the lasing mode, and the stimulated emission rate. The ratio between photon density S to the output power Pf is described in (3.3). The other parameters used have well-known significance [2]. This simple model can be directly implemented with Simulink without any problems.A direct implementation looks like that in Figure1.Figure 1
Construction of rate equation block.Figures1 and 2 show how the rate equations model is constructed in Simulink. with I as input parameter from a signal generator and S, N, and Pf as output parameters. All the parameters in the rate equations can be modified before the simulation starts.Figure 2
Rate equation in Simulink.For the simulation we used the following parameters found in specialized papers [1–5]: λ0=1.502×10-4 cm (laser wavelength), Vact=9×10-11 cm3 (active region volume), Γ=0.44 (optical confinement factor), B=4×10-4 (spontaneous emission factor), g0=3×10-6 cm3/s (gain coefficient), N0=1.2×1018 cm−3 (optical transparency density), τn=3 ns (carrier lifetime), τp=1 ps (photon lifetime), η=0.1 (quantum efficiency), Ne=5.41×1010 cm−3 (equilibrium carrier density), and ε=3.4×10-17 cm3 (gain saturation factor).With the above model different signal forms can be used as input for the quantum well laser. They show that theoretical response of the equations is good in comparison with real results that are obtained in applications. Signals like step function, saw tooth, and sine types are used as input. The results are shown in Figures3 and 4, which illustrate a very fast response at low current level. Low threshold current is the main feature of quantum well lasers, and it is directly shown for the basic form of rate equations. The simulation is not perfect, and this is because of negative solutions for N and S and high power solution of the equations.Figure 3
Output for step input signal of 10 mA amplitude and frequency of 5 MHz.Output for saw tooth input signal of 10 mA amplitude and frequency of 5 MHz (a) and for sinusoidal input signal (b).
(a)(b)


## 4. The Model with Nonlinear Gain Saturation, in SPICE

To obtain nonnegative solutions for the photon density and the density of carriers in terms of a nonnegative current injection, a different version of the standard rate equations must be formulated. In the corresponding generalized equations (3.1) and (3.2), the linear gain-saturation term is replaced by a nonlinear term proposed by Channin or Agrawal as follows [1–3]:
(4.1)dNdt=ηiIqNwVact-Rw(N)-Γcνgrα(N)Φ(S)S,(4.2)dSdt=-Sτp+NwRwβ(N)+NwΓcνgrα(N)Φ(S)S,(4.3)SPf=λτpηcVacthc=ϑ.
Equation (4.1) describes the carrier concentration rate dependence on the injection current, the carrier recombination rate Rw(N) and the stimulated emission rate. To take into account all mechanisms of recombination, we consider: Rw(N)=AN+BN2+CN3, where A, B, and C are, respectively, the unimolecular, radiative, and Auger recombination coefficients.Equation (4.2) shows the photon density rate dependence on photon loss, the rate of coupled recombination into the lasing mode, and the stimulated emission rate. Equation (4.3) relates the output power Pf to the photon density S. The other parameters used have well-known significance [2]. In (4.1) and (4.2), the correlation between the material gain and the carrier density is given by the logarithmic carrier-dependent gain term:
(4.4)α(N)=G0ln⁡(Rw(N)Rw(N0)),
where G0 is the one quantum well gain coefficient and N0 is the optical transparency density while the gain-saturation function can take on one of the following two forms:
(4.5)ϕ-1(S)=11+εΓcS(Channin)orϕ-1(S)=11+εΓcS(Agrawal).
An equivalent circuit model based on (4.1)–(4.3) can be implemented in SPICE as in Figure 5. Unlike models based on the rate equations that use a linear gain-saturation term, this circuit model is applicable for all nonnegative values of injection current.Figure 5
The circuit implementation of the model.Figure6 shows the circuit implementation of the model, using suitable variables transformations [25]: N=Neexp(qV/nkT) and Pf=(m+δ)2, which impose nonnegative solutions for N and Pf. The linear recombination and charge storage in the device are described by diodes D1 and D2 and current sources Ic1 and Ic2. The effects of additional recombination mechanisms and stimulated emission, respectively, on the carrier density are modeled by Br1 and Bs1, according to (4.11) and (4.13).Answer of optical power with variation of the input bias current between 0 and 25 mA (a) and between 0 and 50 mA (b).
(a)(b)The effects of spontaneous and stimulated emission, respectively, are accounted in the circuit byBr2 and Bs2 while Rph and Cph complets the model for the time variation of the photon density. The voltage Bpf describe the optical output power of the laser.The equations for the elements of the circuit in Figure8 according to [2] are as follows: I=IT1+ID1+IC1+Br1+Bs1, where
(4.6)IT1=ID1+IC1,(4.7)2τpdmdt+m=Br2+Bs2,Bpf=(m+δ)2,(4.8)ID1=qNwVactNe2ηiτn[exp⁡(qVnkT)-1],(4.9)ID2=qNwVactNe2ηiτn[exp⁡(qVnkT)-1+2qτpnkTexp⁡(qVnkT)dVdt],(4.10)Ic1=Ic2=qNwVactNe2ηiτn,(4.11)Br1=qNwVactηiRω2(ΘIT1),(4.12)Bs1=λτpNwΓcυgrηiηchcα(ΘIT1)ϕ(ϑ(m+δ)2)ϑ(m+δ)2,(4.13)Br2=NwηcVacthcλϑ(m+δ)Rω2(ΘIT1),(4.14)Bs2=τpNwΓcυgrα(ΘIT1)ϕ(ϑ(m+δ)2)(m+δ)-δ.
The SPICE netlist was implemented in AIM-SPICE. Calculating the parameters in the netlist is time consuming, and if a calculation is wrong, the netlist will fail in SPICE or the result will not be good. Figures 6, 7, and 8 resulting from simulation with simple DC sweep and transitory sweep illustrate the following:(i)
answer of optical power with variation of the input bias current between 0 and 25 mA (a) and between 0 and 50 mA (b) (Figure6);(ii)
transient output power in response between 0 and 100 ps (Figure7);(iii)
transient output power in response between 0 and 5 ns (Figure8).Figure 7
Transient output power in response between 0 and 100 ps.Figure 8
Transient output power in response between 0 and 5 ns. Note: Output power is inV because SPICE cannot output variables in W.


## 5. Numerical Experiment in MathCad

Study of functional characteristics of QWL is possible based on the models described previosly, by implementing MathCad programs to integrate these equations using appropriate algorithms corresponding to different types of pumping signals.The model given by the equations with linear gain-saturation term (3.1)–(3.3) was integrated for a constant current injection Ii=0.05 A. The waveforms corresponding to Ii and to output optical power P are obtained [13].The model given by the equations with nonlinear gain saturation (4.1)–(4.3) was integrated for a constant current injection Ii(t)=0.025 A (Figure 9), sinus current injection (Figure 10):(5.1)Ii(t)=1+sin⁡(5⋅108t),
and rectangular current injection (Figure 11), which for several periods can be written as:
(5.2)Ii(t)=10-8+0.01+0.25Ii(t-2⋅10-9)+0.25Ii(t-7⋅10-9)+0.25Ii(t-12⋅10-9).
Simulations results for these large-signal models have shown excellent agreement with experimental data. A corresponding small-signal model has been derived for the quantum well lasers, to study the modulation properties in the frequency domain [13].Figure 9
The waveforms corresponding toIi and to output optical power P.Figure 10
The waveforms corresponding toIi and to output optical power P.Figure 11
The waveforms corresponding toIi and to output optical power P.


## 6. Conclusions

(1)
In the paper, some modeling methods and techniques for the quantum well lasers have been implemented and validated. A suitable modeling technique is very important for technology designing and analyzing of the integrated optoelectronic circuits.(2)
A method based on the rate equations for using Simulink to describe quantum well lasers was presented. For different signal types like step function, saw tooth, and sinus used as input, the results show a good theoretical response of the equations in comparison with the real results obtained in applications. This method is useful to determine the various regions of operation of the laser diode also.(3)
The SPICE simulation of the circuits shows that modifying one parameter will result in new calculations and new SPICE netlist, which is a rather complicated procedure. Simulation is not limited to SPICE; any all-purpose circuit simulator can be used to get similar results.(4)
Generally, all studied models give satisfactory results for the study of the transitory and dynamic operation at small level of the injection signals. At large injection signal levels the obtained numerical results show the specific limits of each model, according to theoretical analysis.(5)
As a future development, software can be built that integrates circuit simulations and other modeling methods for quantum well lasers to have a tool that models these devices from all points of view.



---

*Source: 736529-2012-03-07.xml*