Monte Carlo simulation of confined fluids of polarizable particles: an efficient iterative treatment of the local field in slab geometry using Ewald summation

We propose a method for treating in Monte Carlo simulations the problem of the induced dipoles for polarizable particle fluids confined in slab geometry and subject to an external field. In order to compute the local field in a reasonable time, a partial update of the induced dipole moments is performed by introducing a cut-off distance, as in bulk systems. This strategy is then combined with a slab adapted 3D-Ewald summation for treating the long-range interactions between the induced dipoles. The method is illustrated by simulations of confined binary mixtures in the canonical and grand canonical ensembles.