An infrared-radio simulation of the extragalactic sky: from the Square Kilometre Array to Herschel

To exploit synergies between the Herschel Space Observatory and next generation radio facilities, we have extended the semi-empirical extragalactic radio continuum simulation of Wilman et al. to the mid- and far-infrared. Here, we describe the assignment of infrared spectral energy distributions (SEDs) to the star-forming galaxies and active galactic nuclei, using Spitzer 24, 70 and 160 μm and SCUBA 850 μm survey results as the main constraints. Star-forming galaxies dominate the source counts, and a model in which their far-infrared–radio correlation and infrared SED assignment procedure are invariant with redshift underpredicts the observed 24 and 70 μm source counts. The 70 μm deficit can be eliminated if the star-forming galaxies undergo stronger luminosity evolution than originally assumed for the radio simulation, a requirement which may be partially ascribed to known non-linearity in the far-infrared–radio correlation at low luminosity if it evolves with redshift. At 24 μm, the shortfall is reduced if the star-forming galaxies develop SEDs with cooler dust and correspondingly stronger polycyclic aromatic hydrocarbon emission features with increasing redshift at a given far-infrared luminosity, but this trend may reverse at z > 1 in order not to overproduce the submillimetre source counts. The resulting model compares favourably with recent Balloon-borne Large Aperture Submillimetre Telescope (BLAST) results, and we have extended the simulation data base to aid the interpretation of Herschel surveys. Such comparisons may also facilitate further model refinement and revised predictions for the Square Kilometre Array and its precursors