“…However, these mechanisms are independent of the QM treatment, and the investigation of such contributions can be effectively performed by coupling our model to other mechanistic treatments, such as polarizable force fields or advanced docking techniques. In addition, we here employ a common, well-established DFT approximation (PBE+D3), which already provides useful information for coarse-grained quantities and trends [21,35], and aptly fits the aim of intensive, high-throughput simulations of several structures, especially in their relaxed positions [45]. Yet, we leave out arguments about the choice of other ab initio levels of theory, which may shed light on more quantitative aspects related to processes beyond the ground-state: reaction coordinates, activation barriers, etc.…”