Support

Support Options

Submit a Support Ticket

 

Questions and Answers

Status: Accepting answers

0 Like

Youngseok Kim

Issue in obtaining solution of Poisson eq. for self-consistent calculation in NEGF

Hi,

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)

dN=n-Nd+((1/beta)*D2*U); dU=(-beta)*(inv(D2-(beta*diag(D))))*dN;U=U+dU;

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

Report abuse

0 Responses

No answers to this question yet. Be the first to answer this question.

Did you know you can earn points for providing good answers?
Learn more about how points are awarded.

nanoHUB.org, a resource for nanoscience and nanotechnology, is supported by the National Science Foundation and other funding agencies. Any opinions, findings, and conclusions or recommendations expressed in this material are those of the author(s) and do not necessarily reflect the views of the National Science Foundation.