Search

We argue that the difference in infrared-to-radio luminosity ratio between local and high-redshift star-forming galaxies reflects {the alternative physical conditions} -- including magnetic field configurations -- of the dominant population of star-forming galaxies in different redshift ranges. We define three galactic types, based on our reference model, with reference to ages of stellar populations. ``Normal'' late-type galaxies dominate the star formation in the nearby Universe; ``starburst'' galaxies take over at higher redshifts, up to z ~1.5; while ``protospheroidal'' galaxies dominate at high redshift. A reanalysis of data from the COSMOS field combined with literature results shows that, for each population, the data are consistent with an almost redshift-independent mean value of the parameter q_IR, which quantifies the infrared-radio correlation. However, we find a hint of an upturn of the mean q_IR at z>~3.5 consistent with the predicted dimming of synchrotron emission due to cooling of relativistic electrons by inverse Compton scattering off the cosmic microwave background. The typical stellar masses increase from normal, to starburst, and to protospheroidal galaxies, accounting for the reported dependence of the mean q_IR on stellar mass. Higher values of q_IR found for high-z strongly lensed dusty galaxies selected at 500 micron might be explained by differential magnification., Comment: 9 pages, 4 figures, accepted for publication in A&A

We present the first deep (72 hours of observations) radio image of the Euclid Deep Field North (EDFN) obtained with the LOw-Frequency ARray (LOFAR) High Band Antenna (HBA) at 144 MHz. The EDFN is the latest addition to the LOFAR Two-Metre Sky Survey (LoTSS) Deep Fields and these observations represent the first data release for this field. The observations produced a 6" resolution image with a central r.m.s. noise of $32\,\mu$Jy\,beam$^{-1}$. A catalogue of $\sim 23,000$ radio sources above a signal-to-noise ratio (SNR) threshold of 5 is extracted from the inner circular 10 deg$^2$ region. We discuss the data analysis and we provide a detailed description of how we derived the catalogue of radio sources and on the issues related to direction-dependent calibration and their effects on the final products. Finally, we derive the radio source counts at 144 MHz in the EDFN using catalogues of mock radio sources to derive the completeness correction factors. The source counts in the EDFN are consistent with those obtained from the first data release of the other LoTSS Deep Fields (ELAIS-N1, Lockman Hole and Bootes), despite the different method adopted to construct the final catalogue and to assess its completeness., Comment: accepted for publication in A&A

Radio-mode feedback is a key ingredient in galaxy formation and evolution models, required to reproduce the observed properties of massive galaxies in the local Universe. We study the cosmic evolution of radio-AGN feedback out to $z\sim2.5$ using a sample of 9485 radio-excess AGN. We combine the evolving radio luminosity functions with a radio luminosity scaling relationship to estimate AGN jet kinetic powers and derive the cosmic evolution of the kinetic luminosity density, $\Omega_{\rm{kin}}$ (i.e. the volume-averaged heating output). Compared to all radio-AGN, low-excitation radio galaxies (LERGs) dominate the feedback activity out to $z\sim2.5$, with both these populations showing a constant heating output of $\Omega_{\rm{kin}} \approx 4-5 \times 10^{32}\,\rm{W\,Mpc^{-3}}$ across $0.5 < z < 2.5$. We compare our observations to predictions from semi-analytical and hydrodynamical simulations, which broadly match the observed evolution in $\Omega_{\rm{kin}}$, although their absolute normalisation varies. Comparison to the Semi-Analytic Galaxy Evolution (SAGE) model suggests that radio-AGN may provide sufficient heating to offset radiative cooling losses, providing evidence for a self-regulated AGN feedback cycle. We integrate the kinetic luminosity density across cosmic time to obtain the kinetic energy density output from AGN jets throughout cosmic history to be $\sim 10^{50}\,\rm{J\,Mpc^{-3}}$. Compared to AGN winds, the kinetic energy density from AGN jets dominates the energy budget at $z \lesssim 2$; this suggests that AGN jets play an important role in AGN feedback across most of cosmic history., Comment: 15 pages, 6 figures, accepted for publication in MNRAS

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., Comment: Accepted for publication in MNRAS. Catalogues available at www.lofar-surveys.org/deepfields.html. 27 pages

By combining high-sensitivity LOFAR 150MHz, uGMRT 400MHz and 1,250MHz, GMRT 610MHz, and VLA 5GHz data in the ELAIS-N1 field, we study the radio spectral properties of radio-detected star-forming galaxies (SFGs) at observer-frame frequencies of 150-5,000MHz. We select ~3,500 SFGs that have both LOFAR 150MHz and GMRT 610MHz detections, and obtain a median two-point spectral index of $\alpha_{150}^{610}=-0.51\pm0.01$. The photometric redshift of these SFGs spans $z=0.01-6.21$. We also measure the two-point radio spectral indices at 150-400-610-1,250MHz and 150-610-5,000MHz respectively for the GMRT 610-MHz-detected SFGs, and find that, on average, the radio spectrum of SFGs is flatter at low frequency than at high frequency. At observer-frame 150-5,000MHz, we find that the radio spectrum slightly steepens with increasing stellar mass. However, we only find that the radio spectrum flattens with increasing optical depth at $V$-band at $\nu<1$GHz. We suggest that spectral ageing due to the energy loss of CR electrons and thermal free-free absorption could be among the possible main physical mechanisms that drive the above two correlations respectively. In addition, both of these mechanisms could physically explain why the radio spectrum is flatter at low frequency than at high frequency., Comment: 17 pages, 12 figures, 1 table, published in MNRAS

We present the detection of 68 sources from the most sensitive radio survey in circular polarisation conducted to date. We use the second data release of the 144 MHz LOFAR Two-metre Sky Survey to produce circularly-polarised maps with median 140 $\mu$Jy beam$^{-1}$ noise and resolution of 20$''$ for $\approx$27% of the northern sky (5634 deg$^{2}$). The leakage of total intensity into circular polarisation is measured to be $\approx$0.06%, and our survey is complete at flux densities $\geq1$ mJy. A detection is considered reliable when the circularly-polarised fraction exceeds 1%. We find the population of circularly-polarised sources is composed of four distinct classes: stellar systems, pulsars, active galactic nuclei, and sources unidentified in the literature. The stellar systems can be further separated into chromospherically-active stars, M dwarfs, and brown dwarfs. Based on the circularly-polarised fraction and lack of an optical counterpart, we show it is possible to infer whether the unidentified sources are likely unknown pulsars or brown dwarfs. By the completion of this survey of the northern sky, we expect to detect 300$\pm$100 circularly-polarised sources., Comment: Accepted for publication in A&A, 15 pages, 8 figures, 1 table. The catalogue will be publicly available at http://lofar-surveys.org/ and via Vizier shortly

