Source classifications, stellar masses and star formation rates are presented for 80,000 radio sources from the first data release of the Low Frequency Array Two-metre Sky Survey (LoTSS) Deep Fields, which represents the widest deep radio survey ever undertaken. Using deep multi-wavelength data spanning from the ultraviolet to the far-infrared, spectral energy distribution (SED) fitting is carried out for all of the LoTSS-Deep host galaxies using four different SED codes, two of which include modelling of the contributions from an active galactic nucleus (AGN). Comparing the results of the four codes, galaxies that host a radiative AGN are identified, and an optimised consensus estimate of the stellar mass and star-formation rate for each galaxy is derived. Those galaxies with an excess of radio emission over that expected from star formation are then identified, and the LoTSS-Deep sources are divided into four classes: star-forming galaxies, radio-quiet AGN, and radio-loud high-excitation and low-excitation AGN. Ninety-five per cent of the sources can be reliably classified, of which more than two-thirds are star-forming galaxies, ranging from normal galaxies in the nearby Universe to highly-starbursting systems at z>4. Star-forming galaxies become the dominant population below 150-MHz flux densities of about 1 mJy, accounting for 90 per cent of sources at a 150-MHz flux density of 100 microJy. Radio-quiet AGN comprise around 10 per cent of the overall population. Results are compared against the predictions of the SKADS and T-RECS radio sky simulations, and improvements to the simulations are suggested.

The far-infrared radio correlation (FIRC) is one of the strongest correlations in astronomy, yet a model that explains this comprehensively does not exist. The new LOFAR all Sky Survey (LoTSS) deep field, ELAIS-N1, allows exploration of this relation in previously unexplored regions of parameter space of radio frequency (150 MHz), luminosity (L150 < 1024.7), redshift (z ∼ 1), and stellar mass M* < 1011.4. We present accurate deblended far-infrared (FIR) flux measurements with robust errors at 24, 100, 160, 250, 350, and 500 μm from Spitzer and the Herschel Space Observatory using XID+. We find that the FIRC has a strong mass dependence, the evolution of which takes the form qTIR(M*) = (2.00 ± 0.01)+(−0.22 ± 0.02)(log(M/M*)−10.05). This matches recent findings in regards to the star formation rate–radio luminosity relation at 150 MHz and results from radio observations in COSMOS at 1.4 GHz with the Jansky Very Large Array (JVLA). Our results provide tighter constraints on the low-redshift end of the FIRC and at lower frequency than the COSMOS observations. In addition, we find a mild evolution with redshift, with a best fit relation qTIR(z) = (1.94 ± 0.01)(1 + z)−0.04 ± 0.01. This evolution is shallower than that suggested by previous results at 150 MHz with the differences explained by the fact that previous studies did not account for the mass dependence. Finally, we present deblended FIR fluxes for 79 609 galaxies across the LoTSS deep fields: Boötes, ELAIS-N1, and Lockman Hole.

Aims. We provide the most accurate estimate yet of the bright end of the infrared (IR) luminosity functions (LFs) and the abundance of hyperluminous IR galaxies (HLIRGs) with IR luminosities >1013L⊙, thanks to the combination of the high sensitivity, angular resolution, and large area of the LOFAR Deep Fields, which probes an unprecedented dynamic range of luminosity and volume.Methods. We cross-match Herschel sources and LOFAR sources in Boötes (8.63 deg2), Lockman Hole (10.28 deg2), and ELAIS-N1 (6.74 deg2) with rms sensitivities of ~32, 22, and 20 μJy beam−1, respectively. We divide the matched samples into “unique” and “multiple” categories. For the multiple matches, we de-blend the Herschel fluxes using the LOFAR positions and the 150-MHz flux densities as priors. We perform spectral energy distribution fitting, combined with multi-wavelength counterpart identifications and photometric redshift estimates, to derive IR luminosities.Results. The depth of the LOFAR data allows us to identify highly complete (~92% completeness) samples of bright Herschel sources with a simple selection based on the 250 μm flux density (45, 40, and 35 mJy in Boötes, Lockman Hole, and ELAIS-N1, respectively). Most of the bright Herschel sources fall into the unique category (i.e. a single LOFAR counterpart). For the multiple matches, there is excellent correspondence between the radio emission and the far-IR emission. We find a good agreement in the IR LFs with a previous study out to z ~ 6 which used de-blended Herschel data. Our sample gives the strongest and cleanest indication to date that the population of HLIRGs has surface densities of ~5 to ~18/deg2 (with variations due to a combination of the applied flux limit and cosmic variance) and an uncertainty of a factor of ≲2. In comparison, the GALFORM semi-analytic model significantly under-predicts the abundance of HLIRGs.

