“…On the other hand, far enough from them (∼ tens of nanometers), the shape of the scalar potential depends on the particular spatial charge distribution. In this regard, a common strategy consists on introducing a long-range part of the Coulomb interaction numerically through the solution of the mesh-dependent mean-field Poisson equation, and a short-range part analytically (the so-called Monte Carlo/Molecular dynamics approach, MCMD, (Gross et al, 1999;Wordelman & Ravaioli, 2000)), or perturbativelly (Barraud et al, 2002). The MCMD approach, however, shows some inconvenients.…”