The rigorous expression to calculate quantum rate constants in the flux correlation function formalism has been well established. For a single molecule with the symmetric double well potential as described in Sec. II A, we assume that the reactant and product region are separated by a dividing surface at x ds i = 0. The population of the ith molecule in the product (right well) and the reactant (left well) region are ob-tained as