$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/c{m}^{3}s)$