As already emphasized, the D3S scheme (like D3 itself) is independent of the choice of damping function to remove double counting and incipient divergence at short range. Since the D3S scheme only slightly affects the elements that compose most of the binding interaction datasets, reoptimizing the functional-dependent parameters, such as s8, sr6 for zero damping, or α1 and α2 for Becke-Johnson damping, should be unnecessary since optimized D3 parameters should also be virtually optimal for the D3S scheme. To validate this expectation, we reoptimized the s8 and sr6 coefficients for the D3(0) and D3S(0) schemes for the B3LYP functional using training sets described in Ref. . We employed a grid approach to find the best combinations of s8 and sr6 coefficients using a product of product of root-meansquare errors as a sorting criterion. The results show no difference between the D3(0) and D3S(0) optimized coefficients with a grid step size of 0.001.