SERENADE II: An ALMA Multi-Band Dust-Continuum Analysis of 28 Galaxies at $5<z<8$ and the Physical Origin of the Dust Temperature Evolution

We present an analysis of ALMA multi-band dust-continuum observations for 28 spectroscopically-confirmed bright Lyman-break galaxies at $5<z<8$. Our sample consists of 11 galaxies at $z\sim6$ newly observed in our ALMA program, which substantially increases the number of $5<z<8$ galaxies with both rest-frame 88 and 158 $\mu{\rm m}$ continuum observations, allowing us to simultaneously measure the IR luminosity and dust temperature for a statistical sample of $z\gtrsim5$ galaxies for the first time. We derive the relationship between the UV slope ($\beta_{\rm UV}$) and infrared excess (IRX) for the $z\sim6$ galaxies, and find a shallower IRX-$\beta_{\rm UV}$ relation compared to the previous results at $z\sim2$--4. Based on the IRX-$\beta_{\rm UV}$ relation consistent with our results and the $\beta_{\rm UV}$-$M_{\rm UV}$ relation including fainter galaxies in the literature, we find a limited contribution of the dust-obscured star formation to the total SFR density, $\sim30\%$ at $z\sim6$. Our measurements of the dust temperature at $z\sim6-7$, $T_{\rm dust}=40.9_{-9.1}^{+10.0}\,{\rm K}$ on average, supports a gentle increase of $T_{\rm dust}$ from $z=0$ to $z\sim6$--7. Using an analytic model with parameters consistent with recent {\it{JWST}} results, we discuss that the observed redshift evolution of the dust temperature can be reproduced by an $\sim0.6\,{\rm dex}$ increase in the gas depletion timescale and $\sim0.4\,{\rm dex}$ decrease of the metallicity. The variety of $T_{\rm dust}$ observed at high redshifts can also be naturally explained by scatters around the star-formation main sequence and average mass-metallicity relation, including an extremely high dust temperature of $T_{\rm dust}>80\,{\rm K}$ observed in a galaxy at $z=8.3$.