Search

We characterise the co-evolution of radio-loud AGN and their galaxies by mapping the dependence of radio-loud AGN activity on stellar mass and star-formation rate (SFR) across cosmic time (out to $z \sim 1.5$). Deep LOFAR radio observations are combined with large galaxy samples to study the incidence of radio-loud AGN across the galaxy population; the AGN are further split into low-excitation radio galaxies (LERGs) and high-excitation radio galaxies (HERGs). We find that LERG activity occurs over a wide range of SFRs, whereas HERGs are typically found in galaxies with ongoing star formation. The LERGs are then split based on their SFRs relative to the main sequence, across redshift. Within quiescent galaxies, LERG activity shows a steep stellar mass dependence with the same normalisation across the past $\sim$ 10 Gyr; this indicates that hot gas fuels LERGs in quiescent galaxies across cosmic time. In massive galaxies ($\log_{10}(M/\rm{M_{\odot}}) \gtrsim 11$), the incidence of LERGs is roughly constant across the galaxy population, suggesting that LERGs in massive galaxies may be fuelled by hot gas regardless of the star-formation activity. At lower masses, however, LERG activity is significantly more enhanced (by a factor of up to 10) in star-forming galaxies compared to quiescent galaxies; this suggests that an additional fuelling mechanism, likely associated with cold gas, may fuel the LERGs in galaxies with higher SFRs. We find that HERGs typically accrete above 1 per cent of the Eddington-scaled accretion rate, and the LERGs typically accrete below this level., Comment: Accepted for publication in MNRAS, 18 pages, 10 figures

We present an analysis of 152,355 radio sources identified in the second data release of the LOFAR Two Metre Sky Survey (LoTSS-DR2) with Sloan Digital Sky Survey (SDSS) spectroscopic redshifts in the range 0.00 < z < 0.57. Using Monte Carlo simulations we determine the reliability of each source exhibiting an excess in radio luminosity relative to that predicted from their Ha emission, and, for a subset of 124,023 sources we combine this measurement with a full BPT analysis. Using these two independent diagnostics we determine the reliability of each source hosting a supermassive black hole of high or low Eddington-scaled accretion rate, and combine the measurements to determine the reliability of sources belonging to each of four physical classes of objects: star forming galaxies (SFGs), radio-quiet active galactic nuclei (RQAGN), and high- or low-excitation radio galaxies (HERGs or emission-line LERGs). The result is a catalogue which enables user-defined samples of radio sources with a reliability threshold suited to their science goal e.g. prioritising purity or completeness. Here we select high-confidence samples of radio sources (>90% reliability) to report: 38,588 radio-excess AGN in the LoTSS DR2 sample (362 HERGs, and 12,648 emission-line LERGs), together with 38,729 SFGs, and 18,726 RQAGN. We validate these results through comparison to literature using independent emission-line measurements, and to widely-adopted WISE photometric selection techniques. While our use of SDSS spectroscopy limits our current analysis to ~4 percent of the LoTSS-DR2 catalogue, our method is directly applicable to data from the forthcoming WEAVE-LOFAR survey which will obtain over a million spectra of 144 MHz selected sources., Comment: Accepted for publication in MNRAS

A fraction of galaxy clusters harbor diffuse radio sources known as radio halos. The currently adopted scenario for their formation is based on second-order Fermi re-acceleration of seed electrons that is driven by merger-driven turbulence in the intra-cluster medium. This mechanism is expected to be inefficient, which implies that a significant fraction of halos should have very steep ($\alpha < -1.5$) energy spectra. We start investigating the potential and current limitations of the combination of the two surveys conducted by LOFAR, LoTSS (144 MHz) and LoLSS (54 MHz), to probe the origin of radio halos. We follow up the 20 radio halos detected in the DR1 of LoTSS, which covers the HETDEX field, with the LoLSS survey, and we study their spectral properties between 54 and 144 MHz. After the removal of compact sources, 9 halos were excluded due to unreliable halo flux density measurements at 54 MHz. Our main finding is that 7 out of 11 ($\sim$ 64%) exhibit an ultra-steep spectrum ($\alpha < -1.5$), which is a key prediction of turbulent re-acceleration models. We also note a tentative trend for more massive systems to host flatter halos, although the currently poor statistics does not allow for a deeper analysis. Our sample suffers from low angular resolution at 54 MHz, which limits the accuracy of the compact-sources subtraction. Nevertheless, this study is the first step towards providing compelling evidence for the existence of a large fraction of radio halos with very steep spectrum, which is a fundamental prediction of turbulent re-acceleration models. In this regard, the forthcoming second data release of LoLSS, along with the integration of LOFAR international stations and the instrumental upgrade to LOFAR2.0, will improve both the statistics and the low-frequency angular resolution, allowing to conclusively determine the origin of radio halos in galaxy clusters., Comment: A&A Letter; 11 pages, 12 figures

Radio observations of a few cool-core galaxy clusters have revealed the presence of diffuse emission on cluster scales, similar to what was found in merging clusters in the form of radio halos. These sources might suggest that a minor merger, while not sufficiently energetic to disrupt the cool core, could still trigger particle acceleration in the intracluster medium on scales of hundreds of kpc. We observed with LOFAR at 144 MHz a sample of twelve cool-core galaxy clusters presenting some level of dynamical disturbances, according to X-ray data. We also performed a systematic search of cold fronts in these clusters, re-analysing archival Chandra data. The clusters PSZ1G139.61+24, A1068 (new detection), MS 1455.0+2232, and RX J1720.1+2638 present diffuse radio emission on a cluster scale. This emission is characterised by a double component: a central mini-halo confined by cold fronts and diffuse emission on larger scales, whose radio power at 144 MHz is comparable to that of radio halos detected in merging systems. The cold fronts in A1068 are a new detection. We also found a candidate plasma depletion layer in this cluster. No sloshing features are found in the other eight clusters. Two of them present a mini-halo, with diffuse radio emission confined to the cluster core. We also found a new candidate mini-halo. Whereas, for the remaining five clusters, we did not detect halo-like emission. For clusters without cluster-scale halos, we derived upper limits to the radio halo power. We found that cluster-scale diffuse radio emission is not present in all cool-core clusters when observed at a low frequency, but it is correlated to the presence of cold fronts. This morphology requires a specific configuration of the merger and so it puts some constraints on the turbulence, which deserves to be investigated in the future with theoretical works., Comment: 34 pages, 30 figures. Accepted for publication in A&A

