Silvaco 仿真二极管反向恢复电流

go atlas
#
mesh
x.mesh loc=0.00 spac=0.5
x.mesh loc=25 spac=0.5

y.mesh loc=-0.5 spac=0.1
y.mesh loc=0 spac=0.5
y.mesh loc=10 spac=0.5
y.mesh loc=35 spac=1
y.mesh loc=40 spac=2
#
region num=1 y.max=0 oxide
region num=2 y.min=0 y.max=40 silicon
#
# #1=anode #2=cathode
electrodenum=2 name=anodex.min=0 length=2.5 y.max=0
electrodenum=2 name=anodex.min=22.5 length=2.5 y.max=0
electrodenum=1 name=anode1x.min=8 length=9 y.max=0
electrode bottomname=cathode

#
dopingreg=2uniform conc=8e14n.type y.min=0 y.max=40
dopingreg=2gauss conc=5e19p.type junction=3 x.min=0 x.max=5 ratio=0.75
dopingreg=2gauss conc=5e19p.type junction=3 x.min=20 x.max=25 ratio=0.75
#
dopinggauss conc=2e20n.type peak=40 junc=35
save outf=mps25_0.str
tonyplot mps25_0.str -mps25_0.set
models srh conmob bgn auger fldmob
impact crowell
material reg=2 taun0=5e-6 taup0=2e-6
contact name=anode1 workf=4.93
solve init
solve
solve vcathode=-0.1
method newton
log outfile=mps251.log
solve vanode=0
solve vcathode=-0.05 vstep=-0.05 vfinal=-0.5 name=cathode
save outf=mps25_1.str
################################################################
go atlas
.begin
#
#Steady-state simulation of circuit with power diode
#
v1 1 0150.
r1 1 21m
l1 2 33uH
adiode3=cathode 4=anode 4=anode1 width=2.e6 infile=mps25_1.str
r2 4 01mg
i1 0 410.
#
.nodeset v(1)=150. v(2)=150. v(3)=150. v(4)=150.5
.numeric vchange=0.1
#
.save outfile=mps25_1_save
.options m2ln print
.end
#
models device=adiode reg=2 conmob fldmob consrh auger bgn
material device=adiode reg=2 taun0=5e-6 taup0=2e-6
#
impact device=adiode reg=2 selb

go atlas
.begin
#
#Reverse recovery of power diode
#
v1 1 0200.
r1 1 21m
l1 2 31uH
adiode3=cathode 4=anode 4=anode1 width=2.e6 infile=mps25_1.str
r2 4 01mg EXP 1mg 1e-3 0. 20ns 10 200
i1 0 410
#
.numeric lte=0.3 toltr=1.e-5 vchange=10.
.options print relpot write=10
#
.log outfile=mps25_1
.load infile=mps25_1_save
.save master=mps25_1
#
.tran 0.1ns 5us
#
.end
models device=adiode reg=2 conmob fldmob consrh auger bgn
material device=adiode reg=2 taun0=5e-6 taup0=2e-6
impact device=adiode reg=2 selb

go atlas
tonyplot mps25_1_tr.log -set mps25_1.set
quit
在仿真反向恢复电流的时候为什么是这样的，电流没有变为0，而是一直维持在反向最大的状态? 大神帮看看那个地方有错。
在仿真反向恢复电流的时候为什么是这样的，电流没有变为0，而是一直维持在反向最大的状态? 大神帮看看那个地方有错。

图像是这样的

我也期待有高手解答呢！