Aide ATLAS SILVACO

Voila le programme

go atlas
set cellWidth=5.000000e+002
set capWidthpercent=8.000000e+000
set divs=1.000000e+001
set contThick=1.000000e-001
set capThick=3.000000e-001
set capDop=1.000000e+020
set windowThick=0.01
set winDop=2.15e17
set emitterThick=0.01
#changed emitDop from 1e16 to 1e20
set emitDop=1e16
set baseThick=3.19467
#changed basDop from 1e16 to 1e20
set baseDop=1e16
set bsfThick=0.03533
set bsfDop=2.15e19
set cellWidthDiv=$cellWidth/$divs
set width3d=100e6/$cellWidth
set capWidth=0.01*$capWidthpercent*$cellWidth/2
set capWidthDiv=$capWidth/($divs/2)
set cellWidthHalf=$cellWidth/2
set bsfLo=0
set bsfHi=$bsfLo-$bsfThick
set bsfDiv=$bsfThick/$divs
set baseLo=$bsfHi
set baseHi=$baseLo-$baseThick
set baseMid=$baseLo-$baseThick/2
set baseDiv=$baseThick/$divs
set emitterLo=$baseHi
set emitterHi=$emitterLo-$emitterThick
set emitterDiv=$emitterThick/$divs
set windowLo=$emitterHi
set windowHi=$windowLo-$windowThick
set windowDiv=$windowThick/$divs
set capLo=$windowHi
set capHi=$capLo-$capThick
#set capDiv=$capThick/$divs
set contLo=$capHi
set contHi=$contLo-$contThick
set contDiv=$contThick/$divs
set lightY=$emitterHi-5
mesh width=$width3d
## X-Mesh
x.mesh loc=-$cellWidthHalf spac=$cellWidthDiv
x.mesh loc=-$capWidth spac=$capWidthDiv
x.mesh loc=$capWidth spac=$capWidthDiv
x.mesh loc=$cellWidthHalf spac=$cellWidthDiv
## Y-Mesh
# Top contact
y.mesh loc=$contHi spac=0
y.mesh loc=$contLo spac=0
# Cap
# Window
y.mesh loc=$windowHi spac=$windowDiv
y.mesh loc=$windowLo spac=$windowDiv
# Emitter
y.mesh loc=$emitterLo spac=$emitterDiv
# Base
y.mesh loc=$baseMid spac=$baseDiv
# BSF
y.mesh loc=$bsfHi spac=$bsfDiv
y.mesh loc=$bsfLo spac=$bsfDiv
## Regions
# Cap
region num=8 material=Vacuum x.min=-$capWidth x.max=$capWidth y.min=$contHi y.max=$contLo
region num=1 material=InGaN x.min=-$capWidth x.max=$capWidth y.min=$capHi y.max=$capLo x.comp=0.20
region num=2 material=Vacuum x.min=-$cellWidthHalf x.max=-$capWidth y.min=$contHi y.max=$capLo
region num=3 material=Vacuum x.min=$capWidth x.max=$cellWidthHalf y.min=$contHi y.max=$capLo
# Window [for Ge cell, use AlGaAs with x.comp=0.2]
region num=4 material=InGaN x.min=-$cellWidthHalf x.max=$cellWidthHalf y.min=$windowHi y.max=$windowLo x.comp=0.20
region num=4 material=AlGaAs x.comp=0.2 x.min=-$cellWidthHalf x.max=$cellWidthHalf y.min=$windowHi y.max=$windowLo
# Emitter
region num=5 material=InGaN x.min=-$cellWidthHalf x.max=$cellWidthHalf y.min=$emitterHi y.max=$emitterLo x.comp=0.20
# Base
region num=6 material=InGaN x.min=-$cellWidthHalf x.max=$cellWidthHalf y.min=$baseHi y.max=$baseLo x.comp=0.20
# BSF
region num=7 material=InGaN x.min=-$cellWidthHalf x.max=$cellWidthHalf y.min=$bsfHi y.max=$bsfLo x.comp=0.20
## Electrodes [for InGaP cell, add cathode (gold) and remove cathode (conductor)]
electrode name=cathode material=Gold x.min=-$capWidth x.max=$capWidth y.min=$contHi y.max=$contLo
#electrode name=cathode x.min=-$cellWidthHalf x.max=$cellWidthHalf y.min=$windowHi y.max=$windowHi
electrode name=anode x.min=-$cellWidthHalf x.max=$cellWidthHalf y.min=$bsfLo y.max=$bsfLo
## Doping [for InGaP cell, uncomment cap doping]
# Cap
doping uniform region=1 n.type conc=$capDop
# Window
doping uniform region=4 n.type conc=$winDop
# Emitter
doping uniform region=5 n.type conc=$emitDop
# Base
doping uniform region=6 p.type conc=$baseDop
# BSF
doping uniform region=7 p.type conc=$bsfDop

