The star formation activity of IllustrisTNG galaxies: main sequence, UVJ diagram, quenched fractions, and systematics

We select galaxies from the IllustrisTNG hydrodynamical simulations ($M_{\rm stars}\gt 10^9 \, {\rm M}_\odot$ at 0 ≤ z ≤ 2) and characterize the shapes and evolutions of their UVJ and star formation rate–stellar mass (SFR–Mstars) diagrams. We quantify the systematic uncertainties related to different criteria to classify star-forming versus quiescent galaxies, different SFR estimates, and by accounting for the star formation measured within different physical apertures. The TNG model returns the observed features of the UVJ diagram at z ≤ 2, with a clear separation between two classes of galaxies. It also returns a tight star-forming main sequence (MS) for $M_{\rm stars}\lt 10^{10.5} \, ({\rm M}_\odot)$ with a ∼0.3 dex scatter at z ∼ 0 in our fiducial choices. If a UVJ-based cut is adopted, the TNG MS exhibits a downwardly bending at stellar masses of about 1010.5−10.7 M⊙. Moreover, the model predicts that ${\sim }80\, (50)$ per cent of 1010.5−11 M⊙ galaxies at z = 0 (z = 2) are quiescent and the numbers of quenched galaxies at intermediate redshifts and high masses are in better agreement with observational estimates than previous models. However, shorter SFR-averaging time-scales imply higher normalizations and scatter of the MS, while smaller apertures lead to underestimating the galaxy SFRs: overall we estimate the inspected systematic uncertainties to sum up to about 0.2−0.3 dex in the locus of the MS and to about 15 percentage points in the fraction of quenched galaxies. While TNG colour distributions are clearly bimodal, this is not the case for the SFR logarithmic distributions in bins of stellar mass (SFR ≳ 10−3 M⊙yr−1). Finally, the slope and z = 0 normalization of the TNG MS are consistent with observational findings; however, the locus of the TNG MS remains lower by about 0.2−0.5 dex at 0.75 ≤ z < 2 than the available observational estimates taken at face value.