MaterialConstants: InAs_100nm.txt

File InAs_100nm.txt, 4.7 KB (added by kim568, 7 years ago)

Prophet codes for InAs

Line 
1#PROPHET script for cylindrical semiconductor InAs nanowire simulations
2# That's Joerg's session
3
4dbase create name=options/cylindrical sval="Y"
5dbase create name=library/physics/ambient/psi/scale rval=1
6dbase create name=library/physics/ambient/psi/SeeAlso sval=/library/physics/nitride/psi
7dbase create name=options/ignoreFPE ival=1
8dbase modify name=library/physics/InAs/eg.0 rval=0.415
9dbase modify name=library/physics/InAs/eg.alpha rval=0.000276
10dbase modify name=library/physics/InAs/eg.beta rval=83
11
12#------------added to avoid a library warning in web-based PROPHET-------------------
13dbase create name=library/physics/ambient/epsilon rval=1.0
14#------------------------------------------------------------------------------------
15
16dbase create name=options/ignoreFPE ival=1
17#define poisson equation for initial values#
18  system name=InAs_poisson
19    +       sysvars=psi
20    +       tmpvars=electrons,holes,eg,nc,nv,ni
21    +       term0=ndiv_fbm.lapflux(psi|psi)@{InAs,ambient,nitride,poly}
22    +       term1=nodal.nscd(electrons,holes,netdope|psi)@{InAs,poly}
23    +       term2=dirichlet.device_dirichlet(netdope,ni,nc,nv|psi)@{InAs/source,InAs/drain,poly/gate}
24    +       term3=constraint.continuity(psi|psi)@{InAs/nitride,poly/nitride}
25    +       func0=eg(tl,netdope|eg)@{InAs}
26    +       func1=10*dos(electrons,holes,tl|nc,nv)@{InAs}
27    +       func2=ni(tl,eg,nc,nv|ni)@{InAs}
28    +       func3=quasiFermi_fd(psi,electrons,ni,nc,tl|electrons)@{InAs}
29    +       func4=quasiFermi_fd(psi,holes,ni,nv,tl|holes)@{InAs}
30    +       init0="tl:${solve/bias/Tcelsius}+${library/physics/celsius}"
31   
32#define the coupled DD-Poisson equations#
33system name=InAs_dd
34    +       sysvars=psi,electrons,holes
35    +       tmpvars=eg,nc,nv,ni,nmob0,pmob0,nin,nip,taun,taup,phin,phip,ec,ev,evac
36    +       term0=dirichlet.device_dirichlet_h(netdope,ni,nc,nv|psi,electrons,holes)@{InAs/source,InAs/drain,poly/gate}
37    +       term1=ndiv_fbm.lapflux(psi|psi)@{InAs,ambient,nitride,poly}
38    +       term2=nodal.nscd(electrons,holes,netdope|psi)@{InAs,poly}
39    +       term3=ndiv_fbm.ncflux(psi,electrons,tl,nmob0,nin,edge|electrons)@{InAs,poly}
40    +       term4=ndiv_fbm.ncflux(psi,holes,tl,pmob0,nip,edge|holes)@{InAs,poly}
41    +       term5=constraint.continuity(psi|psi)@{InAs/nitride,poly/nitride}
42    +       term6=transient.ddt(electrons,holes|electrons,holes)@{InAs,poly}
43    +       term7=nodal.srh_fd(electrons,holes,ni,nin,nip,taun,taup|electrons)@{InAs,poly}
44    +       term8=nodal.srh_fd(electrons,holes,ni,nin,nip,taun,taup|holes)@{InAs,poly}     
45    +       term9=nodal.auger_fd(electrons,holes,nin,nip|electrons)@{InAs}
46    +       term10=nodal.auger_fd(electrons,holes,nin,nip|holes)@{InAs}
47    +       func0=eg(tl,netdope|eg)@{InAs,poly}
48    +       func1=dos(electrons,holes,tl|nc,nv)@{InAs,poly}
49    +       func2=ni(tl,eg,nc,nv|ni)@{InAs,poly}
50    +       func3=mob0(electrons,holes,tl,netdope|nmob0,pmob0)@{InAs,poly}   
51    +       func4=nic(ni,nc,electrons|nin)@{InAs,poly,nitride}
52    +       func5=nic(ni,nv,holes|nip)@{InAs,poly}
53    +       func6=tau(electrons,holes,netdope|taun,taup)@{InAs,poly}
54    +       func7=band(psi,electrons,holes,tl,ni,eg,nc,nv,chi|phin,phip,ec,ev,evac)@{InAs,ambient,poly,nitride}
55    +       init0="tl:${solve/bias/Tcelsius}+${library/physics/celsius}"
56    +       init1="edge:1"
57    +       init2="chi:0"
58   
59#first define the ambient layer
60grid dim=2  mat=ambient
61+yloc=-2,-1,-0.5,0.0,0.5,1,2
62+ydel=0.2,0.1,0.05,0.001,0.05,0.1,0.2
63+xloc=0.065,0.1,0.30
64+xdel=0.0065,0.01,0.03
65
66
67deposit mat=poly start=-1.0 end=1.0 thick=0.0,0.01
68 +xdel=0.0001
69
70deposit mat=nitride start=-1.0 end=1.0 thick=0.0,0.00499
71 +xdel=0.0005
72
73deposit mat=nitride start=-2.0 end=-1.0 thick=0,0.01499
74 +xdel=0.0001 zipper
75
76deposit mat=nitride start=1.0 end=2.0 thick=0,0.01499
77 +xdel=0.0001 zipper
78
79
80#Adding InAs in radial direction
81deposit mat=InAs start=-2 end=2 thick=0,0.001,0.01,0.050
82 +xdel=0.0001,0.001,0.001,0.005
83
84#define electrodes#
85boundary ymin=-2 ymax=-2 xmin=0.00001 xmax=0.05001 name=source
86boundary ymin=2 ymax=2 xmin=0.00001 xmax=0.05001 name=drain
87boundary ymin=-1 ymax=1 xmin=0.055 xmax=0.065 name=gate
88
89#define netdope#
90field set=netdope val=1e18 mat=InAs
91#field set=netdope val=1.0  mat=ambient
92field set=netdope val=1e21 mat=poly
93
94bias initial system=InAs_poisson
95bias system=InAs_dd elec=gate voltage=2.0
96#bias system=InAs_dd elec=drain voltage=0.0
97bias system=InAs_dd elec=drain vstep=0.1 nstep=20
98
99
100
101
102#save I-V data to file in column format"#
103graph iv outfile=FET.iv
104   
105#---------2D plot is slow in web-based Prophet. use 1D slice plot instead -------------------
106#graph quantity=psi ypos=0.0   
107#graph contour color quantity=psi