is a Gaussian wavepacket centered around R 0 = -0.7 a.u. with variance 1/2M ω R , mass M = 550 Da, and frequency ω R = 132.4 cm -1 . The temperature is set to be T = 0 K. For the L-MFE simulation, the corresponding initial condition is ρS (0