go internal
assign name=VB n.value=(0)


go atlas 

mesh width=30 outf=device.str master.out

x.m l=0 s=1
x.m l=40 s=1
y.m l=0 s=0.002
y.m l=-0.1 s=0.02
y.m l=0 s=0.002
y.m l=0.02 s=0.002
y.m l=0.12 s=0.02

region num=1 material=IGZO y.min=0    y.max=0.02
region num=2 material=sio2 y.min=0.02 y.max=0.12
region num=3 material=sio2 y.min=-0.1 y.max=0

elec num=1 name=gate   x.min=5  x.max=35 y.max=0.12 y.min=0.12
elec num=2 name=source x.min=0  x.max=5  y.max=0.02 y.min=0.02
elec num=3 name=drain  x.min=35 x.max=40 y.max=0.02 y.min=0.02
#elec num=4 name=substrate   x.min=5  x.max=35 y.max=-0.1 y.min=-0.1


contact num=1 n.poly
#
# We also define a workfunction for the source and drain that
# is very close to the conduction edge.  In the reference the
# authors observed that without a workfunction the results for
# ohmic boundaries were not significantly different than the
# Schottky model.
#
contact num=2 workf=4.33
contact num=3 workf=4.33
#contact num=4 n.poly


models fermi
#
# Key to the characterization of amorphous materials is the
# definition of the states within the band gap.
#
# a-IGZO channel bulk DOS
defect material=IGZO cont numa=128 numd=64\
      nga=${nga} ega=${ega} wga=${wga}\
      ngd=${ngd} egd=${egd} wgd=${wgd}\
      wta=${wta} nta=${nta} \
	wtd=${wtd} ntd=${ntd} \
	sigtae=1.e-14 sigtah=1.e-14 sigtde=1.e-14 sigtdh=1.e-14 \  
       siggae=1e-14 siggah=1e-14 siggde=1e-14 siggdh=1e-14  tfile=IGZOdefect-nga6e16.dat

# From here we simply extract the Id-Vg characteristic

solve init
solve prev
solve vdrain=0.1
save outf=device_VB=0.str
log outf=IdVg-nga6e16-1=0.log
solve vgate=0 vstep=-0.1 vfinal=-15 name=gate
log off

load inf=device_VB=0.str master
solve prev
log outf=IdVg-nga6e16-2_VB=0.log
solve vgate=0 vstep=0.1 vfinal=15 name=gate
log off

solve name=gate aname=gate ac freq=2000
log outf=CV-CGS-nga6e16-2000-1.log
solve name=gate vgate=0 vstep=-0.05 vstop=-3 aname=gate ac freq=2000
log off

solve init
solve name=gate aname=gate ac freq=2000
log outf=CV-CGS-nga6e16-2000-2.log
solve name=gate vgate=0 vstep=0.05 vfinal=3 aname=gate ac freq=2000
log off

quit