OFET(Also you can know how to use pentacene and poole frenkel mobility here)

 go atlas

mesh space.mult=1

x.m l=0 s=0.5

x.m l=10 s=0.5

x.m l=15 s=0.5

x.m l=20 s=0.5

x.m l=30 s=0.5

y.m l=-0.03 s=0.01

y.m l=0 s=0.005

y.m l=0.015 s=0.003

y.m l=0.03 s=0.005

y.m l=0.0325 s=0.001

y.m l=0.0357 s=0.001

y.m l=0.0557 s=0.0025

region num=1 material=Pentacene x.min=0 x.max=30 y.min=-0.03 y.max=0.03 acceptor=7e17

region num=2 material=Al2O3 x.min=0 x.max=30 y.min=0.03 y.max=0.0357

electrode num=1 name=source x.max=10 y.min=0.0 y.max=0.03 material=Gold

electrode num=2 name=gate x.min=0 x.max=30 y.min=0.0357 y.max=0.0557 material=Aluminium

electrode num=3 name=drain x.min=20 y.min=0.00 y.max=0.03 material=Gold

material material=Pentacene eg300=2.8 nc300=1e21 nv300=1e21 permittivity=4.0

material material=Al2O3 permittivity=4.5

#

mobility region=1 deltaep.pfmob=0.0179 betap.pfmob=7.75e-5 mun=5e-5 mup=0.85

#

models conmob srh fermi print

models region=1 pfmob.p

#

method newton autonr trap

#

contact name=gate workfun=4.6

contact name=source workfun=5.1

contact name=drain workfun=5.1

#

solve init

solve prev

solve vdrain=2.5

log outf=abc0.log

solve vgate=-20 name=gate vstep=1 vfinal=20

tonyplot abc0.log

output val.band con.band band.param

save outf=abc.str

tonyplot abc.str

Tunnel field effect transistor -VTFET

go atlas

mesh auto

#

x.mesh loc=0 spac=0.1

x.mesh loc=.002 spac=.001

x.mesh loc=.004 spac=.001

x.mesh loc=.005 spac=.005

x.mesh loc=.007 spac=.0001

x.mesh loc=.017 spac=.0001

x.mesh loc=.019 spac=.005

x.mesh loc=.020 spac=.001

x.mesh loc=.022 spac=.001

x.mesh loc=.024 spac=0.1

#

y.mesh loc=0 spac=.005

y.mesh loc=.010 spac=.001

y.mesh loc=.018 spac=.001

y.mesh loc=.022 spac=.001

y.mesh loc=.040 spac=.0001

y.mesh loc=.042 spac=.0001

y.mesh loc=.060 spac=.0005

y.mesh loc=.061 spac=.0005

y.mesh loc=.079 spac=.005

#

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

region num=2 material=silicon x.min=.007 x.max=.017 y.min=.001 y.max=.021 

region num=3 material=silicon x.min=.007 x.max=.017 y.min=.021 y.max=.041 

region num=4 material=silicon x.min=.007 x.max=.017 y.min=.041 y.max=.061 

region num=5 material=silicon x.min=0 x.max=.024 y.min=.061 y.max=.081

region num=6 material=hfo2 x.min=0.005 x.max=.007 y.min=.021 y.max=.041

region num=7 material=hfo2 x.min=0.017 x.max=.019 y.min=.021 y.max=.041

#

qtx.mesh loc=0.005 spac=0.01

qtx.mesh loc=0.020 spac=0.001

qtx.mesh loc=0.019 spac=0.001

qtx.mesh loc=0.022 spac=0.01

#

qty.mesh loc=0.002 spac=0.005

qty.mesh loc=0.003 spac=0.0005

qty.mesh loc=0.045 spac=0.0005

qty.mesh loc=0.055 spac=0.0005

#

electrode name=gate x.min=.003 x.max=.005 y.min=.021 y.max=.041

electrode name=gate2 x.min=.019 x.max=.021 y.min=.021 y.max=.041 

electrode name=source x.min=.007 x.max=.017 y.min=0 y.max=0.001 

electrode name=drain x.min=0 x.max=.007 y.min=0.059 y.max=0.061 

electrode name=drain2 x.min=.017 x.max=.024 y.min=0.059 y.max=0.061 

#

material material=silicon me.tunnel=.24 mh.tunnel=.20 

#

doping uniform p.type conc=2e19 reg=2

doping uniform p.type conc=5e16 reg=3

doping uniform n.type conc=5e18 reg=4

#

models bbt.nonlocal qtunn.dir=0 bgn srh ni.fermi conmob cvt print temp=300

#

contact name=gate workfun=4.2

contact name=gate2 workfun=4.2 common=gate

contact name=source 

contact name=drain 

contact name=drain2  common=drain

save outf= abc.str

tonyplot abc.str

#

#

method newton autonr trap maxtrap=10

#solve

solve init

solve prev

solve vdrain=0.1

solve vdrain=0.5

solve vgate=-.3

solve vgate=-.2

log outf=tfet.log 

solve name=gate vgate=-0.2 vstep=0.05 vfinal=2.0

tonyplot tfet.log 

output val.band con.band qfn qfp e.field j.electron j.hole j.conduction j.total 

save outf= abc.str

tonyplot abc.str

quit

result:-