Here, a new implementation for calculating canonical variational rate constants is introduced and tested against methanol dissociation rate constants and the (high pressure) radical recombination reaction, CH 3 þ OH ! CH 3 OH. Reaction paths are described at the CASSCF level, with energy corrections at MRMP2. Molecular properties for both stationary and non-stationary points along the reaction path are used for the calculation of rate constants. Different models for the low vibrational frequencies are discussed: harmonic oscillator and free rotor models. The CH 3 þ OH ! CH 3 OH recombination rate constants are determined from dissociation rate constants and equilibrium constants. Our calculated rate constants agree with previously reported data for these reactions, suggesting the good performance of our implemented code in calculating canonical variational rate constants. The rate expressions are suggested: as: k dis (T) ¼ 8.42 Â 10 þ16 Â e À96.18/RT and k rec (T) ¼ 2.05 Â 10 À10 Â (T/298) À0.24 Â e À0.30/RT , in units s À1 and cm 3 molecule À1 s À1 , respectively, and over the temperature range 300-2000 K, with activation energies in kcal/mol.