The upcoming Square Kilometer Array (SKA) will set a new standard regarding data volume generated by an astronomical instrument, which is likely to challenge widely adopted data-analysis tools that scale inadequately with the data size. The aim of this study is to develop a new source detection and characterization method for massive radio astronomical datasets based on modern deep-learning object detection techniques. For this, we seek to identify the specific strengths and weaknesses of this type of approach when applied to astronomical data. We introduce YOLO-CIANNA, a highly customized deep-learning object detector designed specifically for astronomical datasets. In this paper, we present the method and describe all the elements introduced to address the specific challenges of radio astronomical images. We then demonstrate the capabilities of this method by applying it to simulated 2D continuum images from the SKA observatory Science Data Challenge 1 (SDC1) dataset. Using the SDC1 metric, we improve the challenge-winning score by +139% and the score of the only other post-challenge participation by +61%. Our catalog has a detection purity of 94% while detecting 40 to 60% more sources than previous top-score results, and exhibits strong characterization accuracy. The trained model can also be forced to reach 99% purity in post-process and still detect 10 to 30% more sources than the other top-score methods. It is also computationally efficient, with a peak prediction speed of 500 images of 512x512 pixels per second on a single GPU. YOLO-CIANNA achieves state-of-the-art detection and characterization results on the simulated SDC1 dataset and is expected to transfer well to observational data from SKA precursors., Comment: 42 pages. Custom layout edit of the corresponding in press A&A paper

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

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

Radio mini-halos are clouds of diffuse, low surface brightness synchrotron emission that surround the Brightest Cluster Galaxy (BCG) in massive cool-core galaxy clusters. In this paper, we use third generation calibration (3GC), also called direction-dependent (DD) calibration, and point source subtraction on MeerKAT extragalactic continuum data. We calibrate and image archival MeerKAT L-band observations of a sample of five galaxy clusters (ACO 1413, ACO 1795, ACO 3444, MACS J1115.8+0129, MACS J2140.2-2339). We use the CARACal pipeline for direction-independent (DI) calibration, DDFacet and killMS for 3GC, followed by visibility-plane point source subtraction to image the underlying mini-halo without bias from any embedded sources. Our 3GC process shows a drastic improvement in artefact removal, to the extent that the local noise around severely affected sources was halved and ultimately resulted in a 7\% improvement in global image noise. Thereafter, using these spectrally deconvolved Stokes I continuum images, we directly measure for four mini-halos the flux density, radio power, size and in-band integrated spectra. Further to that, we show the in-band spectral index maps of the mini-halo (with point sources). We present a new mini-halo detection hosted by MACS J2140.2-2339, having flux density $S_{\rm 1.28\,GHz} = 2.61 \pm 0.31$ mJy, average diameter 296 kpc and $\alpha^{\rm 1.5\,GHz}_{\rm 1\,GHz} = 1.21 \pm 0.36$. We also found a $\sim$100 kpc southern extension to the ACO 3444 mini-halo which was not detected in previous VLA L-band observations. Our description of MeerKAT wide-field, wide-band data reduction will be instructive for conducting further mini-halo science., Comment: 16 pages. 10 figures