The effect of Active Galactic Nuclei (AGN) on their host galaxies -- in particular their levels of star formation -- remains one of the key outstanding questions of galaxy evolution. Successful cosmological models of galaxy evolution require a fraction of energy released by an AGN to be redistributed into the interstellar medium to reproduce the observed stellar mass and luminosity function and to prevent the formation of over-massive galaxies. Observations have confirmed that the radio-AGN population is energetically capable of heating and redistributing gas at all phases, however, direct evidence of AGN enhancing or quenching star formation remains rare. With modern, deep radio surveys and large integral field spectroscopy (IFS) surveys, we can detect fainter synchrotron emission from AGN jets and accurately probe the star-forming properties of galaxies, respectively. In this paper, we combine data from the LOw Frequency ARray Two-meter Sky Survey with data from one of the largest optical IFS surveys, Mapping Nearby Galaxies at Apache Point Observatory to probe the star-forming properties of 307 local (z $<$ 0.15) galaxies that host radio-detected AGN (RDAGN). We compare our results to a robust control sample of non-active galaxies that each match the stellar mass, redshift, visual morphology, and inclination of a RDAGN host. We find that RDAGN and control galaxies have broad SFR distributions, typically lie below the star-forming main-sequence, and have negative stellar light-weighted age gradients. These results indicate that AGN selected based on their current activity are not responsible for suppressing their host galaxies' star formation. Rather, our results support the maintenance mode role that radio AGN are expected to have in the local Universe., Comment: Accepted for publication in A&A. 24 pages, 14 figures

The \textit{Herschel} Astrophysical Terahertz Large Area Survey (H-ATLAS), that has covered about 642 sq. deg. in 5 bands from 100 to 500 $\mu\rm m$, allows a blind flux-limited selection of blazars at sub-mm wavelengths. However, blazars constitute a tiny fraction of H-ATLAS sources and therefore identifying them is not a trivial task. Using the data on known blazars detected by the H-ATLAS we have defined a locus for 500\,$\mu$m selected blazars and exploited it to select blazar candidates in the H-ATLAS fields. Candidates and known blazars in the H-ATLAS equatorial and South Galactic Pole fields were followed up with the Australia Telescope Compact Array (ATCA) or with the Karl G. Jansky Very Large Array (VLA), and matched with existing radio- and mm-catalogues to reconstruct the spectral behaviour over at least 6 orders of magnitude in frequency. We identified a selection approach that, combining the information in the sub-mm and radio domains, efficiently singles out genuine blazars. In this way, we identified a sample of 39 blazars brighter than $S_{500\mu\rm m} = 35\,$mJy in the H-ATLAS fields. Tests made cross-matching the H-ATLAS catalogues with large catalogues of blazar candidates indicate that the sample is complete. The derived counts are compared with model predictions finding good consistency with the C2Ex model and with estimates based on ALMA data., Comment: 16 pages, 7 figures, accepted by MNRAS

Feedback from low-excitation radio galaxies (LERGs) plays a key role in the lifecycle of massive galaxies in the local Universe; their evolution, and the impact of these active galactic nuclei on early galaxy evolution, however, remain poorly understood. We use a sample of 10481 LERGs from the first data release of the LOFAR Two-meter Sky Survey Deep Fields, covering $\sim$ 25 deg$^2$, to present the first measurement of the evolution of the radio luminosity function (LF) of LERGs out to $z\sim2.5$; this shows relatively mild evolution. We split the LERGs into those hosted by quiescent and star-forming galaxies, finding a new dominant population of LERGs hosted by star-forming galaxies at high redshifts. The incidence of LERGs in quiescent galaxies shows a steep dependence on stellar-mass out to $z \sim1.5$, consistent with local Universe measurements of accretion occurring from cooling of hot gas haloes. The quiescent-LERGs dominate the LFs at $z<1$, showing a strong decline in space density with redshift, tracing that of the available host galaxies, while there is an increase in the characteristic luminosity. The star-forming LERG LF increases with redshift, such that this population dominates the space densities at most radio-luminosities by $z \sim 1$. The incidence of LERGs in star-forming galaxies shows a much weaker stellar-mass dependence, and increases with redshift, suggesting a different fuelling mechanism compared to their quiescent counterparts, potentially associated with the cold gas supply present in the star-forming galaxies., Comment: 26 pages, 16 figures, accepted for publication in MNRAS. LF data available at: https://github.com/rohitk-10/AGN_LF_Kondapally22 . Author list error corrected in updated version

In this data release from the LOFAR Two-metre Sky Survey (LoTSS) we present 120-168MHz images covering 27% of the northern sky. Our coverage is split into two regions centred at approximately 12h45m +44$^\circ$30' and 1h00m +28$^\circ$00' and spanning 4178 and 1457 square degrees respectively. The images were derived from 3,451hrs (7.6PB) of LOFAR High Band Antenna data which were corrected for the direction-independent instrumental properties as well as direction-dependent ionospheric distortions during extensive, but fully automated, data processing. A catalogue of 4,396,228 radio sources is derived from our total intensity (Stokes I) maps, where the majority of these have never been detected at radio wavelengths before. At 6" resolution, our full bandwidth Stokes I continuum maps with a central frequency of 144MHz have: a median rms sensitivity of 83$\mu$Jy/beam; a flux density scale accuracy of approximately 10%; an astrometric accuracy of 0.2"; and we estimate the point-source completeness to be 90% at a peak brightness of 0.8mJy/beam. By creating three 16MHz bandwidth images across the band we are able to measure the in-band spectral index of many sources, albeit with an error on the derived spectral index of +/-0.2 which is a consequence of our flux-density scale accuracy and small fractional bandwidth. Our circular polarisation (Stokes V) 20" resolution 120-168MHz continuum images have a median rms sensitivity of 95$\mu$Jy/beam, and we estimate a Stokes I to Stokes V leakage of 0.056%. Our linear polarisation (Stokes Q and Stokes U) image cubes consist of 480 x 97.6 kHz wide planes and have a median rms sensitivity per plane of 10.8mJy/beam at 4' and 2.2mJy/beam at 20"; we estimate the Stokes I to Stokes Q/U leakage to be approximately 0.2%. Here we characterise and publicly release our Stokes I, Q, U and V images in addition to the calibrated uv-data., Comment: 23 figures, 1 table and 29 pages. The catalogues, images and uv-data associated with this data release are publicly available via https://lofar-surveys.org/

We investigate the radio-far infrared (FIR) correlation for a sample of $28$ bright high-redshift ($1 \lesssim z \lesssim 4$) star-forming galaxies selected in the FIR from the Herschel-ATLAS fields as candidates to be strongly gravitationally lensed. The radio information comes either from high sensitivity dedicated ATCA observations at $2.1$ GHz or from cross-matches with the FIRST survey at $1.4$ GHz. By taking advantage of source brightness possibly enhanced by lensing magnification, we identify a weak evolution with redshift out to $z\lesssim 4$ of the FIR-to-radio luminosity ratio $q_{\rm FIR}$. We also find that the $q_{\rm FIR}$ parameter as a function of the radio power $L_{1.4\,\rm GHz}$ displays a clear decreasing trend, similarly to what is observed for optically/radio selected lensed quasars found in literature, yet covering a complementary region in the $q_{\rm FIR}-L_{1.4\,\rm GHz}$ diagram. We interpret such a behavior in the framework of an in-situ galaxy formation scenario, as a result of the transition from an early dust-obscured star-forming phase (mainly pinpointed by our FIR selection) to a late radio-loud quasar phase (preferentially sampled by the optical/radio selection)., Comment: 12 pages, 6 figures. Minor changes and references updated to match the published version

