We interpret recent ALMA observations of z > 6 normal star forming galaxies by means of a semi-numerical method, which couples the output of a cosmological hydrodynamical simulation with a chemical evolution model which accounts for the contribution to dust enrichment from supernovae, asymptotic giant branch stars and grain growth in the interstellar medium. We find that while stellar sources dominate the dust mass of small galaxies, the higher level of metal enrichment experienced by galaxies with M star > 10 9 M ⊙ allows efficient grain growth, which provides the dominant contribution to the dust mass. Even assuming maximally efficient supernova dust production, the observed dust mass of the z = 7.5 galaxy A1689-zD1 requires very efficient grain growth. This, in turn, implies that in this galaxy the average density of the cold and dense gas, where grain growth occurs, is comparable to that inferred from observations of QSO host galaxies at similar redshifts. Although plausible, the upper limits on the dust continuum emission of galaxies at 6.5 < z < 7.5 show that these conditions must not apply to the bulk of the high redshift galaxy population.

We attempt to interpret existing data on the evolution of the UV luminosity function and UV colours, β, of galaxies at 5 ≤ z ≤ 8, to improve our understanding of their dust content and ISM properties. To this aim, we post-process the results of a cosmological hydrodynamical simulation with a chemical evolution model, which includes dust formation by supernovae and intermediate mass stars, dust destruction in supernova shocks, and grain growth by accretion of gas-phase elements in dense gas. We find that observations require a steep, Small Magellanic Cloud-like extinction curve and a clumpy dust distribution, where stellar populations younger than 15 Myr are still embedded in their dusty natal clouds. Investigating the scatter in the colour distribution and stellar mass, we find that the observed trends can be explained by the presence of two populations: younger, less massive galaxies where dust enrichment is mainly due to stellar sources, and massive, more chemically evolved ones, where efficient grain growth provides the dominant contribution to the total dust mass. Computing the IR-excess -UV color relation we find that all but the dustiest model galaxies follow a relation shallower than the Meurer et al. (1999) one, usually adopted to correct the observed UV luminosities of high-z galaxies for the effects of dust extinction. As a result, their total star formation rates might have been over-estimated. Our study illustrates the importance to incorporate a proper treatment of dust in simulations of high-z galaxies, and that massive, dusty, UV-faint galaxies might have already appeared at z 7.

Context. Apertif is a phased-array feed system for the Westerbork Synthesis Radio Telescope, providing forty instantaneous beams over 300 MHz of bandwidth. A dedicated survey program utilizing this upgrade started on 1 July 2019, with the last observations taken on 28 February 2022. The imaging survey component provides radio continuum, polarization, and spectral line data. Aims. Public release of data is critical for maximizing the legacy of a survey. Toward that end, we describe the release of data products from the first year of survey operations, through 30 June 2020. In particular, we focus on defining quality control metrics for the processed data products. Methods. The Apertif imaging pipeline, Apercal, automatically produces non-primary beam corrected continuum images, polarization images and cubes, and uncleaned spectral line and dirty beam cubes for each beam of an Apertif imaging observation. For this release, processed data products are considered on a beam-by-beam basis within an observation. We validate the continuum images by using metrics that identify deviations from Gaussian noise in the residual images. If the continuum image passes validation, we release all processed data products for a given beam. We apply further validation to the polarization and line data products and provide flags indicating the quality of those data products. Results. We release all raw observational data from the first year of survey observations, for a total of 221 observations of 160 independent target fields, covering approximately one thousand square degrees of sky. Images and cubes are released on a per beam basis, and 3374 beams (of 7640 considered) are released. The median noise in the continuum images is 41.4 uJy beam −1 , with a slightly lower median noise of 36.9 uJy beam −1 in the Stokes V polarization image. The median angular resolution is 11.6 /sin δ. The median noise for all line cubes, with a spectral resolution of 36.6 kHz, is 1.6 mJy beam −1 , corresponding to a 3-σ H i column density sensitivity of 1.8×10 20 atoms cm −2 over 20 km s −1 (for a median angular resolution of 24 ×15 ). Line cubes at lower frequency have slightly higher noise values, consistent with the global RFI environment and overall Apertif system performance. We also provide primary beam images for each individual Apertif compound beam. The data are made accessible using a Virtual Observatory interface and can be queried using a variety of standard tools.