Investigating the influence of radio-faint active galactic nuclei on the infrared-radio correlation of massive galaxies

Context. It is well known that star-forming galaxies (SFGs) exhibit a tight correlation between their radio and infrared emissions, commonly referred to as the infrared-radio correlation (IRRC). Recent empirical studies have reported a dependence of the IRRC on the galaxy stellar mass, in which more massive galaxies tend to show lower infrared-to-radio ratios (qIR) with respect to less massive galaxies. One possible, yet unexplored, explanation is a residual contamination of the radio emission from active galactic nuclei (AGNs), not captured through "radio-excess" diagnostics. Aims. To investigate this hypothesis, we aim to statistically quantify the contribution of AGN emission to the radio luminosities of SFGs located within the scatter of the IRRC. Methods. Our Very Large Baseline Array (VLBA) AGN-sCAN program has targeted 500 galaxies that follow the qIR distribution of the IRRC, i.e., with no prior evidence for radio-excess AGN emission based on low-resolution (similar to arcsec) VLA radio imaging. Our VLBA 1.4 GHz observations reach a 5 sigma sensitivity limit of 25 mu Jy/beam, corresponding to a radio-brightness temperature of T-b similar to 10(5) K. This classification serves as a robust AGN diagnostic, regardless of the host galaxy's star formation rate. Results. We detect four VLBA sources in the deepest regions, which are also the faintest VLBI-detected AGNs in SFGs to date. The effective AGN detection rate is 9%, when considering a control sample matched in mass and sensitivity, which is in good agreement with the extrapolation of previous radio AGN number counts. Despite the non-negligible AGN flux contamination (similar to 30%) in our individual VLBA detections, we find that the peak of the qIR distribution is completely unaffected by this correction. Although we cannot rule out a high incidence of radio-silent AGNs at (sub)mu Jy levels among the VLBA non-detections, we derive a conservative upper limit of < 0.1 dex of their cumulative impact on the qIR distribution. We conclude that residual AGN contamination from non-radio-excess AGNs is unlikely to be the primary driver of the M-star - dependent IRRC.