我建立了InGaAs/InP的红外光电探测器,其光谱响应应该是在1.0um至1.65um。而我仿真出来的结果为0.6um至1um

求教为什么会出现这种结果?
附上代码:
go atlas
#mesh
mesh space.mult=1.0
x.mesh loc=0.0 spac=0.25
x.mesh loc=5.0 spac=0.25
y.mesh loc=0.0 spac=0.5
y.mesh loc=1.0 spac=0.25
y.mesh loc=3.5 spac=0.25
y.mesh loc=4.0 spac=0.5
y.mesh loc=6.0 spac=1.0
#region
region num=1 material=InP y.max=1.0
region num=2 material=InGaAs y.min=1.0 y.max=3.5 compx.top=0.53 compx.bottom=0.53
region num=3 material=InP y.min=3.5 y.max=4.0
region num=4 material=InP Y.min=4.0
#electrode
electrode num=1 name=anode x.min=0.0 x.max=0.75 y.max=0.0
electrode num=2 name=anode x.min=4.25 x.max=5.0 y.max=0.0
electrode num=3 name=cathode y.min=6.0
#doping
doping uniform region=1 p.type conc=5e18
doping uniform region=2 n.type conc=1e16
doping uniform region=3 n.type conc=5e18
doping uniform region=4 n.type conc=5e11
save outf=ingaas.str
go atlas
init infile=ingaas.str
#material
material material=InP align=0.36 mun=5000 mup=250 augn=9e-31 augp=9e-31 taun0=6e-12 taup0=6e-12 copt=1.2e-10 permittivity=12.5 eg300=1.34
material material=InGaAs align=0.36 nc300=2.1e17 nv300=7.7e18 mun=10000 mup=500 augn=2e-27 augp=2e-27 taun0=7e-6 taup0=7e-6 copt=9.6e-11 permittivity=13.9 eg300=0.75
#model
model srh fldmob conmob auger bbt.std print fermidirac
#light
beam num=1 x.origin=2.5 y.origin=-1.0 angle=90 wavelength=1.3 rays=101
#method
method newton trap carr=2
output band.param con.band val.band opt.intens
#
solve init
solve vcathode=0.05 vstep=0.05 vfinal=0.5 name=cathode
#
log outf=response.log master
solve prev b1=1 lambda=0.7
solve prev b1=1 lambda=0.725
....
solve prev b1=1 lambda=1.75
solve prev b1=1 lambda=1.775
solve prev b1=1 lambda=1.8
log off
save outf=ingaas.str
tonyplot ingaas.str
tonyplot response.log
quit