Operator-Split Method for Variance Reduction in Stochastic Solutions of the Wigner Equation

M. Nedjalkov 1 , E. Atanassov 2 , H. Kosina 1  and S. Selberherr 1
  • 1 Institute for Microelectronics, TU-Vienna
  • 2 CLPP, Bulgarian Academy of Sciences, Sofia, Bulgaria

The standard stochastic approach for simulation of carrier transport described by the Wigner equation introduces hard computational requirements. Averaged values of physical quantities are evaluated by means of numerical trajectories which accumulate statistical weight. The weight can take large positive and negative values which introduce large variance of the calculations. Aiming at variance reduction, we utilize the idea to split the weight so that a part is assigned to the trajectory and a part is left on a phase space grid for future processing. Formally this corresponds to a splitting of the kernel of the Wigner equation into two components. An operator equation is derived which couples the two kernel components and gives an answer how to further process the stored weight. The obtained Monte Carlo algorithm resembles a physical process of generation or annihilation of positive and negative particles. Variance reduction is achieved due to the partial annihilation of particles with opposite sign. Simulation results for a resonant-tunneling diode are presented.

Purchase article
Get instant unlimited access to the article.
Log in
Already have access? Please log in.

Log in with your institution

Journal + Issues

This quarterly published journal presents original articles on the theory and applications of Monte Carlo and Quasi-Monte Carlo methods. Launched in 1995 the journal covers all stochastic numerics topics with emphasis on the theory of Monte Carlo methods and new applications in all branches of science and technology.