-
fast-resolve: Fast Bayesian Radio Interferometric Imaging
Authors:
Jakob Roth,
Philipp Frank,
Hertzog L. Bester,
Oleg M. Smirnov,
Rüdiger Westermann,
Torsten A. Enßlin
Abstract:
Context: Interferometric imaging is algorithmically and computationally challenging as there is no unique inversion from the measurement data back to the sky maps, and the datasets can be very large. Many imaging methods already exist, but most of them focus either on the accuracy or the computational aspect. Aims: This paper aims to reduce the computational complexity of the Bayesian imaging algo…
▽ More
Context: Interferometric imaging is algorithmically and computationally challenging as there is no unique inversion from the measurement data back to the sky maps, and the datasets can be very large. Many imaging methods already exist, but most of them focus either on the accuracy or the computational aspect. Aims: This paper aims to reduce the computational complexity of the Bayesian imaging algorithm resolve, enabling the application of Bayesian imaging for larger datasets. Methods: By combining computational shortcuts of the CLEAN algorithm with the Bayesian imaging algorithm resolve we developed an accurate and fast imaging algorithm which we name fast-resolve. Results: We validate the accuracy of the presented fast-resolve algorithm by comparing it with results from resolve on VLA Cygnus A data. Furthermore, we demonstrate the computational advantages of fast-resolve on a large MeerKAT ESO 137-006 dataset which is computationally out of reach for resolve. Conclusions: The presented algorithm is significantly faster than previous Bayesian imaging algorithms, broadening the applicability of Bayesian interferometric imaging. Specifically for the single channel VLA Cygnus A datasets fast-resolve is about $144$ times faster than resolve. For the MeerKAT dataset with multiple channels the computational speedup of fast-resolve is even larger.
△ Less
Submitted 13 June, 2024;
originally announced June 2024.
-
The RATT PARROT: serendipitous discovery of a peculiarly scintillating pulsar in MeerKAT imaging observations of the Great Saturn-Jupiter Conjunction of 2020. I. Dynamic imaging and data analysis
Authors:
O. M. Smirnov,
B. W. Stappers,
C. Tasse,
H. L. Bester,
H. Bignall,
M. A. Walker,
M. Caleb,
K. M. Rajwade,
S. Buchner,
P. Woudt,
M. Ivchenko,
L. Roth,
J. E. Noordam,
F. Camilo
Abstract:
We report on a radiopolarimetric observation of the Saturn-Jupiter Great Conjunction of 2020 using the MeerKAT L-band system, initially carried out for science verification purposes, which yielded a serendipitous discovery of a pulsar. The radiation belts of Jupiter are very bright and time variable: coupled with the sensitivity of MeerKAT, this necessitated development of dynamic imaging techniqu…
▽ More
We report on a radiopolarimetric observation of the Saturn-Jupiter Great Conjunction of 2020 using the MeerKAT L-band system, initially carried out for science verification purposes, which yielded a serendipitous discovery of a pulsar. The radiation belts of Jupiter are very bright and time variable: coupled with the sensitivity of MeerKAT, this necessitated development of dynamic imaging techniques, reported on in this work. We present a deep radio "movie" revealing Jupiter's rotating magnetosphere, a radio detection of Callisto, and numerous background radio galaxies. We also detect a bright radio transient in close vicinity to Saturn, lasting approximately 45 minutes. Follow-up deep imaging observations confirmed this as a faint compact variable radio source, and yielded detections of pulsed emission by the commensal MeerTRAP search engine, establishing the object's nature as a radio emitting neutron star, designated PSR J2009-2026. A further observation combining deep imaging with the PTUSE pulsar backend measured detailed dynamic spectra for the object. While qualitatively consistent with scintillation, the magnitude of the magnification events and the characteristic timescales are odd. We are tentatively designating this object a pulsar with anomalous refraction recurring on odd timescales (PARROT). As part of this investigation, we present a pipeline for detection of variable sources in imaging data, with dynamic spectra and lightcurves as the products, and compare dynamic spectra obtained from visibility data with those yielded by PTUSE. We discuss MeerKAT's capabilities and prospects for detecting more of such transients and variables.
△ Less
Submitted 20 December, 2023; v1 submitted 19 December, 2023;
originally announced December 2023.
-
MeerKAT follow-up of enigmatic GLEAM 4-Jy (G4Jy) sources
Authors:
Precious K. Sejake,
Sarah V. White,
Ian Heywood,
Kshitij Thorat,
Hertzog L. Bester,
Sphesihle Makhathini,
Bernie Fanaroff
Abstract:
We present the results from studying 140 radio sources in the GLEAM (GaLactic and Extragalactic All-sky MWA [Murchison Widefield Array]) 4-Jy (G4Jy) Sample. These sources were followed-up with MeerKAT to assess their radio morphology and enable host-galaxy identification, as existing radio images of 25 to 45-arcsec resolution do not provide sufficient information. We refer to these sources as the…
▽ More
We present the results from studying 140 radio sources in the GLEAM (GaLactic and Extragalactic All-sky MWA [Murchison Widefield Array]) 4-Jy (G4Jy) Sample. These sources were followed-up with MeerKAT to assess their radio morphology and enable host-galaxy identification, as existing radio images of 25 to 45-arcsec resolution do not provide sufficient information. We refer to these sources as the MeerKAT-2019 subset. The aim is to identify the host galaxy of these sources by visually inspecting the overlays comprising radio data from four surveys (at 150, 200, 843/1400, and 1300 MHz). Our morphological classification and host-galaxy identification relies upon the ~7-arcsec resolution images from MeerKAT (1300 MHz). Through the visual inspection of the overlays, 14 radio sources in the MeerKAT-2019 subset have wide-angle tail (WAT) morphology, 10 are head-tail, and 5 have X-, S-/Z-shaped morphology. Most of the remaining sources have the radio morphology of typical symmetric lobes. Of 140 sources, we find host galaxies for 98 sources, leaving 42 with no identified host galaxy. These 42 sources still have ambiguous identification even with higher resolution images from MeerKAT.
△ Less
Submitted 13 February, 2023; v1 submitted 22 November, 2022;
originally announced November 2022.
-
The 1.28 GHz MeerKAT Galactic Center Mosaic
Authors:
I. Heywood,
I. Rammala,
F. Camilo,
W. D. Cotton,
F. Yusef-Zadeh,
T. D. Abbott,
R. M. Adam,
G. Adams,
M. A. Aldera,
K. M. B. Asad,
E. F. Bauermeister,
T. G. H. Bennett,
H. L. Bester,
W. A. Bode,
D. H. Botha,
A. G. Botha,
L. R. S. Brederode,
S. Buchner,
J. P. Burger,
T. Cheetham,
D. I. L. de Villiers,
M. A. Dikgale-Mahlakoana,
L. J. du Toit,
S. W. P. Esterhuyse,
B. L. Fanaroff
, et al. (86 additional authors not shown)
Abstract:
The inner $\sim$200 pc region of the Galaxy contains a 4 million M$_{\odot}$ supermassive black hole (SMBH), significant quantities of molecular gas, and star formation and cosmic ray energy densities that are roughly two orders of magnitude higher than the corresponding levels in the Galactic disk. At a distance of only 8.2 kpc, the region presents astronomers with a unique opportunity to study a…
▽ More
The inner $\sim$200 pc region of the Galaxy contains a 4 million M$_{\odot}$ supermassive black hole (SMBH), significant quantities of molecular gas, and star formation and cosmic ray energy densities that are roughly two orders of magnitude higher than the corresponding levels in the Galactic disk. At a distance of only 8.2 kpc, the region presents astronomers with a unique opportunity to study a diverse range of energetic astrophysical phenomena, from stellar objects in extreme environments, to the SMBH and star-formation driven feedback processes that are known to influence the evolution of galaxies as a whole. We present a new survey of the Galactic center conducted with the South African MeerKAT radio telescope. Radio imaging offers a view that is unaffected by the large quantities of dust that obscure the region at other wavelengths, and a scene of striking complexity is revealed. We produce total intensity and spectral index mosaics of the region from 20 pointings (144 hours on-target in total), covering 6.5 square degrees with an angular resolution of 4$"$,at a central frequency of 1.28 GHz. Many new features are revealed for the first time due to a combination of MeerKAT's high sensitivity, exceptional $u,v$-plane coverage, and geographical vantage point. We highlight some initial survey results, including new supernova remnant candidates, many new non-thermal filament complexes, and enhanced views of the Radio Arc Bubble, Sgr A and Sgr B regions. This project is a SARAO public legacy survey, and the image products are made available with this article.
△ Less
Submitted 27 January, 2022; v1 submitted 25 January, 2022;
originally announced January 2022.
-
Simulations of primary beam effects on the cosmic bispectrum phase observed with the Hydrogen Epoch of Reionization Array
Authors:
N. Charles,
G. Bernardi,
H. L. Bester,
O. M. Smirnov,
C. Carilli,
P. M. Keller,
N. Kern,
B. Nikolic,
N. Thyagarajan,
E. de Lera Acedo,
N. Fagnoni,
M. G. Santos
Abstract:
The 21~cm transition from neutral Hydrogen promises to be the best observational probe of the Epoch of Reionisation. The main difficulty in measuring the 21 cm signal is the presence of bright foregrounds that require very accurate interferometric calibration. Closure quantities may circumvent the calibration requirements but may be, however, affected by direction dependent effects, particularly a…
▽ More
The 21~cm transition from neutral Hydrogen promises to be the best observational probe of the Epoch of Reionisation. The main difficulty in measuring the 21 cm signal is the presence of bright foregrounds that require very accurate interferometric calibration. Closure quantities may circumvent the calibration requirements but may be, however, affected by direction dependent effects, particularly antenna primary beam responses. This work investigates the impact of antenna primary beams affected by mutual coupling on the closure phase and its power spectrum. Our simulations show that primary beams affected by mutual coupling lead to a leakage of foreground power into the EoR window, which can be up to $\sim4$ orders magnitude higher than the case where no mutual coupling is considered. This leakage is, however, essentially confined at $k < 0.3$~$h$~Mpc$^{-1}$ for triads that include 29~m baselines. The leakage magnitude is more pronounced when bright foregrounds appear in the antenna sidelobes, as expected. Finally, we find that triads that include mutual coupling beams different from each other have power spectra similar to triads that include the same type of mutual coupling beam, indicating that beam-to-beam variation within triads (or visibility pairs) is not the major source of foreground leakage in the EoR window.
△ Less
Submitted 19 May, 2022; v1 submitted 23 December, 2021;
originally announced December 2021.
-
MIGHTEE: total intensity radio continuum imaging and the COSMOS / XMM-LSS Early Science fields
Authors:
I. Heywood,
M. J. Jarvis,
C. L. Hale,
I. H. Whittam,
H. L. Bester,
B. Hugo,
J. S. Kenyon,
M. Prescott,
O. M. Smirnov,
C. Tasse,
J. M. Afonso,
P. N. Best,
J. D. Collier,
R. P. Deane,
B. S. Frank,
M. J. Hardcastle,
K. Knowles,
N. Maddox,
E. J. Murphy,
I. Prandoni,
S. M. Randriamampandry,
M. G. Santos,
S. Sekhar,
F. Tabatabaei,
A. R. Taylor
, et al. (1 additional authors not shown)
Abstract:
MIGHTEE is a galaxy evolution survey using simultaneous radio continuum, spectro-polarimetry, and spectral line observations from the South African MeerKAT telescope. When complete, the survey will image $\sim$20 deg$^{2}$ over the COSMOS, E-CDFS, ELAIS-S1, and XMM-LSS extragalactic deep fields with a central frequency of 1284 MHz. These were selected based on the extensive multiwavelength dataset…
▽ More
MIGHTEE is a galaxy evolution survey using simultaneous radio continuum, spectro-polarimetry, and spectral line observations from the South African MeerKAT telescope. When complete, the survey will image $\sim$20 deg$^{2}$ over the COSMOS, E-CDFS, ELAIS-S1, and XMM-LSS extragalactic deep fields with a central frequency of 1284 MHz. These were selected based on the extensive multiwavelength datasets from numerous existing and forthcoming observational campaigns. Here we describe and validate the data processing strategy for the total intensity continuum aspect of MIGHTEE, using a single deep pointing in COSMOS (1.6 deg$^{2}$) and a three-pointing mosaic in XMM-LSS (3.5 deg$^{2}$). The processing includes the correction of direction-dependent effects, and results in thermal noise levels below 2~$\mathrmμ$Jy beam$^{-1}$ in both fields, limited in the central regions by classical confusion at $\sim$8$''$ angular resolution, and meeting the survey specifications. We also produce images at $\sim$5$''$ resolution that are $\sim$3 times shallower. The resulting image products form the basis of the Early Science continuum data release for MIGHTEE. From these images we extract catalogues containing 9,896 and 20,274 radio components in COSMOS and XMM-LSS respectively. We also process a close-packed mosaic of 14 additional pointings in COSMOS and use these in conjunction with the Early Science pointing to investigate methods for primary beam correction of broadband radio images, an analysis that is of relevance to all full-band MeerKAT continuum observations, and wide field interferometric imaging in general. A public release of the MIGHTEE Early Science continuum data products accompanies this article.
△ Less
Submitted 1 October, 2021;
originally announced October 2021.
-
Solution intervals considered harmful: on the optimality of radio interferometric gain solutions
Authors:
Ulrich Armel Mbou Sob,
Hertzog Landman Bester,
Oleg M. Smirnov,
Jonathan Kenyon,
Cyndie Russeeawon
Abstract:
Solution intervals are often used to improve the signal-to-noise ratio during radio interferometric gain calibration. This work investigates how factors such as the noise level, intrinsic gain variability, degree of model incompleteness, and the presence of radio frequency interference impact the selection of solution intervals for calibration. We perform different interferometric simulations to d…
▽ More
Solution intervals are often used to improve the signal-to-noise ratio during radio interferometric gain calibration. This work investigates how factors such as the noise level, intrinsic gain variability, degree of model incompleteness, and the presence of radio frequency interference impact the selection of solution intervals for calibration. We perform different interferometric simulations to demonstrate how these factors, in combination with the choice of solution intervals, affect calibration and imaging outputs and discuss practical guidelines for choosing optimal solution intervals. Furthermore, we present an algorithm capable of automatically selecting suitable solution intervals during calibration. By applying the algorithm to both simulated and real data, we show that it can successfully choose solution intervals that strike a good balance between capturing intrinsic gain variability and not fitting noise as long as the data are not too inhomogeneously flagged. Furthermore, we elaborate on several practical aspects that emphasize the need to develop regularised calibration algorithms that do not require solution intervals.
△ Less
Submitted 23 March, 2021;
originally announced March 2021.
-
A practical preconditioner for wide-field continuum imaging of radio interferometric data
Authors:
Hertzog L. Bester,
Audrey Repetti,
Simon Perkins,
Oleg M. Smirnov,
Jonathan S. Kenyon
Abstract:
The celebrated CLEAN algorithm has been the cornerstone of deconvolution algorithms in radio interferometry almost since its conception in the 1970s. For all its faults, CLEAN is remarkably fast, robust to calibration artefacts and in its ability to model point sources. We demonstrate how the same assumptions that afford CLEAN its speed can be used to accelerate more sophisticated deconvolution al…
▽ More
The celebrated CLEAN algorithm has been the cornerstone of deconvolution algorithms in radio interferometry almost since its conception in the 1970s. For all its faults, CLEAN is remarkably fast, robust to calibration artefacts and in its ability to model point sources. We demonstrate how the same assumptions that afford CLEAN its speed can be used to accelerate more sophisticated deconvolution algorithms.
△ Less
Submitted 20 January, 2021;
originally announced January 2021.
-
Comparison of classical and Bayesian imaging in radio interferometry
Authors:
Philipp Arras,
Hertzog L. Bester,
Richard A. Perley,
Reimar Leike,
Oleg Smirnov,
Rüdiger Westermann,
Torsten A. Enßlin
Abstract:
CLEAN, the commonly employed imaging algorithm in radio interferometry, suffers from a number of shortcomings: in its basic version it does not have the concept of diffuse flux, and the common practice of convolving the CLEAN components with the CLEAN beam erases the potential for super-resolution; it does not output uncertainty information; it produces images with unphysical negative flux regions…
▽ More
CLEAN, the commonly employed imaging algorithm in radio interferometry, suffers from a number of shortcomings: in its basic version it does not have the concept of diffuse flux, and the common practice of convolving the CLEAN components with the CLEAN beam erases the potential for super-resolution; it does not output uncertainty information; it produces images with unphysical negative flux regions; and its results are highly dependent on the so-called weighting scheme as well as on any human choice of CLEAN masks to guiding the imaging. Here, we present the Bayesian imaging algorithm resolve which solves the above problems and naturally leads to super-resolution. We take a VLA observation of Cygnus~A at four different frequencies and image it with single-scale CLEAN, multi-scale CLEAN and resolve. Alongside the sky brightness distribution resolve estimates a baseline-dependent correction function for the noise budget, the Bayesian equivalent of weighting schemes. We report noise correction factors between 0.4 and 429. The enhancements achieved by resolve come at the cost of higher computational effort.
△ Less
Submitted 25 January, 2021; v1 submitted 26 August, 2020;
originally announced August 2020.
-
Radio Interferometric Calibration Using a Complex Student's t-distribution and Wirtinger Derivatives
Authors:
Ulrich Armel Mbou Sob,
Hertzog Landman Bester,
Oleg Smirnov,
Jonathan Kenyon,
Trienko Grobler
Abstract:
Radio interferometric gain calibration can be biased by incomplete sky models and radio frequency interference, resulting in calibration artefacts that can restrict the dynamic range of the resulting images. It has been suggested that calibration algorithms employing heavy-tailed likelihood functions are less susceptible to this due to their robustness against outliers in the data. We present an a…
▽ More
Radio interferometric gain calibration can be biased by incomplete sky models and radio frequency interference, resulting in calibration artefacts that can restrict the dynamic range of the resulting images. It has been suggested that calibration algorithms employing heavy-tailed likelihood functions are less susceptible to this due to their robustness against outliers in the data. We present an algorithm based on a Student's t-distribution which leverages the framework of complex optimisation and Wirtinger calculus for efficient and robust interferometric gain calibration. We integrate this algorithm as an option in the newly released calibration software package, CubiCal. We demonstrate that the algorithm can mitigate some of the biases introduced by incomplete sky models and radio frequency interference by applying it to both simulated and real data. Our results show significant improvements compared to a conventional least-squares solver which assumes a Gaussian likelihood function. Furthermore, we provide some insight into why the algorithm outperforms the conventional solver, and discuss specific scenarios (for both direction-independent and direction-dependent self-calibration) where this is expected to be the case.
△ Less
Submitted 17 October, 2019;
originally announced October 2019.
-
Testing the Copernican principle with future radio-astronomy observations
Authors:
Hertzog L. Bester,
Julien Larena,
Nigel T. Bishop
Abstract:
We use a direct observational approach to investigate the possibility of testing the Copernican principle with data from upcoming radio surveys. In particular we illustrate the importance of measuring derivatives transverse to the past light-cone when prior knowledge of the value of the cosmological constant is not available.
We use a direct observational approach to investigate the possibility of testing the Copernican principle with data from upcoming radio surveys. In particular we illustrate the importance of measuring derivatives transverse to the past light-cone when prior knowledge of the value of the cosmological constant is not available.
△ Less
Submitted 2 May, 2017;
originally announced May 2017.
-
Numerically reconstructing the geometry of the Universe from data
Authors:
Hertzog L. Bester,
Julien Larena,
Nigel T. Bishop
Abstract:
We give an outline of an algorithm designed to reconstruct the background cosmological metric within the class of spherically symmetric dust universes that may include a cosmological constant. Luminosity and age data are used to derive constraints on the geometry of the universe up to a redshift of $z = 1.75$. It is shown that simple radially inhomogeneous void models that are sometimes used as al…
▽ More
We give an outline of an algorithm designed to reconstruct the background cosmological metric within the class of spherically symmetric dust universes that may include a cosmological constant. Luminosity and age data are used to derive constraints on the geometry of the universe up to a redshift of $z = 1.75$. It is shown that simple radially inhomogeneous void models that are sometimes used as alternative explanations for the apparent acceleration of the late time Universe cannot be ruled out by these data alone.
△ Less
Submitted 28 January, 2016; v1 submitted 27 January, 2016;
originally announced January 2016.
-
Towards the geometry of the universe from data
Authors:
Hertzog L. Bester,
Julien Larena,
Nigel T. Bishop
Abstract:
We present a new algorithm that can reconstruct the full distributions of metric components within the class of spherically symmetric dust universes that may include a cosmological constant. The algorithm is capable of confronting this class of solutions with arbitrary data and opens a new observational window to determine the value of the cosmological constant. In this work we use luminosity and…
▽ More
We present a new algorithm that can reconstruct the full distributions of metric components within the class of spherically symmetric dust universes that may include a cosmological constant. The algorithm is capable of confronting this class of solutions with arbitrary data and opens a new observational window to determine the value of the cosmological constant. In this work we use luminosity and age data to constrain the geometry of the universe up to a redshift of $z = 1.75$. We show that, although current data are perfectly compatible with homogeneous models of the universe, simple radially inhomogeneous void models that are sometimes used as alternative explanations for the apparent acceleration of the late time universe cannot yet be ruled out. In doing so we reconstruct the density of cold dark matter out to $z = 1.75$ and derive constraints on the metric components when the universe was 10.5 Gyr old within a comoving volume of approximately 1 Gpc$^{3}$.
△ Less
Submitted 4 September, 2015; v1 submitted 4 June, 2015;
originally announced June 2015.
-
What's Inside the Cone? Numerically reconstructing the metric from observations
Authors:
Hertzog L. Bester,
Julien Larena,
Petrus J. van der Walt,
Nigel T. Bishop
Abstract:
We investigate the possibility of using Gaussian process regression to smooth data on the current past null-cone for use as the input to a relativistic integration scheme. The algorithm we present is designed to reconstruct the metric of spacetime within the class of spherically symmetric dust universes, with or without a cosmological constant. Assuming that gravity is well described by General Re…
▽ More
We investigate the possibility of using Gaussian process regression to smooth data on the current past null-cone for use as the input to a relativistic integration scheme. The algorithm we present is designed to reconstruct the metric of spacetime within the class of spherically symmetric dust universes, with or without a cosmological constant. Assuming that gravity is well described by General Relativity, we demonstrate how the algorithm can be employed to test the Copernican principle based on currently available observations. It is shown that currently available data is not sufficient for a conclusive result. The intrinsic noise present in realistic data presents a challenge for our smoothing algorithm and we discuss some of its limitations as well as possible extensions to it. We conclude by demonstrating how a direct determination of the cosmological constant is possible using redshift drift data.
△ Less
Submitted 11 February, 2014; v1 submitted 4 December, 2013;
originally announced December 2013.
