During GaN vertical power device simulation, faces convergence problem at breakdown point. I used voltage, current boundary condition and curvetrace method. But can not get stable solution beyond 1000V. I also tuned the methods parameter for very low bias point increment like, method newton climit=1e-9 maxtrap=40 itlimit=40 dvmax=1e8 ir.tol=1e-40 cr.tol=1e-40 ix.tol=1e-40 px.tol=1e-30 pr.tol=1e-45 cx.tol=1e-30.
I also try high simflags for numerical precision like ~160bit, 256bit precision
During 256 bit precision get very long simulation time, although used all available system like "-P all"