Additional information: Use policyThe full-text may be used and/or reproduced, and given to third parties in any format or medium, without prior permission or charge, for personal research or study, educational, or not-for-prot purposes provided that:• a full bibliographic reference is made to the original source • a link is made to the metadata record in DRO • the full-text is not changed in any way The full-text must not be sold in any format or medium without the formal permission of the copyright holders.Please consult the full DRO policy for further details. ABSTRACTWe present a new version of the GALFORM semi-analytical model of galaxy formation. This brings together several previous developments of GALFORM into a single unified model, including a different initial mass function (IMF) in quiescent star formation and in starbursts, feedback from active galactic nuclei supressing gas cooling in massive halos, and a new empirical star formation law in galaxy disks based on their molecular gas content. In addition, we have updated the cosmology, introduced a more accurate treatment of dynamical friction acting on satellite galaxies, and updated the stellar population model. The new model is able to simultaneously explain both the observed evolution of the K-band luminosity function and stellar mass function, and the number counts and redshift distribution of sub-mm galaxies selected at 850µm. This was not previously achieved by a single physical model within the ΛCDM framework, but requires having an IMF in starbursts that is somewhat top-heavy. The new model is tested against a wide variety of observational data covering wavelengths from the far-UV to sub-mm, and redshifts from z = 0 to z = 6, and is found to be generally successful. These observations include the optical and near-IR luminosity functions, HI mass function, fraction of early type galaxies, Tully-Fisher, metallicity-luminosity and size-luminosity relations at z = 0, as well as far-IR number counts, and far-UV luminosity functions at z ∼ 3 − 6. Discrepancies are however found in galaxy sizes and metallicities at low luminosities, and in the abundance of low mass galaxies at high-z, suggesting the need for a more sophisticated model of supernova feedback.

We report the detection of extended Lyα haloes around 145 individual star-forming galaxies at redshifts 3 ≤ z ≤ 6 in the Hubble Ultra Deep Field observed with the Multi-Unit Spectroscopic Explorer (MUSE) at ESO-VLT. Our sample consists of continuum-faint (−15 ≥ M UV ≥ −22) Lyα emitters (LAEs). Using a 2D, two-component (continuum-like and halo) decomposition of Lyα emission assuming circular exponential distributions, we measure scale lengths and luminosities of Lyα haloes. We find that 80% of our objects having reliable Lyα halo measurements show Lyα emission that is significantly more extended than the UV continuum detected by HST (by a factor ≈4 to >20). The median exponential scale length of the Lyα haloes in our sample is ≈4.5 kpc with a few haloes exceeding 10 kpc. By comparing the maximal detected extent of the Lyα emission with the predicted dark matter halo virial radii of simulated galaxies, we show that the detected Lyα emission of our selected sample of Lyα emitters probes a significant portion of the cold circum-galactic medium of these galaxies (>50% in average). This result therefore shows that there must be significant HI reservoirs in the circum-galactic medium and reinforces the idea that Lyα haloes are ubiquitous around high-redshift Lyα emitting galaxies. Our characterization of the Lyα haloes indicates that the majority of the Lyα flux comes from the halo (≈65%) and that their scale lengths seem to be linked to the UV properties of the galaxies (sizes and magnitudes). We do not observe a significant Lyα halo size evolution with redshift, although our sample for z > 5 is very small. We also explore the diversity of the Lyα line profiles in our sample and we find that the Lyα lines cover a large range of full width at half maximum (FWHM) from 118 to 512 km s −1 . While the FWHM does not seem to be correlated to the Lyα scale length, most compact Lyα haloes and those that are not detected with high significance tend to have narrower Lyα profiles (<350 km s −1 ). Finally, we investigate the origin of the extended Lyα emission but we conclude that our data do not allow us to disentangle the possible processes, i.e. scattering from star-forming regions, fluorescence, cooling radiation from cold gas accretion, and emission from satellite galaxies.

Detailed studies of the spectral energy distribution (SED) of normal galaxies have increasingly been used to understand the physical mechanism dominating their integrated emission, mainly owing to the availability of high quality multi-wavelength data from the UV to the far-infrared (FIR). However, systems hosting dust-enshrouded nuclear starbursts and/or an accreting supermassive black hole (an active galactic nucleus or AGN) are especially challenging to study. This is due to the complex interplay between the heating by massive stars and the AGN, the absorption and emission of radiation from dust, as well as the presence of the underlying old stellar population. We used the latest release of CIGALE, a fast state-of-the-art galaxy SED-fitting model relying on energy balance, to study the influence of an AGN in a self consistent manner in estimating both the star formation rate (SFR) and stellar mass in galaxies, as well as to calculate the contribution of the AGN to the power output of the host. Using the semi-analytical galaxy formation model , we created a suite of mock galaxy SEDs using realistic star formation histories (SFH). We also added an AGN of Type-1, Type-2, or intermediate-type whose contribution to the bolometric luminosity can be variable. We performed an SED-fitting of these catalogues with CIGALE, assuming three different SFHs: a single-exponentially-decreasing (1τ-dec), a double-exponentiallydecreasing (2τ-dec), and a delayed SFH. Constraining the overall contribution of an AGN to the total infrared luminosity (frac AGN ) is very challenging for frac AGN < 20%, with uncertainties of ∼5-30% for higher fractions depending on the AGN type, while FIR and sub-mm are essential. The AGN power has an impact on the estimation of M * in Type-1 and intermediate-type AGNs but has no effect on galaxies hosting Type-2 AGNs. We find that in the absence of AGN emission, the best estimates of M * are obtained using the 2τ-dec model but at the expense of realistic ages of the stellar population. The delayed SFH model provides good estimates of M * and SFR, with a maximum offset of 10% as well as better estimates of the age. Our analysis shows that the under-estimation of the SFR increases with frac AGN for Type-1 systems, as well as for low contributions of an intermediate AGN type, but it is quite insensitive to the emission of Type-2 AGNs up to frac AGN ∼ 45%. A lack of sampling the FIR, or sub-mm domain systematically over-estimates the SFR (<20%), independent of the contribution of the AGN. Similarly, the UV emission is critical in accurately retrieving both the M * for Type-1 and intermediate-type AGN and the SFR of all three AGN types. We show that the presence of AGN emission introduces a scatter to the SFR-M * main sequence relation derived from SED-fitting, which is driven by the uncertainties on M * . Finally, we used our mock catalogues to test the popular IR SED-fitting code DIR and show that frac AGN is under-estimated but that the SFR is recovered well for Type-1 and intermedia...