Not yet reviewed by Pith; the record is open.

We perform a stacking analysis of {\it Planck}, {\it AKARI}, Infrared Astronomical Satellite ($IRAS$), Wide-field Infrared Survey Eplorer ($WISE$), and {\it Herschel} images of the largest number of (candidate) protoclusters at $z\sim3.8$ selected from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). Stacking the images of the $179$ candidate protoclusters, the combined infrared (IR) emission of the protocluster galaxies in the observed $12-850~\mu$m wavelength range is successfully detected with $>5\sigma$ significance (at $Planck$). This is the first time that the average IR spectral energy distribution (SED) of a protocluster has been constrained at $z\sim4$. The observed IR SEDs of the protoclusters exhibit significant excess emission in the mid-IR compared to that expected from typical star-forming galaxies (SFGs). They are reproduced well using SED models of intense starburst galaxies with warm/hot dust heated by young stars, or by a population of active galactic nuclei (AGN)/SFG composites. For the pure star-forming model, a total IR (from 8 to 1000 $\mu$m) luminosity of $19.3_{-4.2}^{+0.6}\times10^{13}~L_{\odot}$ and a star formation rate (SFR) of $16.3_{-7.8}^{+1.0}\times10^3~M_{\odot}$ yr$^{-1}$ are found whereas for the AGN/SFG composite model, $5.1_{-2.5}^{+2.5}\times10^{13}~L_{\odot}$ and $2.1^{+6.3}_{-1.7}\times10^3~M_{\odot}$ yr$^{-1}$ are found. Uncertainty remaining in the total SFRs; however, the IR luminosities of the most massive protoclusters are likely to continue increasing up to $z\sim4$. Meanwhile, no significant IR flux excess is observed around optically selected QSOs at similar redshifts, which confirms previous results. Our results suggest that the $z\sim4$ protoclusters trace dense, intensely star-forming environments that may also host obscured AGNs missed by the selection in the optical.