Interestingly, after applying the bond-based correction, the MAE of all three functionals falls to a similar level of errors irrespective of system size. The trend of the uncorrected and the corrected line of the B3LYP-D functional in Figure exhibits a behavior similar to the trend in formation energies of molecules with B3LYP vs G3 methods as studied by Wodrich et al. The fact that our method can capture a similar trend for reaction energy without performing Gn-level calculations is very interesting. The error reduction of all the post-SCF functionals to the same level of MAE after the correction implies that the accuracy of the correction method does not depend on the post-SCF functionals used. Therefore, in frameworks where a large number of energy calculations are needed, e.g. potential energy surface exploration, the combination of cheap functionals with bond-based correction would save computational time and at the same time provide better accuracy than the higher-level DFT methods.