Coupled cluster calculations with singles, doubles, and non-iterative perturbative triples using local approximations (LCCSD(T0)) were performed with the cc-pVQZ 47 basis set. Orbitals were localised according to the Pipek-Mezey 48 scheme. Density fitting approximations were used throughout. NPA and Chelpg 54 charges were used for the electrostatic embedding. Force field parameters for the LJ potentials were taken from the CHARMM Generalised Force Field (CGFF) and the MM3 56 force field. Clusters were generated using a cutoff of 21 Å (40 a.u.), in-cluding 108 background charges at the cluster surface.