“…However, for large chemical systems, such as those encountered in heterogeneous catalysis, where adequate models typically include several hundred atoms, the “chemical accuracy” (4 kJ mol −1 for the energy barriers, one order of magnitude for the pre‐exponential factors) required for useful predictions is not attained. The exponential scaling with the system size of both the electronic structure methods (potential energy surface) and the nuclear motion problem (vibrational partition function) limits the applicable method to density functional theory (DFT) for the potential energy surface and the harmonic approximation for the vibrations 4, 5, 6, 7, 8. Depending on the functional, energy barriers may be in error by 10–20 kJ mol −1 , and the harmonic approximation is most problematic for low‐frequency modes, which are known to make the largest contribution to the partition functions 9…”