Finally, we wanted to understand to what extent our results might depend on our choice of DFT method and basis set. To that end, we tested the method/basis set combination of PW6B95/def2-TZVPD. The PW6B95 functional is a hybrid meta-GGA functional that is considerably newer than B3LYP and, along with SOGGA11-X , was recently found to have the most accurate gas phase dipole moment predictions of any non-double hybrid functional relative to CCSD(T)/CBS values for a benchmark set of 152 molecular species. Indeed, the accuracy of PW6B95 exceeded that of restricted MP2 in this previous study. The def2-TZVPD augmented basis set differs meaningfully from aug-cc-pV(T+d)Z in terms of the philosophy used in its construction and is more minimally augmented. Together with a PCM model of benzene, this combination yields dipole moments that are 1.178 times larger (95% CI: [1.157, 1.207]) than experimental gas phase (Figure ). This is very similar to the B3LYP/aug-cc-pV(T+d)Z result of 1.182, and the same is true for the PCM water results as well. It is interesting to note that two quite different combinations of DFT method and basis set yield such similar predictions for the characteristic overpolarization. Moreover, the variation of data with respect to the regression line (i.e., the weighted residuals) and the fractional overpolarization for this combination with both PCM benzene and PCM water are quite similar to the corresponding B3LYP/aug-cc-pV(T+d)Z results (Figures and). These results lend credence to the notion that many modern density functionals in combination with suitable augmented basis sets could provide a reasonable foundation for these calculations.