where ρSS is the steady-state reduced density matrix (RDM) obtained using the HEOM steady-state solver , which has significant advantages over the Krylov subspace solver or long-time propagation using RK-4 integrator. Details on the procedures of the HEOM steady-state solver can be found in Ref. 62. For all the numerical simulations, the dipole operator is expressed as