2019MNRAS.489.4196L

We combine the SHARK semi-analytic model of galaxy formation with the PROSPECT software tool for spectral energy distribution (SED) generation to study the multiwavelength emission of galaxies from the far-ultraviolet (FUV) to the far-infrared (FIR) at 0 <= z <= 10. We produce a physical model for the attenuation of galaxies across cosmic time by combining a local Universe empirical relation to compute the dust mass of galaxies from their gas metallicity and mass, attenuation curves derived from radiative transfer calculations of galaxies in the EAGLE hydrodynamic simulation suite, and the properties of SHARK galaxies. We are able to produce a wide range of galaxies, from the z = 8 star-forming galaxies with almost no extinction, z = 2 submillimetre galaxies, down to the normal star-forming and red-sequence galaxies at z = 0. Quantitatively, we find that SHARK reproduces the observed (i) z = 0 FUV-to-FIR, (ii) 0 <= z <= 3 rest-frame K-band, and (iii) 0 <= z <= 10 rest-frame FUV luminosity functions, (iv) z <= 8 UV slopes, (v) the FUV-to-FIR number counts (including the widely disputed 850 µm), (vi) redshift distribution of bright 850 µm galaxies, and (vii) the integrated cosmic SED from z = 0 to 1 to an unprecedented level. This is achieved without the need to invoke changes in the stellar initial mass function, dust-to-metal mass ratio, or metal enrichment time-scales. Our model predicts star formation in galaxy discs to dominate in the FUV-to-optical, while bulges dominate at the NIR at all redshifts. The FIR sees a strong evolution in which discs dominate at z <= 1 and starbursts (triggered by both galaxy mergers and disc instabilities, in an even mix) dominate at higher redshifts, even out to z = 10.