General Linear Scaling Implementation of Polarizable Embedding Schemes

A general framework to treat polarizable embedding schemes in the context of QM/MM calculations is presented. Such a framework is completely general, as it allows in principle one to treat any electrostatic distribution and polarization model with minimal modifications and achieves linear scaling in computational cost and memory requirements. The performances and scaling of the new implementation are demonstrated with benchmark calculations on water clusters including up to 1146240 atoms.