where ω max ≫ ω b is the maximum frequency when discretizing the bath frequencies. To sample within the classical regime, we choose ω max = 5 ω b . The bath parameters are ω b = 18 cm -1 (which is 2.2 meV) and λ b = 50 cm -1 (which is 6.2 meV). The characteristic phonon frequency ω b leads to a bath correlation function decay timescale of 300 fs. The initial bath condition ρb (see Eq. 55) is described by its Wigner transform with the following analytic expression,