Going deeper into the dark with COSMOS-Web: JWST unveils the total contribution of Radio-Selected NIRfaint galaxies to the cosmic Star Formation Rate Density

We present the first follow-up with JWST of radio-selected NIRfaint galaxies as part of the COSMOS-Web survey. By selecting galaxies detected at radio frequencies ($S_{\rm 3 GHz}>11.5$ $\mu$Jy; i.e. S/N$>5$) and with faint counterparts at NIR wavelengths (F150W$>26.1$ mag), we collect a sample of 127 likely dusty star-forming galaxies (DSFGs). We estimate their physical properties through SED fitting, compute the first radio luminosity function for these types of sources, and their contribution to the total cosmic star formation rate density. Our analysis confirms that these sources represent a population of highly dust-obscured ($\langle A_{\rm v} \rangle \sim3.5$ mag), massive ($\langle M_\star \rangle \sim10^{10.8}$ M$_\odot$) and star-forming galaxies ($\langle {\rm SFR} \rangle\sim300$ M$_\odot$ yr$^{-1}$) located at $\langle z \rangle\sim3.6$, representing the high-redshift tail of the full distribution of radio sources. Our results also indicate that these galaxies could dominate the bright end of the radio luminosity function and reach a total contribution to the cosmic star formation rate density equal to that estimated only considering NIR-bright sources at $z\sim4.5$. Finally, our analysis further confirms that the radio selection can be employed to collect statistically significant samples of DSFGs, representing a complementary alternative to the other selections based on JWST colors or detection at FIR/(sub)mm wavelengths.