MFIS finfet

go devedit simflags="-3d"

work.area x1=-0.01 y1=-0.01 x2=0.02 y2=0.045

#

#

region reg=1 mat=Silicon color=0xffcc00 pattern=0x4 z1=0.01 z2=0.09 \

polygon="0,0 0.01,0 0.01,0.015 0,0.015"

constr.mesh region=1 default max.height=0.01 max.width=0.02

#

#

region reg=2 mat="Silicon Oxide" color=0xff pattern=0x2 z1=0.01 z2=0.09 \

polygon="0.013,0.015 0.01,0.015 0.01,0 0,0 0,0.015 -0.003,0.015 -0.003,-0.003 0.013,-0.003"

constr.mesh region=2 default max.height=0.001 max.width=0.01

#

#

region reg=3 name=gate mat=PolySilicon elec.id=1 work.func=0 color=0xffff00 pattern=0x5 z1=0.035 z2=0.065 \

polygon="0.0165,-0.0065 -0.0065,-0.0065 -0.0065,0.015 -0.01,0.015 -0.01,-0.01 0.02,-0.01 0.02,0.015 0.0165,0.015"

#polygon="0.013,-0.003 -0.003,-0.003 -0.003,0.015 -0.01,0.015 -0.01,-0.01 0.02,-0.01 0.02,0.015 0.013,0.015"

constr.mesh region=3 default max.height=0.01 max.width=0.02

#

#

region reg=4 mat="Silicon Oxide" color=0xff pattern=0x2 z1=0 z2=0.1 \

polygon="-0.01,0.015 -0.003,0.015 0,0.015 0.01,0.015 0.013,0.015 0.02,0.015 0.02,0.045 -0.01,0.045"

constr.mesh region=4 default max.height=0.01 max.width=0.02

#

#

region reg=5 name=drain mat=Aluminum elec.id=2 work.func=0 color=0xffc8c8 pattern=0x7 z1=0 z2=0.01 \

polygon="0,0 0.01,0 0.01,0.015 0,0.015"

constr.mesh region=5 default max.height=0.01 max.width=0.02

#

#

region reg=6 name=source mat=Aluminum elec.id=3 work.func=0 color=0xffc8c8 pattern=0x7 z1=0.09 z2=0.1 \

polygon="0,0 0.01,0 0.01,0.015 0,0.015"

constr.mesh region=6 default max.height=0.01 max.width=0.02

#

#

region reg=7 name=ferric material=HFO2 work.func=0 color=0xffc8c8 pattern=0x5 z1=0.035 z2=0.065 \

polygon="0.013,-0.003 -0.003,-0.003 -0.003,0.015 -0.0065,0.015 -0.0065,-0.0065 0.0165,-0.0065 0.0165,0.015 0.013,0.015"

constr.mesh region=3 default max.height=0.01 max.width=0.001

#

#

# Set Meshing Parameters

#

base.mesh height=0.01 width=0.005

#

bound.cond !apply max.slope=28 max.ratio=300 rnd.unit=0.0001 line.straightening=1 align.points when=automatic

#

imp.refine min.spacing=0.02 z=0

#

constr.mesh max.angle=90 max.ratio=300 max.height=10000 \

max.width=10000 min.height=0.0001 min.width=0.0001

#

constr.mesh type=Semiconductor default

#

constr.mesh type=Insulator default

#

constr.mesh type=Metal default

#

constr.mesh type=Other default

#

Mesh Mode=MeshBuild

z.plane z=0 spacing=0.1

z.plane z=0.005 spacing=0.1

z.plane z=0.01 spacing=0.1

z.plane z=0.012 spacing=0.1

z.plane z=0.0167 spacing=0.1

z.plane z=0.0233 spacing=0.1

z.plane z=0.0267 spacing=0.1

z.plane z=0.030 spacing=0.1

z.plane z=0.0315 spacing=0.1

z.plane z=0.033 spacing=0.1

