Medici 仿真的隧穿结I-V特性曲线不正确，为什么？

MESH RECT
X.MESH X.MIN=-1 X.MAX=1 N.SPACES=20
Y.MESH L=0 N=1
Y.MESH L=0.5 N=11
REGION Silicon
ELECTR NAME=A LEFT
ELECTR NAME=C RIGHT
PROFILE UNIFORM P-TYPE N.PEAK=5E23 X.MIN=-1 X.MAX=0 OUT.FILE=GRID
PROFILE UNIFORM N-TYPE N.PEAK=5E23 X.MIN=0 X.MAX=1
PLOT.2D GRID TITLE="INITIAL GRID" FILL SCALE
REGRID DOPINGRATIO=2 SMOOTH=1 IN.FILE=GRID
PLOT.2D GRID TITLE="DOPING REGRID" FILL SCALE
MODELS SRH AUGER CCSMOBBTBT INCOMPLE IMPACT.I R.TUNNEL
SYMBOL NEWTON CARRIERS=2
METHOD AUTONR NRCRITER=0.1
SOLVE V(A)=0 V(C)=0
PLOT.1D CONDUC BOTTOM=-2 TOP=2
+ TITLE="BAND"
PLOT.1D VALENC color=4 UNCHAN
PLOT.1D QFN color=2UNCHAN
SOLVE
SYMBOL NEWTON CARRIERS=2
METHOD AUTONR NRCRITER=0.1
SOLVE V(C)=0 ELEC=A V(A)=1 VSTEP=0.01 NSTEP=50
PLOT.1D X.AXIS=V(A) Y.AXIS=I(A)
+ TITLE="I-V"
PLOT.1D CONDUCBOTTOM=-2 TOP=4
+ TITLE="CHANGED BAND"
PLOT.1D VALENCUNCHAN
PLOT.1D QFN COLOR=4 UNCHAN
仿出的I-V特性曲线是错的，想请大家帮忙讨论一下是为什么，还有，画出的BNAD图的纵坐标是电势，而不是势能，我想画的是能带图，请问怎么画?不胜感激~

你应该把图放上来, 5E23这样高的pn结是做不出来的。