“…Looking at final expression ͑B7͒, rewritten here as W k ͑r ជ , t͒ = W mean ͑r ជ , t͒ − V͑r ជ , r ជ k ͓t͔͒, it seems that W k ͑r ជ , t͒ can be computed from a unique mean-field solution of the Poisson equation W mean ͑r ជ , t͒ when subtracting the appropriate two-particle potential V͑r ជ , r ជ k ͓t͔͒. In fact, this is the strategy used in several recent works [28][29][30][31][32] to improve the treatments of the short-range Coulomb interaction in electron device Monte Carlo simulators. However, this strategy is not as general as our procedure because these works need an analytical expression for the two-particle Coulomb interaction V͑r ជ , r ជ k ͓t͔͒.…”