Organic solar cell (with Nanoscale Layers and Films)

silvaco code:-

go atlas

#

mesh auto

x.mesh Loc=0 Spacing=1

x.mesh Loc=1 Spacing=1

y.mesh Loc=0 Spacing=0.1

y.mesh Loc=0.47 Spacing=0.1

#

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

region num=2 x.min=0 x.max=1 y.min=0.20 y.max=0.27 material=ZnO

region num=3 x.min=0 x.max=1 y.min=0.27 y.max=0.37 user.material=P3HT

region num=4 x.min=0 x.max=1 y.min=0.37 y.max=0.47 material=silver

#

electrode num=1 name=cathode x.min=0 x.max=1 y.min=0 y.max=0

electrode num=2 name=anode x.min=0 x.max=1 y.min=0.47 y.max=0.47

#

Material name=P3HT Eg300=1.95 Nc300=2e21 Nv300=2e21 Permittivity=3 Affinity=3.2 Mun=1e-6 Mup=1.2e-3 \

taun0=1e-6 taup0=1e-7 USER.GROUP=semiconductor USER.DEFAULT=silicon sopra=P3HT.nk

material material=ITO sopra=Ito2.nk

material material=silver sopra=Ag.nk

#

doping uniform conc=1E15 p.type reg=3

doping uniform conc=1E21 n.type reg=2

#

Contact name=cathode Workfun=4.3

Contact name=anode Workfun=4.7

#

save outf=p3ht.str

tonyplot p3ht.str

go atlas

mesh infile = p3ht.str

#

models conmob fldmob srh optr print

method newton trap maxtrap=10

output con.band val.band band.par e.mobility e.velocity

beam num=1 am1.5 x.o=.5 y.o=-1 angle=90 wavel.start=0.305 wavel.end=3.455 wavel.num=100 \

power.file=C:/sedatools/examples/solar/solarex09.spec

#

solve init

solve prev

solve b1=1

log outf=mie.log

solve name=anode vanode=0 vstep=0.01 vfinal=0.9

log off

tonyplot mie.log

#

extract init infile="mie.log"

extract name="IV" curve(v."Anode",i."Cathode") outfile="curveIV.dat"

#

extract name="Jsc" y.val from curve(v."Anode", i."Cathode") where x.val=0

extract name="Voc" x.val from curve(v."Anode", i."Cathode") where y.val=0

extract name="Pm" max(curve(v."Anode",(v."Anode"* i."Cathode")))

extract name="Vm" x.val from curve(v."Anode",(v."Anode"* i."Cathode")) \

where y.val=$"Pm"

extract name="Im" $"Pm"/$"Vm"

extract name="FF" $"Pm"/($"Jsc"*$"Voc")

tonyplot curveIV.dat

Quit

result:-

if you need p3ht  then i will write below and save it as P3HT.nk file

12

0.25024 1.55318 0.25475

0.33525 1.48282 0.33525

0.41156 1.23235 0.43828

0.50558 1.72542 0.50558

0.56804 2.24214 0.56354

0.59058 2.17122 0.59509

0.60829 2.1791 0.63083

0.75609 1.8425 0.87265

1.00241 1.75638 1.09192

1.25324 1.7333 1.42325

1.49505 1.71754 1.4268

(updated) perovskite solar cell

Silvaco code:-

go atlas

mesh auto

x.m l=0.0 s=1

x.m l=1.0 s=1

y.m l=0.175 s=0.01

y.m l=0.375 s=0.01

y.m l=0.4 s=.01

y.m l=0.9 s=.1

y.m l=1.4 s=.1

y.m l=1.6 s=.1

region num=1 mat=ITO x.min=0 x.max=1.0 y.min=.175 y.max=0.375

region num=2 user.material=graphene x.min=0 x.max=1.0 y.min=0.375 y.max=0.4

region num=3 user.material=perovskite x.min=0 x.max=1.0 y.min=.4 y.max=.9

region num=4 user.material=spiro x.min=0 x.max=1.0 y.min=.9 y.max=1.6

elec num=1 name=anode x.min=0.0 x.max=1.0 y.min=0.175 y.max=0.175

elec num=2 name=cathode x.min=0 x.max=1.0 y.min=1.6 y.max=1.7 mat=gold

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

material material=ITO sopra=Ito2.nk

material material=graphene eg300=0.25 affinity=4.1 permittivity=25 NC300=2e18 NV300=1.8e18 MUN=1.8e4 MUP=1.4e4 user.group=conductor user.default=sio2

material material=pervoskite eg300=1.55 affinity=3.93 permittivity=100 NC300=3.97e19 NV300=3.97e17 MUN=2 MUP=2 user.group=conductor user.default=sio2

material material=spiro eg300=3.0 affinity=2.45 permittivity=3 NC300=2e18 NV300=2e19 MUN=0.02 MUP=0.02 user.group=conductor user.default=sio2

doping uniform n.type conc=1e15 x.min=0 x.max=1.0 y.min=.375 y.max=.4

doping uniform p.type conc=1e19 x.min=0 x.max=1.0 y.min=.4 y.max=.9

contact name=anode

contact name=cathode workfun=5.2

save outf= struc.str

tonyplot struc.str

#############

go atlas

mesh infile = struc.str

models conmob fldmob srh optr chuang print

method newton autonr trap maxtrap=10

output con.band val.band band.par e.mobility e.velocity

beam num=1 am1.5 x.o=.5 y.o=-.5 angle=90 wavel.start=0.305 wavel.end=3.455 wavel.num=100 \

power.file=C:/sedatools/examples/solar/solarex09.spec

solve init

solve prev

solve b1=0

log outf=mie.log

solve b1=1 beam=1 lambda=0.3 wstep=0.025 wfinal=1.0

extract init infile="mie.log"

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

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

tonyplot -overlay EQE1_0.dat EQE1_1.dat

solve init

solve prev

solve b1=0

log outfile=DARK.log

solve name=anode vanode=-.1 vstep=0.01 vfinal=.2

log off

tonyplot DARK.log

Solve init

solve prev

solve b1=1

log outf=illumination.log

solve name=anode vanode=-.1 vstep=0.01 vfinal=.2

log off

tonyplot -overlay DARK.log illumination.log

#parameters

extract init infile="illumination.log"

extract name="Jsc" max(curve(v."anode", i."cathode"))

extract name="JscmAcm2" $Jsc*1e08*1e03

extract name="Voc" x.val from curve(v."anode",i."cathode") where y.val=0.0

extract name="Pm" max(curve(v."anode", (v."anode" *i."cathode")))

extract name="Vm" x.val from curve(v."anode",(v."anode"*i."cathode") ) \

where y.val=$"Pm"

extract name="Im" $"Pm"/$"Vm"

extract name="FF" ($"Pm"/($"Jsc"*$"Voc"))*100

extract name="Opt_int" max(beam."1")

extract name="Eff" (1e8*$Pm/$Opt_int)*100

extract name="Power" curve(v."anode", (v."anode" * i."anode" *(-1))) outf="Power_a.dat"

tonyplot Power_a.dat

save outf= finalstrucre.str

tonyplot finalstrucre.str

quit

result:-