The Square Kilometre Array Observatory (SKAO) will explore the radio sky to new depths in order to conduct transformational science. SKAO data products made available to astronomers will be correspondingly large and complex, requiring the application of advanced analysis techniques to extract key science findings. To this end, SKAO is conducting a series of Science Data Challenges, each designed to familiarise the scientific community with SKAO data and to drive the development of new analysis techniques. We present the results from Science Data Challenge 2 (SDC2), which invited participants to find and characterise 233245 neutral hydrogen (Hi) sources in a simulated data product representing a 2000~h SKA MID spectral line observation from redshifts 0.25 to 0.5. Through the generous support of eight international supercomputing facilities, participants were able to undertake the Challenge using dedicated computational resources. Alongside the main challenge, `reproducibility awards' were made in recognition of those pipelines which demonstrated Open Science best practice. The Challenge saw over 100 participants develop a range of new and existing techniques, with results that highlight the strengths of multidisciplinary and collaborative effort. The winning strategy -- which combined predictions from two independent machine learning techniques to yield a 20 percent improvement in overall performance -- underscores one of the main Challenge outcomes: that of method complementarity. It is likely that the combination of methods in a so-called ensemble approach will be key to exploiting very large astronomical datasets., Comment: Under review by MNRAS; 28 pages, 16 figures

A Faraday rotation measure (RM) catalogue, or RM Grid, is a valuable resource for the study of cosmic magnetism. Using the second data release (DR2) from the LOFAR Two-metre Sky Survey (LoTSS), we have produced a catalogue of 2461 extragalactic high-precision RM values across 5720 deg$^{2}$ of sky (corresponding to a polarized source areal number density of $\sim$0.43 deg$^{-2}$). The linear polarization and RM properties were derived using RM synthesis from the Stokes $Q$ and $U$ channel images at an angular resolution of 20'' across a frequency range of 120 to 168 MHz with a channel bandwidth of 97.6 kHz. The fraction of total intensity sources ($>1$ mJy beam$^{-1}$) found to be polarized was $\sim$0.2%. The median detection threshold was 0.6 mJy beam$^{-1}$ ($8\sigma_{QU}$), with a median RM uncertainty of 0.06 rad m$^{-2}$ (although a systematic uncertainty of up to 0.3 rad m$^{-2}$ is possible, after the ionosphere RM correction). The median degree of polarization of the detected sources is 1.8%, with a range of 0.05% to 31%. Comparisons with cm-wavelength RMs indicate minimal amounts of Faraday complexity in the LoTSS detections, making them ideal sources for RM Grid studies. Host galaxy identifications were obtained for 88% of the sources, along with redshifts for 79% (both photometric and spectroscopic), with the median redshift being 0.6. The focus of the current catalogue was on reliability rather than completeness, and we expect future versions of the LoTSS RM Grid to have a higher areal number density. In addition, 25 pulsars were identified, mainly through their high degrees of linear polarization., Comment: 20 pages, 23 figures, 3 tables. Accepted for publication in MNRAS. Go to https://lofar-mksp.org/data/ to download the catalogue and associated data products

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

Abell 3266 is a massive and complex merging galaxy cluster that exhibits significant substructure. We present new, highly sensitive radio continuum observations of Abell 3266 performed with the Australian Square Kilometre Array Pathfinder (0.8$-$1.1 GHz) and the Australia Telescope Compact Array (1.1$-$3.1 GHz). These deep observations provide new insights into recently-reported diffuse non-thermal phenomena associated with the intracluster medium, including a 'wrong-way' relic, a fossil plasma source, and an as-yet unclassified central diffuse ridge, which we reveal comprises the brightest part of a large-scale radio halo detected here for the first time. The 'wrong-way' relic is highly atypical of its kind: it exhibits many classical signatures of a shock-related radio relic, while at the same time exhibiting strong spectral steepening. While radio relics are generally consistent with a quasi-stationary shock scenario, the 'wrong-way' relic is not. We study the spectral properties of the fossil plasma source; it exhibits an ultra-steep and highly curved radio spectrum, indicating an extremely aged electron population. The larger-scale radio halo fills much of the cluster centre, and presents a strong connection between the thermal and non-thermal components of the intracluster medium, along with evidence of substructure. Whether the central diffuse ridge is simply a brighter component of the halo, or a mini-halo, remains an open question. Finally, we study the morphological and spectral properties of the multiple complex radio galaxies in this cluster in unprecedented detail, tracing their evolutionary history., Comment: Main article contains 26 pages, 15 figures, 6 tables. Supplementary material contains figures illustrating the improvement gain by applying third-generation calibration and imaging techniques to ATCA and ASKAP data. Manuscript published Open Access in MNRAS on August 1st

We report the discovery of a ring of low surface brightness radio emission around the Calvera pulsar, a high Galactic latitude, isolated neutron star, in the LOFAR Two-metre Sky Survey (LoTSS). It is centered at $\alpha=14\mathrm{h}11\mathrm{m}12.6\mathrm{s}$, $\delta=+79^\mathrm{o}23'15"$, has inner and outer radii of $14.2'$ and $28.4'$, and an integrated flux density at 144 MHz of $1.08\pm0.15$ Jy. The ring center is offset by $4.9'$ from the location of the Calvera pulsar. H$\alpha$ observations with the Isaac Newton Telescope show no coincident optical emission, but do show a small ($\sim20"$) optical structure internal to the ring. We consider three possible interpretations for the ring: that it is an H~II region, a supernova remnant (SNR), or an Odd Radio Circle (ORC). The positional coincidence of the ring, the pulsar, and an X-ray-emitting non-equilibrium ionisation plasma previously detected, lead us to prefer the SNR interpretation. If the source is indeed a SNR and its association with the Calvera pulsar is confirmed, then Calvera's SNR, or G118.4+37.0, will be one of few SNRs in the Galactic halo., Comment: Accepted in A&A

Context. Cosmic rays and magnetic fields are key ingredients in galaxy evolution, regulating both stellar feedback and star formation. Their properties can be studied with low-frequency radio continuum observations, free from thermal contamination. Aims. We define a sample of 76 nearby (< 30 Mpc) galaxies, with rich ancillary data in the radio continuum and infrared from the CHANG-ES and KINGFISH surveys, which will be observed with the LOFAR Two-metre Sky Survey (LoTSS) at 144 MHz. Methods. We present maps for 45 of them as part of the LoTSS data release 2 (LoTSS-DR2), where we measure integrated flux densities and study integrated and spatially resolved radio spectral indices. We investigate the radio-SFR relation, using star-formation rates (SFR) from total infrared and H $\alpha$ + 24-$\mu$m emission. Results. The radio-SFR relation at 144 MHz is clearly super-linear with $L_{144} \propto SFR^{1.4-1.5}$. The mean integrated radio spectral index between 144 and $\approx$1400 MHz is $\langle \alpha\rangle = -0.56 \pm 0.14$, in agreement with the injection spectral index for cosmic ray electrons (CRE). However, the radio spectral index maps show a variation of spectral indices with flatter spectra associated with star-forming regions and steeper spectra in galaxy outskirts and, in particular, in extra-planar regions. We found that galaxies with high star-formation rates (SFR) have steeper radio spectra; we find similar correlations with galaxy size, mass, and rotation speed. Conclusions. Galaxies that are larger and more massive are better electron calorimeters, meaning that the CRE lose a higher fraction of their energy within the galaxies. This explains the super-linear radio-SFR relation, with more massive, star-forming galaxies being radio bright. We propose a semi-calorimetric radio-SFR relation, which employs the galaxy mass as a proxy for the calorimetric efficiency., Comment: Accepted to Astronomy & Astrophysics. 16 pages, 8 figures, 3 tables + 45-page image atlas

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 study the properties and the origin of the radio emission in the most luminous early-type galaxies (ETGs) in the nearby Universe (MK<-25, recession velocity < 7,500 km/s) as seen by the 150 MHz Low-Frequency ARray (LOFAR) observations. LOFAR images are available for 188 of these giant ETGs (gETGs) and 146 (78%) of them are detected above a typical luminosity of ~10E21 W/Hz. They show a large spread in power, reaching up to ~10E26 W/Hz. We confirm a positive link between the stellar luminosity of gETGs and their median radio power, the detection rate, and the fraction of extended sources. About two-thirds (91) of the detected gETGs are unresolved, with sizes <4 kpc, confirming the prevalence of compact radio sources in local sources. Forty-six gETGs show extended emission on scales ranging from 4 to 340 kpc, at least 80% of which have a FRI class morphology. Based on the morphology and spectral index of the extended sources, ~30% of them might be remnant or restarted sources but further studies are needed to confirm this. Optical spectroscopy (available for 44 gETGs) indicates that for seven of them the nuclear gas is ionized by young stars suggesting a contribution to their radio emission from star forming regions. Their radio luminosities correspond to a star formation rate (SFR) in the range 0.1-8 Msun/yr and a median specific SFR of 0.8x10E-12 yr-1. The gas flowing toward the center of gETGs can accrete onto the supermassive black hole but also stall at larger radii and form new stars, an indication that feedback does not completely quench star formation. The most luminous gETGs (25 galaxies with MK < -25.8) are all detected at 150 MHz however they are not all currently turned on: at least four of them are remnant sources and at least one is likely powered by star formation., Comment: Pre-proofs version - Accepted for publication in A&A

High quality low-frequency radio surveys have the promise of advancing our understanding of many important topics in astrophysics, including the life cycle of active galactic nuclei (AGN), particle acceleration processes in jets, the history of star formation, and exoplanet magnetospheres. Currently leading low-frequency surveys reach an angular resolution of a few arcseconds. However, this resolution is not yet sufficient to study the more compact and distant sources in detail. Sub-arcsecond resolution is therefore the next milestone in advancing these fields. The biggest challenge at low radio frequencies is the ionosphere. If not adequately corrected for, ionospheric seeing blurs the images to arcsecond or even arcminute scales. Additionally, the required image size to map the degree-scale field of view of low-frequency radio telescopes at this resolution is far greater than what typical soft- and hardware is currently capable of handling. Here we present for the first time (to the best of our knowledge) widefield sub-arcsecond imaging at low radio frequencies. We derive ionospheric corrections in a few dozen individual directions and apply those during imaging efficiently using a recently developed imaging algorithm (arXiv:1407.1943, arXiv:1909.07226). We demonstrate our method by applying it to an eight hour observation of the International LOw Frequency ARray (LOFAR) Telescope (ILT) (arXiv:1305.3550). Doing so we have made a sensitive $7.4\ \mathrm{deg}^2$ $144\ \mathrm{MHz}$ map at a resolution of $0.3''$ reaching $25\ \mu\mathrm{Jy\ beam}^{-1}$ near the phase centre. The estimated $250,000$ core hours used to produce this image, fit comfortably in the budget of available computing facilities. This result will enable future mapping of the entire northern low-frequency sky at sub-arcsecond resolution., Comment: 19 pages, 4 figures; This preprint has not undergone peer review or any post-submission improvements or corrections. The Version of Record of this article is published in Nature Astronomy, and is available online at https://doi.org/10.1038/s41550-021-01573-z

The inner $\sim$200 pc region of the Galaxy contains a 4 million M$_{\odot}$ supermassive black hole (SMBH), significant quantities of molecular gas, and star formation and cosmic ray energy densities that are roughly two orders of magnitude higher than the corresponding levels in the Galactic disk. At a distance of only 8.2 kpc, the region presents astronomers with a unique opportunity to study a diverse range of energetic astrophysical phenomena, from stellar objects in extreme environments, to the SMBH and star-formation driven feedback processes that are known to influence the evolution of galaxies as a whole. We present a new survey of the Galactic center conducted with the South African MeerKAT radio telescope. Radio imaging offers a view that is unaffected by the large quantities of dust that obscure the region at other wavelengths, and a scene of striking complexity is revealed. We produce total intensity and spectral index mosaics of the region from 20 pointings (144 hours on-target in total), covering 6.5 square degrees with an angular resolution of 4$"$,at a central frequency of 1.28 GHz. Many new features are revealed for the first time due to a combination of MeerKAT's high sensitivity, exceptional $u,v$-plane coverage, and geographical vantage point. We highlight some initial survey results, including new supernova remnant candidates, many new non-thermal filament complexes, and enhanced views of the Radio Arc Bubble, Sgr A and Sgr B regions. This project is a SARAO public legacy survey, and the image products are made available with this article., Comment: 29 pages, 15 figures, 3 tables, accepted for publication in the The Astrophysical Journal. Replacement arXiv version with higher quality figures

The Shapley Supercluster ($\langle z \rangle\approx0.048$) contains several tens of gravitationally bound clusters and groups, making it it is an ideal subject for radio studies of cluster mergers. We used new high sensitivity radio observations to investigate the less energetic events of mass assembly in the Shapley Supercluster from supercluster down to galactic scales. We created total intensity images of the full region between A 3558 and A 3562, from $\sim 230$ to $\sim 1650$ MHz, using ASKAP, MeerKAT and the GMRT, with sensitivities ranging from $\sim 6$ to $\sim 100$ $\mu$Jy beam$^{-1}$. We performed a detailed morphological and spectral study of the extended emission features, complemented with ESO-VST optical imaging and X-ray data from XMM-Newton. We report the first GHz frequency detection of extremely low brightness intercluster diffuse emission on a $\sim 1$ Mpc scale connecting a cluster and a group, namely: A 3562 and the group SC 1329--313. It is morphologically similar to the X-ray emission in the region. We also found (1) a radio tail generated by ram pressure stripping in the galaxy SOS 61086 in SC 1329-313; (2) a head-tail radio galaxy, whose tail is broken and culminates in a misaligned bar; (3) ultrasteep diffuse emission at the centre of A 3558. Finally (4), we confirm the ultra-steep spectrum nature of the radio halo in A 3562. Our study strongly supports the scenario of a flyby of SC 1329-313 north of A 3562 into the supercluster core. [abridged...], Comment: 23 pages, 17 figures, A&A accepted

MeerKAT's large number of antennas, spanning 8 km with a densely packed 1 km core, create a powerful instrument for wide-area surveys, with high sensitivity over a wide range of angular scales. The MeerKAT Galaxy Cluster Legacy Survey (MGCLS) is a programme of long-track MeerKAT L-band (900-1670 MHz) observations of 115 galaxy clusters, observed for $\sim$6-10 hours each in full polarisation. The first legacy product data release (DR1), made available with this paper, includes the MeerKAT visibilities, basic image cubes at $\sim$8" resolution, and enhanced spectral and polarisation image cubes at $\sim$8" and 15" resolutions. Typical sensitivities for the full-resolution MGCLS image products are $\sim$3-5 {\mu}Jy/beam. The basic cubes are full-field and span 4 deg^2. The enhanced products consist of the inner 1.44 deg^2 field of view, corrected for the primary beam. The survey is fully sensitive to structures up to $\sim$10' scales and the wide bandwidth allows spectral and Faraday rotation mapping. HI mapping at 209 kHz resolution can be done at $0

Active Galactic Nuclei (AGN) inject large amounts of energy into their host galaxies and surrounding environment, shaping their properties and evolution. In particular, AGN jets inflate cosmic-ray lobes, which can rise buoyantly as light `bubbles' in the surrounding medium, displacing and heating the encountered thermal gas and thus halting its spontaneous cooling. These bubbles have been identified in a wide range of systems. However, due to the short synchrotron lifetime of electrons, the most advanced phases of their evolution have remained observationally unconstrained, preventing us to fully understand their coupling with the external medium, and thus AGN feedback. Simple subsonic hydrodynamic models predict that the pressure gradients, naturally present around the buoyantly rising bubbles, transform them into toroidal structures, resembling mushroom clouds in a stratified atmosphere. The way and timescales on which these tori will eventually disrupt depend on various factors including magnetic fields and plasma viscosity. Here we report LOFAR observations below 200 MHz, sensitive to the oldest radio-emitting particles, showing the late evolution of multiple generations of cosmic-ray AGN bubbles in a galaxy group with unprecedented level of detail. The bubbles' buoyancy power can efficiently offset the radiative cooling of the intragroup medium. However, the bubbles have still not thoroughly mixed with the thermal gas, after hundreds of million years, likely under the action of magnetic fields., Comment: 34 pages, 9 figures, 5 tables, published in Nature Astronomy

