Simulating nanoscale dielectric response
  We introduce a constrained energy functional to describe dielectric response. We demonstrate that the local functional is a generalization of the long ranged Marcus energy. Our re-formulation is used to implement a cluster Monte Carlo algorithm for the simulation of dielectric media. The algorithm avoids solving the Poisson equation and remains efficient in the presence of spatial heterogeneity, nonlinearity and scale dependent dielectric properties. 
A.C. Maggs, R. Everaers