InGaN solar cell

 go atlas 

set temp=400

set ptype=1e15

set ntype=5e19

set thp=0.3

set thn=0.3

mesh auto

#meshing============================================================

x.m loc=0 spac=.1

x.m loc=1 spac=.1

#

y.m loc=0 spac=.2

y.m loc=0.2 spac=.1

y.m loc=0.2+$thp spac=.1

y.m loc=0.2+$thp+$thn spac=.1

#region============================================================

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

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

region num=3 material=InGaN x.min=0 x.max=1 y.min=0.2 y.max=0.2+$thp acceptor=$ptype

region num=4 material=InGaN x.min=0 x.max=1 y.min=0.2+$thp y.max=0.2+$thp+$thn donor=$ntype

#electrode============================================================

electrode num=1 name=anode x.min=0 x.max=1 y.min=0 y.max=0.2 mat=ITO

electrode num=2 name=cathode x.min=0 x.max=1 y.min=0.2+$thp+$thn y.max=0.2+$thp+$thn material=InGaN 

 

#material============================================================

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

material material=ITO eg300=3.5 affinity=2.3 permittivity=9  NC300=5e17 NV300=5e17 MUN=30 MUP=5 sopra=Ito2.nk taun0=1e-9 taup0=1e-9

material mat=InGaN EG300=1.39 PERMITTIVITY=12.8 AFFINITY=5.4 MUN0=1250 MUP0=650  NC300=5e17 NV300=5e17 taup0=1e-9 taun0=1e-9 

#

#beam============================================================

beam num=1 am1.5g x.origin=0.0 y.origin=-2.0 angle=90.0 tr.mat  

contact name=anode SURF.REC VSURFN=1e3 VSURFP=1e3 

contact name=cathode SURF.REC VSURFN=1e3 VSURFP=1e3 

models srh fermi ni.fermi auger optr print 

method newton autonr trap maxtrap=10 itlimit=40

#

solve init 

solve b1=0.001

solve b1=0.01

solve b1=0.1

solve b1=1

#

log outf=cd.log

solve lambda=0.3 wstep=0.01 wfinal=1.1

log off

extract init infile="cd.log"

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

tonyplot cd.dat

quit