We also consider the phase correction proposed by Shenvi and co-workers, which introduces a new effective Hamiltonian to correct the relative phase between Gaussian WP components on different PESs. For two-level systems, the relative phase caused by the distinct motions of Gaussian WP components can be introduced by replacing the diagonal elements of the adiabatic Hamiltonian, 71