A combined quantum-classical dynamics method for calculating thermal rate constants of chemical reactions in solution