Diffuse radio emission at the centre of galaxy clusters has been observed both in merging clusters on scales of Mpc, called giant radio haloes, and in relaxed systems with a cool-core on smaller scales, named mini haloes. Giant radio haloes and mini haloes are thought to be distinct classes of sources. However, recent observations have revealed the presence of diffuse radio emission on Mpc scales in clusters that do not show strong dynamical activity. RX J1720.1+2638 is a cool-core cluster, presenting both a bright central mini halo and a fainter diffuse, steep-spectrum emission extending beyond the cluster core that resembles giant radio halo emission. In this paper, we present new observations performed with the LOFAR Low Band Antennas (LBA) at 54 MHz. These observations, combined with data at higher frequencies, allow us to constrain the spectral properties of the radio emission. The large-scale emission presents an ultra-steep spectrum with $\alpha_{54}^{144}\sim3.2$. The radio emission inside and outside the cluster core have strictly different properties, as there is a net change in spectral index and they follow different radio-X-ray surface brightness correlations. We argue that the large-scale diffuse emission is generated by particles re-acceleration after a minor merger. While for the central mini halo we suggest that it could be generated by secondary electrons and positrons from hadronic interactions of relativistic nuclei with the dense cool-core gas, as an alternative to re-acceleration models., Comment: 13 pages, 10 figures; accepted for publication in MNRAS

