Textured solar cell (Spectral response)

 code:---

go athena

#

line x loc=0.00 spac=0.1

line x loc=0.25 spac=0.1

line x loc=0.50 spac=0.1

line x loc=0.75 spac=0.1

line x loc=1.0 spac=0.1

line x loc=1.25 spac=0.1

line x loc=1.5  spac=0.1

line x loc=1.75  spac=0.1

line x loc=2.0  spac=0.1

line x loc=2.125  spac=0.1

#

line y loc=0.00 spac=0.1

line y loc=4.00 spac=0.1

line y loc=8.00 spac=0.1

line y loc=12.00 spac=0.5

line y loc=20.00 spac=2.0 

line y loc=24.00 spac=2.0

# doped with phosphorous

init silicon c.phosphor=1.0e15 two.d orientation=100

#ecthing process

etch silicon start x=0 y=0

etch cont x=0.25 y=0

etch cont x=0.5 y=4

etch cont x=0.75 y=4

etch cont x=1 y=0

etch cont x=1.25 y=0

etch cont x=1.5 y=4

etch cont x=1.75 y=4

etch done x=2.0 y=0

struct mirror right

save outf=abc.str

tonyplot abc.str

go atlas

mesh infile=abc.str

electrode num=1 name=anode x.min=0.0 x.max=.25 top mat=silver

electrode num=2 name=anode x.min=1 x.max=1.25 top mat=silver

electrode num=3 name=anode x.min=2 x.max=2.25 top mat=silver

electrode num=4 name=anode x.min=3 x.max=3.25 top mat=silver

electrode num=5 name=anode x.min=4 x.max=4.25 top mat=silver

electrode num=6 name=cathode y.min=24 y.max=24  mat=silver

material material=silicon mun0=1000 mup0=500 eg300=2.0 taun0=1e-6 taup0=1e-6 sopra=Si100_2.nk

beam num=1 am1.5g x.origin=0 y.origin=-1.0 angle=90 wavel.start=0.4 wavel.end=1 wavel.num=100 

mat taun0=1e-6 taup0=1e-6

output val.band con.band opt.int

models conmob fldmob srh auger print

method  newton autonr trap

contact name=anode workfun=5.1

contact name=cathode 

solve init

solve b1=0.001

solve b1=0.01

solve b1=0.1

solve b1=1

log outf=iv.log

solve b1=1 lambda=0.3 wstep=0.01 wfinal=1.0 

log off

tonyplot  iv.log

extract init infile="iv.log"

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

tonyplot EQE01_1.dat

save outf=abc.str

tonyplot abc.str

result:---