QM/MM boundaries across covalent bonds: A frozen localized molecular orbital-based approach for the effective fragment potential method

A computational methodology to treat the covalent boundary between QM and M M regions in the Effective Fragment Potential (EFP) method, by defining a b uffer region consisting of frozen localized molecular orbitals (LMOs), is i ntroduced. The implementation of energy, gradient, and EFP parameter evalua tions in the presence of frozen LMOs is discussed. The magnitude and source of errors introduced by various choices of buffer region is studied for th e proton affinities of lysine and the tripeptide glycine-lysine-glycine. It is shown that by reasonable choice of the frozen density buffer region the proton affinity error can be consistently decreased to less than 0.5 kcal/ mol compared to the full ab initio calculation.