We employed here the latter, less costly procedure as our tests on the model system indicated that both procedures obtain comparable interaction energies within chemical accuracy (see section S1.1, table ). The short-range attenuation parameter b of the studied functionals was taken from previous works benchmarking the S22 and S66 datasets. On the other hand, as the optimization of the long-range parameter C yield only slight improvement on the aforementioned datasets, its original value (C = 0.0093) is hence used here.