Showing 1–50 of 352 results for author: Hardcastle, M J

J0011+3217: A peculiar radio galaxy with a one-sided secondary lobe and misaligned giant primary lobes

Authors: Shobha Kumari, Sabyasachi Pal, Martin J. Hardcastle, Maya A. Horton

Abstract: From the LOFAR Two-metre Sky Survey second data release (LoTSS DR2) at 144 MHz, we identified a peculiar radio galaxy, J0011+3217. It has a large, one-sided diffuse secondary wing that stretches up to 0.85 Mpc (roughly 85\% of the size of the primary lobe). The linear size of the primary lobe of the galaxy is 0.99 Mpc. This peculiar source is a giant radio galaxy with a misaligned primary lobe. Th… ▽ More From the LOFAR Two-metre Sky Survey second data release (LoTSS DR2) at 144 MHz, we identified a peculiar radio galaxy, J0011+3217. It has a large, one-sided diffuse secondary wing that stretches up to 0.85 Mpc (roughly 85\% of the size of the primary lobe). The linear size of the primary lobe of the galaxy is 0.99 Mpc. This peculiar source is a giant radio galaxy with a misaligned primary lobe. There is an optical galaxy 16 kpc (7 arcsec) from the host active galactic nucleus of J0011+3217. J0011+3217 has a radio luminosity of $1.65\times 10^{26}$ W Hz$^{-1}$ at 144 MHz with a spectral index of $-0.80$ between 144 and 607 MHz. J0011+3217 is located 1.2 Mpc from the centre of the Abell 7 cluster. The Abell 7 cluster has a redshift of 0.104 and a mass ($M_{500}$) of 3.71 $\times 10^{14}$ M$_\odot$. The cluster is associated with strong X-ray emission. We studied the X-ray emission around the cluster and from the region surrounding J0011+3217 using an XMM-Newton image of J0011+3217, and we analysed the velocity structure and spatial distribution of galaxies in the cluster, showing that J0011+3217 inhabits an offset group of galaxies that are moving with respect to Abell 7. The off-axis distortion, or bending, of the primary lobe of J0011+3217 in the outer edges has a strong effect on the relative motion of the surrounding medium; this in turn causes the bending of the jets in the opposite direction of the cluster (like wide-angle tailed sources). We suggest that the morphology of J0011+3217 is influenced by ram pressure created by the Abell 7 cluster, highlighting the complex interactions between the source and the surrounding cluster environment. △ Less

Submitted 9 July, 2024; v1 submitted 21 June, 2024; originally announced June 2024.

Comments: Accepted for publication in Astronomy and Astrophysics

The LOFAR Two-metre Sky Survey: The nature of the faint source population and SFR-radio luminosity relation using Prospector

Authors: Soumyadeep Das, Daniel J. B. Smith, Paul Haskell, Martin J. Hardcastle, Philip N. Best, Kenneth J. Duncan, Marina I. Arnaudova, Shravya Shenoy, Rohit Kondapally, Rachel K. Cochrane, Alyssa B. Drake, Gülay Gürkan, Katarzyna Małek, Leah K. Morabito, Isabella Prandoni

Abstract: Spectral energy distribution (SED) fitting has been extensively used to determine the nature of the faint radio source population. Recent efforts have combined fits from multiple SED-fitting codes to account for the host galaxy and any active nucleus that may be present. We show that it is possible to produce similar-quality classifications using a single energy-balance SED fitting code, Prospecto… ▽ More Spectral energy distribution (SED) fitting has been extensively used to determine the nature of the faint radio source population. Recent efforts have combined fits from multiple SED-fitting codes to account for the host galaxy and any active nucleus that may be present. We show that it is possible to produce similar-quality classifications using a single energy-balance SED fitting code, Prospector, to model up to 26 bands of UV$-$far-infrared aperture-matched photometry for $\sim$31,000 sources in the ELAIS-N1 field from the LOFAR Two-Metre Sky Survey (LoTSS) Deep fields first data release. One of a new generation of SED-fitting codes, Prospector accounts for potential contributions from radiative active galactic nuclei (AGN) when estimating galaxy properties, including star formation rates (SFRs) derived using non-parametric star formation histories. Combining this information with radio luminosities, we classify 92 per cent of the radio sources as a star-forming galaxy, high/low-excitation radio galaxy, or radio-quiet AGN and study the population demographics as a function of 150 MHz flux density, luminosity, SFR, stellar mass, redshift and apparent $r$-band magnitude. Finally, we use Prospector SED fits to investigate the SFR$-$150 MHz luminosity relation for a sample of $\sim$$133,000~3.6~μ$m-selected $z<1$ sources, finding that the stellar mass dependence is significantly weaker than previously reported, and may disappear altogether at $\log_{10} (\mathrm{SFR}/M_\odot~\mathrm{yr}^{-1}) > 0.5$. This approach makes it significantly easier to classify radio sources from LoTSS and elsewhere, and may have important implications for future studies of star-forming galaxies at radio wavelengths. △ Less

Submitted 2 May, 2024; originally announced May 2024.

Comments: 20 pages, 13 figures, accepted for publication in MNRAS

Constraining the giant radio galaxy population with machine learning and Bayesian inference

Authors: Rafaël I. J. Mostert, Martijn S. S. L. Oei, B. Barkus, Lara Alegre, Martin J. Hardcastle, Kenneth J. Duncan, Huub J. A. Röttgering, Reinout J. van Weeren, Maya Horton

Abstract: Large-scale sky surveys at low frequencies, like the LOFAR Two-metre Sky Survey (LoTSS), allow for the detection and characterisation of unprecedented numbers of giant radio galaxies (GRGs, or 'giants'). In this work, by automating the creation of radio--optical catalogues, we aim to significantly expand the census of known giants. We then combine this sample with a forward model to constrain GRG… ▽ More Large-scale sky surveys at low frequencies, like the LOFAR Two-metre Sky Survey (LoTSS), allow for the detection and characterisation of unprecedented numbers of giant radio galaxies (GRGs, or 'giants'). In this work, by automating the creation of radio--optical catalogues, we aim to significantly expand the census of known giants. We then combine this sample with a forward model to constrain GRG properties of cosmological interest. In particular, we automate radio source component association through machine learning and optical host identification for resolved radio sources. We create a radio--optical catalogue for the full LoTSS Data Release 2 (DR2) and select all possible giants. We combine our candidates with an existing catalogue of LoTSS DR2 crowd-sourced GRG candidates and visually confirm or reject them. To infer intrinsic GRG properties from GRG observations, we develop further a population-based forward model that takes into account selection effects and constrain its parameters using Bayesian inference. We confirm 5,647 previously unknown giants from the crowd-sourced catalogue and 2,597 previously unknown giants from the ML-driven catalogue. Our confirmations and discoveries bring the total number of known giants to at least 11,585. We predict a comoving GRG number density $n_\mathrm{GRG} = 13 \pm 10\ (100\ \mathrm{Mpc})^{-3}$, close to a recent estimate of the number density of luminous non-giant radio galaxies. We derive a current-day GRG lobe volume-filling fraction $V_\mathrm{GRG-CW}(z = 0) = 1.4 \pm 1.1 \cdot 10^{-5}$ in clusters and filaments of the Cosmic Web. Our analysis suggests that giants are more common than previously thought. Moreover, tentative results imply that it is possible that magnetic fields once contained in giants pervade a significant ($\gtrsim 10\%$) fraction of today's Cosmic Web. △ Less

Submitted 30 April, 2024; originally announced May 2024.

Comments: 25 pages, 15 figures, 3 tables, accepted for publication in Astronomy & Astrophysics

Luminous giants populate the dense Cosmic Web: The radio luminosity-environmental density relation for radio galaxies in action

Authors: Martijn S. S. L. Oei, Reinout J. van Weeren, Martin J. Hardcastle, Aivin R. D. J. G. I. B. Gast, Florent Leclercq, Huub J. A. Röttgering, Pratik Dabhade, Tim W. Shimwell, Andrea Botteon

Abstract: Giant radio galaxies (GRGs, giant RGs, or giants) are megaparsec-scale, jet-driven outflows from accretion disks of supermassive black holes, and represent the most extreme pathway by which galaxies can impact the Cosmic Web around them. A long-standing but unresolved question is why giants are so much larger than other radio galaxies. It has been proposed that, in addition to having higher jet po… ▽ More Giant radio galaxies (GRGs, giant RGs, or giants) are megaparsec-scale, jet-driven outflows from accretion disks of supermassive black holes, and represent the most extreme pathway by which galaxies can impact the Cosmic Web around them. A long-standing but unresolved question is why giants are so much larger than other radio galaxies. It has been proposed that, in addition to having higher jet powers than most RGs, giants might live in especially low-density Cosmic Web environments. In this work, we aim to test this hypothesis by pinpointing Local Universe giants and other RGs in physically principled, Bayesian large-scale structure reconstructions. More specifically, we localised a LOFAR Two-metre Sky Survey (LoTSS) DR2-dominated sample of luminous ($l_ν(ν= 150\ \mathrm{MHz}) \geq 10^{24}\ \mathrm{W\ Hz^{-1}}$) giants and a control sample of LoTSS DR1 RGs, both with spectroscopic redshifts up to $z_\mathrm{max} = 0.16$, in the BORG SDSS Cosmic Web reconstructions. We measured the Cosmic Web density for each RG; for the control sample, we then quantified the relation between RG radio luminosity and Cosmic Web density. With the BORG SDSS tidal tensor, we also measured for each RG whether the gravitational dynamics of its Cosmic Web environment resemble those of clusters, filaments, sheets, or voids. Luminous giants populate large-scale environments that tend to be denser than those of general RGs. This shows that -- at least at high jet powers -- low-density environments are no prerequisite for giant growth. This result is corroborated by gravitational dynamics classification and a cluster catalogue crossmatching analysis. This work presents more than a thousand inferred megaparsec-scale densities around radio galaxies. Our findings are consistent with the view that giants are regular, rather than mechanistically special, members of the radio galaxy population. △ Less

Submitted 27 April, 2024; originally announced April 2024.

Comments: 28 pages, 21 figures, 4 tables, accepted for publication in Astronomy & Astrophysics

The nature of compact radio-loud AGN: a systematic look at the LOFAR AGN population

Authors: J. Chilufya, M. J. Hardcastle, J. C. S. Pierce, J. H. Croston, B. Mingo, X. Zheng, R. D. Baldi, H. J. A. Röttgering

