and DFT-D4 approaches was accounted for by the third-order Axilrod-Teller-Muto (ATM) term. In DFT-D4, partial atomic charges were determined using the classical electronegativity equilibration procedure. We note that, the rational BJ damping function was applied in all DFT-D4 calculations even though, as suggested by the author of DFT-D4 paper, the BJ abbreviation is not included in the name of the method.