Simulations of Coulomb systems confined by polarizable surfaces using periodic Green functions

We present an efficient approach for simulating Coulomb systems confined by planar polarizable surfaces. The method is based on the solution of Poisson equation using periodic Green functions. It is shown that the electrostatic energy arising from surface polarization can be decoupled from the energy of periodic replicas. This allows us to combine an efficient Ewald summation method for the replicas with the polarization contribution calculated using Green function techniques. We apply the method to calculate density profiles of ions confined between charged dielectric and metal interfaces.