MIGHTEE is a galaxy evolution survey using simultaneous radio continuum, spectro-polarimetry, and spectral line observations from the South African MeerKAT telescope. When complete, the survey will image $\sim$20 deg$^{2}$ over the COSMOS, E-CDFS, ELAIS-S1, and XMM-LSS extragalactic deep fields with a central frequency of 1284 MHz. These were selected based on the extensive multiwavelength datasets from numerous existing and forthcoming observational campaigns. Here we describe and validate the data processing strategy for the total intensity continuum aspect of MIGHTEE, using a single deep pointing in COSMOS (1.6 deg$^{2}$) and a three-pointing mosaic in XMM-LSS (3.5 deg$^{2}$). The processing includes the correction of direction-dependent effects, and results in thermal noise levels below 2~$\mathrm{\mu}$Jy beam$^{-1}$ in both fields, limited in the central regions by classical confusion at $\sim$8$''$ angular resolution, and meeting the survey specifications. We also produce images at $\sim$5$''$ resolution that are $\sim$3 times shallower. The resulting image products form the basis of the Early Science continuum data release for MIGHTEE. From these images we extract catalogues containing 9,896 and 20,274 radio components in COSMOS and XMM-LSS respectively. We also process a close-packed mosaic of 14 additional pointings in COSMOS and use these in conjunction with the Early Science pointing to investigate methods for primary beam correction of broadband radio images, an analysis that is of relevance to all full-band MeerKAT continuum observations, and wide field interferometric imaging in general. A public release of the MIGHTEE Early Science continuum data products accompanies this article., Comment: 21 pages, 13 figures, 3 tables, to be published in MNRAS

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

Context. Radio halos are megaparsec-scale diffuse radio sources{ mostly} located at the centres of merging galaxy clusters. The common mechanism invoked to explain their origin is the re-acceleration of relativistic particles caused by large-scale turbulence. Aims. Current re-acceleration models predict that a significant number of halos at high redshift should be characterised by very steep spectra ($\alpha<-1.5$) because of increasing inverse Compton energy losses. In this paper, we investigate the spectral index properties of a sample of nine clusters selected from the second Planck Sunyaev-Zel'dovich catalogue showing diffuse radio emission with the Low Frequency Array (LOFAR) in the 120-168 MHz band. This is the first time that radio halos discovered at low frequencies are followed up at higher frequencies. Methods. We analysed upgraded Giant Metrewave Radio Telescope (uGMRT) observations in Bands 3 and 4, that is, 250-500 and 550-900 MHz respectively. These observations were combined with existing LOFAR data to obtain information on the spectral properties of the diffuse radio emission. Results. We find diffuse radio emission in the uGMRT observations for five of the nine high-$z$ radio halos previously discovered with LOFAR. For those, we measure spectral indices in the range of $-1$ to $-1.4$. For the uGMRT non-detections, we estimated that the halos should have a spectral index steeper than $-1.5$. We also confirm the presence of one candidate relic. Conclusions. Despite the small number of clusters, we find evidence that about half of the massive and merging clusters at high redshift host radio halos with a very steep spectrum. This is in line with theoretical predictions, although larger statistical samples are necessary to test models., Comment: 23 pages, 22 figures, 4 tables; accepted for publication in A&A

[abridged] The International LOFAR Telescope is an interferometer with stations spread across Europe. With baselines of up to ~2,000 km, LOFAR has the unique capability of achieving sub-arcsecond resolution at frequencies below 200 MHz, although this is technically and logistically challenging. Here we present a calibration strategy that builds on previous high-resolution work with LOFAR. We give an overview of the calibration strategy and discuss the special challenges inherent to enacting high-resolution imaging with LOFAR, and describe the pipeline, which is publicly available, in detail. We demonstrate the calibration strategy by using the pipeline on P205+55, a typical LOFAR Two-metre Sky Survey (LoTSS) pointing. We perform in-field delay calibration, solution referencing to other calibrators, self-calibration, and imaging of example directions of interest in the field. For this specific field and these ionospheric conditions, dispersive delay solutions can be transferred between calibrators up to ~1.5 degrees away, while phase solution transferral works well over 1 degree. We demonstrate a check of the astrometry and flux density scale. Imaging in 17 directions, the restoring beam is typically 0.3" x 0.2" although this varies slightly over the entire 5 square degree field of view. We achieve ~80 to 300 $\mu$Jy/bm image rms noise, which is dependent on the distance from the phase centre; typical values are ~90 $\mu$Jy/bm for the 8 hour observation with 48 MHz of bandwidth. Seventy percent of processed sources are detected, and from this we estimate that we should be able to image ~900 sources per LoTSS pointing. This equates to ~3 million sources in the northern sky, which LoTSS will entirely cover in the next several years. Future optimisation of the calibration strategy for efficient post-processing of LoTSS at high resolution (LoTSS-HR) makes this estimate a lower limit., Comment: Accepted to a special issue of A&A on sub-arcsecond imaging with LOFAR. 24 pages, 16 figures

The prominent radio source Hercules A features complex structures in its radio lobes. Although it is one of the most comprehensively studied sources in the radio sky, the origin of the ring structures in the Hercules A radio lobes remains an open question. We present the first sub-arcsecond angular resolution images at low frequencies (<300 MHz) of Hercules A, made with the International LOFAR Telescope. With the addition of data from the Karl G. Jansky Very Large Array, we mapped the structure of the lobes from 144 MHz to 7 GHz. We explore the origin of the rings within the lobes of Hercules A, and test whether their properties are best described by a shock model, where shock waves are produced by the jet propagating in the radio lobe, or by an inner-lobe model, where the rings are formed by decelerated jetted plasma. From spectral index mapping our large frequency coverage reveals that the curvature of the different ring spectra increases with distance away from the central active galactic nucleus. We demonstrate that the spectral shape of the rings is consistent with synchrotron aging, which speaks in favor of an inner-lobe model where the rings are formed from the deposition of material from past periods of intermittent core activity., Comment: Accepted for a special issue of A&A on sub-arcsecond imaging with LOFAR

In this paper we present a large sample of jellyfish galaxies in low redshift clusters (z<0.05), identified through 120-168 MHz radio continuum from the LOFAR Two-metre Sky Survey (LoTSS). From a parent sample of 29 X-ray-detected SDSS galaxy clusters and their spectroscopic members, we visually identify 95 star-forming, LoTSS jellyfish galaxies with 144 MHz radio tails. Star formation rates (SFRs) and stellar masses are obtained for all galaxies from SED fits. For each jellyfish galaxy we determine the tail orientation with respect to the cluster centre and quantify the prominence of the radio tails with the 144 MHz shape asymmetry. After carefully accounting for redshift-dependent selection effects, we find that the frequency of jellyfish galaxies is relatively constant from cluster to cluster. LoTSS jellyfish galaxies are preferentially found at small clustercentric radius and large velocity offsets within their host clusters and have radio tails that are oriented away from the cluster centre. These galaxies also show enhanced star formation, relative to both 'normal' cluster galaxies and isolated field galaxies, but generally fall within the scatter of the L144MHz - SFR relation. The properties of the LoTSS jellyfish galaxies identified in this work are fully consistent with expectations from ram pressure stripping. This large sample of jellyfish galaxies will be valuable for further constraining ram pressure stripping and star formation quenching in nearby galaxy clusters. We show that LOFAR is a powerful instrument for identifying ram pressure stripped galaxies across extremely wide fields. Moving forward we will push the search for jellyfish galaxies beyond this initial cluster sample, including a comprehensive survey of the galaxy group regime., Comment: 15 pages plus appendix, 10 figures, accepted for publication in A&A

