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

#go atlas
go atlas simflags="-P 4"
#go atlas simflags="-V 5.30.0.R -P 4"

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

x.m l=0 s=1
x.m l=60 s=1
#SiO2
y.m l=0 s=0.02
y.m l=-0.15 s=0.02
#IGZO
y.m l=0 s=0.002
y.m l=0.03 s=0.002

set Eg={{Eg}}

## DOS based on the reference paper
# acceptor-like tail state
set nta={{nta}}
set wta={{wta}}
# acceptor-like Gaussian state
set nga={{nga}}
set wga={{wga}}
# donor-like tail state
set ntd={{ntd}}
set wtd={{wtd}}
# donor-like Gaussian state
set ngd={{ngd}}
set egd={{egd}}
set wgd={{wgd}}
#drift mobility
set mun={{mun}}

# IGZO channel doping
set cdp={{cdp}}

set T={{T}}

set dosfile="{{defectaccfile}}"
set IdVglog="{{atlaslogfile}}"

region num=1 user.material=a-IGZO y.min=0    y.max=0.03
region num=2 material=sio2 y.min=-0.15 y.max=0

elec num=1 name=gate   x.min=0  x.max=60 y.max=-0.15 y.min=-0.15
elec num=2 name=source x.min=0  x.max=5  y.max=0.03 y.min=0.03
elec num=3 name=drain  x.min=55 x.max=60 y.max=0.03 y.min=0.03

material material=a-IGZO user.default=Silicon  user.group=semiconductor \
                         affinity=4.3 eg300=$Eg nc300=5e18 nv300=5e19 \
                         permittivity=13.0 \
                         mun=$mun mup=0.1
doping uniform n.type concentration=$cdp region=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

# models
models srh fermi temp=$T print
#
# Key to the characterization of amorphous materials is the
# definition of the states within the band gap.

#
# a-IGZO channel bulk DOS
defects reg=1 cont numa=128 numd=128 \
	nga=$nga ega=0.0 wga=$wga \
      	ngd=$ngd egd=$egd wgd=$wgd \
      	wta=$wta nta=$nta \
      	wtd=$wtd ntd=$ntd \
      	sigtae=1e-17 sigtah=1e-15 sigtde=1e-15 sigtdh=1e-17 \
      	siggae=2e-16 siggah=2e-15 siggde=2e-15 siggdh=2e-16 \
      	tfile="$dosfile"

## Interface properties
# energy of IF acceptor-like state
#set aelev=0.5
# density of IF acceptor-like state
#set adens=1e18
# energy of IF donor-like state
#set delev=3.15
# density of IF donor-like state
#set ddens=1e20

# interface properties
#inttrap e.level=$aelev acceptor density=$adens s.i degen=1 \ 
#	y.min=0.029 y.max=0.03 \
#	sign=1e-16 sigp=1e-14
#inttrap e.level=$delev donor density=$ddens s.i degen=1 \ 
#	y.min=0.029 y.max=0.03 \
#	sign=1e-16 sigp=1e-14

#
# Id-Vg simulation
{{IdVgSimulation}}