For small QM regions, the reaction energy shows large oscillations, while starting from QM Overall, our results confirm the slow convergence of the reaction energy with increasing size of the QM region found in earlier studies and emphasize the need for systematic protocols for the construction and selection of the QM region. We note that Jindal and Warshel 14