We report extremely faint 144 MHz radio emission from two gravitationally lensed quasars, SDSS J1004+4112 (z = 1.730) and SDSS J2222+2745 (z = 2.803), using the LOFAR Two Metre Sky Survey (LoTSS) data release 2. After correcting for the lensing magnifications, the two objects have intrinsic flux-densities of 13+/-2 and 58+/-6 uJy, respectively, corresponding to 144 MHz rest-frame luminosities of 10^(23.2+/-0.2) and 10^(24.42+/-0.05) W / Hz, respectively. In the case of SDSS J1004+4112, the intrinsic flux density is close to the confusion limit of LoTSS, making this radio source the faintest to be detected thus far at low frequencies, and the lowest luminosity known at z > 0.65. Under the assumption that all of the radio emission is due to star-formation processes, the quasar host galaxies are predicted to have star-formation rates of 5.5^(+1.8)_(-1.4) and 73^(+34)_(-22) M / yr, respectively. Further multi-wavelength observations at higher angular resolution will be needed to determine if any of the detected radio emission is due to weak jets associated with the quasars., Comment: 5 pages, 1 figure, accepted for publication in MNRAS Letters

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

Thermal gas in the center of galaxy clusters can show substantial motions that generate surface-brightness and temperature discontinuities known as cold fronts. The motions may be triggered by minor or off-axis mergers that preserve the cool core of the system. The dynamics of the thermal gas can also generate radio emission from the intra-cluster medium (ICM) and impact the evolution of clusters' radio sources. We aim to study the central region of Abell 1775, a system in an ambiguous dynamical state at $z=0.072$ which is known to host an extended head-tail radio galaxy, with the goal of investigating the connection between thermal and nonthermal components in its center. We made use of a deep (100 ks) Chandra observation accompanied by LOFAR 144 MHz, GMRT 235 MHz and 610 MHz, and VLA 1.4 GHz radio data. We find a spiral-like pattern in the X-ray surface brightness that is mirrored in the temperature and pseudo-entropy maps. Additionally, we characterize an arc-shaped cold front in the ICM. We interpret these features in the context of a slingshot gas tail scenario. The structure of the head-tail radio galaxy "breaks" at the position of the cold front, showing an extension that is detected only at low frequencies, likely due to its steep and curved spectrum. We speculate that particle reacceleration is occurring in the outer region of this tail, which in total covers a projected size of $\sim800$ kpc. We also report the discovery of revived fossil plasma with ultra-steep spectrum radio emission in the cluster core together with a central diffuse radio source that is bounded by the arc-shaped cold front. The results reported in this work demonstrate the interplay between thermal and nonthermal components in the cluster center and the presence of ongoing particle reacceleration in the ICM on different scales., Comment: 17 pages, 17 figures, 2 tables (including appendix); Updated to match the accepted version in A&A

We report the discovery of two new giant radio galaxies (GRGs) using the MeerKAT International GHz Tiered Extragalactic Exploration (MIGHTEE) survey. Both GRGs were found within a 1 deg^2 region inside the COSMOS field. They have redshifts of z=0.1656 and z=0.3363 and physical sizes of 2.4Mpc and 2.0Mpc, respectively. Only the cores of these GRGs were clearly visible in previous high resolution VLA observations, since the diffuse emission of the lobes was resolved out. However, the excellent sensitivity and uv coverage of the new MeerKAT telescope allowed this diffuse emission to be detected. The GRGs occupy a unpopulated region of radio power - size parameter space. Based on a recent estimate of the GRG number density, the probability of finding two or more GRGs with such large sizes at z<0.4 in a ~1deg^2 field is only 2.7x10^-6, assuming Poisson statistics. This supports the hypothesis that the prevalence of GRGs has been significantly underestimated in the past due to limited sensitivity to low surface brightness emission. The two GRGs presented here may be the first of a new population to be revealed through surveys like MIGHTEE which provide exquisite sensitivity to diffuse, extended emission., Comment: 15 pages; 8 figures, accepted for publication in MNRAS on 2020 Dec 9

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

In recent years, the level of the extragalactic radio background has become a point of considerable interest, with some lines of argument pointing to an entirely new cosmological synchrotron background. The contribution of the known discrete source population to the sky temperature is key to this discussion. Because of the steep spectral index of the excess over the cosmic microwave background, it is best studied at low frequencies where the signal is strongest. The Low-Frequency Array (LOFAR) wide and deep sky surveys give us the best constraints yet on the contribution of discrete extragalactic sources at 144 MHz, and in particular allow us to include contributions from diffuse, low-surface-brightness emission that could not be fully accounted for in previous work. We show that, even with these new data, known sources can still only account for around a quarter of the estimated extragalactic sky temperature at LOFAR frequencies., Comment: 8 pages, 4 figures. This paper is part of the 1st data release of the LoTSS Deep Fields

We present deep polarimetric observations of the European Large Area ISO Survey-North 1 (ELAIS-N1) field using the Low Frequency Array (LOFAR) at 114.9-177.4 MHz. The ELAIS-N1 field is part of the LOFAR Two-metre Sky Survey deep fields data release I. For six eight-hour observing epochs, we align the polarization angles and stack the 20"-resolution Stokes $Q$, $U$-parameter data cubes. This produces a 16 deg$^2$ image with 1$\sigma_{\rm QU}$ sensitivity of 26 $\mu$Jy/beam in the central area. In this paper, we demonstrate the feasibility of the stacking technique, and we generate a catalog of polarized sources in ELAIS-N1 and their associated Faraday rotation measures (RMs). While in a single-epoch observation we detect three polarized sources, this number increases by a factor of about three when we consider the stacked data, with a total of ten sources. This yields a surface density of polarized sources of one per 1.6 deg$^2$. The Stokes $I$ images of three of the ten detected polarized sources have morphologies resembling those of FR I radio galaxies. This represents a greater fraction of this type of source than previously found, which suggests that more sensitive observations may help with their detection., Comment: This paper is part of the 1st data release of the LoTSS Deep Fields. Accepted for publication in A&A. 14 pages, 9 figures, 7 tables

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

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 effects of feedback from high luminosity radio-loud AGN have been extensively discussed in the literature, but feedback from low-luminosity radio-loud AGN is less well understood. The advent of high sensitivity, high angular resolution, large field of view telescopes such as LOFAR is now allowing wide-area studies of such faint sources for the first time. Using the first data release of the LOFAR Two Metre Sky Survey (LoTSS) we report on our discovery of a population of 195 radio galaxies with 150 MHz luminosities between $3\times10^{22}$ and $1.5\times10^{25}\text{ W Hz}^{-1}$ and total radio emission no larger than 80 kpc. These objects, which we term galaxy-scale jets (GSJ), are small enough to be directly influencing the evolution of the host on galaxy scales. We report upon the typical host properties of our sample, finding that 9 per cent are hosted by spirals with the remainder being hosted by elliptical galaxies. Two of the spiral-hosted GSJ are highly unusual with low radio luminosities and FRII-like morphology. The host properties of our GSJ show that they are ordinary AGN observed at a stage in their life shortly after the radio emission has expanded beyond the central regions of the host. Based on our estimates, we find that about half of our GSJ have internal radio lobe energy within an order of magnitude of the ISM energy so that, even ignoring any possible shocks, GSJ are energetically capable of affecting the evolution of the host. The current sample of GSJ will grow in size with future releases of LoTSS and can also form the basis for further studies of feedback from low-luminosity radio sources., Comment: 23 pages, accepted by MNRAS

