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

Low-frequency radio observations show an increasing number of radio galaxies located in galaxy clusters that display peculiar morphologies and spectral profiles. This is the result of the dynamical interaction of the galaxy with the surrounding medium. Studying this phenomenon is key to understanding the evolution of low-energy relativistic particles in the intracluster medium. We present a multi-frequency study of the three head-tail (HT) radio galaxies and the radio halo in the galaxy cluster ZwCl0634.1+4747. We make use of observations at four frequencies performed with LOFAR LBA (53 MHz), HBA (144 MHz), GMRT (323 MHz) and VLA (1518 MHz) data. The use of extremely low radio frequency observations, such as LOFAR at 53 and 144 MHz, allowed us to detect the extension of the tails up to a distance of ~ 1 Mpc. We extracted spectral profiles along the tails in order to identify possible departures from a pure ageing model, such as the Jaffe-Perola (JP) model, which only involves synchrotron and inverse-Compton losses. We found clear evidence of departures from this simple ageing model, such as surface brightness enhancement and spectral flattening along all of the tails. This can be interpreted as the consequence of particle re-acceleration along the tails. Possible explanations for this behaviour include the interaction between a shock and the radio tails or a turbulence-driven re-acceleration mechanism. We show that the latter scenario is able to reproduce the characteristic features that we observed in our profiles.

We report on a radiopolarimetric observation of the Saturn-Jupiter Great Conjunction of 2020 using the MeerKAT L-band system, initially carried out for science verification purposes, which yielded a serendipitous discovery of a pulsar. The radiation belts of Jupiter are very bright and time variable: coupled with the sensitivity of MeerKAT, this necessitated development of dynamic imaging techniques, reported on in this work. We present a deep radio "movie" revealing Jupiter's rotating magnetosphere, a radio detection of Callisto, and numerous background radio galaxies. We also detect a bright radio transient in close vicinity to Saturn, lasting approximately 45 minutes. Follow-up deep imaging observations confirmed this as a faint compact variable radio source, and yielded detections of pulsed emission by the commensal MeerTRAP search engine, establishing the object's nature as a radio emitting neutron star, designated PSR J2009-2026. A further observation combining deep imaging with the PTUSE pulsar backend measured detailed dynamic spectra for the object. While qualitatively consistent with scintillation, the magnitude of the magnification events and the characteristic timescales are odd. We are tentatively designating this object a pulsar with anomalous refraction recurring on odd timescales (PARROT). As part of this investigation, we present a pipeline for detection of variable sources in imaging data, with dynamic spectra and lightcurves as the products, and compare dynamic spectra obtained from visibility data with those yielded by PTUSE. We discuss MeerKAT's capabilities and prospects for detecting more of such transients and variables., Comment: 21 pages, 17 figures, submitted to MNRAS

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

We present the first wide area (2.5 x 2.5 deg^2) LOFAR High Band Antenna image at a resolution of 1.2'' x 2'' with a median noise of approximately 80 microJy per beam. It was made from an 8-hour International LOFAR Telescope (ILT) observation of the ELAIS-N1 field at frequencies ranging from 120 to 168 MHz with the most up-to-date ILT imaging methods. This intermediate resolution falls between the highest possible resolution (0.3'') achievable by using all International LOFAR Telescope (ILT) baselines and the standard 6-arcsecond resolution in the LoTSS (LOFAR Two-meter Sky Survey) image products utilising the LOFAR Dutch baselines only. This is the first demonstration of the feasibility of imaging using the ILT at a resolution of around 1'', which provides unique information on source morphology at scales that fall below the surface brightness limits at higher resolutions. The total calibration and imaging computational time is approximately 52,000 core hours, nearly 5 times more than required to produce a 6'' resolution image. We also present a radio source catalogue containing 2263 sources detected over the 2.5 x 2.5 deg^2 image of the ELAIS-N1 field, with a peak intensity threshold of 5.5 sigma. The catalogue has been cross-matched with the LoTSS deep ELAIS-N1 field radio catalogue, and its flux density and positional accuracy have been investigated and corrected accordingly. We find that approximately 80% of sources which we expect to be detectable based on their peak brightness in the LoTSS 6'' resolution image are detected in this image, which is approximately a factor of two higher than for 0.3'' resolution imaging in the Lockman Hole, implying there is a wealth of information on these intermediate scales., Comment: 20 pages, 19 figures, revised manuscript submitted to 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

