A directab initiodynamics approach for calculating thermal rate constants using variational transition state theory and multidimensional semiclassical tunneling methods. An application to the CH4+H↔CH3+H2reaction