Short-range and long-range Coulomb interactions for 3D Monte Carlo device simulation with discrete impurity distribution