z.plane z=0.035 spacing=0.1

z.plane z=0.041 spacing=0.1

z.plane z=0.047 spacing=0.1

z.plane z=0.053 spacing=0.1

z.plane z=0.059 spacing=0.1

z.plane z=0.065 spacing=0.1

z.plane z=0.067 spacing=0.1

z.plane z=0.0685 spacing=0.1

z.plane z=0.070 spacing=0.1

z.plane z=0.0733 spacing=0.1

z.plane z=0.0767 spacing=0.1

z.plane z=0.0833 spacing=0.1

z.plane z=0.088 spacing=0.1

z.plane z=0.090 spacing=0.1

z.plane z=0.095 spacing=0.1

z.plane z=0.100 spacing=0.1

z.plane max.spacing=1000000 max.ratio=1.5

structure outf=ncfinfet_0.str

tonyplot3d ncfinfet_0.str

go atlas 

electrode name=bulk bottom

doping num=1 n.type uniform conc=1e19 reg=1

material material=Silicon eg300=1.1 affinity=4.05 permitti=11.9 

material material=SiO2  mun0=20.0 mup0=1.0e-5

#

MATERIAL MATERIAL=HFO2 FERRO.EPS=30.25 FERRO.PR=9.0e-6 FERRO.PS=9.5e-6 FERRO.EC=1.1e6

#

contact   name=gate p.poly workfunction=5.1

contact   name=drain workfunction=4.1

contact   name=source workfunction=4.1

MODELS REGION=4 FERRO

method autonr gummel maxtrap=5

solve init

solve vdrain=0.01

solve vdrain=0.05

solve vdrain=0.1

log outf=ncfinfet.log

solve  vgate=0 name=gate vstep=0.05 vfinal=1.1

log off

output  band.temp con.band val.band band.par

save outf=ncfinfet_1.str master

tonyplot ncfinfet.log 

tonyplot3d ncfinfet_1.str 

extract name="subvt" 1.0/slope(maxslope(curve(v."gate",log10(abs(i."drain")))))

extract name="Ioff_n" y.val from curve(v."gate",i."drain") where x.val=0

extract name="Ion_n" y.val from curve(v."gate",i."drain") where x.val=1.0

quit

                       

i will upload more code in future

InGaN solar cell

 go atlas 

set temp=400

set ptype=1e15

set ntype=5e19

set thp=0.3

set thn=0.3

mesh auto

#meshing============================================================

x.m loc=0 spac=.1

x.m loc=1 spac=.1

#

y.m loc=0 spac=.2

y.m loc=0.2 spac=.1

y.m loc=0.2+$thp spac=.1

y.m loc=0.2+$thp+$thn spac=.1

#region============================================================

region num=1 mat=air x.min=0 y.min=0

region num=2 mat=ITO x.min=0 x.max=1 y.min=0 y.max=0.2

region num=3 material=InGaN x.min=0 x.max=1 y.min=0.2 y.max=0.2+$thp acceptor=$ptype

region num=4 material=InGaN x.min=0 x.max=1 y.min=0.2+$thp y.max=0.2+$thp+$thn donor=$ntype

#electrode============================================================

electrode num=1 name=anode x.min=0 x.max=1 y.min=0 y.max=0.2 mat=ITO

electrode num=2 name=cathode x.min=0 x.max=1 y.min=0.2+$thp+$thn y.max=0.2+$thp+$thn material=InGaN 

 

#material============================================================

material material=air real.index=1 imag.index=0

material material=ITO eg300=3.5 affinity=2.3 permittivity=9  NC300=5e17 NV300=5e17 MUN=30 MUP=5 sopra=Ito2.nk taun0=1e-9 taup0=1e-9

material mat=InGaN EG300=1.39 PERMITTIVITY=12.8 AFFINITY=5.4 MUN0=1250 MUP0=650  NC300=5e17 NV300=5e17 taup0=1e-9 taun0=1e-9 

