Here β = 1 k B T is the inverse temperature and J 0 (ω) is the solvent spectral density. The quantity J 0 (ω) is well-described by a Drude-Lorentz functional form (see main text) which, in turn, is fully determined by the solvent correlation time 1 Λ and reorganization energy λ 0 . The Raman cross section has an additional vibrational lineshape d(τ ) which is usually described through a normalized Lorentzian.