I am simulating a PIN device in Sentaurus Workbench (TCAD), and my simulation results show a huge electric field peak. (I apply a negative DC voltage of 1300V to see if the PIN rectifier can withstand it).
I have been trying to model the real device better (try to get a reduced electric field distribution in my simulation) by modifying the resistance of certain regions (ion implanted regions) but my simulation has problems converging whenever I change any parameters (the resistance of regions or if I add interface charge to it). I can add and modify the slightest thing and then it will not converge.
Would you have any suggestions for me as to how to make my solutions converge? I am using Quasistationary Newton convergence algorithm. Using simple Poisson coupling. As of now, I am playing around with the maximum step size and initial step size in the SDEVICE code, but that does not seem to work.
Thanks a lot!