Continuous dielectric permittivity II: An Iterative Method for   Calculating the Polar Component of the Molecular Solvation Gibbs Energy Under   a Smooth Change in the Dielectric Permittivity of a Solution
  An iterative method for calculating the polar component of the solvation Gibbs energy under a smooth change in dielectric permittivity, both between a substrate and a solvent and in a solvent is formulated on the basis of a previously developed model. The method is developed in the approximation of the local relationship D = \eps (r) E between the displacement vectors D and the electric field intensity E. 
F.V. Grigoriev, Oleg Kupervasser, I.P. Kikot'
physics.chem-ph