Aims. We investigate the radio properties of a sample of 850 μm-selected sources from the SCUBA-2 Cosmology Legacy Survey (S2CLS) using new deep, low-frequency radio imaging of the Lockman Hole field from the Low Frequency Array. This sample consists of 53 sources, 41 of which are detected at >5σ at 150 MHz.Methods. Combining these data with additional observations at 324 MHz, 610 MHz, and 1.4 GHz from the Giant Metrewave Radio Telescope and the Jansky Very Large Array, we find a variety of radio spectral shapes and luminosities (L1.4 GHz ranging from ~4 × 1023−1 × 1025) within our sample despite their similarly bright submillimetre flux densities (>4 mJy). We characterise their spectral shapes in terms of multi-band radio spectral indices. Finding strong spectral flattening at low frequencies in ~20% of sources, we investigate the differences between sources with extremely flat low-frequency spectra and those with ‘normal’ radio spectral indices (α > −0.25).Results. As there are no other statistically significant differences between the two subgroups of our sample as split by the radio spectral index, we suggest that any differences are undetectable in galaxy-averaged properties that we can observe with our unresolved images, and likely relate to galaxy properties that we cannot resolve, on scales ≲1 kpc. We attribute the observed spectral flattening in the radio to free–free absorption, proposing that those sources with significant low-frequency spectral flattening have a clumpy distribution of star-forming gas. We estimate an average spatial extent of absorbing material of at most several hundred parsecs to produce the levels of absorption observed in the radio spectra. This estimate is consistent with the highest-resolution observations of submillimetre galaxies in the literature, which find examples of non-uniform dust distributions on scales of ~100 pc, with evidence for clumps and knots in the interstellar medium. Additionally, we find two bright (>6 mJy) S2CLS sources undetected at all other wavelengths. We speculate that these objects may be very high redshift sources, likely residing at z > 4.

We present the Herschel Extragalactic Legacy Project (HELP). This project collates, curates, homogenizes, and creates derived data products for most of the premium multiwavelength extragalactic data sets. The sky boundaries for the first data release cover 1270 deg(2) defined by the Herschel SPIRE extragalactic survey fields; notably the Herschel Multi-tiered Extragalactic Survey (HerMES) and the Herschel Atlas survey (H-ATLAS). Here, we describe the motivation and principal elements in the design of the project. Guiding principles are transparent or 'open' methodologies with care for reproducibility and identification of provenance. A key element of the design focuses around the homogenization of calibration, meta data, and the provision of information required to define the selection of the data for statistical analysis. We apply probabilistic methods that extract information directly from the images at long wavelengths, exploiting the prior information available at shorter wavelengths and providing full posterior distributions rather than maximum-likelihood estimates and associated uncertainties as in traditional catalogues. With this project definition paper, we provide full access to the first data release of HELP; Data Release 1 (DR1), including a monolithic map of the largest SPIRE extragalactic field at 385 deg(2) and 18 million measurements of PACS and SPIRE fluxes. We also provide tools to access and analyse the full HELP data base. This new data set includes far-infrared photometry, photometric redshifts, and derived physical properties estimated from modelling the spectral energy distributions over the full HELP sky. All the software and data presented is publicly available.

We have exploited LOFAR deep observations of the Lockman Hole field at 150 MHz to investigate the relation between the radio luminosity of star-forming galaxies (SFGs) and their star formation rates (SFRs), as well as its dependence on stellar mass and redshift. The adopted source classification, SFRs and stellar masses are consensus estimates based on a combination of four different SED fitting methods. We note a flattening of radio spectra of a substantial minority of sources below $\sim 1.4 $ GHz. Such sources have thus a lower radio-loudness level at 150 MHz than expected from extrapolations from 1.4 GHz using the average spectral index. We found a weak trend towards a lower $\hbox{SFR}/L_{150 \rm MHz}$ ratio for higher stellar mass, $M_\star$. We argue that such a trend may account for most of the apparent redshift evolution of the $L_{150 \rm MHz}/\hbox{SFR}$ ratio, in line with previous work. Our data indicate a weaker evolution than found by some previous analyses. We also find a weaker evolution with redshift of the specific star formation rate than found by several (but not all) previous studies. Our radio selection provides a view of the distribution of galaxies in the $\hbox{SFR}$-$M_\star$ plane complementary to that of optical/near-IR selection. It suggests a higher uniformity of the star formation history of galaxies than implied by some analyses of optical/near-IR data. We have derived luminosity functions at 150 MHz of both SFGs and radio-quiet (RQ) AGN at various redshifts. Our results are in very good agreement with the T-RECS simulations and with literature estimates. We also present explicit estimates of SFR functions of SFGs and RQ AGN at several redshifts derived from our radio survey data.

