“…In geometry optimizations and frequency calculations, the high level (HL) is typically described using a hybrid DFT functional, while the low level (LL) is treated by semi-empirical methods, e.g., AM1 [83], MNDO [84] or PM3 [85], or a molecular mechanics method, e.g., the universal force field (UFF) [86]. Some authors allow the entire cluster model to relax during geometry optimizations, keeping the terminating hydrogens fixed in space to prevent the cluster from collapsing [87,88], while others only relax the RR containing the active center [52,[89][90][91]. Following geometry optimizations at the ONIOM level, single-point energy refinements are sometimes performed on the entire cluster to minimize potential artifacts arising from the specific partitioning of the system [76,87,88,92].…”