#

#beam============================================================

beam num=1 am1.5g x.origin=0.0 y.origin=-2.0 angle=90.0 tr.mat  

contact name=anode SURF.REC VSURFN=1e3 VSURFP=1e3 

contact name=cathode SURF.REC VSURFN=1e3 VSURFP=1e3 

models srh fermi ni.fermi auger optr print 

method newton autonr trap maxtrap=10 itlimit=40

#

solve init 

solve b1=0.001

solve b1=0.01

solve b1=0.1

solve b1=1

#

log outf=cd.log

solve lambda=0.3 wstep=0.01 wfinal=1.1

log off

extract init infile="cd.log"

extract name="EQE2" curve(elect."optical wavelength",-(i."anode")/elect."source photo current"*) outf="EQE01_1.dat" 

tonyplot cd.dat

quit

NC-FINFET

 code:-

go devedit simflags="-3d"

work.area x1=-0.0135 y1=-0.0135 x2=0.0235 y2=0.045

#

#

region reg=1 mat=Silicon color=0xffcc00 pattern=0x4 z1=0.01 z2=0.09 \

polygon="0,0 0.01,0 0.01,0.015 0,0.015"

constr.mesh region=1 default max.height=0.01 max.width=0.02

#

#

region reg=2 name=gate mat=Polysilicon elec.id=1 work.func=0 color=0xffc8c8 pattern=0x7 z1=0.035 z2=0.065 \

polygon="-0.0135,0.015 -0.01,0.015 -0.01,-0.01 0.02,-0.01 0.02,0.015 0.0235,0.015 0.0235,-0.0135 -0.0135,-0.0135 -0.0135,0.015"

constr.mesh region=2 default max.height=0.01 max.width=0.02

#

#

region reg=3 mat="Silicon Oxide" color=0xff pattern=0x2 z1=0.01 z2=0.09 \

polygon="0.013,0.015 0.01,0.015 0.01,0 0,0 0,0.015 -0.003,0.015 -0.003,-0.003 0.013,-0.003"

constr.mesh region=3 default max.height=0.001 max.width=0.001

#

#

region reg=4 name=ferric  mat=HfO2  work.func=0 color=0xffff00 pattern=0x5 z1=0.035 z2=0.065 \

polygon="0.0165,-0.0065 -0.0065,-0.0065 -0.0065,0.015 -0.01,0.015 -0.01,-0.01 0.02,-0.01 0.02,0.015 0.0165,0.015"

constr.mesh region=4 default max.height=0.01 max.width=0.02

#

#

region reg=5 mat="Silicon Oxide" color=0xff pattern=0x2 z1=0 z2=0.1 \

polygon="-0.0135,0.015 -0.003,0.015 0,0.015 0.01,0.015 0.013,0.015 0.0235,0.015 0.0235,0.045 -0.0135,0.045"

constr.mesh region=5 default max.height=0.01 max.width=0.02

#

#

region reg=6 name=drain mat=Aluminum elec.id=2 work.func=0 color=0xffc8c8 pattern=0x7 z1=0 z2=0.01 \

polygon="0,0 0.01,0 0.01,0.015 0,0.015"

constr.mesh region=6 default max.height=0.01 max.width=0.02

#

#

region reg=7 name=source mat=Aluminum elec.id=3 work.func=0 color=0xffc8c8 pattern=0x7 z1=0.09 z2=0.1 \

polygon="0,0 0.01,0 0.01,0.015 0,0.015"

constr.mesh region=7 default max.height=0.01 max.width=0.02

#

#

region reg=8  material=Polysilicon work.func=0 color=0xffc8c8 pattern=0x5 z1=0.035 z2=0.065 \

polygon="0.013,-0.003 -0.003,-0.003 -0.003,0.015 -0.0065,0.015 -0.0065,-0.0065 0.0165,-0.0065 0.0165,0.015 0.013,0.015"

