The most significant contributor to this difference is probably the result of a change in the way force constants are analytically evaluated for solvation models, in which the continuum cavity in G16 is now evaluated over a smoothed rather than the original tessellated -and hence discontinuous -grid used in G03. The smoothed cavity approach will be adopted for all the systems for which kinetic isotope effects are computed in this work.