Assuming an initial Gaussian density of the form of Eq. ( ), the initial phase is given by S(x, t = 0) = p 0 (xx c ), where x c < 0, and p 0 > 0 is the initial momentum, so that the wavepacket is incident from the left side of the PES barrier. Unless explicitly stated, we use atomic units henceforth.