The B3LYP/6-31G(df,3p) model for the calculation of deuterium nuclear quadrupole coupling constants (nqcc's) is shown to yield results as accurate as calculations previously performed at the MP4 level of theory. For 25 molecules, ranging from HD and DF to pyridine and fluorobenzene, the rms difference between the B3LYP nqcc's and the experimental nqcc's is 3.2 kHz (2.7%). For benzene, our calculations suggest that the experimental χbband χccof S. Jans-Bürli, M. Oldani, and A. Bauder, 1989.Mol. Phys.,68, 1111–1123) have been incorrectly assigned with respect to inertia axes and should be reversed. For borane carbonyl and nitric acid, it is shown that nqcc calculations using hydrogen bond lengths given by MP2/6-311 + G(d,p) optimizations in combination with the heavy atom experimental structures significantly improve agreement with the experimental nqcc's.