will be used for the calculation of the total energy. The last term is the van der Waals (vdW) interaction between QM and MM regions. Here we adopt the distance-dependent variance of the classical Lennard-Jones (LJ) 12-6 potential according to Jensen et al. The expression is presented in Section S1 of the supporting information (SI).