Abstract:[Abridged] In recent years, conflicting results have provided an uncertain view of the dust-attenuated properties of $z>4$ star-forming galaxies (SFGs). To solve this, we used the deepest data publicly available in COSMOS to build a mass-complete ($>10^{9.5}\,M_{\odot}$) sample of SFGs at $4<z<5$ and measured their dust-attenuated properties by stacking all archival ALMA band 6 and 7 observations available. Combining this information with their rest-frame ultraviolet emission from the COSMOS2020 catalog, we constrained the IRX ($\equiv L_{\rm IR}/L_{\rm UV}$)--$\beta_{\rm UV}$, IRX--$M_\ast$, and SFR--$M_\ast$ relations at $z\sim4.5$. Finally, using these relations and the stellar mass function of SFGs at $z\sim4.5$, we inferred the unattenuated and dust-attenuated SFRD at this epoch. SFGs at $z\sim4.5$ follow an IRX--$\beta_{\rm UV}$ relation that is consistent with that of local starbursts, while they follow a steeper IRX--$M_\ast$ relation than observed locally. The grain properties of dust in these SFGs seems thus similar to those in local starbursts but its mass and geometry result in lower attenuation in low-mass SFGs. SFGs at $z\sim4.5$ lie on a linear SFR--$M_\ast$ relation, whose normalization varies by 0.3 dex, when we exclude or include from our stacks the ALMA primary targets. The cosmic SFRD$(>M_\ast)$ converges at $M_\ast<10^{9}\,M_\odot$ and is dominated by SFGs with $M_\ast\sim10^{9.5-10.5}\,M_\odot$. The fraction of the cosmic SFRD that is attenuated by dust, ${\rm SFRD}_{\rm IR}(>M_\ast)/ {\rm SFRD}(>M_\ast)$, is $90\pm4\%$ for $M_\ast\,=\,10^{10}\,M_\odot$, $68\pm10\%$ for $M_\ast=10^{8.9}\,M_\odot$ (i.e., $0.03\times M^\star$; $M^\star$ being the characteristic stellar mass of SFGs) and this value converges to $60\pm10\%$ for $M_\ast=10^{8}\,M_\odot$. Even at this early epoch, the fraction of the cosmic SFRD that is attenuated by dust remains thus significant.

Dust Attenuation in Uv-Selected Starbursts at High Redshift and Their Local Counterparts: Implications for the Cosmic Star Formation Rate Density

Dust attenuation, dust content and geometry of star-forming galaxies

The VANDELS survey: Dust attenuation in star-forming galaxies at $\mathbf{z=3-4}$

The HDUV Survey: A Revised Assessment of the Relationship between UV Slope and Dust Attenuation for High-redshift Galaxies

Dust in High-Redshift Galaxies: Reconciling UV Attenuation and IR Emission

Dust attenuation in the restframe ultraviolet: constraints from star-forming galaxies at z~1

The Ultraviolet Spectroscopic Properties of Local Starbursts: Implications at High-Redshift

A^3COSMOS: Measuring the cosmic dust-attenuated star formation rate density at 4 < z < 5

A$^3$COSMOS: Measuring the cosmic dust-attenuated star formation rate density at $4 < z < 5$

Interpreting the observed UV continuum slopes of high-redshift galaxies

The blue UV slopes of z~4 Lyman break galaxies: implications for the corrected star formation rate density

The slippery slope of dust attenuation curves: Correlation of dust attenuation laws with star-to-dust compactness up to z = 4

Star Formation from DLA Gas in the Outskirts of Lyman Break Galaxies at z~3

Mapping dust attenuation and the 2175 Å bump at kpc scales in nearby galaxies

On the Dust properties of the UV galaxies in the redshift range $z \sim 0.6-1.2$

Effects of dust sources on dust attenuation properties in IllustrisTNG galaxies at $z\sim 7$

The Ultraviolet Slopes of Early Universe Galaxies: The Impact of Bursty Star Formation, Dust, and Nebular Continuum Emission

THE SPACE DENSITY OF LUMINOUS DUSTY STAR-FORMING GALAXIES AT Z > 4: SCUBA-2 AND LABOCA IMAGING OF ULTRARED GALAXIES FROM HERSCHEL-ATLAS

Universal Selftrapping in Nonlinear Tight-binding Lattices

Herschel PEP/HerMES: the redshift evolution (0 ≤ z ≤ 4) of dust attenuation and of the total (UV+IR) star formation rate density

On the dust properties of the UV galaxies in the redshift range z ∼ 0.6–1.2