In concordance with the considerable variation of data about the regression line displayed in Figures and, the weighted residuals of the HF/6-31G* method show the broadest distribution of all four method/basis set combinations that we examined. In particular, the maximum and median absolute weighted residuals for HF/6-31G* are 0.261 and 0.052, respectively (Figure ). For comparison, the maximum and median absolute weighted residuals for HF/aug-cc-pV(T+d)Z are 0.135 and 0.040, respectively (Figure ). This complements our earlier observation regarding the 95% CI for the slope estimate for this method vs. HF/6-31G* and suggests that the use of a more sophisticated basis set is indeed useful in improving the consistency of the overpolarization predictions, regardless of the underlying electronic structure method employed. The predictions of B3LYP/aug-cc-pV(T+d)Z with PCM benzene are even more consistent, however, with maximum and median absolute weighted residuals of 0.104 and 0.036, respectively (Figure ). In other words, the maximum absolute weighted residual for this combination is approximately 2.5 times smaller than for HF/6-31G*, and the median absolute weighted residual is 30% smaller. There is a reversal of this trend for B3LYP/aug-cc-pV(T+d)Z