Harmonic restraints also lead to more rapid convergence of the ZPE correction, which is converged even for model III. This is a useful observation because DFT analytic frequency calculations incur a cubic-scaling memory bottleneck for solution of the coupled-perturbed equations. With harmonic restraints, it may be possible to use smaller models to estimate the thermal vibrational contribution to ∆ rxn G.