## Material properties
# Opaque contact [comment out for InGaP cell]
material region=8 real.index=1.2 imag.index=1.8
# Vacuum (for zero reflection) [change to match window material (InGaP use Vacuum_AlInP)]
# [for InGaP cell, comment out region 1]
#material region=1 index.file=Vacuum_InGaP.opt
#material region=2 index.file=VacuumIn20Ga80N.opt
#material region=3 index.file=VacuumIn20Ga80N.opt
#InGaN
#material material=InGaN EG300=2.6612 index.file=In20Ga80N.opt
# Gold
material material=Gold real.index=1.2 imag.index=1.8
## Models [InGaP cell, 1; GaAs cell, 5&6; InGaNAs cell, 7]
models region=1 CONMOB
MODELS CHUANG CONMOB FLDMOB SRH OPTR PRINT
## Light beams [GaAs b1,0.55-0.75,200 b2,0.75-0.88,65] 0.12-2.7,50 [630,825]
beam num=1 x.origin=0 y.origin=$lightY angle=90 back.refl \
wavel.start=0.12 wavel.end=2.4 wavel.num=50
struct outfile=SingleCell_webf.str
#tonyplot SingleCell_webf.str
solve init
method gummel newton maxtraps=10 itlimit=25
solve b1=0.9
## Getting Isc for I-V curve points
method newton maxtraps=10 itlimit=100
solve b1=0.95
extract name="isc" max(i."cathode")
set isc=$isc*$width3d
set i1=$isc/10
set i2=$i1+$isc/10
set i3=$i2+$isc/10
set i4=$i3+$isc/10
set i5=$i4+$isc/10
set i6=$i5+$isc/20
set i7=$i6+$isc/20
set i8=$i7+$isc/20
set i9=$i8+$isc/20
set i10=$i9+$isc/20
set i11=$i10+$isc/40
set i12=$i11+$isc/40
set i13=$i12+$isc/40
set i14=$i13+$isc/40
set i15=$i14+$isc/40
set i16=$i15+$isc/80
set i17=$i16+$isc/80
set i18=$i17+$isc/80
set i19=$i18+$isc/80
set i20=$i19+$isc/80
set i21=$i20+$isc/80
set i22=$i21+$isc/80
set i23=$i22+$isc/80
set i24=$i23+$isc/80
set i25=$i24+$isc/80-0.00001
##
log outfile=In20Ga80N.log
method newton maxtraps=10 itlimit=100
solve b1=0.95
contact name=anode current
method newton maxtraps=10 itlimit=100
## Pmax points [InGaP 18-25; GaAs 15-25; InGaNAs 13-25; Ge 11-25]
solve ianode=-$i25 b1=0.95
solve ianode=-$i24 b1=0.95
solve ianode=-$i23 b1=0.95
solve ianode=-$i22 b1=0.95
solve ianode=-$i21 b1=0.95
solve ianode=-$i20 b1=0.95
solve ianode=-$i19 b1=0.95
solve ianode=-$i18 b1=0.95
solve ianode=-$i17 b1=0.95
solve ianode=-$i16 b1=0.95
solve ianode=-$i15 b1=0.95
solve ianode=-$i14 b1=0.95
solve ianode=-$i13 b1=0.95
solve ianode=-$i12 b1=0.95
solve ianode=-$i11 b1=0.95
solve ianode=-$i10 b1=0.95
solve ianode=-$i9 b1=0.95
solve ianode=-$i8 b1=0.95
solve ianode=-$i7 b1=0.95
solve ianode=-$i6 b1=0.95
solve ianode=-$i5 b1=0.95
solve ianode=-$i4 b1=0.95
solve ianode=-$i3 b1=0.95
solve ianode=-$i2 b1=0.95
solve ianode=-$i1 b1=0.95
solve ianode=0 b1=0.95
log off
extract name="iv" curve(v."anode", i."cathode") outfile="IVcurveIn20Ga80N.dat"
tonyplot IVcurveIn20Ga80N.dat
log outfil