Radio-loud active galaxies have two accretion modes [radiatively inefficient (RI) and radiatively efficient (RE)], with distinct optical and infrared signatures, and two jet dynamical behaviours, which in arcsec- to arcmin-resolution radio surveys manifest primarily as centre- or edge-brightened structures [Fanaroff-Riley (FR) class I and II]. The nature of the relationship between accretion mode and radio morphology (FR class) has been the subject of long debate. We present a comprehensive investigation of this relationship for a sample of 286 well-resolved radio galaxies in the LOFAR Two-metre Sky Survey Deep Fields (LoTSS-Deep) first data release, for which robust morphological and accretion mode classifications have been made. We find that two-thirds of luminous FRII radio galaxies are RI, and identify no significant differences in the visual appearance or source dynamic range (peak/mean surface brightness) of the RI and RE FRIIs, demonstrating that both RI and RE systems can produce FRII structures. We also find a significant population of low-luminosity FRIIs (predominantly RI), supporting our earlier conclusion that FRII radio structures can be produced at all radio luminosities. We demonstrate that in the luminosity range where both morphologies are present, the probability of producing FRI or FRII radio morphology is directly linked to stellar mass, while across all morphologies and luminosities, RE accretion occurs in systems with high specific star formation rate, presumably because this traces fuel availability. In summary, the relationship between accretion mode and radio morphology is very indirect, with host-galaxy environment controlling these two key parameters in different ways., Comment: 22 pages, 17 figures, 7 tables. Accepted for publication in MNRAS. Update with minor changes to match published version

The identification of an emission line is unambiguous when multiple spectral features are clearly visible in the same spectrum. However, in many cases, only one line is detected, making it difficult to correctly determine the redshift. We developed a freely available unsupervised machine-learning algorithm based on unbiased topology (UMLAUT) that can be used in a very wide variety of contexts, including the identification of single emission lines. To this purpose, the algorithm combines different sources of information, such as the apparent magnitude, size and color of the emitting source, and the equivalent width and wavelength of the detected line. In each specific case, the algorithm automatically identifies the most relevant ones (i.e., those able to minimize the dispersion associated with the output parameter). The outputs can be easily integrated into different algorithms, allowing us to combine supervised and unsupervised techniques and increasing the overall accuracy. We tested our software on WISP (WFC3 IR Spectroscopic Parallel) survey data. WISP represents one of the closest existing analogs to the near-IR spectroscopic surveys that are going to be performed by the future Euclid and Roman missions. These missions will investigate the large-scale structure of the universe by surveying a large portion of the extragalactic sky in near-IR slitless spectroscopy, detecting a relevant fraction of single emission lines. In our tests, UMLAUT correctly identifies real lines in 83.2% of the cases. The accuracy is slightly higher (84.4%) when combining our unsupervised approach with a supervised approach we previously developed.

We present the first sub-mJy ($\approx0.7$ mJy beam$^{-1}$) survey to be completed below 100 MHz, which is over an order of magnitude deeper than previously achieved for widefield imaging of any field at these low frequencies. The high resolution ($15 \times 15$ arcsec) image of the Bo\"otes field at 34-75 MHz is made from 56 hours of observation with the LOw Frequency ARray (LOFAR) Low Band Antenna (LBA) system. The observations and data reduction, including direction-dependent calibration, are described here. We present a radio source catalogue containing 1,948 sources detected over an area of $23.6$ deg$^2$, with a peak flux density threshold of $5\sigma$. Using existing datasets, we characterise the astrometric and flux density uncertainties, finding a positional uncertainty of $\sim1.2$ arcsec and a flux density scale uncertainty of about 5 per cent. Using the available deep 144-MHz data, we identified 144-MHz counterparts to all the 54-MHz sources, and produced a matched catalogue within the deep optical coverage area containing 829 sources. We calculate the Euclidean-normalised differential source counts and investigate the low-frequency radio source spectral indices between 54 and 144 MHz, both of which show a general flattening in the radio spectral indices for lower flux density sources, from $\sim-0.75$ at 144-MHz flux densities between 100-1000 mJy to $\sim-0.5$ at 144-MHz flux densities between 5-10 mJy, due to a growing population of star forming galaxies and compact core-dominated AGN., Comment: 15 pages, 17 figures. Accepted for publication by A&A. Data available online at lofar-surveys.org

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., Comment: 18 pages, 9 figures, 2 tables, A&A accepted

The Additional Representative Images for Legacy (ARI-L) project is a European Development project for ALMA Upgrade approved by the Joint ALMA Observatory (JAO) and the European Southern Observatory (ESO), started in June 2019. It aims to increase the legacy value of the ALMA Science Archive (ASA) by bringing the reduction level of ALMA data from Cycles 2-4 close to that of data from more recent Cycles processed for imaging with the ALMA Pipeline. As of mid-2021 more than 150000 images have been returned to the ASA for public use. At its completion in 2022, the project will have provided enhanced products for at least 70% of the observational data from Cycles 2-4 processable with the ALMA Pipeline. In this paper we present the project rationale, its implementation, and the new opportunities offered to ASA users by the ARI-L products. The ARI-L cubes and images complement the much limited number of archival image products generated during the data quality assurance stages (QA2), which cover only a small fraction of the available data for those Cycles. ARI-L imaging products are highly relevant for many science cases and significantly enhance the possibilities for exploiting archival data. Indeed, ARI-L products facilitate archive access and data usage for science purposes even for non-expert data miners, provide a homogeneous view of all data for better dataset comparisons and download selections, make the archive more accessible to visualization and analysis tools, and enable the generation of preview images and plots similar to those possible for subsequent Cycles., Comment: 15 pages. Accepted for publication in PASP

We make use of multi-wavelength data of a large hyperluminous infrared (HLIRG) sample to derive their main physical properties, e.g., stellar mass, star-formation rate (SFR), volume density, contribution to the cosmic stellar mass density and to the cosmic SFR density. We also study the black hole (BH) growth rate and its relationship with the SFR of the host galaxy. We select 526 HLIRGs in three deep fields (Bo$\"o$tes, Lockman-Hole, ELAIS-N1) and adopt 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 energy balance principle. We use two different active galactic nucleus (AGN) models in CIGALE and three AGN models in CYGNUS to compare the results estimated using different SED fitting codes and different AGN models. The stellar mass, total IR luminosity and AGN luminosity agree well between 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 previous contradictory results could partly be due to different choices of methods. HLIRGs are ultra-massive galaxies with 99% of them having stellar masses larger than $10^{11} M_{\odot}$. Our results reveal a higher space density of ultra-massive galaxies than found in previous surveys or predicted by 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 which SFR decreases as AGN luminosity increases when using CYGNUS estimates, possibly implying quenching by AGN, while this trend is much weaker when using CIGALE estimates. This difference is also influenced by the dispersion between SlFR estimates obtained by the two codes., Comment: 21 pages, 15 figures. Accepted in A&A

We investigate the radio properties of a sample of 53 sources selected at 850 $\mu$m from the SCUBA-2 Cosmology Legacy Survey using new deep, low-frequency radio imaging of the Lockman Hole field from the Low Frequency Array. Combining these data with additional radio observations from the GMRT and the JVLA, we find a variety of radio spectral shapes and luminosities within our sample despite their similarly bright submillimetre flux densities. 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. 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 $\lesssim$ 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) submm sources undetected at all other wavelengths. We speculate that these objects may be very high redshift sources, likely residing at z > 4., Comment: 15 pages, 10 figures, accepted by A&A

