G2) where m ν is the nuclear masse set to be 1, ∆ is the diabatic electronic coupling set to be 1, and ε is the energy bias between two electronic states. The F nuclear modes each have their frequency ω ν and vibronic coupling coefficient c ν , determined through the spectral density of the bath J(ω). Here, we use an Ohmic spectral density