$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}P{3}_{1}P{4}_1$
${\beta }_2{R}_{th,2}P{3}_{2}P{4}_2$
${\beta }_3{R}_{th,3}P{3}_{3}P{4}_3$
${\beta }_4{R}_{th,4}P{3}_{4}P{4}_4$
${\beta }_5{R}_{th,5}P{3}_{5}P{4}_5$
.....
..... total N layers (regions)
${\beta }_N{R}_{th,N}P{3}_{N}P{4}_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. It might 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/c{m}^{3}s)$