6533b828fe1ef96bd12884cf

RESEARCH PRODUCT

Iterative integral equation methods for structural coarse-graining

Marvin P. BernhardtNico F. A. Van Der VegtMartin Hanke

subject

Quantitative Biology::BiomoleculesMonte Carlo methodGeneral Physics and AstronomyInverseRadial distribution functionIntegral equationInversion (discrete mathematics)symbols.namesakeBoltzmann constantConvergence (routing)symbolsApplied mathematicsGranularityPhysical and Theoretical ChemistryMathematics

description

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.

https://doi.org/10.1063/5.0038633