Broadband spectral energy distributions of SDSS-selected quasars and of their host galaxies: intense activity at the edge of the quenching

We present a multiwavelength reconstruction, from the near-ultraviolet (near-UV) to the far-infrared (far-IR), of the mean spectral energy distribution (SED) of optically selected quasars (QSOs) from the Sloan Digital Sky Survey (SDSS) over the redshift range $1 \le z \le 5$. We achieve this goal by combining photometric information from SDSS, UKIDSS, and WISE surveys as well as complementing with a stacking analysis of Herschel, AKARI, and Planck maps at the location of the QSOs. The near-UV and optical parts of the reconstructed mean rest-frame SED are similar to those found in other studies. However, the SED shows an excess at 1-2 $\mu$m (when compared to the aforementioned SEDs normalized in the near-UV) and a prominent bump around 4-6 $\mu$m, followed by a decrease out to $\sim 20 \,\mu$m and a subsequent far-IR bump. We perform SED fitting to the composite SEDs in order to estimate the average active galactic nuclei (AGN) luminosity $L_{\rm AGN}$ and star formation rate (SFR) as function of cosmic time, finding typical $L_{\rm AGN} \sim 10^{46} - 10^{47}$ erg/s and SFR $\sim 50 - 1000\, M_{\odot}/$yr. We develop mid-IR based criteria to split the QSO sample, finding that these allow us to move along the average relationship in the $L_{\rm AGN}-$SFR plane toward increasing luminosities. We also find the detection in the far-IR band to be an effective criterion to select objects where the star formation is on the verge of being quenched by AGN activity.