-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
Filamentary Baryons and Where to Find Them: A forecast of synchrotron radiation from merger and accretion shocks in the local Cosmic Web
Authors:
Martijn S. S. L. Oei,
Reinout J. van Weeren,
Franco Vazza,
Florent Leclercq,
Akshatha Gopinath,
Huub J. A. Röttgering
Abstract:
We generate probabilistic predictions of the low-redshift ($z < 0.2$) synchrotron Cosmic Web for half of the Northern Sky. In particular, we predict the contribution to the specific intensity function at $ν_\mathrm{obs} = 150\ \mathrm{MHz}$ from merger shocks in clusters and accretion shocks in filaments, both of which arise during large-scale structure formation. We assume a primordial magnetogen…
▽ More
We generate probabilistic predictions of the low-redshift ($z < 0.2$) synchrotron Cosmic Web for half of the Northern Sky. In particular, we predict the contribution to the specific intensity function at $ν_\mathrm{obs} = 150\ \mathrm{MHz}$ from merger shocks in clusters and accretion shocks in filaments, both of which arise during large-scale structure formation. We assume a primordial magnetogenesis scenario, but our method is general enough to allow for an exploration of alternative magnetogenesis scenarios - and even alternative radiation mechanisms and spectral windows. In the future, by comparing different predictions, one could infer the most plausible physical model from a synchrotron Cosmic Web detection. Our method combines Bayesian large-scale structure reconstructions, snapshots of an MHD cosmological simulation, a Gaussian random field, and a ray tracing approach. The results help to select targets for deep observations and can be used in actual detection experiments. We highlight predictions for the Hercules Cluster, the Coma Cluster, Abell 2199, Abell 2255, the Lockman Hole and the Ursa Major Supercluster. At degree-scale resolution, the median specific intensity reaches $m_{I_ν} \sim 10^{-1}\ ξ_e\ \mathrm{Jy\ deg^{-2}}$, where $ξ_e$ is the electron acceleration efficiency.
△ Less
Submitted 10 March, 2022;
originally announced March 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.
-
A probabilistic approach to direction-dependent ionospheric calibration
Authors:
J. G. Albert,
M. S. S. L. Oei,
R. J. van Weeren,
H. T. Intema,
H. J. A. Röttgering
Abstract:
Calibrating for direction-dependent ionospheric distortions in visibility data is one of the main technical challenges that must be overcome to advance low-frequency radio astronomy. In this paper, we propose a novel probabilistic, tomographic approach that utilises Gaussian processes to calibrate direction-dependent ionospheric phase distortions in low-frequency interferometric data. We suggest t…
▽ More
Calibrating for direction-dependent ionospheric distortions in visibility data is one of the main technical challenges that must be overcome to advance low-frequency radio astronomy. In this paper, we propose a novel probabilistic, tomographic approach that utilises Gaussian processes to calibrate direction-dependent ionospheric phase distortions in low-frequency interferometric data. We suggest that the ionospheric free electron density can be modelled to good approximation by a Gaussian process restricted to a thick single layer, and show that under this assumption the differential total electron content must also be a Gaussian process. We perform a comparison with a number of other widely successful Gaussian processes on simulated differential total electron contents over a wide range of experimental conditions, and find that, in all experimental conditions, our model is better able to represent observed data and generalise to unseen data. The mean equivalent source shift imposed by our predictive errors are half as large as the best competitor model's. We find that it is possible to partially constrain the ionosphere's hyperparameters from sparse-and-noisy observed data. Our model provides an alternative explanation for observed phase structure functions deviating from Kolmogorov's 5/3 turbulence, turnover at high baselines, and diffractive scale anisotropy. We show that our model implicitly cheaply performs tomography of the free electron density. Moreover, we find that even a fast, low-resolution approximation to our model yields better results than the best alternative Gaussian process, implying that the geometric coupling between directions and antennae is a powerful prior that should not be ignored.
△ Less
Submitted 24 October, 2019; v1 submitted 10 October, 2019;
originally announced October 2019.