檢視 $RandomizedR 的原始碼

$RandomizedR is a special command to setup the random alloy structures. The format is

 $RandomizedR
 Region_Number_used Random_type location_N_used seeding_num
 location_1 location_2 ...... location_N
 Region_1 method P1 P2 P3 P4 P5  .... 
 Region_2 method P1 P2 P3 P4 P5  .... 
 ...
 Region_N method  P1 P2 P3 P4 P5  .... 


 Method=0 fully random number
 Method=1 P1=dx(um), P2=average_composition, P3=sigma, P4=indium_low P5=delta- In%
 Method=2 P1=dx(um), P2=average_composition, P3=sigma, P4=indium_low P5=delta- In%
 Method=3 P1=dx(um), P2=average_composition, P3=sigma, P4=indium_low P5=delta- In%

If method ==11 12 13, the format would be 

 $RandomizedR
 Region_Number_used Random_type location_N_used seeding_num
 location_1 location_2 ...... location_N
 Region_1 method Region_ini  Region_fin P1 P2 P3 P4 P5  .... 
 Region_2 method Region_ini  Region_fin P1 P2 P3 P4 P5  .... 
 ...
 Region_N method Region_ini  Region_fin  P1 P2 P3 P4 P5  .... 

 Method=11 use_region_min use_region_max  P1=dx(um), P2=average_composition, P3=sigma, P4=indium_low P5=delta- In%
 Method=12 use_region_min use_region_max  P1=dx(um), P2=average_composition, P3=sigma, P4=indium_low P5=delta- In%
 Method=13 use_region_min use_region_max  P1=dx(um), P2=average_composition, P3=sigma, P4=indium_low P5=delta- In%

 The difference between 1(11) and 2(12) is that Method=2 (12) will not replace the trap-related coefficients, nonradiative lifetime, doping, impurity, and generation of the original region. 
 Method==2(12)  use the original trap-related coefficients, nonradiative lifetime, doping, impurity, and generation of the original region
 Method==3(13)  use the original polarization, trap-related coefficients, nonradiative lifetime, doping, impurity, and generation of the original region
 
 The idea is to generate virtual regions. The region number is "location_N_used" 
 The second line defines the ID of virtual regions (location_N_used). For example location_N_used=10
 13 14 15 16 17 18 19 20 21 22 23    means these regions are virtual region
 P1 is the lattice size. typically 2.83e-4 um. 
 P2=average_com of this QW . For example 0.15 is 15% average indium composition
 P3=sigma  -> broaden factor of averaging local indium composition. Typically 2-4 lattice size. Unit is um. 
 P4=indium_low  define indium composition of virtual region 1. 
 P5=delta-In%   defines the composition difference of each virtual region. It should be linear.
 
 Note that no matter which method is used, the mobility, B, and C will use the original region setting.


Example:


 $RandomizedR
 6 1 20 -10       #total 20 regions as virtual region
 13 14 15 16 17 18 19 20 21 22 23 .... 32     # The region id of 20 virtual region 
 3 0                                       （QW location)
 4 0
 5 0
 6 0
 7  1  2.83e-4  0.15  9.8e-4 0.01 0.02
 8  2  2.83e-4  0.15  9.8e-4 0.01 0.02         # use the original lifetime, doping, , impurity, traps, and generation of the original region.
 9  3  2.83e-4  0.15  9.8e-4 0.01 0.02         # use the original polarization, lifetime, doping, , impurity, traps, and generation of the original region.
 10 11  1  10   2.83e-4  0.15  9.8e-4 0.01 0.02   # use virtual region 1-10 as 1% to 19%  if composition smaller than 1%, it would be 1%. If > 19% it would be 19%
 10 11  11 20   2.83e-4  0.40  9.8e-4 0.21 0.02    # use virtual region 11-20 as 21% to 39%  if composition smaller than 21%, it would be 21%. If > 39% it would be 39%