Exact Correction of the Self-Force Problem in Monte Carlo Device Simulation

The self-force is a specific problem of self-consistent Monte Carlo-Poisson simulation resulting in an un-physical field component acting on a particle coming from the particle itself (the self-force). Several approaches have been proposed in literature to mitigate this problem, but all of them suffer to some extent of approximations and/or limitations. In this paper we propose a new and mathematically exact correction of the self-force problem based on a numerical approach. Although computationally expensive, it has no restriction and can be always applied. The new method has been tested on the difficult problem of plasma oscillation simulation providing the expected plasma energy from theory. Moreover, the same mathematical framework introduced here for the self-force correction can be readily applied also for the exact calculation of the reference force in the Particle-Particle-Particle-Mesh (P3M) method. The accuracy of such approach to P3M method is demonstrated by simulating the bulk low field mobility dependence on doping concentration.