Since Ĥa,jk is the only term that contributes to the population transfer between dark states, such transition is equivalent to a population transfer between unshifted harmonic oscillators, only involving contributions from off-diagonal system-bath coupling. The correlation function is given by (for details of the derivation, see Sec. IIA of the supplementary material), and the final expression is