The Low Frequency Array (LOFAR) is the only existing radio interferometer able to observe at ultra-low frequencies (<100 MHz) with high resolution (<15") and high sensitivity (<1 mJy/beam). To exploit these capabilities, the LOFAR Surveys Key Science Project is using the LOFAR Low Band Antenna (LBA) to carry out a sensitive wide-area survey at 41-66 MHz named the LOFAR LBA Sky Survey (LoLSS). LoLSS is covering the whole northern sky above declination 24 deg with a resolution of 15" and a sensitivity of 1-2 mJy/beam (1 sigma) depending on declination, field properties, and observing conditions. Here we present the first data release. An automated pipeline was used to reduce the 95 fields included in this data release. The data reduction procedures developed for this project have general application and are currently being used to process LOFAR LBA interferometric observations. Compared to the preliminary release, direction-dependent errors have been corrected for during the calibration process. This results in a typical sensitivity of 1.55 mJy/beam at the target resolution of 15". The first data release of the LOFAR LBA Sky Survey covers 650 sqdeg in the HETDEX spring field. The resultant data products released to the community include mosaic images (I and V Stokes) of the region, and a catalogue of 42463 detected sources and related Gaussian components used to describe sources' morphologies. Separate catalogues for 6 in-band frequencies are also released. The first data release of LoLSS shows that, despite the influences of the ionosphere, LOFAR can conduct large-scale surveys in the frequency window 42-66 MHz with unprecedentedly high sensitivity and resolution. The data can be used to derive unique information on the low-frequency spectral properties of many thousands of sources with a wide range of applications in extragalactic and galactic astronomy., Comment: 20 pages, 22 figures, images and catalogues available at https://www.lofar-surveys.org/lolss.html

Extended radio sources in the sky require a dense sampling of short baselines to be properly imaged by interferometers. This problem arises in many areas of radio astronomy, such as in the study of galaxy clusters, which may host Mpc-scale diffuse synchrotron sources in the form of radio halos. In clusters where no radio halos are detected, owing to intrinsic absence of emission or extrinsic (instrumental and/or observational) effects, it is possible to determine upper limits. We consider a sample of Planck galaxy clusters from the Second Data Release of the LOFAR Two Meter Sky Survey (LoTSS-DR2) where no radio halos are detected. We use this sample to test the capabilities of LOFAR to recover diffuse extended emission and derive upper limits. Through the injection technique, we simulate radio halos with various surface brightness profiles. We then predict the corresponding visibilities and image them along with the real visibilities. This method allows us to test the fraction of flux density losses owing to inadequate uv-coverage and obtain thresholds at which the mock emission becomes undetectable by visual inspection. The dense uv-coverage of LOFAR at short spacings allows to recover $\gtrsim90\%$ of the flux density of targets with sizes up to $\sim 15'$. We find a relation that provides upper limits based on the image noise and extent (in terms of number of beams) of the mock halo. This relation can be safely adopted to obtain upper limits without injecting when artifacts introduced by the subtraction of the discrete sources are negligible in the central region of the cluster. Otherwise, the injection process and visual inspection of the images are necessary to determine more reliable limits. Through these methods, we obtain upper limits for 75 clusters to be exploited in ongoing statistical studies., Comment: 20 pages (including Appendices). 21 Figures. 4 Tables. Accepted for publication by A&A. For the project Webpage see https://lofar-surveys.org/planck_dr2.html

Context. It is well-established that shock waves in the intracluster medium launched by galaxy cluster mergers can produce synchrotron emission, which is visible to us at radio frequencies as radio relics. However, the particle acceleration mechanism producing these relics is still not fully understood. It is also unclear how relics relate to radio halos, which trace merger-induced turbulence in the intracluster medium. Aims. We aim to perform the first statistical analysis of radio relics in a mass-selected sample of galaxy clusters, using homogeneous observations. Methods. We analysed all relics observed by the Low Frequency Array Two Metre Sky Survey Data Release 2 (LoTSS DR2) at 144 MHz, hosted by galaxy clusters in the second Planck catalogue of SZ sources (PSZ2). We measured and compared the relic properties in a uniform, unbiased way. In particular, we developed a method to describe the characteristic downstream width in a statistical manner. Additionally, we searched for differences between radio relic-hosting clusters with and without radio halos. Results. We find that, in our sample, $\sim$ 10% of galaxy clusters host at least one radio relic. We confirm previous findings, at higher frequencies, of a correlation between the relic-cluster centre distance and the longest linear size, as well as the radio relic power and cluster mass. However, our findings suggest that we are still missing a population of low-power relics. We also find that relics are wider than theoretically expected, even with optimistic downstream conditions. Finally, we do not find evidence of a single property that separates relic-hosting clusters with and without radio halos., Comment: 28 pages, 17 figures. Accepted for publication in A&A

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

We examine the particle acceleration mechanism in the Mpc-scale bridge between Abell 399 and Abell 401 and assess in particular if the synchrotron emission originates from first-order or second-order Fermi re-acceleration. We use deep (~40 hours) LOw-Frequency ARray (LOFAR) observations from Abell 399 and Abell 401 and apply improved direction-dependent calibration to produce deep radio images at three different resolutions at 144 MHz. With a point-to-point analysis we find in the bridge trends between the radio emission from our new maps and X-ray emission from an XMM Newton observation. By analyzing our observations and results, we argue that second-order Fermi re-acceleration is currently the most favoured process to explain the emission from the radio bridge, where past AGN activity may be responsible for the supply of fossil plasma needed for in-situ re-acceleration. The radio halos from Abell 401 and Abell 399 are also consistent with a second-order Fermi re-acceleration model.

The central regions of galaxy clusters are permeated by magnetic fields and filled with relativistic electrons. When clusters merge, the magnetic fields are amplified and relativistic electrons are re-accelerated by turbulence in the intra cluster medium. These electrons reach energies of 1 -- 10 GeV and, in the presence of magnetic fields, produce diffuse radio halos that typically cover an area of ~1 square Mpc. Here we report observations of four clusters whose radio halos are embedded in much more extended, diffuse radio emission, filling a volume 30 times larger than that of radio halos. The emissivity in these larger features is about 20 times lower than the emissivity in radio halos. We conclude that relativistic electrons and magnetic fields extend far beyond radio halos, and that the physical conditions in the outer regions of the clusters are quite different from those in the radio halos., Comment: 21 pages, 8 figures, Nature in press

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

We present LOw Frequency ARray observations of the Coma cluster field at 144\,MHz. The cluster hosts one of the most famous radio halos, a relic, and a low surface-brightness bridge. We detect new features that allow us to make a step forward in the understanding of particle acceleration in clusters. The radio halo extends for more than 2 Mpc, which is the largest extent ever reported. To the North-East of the cluster, beyond the Coma virial radius, we discover an arc-like radio source that could trace particles accelerated by an accretion shock. To the West of the halo, coincident with a shock detected in the X-rays, we confirm the presence of a radio front, with different spectral properties with respect to the rest of the halo. We detect a radial steepening of the radio halo spectral index between 144 MHz and 342 MHz, at $\sim 30^{\prime}$ from the cluster centre, that may indicate a non constant re-acceleration time throughout the volume. We also detect a mild steepening of the spectral index towards the cluster centre. For the first time, a radial change in the slope of the radio-X-ray correlation is found, and we show that such a change could indicate an increasing fraction of cosmic ray versus thermal energy density in the cluster outskirts. Finally, we investigate the origin of the emission between the relic and the source NGC 4789, and we argue that NGC4789 could have crossed the shock originating the radio emission visible between its tail and the relic., Comment: submitted to ApJ

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 empirical relationship between the non-thermal 5GHz radio luminosity and the soft X-ray luminosity of active stellar coronae, canonically called the G\"udel-Benz relationship (G\"udel & Benz 1993), has been a cornerstone of stellar radio astronomy as it explicitly ties the radio emission to the coronal heating mechanisms. The relationship extends from microflares on the Sun to the coronae of the most active stars suggesting that active coronae are heated by a flare-like process (Benz & G\"udel 1994). The relationship is thought to originate from a consistent partition of the available flare energy into relativistic charges, that emit in the radio-band via the incoherent gyrosynchrotron mechanism, and heating of the bulk coronal plasma, that emits in the X-ray band via the Bremsstrahlung mechanism. Consequently, coherent emission from stellar and sub-stellar objects is not expected to adhere to this empirical relationship, as is observed in ultracool dwarf stars and brown dwarfs. Here we report a population of radio-detected chromospherically active stars that surprisingly follows the G\"udel-Benz relationship despite their radio emission being classified as coherent emission by virtue of its high circularly polarised fraction and high brightness temperature. Our results prompt a re-examination of the physics behind the G\"udel-Benz relationship, its implication for the mechanism of coronal heating and particle acceleration in active stars and the phenomenological connection between solar and stellar flares., Comment: Under review, ApJL

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

Coherent low-frequency ($\lesssim 200$ MHz) radio emission from stars encodes the conditions of the outer corona, mass-ejection events, and space weather. Previous low-frequency searches for radio emitting stellar systems have lacked the sensitivity to detect the general population, instead largely focusing on targeted studies of anomalously active stars. Here we present 19 detections of coherent radio emission associated with known M~dwarfs from a blind flux-limited low-frequency survey. Our detections show that coherent radio emission is ubiquitous across the M~dwarf main sequence, and that the radio luminosity is independent of known coronal and chromospheric activity indicators. While plasma emission can generate the low-frequency emission from the most chromospherically active stars of our sample, the origin of the radio emission from the most quiescent sources is yet to be ascertained. Large-scale analogues of the magnetospheric processes seen in gas-giant planets likely drive the radio emission associated with these quiescent stars. The slowest-rotating stars of this sample are candidate systems to search for star-planet interaction signatures., Comment: Accepted to Nature Astronomy; 31 pages, 9 figures, and 2 tables

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 recent years, thanks to new facilities such as LOFAR capable of sensitive observations, much work has been done on the detection of stellar radio emission at low frequencies. Such emission has commonly been shown to be coherent emission, generally attributed to electron-cyclotron maser emission, and has usually been detected from main-sequence M dwarfs. Here we report the first detection of coherent emission at low frequencies from T Tauri stars, which are known to be associated with high levels of stellar activity. Using LOFAR, we have detected several bright radio bursts at 150 MHz from two weak-line T Tauri stars: KPNO-Tau 14 and LkCa 4. All of the bursts have high brightness temperatures ($10^{13} - 10^{14}\ \mathrm{K}$) and high circular polarization fractions (60 - 90 \%), indicating that they must be due to a coherent emission mechanism. This could be either plasma emission or electron-cyclotron maser (ECM) emission. Due to the exceptionally high brightness temperatures seen in at least one of the bursts ($\geq 10^{14}\ \mathrm{K}$), and the high circular polarization levels, it seems unlikely that plasma emission could be the source and so ECM is favoured as the most likely emission mechanism. Assuming this is the case, the required magnetic field in the emission regions would be 40 - 70 G. We determine that the most likely method of generating ECM emission is plasma co-rotation breakdown in the stellar magnetosphere. There remains the possibility, however, it could be due to an interaction with an orbiting exoplanet., Comment: 16 pages, 9 figures, accepted for publication in A&A

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

LOFAR is the only radio telescope that is presently capable of high-sensitivity, high-resolution (<1 mJy/b and <15") observations at ultra-low frequencies (<100 MHz). To utilise these capabilities, the LOFAR Surveys Key Science Project is undertaking a large survey to cover the entire northern sky with Low Band Antenna (LBA) observations. The LOFAR LBA Sky Survey (LoLSS) aims to cover the entire northern sky with 3170 pointings in the frequency range 42-66 MHz, at a resolution of 15 arcsec and at a sensitivity of 1 mJy/b. Here we outline the survey strategy, the observational status, the current calibration techniques, and briefly describe several scientific motivations. We also describe the preliminary public data release. The preliminary images were produced using a fully automated pipeline that aims to correct all direction-independent effects in the data. Whilst the direction-dependent effects, such as those from the ionosphere, are not yet corrected, the images presented in this work are still 10 times more sensitive than previous surveys available at these low frequencies. The preliminary data release covers 740 sqdeg around the HETDEX spring field region at a resolution of 47" with a median noise level of 5 mJy/b. The images and the catalogue with 25,247 sources are publicly released. We demonstrate that the system is capable of reaching an rms noise of 1 mJy/b and the resolution of 15" once direction-dependent effects are corrected for. LoLSS will provide the ultra-low-frequency information for hundreds of thousands of radio sources, providing critical spectral information and producing a unique dataset that can be used for a wide range of science topics such as: the search for high redshift galaxies and quasars, the study of the magnetosphere of exoplanets, and the detection of the oldest populations of cosmic-rays in galaxies, clusters of galaxies, and from AGN activity., Comment: 19 pages, 14 figures, Accepted A&A, catalogue and images on www.lofar-surveys.org

Recently detected coherent low-frequency radio emission from M dwarf systems shares phenomenological similarities with emission produced by magnetospheric processes from the gas giant planets of our Solar System. Such beamed electron-cyclotron maser emission can be driven by a star-planet interaction or a breakdown in co-rotation between a rotating plasma disk and a stellar magnetosphere. Both models suggest that the radio emission could be periodic. Here we present the longest low-frequency interferometric monitoring campaign of an M dwarf system, composed of twenty-one $\approx$8 hour epochs taken in two series of observing blocks separated by a year. We achieved a total on-source time of 6.5 days. We show that the M dwarf binary CR Draconis has a low-frequency 3$\sigma$ detection rate of 90$^{+5}_{-8}$% when a noise floor of $\approx$0.1 mJy is reached, with a median flux density of 0.92 mJy, consistent circularly polarised handedness, and a median circularly polarised fraction of 66%. We resolve three bright radio bursts in dynamic spectra, revealing the brightest is elliptically polarised, confined to 4 MHz of bandwidth centred on 170 MHz, and reaches a flux density of 205 mJy. The burst structure is mottled, indicating it consists of unresolved sub-bursts. Such a structure shares a striking resemblance with the low-frequency emission from Jupiter. We suggest the near-constant detection of high brightness temperature, highly-circularly-polarised radiation that has a consistent circular polarisation handedness implies the emission is produced via the electron-cyclotron maser instability. Optical photometric data reveal the system has a rotation period of 1.984$\pm$0.003 days. We observe no periodicity in the radio data, but the sampling of our radio observations produces a window function that would hide the near two-day signal., Comment: Accepted for publication in A&A, 16 pages, 7 figures, 2 tables

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

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