Lack of convergence - thermodynamic symulation is Sentaurus Device

File {
   * input files:
   Grid="3D_N_r40o4_s4e+17_p5e+18_msh.tdr"
   Parameter="pp51_des.par"
   * output files:
   Plot=   "n51_des.tdr"
   Current="n51_des.plt"
   Output= "n51_des.log"

   NewtonPlot = "n51_des_newton.tdr"
}


Electrode {
	{ Name="Source" 		Voltage=0.0 Resistor= 40 }
	{ Name="Drain"  		Voltage=0.0 Resistor= 40 }
	{ Name="Gate1" 		Voltage=0.0 }
	{ Name="Gate2" 		Voltage=0.0 }
	{ Name="SourceBot" 	Voltage=0.0 }
	{ Name="DrainBot"  	Voltage=0.0 }
	{ Name="Gate1Bot" 	Voltage=0.0 }
	{ Name="Gate2Bot" 	Voltage=0.0 }
}

Thermode {
	{ Name = "Source"	Temperature = 300 SurfaceResistance=2e-4 }
	{ Name = "Drain"	Temperature = 300 SurfaceResistance=2e-4 }
	{ Name = "Gate1"	Temperature = 300 SurfaceResistance=2e-4 }
	{ Name = "Gate2"	Temperature = 300 SurfaceResistance=2e-4 }
}

Physics{
   Thermodynamic
   RecGenHeat
   Hydrodynamic( eTemperature
 )  
   EffectiveIntrinsicDensity(BandGapNarrowing ( OldSlotboom ))
}

Physics(Material="Silicon"){
   eQuantumPotential

   Mobility(
	   DopingDep
      PhuMob
      eHighFieldSaturation
      Enormal
   )
	Recombination(
    	SRH( DopingDep TempDependence )
  )

}

Insert= "PlotSection_des.cmd"
Insert= "MathSection_des.cmd"

Solve {
*- Creating initial guess:
   Coupled(Iterations= 100 LineSearchDamping= 1e-4){ Poisson eQuantumPotential
 } 
   Coupled { Poisson eQuantumPotential Electron
 }
   Coupled { Poisson eQuantumPotential Electron eTemperature Temperature
 }

*- Ramp to drain to Vd
   Quasistationary( 
      InitialStep= 1e-2 Increment= 1.35 
      MinStep= 1e-5 MaxStep= 0.2 
      Goal { Name="Drain" Voltage=1.25
} 
   ){ Coupled { Poisson eQuantumPotential Electron eTemperature Temperature
 } } 

 
*- Vg sweep 
   NewCurrentFile="IdVg_" 
   Quasistationary( 
      DoZero 
      InitialStep= 1e-3 Increment= 1.5 
      MinStep= 1e-6 MaxStep= 0.04 
      Goal { Name="Gate1" Voltage=1.25
 } 
   ){ Coupled { Poisson eQuantumPotential Electron eTemperature Temperature
 } 
      CurrentPlot( Time=(Range=(0 1) Intervals= 30)  )
   }
}

Math {
   Extrapolate
   Notdamped= 100
   Iterations= 25
   ExitOnFailure

   Number_Of_Threads= maximum
	
	Method=ILS
	
   *Cnovergence
   ExtendedPrecision
   #IncompleteNewton
   CNormPrint
	
}

Plot{
*--Density and Currents, etc
   eDensity hDensity
   TotalCurrent/Vector eCurrent/Vector hCurrent/Vector
   eMobility hMobility
   eVelocity hVelocity
   eQuasiFermi hQuasiFermi

*--Temperature 
   eTemperature Temperature hTemperature
   RecombinationHeat eJouleHeat hJouleHeat
   lHeatFlux
   TotalHeat

*--Fields and charges
   ElectricField/Vector Potential SpaceCharge

*--Doping Profiles
   Doping DonorConcentration AcceptorConcentration

*--Generation/Recombination
   SRH Band2Band * Auger
   * AvalancheGeneration eAvalancheGeneration hAvalancheGeneration

*--Driving forces
   eGradQuasiFermi/Vector hGradQuasiFermi/Vector
   eEparallel hEparallel eENormal hENormal

*--Band structure/Composition
   BandGap 
   BandGapNarrowing
   Affinity
   ConductionBand ValenceBand
   eQuantumPotential hQuantumPotential

*--Gate Tunneling
   * eBarrierTunneling hBarrierTunneling  BarrierTunneling
   * eDirectTunnel hDirectTunnel
}