The reparameterized D3S scheme (see the main text for details on the scheme definition) was implemented and tested in the Q-Chem 6.2 code. Q-Chem was also used to optimize the problematic cases where unphysical local minima appeared with the standard D3 schemes. These geometries were reoptimized at the B3LYP-D3S(BJ)/def2-SVP level of theory. A modified version of the stand-alone D3 code distributed by the Grimme group is also available for download to enable other codes to be readily updated to use the D3S method, or to validate independent implementations. A set of 178 atom-molecule and molecule-molecule pairs, including main group chemistry species ranging from organic hydrocarbons to molecules like SiH4 and SF6, was used for comparison with experimentally obtained molecular C6 coefficients. Twenty-three molecules in this set were optimized at the B3LYP/6-31G(2df,2p) level of theory, similar to the original paper. All structures were considered in the singlet spin state, except for the O2 molecule, which was set to a triplet spin state. Combinations between these molecules, as well as with nine single atoms (H, He, Li, N, O, Ne, Ar, Kr, Xe), resulted in 178 different molecular pairs. The molecular C6 coefficients were computed by summing atomic C6 AB coefficients between all intermolecular atomic pairs. A table with all the computed molecular C6 coefficients is provided in the ESI (Table ). We note, that molecular pair C6 coefficients for D3 scheme are independent of the damping function and only depend on how the C6 AB are defined (i.e. in our case results for D3 and D3S will be different).