and r c,a in the latter is the covalent radius of the ath MM atom. The modified interaction tensor is introduced in order to mitigate divergences at short distances thus avoiding the electron spill-out problem which is pervasive in electrostatic embedding QM/MM due to the lack of Pauli repulsion between the QM and MM subsystems. The modified form of the interaction potential was introduced and validated in Ref. 4.