Derivation and implementation of multi-dimensional transition state theory to accurately determine the rate constants of rare events