go atlas simflags="-p 4"

mesh three.d

#

x.m l=-0.2 spac=0.02

x.m l=-0.02 spac=0.005

x.m l=0.02 spac=0.005

x.m l=0.2 spac=0.02

#

y.m l=-0.275 spac=0.02

y.m l=-0.135 spac=0.04

y.m l=-0.045 spac=0.005

y.m l=0.045 spac=0.005

y.m l=0.135 spac=0.04

y.m l=0.275 spac=0.02

#

#

z.m l=-0.1 spac=0.03

z.m l=-0.005 spac=0.001

z.m l=0.06 spac=0.005

z.m l=0.21 spac=0.2

z.m l=0.5 spac=0.2

#

region num=1 material= silicon x.min=-0.2 x.max=0.2 y.min=-0.275 y.max=0.275 \

z.min=0.21 z.max=0.5

region num=2 material= sio2 x.min=-0.2 x.max=0.2 y.min=-0.275 y.max=0.275 \

z.min=0.06 z.max=0.21

region num=3 material= poly x.min=-0.2 x.max=0.2 y.min=-0.045 y.max=0.045 \

z.min=-0.1 z.max=0.06

region num=4 material= hfo2 x.min=-0.025 x.max=0.025 y.min=-0.045 y.max=0.045 \

z.min=-0.005 z.max=0.06

region num=5 material= silicon x.min=-0.02 x.max=0.02 y.min=-0.135 y.max=0.135 \

z.min=0 z.max=0.06

region num=6 material= silicon x.min=-0.02 x.max=0.02 y.min=-0.275 y.max=-0.135 \

z.min=0 z.max=0.06

region num=7 material= silicon x.min=-0.02 x.max=0.02 y.min=0.135 y.max=0.275 \

z.min=0 z.max=0.06

#

elec num=1 x.min=-0.2 x.max=0.2 y.min=-0.045 y.max=0.045 z.min=-0.1 z.max=-0.1 name=gate

elec num=2 x.min=-0.02 x.max=0.02 y.min=-0.275 y.max=-0.045 z.min=0 z.max=0 name=source

elec num=3 x.min=-0.02 x.max=0.02 y.min=0.045 y.max=0.275 z.min=0 z.max=0 name=drain

#

doping uniform p.type conc=1.e17 x.min=-0.2 x.max=0.2 y.min=-0.275 y.max=0.275 z.min=0.21 z.max=0.5

doping uniform p.type conc=1.e17 x.min=-0.02 x.max=0.02 y.min=-0.135 y.max=0.135 z.min=0 z.max=0.06

doping uniform n.type conc=3.e22 x.min=-0.02 x.max=0.02 y.min=-0.275 y.max=-0.135 z.min=0 z.max=0.06

doping uniform n.type conc=3.e22 x.min=-0.02 x.max=0.02 y.min=0.135 y.max=0.275 z.min=0 z.max=0.06

save outf=fin_0.str

interface qf=3e10

contact name=gate n.poly

models conmob srh auger bgn fldmob print

solve init

method newton trap

solve vdrain=0.05

log outf=fin_0.log

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

extract name="vt" xintercept(maxslope(curve(v."gate",i."drain"))) \

- ave(v."drain")/2.0

quit