In the present paper, kinetic isotope effects of the title reaction are studied with canonical variational transition state theory on the modified Wang Bowman (MWB) potential energy surface (PES) (Chem Phys Lett 2005, 409, 249) and the ab initio calculations at the quadratic configuration interaction (QCISD (T, full))/aug-cc-pVTZ//QCISD (full)/cc-pVTZ level. The calculated rate constants for the isotopic variants of this title reaction on the MWB PES have good agreement with those of the present ab initio calculations over the temperature range of 20-5000 K for the forward reactions and 800-5000 K for the reverse reactions, respectively. In particular, the forward rate constants for the title reaction and its isotopically substituted reactions have negative temperature dependences at about 40 K. Rate expressions are presented for all the studied reactions. C