Iterative integral equation methods for structural coarse-graining

In this paper, new Newton and Gauss-Newton methods for iterative coarse-graining based on integral equation theory are evaluated and extended. In these methods, the potential update is calculated from the current and target radial distribution function, similar to iterative Boltzmann inversion, but gives a potential update of quality comparable with inverse Monte Carlo. This works well for the coarse-graining of molecules to single beads, which we demonstrate for water. We also extend the methods to systems that include coarse-grained bonded interactions and examine their convergence behavior. Finally, using the Gauss-Newton method with constraints, we derive a model for single bead methanol in implicit water, which matches the osmotic pressure of the atomistic reference. An implementation of all new methods is provided for the open-source VOTCA package.