The leading contributions to reaction kinetics from vibrational anharmonicity are of three types: anharmonic ZPE due mainly to high-frequency modes, multiple-structure anharmonicity due to multiple conformers, and torsional potential anharmonicity. For polyatomic reactions, the ZPE is large, and therefore a small change in the ZPE can have a large effect on the predicted rate constants. For inclusion of anharmonic ZPE in CTST and CVT/SCT, a practical method is harmonic analysis with reactant frequencies scaled by a factor optimized against stable molecules and with scale factors for transition structures determined by HDCPT2. Although HDCPT2 can be problematic for large-amplitude vibrations, the ZPE is dominated by high-frequency modes for which it should be acceptable. However, this issue does not arise in the present CNEO-TST calculations because the ZPE is included via the CNEO-DFT calculation, which does not make the harmonic approximation.