Abstract: We investigate the nature of low-luminosity radio-loud active galactic nuclei (RLAGN) selected from the LOFAR Two-metre Sky Survey (LoTSS) first data release (DR1). Using optical, mid-infrared, and radio data, we have conservatively selected 55 radiative AGN candidates from DR1 within the redshift range $0.03<z<0.1$. We show using high-frequency {\it Karl G. Jansky} Very Large Array (VLA) observat… ▽ More We investigate the nature of low-luminosity radio-loud active galactic nuclei (RLAGN) selected from the LOFAR Two-metre Sky Survey (LoTSS) first data release (DR1). Using optical, mid-infrared, and radio data, we have conservatively selected 55 radiative AGN candidates from DR1 within the redshift range $0.03<z<0.1$. We show using high-frequency {\it Karl G. Jansky} Very Large Array (VLA) observations that 10 out of 55 objects show radio emission on scales $>$$1-3$ kpc, 42 are compact at the limiting resolution of 0.35 arcsec (taking an upper limit on the projected physical size, this corresponds to less than 1 kpc), and three are undetected. The extended objects display a wide range of radio morphologies: two-jet (5), one-jet (4), and double-lobed (1). We present the radio spectra of all detected radio sources which range from steep to flat/inverted and span the range seen for other compact radio sources such as compact symmetric objects (CSOs), compact steep spectrum (CSS) sources, and gigahertz peaked-spectrum (GPS) sources. Assuming synchrotron self-absorption (SSA) for flat/inverted radio spectrum sources, we predict small physical sizes for compact objects to range between $2-53$ pc. Alternatively, using free-free absorption (FFA) models, we have estimated the free electron column depth for all compact objects, assuming a homogeneous absorber. We find that these objects do not occupy a special position on the power/linear size ($P-D$) diagram but some share a region with radio-quiet quasars (RQQs) and so-called `FR0' sources in terms of radio luminosity and linear size. △ Less

Submitted 29 February, 2024; originally announced February 2024.

Comments: Accepted for publication in the Monthly Notices of the Royal Astronomy Society (MNRAS) journal

The LOFAR-eFEDS survey: The incidence of radio and X-ray AGN and the disk-jet connection

Authors: Z. Igo, A. Merloni, D. Hoang, J. Buchner, T. Liu, M. Salvato, R. Arcodia, S. Bellstedt, M. Brüggen, J. H. Croston, F. de Gasperin, A. Georgakakis, M. J. Hardcastle, K. Nandra, Q. Ni, T. Pasini, T. Shimwell, J. Wolf

Abstract: Radio jets are present in a diverse sample of AGN. However, the mechanisms of jet powering are not fully understood, and it is yet unclear to what extent they obey mass-invariant scaling relations, similar to those found for the triggering and fuelling of X-ray selected AGN. We study the incidence of eROSITA/eFEDS X-ray and LOFAR radio AGN as a function of several stellar mass normalised AGN power… ▽ More Radio jets are present in a diverse sample of AGN. However, the mechanisms of jet powering are not fully understood, and it is yet unclear to what extent they obey mass-invariant scaling relations, similar to those found for the triggering and fuelling of X-ray selected AGN. We study the incidence of eROSITA/eFEDS X-ray and LOFAR radio AGN as a function of several stellar mass normalised AGN power indicators. A new sample of radio AGN from the LOFAR-eFEDS survey is defined and we publicly release this catalogue, including host galaxy counterparts from the Legacy Survey DR9, LOFAR radio morphologies and host galaxy properties from the complete, spectroscopic (z<0.4) GAMA09 survey. The fraction of GAMA09 galaxies hosting radio, X-ray and both radio and X-ray AGN are calculated as a function of the specific black hole kinetic ($λ_{\rm Jet}$) and radiative ($λ_{\rm Edd}$) power. The incidence of eFEDS X-ray AGN as a function of $λ_{\rm Edd}$ shows the same mass-invariance as found in past studies. Meanwhile, radio AGN, regardless of their morphology, are more likely to be hosted in more massive galaxies, at all $λ_{\rm Jet}$. Across the stellar mass range, the compact radio AGN incidence follows the same power-law distribution, showing that it is not only high mass galaxies that host high power radio AGN and vice versa. On the other hand, the incidence of compact and complex radio AGN is boosted at the highest jet powers, diverging from a simple power-law. Interestingly, this increased incidence cannot be explained by more powerful radio AGN lying in more dense environments which could naturally boost their radio luminosity. Overall, we show that statistical incidence studies are a powerful method to probe disk-jet coupling for different AGN accretion modes, although future work on a more reliable determination of jet power for diverse samples of radio AGN is needed. △ Less

Submitted 26 February, 2024; originally announced February 2024.

Comments: 30 pages, 32 figures, accepted for publication in A&A

Journal ref: A&A 686, A43 (2024)

Exploring the radio-loudness of SDSS quasars with spectral stacking

Authors: M. I. Arnaudova, D. J. B. Smith, M. J. Hardcastle, S. Das, A. Drake, K. Duncan, G. Gürkan, M. Magliocchetti, L. K. Morabito, J. W. Petley, S. Shenoy, C. Tasse

Abstract: We use new 144 MHz observations over 5634 deg$^2$ from the LOFAR Two-metre Sky Survey (LoTSS) to compile the largest sample of uniformly-selected, spectroscopically-confirmed quasars from the 14th data release of the Sloan Digital Sky Survey (SDSS-DR14). Using the classical definition of radio-loudness, $R=\log(L_{\rm{1.4GHz}}/L_{i})$, we identify 3,697 radio-loud (RL) and 111,132 radio-quiet (RQ)… ▽ More We use new 144 MHz observations over 5634 deg$^2$ from the LOFAR Two-metre Sky Survey (LoTSS) to compile the largest sample of uniformly-selected, spectroscopically-confirmed quasars from the 14th data release of the Sloan Digital Sky Survey (SDSS-DR14). Using the classical definition of radio-loudness, $R=\log(L_{\rm{1.4GHz}}/L_{i})$, we identify 3,697 radio-loud (RL) and 111,132 radio-quiet (RQ) sources at $0.6<z<3.4$. To study their properties, we develop a new rest-frame spectral stacking algorithm, designed with forthcoming massively-multiplexed spectroscopic surveys in mind, and use it to create high signal-to-noise composite spectra of each class, matched in redshift and absolute $i$-band magnitude. We show that RL quasars have redder continuum and enhanced [OII] emission than their RQ counterparts. These results persist when additionally matching in black hole mass, suggesting that this parameter is not the defining factor in making a QSO radio-loud. We find that these features are not gradually varying as a function of radio-loudness but are maintained even when probing deeper into the RQ population, indicating that a clear-cut division in radio-loudness is not apparent. Upon examining the star formation rates (SFRs) inferred from the [OII] emission line, with the contribution from AGN removed using the [NeV] line, we find that RL quasars have a significant excess of star-formation relative to RQ quasars out to $z=1.9$ at least. Given our findings, we suggest that radio-loud sources either preferably reside in gas-rich systems with rapidly-spinning black holes, or represent an earlier obscured phase of QSO evolution. △ Less

Submitted 16 January, 2024; originally announced January 2024.

Comments: 22 pages, 17 figures, accepted for publication in MNRAS

Ubiquitous radio emission in quasars: predominant AGN origin and a connection to jets, dust and winds

Authors: G. Calistro Rivera, D. M. Alexander, C. M. Harrison, V. A. Fawcett, P. N. Best, W. L. Williams, M. J. Hardcastle, D. J. Rosario, D. J. B. Smith, M. I. Arnaudova, E. Escott, G. Gürkan, R. Kondapally, G. Miley, L. K. Morabito, J. Petley, I. Prandoni, H. J. A. Röttgering, B. -H. Yue

Abstract: We present a comprehensive study of the physical origin of radio emission in optical quasars at redshifts z < 2.5. We focus particularly on the associations between compact radio emission, dust reddening, and outflows identified in our earlier work. Leveraging the deepest low-frequency radio data available to date (LoTSS Deep DR1), we achieve radio detection fractions of up to 94%, demonstrating t… ▽ More We present a comprehensive study of the physical origin of radio emission in optical quasars at redshifts z < 2.5. We focus particularly on the associations between compact radio emission, dust reddening, and outflows identified in our earlier work. Leveraging the deepest low-frequency radio data available to date (LoTSS Deep DR1), we achieve radio detection fractions of up to 94%, demonstrating the virtual ubiquity of radio emission in quasars, and a continuous distribution in radio loudness. Through our analysis of radio properties, combined with spectral energy distribution modeling of multiwavelength photometry, we establish that the primary source of radio emission in quasars is the AGN, rather than star formation. Modeling the dust reddening of the accretion disk emission shows a continuous increase in radio detection in quasars as a function of the reddening parameter E(B-V), suggesting a causal link between radio emission and dust reddening. Confirming previous findings, we observe that the radio excess in red quasars is most pronounced for sources with compact radio morphologies and intermediate radio loudness. We find a significant increase in [Oiii] and Civ outflow velocities for red quasars not seen in our control sample, with particularly powerful [Oiii] winds in those around the radio-quiet/radio-loud threshold. Based on the combined characterisation of radio, reddening, and wind properties in our sample, we favor a model in which the compact radio emission observed in quasars originates in compact radio jets and their interaction with a dusty, circumnuclear environment. Our results align with the theory that jet-induced winds and shocks resulting from this interaction are the origin of the enhanced radio emission in red quasars. Further investigation of this model is crucial for advancing our understanding of quasar feedback mechanisms and their role in galaxy evolution. △ Less

Submitted 15 December, 2023; originally announced December 2023.

Comments: 11 pages, 8 Figures, submitted to A&A

LOFAR HBA Observations of the Euclid Deep Field North (EDFN)

Authors: M. Bondi, R. Scaramella, G. Zamorani, P. Ciliegi, F. Vitello, M. Arias, P. N. Best, M. Bonato, A. Botteon, M. Brienza, G. Brunetti, M. J. Hardcastle, M. Magliocchetti, F. Massaro, L. K. Morabito, L. Pentericci, I. Prandoni, H. J. A. Röttgering, T. W. Shimwell, C. Tasse, R. J. van Weeren, G. J. White

Abstract: 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 ce… ▽ More 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\,μ$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. △ Less

Submitted 11 December, 2023; originally announced December 2023.

Comments: accepted for publication in A&A

Numerical modelling of the lobes of radio galaxies -- Paper V: Universal Pressure Profile cluster atmospheres

Authors: Michael Stimpson, Martin J. Hardcastle, Martin G. H. Krause

Abstract: We present relativistic magnetohydrodynamic modelling of jets running into hydrostatic, spherically symmetric cluster atmospheres. For the first time in a numerical simulation, we present model cluster atmospheres based upon the Universal Pressure Profile (UPP), incorporating a temperature profile for a typical self-similar atmosphere described by only one parameter - $M_{500}$. We explore a compr… ▽ More We present relativistic magnetohydrodynamic modelling of jets running into hydrostatic, spherically symmetric cluster atmospheres. For the first time in a numerical simulation, we present model cluster atmospheres based upon the Universal Pressure Profile (UPP), incorporating a temperature profile for a typical self-similar atmosphere described by only one parameter - $M_{500}$. We explore a comprehensive range of realistic atmospheres and jet powers and derive dynamic, energetic and polarimetric data which provide insight into what we should expect of future high-resolution studies of AGN outflows. From the simulated synchrotron emission maps which include Doppler beaming we find sidedness distributions that agree well with observations. We replicated a number of findings from our previous work, such as higher power jets inflating larger aspect-ratio lobes and the cluster environment impacting the distribution of energy between the lobe and shocked regions. Comparing UPP and $β$-profiles we find that the cluster model chosen results in a different morphology for the resultant lobes with the UPP more able to clear lobe material from the core; and that these different atmospheres influence the ratio between the various forms of energy in the fully developed lobes. This work also highlights the key role played by Kelvin-Helmholtz (KH) instabilities in the formation of realistic lobe aspect-ratios. Our simulations point to the need for additional lobe-widening mechanisms at high jet powers, for example jet precession. Given that the UPP is our most representative general cluster atmosphere, these numerical simulations represent the most realistic models yet for spherically symmetric atmospheres. △ Less

Submitted 4 December, 2023; v1 submitted 1 October, 2023; originally announced October 2023.

Comments: 20 pages, 18 figures, 2 tables, accepted for publication in MNRAS, minor correction to references in section 4.2

No strong radio absorption detected in the low-frequency spectra of radio-loud quasars at z > 5.6

Authors: A. J. Gloudemans, A. Saxena, H. Intema, J. R. Callingham, K. J. Duncan, H. J. A. Rottgering, S. Belladitta, M. J. Hardcastle, Y. Harikane, C. Spingola

Abstract: We present the low-frequency radio spectra of 9 high-redshift quasars at $5.6 \leq z \leq 6.6$ using the Giant Metre Radio Telescope band-3, -4, and -5 observations ($\sim$300-1200 MHz), archival Low Frequency Array (LOFAR; 144 MHz), and Very Large Array (VLA; 1.4 and 3 GHz) data. Five of the quasars in our sample have been discovered recently, representing some of the highest redshift radio brigh… ▽ More We present the low-frequency radio spectra of 9 high-redshift quasars at $5.6 \leq z \leq 6.6$ using the Giant Metre Radio Telescope band-3, -4, and -5 observations ($\sim$300-1200 MHz), archival Low Frequency Array (LOFAR; 144 MHz), and Very Large Array (VLA; 1.4 and 3 GHz) data. Five of the quasars in our sample have been discovered recently, representing some of the highest redshift radio bright quasars known at low-frequencies. We model their radio spectra to study their radio emission mechanism and age of the radio jets by constraining the spectral turnover caused by synchrotron self-absorption (SSA) or free-free absorption (FFA). Besides J0309+2717, a blazar at $z=6.1$, our quasars show no sign of a spectral flattening between 144 MHz and a few GHz, indicating there is no strong SSA or FFA absorption in the observed frequency range. However, we find a wide range of spectral indices between $-1.6$ and $0.05$, including the discovery of 3 potential ultra-steep spectrum quasars. Using further archival VLBA data, we confirm that the radio SED of the blazar J0309+2717 likely turns over at a rest-frame frequency of 0.6-2.3 GHz (90-330 MHz observed frame), with a high-frequency break indicative of radiative ageing of the electron population in the radio lobes. Ultra-low frequency data below 50 MHz are necessary to constrain the absorption mechanism for J0309+2717 and the turnover frequencies for the other high-$z$ quasars in our sample. A relation between linear radio jet size and turnover frequency has been established at low redshifts. If this relation were to hold at high redshifts, the limits on the turnover frequency of our sample suggest the radio jet sizes must be more extended than the typical sizes observed in other radio-bright quasars at similar redshift. To confirm this deep radio follow-up observations with high spatial resolution are required. △ Less

Submitted 19 September, 2023; v1 submitted 7 September, 2023; originally announced September 2023.

Comments: 15 pages, 7 figures, accepted for publication in A&A, minor correction Tab. 1

Journal ref: A&A 678, A128 (2023)

Plausible association of distant late M dwarfs with low-frequency radio emission

Authors: A. J. Gloudemans, J. R. Callingham, K. J. Duncan, A. Saxena, Y. Harikane, G. J. Hill, G. R. Zeimann, H. J. A. Rottgering, M. J. Hardcastle, J. S. Pineda, T. W. Shimwell, D. J. B. Smith, J. D. Wagenveld

Abstract: We present the serendipitous discovery of 8 distant ($>$ 50 pc) late M dwarfs with plausible associated radio emission at 144 MHz. The M dwarf nature of our sources has been confirmed with optical spectroscopy performed using HET/LRS2 and Subaru/FOCAS, and their radio flux densities are within the range of 0.5-1.0 mJy at 144 MHz. Considering the radio-optical source separation and source densities… ▽ More We present the serendipitous discovery of 8 distant ($>$ 50 pc) late M dwarfs with plausible associated radio emission at 144 MHz. The M dwarf nature of our sources has been confirmed with optical spectroscopy performed using HET/LRS2 and Subaru/FOCAS, and their radio flux densities are within the range of 0.5-1.0 mJy at 144 MHz. Considering the radio-optical source separation and source densities of the parent catalogues, we suggest that it is statistically probable the M dwarfs are associated with the radio emission. However, it remains plausible that for some of the sources the radio emission originates from an optically faint and red galaxy hiding behind the M dwarf. The isotropic radio luminosities ($\sim10^{17-18}$ erg s$^{-1}$ Hz$^{-1}$) of the M dwarfs suggest that if the association is real, the radio emission is likely driven by a coherent emission process produced via plasma or electron-cyclotron maser instability processes, which is potentially caused by binary interaction. Long term monitoring in the radio and high-resolution radio follow-up observations are necessary to search for any variability and pinpoint the radio emission to determine whether our tentative conclusion that these ultracool dwarfs are radio emitting is correct. If the low-frequency radio emission is conclusively associated with the M dwarfs, this would reveal a new population of optically faint and distant ($>$ 50 pc) radio emitting M dwarfs. △ Less

Submitted 4 September, 2023; originally announced September 2023.

Comments: 10 pages, 5 figures, accepted for publication in A&A

Journal ref: A&A 678, A161 (2023)

The LOFAR Two-Metre Sky Survey (LoTSS): VI. Optical identifications for the second data release

Authors: M. J. Hardcastle, M. A. Horton, W. L. Williams, K. J. Duncan, L. Alegre, B. Barkus, J. H. Croston, H. Dickinson, E. Osinga, H. J. A. Röttgering, J. Sabater, T. W. Shimwell, D. J. B. Smith, P. N. Best, A. Botteon, M. Brüggen, A. Drabent, F. de Gasperin, G. Gürkan, M. Hajduk, C. L. Hale, M. Hoeft, M. Jamrozy, M. Kunert-Bajraszewska, R. Kondapally , et al. (27 additional authors not shown)

Abstract: The second data release of the LOFAR Two-Metre Sky Survey (LoTSS) covers 27% of the northern sky, with a total area of $\sim 5,700$ deg$^2$. The high angular resolution of LOFAR with Dutch baselines (6 arcsec) allows us to carry out optical identifications of a large fraction of the detected radio sources without further radio followup; however, the process is made more challenging by the many ext… ▽ More The second data release of the LOFAR Two-Metre Sky Survey (LoTSS) covers 27% of the northern sky, with a total area of $\sim 5,700$ deg$^2$. The high angular resolution of LOFAR with Dutch baselines (6 arcsec) allows us to carry out optical identifications of a large fraction of the detected radio sources without further radio followup; however, the process is made more challenging by the many extended radio sources found in LOFAR images as a result of its excellent sensitivity to extended structure. In this paper we present source associations and identifications for sources in the second data release based on optical and near-infrared data, using a combination of a likelihood-ratio cross-match method developed for our first data release, our citizen science project Radio Galaxy Zoo: LOFAR, and new approaches to algorithmic optical identification, together with extensive visual inspection by astronomers. We also present spectroscopic or photometric redshifts for a large fraction of the optical identifications. In total 4,116,934 radio sources lie in the area with good optical data, of which 85% have an optical or infrared identification and 58% have a good redshift estimate. We demonstrate the quality of the dataset by comparing it with earlier optically identified radio surveys. This is by far the largest ever optically identified radio catalogue, and will permit robust statistical studies of star-forming and radio-loud active galaxies. △ Less

Submitted 31 August, 2023; originally announced September 2023.

Comments: 29 pages. Accepted by A&A; data products available at https://lofar-surveys.org/dr2_release.html

Journal ref: A&A 678, A151 (2023)

Cosmic evolution of radio-AGN feedback: confronting models with data

Authors: R. Kondapally, P. N. Best, M. Raouf, N. L. Thomas, R. Davé, S. S. Shabala, H. J. A. Röttgering, M. J. Hardcastle, M. Bonato, R. K. Cochrane, K. Małek, L. K. Morabito, I. Prandoni, D. J. B. Smith

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

Submitted 20 June, 2023; originally announced June 2023.

Comments: 15 pages, 6 figures, accepted for publication in MNRAS

ViCTORIA project: The LOFAR HBA Virgo Cluster Survey

Authors: H. W. Edler, F. de Gasperin, T. W. Shimwell, M. J. Hardcastle, A. Boselli, V. Heesen, H. McCall, D. J. Bomans, M. Brüggen, E. Bulbul, K. T. Chŷzy, A. Ignesti, A. Merloni, F. Pacaud, T. H. Reiprich, I. D. Roberts, H. J. A. Rottgering, R. J. van Weeren

Abstract: The Virgo cluster is the nearest massive galaxy cluster and thus a prime target to study astrophysical processes in dense large-scale environments. In the radio band, we can probe the non-thermal components of the inter-stellar medium (ISM), intracluster medium (ICM) and of active galactic nuclei (AGN). With the ViCTORIA (Virgo Cluster multi-Telescope Observations in Radio of Interacting galaxies… ▽ More The Virgo cluster is the nearest massive galaxy cluster and thus a prime target to study astrophysical processes in dense large-scale environments. In the radio band, we can probe the non-thermal components of the inter-stellar medium (ISM), intracluster medium (ICM) and of active galactic nuclei (AGN). With the ViCTORIA (Virgo Cluster multi-Telescope Observations in Radio of Interacting galaxies and AGN) project, we are carrying out multiple wide-field surveys of the Virgo cluster at different frequencies. We aim to investigate the impact of the environment on the evolution of galaxies and the contribution of AGN to the ICM-heating, from the inner cluster regions out to beyond the virial radius. We present a survey of the cluster at 120-168 MHz using LOFAR. We image a 132 deg$^2$ region of the cluster, reaching an order of magnitude greater sensitivity than existing wide-field radio surveys of this field at three times higher spatial resolution compared to other low-frequency observations. We developed a tailored data processing strategy to subtract the bright central radio galaxy M87 from the data. This allowed us to correct for the systematic effects due to ionospheric variation as a function of time and direction. In the final mosaic with a resolution of 9"x5", we reach a median noise level of 140 $μ$Jy/beam inside the virial radius and 280 $μ$Jy/beam for the full area. We detect 112 Virgo member galaxies and 114 background galaxies. In at least 18 cases, the radio morphology of the cluster member galaxies shows clear signs of ram-pressure stripping. This includes three previously unreported candidates. In addition, we reveal for the first time 150 kpc long tails from a previous epoch of AGN activity for NGC 4472 (M 49). While no cluster-scale diffuse radio sources are discovered, we find the presence of an extended radio signature of the W$'$-group. △ Less

Submitted 7 June, 2023; originally announced June 2023.

Comments: 18 pages, 12 figures, accepted for publication in A&A

Journal ref: A&A 676, A24 (2023)

The LOFAR Two-metre Sky Survey: the radio view of the cosmic star formation history

Authors: R. K. Cochrane, R. Kondapally, P. N. Best, J. Sabater, K. J. Duncan, D. J. B. Smith, M. J. Hardcastle, H. J. A. Röttgering, I. Prandoni, P. Haskell, G. Gürkan, G. K. Miley

Abstract: We present a detailed study of the cosmic star formation history over $90$ per cent of cosmic time ($0\lesssim z\lesssim4$), using deep, radio continuum observations that probe star formation activity independent of dust. The Low Frequency Array Two Metre Sky Survey has imaged three well-studied extragalactic fields, Elais-N1, Boötes and the Lockman Hole, reaching $\sim20\,μ\rm{Jy/beam}$ rms sensi… ▽ More We present a detailed study of the cosmic star formation history over $90$ per cent of cosmic time ($0\lesssim z\lesssim4$), using deep, radio continuum observations that probe star formation activity independent of dust. The Low Frequency Array Two Metre Sky Survey has imaged three well-studied extragalactic fields, Elais-N1, Boötes and the Lockman Hole, reaching $\sim20\,μ\rm{Jy/beam}$ rms sensitivity at $150\,\rm{MHz}$. The availability of high-quality ancillary data from ultraviolet to far-infrared wavelengths has enabled accurate photometric redshifts and the robust separation of radio-bright AGN from their star-forming counterparts. We capitalise on this unique combination of deep, wide fields and robustly-selected star-forming galaxies to construct radio luminosity functions and derive the cosmic star formation rate density. We carefully constrain and correct for scatter in the $L_{150\,\rm{MHz}}-\rm{SFR}$ relation, which we find to be $\sim0.3\,\rm{dex}$. Our derived star formation rate density lies between previous measurements at all redshifts studied. We derive higher star formation rate densities between $z\sim0$ and $z\sim3$ than are typically inferred from short wavelength emission; at earlier times, this discrepancy is reduced. Our measurements are generally in good agreement with far-infrared and radio-based studies, with small offsets resulting from differing star formation rate calibrations. △ Less

Submitted 24 May, 2023; originally announced May 2023.

Comments: 22 pages, 14 figures; accepted for publication in MNRAS

Spectral age distribution for radio-loud active galaxies in the XMM-LSS field

Authors: Siddhant Pinjarkar, Martin J. Hardcastle, Jeremy J. Harwood, Dharam V. Lal, Peter W. Hatfield, Matt J. Jarvis, Zara Randriamanakoto, Imogen H. Whittam

Abstract: Jets of energetic particles, as seen in FR type-I and FR type-II sources, ejected from the center of Radio-Loud AGN affect the sources surrounding intracluster medium/intergalactic medium. Placing constraints on the age of such sources is important in order to measure the jet powers and determine the effects on feedback. To evaluate the age of these sources using spectral age models, we require hi… ▽ More Jets of energetic particles, as seen in FR type-I and FR type-II sources, ejected from the center of Radio-Loud AGN affect the sources surrounding intracluster medium/intergalactic medium. Placing constraints on the age of such sources is important in order to measure the jet powers and determine the effects on feedback. To evaluate the age of these sources using spectral age models, we require high-resolution multi-wavelength data. The new sensitive and high-resolution MIGHTEE survey of the XMM-LSS field along with data from the Low Frequency Array (LOFAR) and the Giant Metrewave Radio Telescope (GMRT) provide data taken at different frequencies with similar resolution, which enables us to determine the spectral age distribution for radio loud AGN in the survey field. In this study we present a sample of 28 radio galaxies with their best fitting spectral age distribution analyzed using the Jaffe-Perola (JP) model on a pixel-by-pixel basis. Fits are generally good and objects in our sample show maximum ages within the range of 2.8 Myr to 115 Myr with a median of 8.71 Myr. High-resolution maps over a range of frequencies are required to observe detailed age distributions for small sources and high-sensitivity maps will be needed in order to observe fainter extended emission. We do not observe any correlation between the total physical size of the sources and their age and we speculate both dynamical models and the approach to spectral age analysis may need some modification to account for our observations. △ Less

Submitted 12 May, 2023; originally announced May 2023.

Comments: 20 pages, 9 figures

The LOFAR Two-metre Sky Survey Deep Fields Data Release 1: V. Survey description, source classifications and host galaxy properties

Authors: P. N. Best, R. Kondapally, W. L. Williams, R. K. Cochrane, K. J. Duncan, C. L. Hale, P. Haskell, K. Malek, I. McCheyne, D. J. B. Smith, L. Wang, A. Botteon, M. Bonato, M. Bondi, G. Calistro Rivera, F. Gao, G. Gurkan, M. J. Hardcastle, M. J. Jarvis, B. Mingo, H. Miraghaei, L. K. Morabito, D. Nisbet, I. Prandoni, H. J. A. Rottgering , et al. (4 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 9 May, 2023; originally announced May 2023.

Comments: Accepted for publication in MNRAS. Catalogues available at www.lofar-surveys.org/deepfields.html. 27 pages

Finding AGN remnant candidates based on radio morphology with machine learning

Authors: Rafael I. J. Mostert, Raffaella Morganti, Marisa Brienza, Kenneth J. Duncan, Martijn S. S. L. Oei, Huub J. A. Rottgering, Lara Alegre, Martin J. Hardcastle, Nika Jurlin

Abstract: Remnant radio galaxies represent the dying phase of radio-loud active galactic nuclei (AGN). Large samples of remnant radio galaxies are important for quantifying the radio galaxy life cycle. The remnants of radio-loud AGN can be identified in radio sky surveys based on their spectral index, or, complementary, through visual inspection based on their radio morphology. However, this is extremely ti… ▽ More Remnant radio galaxies represent the dying phase of radio-loud active galactic nuclei (AGN). Large samples of remnant radio galaxies are important for quantifying the radio galaxy life cycle. The remnants of radio-loud AGN can be identified in radio sky surveys based on their spectral index, or, complementary, through visual inspection based on their radio morphology. However, this is extremely time-consuming when applied to the new large and sensitive radio surveys. Here we aim to reduce the amount of visual inspection required to find AGN remnants based on their morphology, through supervised machine learning trained on an existing sample of remnant candidates. For a dataset of 4107 radio sources, with angular sizes larger than 60 arcsec, from the LOw Frequency ARray (LOFAR) Two-Metre Sky Survey second data release (LoTSS-DR2), we started with 151 radio sources that were visually classified as 'AGN remnant candidate'. We derived a wide range of morphological features for all radio sources from their corresponding Stokes-I images: from simple source catalogue-derived properties, to clustered Haralick-features, and self-organising map (SOM) derived morphological features. We trained a random forest classifier to separate the 'AGN remnant candidates' from the not yet inspected sources. The SOM-derived features and the total to peak flux ratio of a source are shown to be most salient to the classifier. We estimate that $31\pm5\%$ of sources with positive predictions from our classifier will be labelled 'AGN remnant candidates' upon visual inspection, while we estimate the upper bound of the $95\%$ confidence interval for 'AGN remnant candidates' in the negative predictions at $8\%$. Visual inspection of just the positive predictions reduces the number of radio sources requiring visual inspection by $73\%$. △ Less

Submitted 12 April, 2023; originally announced April 2023.

Comments: 23 pages; accepted for publication in A&A

Journal ref: A&A 674, A208 (2023)

SKA Science Data Challenge 2: analysis and results

Authors: P. Hartley, A. Bonaldi, R. Braun, J. N. H. S. Aditya, S. Aicardi, L. Alegre, A. Chakraborty, X. Chen, S. Choudhuri, A. O. Clarke, J. Coles, J. S. Collinson, D. Cornu, L. Darriba, M. Delli Veneri, J. Forbrich, B. Fraga, A. Galan, J. Garrido, F. Gubanov, H. Håkansson, M. J. Hardcastle, C. Heneka, D. Herranz, K. M. Hess , et al. (83 additional authors not shown)

Abstract: 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 t… ▽ More 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. △ Less

Submitted 14 March, 2023; originally announced March 2023.

Comments: Under review by MNRAS; 28 pages, 16 figures

New mechanisms for forming multiple hotspots in radio jets

Authors: Maya A. Horton, Martin G. H. Krause, Martin J. Hardcastle

Abstract: Hotspots of radio galaxies are regions of shock-driven particle acceleration. Multiple hotspots have long been identified as potential indicators of jet movement or precession. Two frequent explanations describe a secondary hotspot as either the location of a prior jet termination point, or a deflected backflow-driven shock: the so-called Dentist's Drill and Splatter Spot models. We created high-r… ▽ More Hotspots of radio galaxies are regions of shock-driven particle acceleration. Multiple hotspots have long been identified as potential indicators of jet movement or precession. Two frequent explanations describe a secondary hotspot as either the location of a prior jet termination point, or a deflected backflow-driven shock: the so-called Dentist's Drill and Splatter Spot models. We created high-resolution simulations of precessing jets with a range of parameters. In addition to the existing mechanisms, our results show three additional mechanisms for multiple hotspot formation: (1) the splitting of a large terminal hotspots into passive and active components; (2) jet stream splitting resulting in two active hotspots; (3) dynamic multiple hotspot complexes that form as a result of jet termination in a turbulent cocoon, linked here to rapid precession. We show that these distinct types of multiple hotspots are difficult to differentiate in synthetic radio maps, particularly hotspot complexes which can easily be mistaken for the jet itself. We discuss the implication for hypothesised binary supermassive black hole systems where jet precession is a key component of the morphology, and show a selection of potential precession candidates found using the LOFAR Two-Metre Sky Survey Data Release 2 (LoTSS DR2). △ Less

Submitted 27 February, 2023; originally announced February 2023.

Comments: 14 pages, 10 figures

A low frequency sub-arcsecond view of powerful radio galaxies in rich-cluster environments: 3C 34 and 3C 320

Authors: V. H. Mahatma, A. Basu, M. J. Hardcastle, L. K. Morabito, R. J. van Weeren

Abstract: Models of radio galaxy physics have been primarily based on high frequency ($\geqslant$1 GHz) observations of their jets, hotspots, and lobes. Without highly resolved low frequency observations, which provide information on older plasma, our understanding of the dynamics of radio galaxies and their interaction with their environment is limited. Here, we present the first sub-arcsecond (0.3") resol… ▽ More Models of radio galaxy physics have been primarily based on high frequency ($\geqslant$1 GHz) observations of their jets, hotspots, and lobes. Without highly resolved low frequency observations, which provide information on older plasma, our understanding of the dynamics of radio galaxies and their interaction with their environment is limited. Here, we present the first sub-arcsecond (0.3") resolution images at 144 MHz of two powerful radio galaxies situated in rich cluster environments, namely 3C 34 and 3C 320, using the International Low Frequency Array Telescope. We detect for the first time at low frequencies a plethora of structures in these objects, including strikingly large filaments across the base of the lobes in both sources, which are spatially associated with dense regions in the ambient medium. For 3C 34 we report a spectral flattening in the region of the central filament, suggesting that the origin of the filaments is related to the presence of large-scale ordered magnetic fields. We also report periodic total intensity and spectral index banding of diffuse emission in the eastern lobe, seen for the first time in radio galaxy lobes. The hotspot complexes are resolved into multiple fragments of varying structure and spectral index; we discuss the implications for particle acceleration and jet termination models. We find at most smooth gradients in the spectral behaviour of the hotspot structure suggesting that particle acceleration, if present, may be occurring throughout the complex, in contrast to simple models, but different jet termination models may apply to both sources. △ Less

Submitted 2 February, 2023; originally announced February 2023.

Comments: 16 pages, 7 figures, 2 tables. Accepted for publication in Monthly Notices of the Royal Astronomical Society, 2023 February 01

The LOFAR LBA Sky Survey II. First data release

Authors: F. de Gasperin, H. W. Edler, W. L. Williams, J. R. Callingham, B. Asabere, M. Bruggen, G. Brunetti, T. J. Dijkema, M. J. Hardcastle, M. Iacobelli, A. Offringa, M. J. Norden, H. J. A. Rottgering, T. Shimwell, R. J. van Weeren, C. Tasse, D. J. Bomans, A. Bonafede, A. Botteon, R. Cassano, K. T. Chyzy, V. Cuciti, K. L. Emig, M. Kadler, G. Miley , et al. (5 additional authors not shown)

Abstract: 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). LoLS… ▽ More 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. △ Less

Submitted 30 January, 2023; originally announced January 2023.

Comments: 20 pages, 22 figures, images and catalogues available at https://www.lofar-surveys.org/lolss.html

Journal ref: A&A 673, A165 (2023)

The Faraday Rotation Measure Grid of the LOFAR Two-metre Sky Survey: Data Release 2

Authors: S. P. O'Sullivan, T. W. Shimwell, M. J. Hardcastle, C. Tasse, G. Heald, E. Carretti, M. Brüggen, V. Vacca, C. Sobey, C. L. Van Eck, C. Horellou, R. Beck, M. Bilicki, S. Bourke, A. Botteon, J. H. Croston, A. Drabent, K. Duncan, V. Heesen, S. Ideguchi, M. Kirwan, L. Lawlor, B. Mingo, B. Nikiel-Wroczyński, J. Piotrowska , et al. (2 additional authors not shown)

Abstract: 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 polar… ▽ More 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σ_{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. △ Less

Submitted 18 January, 2023; originally announced January 2023.

Comments: 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

CosmoDRAGoN simulations -- I. Dynamics and observable signatures of radio jets in cosmological environments

Authors: Patrick M. Yates-Jones, Stanislav S. Shabala, Chris Power, Martin G. H. Krause, Martin J. Hardcastle, Elena A. N. Mohd Noh Velastín, Georgia S. C. Stewart

Abstract: We present the Cosmological Double Radio Active Galactic Nuclei (CosmoDRAGoN) project: a large suite of simulated AGN jets in cosmological environments. These environments sample the intra-cluster media of galaxy clusters that form in cosmological smooth particle hydrodynamics (SPH) simulations, which we then use as inputs for grid-based hydrodynamic simulations of radio jets. Initially conical je… ▽ More We present the Cosmological Double Radio Active Galactic Nuclei (CosmoDRAGoN) project: a large suite of simulated AGN jets in cosmological environments. These environments sample the intra-cluster media of galaxy clusters that form in cosmological smooth particle hydrodynamics (SPH) simulations, which we then use as inputs for grid-based hydrodynamic simulations of radio jets. Initially conical jets are injected with a range of jet powers, speeds (both relativistic and non-relativistic), and opening angles; we follow their collimation and propagation on scales of tens to hundreds of kiloparsecs, and calculate spatially-resolved synthetic radio spectra in post-processing. In this paper, we present a technical overview of the project, and key early science results from six representative simulations which produce radio sources with both core- (Fanaroff-Riley Type I) and edge-brightened (Fanaroff-Riley Type II) radio morphologies. Our simulations highlight the importance of accurate representation of both jets and environments for radio morphology, radio spectra, and feedback the jets provide to their surroundings. △ Less

Submitted 28 February, 2023; v1 submitted 20 December, 2022; originally announced December 2022.

Comments: 25 pages, 20 figures. Accepted for publication in PASA

V-LoTSS: The Circularly-Polarised LOFAR Two-metre Sky Survey

Authors: J. R. Callingham, T. W. Shimwell, H. K. Vedantham, C. G. Bassa, S. P. O'Sullivan, T. W. H. Yiu, S. Bloot, P. N. Best, M. J. Hardcastle, M. Haverkorn, R. D. Kavanagh, L. Lamy, B. J. S. Pope, H. J. A. Röttgering, D. J. Schwarz, C. Tasse, R. J. van Weeren, G. J. White, P. Zarka, D. J. Bomans, A. Bonafede, M. Bonato, A. Botteon, M. Bruggen, K. T. Chyży , et al. (22 additional authors not shown)

Abstract: 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 $μ$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 polar… ▽ More 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 $μ$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. △ Less

Submitted 19 December, 2022; originally announced December 2022.

Comments: 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

Journal ref: A&A 670, A124 (2023)

The wide-field, multiplexed, spectroscopic facility WEAVE: Survey design, overview, and simulated implementation

Authors: Shoko Jin, Scott C. Trager, Gavin B. Dalton, J. Alfonso L. Aguerri, J. E. Drew, Jesús Falcón-Barroso, Boris T. Gänsicke, Vanessa Hill, Angela Iovino, Matthew M. Pieri, Bianca M. Poggianti, D. J. B. Smith, Antonella Vallenari, Don Carlos Abrams, David S. Aguado, Teresa Antoja, Alfonso Aragón-Salamanca, Yago Ascasibar, Carine Babusiaux, Marc Balcells, R. Barrena, Giuseppina Battaglia, Vasily Belokurov, Thomas Bensby, Piercarlo Bonifacio , et al. (190 additional authors not shown)

Abstract: WEAVE, the new wide-field, massively multiplexed spectroscopic survey facility for the William Herschel Telescope, will see first light in late 2022. WEAVE comprises a new 2-degree field-of-view prime-focus corrector system, a nearly 1000-multiplex fibre positioner, 20 individually deployable 'mini' integral field units (IFUs), and a single large IFU. These fibre systems feed a dual-beam spectrogr… ▽ More WEAVE, the new wide-field, massively multiplexed spectroscopic survey facility for the William Herschel Telescope, will see first light in late 2022. WEAVE comprises a new 2-degree field-of-view prime-focus corrector system, a nearly 1000-multiplex fibre positioner, 20 individually deployable 'mini' integral field units (IFUs), and a single large IFU. These fibre systems feed a dual-beam spectrograph covering the wavelength range 366$-$959\,nm at $R\sim5000$, or two shorter ranges at $R\sim20\,000$. After summarising the design and implementation of WEAVE and its data systems, we present the organisation, science drivers and design of a five- to seven-year programme of eight individual surveys to: (i) study our Galaxy's origins by completing Gaia's phase-space information, providing metallicities to its limiting magnitude for $\sim$3 million stars and detailed abundances for $\sim1.5$ million brighter field and open-cluster stars; (ii) survey $\sim0.4$ million Galactic-plane OBA stars, young stellar objects and nearby gas to understand the evolution of young stars and their environments; (iii) perform an extensive spectral survey of white dwarfs; (iv) survey $\sim400$ neutral-hydrogen-selected galaxies with the IFUs; (v) study properties and kinematics of stellar populations and ionised gas in $z<0.5$ cluster galaxies; (vi) survey stellar populations and kinematics in $\sim25\,000$ field galaxies at $0.3\lesssim z \lesssim 0.7$; (vii) study the cosmic evolution of accretion and star formation using $>1$ million spectra of LOFAR-selected radio sources; (viii) trace structures using intergalactic/circumgalactic gas at $z>2$. Finally, we describe the WEAVE Operational Rehearsals using the WEAVE Simulator. △ Less

Submitted 31 October, 2023; v1 submitted 7 December, 2022; originally announced December 2022.

Comments: 41 pages, 27 figures, accepted for publication by MNRAS; updated version including information on individual grants in a revised Acknowledgements section, corrections to the affiliation list, and an updated references list

Extragalactic Peaked-Spectrum Radio Sources at Low-Frequencies are Young Radio Galaxies

Authors: M. M. Slob, J. R. Callingham, H. J. A. Röttgering, W. L. Williams, K. J. Duncan, F. de Gasperin, M. J. Hardcastle, G. K. Miley

Abstract: We present a sample of 373 peaked-spectrum (PS) sources with spectral peaks around 150MHz, selected using a subset of two LOFAR all-sky surveys, the LOFAR Two Meter Sky Survey and the LOFAR LBA Sky Survey. These surveys are the most sensitive low-frequency widefield surveys to date, allowing us to select low-luminosity PS sources. Our sample increases the number of known PS sources in our survey a… ▽ More We present a sample of 373 peaked-spectrum (PS) sources with spectral peaks around 150MHz, selected using a subset of two LOFAR all-sky surveys, the LOFAR Two Meter Sky Survey and the LOFAR LBA Sky Survey. These surveys are the most sensitive low-frequency widefield surveys to date, allowing us to select low-luminosity PS sources. Our sample increases the number of known PS sources in our survey area by a factor 50. The 5GHz luminosity distribution of our PS sample shows we sample the lowest luminosity PS sources to-date by nearly an order of magnitude. Since high-frequency PS sources and compact steep-spectrum sources are hypothesised to be the precursors to large radio galaxies, we investigate whether this is also the case for our sample of low-frequency PS sources. Using optical line emission criteria, we find that our PS sources are predominately high-excitation radio galaxies instead of low-excitation radio galaxies, corresponding to a quickly evolving population. We compute the radio source counts of our PS sample, and find they are scaled down by a factor of $\sim$40 compared to a general sample of radio-loud active galactic nuclei (AGN). This implies that the lifetimes of PS sources are 40 times shorter than large scale radio galaxies, if their luminosity functions are identical. To investigate this, we compute the first radio luminosity function for a homogeneously-selected PS sample. We find that for 144MHz luminosities $\gtrsim 10^{25}$W Hz$^{-1}$, the PS luminosity function has the same shape as an unresolved radio-loud AGN population but shifted down by a factor of $\sim$10. We interpret this as strong evidence that these high-luminosity PS sources evolve into large-scale radio-loud AGN. For local, low-luminosity PS sources, there is a surplus of PS sources, which we hypothesise to be the addition of frustrated PS sources that do not evolve into large-scale AGN. △ Less

Submitted 29 October, 2022; originally announced October 2022.

Comments: 21 pages, accepted in Astronomy & Astrophysics

Journal ref: A&A 668, A186 (2022)

Measuring the giant radio galaxy length distribution with the LoTSS

Authors: Martijn S. S. L. Oei, Reinout J. van Weeren, Aivin R. D. J. G. I. B. Gast, Andrea Botteon, Martin J. Hardcastle, Pratik Dabhade, Tim W. Shimwell, Huub J. A. Röttgering, Alexander Drabent

Abstract: Radio galaxies are luminous structures created by the jets of supermassive black holes, and consist of atomic nuclei, relativistic electrons, and magnetic fields. In exceptional cases, radio galaxies attain cosmological, megaparsec extents - and thus turn into giants. Giants embody the most extreme known mechanism through which galaxies can impact the Cosmic Web around them. The triggers of giant… ▽ More Radio galaxies are luminous structures created by the jets of supermassive black holes, and consist of atomic nuclei, relativistic electrons, and magnetic fields. In exceptional cases, radio galaxies attain cosmological, megaparsec extents - and thus turn into giants. Giants embody the most extreme known mechanism through which galaxies can impact the Cosmic Web around them. The triggers of giant growth remain a mystery. Excitingly, new sensitive low-frequency sky surveys hold promise to change this situation. In this work, we perform a precision measurement of the distribution of giant growth's central dynamical quantity: total length. We first construct a statistical geometric framework for radio galaxies that is both rigorous and practical. We then search the LOFAR Two-metre Sky Survey for giants, discovering 2060 previously unknown specimina: more than have been found in all preceding literature combined. Spectacular discoveries include the longest giant hosted by an elliptical galaxy, the longest giant hosted by a spiral galaxy, and 13 giants with an angular length larger than that of the full Moon. By combining theory and observations - and carefully forward modelling selection effects - we infer that giant radio galaxy lengths are well described by a Pareto distribution with tail index $-3.5 \pm 0.5$. This finding is a new observational constraint for models and simulations of radio galaxy growth. In addition, for the first time, we determine the comoving number density of giants, $5 \pm 2\ (100\ \mathrm{Mpc})^{-3}$, and the volume-filling fraction of giant radio galaxy lobes in clusters and filaments, $5\substack{+8\\-2}\cdot 10^{-6}$. We conclude that giants are truly rare - not only from an observational perspective, but also from a cosmological one. At any moment in time, most clusters and filaments - the building blocks of the modern Cosmic Web - do not harbour giants. △ Less

Submitted 1 February, 2023; v1 submitted 18 October, 2022; originally announced October 2022.

Comments: 41 pages, 28 figures, 5 tables, published in Astronomy & Astrophysics

Journal ref: A&A 672, A163 (2023)

An intergalactic medium temperature from a giant radio galaxy

Authors: Martijn S. S. L. Oei, Reinout J. van Weeren, Martin J. Hardcastle, Franco Vazza, Tim W. Shimwell, Florent Leclercq, Marcus Brüggen, Huub J. A. Röttgering

Abstract: The warm-hot intergalactic medium (warm-hot IGM, or WHIM) pervades the filaments of the Cosmic Web and harbours half of the Universe's baryons. The WHIM's thermodynamic properties are notoriously hard to measure. Here we estimate a galaxy group - WHIM boundary temperature using a new method. In particular, we use a radio image of the giant radio galaxy (giant RG, or GRG) created by NGC 6185, a mas… ▽ More The warm-hot intergalactic medium (warm-hot IGM, or WHIM) pervades the filaments of the Cosmic Web and harbours half of the Universe's baryons. The WHIM's thermodynamic properties are notoriously hard to measure. Here we estimate a galaxy group - WHIM boundary temperature using a new method. In particular, we use a radio image of the giant radio galaxy (giant RG, or GRG) created by NGC 6185, a massive nearby spiral. We analyse this extraordinary object with a Bayesian 3D lobe model and deduce an equipartition pressure $P_\mathrm{eq} = 6 \cdot 10^{-16}\ \mathrm{Pa}$ -- among the lowest found in RGs yet. Using an X-ray-based statistical conversion for Fanaroff-Riley II RGs, we find a true lobe pressure $P = 1.5\substack{+1.7\\-0.4}\cdot 10^{-15}\ \mathrm{Pa}$. Cosmic Web reconstructions, group catalogues, and MHD simulations furthermore imply an $\mathrm{Mpc}$-scale IGM density $1 + δ_\mathrm{IGM} = 40\substack{+30\\-10}$. The buoyantly rising lobes are crushed by the IGM at their inner side, where an approximate balance between IGM and lobe pressure occurs: $P_\mathrm{IGM} \approx P$. The ideal gas law then suggests an IGM temperature $T_\mathrm{IGM} = 11\substack{+12\\-5} \cdot 10^6\ \mathrm{K}$, or $k_\mathrm{B}T_\mathrm{IGM} = 0.9\substack{+1.0\\-0.4}\ \mathrm{keV}$, at the virial radius -- consistent with X-ray-derived temperatures of similarly massive groups. Interestingly, the method is not performing at its limit: in principle, estimates $T_\mathrm{IGM} \sim 4 \cdot 10^6\ \mathrm{K}$ are already possible -- rivalling the lowest X-ray measurements available. The technique's future scope extends from galaxy group outskirts to the WHIM. In conclusion, we demonstrate that observations of GRGs in Cosmic Web filaments are finally sensitive enough to probe the thermodynamics of galaxy groups and beyond. △ Less

Submitted 18 October, 2022; originally announced October 2022.

Comments: 17 pages, 14 figures, 5 tables, accepted for publication in Monthly Notices of the Royal Astronomical Society

Discovery of 24 radio-bright quasars at $4.9 \leq z \leq6.6$ using low-frequency radio observations

Authors: A. J. Gloudemans, K. J. Duncan, A. Saxena, Y. Harikane, G. J. Hill, G. R. Zeimann, H. J. A. Rottgering, D. Yang, P. N. Best, E. Banados, A. Drabent, M. J. Hardcastle, J. F. Hennawi, G. Lansbury, M. Magliocchetti, G. K. Miley, R. Nanni, T. W. Shimwell, D. J. B. Smith, B. P. Venemans, J. D. Wagenveld

Abstract: High redshift quasars ($z>5$) that also shine brightly at radio wavelengths are unique signposts of supermassive black hole activity in the early universe. However, bright radio sources at $z\ge5$ are extremely rare and therefore we have started a campaign to search for new high-$z$ quasars by combining an optical dropout selection driven by the $g$, $r$, and $z$ bands from the Dark Energy Spectro… ▽ More High redshift quasars ($z>5$) that also shine brightly at radio wavelengths are unique signposts of supermassive black hole activity in the early universe. However, bright radio sources at $z\ge5$ are extremely rare and therefore we have started a campaign to search for new high-$z$ quasars by combining an optical dropout selection driven by the $g$, $r$, and $z$ bands from the Dark Energy Spectroscopic Instrument (DESI) Legacy Imaging Surveys with low-frequency radio observations from the LOFAR Two-metre Sky Survey (LoTSS). Currently, LoTSS covers a large fraction of the northern sky (5720 deg$^2$) to such a depth (median noise level of 83 $μ$Jy beam$^{-1}$) that about 30% of the general quasar population is detected $-$ which is a factor of 5-10 more than previous large sky radio surveys such as NVSS and FIRST, respectively. In this paper, we present the discovery of 20 new quasars (and the independent confirmation of 4) between $4.9\leq z\leq 6.6$. Out of the 24 quasars, 21 satisfy the traditional radio-loudness criterion of $R=f_{5\text{GHz}}/f_{4400A} > 10$, with the full sample spanning $R\sim$6-1000, thereby more than doubling the sample of known radio-loud quasars at $z \ge 5$. Our radio detection requirement strongly decreases the contamination of stellar sources and allows one to select these quasars in a broad redshift range. Despite selecting our quasar candidates using fewer and less conservative colour restrictions, both the optical and near-infrared colours, Ly$α$ emission line properties, and dust reddening, $E(B-V)$, measurements of our quasar sample do not deviate from the known radio-quiet quasar population, suggesting similar optical quasar properties of the radio-loud and radio-quiet quasar population at high-$z$. Our campaign demonstrates the potential for discovering new high-$z$ quasar populations through next generation radio continuum surveys. △ Less

Submitted 21 October, 2022; v1 submitted 4 October, 2022; originally announced October 2022.

Comments: 19 pages, 11 figures, Tab A1 corrected, accepted for publication in A&A

Radio source-component association for the LOFAR Two-metre Sky Survey with region-based convolutional neural networks

Authors: Rafaël I. J. Mostert, Kenneth J. Duncan, Lara Alegre, Huub J. A. Röttgering, Wendy L. Williams, Philip N. Best, Martin J. Hardcastle, Raffaella Morganti

Abstract: Radio loud active galactic nuclei (RLAGNs) are often morphologically complex objects that can consist of multiple, spatially separated, components. Astronomers often rely on visual inspection to resolve radio component association. However, applying visual inspection to all the hundreds of thousands of well-resolved RLAGNs that appear in the images from the Low Frequency Array (LOFAR) Two-metre Sk… ▽ More Radio loud active galactic nuclei (RLAGNs) are often morphologically complex objects that can consist of multiple, spatially separated, components. Astronomers often rely on visual inspection to resolve radio component association. However, applying visual inspection to all the hundreds of thousands of well-resolved RLAGNs that appear in the images from the Low Frequency Array (LOFAR) Two-metre Sky Survey (LoTSS) at $144$ MHz, is a daunting, time-consuming process, even with extensive manpower. Using a machine learning approach, we aim to automate the radio component association of large ($> 15$ arcsec) radio components. We turned the association problem into a classification problem and trained an adapted Fast region-based convolutional neural network to mimic the expert annotations from the first LoTSS data release. We implemented a rotation data augmentation to reduce overfitting and simplify the component association by removing unresolved radio sources that are likely unrelated to the large and bright radio components that we consider using predictions from an existing gradient boosting classifier. For large ($> 15$ arcsec) and bright ($> 10$ mJy) radio components in the LoTSS first data release, our model provides the same associations for $85.3\%\pm0.6$ of the cases as those derived when astronomers perform the association manually. When the association is done through public crowd-sourced efforts, a result similar to that of our model is attained. Our method is able to efficiently carry out manual radio-component association for huge radio surveys and can serve as a basis for either automated radio morphology classification or automated optical host identification. This opens up an avenue to study the completeness and reliability of samples of radio sources with extended, complex morphologies. △ Less

Submitted 28 September, 2022; originally announced September 2022.

Comments: 22 pages; accepted for publication in A&A

Journal ref: A&A 668, A28 (2022)

A numerical study of the interplay between Fermi acceleration mechanisms in radio lobes of FR-II radio galaxies

Authors: Sayan Kundu, Bhargav Vaidya, Andrea Mignone, Martin J. Hardcastle

Abstract: Context: Radio-loud AGNs are thought to possess various sites of particle acceleration, which gives rise to the observed non-thermal spectra. Stochastic turbulent acceleration (STA) and diffusive shock acceleration (DSA) are commonly cited as potential sources of high-energy particles in weakly magnetized environments. Together, these acceleration processes and various radiative losses determine t… ▽ More Context: Radio-loud AGNs are thought to possess various sites of particle acceleration, which gives rise to the observed non-thermal spectra. Stochastic turbulent acceleration (STA) and diffusive shock acceleration (DSA) are commonly cited as potential sources of high-energy particles in weakly magnetized environments. Together, these acceleration processes and various radiative losses determine the emission characteristics of these extra-galactic radio sources. Aims: The purpose of this research is to investigate the dynamical interplay between the STA and DSA in the radio lobes of FR-II radio galaxies, as well as the manner in which these acceleration mechanisms, along with a variety of radiative losses, collectively shape the emission features seen in these extra-galactic sources. Methods: A phenomenologically motivated model of STA is considered and subsequently employed on a magneto-hydrodynamically simulated radio lobe through a novel hybrid Eulerian-Lagrangian framework. Results: STA gives rise to a curved particle spectrum that is morphologically different from the usual shock-accelerated spectrum. As a consequence of this structural difference in the underlying particle energy spectrum, various multi-wavelength features arise in the spectral energy distribution of the radio lobe. Additionally, we observe enhanced diffuse X-ray emission from radio lobes for cases where STA is taken into account in addition to DSA. △ Less

Submitted 26 August, 2022; originally announced August 2022.

Comments: 15 pages, 10 figures, Accepted for publication in A&A. Comments are welcome

Journal ref: A&A 667, A138 (2022)

Possible discovery of Calvera's supernova remnant

Authors: M. Arias, A. Botteon, C. G. Bassa, S. van der Jagt, R. J. van Weeren, S. P. O'Sullivan, Q. Bosschaart, R. S. Dullaart, M. J. Hardcastle, J. W. T. Hessels, T. Shimwell, M. M. Slob, J. A. Sturm, C. Tasse, N. C. M. A. Theijssen, J. Vink

Abstract: 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 $α=14\mathrm{h}11\mathrm{m}12.6\mathrm{s}$, $δ=+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 r… ▽ More 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 $α=14\mathrm{h}11\mathrm{m}12.6\mathrm{s}$, $δ=+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$α$ 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. △ Less

Submitted 28 July, 2022; originally announced July 2022.

Comments: Accepted in A&A

A machine learning classifier for LOFAR radio galaxy cross-matching techniques

Authors: Lara Alegre, Jose Sabater, Philip Best, Rafaël I. J. Mostert, Wendy L. Williams, Gülay Gürkan, Martin J. Hardcastle, Rohit Kondapally, Tim W. Shimwell, Daniel J. B. Smith

Abstract: New-generation radio telescopes like LOFAR are conducting extensive sky surveys, detecting millions of sources. To maximise the scientific value of these surveys, radio source components must be properly associated into physical sources before being cross-matched with their optical/infrared counterparts. In this paper, we use machine learning to identify those radio sources for which either source… ▽ More New-generation radio telescopes like LOFAR are conducting extensive sky surveys, detecting millions of sources. To maximise the scientific value of these surveys, radio source components must be properly associated into physical sources before being cross-matched with their optical/infrared counterparts. In this paper, we use machine learning to identify those radio sources for which either source association is required or statistical cross-matching to optical/infrared catalogues is unreliable. We train a binary classifier using manual annotations from the LOFAR Two-metre Sky Survey (LoTSS). We find that, compared to a classification model based on just the radio source parameters, the addition of features of the nearest-neighbour radio sources, the potential optical host galaxy, and the radio source composition in terms of Gaussian components, all improve model performance. Our best model, a gradient boosting classifier, achieves an accuracy of 95 per cent on a balanced dataset and 96 per cent on the whole (unbalanced) sample after optimising the classification threshold. Unsurprisingly, the classifier performs best on small, unresolved radio sources, reaching almost 99 per cent accuracy for sources smaller than 15 arcsec, but still achieves 70 per cent accuracy on resolved sources. It flags 68 per cent more sources than required as needing visual inspection, but this is still fewer than the manually-developed decision tree used in LoTSS, while also having a lower rate of wrongly accepted sources for statistical analysis. The results have an immediate practical application for cross-matching the next LoTSS data releases and can be generalised to other radio surveys. △ Less

Submitted 4 July, 2022; originally announced July 2022.

Comments: Accepted for publication in MNRAS. 24 pages, 13 Figures. Submission includes a data table which is also available at https://drive.google.com/file/d/1mwKWJW6oCqh-joNR1uV3PPPYDdYz0n_J/view?usp=sharing

Relic jet activity in "Hanny's Voorwerp" revealed by the LOFAR Two metre Sky Survey

Authors: D. J. B. Smith, M. G. Krause, M. J. Hardcastle, A. B. Drake

Abstract: We report new observations of "Hanny's Voorwerp" (hereafter HV) taken from the second data release of the LOFAR Two-metre Sky Survey (LoTSS). HV is a highly-ionised region in the environs of the galaxy IC2497, first discovered by the Galaxy Zoo project. The new 150MHz observations are considered in the context of existing multi-frequency radio data and archival narrow-band imaging from the Hubble… ▽ More We report new observations of "Hanny's Voorwerp" (hereafter HV) taken from the second data release of the LOFAR Two-metre Sky Survey (LoTSS). HV is a highly-ionised region in the environs of the galaxy IC2497, first discovered by the Galaxy Zoo project. The new 150MHz observations are considered in the context of existing multi-frequency radio data and archival narrow-band imaging from the Hubble Space Telescope, centred on the [Oiii] emission line. The combined sensitivity and spatial resolution of the LoTSS data -- which far exceed what was previously available at radio frequencies -- reveal clear evidence for large-scale extended emission emanating from the nucleus of IC2497. The radio jet appears to have punched a hole in the neutral gas halo, in a region co-located with HV. The new 150MHz data, alongside newly-processed archival 1.64GHz eVLA data, reveal that the extended emission has a steep spectrum, implying an age $>10^8$yr. The jet supplying the extended 150MHz structure must have "turned off" long before the change in X-ray luminosity reported in recent works. In this picture, a combination of jet activity and the influence of the radiatively efficient active galactic nucleus are responsible for the unusual appearance of HV. △ Less

Submitted 13 June, 2022; originally announced June 2022.

Comments: 7 pages, 4 figures, accepted for publication in MNRAS

Star Formation and AGN Feedback in the Local Universe: Combining LOFAR and MaNGA

Authors: C. R. Mulcahey, S. K. Leslie, T. M. Jackson, J. E. Young, I. Prandoni, M. J. Hardcastle, N. Roy, K. Małek, M. Magliocchetti, M. Bonato, H. J. A. Röttgering, A. Drabent

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

Submitted 2 June, 2022; originally announced June 2022.

Comments: Accepted for publication in A&A. 24 pages, 14 figures

Journal ref: A&A 665, A144 (2022)

Cosmic evolution of low-excitation radio galaxies in the LOFAR Two-meter Sky Survey Deep Fields

Authors: R. Kondapally, P. N. Best, R. K. Cochrane, J. Sabater, K. J. Duncan, M. J. Hardcastle, P. Haskell, B. Mingo, H. J. A. Röttgering, D. J. B. Smith, W. L. Williams, M. Bonato, G. Calistro Rivera, F. Gao, C. L. Hale, K. Małek, G. K. Miley, I. Prandoni, L. Wang

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

Submitted 22 April, 2022; v1 submitted 15 April, 2022; originally announced April 2022.

Comments: 26 pages, 16 figures, accepted for publication in MNRAS. LF data available at: https://github.com/rohitk-10/AGN_LF_Kondapally22 . Author list error corrected in updated version

Nearby galaxies in LoTSS-DR2: insights into the non-linearity of the radio-SFR relation

Authors: V. Heesen, M. Staffehl, A. Basu, R. Beck, M. Stein, F. S. Tabatabaei, M. J. Hardcastle, K. T. Chyży, T. W. Shimwell, B. Adebahr, R. Beswick, D. J. Bomans, A. Botteon, E. Brinks, M. Brüggen, R. -J. Dettmar, A. Drabent, F. de Gasperin, G. Gürkan, G. H. Heald, C. Horellou, B. Nikiel-Wroczynski, R. Paladino, J. Piotrowska, H. J. A. Röttgering , et al. (2 additional authors not shown)

Abstract: 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… ▽ More 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 $α$ + 24-$μ$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 α\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. △ Less

Submitted 1 April, 2022; originally announced April 2022.

Comments: Accepted to Astronomy & Astrophysics. 16 pages, 8 figures, 3 tables + 45-page image atlas

Journal ref: A&A 664, A83 (2022)

Faraday tomography of LoTSS-DR2 data: I. Faraday moments in the high-latitude outer Galaxy and revealing Loop III in polarisation

Authors: Ana Erceg, Vibor Jelić, Marijke Haverkorn, Andrea Bracco, Timothy W. Shimwell, Cyril Tasse, John M. Dickey, Lana Ceraj, Alexander Drabent, Martin J. Hardcastle, Luka Turić

Abstract: Observations of synchrotron emission at low radio frequencies reveal a labyrinth of polarised Galactic structures. However, the explanation for the wealth of structures remains uncertain due to the complex interactions between the interstellar medium and the magnetic field. A multi-tracer approach to the analysis of large sky areas is needed. This paper aims to use polarimetric images from the LOF… ▽ More Observations of synchrotron emission at low radio frequencies reveal a labyrinth of polarised Galactic structures. However, the explanation for the wealth of structures remains uncertain due to the complex interactions between the interstellar medium and the magnetic field. A multi-tracer approach to the analysis of large sky areas is needed. This paper aims to use polarimetric images from the LOFAR Two metre Sky Survey (LoTSS) to produce the biggest mosaic of polarised emission in the northern sky at low radio frequencies (150 MHz) to date. The large area this mosaic covers allows for detailed morphological and statistical studies of polarised structures in the high-latitude outer Galaxy, including the well-known Loop III region. We produced a 3100 square degree Faraday tomographic cube using a rotation measure synthesis tool. We calculated the statistical moments of Faraday spectra and compared them with data sets at higher frequencies (1.4 GHz) and with a map of a rotation measure derived from extragalactic sources. The mosaic is dominated by polarised emission connected to Loop III. Additionally, the mosaic reveals an abundance of other morphological structures, mainly {narrow and extended} depolarisation canals, which are found to be ubiquitous. We find a correlation between the map of an extragalactic rotation measure and the LoTSS first Faraday moment image. The ratio of the two deviates from a simple model of a Burn slab (Burn 1966) along the line of sight, which highlights the high level of complexity in the magnetoionic medium that can be studied at these frequencies. △ Less

Submitted 4 July, 2022; v1 submitted 2 March, 2022; originally announced March 2022.

Comments: 20 pages, 25 figures, accepted for publication in A&A

Journal ref: A&A 663, A7 (2022)

The LOFAR Two-metre Sky Survey -- V. Second data release

Authors: T. W. Shimwell, M. J. Hardcastle, C. Tasse, P. N. Best, H. J. A. Röttgering, W. L. Williams, A. Botteon, A. Drabent, A. Mechev, A. Shulevski, R. J. van Weeren, L. Bester, M. Brüggen, G. Brunetti, J. R. Callingham, K. T. Chyży, J. E. Conway, T. J. Dijkema, K. Duncan, F. de Gasperin, C. L. Hale, M. Haverkorn, B. Hugo, N. Jackson, M. Mevius , et al. (81 additional authors not shown)

Abstract: 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 th… ▽ More 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$μ$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$μ$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. △ Less

Submitted 23 February, 2022; originally announced February 2022.

Comments: 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/

The Planck clusters in the LOFAR sky. I. LoTSS-DR2: new detections and sample overview

Authors: A. Botteon, T. W. Shimwell, R. Cassano, V. Cuciti, X. Zhang, L. Bruno, L. Camillini, R. Natale, A. Jones, F. Gastaldello, A. Simionescu, M. Rossetti, H. Akamatsu, R. J. van Weeren, G. Brunetti, M. Brüggen, C. Groenveld, D. N. Hoang, M. J. Hardcastle, A. Ignesti, G. Di Gennaro, A. Bonafede, A. Drabent, H. J. A. Röttgering, M. Hoeft , et al. (1 additional authors not shown)

Abstract: Relativistic electrons and magnetic fields permeate the intra-cluster medium (ICM) and manifest themselves as diffuse sources of synchrotron emission observable at radio wavelengths, namely radio halos and radio relics. Although there is broad consensus that the formation of these sources is connected to turbulence and shocks in the ICM, the details of the required particle acceleration, the stren… ▽ More Relativistic electrons and magnetic fields permeate the intra-cluster medium (ICM) and manifest themselves as diffuse sources of synchrotron emission observable at radio wavelengths, namely radio halos and radio relics. Although there is broad consensus that the formation of these sources is connected to turbulence and shocks in the ICM, the details of the required particle acceleration, the strength and morphology of the magnetic field in the cluster volume, and the influence of other sources of high-energy particles are poorly known. Sufficiently large samples of radio halos and relics, which would allow us to examine the variation among the source population and pinpoint their commonalities and differences, are still missing. At present, large numbers of these sources are easiest to detect at low radio frequencies, where they shine brightly. We examined the low-frequency radio emission from all 309 clusters in the second catalog of Planck Sunyaev Zel'dovich detected sources that lie within the 5634 deg$^2$ covered by the Second Data Release of the LOFAR Two-meter Sky Survey (LoTSS-DR2). We produced LOFAR images at different resolutions, with and without discrete sources subtracted, and created overlays with optical and X-ray images before classifying the diffuse sources in the ICM, guided by a decision tree. Overall, we found 83 clusters that host a radio halo and 26 that host one or more radio relics (including candidates). About half of them are new discoveries. The detection rate of clusters hosting a radio halo and one or more relics in our sample is $30\pm11$% and $10\pm6$%, respectively. Extrapolating these numbers, we anticipate that once LoTSS covers the entire northern sky it will provide the detection of $251\pm92$ clusters hosting a halo and $83\pm50$ clusters hosting at least one relic from Planck clusters alone. △ Less

Submitted 15 April, 2022; v1 submitted 23 February, 2022; originally announced February 2022.

Comments: Matched to the A&A published version. Abstract abridged to meet arXiv requirements. All FITS images and tables produced in this work are publicly available on the project website https://lofar-surveys.org/planck_dr2.html

Journal ref: A&A 660, A78 (2022)

The discovery of a radio galaxy of at least 5 Mpc

Authors: Martijn S. S. L. Oei, Reinout J. van Weeren, Martin J. Hardcastle, Andrea Botteon, Tim W. Shimwell, Pratik Dabhade, Aivin R. D. J. G. I. B. Gast, Huub J. A. Röttgering, Marcus Brüggen, Cyril Tasse, Wendy L. Williams, Aleksandar Shulevski

Abstract: We discover what is in projection the largest known structure of galactic origin: a giant radio galaxy with a projected proper length of $4.99 \pm 0.04\ \mathrm{Mpc}$. The source, named Alcyoneus, was first identified in low-resolution LOFAR Two-metre Sky Survey images from which angularly compact sources had been removed. Being an extreme example in its class, Alcyoneus could shed light on the ma… ▽ More We discover what is in projection the largest known structure of galactic origin: a giant radio galaxy with a projected proper length of $4.99 \pm 0.04\ \mathrm{Mpc}$. The source, named Alcyoneus, was first identified in low-resolution LOFAR Two-metre Sky Survey images from which angularly compact sources had been removed. Being an extreme example in its class, Alcyoneus could shed light on the main mechanisms that drive radio galaxy growth. We find that - beyond geometry - Alcyoneus and its host galaxy appear suspiciously ordinary: the total low-frequency luminosity density, stellar mass and supermassive black hole mass are all lower than, though similar to, those of the medial giant radio galaxy (percentiles $45 \pm 3\%$, $25 \pm 9 \%$ and $23 \pm 11 \%$, respectively). The source resides in a filament of the Cosmic Web, with which it might have significant thermodynamic interaction. At $5 \cdot 10^{-16}\ \mathrm{Pa}$, the pressures in the lobes are the lowest hitherto found, and Alcyoneus therefore represents one of the most promising radio galaxies yet to probe the warm-hot intergalactic medium. △ Less

Submitted 10 February, 2022; originally announced February 2022.

Comments: 18 pages, 14 figures, 3 tables, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 660, A2 (2022)

Radio AGN in nearby dwarf galaxies: the important role of AGN in dwarf-galaxy evolution

Authors: F. Davis, S. Kaviraj, M. J. Hardcastle, G. Martin, R. A. Jackson, K. Kraljic, K. Malek, S. Peirani, D. J. B. Smith, M. Volonteri, L. Wang

Abstract: We combine deep optical and radio data, from the Hyper Suprime-Cam and the Low-Frequency Array (LOFAR) respectively, to study 78 radio AGN in nearby (z<0.5) dwarf galaxies. Comparison to a control sample, matched in stellar mass and redshift, indicates that the AGN and controls reside in similar environments, show similar star-formation rates (which trace gas availability) and exhibit a comparable… ▽ More We combine deep optical and radio data, from the Hyper Suprime-Cam and the Low-Frequency Array (LOFAR) respectively, to study 78 radio AGN in nearby (z<0.5) dwarf galaxies. Comparison to a control sample, matched in stellar mass and redshift, indicates that the AGN and controls reside in similar environments, show similar star-formation rates (which trace gas availability) and exhibit a comparable incidence of tidal features (which indicate recent interactions). We explore the AGN properties by combining the predicted gas conditions in dwarfs from a cosmological hydrodynamical simulation with a Monte-Carlo suite of simulated radio sources, based on a semi-analytical model for radio-galaxy evolution. In the subset of LOFAR-detectable simulated sources, which have a similar distribution of radio luminosities as our observed AGN, the median jet powers, ages and accretion rates are $\sim 10^{35}$ W, $\sim 5$ Myr and $\sim 10^{-3.4}$ M$_{\odot}$ yr$^{-1}$ respectively. The median mechanical energy output of these sources is $\sim100$ times larger than the median binding energy expected in dwarf gas reservoirs, making AGN feedback plausible. Since special circumstances (in terms of environment, gas availability and interactions) are not necessary for the presence of AGN, and the central gas masses are predicted to be an order of magnitude larger than that required to fuel the AGN, AGN triggering in dwarfs is likely to be stochastic and a common phenomenon. Together with the plausibility of energetic feedback, this suggests that AGN could be important drivers of dwarf-galaxy evolution, as is the case in massive galaxies. △ Less

Submitted 24 January, 2022; originally announced January 2022.

Comments: Accepted for publication in MNRAS

The galaxy group NGC 507: newly detected AGN remnant plasma transported by sloshing

Authors: M. Brienza, L. Lovisari, K. Rajpurohit, A. Bonafede, F. Gastaldello, M. Murgia, F. Vazza, E. Bonnassieux, A. Botteon, G. Brunetti, A. Drabent, M. J. Hardcastle, T. Pasini, C. J. Riseley, H. J. A. Röttgering, T. Shimwell, A. Simionescu, R. J. van Weeren

Abstract: Jets from active galactic nuclei (AGN) are known to recurrently enrich their surrounding medium with mildly-relativistic particles and magnetic fields. Here, we present a detailed multi-frequency analysis of the nearby (z=0.01646) galaxy group NGC 507. In particular, we present new high-sensitivity and high spatial resolution radio images in the frequency range 144-675 MHz obtained using LOFAR and… ▽ More Jets from active galactic nuclei (AGN) are known to recurrently enrich their surrounding medium with mildly-relativistic particles and magnetic fields. Here, we present a detailed multi-frequency analysis of the nearby (z=0.01646) galaxy group NGC 507. In particular, we present new high-sensitivity and high spatial resolution radio images in the frequency range 144-675 MHz obtained using LOFAR and uGMRT observations. These reveal the presence of previously undetected diffuse radio emission with complex, filamentary morphology, likely related to a previous outburst of the central galaxy. Based on spectral ageing considerations, we derived that the plasma was first injected by the AGN 240-380 Myr ago and is now cooling. Our analysis of deep archival XMM-Newton data confirms that the system is dynamically disturbed, as previously suggested. We detect two discontinuities in the X-ray surface brightness distribution (in East and South direction) tracing a spiral pattern, which we interpret as cold fronts produced by sloshing motions. The remarkable spatial coincidence observed between the newly-detected arc-like radio filament and the southern concave X-ray discontinuity strongly suggests that the remnant plasma has been displaced by the sloshing motions on large scales. Overall, NGC 507 represents one of the clearest examples known to date in which a direct interaction between old AGN remnant plasma and the external medium is observed in a galaxy group. Our results are consistent with simulations, which suggest that filamentary emission can be created by the cluster/group weather disrupting AGN lobes and spreading their relativistic content into the surrounding medium. △ Less

Submitted 19 August, 2022; v1 submitted 12 January, 2022; originally announced January 2022.

Comments: 21 pages, 15 figures, 3 tables, A&A accepted

Journal ref: A&A 661, A92 (2022)

Accretion mode versus radio morphology in the LOFAR Deep Fields

Authors: B. Mingo, J. H. Croston, P. N. Best, K. J. Duncan, M. J. Hardcastle, R. Kondapally, I. Prandoni, J. Sabater, T. W. Shimwell, W. L. Williams, R. D. Baldi, M. Bonato, M. Bondi, P. Dabhade, G. Gürkan, J. Ineson, M. Magliocchetti, G. Miley, J. C. S. Pierce, H. J. A. Röttgering

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

Submitted 17 February, 2022; v1 submitted 12 January, 2022; originally announced January 2022.

Comments: 22 pages, 17 figures, 7 tables. Accepted for publication in MNRAS. Update with minor changes to match published version

A snapshot of the oldest AGN feedback phases

Authors: M. Brienza, T. W. Shimwell, F. de Gasperin, I. Bikmaev, A. Bonafede, A. Botteon, M. Brüggen, G. Brunetti, R. Burenin, A. Capetti, E. Churazov, M. J. Hardcastle, I. Khabibullin, N. Lyskova, H. J. A. Röttgering, R. Sunyaev, R. J. van Weeren, F. Gastaldello, S. Mandal, S. Purser, A. Simionescu, C. Tasse

Abstract: 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 identifi… ▽ More 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. △ Less

Submitted 18 October, 2021; originally announced October 2021.

Comments: 34 pages, 9 figures, 5 tables, published in Nature Astronomy

Low frequency radio properties of the $z>5$ quasar population

Authors: A. J. Gloudemans, K. J. Duncan, H. J. A. Röttgering, T. W. Shimwell, B. P. Venemans, P. N. Best, M. Brüggen, G. Calistro Rivera, A. Drabent, M. J. Hardcastle, G. K. Miley, D. J. Schwarz, A. Saxena, D. J. B. Smith, W. L. Williams

Abstract: Optically luminous quasars at $z > 5$ are important probes of super-massive black hole (SMBH) formation. With new and future radio facilities, the discovery of the brightest low-frequency radio sources in this epoch would be an important new probe of cosmic reionization through 21-cm absorption experiments. In this work, we systematically study the low-frequency radio properties of a sample of 115… ▽ More Optically luminous quasars at $z > 5$ are important probes of super-massive black hole (SMBH) formation. With new and future radio facilities, the discovery of the brightest low-frequency radio sources in this epoch would be an important new probe of cosmic reionization through 21-cm absorption experiments. In this work, we systematically study the low-frequency radio properties of a sample of 115 known spectroscopically confirmed $z>5$ quasars using the second data release of the Low Frequency Array (LOFAR) Two Metre Sky survey (LoTSS-DR2), reaching noise levels of $\sim$80 $μ$Jy beam$^{-1}$ (at 144 MHz) over an area of $\sim5720$ deg$^2$. We find that 41 sources (36%) are detected in LoTSS-DR2 at $>2 σ$ significance and we explore the evolution of their radio properties (power, spectral index, and radio loudness) as a function of redshift and rest-frame ultra-violet properties. We obtain a median spectral index of $-0.29^{+0.10}_{-0.09}$ by stacking 93 quasars using LoTSS-DR2 and Faint Images of the Radio Sky at Twenty Centimetres (FIRST) data at 1.4 GHz, in line with observations of quasars at $z<3$. We compare the radio loudness of the high-$z$ quasar sample to a lower-$z$ quasar sample at $z\sim2$ and find that the two radio loudness distributions are consistent with no evolution, although the low number of high-z quasars means that we cannot rule out weak evolution. Furthermore, we make a first order empirical estimate of the $z=6$ quasar radio luminosity function, which is used to derive the expected number of high-$z$ sources that will be detected in the completed LoTSS survey. This work highlights the fact that new deep radio observations can be a valuable tool in selecting high-$z$ quasar candidates for follow-up spectroscopic observations by decreasing contamination of stellar dwarfs and reducing possible selection biases introduced by strict colour cuts. △ Less

Submitted 12 October, 2021; originally announced October 2021.

Comments: 11 pages, 10 figures, accepted for publication in A&A

The Application of Ridgelines in Extended Radio Source Cross-Identification

Authors: B. Barkus, J. H. Croston, J. Piotrowska, B. Mingo, P. N. Best, M. J. Hardcastle, R. I. J. Mostert, H. J. A. Röttgering, J. Sabater, B. Webster, W. L. Williams

Abstract: Extended radio sources are an important minority population in modern deep radio surveys, because they enable detailed investigation of the physics governing radio-emitting regions such as active galaxies and their environments. Cross-identification of radio sources with optical host galaxies is challenging for this extended population, due to their morphological complexity and multiple potential… ▽ More Extended radio sources are an important minority population in modern deep radio surveys, because they enable detailed investigation of the physics governing radio-emitting regions such as active galaxies and their environments. Cross-identification of radio sources with optical host galaxies is challenging for this extended population, due to their morphological complexity and multiple potential counterparts. In the first data release of the Low-frequency array (LOFAR) Two-metre Sky Survey (LoTSS DR1) the automated likelihood ratio for compact sources was supplemented by a citizen science visual identification process for extended sources. In this paper we present a novel method for automating the host identification of extended sources by using ridgelines, which trace the assumed direction of fluid-flow through the points of highest flux density. Applying a new code, RL-Xid, to LoTSS DR1, we demonstrate that ridgelines are versatile; by providing information about spatial structure and brightness distributions, they can be used both for optical host identification and morphological studies in radio surveys. RL-Xid draws ridgelines for 85 per cent of sources brighter than 10 mJy and larger than 15 arcsec, with an improved performance of 96 per cent for the subset >30 mJy and >60 arcsec. Using a sample of sources with known hosts from LoTSS DR1, we demonstrate that RL-Xid successfully identifies the host for 98 per cent of the sources with successfully drawn ridgelines, and performs at a comparable level to visual identification via citizen science. We also demonstrate that ridgeline brightness profiles provide a promising automated technique for morphological classification. △ Less

Submitted 12 October, 2021; v1 submitted 11 October, 2021; originally announced October 2021.

Comments: 16 pages, 7 figures, accepted for publication in MNRAS

Investigating the spectra and physical nature of galaxy scale jets

Authors: B. Webster, J. H. Croston, J. J. Harwood, R. D. Baldi, M. J. Hardcastle, B. Mingo, H. J. A. Rottgering

Abstract: Previous studies have shown that physically small, low-luminosity radio galaxies, which we refer to as galaxy scale jets (GSJ), could potentially have a significant effect upon the host galaxy's evolution. Using 6 arcsec resolution images taken from the first release of the LOFAR Two Metre Sky Survey (LoTSS DR1), we identified a representative sample of nine potential GSJ for which we obtained hig… ▽ More Previous studies have shown that physically small, low-luminosity radio galaxies, which we refer to as galaxy scale jets (GSJ), could potentially have a significant effect upon the host galaxy's evolution. Using 6 arcsec resolution images taken from the first release of the LOFAR Two Metre Sky Survey (LoTSS DR1), we identified a representative sample of nine potential GSJ for which we obtained high-resolution, 2-4 GHz data using the Karl G. Jansky Very Large Array (VLA). Using these data we aim to verify the GSJ nature of these sources as well as investigating the potential role of feedback. Our VLA images reveal a diversity of structures, confirm the hosts for four of the sources and find that a fifth is the first known example of a galaxy-scale remnant showing that some radio galaxies never grow beyond the GSJ stage. We also derive spectral ages and the first estimates of the lobe expansion speeds of GSJ. We find our GSJ have maximum spectral ages of 60 Myr with most between about 5 and 20 Myr, consistent with being located along an evolutionary path joining compact sources and larger radio galaxies. We find lobe advance speeds a few times the local sound speed, with most GSJ predicted to be driving strong shocks into their environment and having a significant impact upon the host's evolution. Our discovery of a remnant GSJ, which will eventually transfer all of its energy directly into the local environment, represents an important and previously hidden aspect of AGN life cycles. △ Less

Submitted 8 October, 2021; originally announced October 2021.

Comments: 19 pages, 12 figures, accepted by MNRAS