$activateStimEL

$activateStimEL is a command used to enable the calculation of stimulate Emission.

The stimulate emission coefficient is calculated by:

$B_{stim} = B_{0} \times \beta \times \frac{R_{rad}}{R_{th}-R_{rad}}$
$R_{stim} = B_{stim} \times n \times p$
$B_{0}:$ is the radiative recombination coefficient
$R_{th}:$ is the threshold recombination rate
$R_{rad}:$ is the spontaneous radiative recombination rates decided by  $B_{0}\times n \times p$
$\beta:$ is the fitting coefficient. 

The way to activate the command is as following

$activateStimEL
Method QW_number 1st_QW_layerNum 2nd_QW_layerNum 3rd_QW_layerNum ...
$\beta_{1}  R_{th,1}$ , P3_1, P4_1
$\beta_{2}  R_{th,2}$ , P3_2, P4_2
$\beta_{3}  R_{th,3}$ , P3_3, P4_3
$\beta_{4}  R_{th,4}$ , P3_4, P4_4
$\beta_{5}  R_{th,5}$ , P3_5, P4_5
.....
..... total N layers (regions)
$\beta_{N}  R_{th,N}$ , P3_N, P4_N

Method: 1 -> Given $\beta and R_{th}$ of each layer. If this is layer is not QW, set $\beta=0.0~~,and~~R_{th}=0.0$
Method: 2 -> defined for future expansion
P3, P4 : Not used, may be used when Method 2 is activate in the future

For example:

$activateStimEL
1 1 2 
0.000 5.0e27 0.00 0.00
1.000d-4 5.0e27 0.00 0.00
0.000 5.0e27 0.00 0.00
0.000 5.0e27 0.00 0.00

In this command, use method 1. There is 1 QW which is the 2nd layer. The $\beta$ and $R_{th}$ of the 2nd layer is 1.0e-4 and 5.0e27 $(1/cm^{3}s)$