Various studies have laid claim to finding an alignment of the polarization vectors or radio jets of active galactic nuclei (AGN) over large distances, but these results have proven controversial and so far, there is no clear explanation for this observed alignment. To investigate this case further, we tested the hypothesis that the position angles of radio galaxies are randomly oriented in the sky by using data from the Low-Frequency Array (LOFAR) Two-metre Sky Survey (LoTSS). A sample of 7,555 double-lobed radio galaxies was extracted from the list of 318,520 radio sources in the first data release of LoTSS at 150 MHz. We performed statistical tests for uniformity of the two-dimensional (2D) orientations for the complete 7,555 source sample. We also tested the orientation uniformity in three dimensions (3D) for the 4,212 source sub-sample with photometric or spectroscopic redshifts. Our sample shows a significant deviation from uniformity (p-value < $10^{-5}$) in the 2D analysis at angular scales of about four degrees, mainly caused by sources with the largest flux densities. No significant alignment was found in the 3D analysis. Although the 3D analysis has access to fewer sources and suffers from uncertainties in the photometric redshift, the lack of alignment in 3D points towards the cause of the observed effect being unknown systematics or biases that predominantly affect the brightest sources, although this has yet to be demonstrated irrefutably and should be the subject of subsequent studies., Comment: 10 pages, 13 figures. Accepted for publication in Astronomy & Astrophysics

We present LOFAR observations at 150 MHz of the borderline FRI/FRII giant radio galaxy NGC 6251. This paper presents the most sensitive and highest-resolution images of NGC 6251 at these frequencies to date, revealing for the first time a low-surface-brightness extension to the northern lobe, and a possible backflow associated with the southern lobe. The integrated spectra of components of NGC 6251 are consistent with previous measurements at higher frequencies, similar to results from other LOFAR studies of nearby radio galaxies. We find the outer structures of NGC 6251 to be either at equipartition or slightly electron dominated, similar to those of FRII sources rather than FRIs; but this conclusion remains tentative because of uncertainties associated with the geometry and the extrapolation of X-ray measurements to determine the external pressure distribution on the scale of the outer lobes. We place lower limits on the ages of the extension of the northern lobe and the backflow of the southern lobe of $t \gtrsim 250$ Myr and $t \gtrsim 210$ Myr respectively. We present the first detection of polarisation at 150 MHz in NGC 6251. Taking advantage of the high Faraday resolution of LOFAR, we place an upper limit on the magnetic field in the group of $B < 0.2 (\Lambda_B / 10 {\rm kpc})^{-0.5} \mu$G for a coherence scale of $\Lambda_B < 60 {\rm kpc}$ and $B < 13 \mu$G for $\Lambda_B = 240$ kpc., Comment: 19 pages, 17 figures; accepted by MNRAS

Radio galaxies are known to go through cycles of activity, where phases of apparent quiescence can be followed by repeated activity of the central supermassive black hole. A better understanding of this cycle is crucial for ascertaining the energetic impact that the jets have on the host galaxy, but little is known about it. We used deep LOFAR images at 150 MHz of the Lockman Hole extragalactic field to select a sample of 158 radio sources with sizes $> 60^{\prime\prime}$ in different phases of their jet life cycle. Using a variety of criteria (e.g. core prominence combined with low-surface brightness of the extended emission and steep spectrum of the central region) we selected a subsample of candidate restarted radio galaxies representing between 13% and 15% of the 158 sources of the main sample. We compare their properties to the rest of the sample, which consists of remnant candidates and active radio galaxies. Optical identifications and characterisations of the host galaxies indicate similar properties for candidate restarted, remnant, and active radio galaxies, suggesting that they all come from the same parent population. The fraction of restarted radio galaxies is slightly higher with respect to remnants, suggesting that the restarted phase can often follow after a relatively short remnant phase (the duration of the remnant phase being a few times 10$^{7}$ years). This confirms that the remnant and restarted phases are integral parts of the life cycle of massive elliptical galaxies. A preliminary investigation does not suggest a strong dependence of this cycle on the environment surrounding any given galaxy., Comment: 16 pages, 9 figures, 4 tables, A&A accepted, companion paper to Shabala et al., also in today arXiv posting

Measuring the properties of extragalactic magnetic fields through the effect of Faraday rotation provides a means to understand the origin and evolution of cosmic magnetism. Here we use data from the LOFAR Two-Metre Sky Survey (LoTSS) to calculate the Faraday rotation measure (RM) of close pairs of extragalactic radio sources. By considering the RM difference ($\Delta$RM) between physical pairs (e.g. double-lobed radio galaxies) and non-physical pairs (i.e. close projected sources on the sky), we statistically isolate the contribution of extragalactic magnetic fields to $\Delta$RM along the line of sight between non-physical pairs. From our analysis, we find no significant difference between the $\Delta$RM distributions of the physical and non-physical pairs, limiting the excess Faraday rotation contribution to $< 1.9$ rad/m$^2$ ($\sim$$95\%$ confidence). We use this limit with a simple model of an inhomogeneous universe to place an upper limit of 4 nG on the cosmological co-moving magnetic field strength on Mpc scales. We also compare the RM data with a more realistic suite of cosmological MHD simulations, that explore different magnetogenesis scenarios. Both magnetization of the large scale structure by astrophysical processes such as galactic and AGN outflows, and simple primordial scenarios with seed magnetic field strengths $< 0.5$ nG cannot be rejected by the current data; while stronger primordial fields or models with dynamo amplification in filaments are disfavoured., Comment: 13 pages, 8 figures, 2 tables. Accepted to MNRAS. Full version of Table 2 available at https://www.dropbox.com/s/axqrl8sbmj1f2c1/RMpairs_Table_2.dat?dl=0

