During QCT evolution process and in the final products, the vibrational energy of the system may be lower than the zero point vibrational energy (ZPE). The so called hard-ZPE constraint is applied. Namely, if the vibrational energy of each product is less than its ZPE, then this trajectory is not considered in the final statistics. Further, anharmonicity effects may affect the ZPE. To this end, the scale factor, 0.96. was obtained from the two reactants, OH and CH3OH, and the corresponding products, CH2OH, H2O, and CH3O, according to the formula, 96 demonstrate that the channel of CH2OH + H2O always dominates over this temperature range. As shown, the hard-aZPE results are consistent with and slightly higher than the theoretical results of Ocana et al., Shannon et al., at 300 ~ 400 K and 425 ~ 900 K, respectively, and the experimental results by Liu et al. at 900 ~ 1000 K. To explore the reaction mechanism, several movies of trajectories are provided in SI. The Ring-polymer Molecular dynamics method (RPMD) or other methods, which can well describe the quantum effects, could be used in the future studies of kinetics below 300 K or even at 10 K on this newly fitted PIP-NN PES to explore the temperature dependences of rate coefficients.