Issue in obtaining solution of Poisson eq. for self-consistent calculation in NEGF
I’ve been working on an exercise matlab code posted by Prof. S. Datta. ( https://nanohub.org/resources/19564 ) In that example, the code performs 1D simple self-consistent calculation btw NEGF and Poisson eq. It uses Newton-Raphson method to get a solution of Poisson eq. as below and it converges smoothly after few iterations. (refer Appendix B in linked thesis for detail of Newton-Raphson method: https://engineering.purdue.edu/gekcogrp/publications/theses/PhD_11_2011_Sunhee_Lee_PhD_Thesis_main.pdf)
However, if I get a solution from direct inversion as below, it doesn’t converges at all.
N=n-Nd; U =(-beta)*inv(D2)*N;
I found that in a simple case like PN junction, direct inversion works well. Why it is not working in this case? In what case I need to work with more numerically stable method such as Newton-Raphson method?
I will be appreciated if you can help me.
Thanks. Youngseok Kim