silvaco mixedmode 加入光照 问题

本人在silvaco mixedmode 里面仿真一个加入光照情况下的光电二极管，即时窗口里面跑到显示performing 2D ray trace 就跑不动了。不知道和原因，求大神指教，附上带代码。
###PD_tes
go atlas
#
# DC point
#
.begin
v1 0 2 0
adiode 2=cathode 1=anode width =10 infile =diode01.str
r1 1 0 10
o1 1 1
.nodeset v(1)=0 v(2)=0
.numeric lte=1e-3toldc=1e-4 toltr=1e-3 vchange =0.001 imaxdc=50000
.options fulln print debug
.log outfile =ph_diode
.save outfile =ph_diode_dc_point
.end
material material =alu imag.index=1000
material material =silicon taun0=1e-6 taup0=1e-6
#beam num=1 x.origin=1.6 y.origin=-2.0 angle=90.0 power.file=solarex01.spec
beamnum=1 x.origin =1.6 y.origin=-2 angle=90 wavelength =0.3
models conmob fldmob consrh
method gummel newton atrap=0.0001maxtrap=10
save outfile =ph_diode.str
go atlas
tonyplot ph_diode.str
.begin
v1 0 2 0
adiode 2=cathode 1=anode width =10 infile =diode01.str
r1 1 0 10
o1 1 1
.nodeset v(1)=0 v(2)=0
.numeric lte=1e-3toldc=1e-4 toltr=1e-3 vchange =0.001 imaxdc=50000
.options fulln print debug
.load infile=ph_diode_dc_point
.log outfile =ph_diode
.save master =ph_diode
.dc v1 0 5 0.25
.end
material material =alu imag.index=1000
material material =silicon taun0=1e-6 taup0=1e-6
#beam num=1 x.origin=1.6 y.origin=-2.0 angle=90.0 power.file=solarex01.spec
beamnum=1 x.origin =1.6 y.origin=-2 angle=90 wavelength =0.3
models conmob fldmob consrh
method gummel newton atrap=0.0001maxtrap=10
go atlas
tonyplotph_diode_dc_1.log
.begin
v1 0 2 0 pulse 0 3 0 2ns 10ns 20us 100us
adiode 2=cathode 1=anode width =10 infile =diode01.str
r1 1 0 10
o1 1 1 pulse 1 4 0 2ns 10ns 50us 100us
.nodeset v(1)=0 v(2)=0

.options fulln print debug
.load infile=ph_diode_dc_point
.log outfile =ph_diode

.tran 1ns 100us
.end
material material =alu imag.index=1000
material material =silicon taun0=1e-6 taup0=1e-6
#beam num=1 x.origin=1.6 y.origin=-2.0 angle=90.0 power.file=solarex01.spec
beamnum=1 x.origin =1.6 y.origin=-2 angle=90 wavelength =0.3
models conmob fldmob consrh
method gummel newton atrap=0.0001maxtrap=10
go atlas

tonyplotph_diode_tr.log

帖子发出去，基本都会沉啊，不给力