做了一个SIC二极管的例子

shg841100

导通特性能出来可是击穿特性曲线出不来！
File{
     grid   ="n1_mdr.grd"
     doping ="n1_mdr.dat"
     plot   ="n2_des.dat"
     current="n2_des.plt"
     output ="n2_des.log"
     parameter = "pp2_des.par"
}
Electrode {
   { Name="A"  Voltage= 0.0 Resist= 1e13 }
   { Name="K" Voltage= 0.0 }
}

Physics {
   Recombination (
      SRH(DopingDependence TempDependence)
      Auger
      Avalanche (Okuto Eparallel)
   )
   OptBeam(
      WaveInt= 1e12    # [1/cm/s]
      SemAbs= 1e0     # [1/cm]
      SemSurf= -0.4 # [cm]
      SemWindow= ( 0.0 90.e-4 ) # [cm]
   )
   Mobility (
      DopingDependence
      HighFieldSaturation
   )
   EffectiveIntrinsicDensity ( oldSlotboom NoFermi )
   Temperature= 300
}
Plot {
   eDensity hDensity
   eCurrent hCurrent
   ElectricField
   eQuasiFermi hQuasiFermi
   egradQuasiFermi hgradQuasiFermi
   Potential Doping SpaceCharge
   SRH Auger
   AvalancheGeneration
   eAvalanche hAvalanche
   eMobility hMobility
   DonorConcentration AcceptorConcentration
   Doping
   eVelocity hVelocity
   BarrierTunneling
   ConductionBandEnergy ValenceBandEnergy BandGap
   eQuasiFermi hQuasiFermi
   NonLocal
}
Math {
   Extrapolate
   Digits= 7
   ErrEff(electron)= 1e-12
   ErrEff(hole)= 1e-12
   RHSmax= 1e30
   RHSmin= 1e-30
   CdensityMin= 1.e-30
   Notdamped= 20
   Iterations= 15
   RecBoxIntegr
}
Solve {
*- Creating initial guess:
   Coupled(Iterations=100 LineSearchDamping=1e-4){ Poisson }
   Coupled { Poisson Electron }
*- Ramp to p_side  
   Quasistationary(
      InitialStep=1e-1 Increment=1.35
      MinStep=1e-5 MaxStep=0.05
      Goal { Name="K" Voltage=700.0 }
   ){ Coupled { Poisson Electron } }

      CurrentPlot( Time=(Range=(0 1) Intervals=20)  )
   }