SEMPER I. Radio Predictions for Star-Forming Galaxies at $0<z<5$

Abstract: [Abridged] SFGs are the dominant population in the faint radio sky, corresponding to flux densities at 1.4 GHz $< 0.1$ mJy. A panchromatic approach is essential for selecting SFGs in the radio band and understanding star formation processes over cosmic time. Semi-empirical models are valuable tools to effectively study galaxy formation and evolution, relying on minimal assumptions and exploiting empirical relations between galaxy properties and enabling us to take full advantage of the recent progress in radio and optical/near-infrared (NIR) observations. In this paper, we develop the Semi-EMPirical model for Extragalactic Radio emission (SEMPER) to predict radio luminosity functions and number counts at 1.4 GHz and 150 MHz for SFGs. SEMPER is based on state-of-the-art empirical relations and combines the redshift-dependent galaxy stellar mass functions obtained from the recent COSMOS2020 catalogue, which exploits deep near-infrared observations, with up-to-date observed scaling relations, such as the galaxy main sequence and the mass-dependent far-infrared/radio correlation across cosmic time. Our luminosity functions are compared with recent observational determinations from several radio telescopes, along with previous semi-empirical models and simulations. Our semi-empirical model successfully reproduces the observed luminosity functions at 1.4 GHz and 150 MHz up to $z\sim 5$ and the most recent number count statistics from radio observations in the LoTSS deep fields. Our model, based on galaxies selected in the NIR, naturally predicts the presence of radio-selected massive and/or dust-obscured galaxies already in place at high redshift ($z\gtrsim3.5$), as suggested by recent results from JWST. Our predictions offer an excellent benchmark for upcoming updates from JWST and future ultra-deep radio surveys planned with the SKA and its precursors.