constr.mesh region=8 default max.height=0.01 max.width=0.02

#

#

# Set Meshing Parameters

#

base.mesh height=0.01 width=0.005

#

bound.cond !apply max.slope=28 max.ratio=300 rnd.unit=0.0001 line.straightening=1 align.points when=automatic

#

imp.refine min.spacing=0.02 z=0

#

constr.mesh max.angle=90 max.ratio=300 max.height=10000 \

max.width=10000 min.height=0.0001 min.width=0.0001

#

constr.mesh type=Semiconductor default

#

constr.mesh type=Insulator default

#

constr.mesh type=Metal default

#

constr.mesh type=Other default

#

constr.mesh region=1 default max.height=0.01 max.width=0.02

#

constr.mesh region=2 default max.height=0.001 max.width=0.001

#

constr.mesh region=3 default max.height=0.01 max.width=0.02

#

constr.mesh region=4 default max.height=0.01 max.width=0.02

#

constr.mesh region=5 default max.height=0.01 max.width=0.02

#

constr.mesh region=6 default max.height=0.01 max.width=0.02

Mesh Mode=MeshBuild

z.plane z=0 spacing=0.1

z.plane z=0.005 spacing=0.1

z.plane z=0.01 spacing=0.1

z.plane z=0.012 spacing=0.1

z.plane z=0.0167 spacing=0.1

z.plane z=0.0233 spacing=0.1

z.plane z=0.0267 spacing=0.1

z.plane z=0.030 spacing=0.1

z.plane z=0.0315 spacing=0.1

z.plane z=0.033 spacing=0.1

z.plane z=0.035 spacing=0.1

z.plane z=0.041 spacing=0.1

z.plane z=0.047 spacing=0.1

z.plane z=0.053 spacing=0.1

z.plane z=0.059 spacing=0.1

z.plane z=0.065 spacing=0.1

z.plane z=0.067 spacing=0.1

z.plane z=0.0685 spacing=0.1

z.plane z=0.070 spacing=0.1

z.plane z=0.0733 spacing=0.1

z.plane z=0.0767 spacing=0.1

z.plane z=0.0833 spacing=0.1

z.plane z=0.088 spacing=0.1

z.plane z=0.090 spacing=0.1

z.plane z=0.095 spacing=0.1

z.plane z=0.100 spacing=0.1

z.plane max.spacing=1000000 max.ratio=1.5

structure outf=ncfinfet_0.str

go atlas 

electrode name=bulk bottom

doping num=1 n.type uniform conc=1e18 reg=6

doping num=2 n.type uniform conc=1e18 reg=7

material material=Silicon eg300=1.1245 affinity=4.05 permitti=11.9 

material material=HFO2 ferro.ps=0.5e-6 ferro.pr=0.4e-6 ferro.ec=100000.0 ferro.eps=10.0

material material=Oxide eg300=8.05 affinity=1.00 permittivity=3.9  nc300=2.8e19 nv300=1.04e19

#

contact name=gate workf=5.1

model fermi bqp.n srh ferro hcte.el  evsatmod=0 fldmob print

method newton dvlimit=1.0 maxtrap=6 

solve init

solve vdrain=0.01

solve vdrain=0.05

solve vdrain=0.1

log outf=ncfinfet.log

solve vgate=0.0 name=gate vstep=0.1 vfinal=1.0

log off

output p.quantum band.temp con.band val.band band.par

save outf=ncfinfet_1.str master

tonyplot ncfinfet.log 

tonyplot3d ncfinfet_1.str -set ncfinfet_0.set

tonyplot3d ncfinfet_1.str -set ncfinfet_1.set

extract name="subvt" 1.0/slope(maxslope(curve(v."gate",log10(abs(i."drain")))))

extract name="Ioff_n" y.val from curve(v."gate",i."drain") where x.val=0

extract name="Ion_n" y.val from curve(v."gate",i."drain") where x.val=1

quit

result:-