Reply to comment on “A computational study of the structures of van der Waals and hydrogen-bonded complexes of ethene and ethyne”

We have previously found that current density functionals do not predict bond energies, normal mode frequencies and distortion constants consistent with MP2 calculations for weakly bonded complexes. Jeong and Han report calculations performed using extended basis sets for two of our complexes. Their results underscore our previous conclusions and show that it is necessary to use rather large basis sets in order to compare the relative accuracies of BLYP and B3LYP for weakly bonded systems.