Fundamental differences in the radio properties of red quasars (QSOs), as compared to blue QSOs, have been recently discovered, positioning them as a potential key population in the evolution of galaxies and black holes across cosmic time. To elucidate their nature, we exploited a rich compilation of photometry and spectroscopic data to model their spectral energy distributions (SEDs) from the UV to the FIR and characterise their emission-line properties. Following a systematic comparison approach, we infer the AGN accretion, obscuration, and host galaxy properties in a sample of ~1800 QSOs at 0.21000 km/s) in red QSOs. We find that red QSOs that exhibit evidence for high-velocity winds present a stronger signature of the infrared excess, suggesting a causal connection between reddening and the presence of hot dust in QSO winds. We propose that dusty winds at nuclear scales are potentially the physical ingredient responsible for the colours in red QSOs, as well as a key parameter for the regulation of accretion material in the nucleus., Comment: Accepted for publication in Astronomy & Astrophysics. 22 pages, 16 figures. Abstract abridged

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 > 10^13 L_solar, 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. We cross-match Herschel sources and LOFAR sources in Bootes (8.63 deg^2), Lockman Hole (10.28 deg^2), and ELAIS-N1 (6.74 deg^2) with rms sensitivities of around 32, 22, and 20 mJy per beam, 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 (SED) fitting, combined with multi-wavelength counterpart identifications and photometric redshift estimates, to derive IR luminosities. The depth of the LOFAR data allows us to identify highly complete (around 92% completeness) samples of bright Herschel sources with a simple selection based on the 250 micron flux density (45, 40, and 35 mJy in Bootes, 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 around 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 around 5 to 18 / deg^2 (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., Comment: 16 pages, 11 figures, accepted for publication in A&A as part of the LOFAR Deep Fields Paper Splash

With the advent of new generation low-frequency telescopes, such as the LOw Frequency ARray (LOFAR), and improved calibration techniques, we have now started to unveil the sub GHz radio sky with unprecedented depth and sensitivity. The LOFAR Two Meter Sky Survey (LoTSS) is an ongoing project in which the whole northern radio sky will be observed at 150 MHz with a sensitivity better than 100 $\mu$Jy beam$^{-1}$ at a resolution of \asec{6}. Additionally, deeper observations are planned to cover smaller areas with higher sensitivity. The Lockman Hole, the Bo\"otes and the Elais-N1 regions are among the most well known northern extra-galactic fields, and the deepest of the LoTSS Deep Fields so far. We exploit these deep observations to derive the deepest radio source counts at 150~MHz to date. Our counts are in broad agreement with those from the literature, and show the well known upturn at $\leq$ few mJy, mainly associated with the emergence of the star-forming galaxy population. More interestingly, our counts show for the first time a very pronounced drop around S$\sim$2 mJy, which results in a prominent `bump' at sub-mJy flux densities. Such a feature was not observed in previous counts' determinations (neither at 150 MHz nor at higher frequency). While sample variance can play a role in explaining the observed discrepancies, we believe this is mostly the result of a careful analysis aimed at deblending confused sources and removing spurious sources and artifacts from the radio catalogues. This `drop and bump' feature cannot be reproduced by any of the existing state-of-the-art evolutionary models, and appears to be associated with a deficiency of AGN at intermediate redshift ($1

The Low Frequency Array (LOFAR) is an ideal instrument to conduct deep extragalactic surveys. It has a large field of view and is sensitive to large scale and compact emission. It is, however, very challenging to synthesize thermal noise limited maps at full resolution, mainly because of the complexity of the low-frequency sky and the direction dependent effects (phased array beams and ionosphere). In this first paper of a series we present a new calibration and imaging pipeline that aims at producing high fidelity, high dynamic range images with LOFAR High Band Antenna data, while being computationally efficient and robust against the absorption of unmodeled radio emission. We apply this calibration and imaging strategy to synthesize deep images of the Bootes and LH fields at 150 MHz, totaling $\sim80$ and $\sim100$ hours of integration respectively and reaching unprecedented noise levels at these low frequencies of $\lesssim30$ and $\lesssim23$ $\mu$Jy/beam in the inner $\sim3$ deg$^2$. This approach is also being used to reduce the LoTSS-wide data for the second data release., Comment: This paper is part of the 1st data release of the LoTSS Deep Fields. Electronic data catalogues will be made available on journal publication

Low-frequency radio observations are revealing an increasing number of diffuse synchrotron sources from galaxy clusters, dominantly in the form of radio halos or radio relics. The existence of this diffuse synchrotron emission indicates the presence of relativistic particles and magnetic fields. It is still an open question what mechanisms exactly are responsible for the population of relativistic electrons driving this synchrotron emission. The LOFAR Two-metre Sky Survey Deep Fields offer a unique view of this problem. Reaching noise levels below 30 $\mu$Jy/beam, these are the deepest images made at the low frequency of 144 MHz. This paper presents a search for diffuse emission in galaxy clusters in the first data release of the LOFAR Deep Fields. We detect a new high-redshift radio halo with a flux density of $8.9 \pm 1.0$ mJy and corresponding luminosity of $P_{144\mathrm{MHz}}=(3.6 \pm 0.6)\times10^{25}$ W Hz$^{-1}$ in an X-ray detected cluster at $z=0.77$ with a mass estimate of $M_{500} = 3.3_{-1.7}^{+1.1} \times 10^{14} M_\odot.$ Deep upper limits are placed on clusters with non-detections. We compare the results to the correlation between halo luminosity and cluster mass derived for radio halos found in the literature. This study is one of few to find diffuse emission in low mass ($M_{500} < 5\times10^{14} M_\odot$) systems and shows that deep low-frequency observations of galaxy clusters are fundamental for opening up a new part of parameter space in the study of non-thermal phenomena in galaxy clusters., Comment: This paper is part of the 1st data release of the LoTSS Deep Fields. Re-submitted paper (reply to referee)

The LOFAR Two-metre Sky Survey (LoTSS) will cover the full northern sky and, additionally, aims to observe the LoTSS deep fields to a noise level of ~10 microJy/bm over several tens of square degrees in areas that have the most extensive ancillary data. This paper presents the ELAIS-N1 deep field, the deepest of the LoTSS deep fields to date. With an effective observing time of 163.7 hours, it reaches a root mean square (RMS) noise level below 20 microJy/bm in the central region (and below 30 microJy/bm over 10 square degrees). The resolution is 6 arcsecs and 84862 radio sources were detected in the full area (68 sq. deg.) with 74127 sources in the highest quality area at less than 3 degrees from the pointing centre. The observation reaches a sky density of more than 5000 sources per sq. deg. in the central ~5 sq. deg. region. We present the calibration procedure, which addresses the special configuration of some observations and the extended bandwidth covered (115 to 177 MHz; central frequency 146.2 MHz) compared to standard LoTSS. We also describe the methods used to calibrate the flux density scale using cross-matching with sources detected by other radio surveys in the literature. We find the flux density uncertainty related to the flux density scale to be ~6.5%. By studying the variations of the flux density measurements between different epochs, we show that relative flux density calibration is reliable out to about a 3 degree radius, but that additional flux density uncertainty is present for all sources at about the 3 per cent level; this is likely to be associated with residual calibration errors, and is shown to be more significant in datasets with poorer ionosphere conditions. We also provide intra-band spectral indices, which can be useful to detect sources with unusual spectral properties. The final uncertainty in the flux densities is estimated to be ~10% for ELAIS-N1., Comment: Accepted for publication in A&A. This paper is part of the 1st data release of the LoTSS Deep Fields. Electronic data catalogues will be made available on journal publication. 21 pages, 17 figures

Lyman alpha emitters (LAEs) in the Epoch of Reionization (EoR) offer valuable probes of early galaxy evolution and the process of reionization; 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 ELAIS-N1 field to search for radio galaxies in the EoR and study the low-frequency radio properties of $z\gtrsim5.7$ LAE emitters. Our LOFAR observations reach an unprecedented noise level of $\sim20\,\mu$Jy beam$^{-1}$ at 150MHz, and we detect five candidate LAEs at $>5\sigma$ significance. Based on detailed spectral energy distribution modelling of independent multi-wavelength observations, we conclude that these sources are likely [OII] emitters at $z=1.47$, yielding no reliable $z\gtrsim5.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 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 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$\sigma$ upper limits on radio luminosity of 8.2$\times$10$^{23}$ and 8.7$\times$10$^{23}$ 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$_{\odot}$ yr$^{-1}$., Comment: 17 pages, 11 figures, accepted for publication in A&A. The abstract is abridged. This paper is part of the 1st data release of the LoTSS Deep Fields

The Low Frequency Array (LOFAR) Two-metre Sky Survey (LoTSS) is a sensitive, high-resolution 120-168 MHz survey split across multiple tiers over the northern sky. The first LoTSS Deep Fields data release consists of deep radio continuum imaging at 150 MHz of the Bo\"{o}tes, European Large Area Infrared Space Observatory Survey-North 1 (ELAIS-N1), and Lockman Hole fields, down to rms sensitivities of $\sim$32, 20, and 22 $\mu$Jy beam$^{-1}$, respectively. In this paper we present consistent photometric redshift (photo-$z$) estimates for the optical source catalogues in all three fields - totalling over 7 million sources ($\sim5$ million after limiting to regions with the best photometric coverage). Our photo-$z$ estimation uses a hybrid methodology that combines template fitting and machine learning and is optimised to produce the best possible performance for the radio continuum selected sources and the wider optical source population. Comparing our results with spectroscopic redshift samples, we find a robust scatter ranging from 1.6 to 2% for galaxies and 6.4 to 7% for identified optical, infrared, or X-ray selected active galactic nuclei (AGN). Our estimated outlier fractions ($\left | z_{\text{phot}} - z_{\text{spec}} \right | / (1+z_{\text{spec}}) > 0.15$) for the corresponding subsets range from 1.5 to 1.8% and 18 to 22%, respectively. Replicating trends seen in analyses of previous wide-area radio surveys, we find no strong trend in photo-$z$ quality as a function of radio luminosity for a fixed redshift. We exploit the broad wavelength coverage available within each field to produce galaxy stellar mass estimates for all optical sources at $z < 1.5$. Stellar mass functions derived for each field are used to validate our mass estimates, with the resulting estimates in good agreement between each field and with published results from the literature., Comment: 23 pages, 14 figures. Accepted for publication in A&A. This paper is part of the 1st data release of the LoTSS Deep Fields - Electronic data catalogues will be made available on journal publication

In this paper, we investigate the relationship between 150MHz luminosity and star formation rate (the SFR-L150 relation) using 150MHz measurements for a near-infrared selected sample of 118,517 $z<1$ galaxies. New radio survey data offer compelling advantages for studying star formation in galaxies, with huge increases in sensitivity, survey speed and resolution over previous generation surveys, and remaining impervious to extinction. The LOFAR Surveys Key Science Project is transforming our understanding of the low-frequency radio sky, with the 150MHz data over the ELAIS-N1 field reaching an RMS sensitivity of 20uJy/beam over 10 deg$^2$ at 6" resolution. All of the galaxies studied have SFR and stellar mass estimates derived from energy balance SED fitting, using redshifts and aperture-matched forced photometry from the LOFAR Two-metre Sky Survey (LoTSS) deep fields data release. The impact of active galactic nuclei is minimised by leveraging the deep ancillary data alongside outlier-resistant median-likelihood methods. We find a linear and non-evolving SFR-L150 relation, apparently consistent with expectations based on calorimetric arguments, down to the lowest SFRs. However, we also recover compelling evidence for stellar mass dependence in line with previous work on this topic, in the sense that higher mass galaxies have a larger 150MHz luminosity at a given SFR, suggesting that the overall agreement with calorimetric arguments may be a coincidence. We conclude that in the absence of AGN, 150MHz observations can be used to measure accurate galaxy SFRs out to $z=1$ at least, but it is necessary to account for stellar mass in order to obtain 150MHz-derived SFRs accurate to <0.5 dex. Our best-fit relation is $\log_{10} (L_\mathrm{150 MHz} / W\,Hz^{-1}) = (0.90\pm 0.01) \log_{10}(\psi/M_\odot\,\mathrm{yr}^{-1}) + (0.33 \pm 0.04) \log_{10} (M/10^{10}M_\odot) + 22.22 \pm 0.02$. (Abridged), Comment: 17 pages, 17 figures, accepted for publication in A&A. This paper is part of the 1st data release of the LoTSS Deep Fields

We present the source associations, cross-identifications, and multi-wavelength properties of the faint radio source population detected in the deep tier of the LOFAR Two Metre Sky Survey (LoTSS): the LoTSS Deep Fields. The first LoTSS Deep Fields data release consists of deep radio imaging at 150~MHz of the ELAIS-N1, Lockman Hole, and Bo\"{o}tes fields, down to RMS sensitives of around 20, 22, and 32$~\mu$Jy\,beam$^{-1}$, respectively. These fields are some of the best studied extra-galactic fields in the northern sky, with existing deep, wide-area panchromatic photometry from X-ray to infrared wavelengths, covering a total of $\approx$~26~\mbox{deg$^{2}$}. We first generated improved multi-wavelength catalogues in ELAIS-N1 and Lockman Hole; combined with the existing catalogue for Bo\"{o}tes, we present forced, matched aperture photometry for over 7.2 million sources across the three fields. We identified multi-wavelength counterparts to the radio detected sources, using a combination of the Likelihood Ratio method and visual classification, which greatly enhances the scientific potential of radio surveys and allows for the characterisation of the photometric redshifts and the physical properties of the host galaxies. The final radio-optical cross-matched catalogue consists of 81\,951 radio-detected sources, with counterparts identified and multi-wavelength properties presented for 79\,820 ($>$97\%) sources. We also examine the properties of the host galaxies, and through stacking analysis find that the radio population with no identified counterpart is likely dominated by AGN at $z\sim3-4$. This dataset contains one of the largest samples of radio-selected star-forming galaxies and active galactic nuclei (AGN) at these depths, making it ideal for studying the history of star-formation, and the evolution of galaxies and AGN across cosmic time., Comment: 34 pages, 14 figures and 14 tables. Accepted for publication in A&A special issue. This paper is part of the 1st data release of the LoTSS Deep Fields. Electronic data catalogues will be made available on journal publication

The morphology of quiescent galaxies has been found to be correlated with the activity of their central super massive black hole. In this work, we use data from the first data release of the LOFAR Two$-$Metre Sky Survey (LoTSS DR1) and the Sloan Digital Sky Survey Data Release 7 (SDSS DR7) to select more than 15 000 quiescent galaxies at $z<0.3$ to investigate the connection between radio-loud active galactic nuclei (RLAGNs) and the morphology of their host galaxy. Taking advantage of the depth of LoTSS, we find that the fraction of RLAGNs with $L_{\rm 150\,MHz}>10^{21}\rm\,W\,Hz^{-1}$ at fixed stellar mass, velocity dispersion, or surface mass density does not depend on the galaxy projected axis ratio ($q$). However, the high-power ($L_{\rm 150\,MHz}>10^{23}\rm\,W\,Hz^{-1}$) RLAGNs are more likely to be found in massive, round galaxies, while the low- and intermediate-power ($L_{\rm 150\,MHz}\leq10^{23}\rm\,W\,Hz^{-1}$) RLAGNs have similar distributions of $q$ to non-RLAGN galaxies. We argue that our results support the picture that high-power RLAGNs are more easily triggered in galaxies with a merger-rich history, while low-power RLAGNs can be triggered in galaxies growing mainly via secular processes. Our work also supports the idea that the low-luminosity RLAGN may be sufficient for maintenance-mode feedback in low-mass quiescent galaxies with disc-like morphology, which is based on a simple extrapolation from the observed energy balance between cooling and RLAGN-induced cavities in massive clusters. We find no significant difference between the $q$ distributions of RLAGNs likely to be found in clusters and those likely not found in clusters after controlling the radio luminosity and stellar mass of the two samples, indicating that the environment does not significantly influence the morphology--RLAGN correlation., Comment: 14 pages, 8 figures, accepted for publication in A&A

As the largest radio telescope in the world, the Square Kilometre Array (SKA) will lead the next generation of radio astronomy. The feats of engineering required to construct the telescope array will be matched only by the techniques developed to exploit the rich scientific value of the data. To drive forward the development of efficient and accurate analysis methods, we are designing a series of data challenges that will provide the scientific community with high-quality datasets for testing and evaluating new techniques. In this paper we present a description and results from the first such Science Data Challenge (SDC1). Based on SKA MID continuum simulated observations and covering three frequencies (560 MHz, 1400MHz and 9200 MHz) at three depths (8 h, 100 h and 1000 h), SDC1 asked participants to apply source detection, characterization and classification methods to simulated data. The challenge opened in November 2018, with nine teams submitting results by the deadline of April 2019. In this work we analyse the results for 8 of those teams, showcasing the variety of approaches that can be successfully used to find, characterise and classify sources in a deep, crowded field. The results also demonstrate the importance of building domain knowledge and expertise on this kind of analysis to obtain the best performance. As high-resolution observations begin revealing the true complexity of the sky, one of the outstanding challenges emerging from this analysis is the ability to deal with highly resolved and complex sources as effectively as the unresolved source population., Comment: 18 pages, 11 figures, accepted by MNRAS

Future surveys focusing on understanding the nature of dark energy (e.g., Euclid and WFIRST) will cover large fractions of the extragalactic sky in near-IR slitless spectroscopy. These surveys will detect a large number of galaxies that will have only one emission line in the covered spectral range. In order to maximize the scientific return of these missions, it is imperative that single emission lines are correctly identified. Using a supervised machine-learning approach, we classified a sample of single emission lines extracted from the WFC3 IR Spectroscopic Parallel survey (WISP), one of the closest existing analogs to future slitless surveys. Our automatic software integrates a SED fitting strategy with additional independent sources of information. We calibrated it and tested it on a "gold" sample of securely identified objects with multiple lines detected. The algorithm correctly classifies real emission lines with an accuracy of 82.6%, whereas the accuracy of the SED fitting technique alone is low (~50%) due to the limited amount of photometric data available (<=6 bands). While not specifically designed for the Euclid and WFIRST surveys, the algorithm represents an important precursor of similar algorithms to be used in these future missions.

Context:Blazars are the rarest and most powerful active galactic nuclei, playing a crucial and growing role in today multi-frequency and multi-messenger astrophysics. Current blazar catalogs, however, are incomplete and particularly depleted at low Galactic latitudes. Aims: We aim at augmenting the current blazar census to build a catalog of blazar candidates with homogeneous sky coverage that can provide candidate counterparts to unassociated gamma-ray sources, sources of high-energy neutrino emission, and ultra-high energy cosmic rays. Methods: Starting from the ALMA Calibrator Catalog we built a catalog of 1580 blazar candidates (ALMA Blazar Candidates, ABC) for which we collect multi-wavelength information. We also compared ABC sources with existing blazar catalogs. Results: The ABC catalogue fills the lack of low Galactic latitude sources in current blazar catalogues. ABC sources are significantly dimmer than known blazars in Gaia g band, and they appear bluer in SDSS and WISE colors. The majority of ABC sources (~ 90%) have optical spectra that classify them as QSO, while the remaining sources resulted galactic objects. ABC sources are similar in X-rays to known blazar, while in gamma-rays they are on average dimmer and softer, indicating a significant contribution of FSRQ sources. Making use of WISE colours, we classified 715 ABC sources as candidate gamma-ray blazar of different classes. Conclusions: We built a new catalogue of 1580 candidate blazars with a rich multi-wavelength data-set, filling the lack of low Galactic latitude sources in current blazar catalogues. This will be particularly important to identify the source population of high energy neutrinos or ultra-high energy cosmic rays. The data collected by the upcoming LSST surveys will provide a key tool to investigate the possible blazar nature of these sources., Comment: 53 pages, 32 figures, 6 tables. Accepted for publication on A&A

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-um 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., SFR_Halpha based on the Halpha line luminosity, SFR_60 based on the rest-frame 60-um luminosity and SFR_IR based on LIR, with a scatter of 0.3 dex. Our best-fit relations between L150 and these SFR tracers are, log L150 (Lsun) = 1.35(0.06) x log SFR_Halpha (Msun/yr) + 3.20(0.06), log L150 (Lsun) = 1.31(0.05) x log SFR_60 (Msun/yr) + 3.14(0.06), and log L150 (Lsun) = 1.37(0.05) x log SFR_IR (Msun/yr) + 3.09(0.05), which show excellent agreement with each other., Comment: 13 pages, 13 figures, accepted for publication in A&A

The Probe of Inflation and Cosmic Origins (PICO) is a proposed probe-scale space mission consisting of an imaging polarimeter operating in frequency bands between 20 and 800 GHz. We describe the science achievable by PICO, which has sensitivity equivalent to more than 3300 Planck missions, the technical implementation, the schedule and cost., Comment: APC White Paper submitted to the Astro2020 decadal panel; 10 page version of the 50 page mission study report arXiv:1902.10541

Planck, SPT and ACT surveys have clearly demonstrated that Cosmic Microwave Background (CMB) experiments, while optimised for cosmological measurements, have made important contributions to the field of extragalactic astrophysics in the last decade. Future CMB experiments have the potential to make even greater contributions. One example is the detection of high-z galaxies with extreme gravitational amplifications. The combination of flux boosting and of stretching of the images has allowed the investigation of the structure of galaxies at z ~3 with the astounding spatial resolution of about 60 pc. Another example is the detection of proto-clusters of dusty galaxies at high z when they may not yet possess the hot intergalactic medium allowing their detection in X-rays or via the Sunyaev-Zeldovich effect. Next generation CMB experiments, like PICO, CORE, CMB-Bharat from space and Simons Observatory and CMB-S4 from the ground, will discover several thousands of strongly lensed galaxies out to z~6 or more and of galaxy proto-clusters caught in the phase when their member galaxies where forming the bulk of their {stars. They will also detect tens of thousands of local dusty galaxies and thousands of radio sources at least up to z~5. Moreover they will measure the polarized emission of thousands of radio sources and of dusty galaxies at mm/sub-mm wavelengths., Comment: 29 pages, 10 figures. Accepted for publication in Frontiers in Astronomy and Space Sciences, section Cosmology

We present a new catalogue of ALMA observations of 3,364 bright, compact radio sources, mostly blazars, used as calibrators. These sources were observed between May 2011 and July 2018, for a total of 47,115 pointings in different bands and epochs. We have exploited the ALMA data to validate the photometry given in the new Planck Multi-frequency Catalogue of Non-thermal sources (PCNT), for which an external validation was not possible so far. We have also assessed the positional accuracy of Planck catalogues and the PCNT completeness limits, finding them to be consistent with those of the Second Planck Catalogue of Compact Sources. The ALMA continuum spectra have allowed us to extrapolate the observed radio source counts at 100 GHz to the effective frequencies of ALMA bands 4, 6, 7, 8 and 9 (145, 233, 285, 467 and 673 GHz, respectively), where direct measurements are scanty, especially at the 3 highest frequencies. The results agree with the predictions of the Tucci et al. model C2Ex, while the model C2Co is disfavoured., Comment: 9 pages, 4 figures, 4 tables, accepted for publication in MNRAS

We present a catalogue of ALMA flux density measurements of 754 calibrators observed between August 2012 and September 2017, for a total of 16,263 observations in different bands and epochs. The flux densities were measured reprocessing the ALMA images generated in the framework of the ALMACAL project, with a new code developed by the Italian node of the European ALMA Regional Centre. A search in the online databases yielded redshift measurements for 589 sources ($\sim$78 per cent of the total). Almost all sources are flat-spectrum, based on their low-frequency spectral index, and have properties consistent with being blazars of different types. To illustrate the properties of the sample we show the redshift and flux density distributions as well as the distributions of the number of observations of individual sources and of time spans in the source frame for sources observed in bands 3 (84$-$116 GHz) and 6 (211$-$275 GHz). As examples of the scientific investigations allowed by the catalogue we briefly discuss the variability properties of our sources in ALMA bands 3 and 6 and the frequency spectra between the effective frequencies of these bands. We find that the median variability index steadily increases with the source-frame time lag increasing from 100 to 800 days, and that the frequency spectra of BL Lacs are significantly flatter than those of flat-spectrum radio quasars. We also show the global spectral energy distributions of our sources over 17 orders of magnitude in frequency., Comment: 8 pages, 9 figures, 2 tables, MNRAS accepted

For a sample of star forming galaxies in the redshift interval 0.15$<$z$<$0.3, we study how both the relative strength of the AGN infra-red emission, compared to that due to the star formation (SF), and the numerical fraction of AGNs, change as a function of the total stellar mass of the hosting galaxy group (M$^{*}_{\mathrm{group}}$), between $10^{10.25}$ and $10^{11.9}$M$_{\odot}$. Using a multi-component SED fitting analysis, we separate the contribution of stars, AGN torus and star formation to the total emission at different wavelengths. This technique is applied to a new multi-wavelength data-set in the SIMES field (23 not redundant photometric bands), spanning the wavelength range from the UV (GALEX) to the far-IR (Herschel) and including crucial AKARI and WISE mid-IR observations (4.5 \mu m$<\lambda<$24 \mu m), where the BH thermal emission is stronger. This new photometric catalog, that includes our best photo-z estimates, is released through the NASA/IPAC Infrared Science Archive (IRSA). Groups are identified through a friends of friends algorithm ($\sim$62% purity, $\sim$51% completeness). We identified a total of 45 galaxies requiring an AGN emission component, 35 of which in groups and 10 in the field. We find BHAR$\propto ($M$^{*}_{\mathrm{group}})^{1.21\pm0.27}$ and (BHAR/SFR)$\propto ($M$^{*}_{\mathrm{group}})^{1.04\pm0.24}$ while, in the same range of M$^{*}_{\mathrm{group}}$, we do not observe any sensible change in the numerical fraction of AGNs. Our results indicate that the nuclear activity (i.e. the BHAR and the BHAR/SFR ratio) is enhanced when galaxies are located in more massive and richer groups., Comment: 31 pages, 23 figures

We present an analysis of the main systematic effects that could impact the measurement of CMB polarization with the proposed CORE space mission. We employ timeline-to-map simulations to verify that the CORE instrumental set-up and scanning strategy allow us to measure sky polarization to a level of accuracy adequate to the mission science goals. We also show how the CORE observations can be processed to mitigate the level of contamination by potentially worrying systematics, including intensity-to-polarization leakage due to bandpass mismatch, asymmetric main beams, pointing errors and correlated noise. We use analysis techniques that are well validated on data from current missions such as Planck to demonstrate how the residual contamination of the measurements by these effects can be brought to a level low enough not to hamper the scientific capability of the mission, nor significantly increase the overall error budget. We also present a prototype of the CORE photometric calibration pipeline, based on that used for Planck, and discuss its robustness to systematics, showing how CORE can achieve its calibration requirements. While a fine-grained assessment of the impact of systematics requires a level of knowledge of the system that can only be achieved in a future study phase, the analysis presented here strongly suggests that the main areas of concern for the CORE mission can be addressed using existing knowledge, techniques and algorithms., Comment: 54 pages, 26 figures, 3 tables

Future observations of cosmic microwave background (CMB) polarisation have the potential to answer some of the most fundamental questions of modern physics and cosmology. In this paper, we list the requirements for a future CMB polarisation survey addressing these scientific objectives, and discuss the design drivers of the CORE space mission proposed to ESA in answer to the "M5" call for a medium-sized mission. The rationale and options, and the methodologies used to assess the mission's performance, are of interest to other future CMB mission design studies. CORE is designed as a near-ultimate CMB polarisation mission which, for optimal complementarity with ground-based observations, will perform the observations that are known to be essential to CMB polarisation scienceand cannot be obtained by any other means than a dedicated space mission., Comment: 79 pages, 14 figures

Observational investigations of the abundance of massive precursors of local galaxy clusters ("proto-clusters") allow us to test the growth of density perturbations, to constrain cosmological parameters that control it, to test the theory of non-linear collapse and how the galaxy formation takes place in dense environments. The Planck collaboration has recently published a catalogue of >~ 2000 cold extra-galactic sub-millimeter sources, i.e. with colours indicative of z >~ 2, almost all of which appear to be over-densities of star-forming galaxies. They are thus considered as proto-cluster candidates. Their number densities (or their flux densities) are far in excess of expectations from the standard scenario for the evolution of large-scale structure. Simulations based on a physically motivated galaxy evolution model show that essentially all cold peaks brighter than S_{545GHz} = 500 mJy found in Planck maps after having removed the Galactic dust emission can be interpreted as positive Poisson fluctuations of the number of high-z dusty proto-clusters within the same Planck beam, rather then being individual clumps of physically bound galaxies. This conclusion does not change if an empirical fit to the luminosity function of dusty galaxies is used instead of the physical model. The simulations accurately reproduce the statistic of the Planck detections and yield distributions of sizes and ellipticities in qualitative agreement with observations. The redshift distribution of the brightest proto-clusters contributing to the cold peaks has a broad maximum at 1.5 <~ z <~ 3. Therefore follow-up of Planck proto-cluster candidates will provide key information on the high-z evolution of large scale structure., Comment: 10 pages, 8 figures, accepted for publication in MNRAS

We describe a space-borne, multi-band, multi-beam polarimeter aiming at a precise and accurate measurement of the polarization of the Cosmic Microwave Background. The instrument is optimized to be compatible with the strict budget requirements of a medium-size space mission within the Cosmic Vision Programme of the European Space Agency. The instrument has no moving parts, and uses arrays of diffraction-limited Kinetic Inductance Detectors to cover the frequency range from 60 GHz to 600 GHz in 19 wide bands, in the focal plane of a 1.2 m aperture telescope cooled at 40 K, allowing for an accurate extraction of the CMB signal from polarized foreground emission. The projected CMB polarization survey sensitivity of this instrument, after foregrounds removal, is 1.7 {\mu}K$\cdot$arcmin. The design is robust enough to allow, if needed, a downscoped version of the instrument covering the 100 GHz to 600 GHz range with a 0.8 m aperture telescope cooled at 85 K, with a projected CMB polarization survey sensitivity of 3.2 {\mu}K$\cdot$arcmin., Comment: 43 pages

We discuss the effects on the CMB, CIB, and thermal SZ effect due to the peculiar motion of an observer with respect to the CMB rest frame, which induces boosting effects. We investigate the scientific perspectives opened by future CMB space missions, focussing on the CORE proposal. The improvements in sensitivity offered by a mission like CORE, together with its high resolution over a wide frequency range, will provide a more accurate estimate of the CMB dipole. The extension of boosting effects to polarization and cross-correlations will enable a more robust determination of purely velocity-driven effects that are not degenerate with the intrinsic CMB dipole, allowing us to achieve a S/N ratio of 13; this improves on the Planck detection and essentially equals that of an ideal cosmic-variance-limited experiment up to a multipole l of 2000. Precise inter-frequency calibration will offer the opportunity to constrain or even detect CMB spectral distortions, particularly from the cosmological reionization, because of the frequency dependence of the dipole spectrum, without resorting to precise absolute calibration. The expected improvement with respect to COBE-FIRAS in the recovery of distortion parameters (in principle, a factor of several hundred for an ideal experiment with the CORE configuration) ranges from a factor of several up to about 50, depending on the quality of foreground removal and relative calibration. Even for 1% accuracy in both foreground removal and relative calibration at an angular scale of 1 deg, we find that dipole analyses for a mission like CORE will be able to improve the recovery of the CIB spectrum amplitude by a factor of 17 in comparison with current results based on FIRAS. In addition to the scientific potential of a mission like CORE for these analyses, synergies with other planned and ongoing projects are also discussed., Comment: 61+5 pages, 17 figures, 25 tables, 8 sections, 5 appendices. In press on JCAP - Version 3 - Minor changes, affiliations fixed, references updated - version in line with corrected proofs

We demonstrate that, for the baseline design of the CORE satellite mission, the polarized foregrounds can be controlled at the level required to allow the detection of the primordial cosmic microwave background (CMB) $B$-mode polarization with the desired accuracy at both reionization and recombination scales, for tensor-to-scalar ratio values of ${r\gtrsim 5\times 10^{-3}}$. We consider detailed sky simulations based on state-of-the-art CMB observations that consist of CMB polarization with $\tau=0.055$ and tensor-to-scalar values ranging from $r=10^{-2}$ to $10^{-3}$, Galactic synchrotron, and thermal dust polarization with variable spectral indices over the sky, polarized anomalous microwave emission, polarized infrared and radio sources, and gravitational lensing effects. Using both parametric and blind approaches, we perform full component separation and likelihood analysis of the simulations, allowing us to quantify both uncertainties and biases on the reconstructed primordial $B$-modes. Under the assumption of perfect control of lensing effects, CORE would measure an unbiased estimate of $r=\left(5 \pm 0.4\right)\times 10^{-3}$ after foreground cleaning. In the presence of both gravitational lensing effects and astrophysical foregrounds, the significance of the detection is lowered, with CORE achieving a $4\sigma$-measurement of $r=5\times 10^{-3}$ after foreground cleaning and $60$% delensing. For lower tensor-to-scalar ratios ($r=10^{-3}$) the overall uncertainty on $r$ is dominated by foreground residuals, not by the 40% residual of lensing cosmic variance. Moreover, the residual contribution of unprocessed polarized point-sources can be the dominant foreground contamination to primordial B-modes at this $r$ level, even on relatively large angular scales, $\ell \sim 50$. Finally, we report two sources of potential bias for the detection of the primordial $B$-modes.[abridged], Comment: 87 pages, 32 figures, 4 tables, expanded abstract. Updated to match version accepted by JCAP

We examine the cosmological constraints that can be achieved with a galaxy cluster survey with the future CORE space mission. Using realistic simulations of the millimeter sky, produced with the latest version of the Planck Sky Model, we characterize the CORE cluster catalogues as a function of the main mission performance parameters. We pay particular attention to telescope size, key to improved angular resolution, and discuss the comparison and the complementarity of CORE with ambitious future ground-based CMB experiments that could be deployed in the next decade. A possible CORE mission concept with a 150 cm diameter primary mirror can detect of the order of 50,000 clusters through the thermal Sunyaev-Zeldovich effect (SZE). The total yield increases (decreases) by 25% when increasing (decreasing) the mirror diameter by 30 cm. The 150 cm telescope configuration will detect the most massive clusters ($>10^{14}\, M_\odot$) at redshift $z>1.5$ over the whole sky, although the exact number above this redshift is tied to the uncertain evolution of the cluster SZE flux-mass relation; assuming self-similar evolution, CORE will detect $\sim 500$ clusters at redshift $z>1.5$. This changes to 800 (200) when increasing (decreasing) the mirror size by 30 cm. CORE will be able to measure individual cluster halo masses through lensing of the cosmic microwave background anisotropies with a 1-$\sigma$ sensitivity of $4\times10^{14} M_\odot$, for a 120 cm aperture telescope, and $10^{14} M_\odot$ for a 180 cm one. [abridged], Comment: 35 pages, 15 figures, to be submitted to JCAP