Modeling Quantum Well Lasers

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 goodmodeling behavior. Numerical simulation inMathCad 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 eachmodel, 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.


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 .

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.approach for the induced polarization.However, the rate equation could also be derived by considering physical phenomena 1, 3, 5 : 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 P f 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 Figure 1.Figures 1 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 P f as output parameters.All the parameters in the rate equations can be modified before the simulation starts.
For the simulation we used the following parameters found in specialized papers 1-5 : λ 0 1.502 × 10 −4 cm laser wavelength , V act 9 × 10 −11 cm 3 active region volume , 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 Figures 3 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.

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 gainsaturation term is replaced by a nonlinear term proposed by Channin or Agrawal as follows 1-3 : Equation 4.1 describes the carrier concentration rate dependence on the injection current, the carrier recombination rate R w N and the stimulated emission rate.To take into account all mechanisms of recombination, we consider: R w N AN BN 2 CN 3 , 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 P f to the photon density S. The other parameters used have wellknown 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: where G 0 is the one quantum well gain coefficient and N 0 is the optical transparency density while the gain-saturation function can take on one of the following two forms: 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 6 shows the circuit implementation of the model, using suitable variables transformations 25 : N N e exp qV/nkT and P f m δ 2 , which impose nonnegative solutions for N and P f .The linear recombination and charge storage in the device are described by diodes D1 and D2 and current sources I c1 and I c2 .The effects of additional recombination mechanisms and stimulated emission, respectively, on the carrier density are modeled by B r1 and B s1 , according to 4.11 and 4.13 .
The effects of spontaneous and stimulated emission, respectively, are accounted in the circuit by B r2 and B s2 while R ph and C ph complets the model for the time variation of the photon density.The voltage B pf describe the optical output power of the laser.
The equations for the elements of the circuit in Figure 8 according to 2 are as follows:  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 Figure 6 ; ii transient output power in response between 0 and 100 ps Figure 7 ; iii transient output power in response between 0 and 5 ns Figure 8 .

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 I i 0.05 A. The waveforms corresponding to I i 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 I i t 0.025 A Figure 9 , sinus current injection Figure 10 : and rectangular current injection Figure 11 , which for several periods can be written as: 5.2   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 .

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.

Figure 1 :
Figure 1: Construction of rate equation block.

Figure 5 :
Figure 5: The circuit implementation of the model.

Figure 6 :
Figure 6: Answer of optical power with variation of the input bias current between 0 and 25 mA a and between 0 and 50 mA b .

Figure 7 :
Figure 7: Transient output power in response between 0 and 100 ps.

Figure 8 :
Figure 8: Transient output power in response between 0 and 5 ns.Note: Output power is in V because SPICE cannot output variables in W.

Figure 9 :
Figure 9: The waveforms corresponding to I i and to output optical power P .

Figure 10 :Figure 11 :
Figure 10: The waveforms corresponding to I i and to output optical power P .