Context. Hyperluminous infrared galaxies (HLIRGs) are shown to have been more abundant in early epochs. The small samples used in earlier studies are not sufficient to draw robust statistical conclusions regarding the physical properties and the power sources of these extreme infrared (IR) bright galaxies. Aims. We make use of multi-wavelength data of a large HLIRG sample to derive the main physical properties, such as stellar mass, star formation rate (SFR), volume density, and the contribution to the cosmic stellar mass density and the cosmic SFR density. We also study the black hole (BH) growth rate and its relationship with the SFR of the host galaxy. Methods. We selected 526 HLIRGs in three deep fields (Boötes, Lockman-Hole, and ELAIS-N1) and adopted two spectral energy distribution (SED) fitting codes: CIGALE, which assumes energy balance, and CYGNUS, which is based on radiative transfer models and does not adopt an energy balance principle. We used two different active galactic nucleus (AGN) models in CIGALE and three AGN models in CYGNUS to compare results that were estimated using different SED fitting codes and a range of AGN models. Results. The stellar mass, total IR luminosity, and AGN luminosity agree well among different models, with a typical median offset of 0.1 dex. The SFR estimates show the largest dispersions (up to 0.5 dex). This dispersion has an impact on the subsequent analysis, which may suggest that the previous contradictory results could partly have been due to the different choices in methods. HLIRGs are ultra-massive galaxies, with 99% of them having stellar masses larger than 10 11 M ⊙ . Our results reveal a higher space density of ultra-massive galaxies than what was found by previous surveys or predicted via simulations. We find that HLIRGs contribute more to the cosmic SFR density as redshift increases. In terms of BH growth, the two SED fitting methods provide different results. We can see a clear trend in whereby SFR decreases as AGN luminosity increases when using CYGNUS estimates. This may possibly imply quenching by AGN in this case, whereas this trend is much weaker when using CIGALE estimates. This difference is also influenced by the dispersion between SFR estimates obtained by the two codes.

Lyman alpha emitters (LAEs) in the Epoch of Reionization (EoR) offer valuable probes of both early galaxy evolution and the process of reionization itself; however, the exact evolution of their abundance and the nature of their emission remain open questions. We combine samples of 229 and 349 LAE candidates at z = 5.7 and z = 6.6, respectively, from the SILVERRUSH narrowband survey with deep Low Frequency Array (LOFAR) radio continuum observations in the European Large Area Infrared Space Observatory Survey-North 1 (ELAIS-N1) field to search for radio galaxies in the EoR and study the low-frequency radio properties of z ≳ 5.7 LAE emitters. Our LOFAR observations reach an unprecedented noise level of ~20 μJy beam−1 at 150 MHz, and we detect five candidate LAEs at >5σ significance. Based on detailed spectral energy distribution modelling of independent multi-wavelength observations in the field, we conclude that these sources are likely [OII] emitters at z = 1.47, yielding no reliable z ≳ 5.7 radio galaxy candidates. We examine the 111 z = 5.7 and z = 6.6 LAE candidates from our panchromatic photometry catalogue that are undetected by LOFAR, finding contamination rates of 81–92% for the z = 5.7 and z = 6.6 subset of the LAE candidate samples. This subset of the full sample is biased towards brighter magnitudes and redder near-infrared colours. The contamination rates of the full sample will therefore likely be lower than the reported values. Contamination of these optically bright LAE samples by likely [OII] emitters is lowered significantly through constraints on the near-infrared colours, highlighting the need for infrared observations to robustly identify bright LAEs in narrowband surveys. Finally, the stacking of radio continuum observations for the robust LAE samples yields 2σ upper limits on radio luminosity of 8.2 × 1023 and 8.7 × 1023 W Hz−1 at z = 5.7 and 6.6, respectively, corresponding to limits on their median star-formation rates of <53 and <56 M⊙ yr−1.

Aims. We aim to study the far-infrared radio correlation (FIRC) at 150 MHz in the local Universe (at a median redshift ⟨z⟩∼0.05) and improve the use of the rest-frame 150 MHz luminosity, L150, as a star-formation rate (SFR) tracer, which is unaffected by dust extinction.Methods. We cross-match the 60 μm selected Revised IRAS Faint Source Survey Redshift (RIFSCz) catalogue and the 150 MHz selected LOFAR value-added source catalogue in the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) Spring Field. We estimate L150 for the cross-matched sources and compare it with the total infrared (IR) luminosity, LIR, and various SFR tracers.Results. We find a tight linear correlation between log L150 and log LIR for star-forming galaxies, with a slope of 1.37. The median qIR value (defined as the logarithm of the LIR to L150 ratio) and its rms scatter of our main sample are 2.14 and 0.34, respectively. We also find that log L150 correlates tightly with the logarithm of SFR derived from three different tracers, i.e., SFRHα based on the Hα line luminosity, SFR60 based on the rest-frame 60 μm luminosity and SFRIR based on LIR, with a scatter of 0.3 dex. Our best-fit relations between L150 and these SFR tracers are, log L150 (L⊙) = 1.35(±0.06) × log SFRHα (M⊙ yr−1) + 3.20(±0.06), log L150 (L⊙) = 1.31(±0.05) × log SFR60 (M⊙ yr−1) + 3.14(±0.06), and log L150 (L⊙) = 1.37 (±0.05) × log SFRIR (M⊙ yr−1) + 3.09(±0.05), which show excellent agreement with each other.