We present the confusion-limited 1.28 GHz MeerKAT DEEP2 image covering one $\approx 68'$ FWHM primary beam area with $7.6''$ FWHM resolution and $0.55 \pm 0.01$ $\mu$Jy/beam rms noise. Its J2000 center position $\alpha=04^h 13^m 26.4^s$, $\delta=-80^\circ 00' 00''$ was selected to minimize artifacts caused by bright sources. We introduce the new 64-element MeerKAT array and describe commissioning observations to measure the primary beam attenuation pattern, estimate telescope pointing errors, and pinpoint $(u,v)$ coordinate errors caused by offsets in frequency or time. We constructed a 1.4 GHz differential source count by combining a power-law count fit to the DEEP2 confusion $P(D)$ distribution from $0.25$ to $10$ $\mu$Jy with counts of individual DEEP2 sources between $10$ $\mu$Jy and $2.5$ mJy. Most sources fainter than $S \sim 100$ $\mu$Jy are distant star-forming galaxies obeying the FIR/radio correlation, and sources stronger than $0.25$ $\mu$Jy account for $\sim93\%$ of the radio background produced by star-forming galaxies. For the first time, the DEEP2 source count has reached the depth needed to reveal the majority of the star formation history of the universe. A pure luminosity evolution of the 1.4 GHz local luminosity function consistent with the Madau & Dickinson (2014) model for the evolution of star-forming galaxies based on UV and infrared data underpredicts our 1.4 GHz source count in the range $-5 \lesssim \log[S(\mathrm{Jy})] \lesssim -4$., Comment: 20 pages, 18 figures. Accepted for publication in ApJ

It is well established that particle acceleration by shocks and turbulence in the intra-cluster medium can produce cluster-scale synchrotron emitting sources. However, the detailed physics of these particle acceleration processes is still not well understood. One of the main open questions is the role of fossil relativistic electrons that have been deposited in the intra-cluster medium by radio galaxies. These synchrotron-emitting electrons are very difficult to study, as their radiative life time is only tens of Myrs at GHz frequencies, and are therefore a relatively unexplored population. Despite the typical steep radio spectrum due to synchrotron losses, these fossil electrons are barely visible even at radio frequencies well below a GHz. However, when a pocket of fossil radio plasma is compressed, it boosts the visibility at sub-GHz frequencies, creating so-called radio phoenices. This compression can be the result of bulk motion and shocks in the ICM due to merger activity. In this paper, we demonstrate the discovery potential of low frequency radio sky surveys to find and study revived fossil plasma sources in galaxy clusters. We used the 150~MHz TGSS and 1.4 GHz NVSS sky surveys to identify candidate radio phoenices. A subset of three candidates were studied in detail using deep multi-band radio observations (LOFAR and GMRT), X-ray (\textit{Chandra} or \textit{XMM-Newton}) and archival optical observations. Two of the three sources are new discoveries. Using these observations, we identified common observational properties (radio morphology, ultra-steep spectrum, X-ray luminosity, dynamical state) that will enable us to identify this class of sources more easily, and help to understand the physical origin of these sources., Comment: 12 pages, 5 figures; Accepted for publication in Astronomy and Astrophysics

The Galactic Centre contains a supermassive black hole with a mass of 4 million suns within an environment that differs markedly from that of the Galactic disk. While the black hole is essentially quiescent in the broader context of active galactic nuclei, X-ray observations have provided evidence for energetic outbursts from its surroundings. Also, while the levels of star formation in the Galactic Centre have been approximately constant over the last few hundred Myr, there is evidence of elevated short-duration bursts, strongly influenced by interaction of the black hole with the enhanced gas density present within the ring-like Central Molecular Zone at Galactic longitude |l| < 0.7 degrees and latitude |b| < 0.2 degrees. The inner 200 pc region is characterized by large amounts of warm molecular gas, a high cosmic ray ionization rate, unusual gas chemistry, enhanced synchrotron emission, and a multitude of radio-emitting magnetised filaments, the origin of which has not been established. Here we report radio imaging that reveals bipolar bubbles spanning 1 degree x 3 degrees (140 parsecs x 430 parsecs), extending above and below the Galactic plane and apparently associated with the Galactic Centre. The structure is edge-brightened and bounded, with symmetry implying creation by an energetic event in the Galactic Centre. We estimate the age of the bubbles to be a few million years, with a total energy of 7 x 10^52 ergs. We postulate that the progenitor event was a major contributor to the increased cosmic-ray density in the Galactic Centre, and is in turn the principal source of the relativistic particles required to power the synchrotron emission of the radio filaments within and in the vicinity of the bubble cavities., Comment: 25 pages, 3 figures, 3 tables, authors' version of a Letter published in Nature on 11 September 2019

In this Letter, we report the discovery of a radio halo in the high-redshift galaxy cluster PSZ2 G099.86+58.45 ($z=0.616$) with the LOw Frequency ARray (LOFAR) at 120-168 MHz. This is one of the most distant radio halos discovered so far. The diffuse emission extends over $\sim$ 1 Mpc and has a morphology similar to that of the X-ray emission as revealed by XMM-Newton data. The halo is very faint at higher frequencies and is barely detected by follow-up 1-2 GHz Karl G.~Jansky Very Large Array (JVLA) observations, which enable us to constrain the radio spectral index to be $\alpha\leq 1.5-1.6$, i.e.; with properties between canonical and ultra-steep spectrum radio halos. Radio halos are currently explained as synchrotron radiation from relativistic electrons that are re-accelerated in the intra-cluster medium (ICM) by turbulence driven by energetic mergers. We show that in such a framework radio halos are expected to be relatively common at $\sim150$ MHz ($\sim30-60\%$) in clusters with mass and redshift similar to PSZ2 G099.86+58.45; however, at least 2/3 of these radio halos should have steep spectrum and thus be very faint above $\sim 1$ GHz frequencies. Furthermore, since the luminosity of radio halos at high redshift depends strongly on the magnetic field strength in the hosting clusters, future LOFAR observations will also provide vital information on the origin and amplification of magnetic fields in galaxy clusters., Comment: 6 pages, 3 figures, accepted (2019, July 17) for publication in ApJ Letter

We present new 144-MHz LOFAR observations of the prototypical `X-shaped' radio galaxy NGC 326, which show that the formerly known wings of the radio lobes extend smoothly into a large-scale, complex radio structure. We argue that this structure is most likely the result of hydrodynamical effects in an ongoing group or cluster merger, for which pre-existing X-ray and optical data provide independent evidence. The large-scale radio structure is hard to explain purely in terms of jet reorientation due to the merger of binary black holes, a previously proposed explanation for the inner structure of NGC 326. For this reason, we suggest that the simplest model is one in which the merger-related hydrodynamical processes account for all the source structure, though we do not rule out the possibility that a black hole merger has occurred. Inference of the black hole-black hole merger rate from observations of X-shaped sources should be carried out with caution in the absence of deep, sensitive low-frequency observations. Some X-shaped sources may be signposts of cluster merger activity, and it would be useful to investigate the environments of these objects more generally., Comment: 7 pages, 2 colour figures. Accepted by MNRAS