-
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.
-
A Spectropolarimetric Study of Pictor A Radio Galaxy with MeerKAT
Authors:
Lexy A. L. Andati,
Lerato M. Baidoo,
Athanaseus J. T. Ramaila,
Oleg M. Smirnov,
Sphesihle Makhathini,
Richard A. Perley
Abstract:
We present the results of a polarimetric study from our new high-sensitivity L-band (0.8--1.7 GHz) observation of Pictor A with the MeerKAT radio telescope. We confirm the presence of the radio jet extending from the nucleus to the western hotspot of this source. Additionally, we show the radio emission expected to be coincident with previously observed X-ray emission in the radio lobes, confirmin…
▽ More
We present the results of a polarimetric study from our new high-sensitivity L-band (0.8--1.7 GHz) observation of Pictor A with the MeerKAT radio telescope. We confirm the presence of the radio jet extending from the nucleus to the western hotspot of this source. Additionally, we show the radio emission expected to be coincident with previously observed X-ray emission in the radio lobes, confirming that the emission mechanism is of inverse Compton origin, as suggested by a previous study. Our spectropolarimetric analysis using the RM-Synthesis technique reveals a relatively uniform mean RM distribution across the lobes of Pictor A, with most lines-of-sight exhibiting single-peaked Faraday spectra. However, a number of the lines-of-sight exhibit single peaked spectra with a wide base or multiple peaks, suggesting the presence of multiple Faraday components or a Faraday thick structure along Pictor A's lines-of-sight. We also confirm the asymmetry in RM variability and depolarization between the two lobes of this source which were reported in a previous study.
△ Less
Submitted 27 February, 2024;
originally announced February 2024.
-
The VLBA CANDELS GOODS-North Survey. I -- Survey Design, Processing, Data Products, and Source Counts
Authors:
Roger P. Deane,
Jack F. Radcliffe,
Ann Njeri,
Alexander Akoto-Danso,
Gianni Bernardi,
Oleg M. Smirnov,
Rob Beswick,
Michael A. Garrett,
Matt J. Jarvis,
Imogen H. Whittam,
Stephen Bourke,
Zsolt Paragi
Abstract:
The past decade has seen significant advances in wide-field cm-wave very long baseline interferometry (VLBI), which is timely given the wide-area, synoptic survey-driven strategy of major facilities across the electromagnetic spectrum. While wide-field VLBI poses significant post-processing challenges that can severely curtail its potential scientific yield, many developments in the km-scale conne…
▽ More
The past decade has seen significant advances in wide-field cm-wave very long baseline interferometry (VLBI), which is timely given the wide-area, synoptic survey-driven strategy of major facilities across the electromagnetic spectrum. While wide-field VLBI poses significant post-processing challenges that can severely curtail its potential scientific yield, many developments in the km-scale connected-element interferometer sphere are directly applicable to addressing these. Here we present the design, processing, data products, and source counts from a deep (11 $μ$Jy beam$^{-1}$), quasi-uniform sensitivity, contiguous wide-field (160 arcmin$^2$) 1.6 GHz VLBI survey of the CANDELS GOODS-North field. This is one of the best-studied extragalactic fields at milli-arcsecond resolution and, therefore, is well-suited as a comparative study for our Tera-pixel VLBI image. The derived VLBI source counts show consistency with those measured in the COSMOS field, which broadly traces the AGN population detected in arcsecond-scale radio surveys. However, there is a distinctive flattening in the $ S_{\rm 1.4GHz}\sim$100-500 $μ$Jy flux density range, which suggests a transition in the population of compact faint radio sources, qualitatively consistent with the excess source counts at 15 GHz that is argued to be an unmodelled population of radio cores. This survey approach will assist in deriving robust VLBI source counts and broadening the discovery space for future wide-field VLBI surveys, including VLBI with the Square Kilometre Array, which will include new large field-of-view antennas on the African continent at $\gtrsim$1000~km baselines. In addition, it may be useful in the design of both monitoring and/or rapidly triggered VLBI transient programmes.
△ Less
Submitted 22 January, 2024;
originally announced January 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.
-
PolarVis: Towards Web-based Polarimetric Analysis
Authors:
Lexy A. L. Andati,
O. M. Smirnov,
S. Makhathini,
L. M. Sebokolodi
Abstract:
Astronomers performing polarimetric analysis on astronomical images often have to manually identify locations on their objects of interest, such as galaxies, which exhibit the influence of magnetic forces due to interaction with their environments or inherent processes. These locations are known as Lines of Sight (LoS). Analysing the various lines of sight can provide insight into the electromagne…
▽ More
Astronomers performing polarimetric analysis on astronomical images often have to manually identify locations on their objects of interest, such as galaxies, which exhibit the influence of magnetic forces due to interaction with their environments or inherent processes. These locations are known as Lines of Sight (LoS). Analysing the various lines of sight can provide insight into the electromagnetic nature of the astrophysical object in question and its surroundings. For each LoS, astronomers generate diagnostic plots to map out the variation of the corresponding electromagnetic field, such as those of fractional polarisation and Faraday spectra. However, associating the different LoS diagnostic plots to their positions on an astronomical image requires alternating between the plots and the images. As a result, determining whether the location of the LoS influences its magnetic field variation by analysing its diagnostic plots becomes arduous due to the absence of a direct way of linking the two. PolarVis is an effort towards allowing an almost instant view of the interactive diagnostic plots corresponding to a given line of sight at the click of a button on that line of sight on the image, using an interactive web-based FITS viewer -- JS9.
△ Less
Submitted 12 December, 2023;
originally announced December 2023.
-
The SARAO MeerKAT 1.3 GHz Galactic Plane Survey
Authors:
S. Goedhart,
W. D. Cotton,
F. Camilo,
M. A. Thompson,
G. Umana,
M. Bietenholz,
P. A. Woudt,
L. D. Anderson,
C. Bordiu,
D. A. H. Buckley,
C. S. Buemi,
F. Bufano,
F. Cavallaro,
H. Chen,
J. O. Chibueze,
D. Egbo,
B. S. Frank,
M. G. Hoare,
A. Ingallinera,
T. Irabor,
R. C. Kraan-Korteweg,
S. Kurapati,
P. Leto,
S. Loru,
M. Mutale
, et al. (105 additional authors not shown)
Abstract:
We present the SARAO MeerKAT Galactic Plane Survey (SMGPS), a 1.3 GHz continuum survey of almost half of the Galactic Plane (251°$\le l \le$ 358°and 2°$\le l \le$ 61°at $|b| \le 1.5°$). SMGPS is the largest, most sensitive and highest angular resolution 1 GHz survey of the Plane yet carried out, with an angular resolution of 8" and a broadband RMS sensitivity of $\sim$10--20 $μ$ Jy/beam. Here we d…
▽ More
We present the SARAO MeerKAT Galactic Plane Survey (SMGPS), a 1.3 GHz continuum survey of almost half of the Galactic Plane (251°$\le l \le$ 358°and 2°$\le l \le$ 61°at $|b| \le 1.5°$). SMGPS is the largest, most sensitive and highest angular resolution 1 GHz survey of the Plane yet carried out, with an angular resolution of 8" and a broadband RMS sensitivity of $\sim$10--20 $μ$ Jy/beam. Here we describe the first publicly available data release from SMGPS which comprises data cubes of frequency-resolved images over 908--1656 MHz, power law fits to the images, and broadband zeroth moment integrated intensity images. A thorough assessment of the data quality and guidance for future usage of the data products are given. Finally, we discuss the tremendous potential of SMGPS by showcasing highlights of the Galactic and extragalactic science that it permits. These highlights include the discovery of a new population of non-thermal radio filaments; identification of new candidate supernova remnants, pulsar wind nebulae and planetary nebulae; improved radio/mid-IR classification of rare Luminous Blue Variables and discovery of associated extended radio nebulae; new radio stars identified by Bayesian cross-matching techniques; the realisation that many of the largest radio-quiet WISE HII region candidates are not true HII regions; and a large sample of previously undiscovered background HI galaxies in the Zone of Avoidance.
△ Less
Submitted 2 May, 2024; v1 submitted 12 December, 2023;
originally announced December 2023.
-
A Multi-Frequency View of the Radio Phoenix in the Abell 85 Cluster
Authors:
Ramij Raja,
Majidul Rahaman,
Abhirup Datta,
Oleg M. Smirnov
Abstract:
Radio phoenices are complex and filamentary diffuse radio sources found in both merging and relaxed clusters. The formation of these sources was proposed to be adiabatic compression of old Active Galactic Nucleus (AGN) plasma in shock waves. Most of the previous spectral studies of these sources were limited to integrated spectral indices, which were found to be very steep as well as show a curved…
▽ More
Radio phoenices are complex and filamentary diffuse radio sources found in both merging and relaxed clusters. The formation of these sources was proposed to be adiabatic compression of old Active Galactic Nucleus (AGN) plasma in shock waves. Most of the previous spectral studies of these sources were limited to integrated spectral indices, which were found to be very steep as well as show a curved spectrum. Here, we have performed a multi-frequency investigation of the radio phoenix in the Abell 85 cluster. Owing to the sensitive high-resolution observations, we found some of the finer filamentary structures that were previously undetected. We produced resolved spectral index maps of the radio phoenix between 148, 323, 700, and 1280 MHz. The orientation of the filaments, as well as the gradient across the spectral index maps suggest the possible direction of the shock motion from northeast to southwest. The integrated spectral index of the radio phoenix was found to be very steep with a break at around 700 MHz, indicating the re-energization of fossil electrons being recent. Furthermore, the spectral index of the filaments was found to be less steep compared to the non-filamentary regions, implying greater energy injection in the filaments. The observed features in the radio phoenix in the Abell 85 cluster seem to be in support of the adiabatic shock compression mechanism.
△ Less
Submitted 25 September, 2023;
originally announced September 2023.
-
Radio multifrequency observations of Abell~781 with the WSRT
Authors:
B. V. Hugo,
G. Bernardi,
O. M. Smirnov,
D. Dallacasa,
T. Venturi,
M. Murgia,
R. F. Pizzo
Abstract:
The `Main' galaxy cluster in the Abell 781 system is undergoing a significant merger and accretion process with peripheral emission to the north and southeastern flanks of the merging structure. Here we present a full polarimetric study of this field, using radio interferometric data taken at 21 and 92 cm with the Westerbork Synthesis Radio Telescope, to a sensitivity better than any 21 cm (L-band…
▽ More
The `Main' galaxy cluster in the Abell 781 system is undergoing a significant merger and accretion process with peripheral emission to the north and southeastern flanks of the merging structure. Here we present a full polarimetric study of this field, using radio interferometric data taken at 21 and 92 cm with the Westerbork Synthesis Radio Telescope, to a sensitivity better than any 21 cm (L-band) observation to date. We detect evidence of extended low-level emission of 1.9 mJy associated with the Main cluster at 21 cm, although this detection necessitates further follow-up by modern instruments due to the limited resolution of the Westerbork Synthesis Radio Telescope. Our polarimetric study indicates that, most likely, the peripheral emission associated with this cluster is not a radio relic.
△ Less
Submitted 18 September, 2023;
originally announced September 2023.
-
A radio bridge connecting the minihalo and phoenix in the Abell 85 cluster
Authors:
Ramij Raja,
Majidul Rahaman,
Abhirup Datta,
Oleg M. Smirnov
Abstract:
Galaxy clusters are located at the nodes of cosmic filaments and therefore host a lot of hydrodynamical activity. However, cool core clusters are considered to be relatively relaxed systems without much merging activity. The Abell 85 cluster is a unique example where the cluster hosts both a cool core and multiple ongoing merging processes. In this work, we used 700 MHz uGMRT as well as MeerKAT L-…
▽ More
Galaxy clusters are located at the nodes of cosmic filaments and therefore host a lot of hydrodynamical activity. However, cool core clusters are considered to be relatively relaxed systems without much merging activity. The Abell 85 cluster is a unique example where the cluster hosts both a cool core and multiple ongoing merging processes. In this work, we used 700 MHz uGMRT as well as MeerKAT L-band observations, carried out as part of the MGCLS, of the Abell 85. We reconfirm the presence of a minihalo in the cluster centre at 700MHz that was recently discovered in MGCLS. Furthermore, we discovered a radio bridge connecting the central minihalo and the peripheral radio phoenix. The mean surface brightness, size and flux density of the bridge at 700 MHz is found to be $\sim 0.14\ μ$Jy/arcsec$^2$, $\sim 220$ kpc and $\sim 4.88$ mJy, respectively, with a spectral index of $α_{700}^{1.28} = -0.92$. Although the origin of the seed relativistic electrons is still unknown, turbulent re-acceleration caused by both the spiralling sloshing gas in the intracluster medium (ICM) and the post-shock turbulence from the outgoing merging shock associated with the phoenix formation may be responsible for the bridge.
△ Less
Submitted 17 August, 2023;
originally announced August 2023.
-
MIGHTEE-\HI: Possible interactions with the galaxy NGC~895
Authors:
Brenda Namumba,
Javier Román,
Jesus Falcon Barroso,
Johan H. Knapen,
Ianjamasimanana Roger,
Elizabeth Naluminsa,
Gyula I. G. Jozsa,
Marie Korsaga,
Natasha Maddox,
Brad Frank,
Sinenhlanhla Sikhosana,
Samuel Legodi,
Claude Carignan,
Anastasia A. Ponomareva,
Tom Jarrett,
Danielle Lucero,
Oleg M. Smirnov,
Thijs van der Hulst,
D. J. Pisano,
kasia Malek,
Lucia Marchetti,
Mattia Vaccari,
Matt Jarvis,
Maarten Baes,
Martin Meyer
, et al. (7 additional authors not shown)
Abstract:
The transformation and evolution of a galaxy is strongly influenced by interactions with its environment. Neutral hydrogen (\HI) is an excellent way to trace these interactions. Here, we present \HI\ observations of the spiral galaxy NGC~895, which was previously thought to be isolated. High-sensitivity \HI\ observations from the MeerKAT large survey project MIGHTEE reveal possible interaction fea…
▽ More
The transformation and evolution of a galaxy is strongly influenced by interactions with its environment. Neutral hydrogen (\HI) is an excellent way to trace these interactions. Here, we present \HI\ observations of the spiral galaxy NGC~895, which was previously thought to be isolated. High-sensitivity \HI\ observations from the MeerKAT large survey project MIGHTEE reveal possible interaction features, such as extended spiral arms, and the two newly discovered \HI\ companions, that drive us to change the narrative that it is an isolated galaxy. We combine these observations with deep optical images from the Hyper Suprime Camera to show an absence of tidal debris between NGC 895 and its companions. We do find an excess of light in the outer parts of the companion galaxy MGTH$\_$J022138.1-052631 which could be an indication of external perturbation and thus possible sign of interactions. Our analysis shows that NGC~895 is an actively star-forming galaxy with a SFR of $\mathrm{1.75 \pm 0.09 [M_{\odot}/yr]}$, a value typical for high stellar mass galaxies on the star forming main sequence. It is reasonable to state that different mechanisms may have contributed to the observed features in NGC~895 and this emphasizes the need to revisit the target with more detailed observations. Our work shows the high potential and synergy of using state-of-the-art data in both \HI\ and optical to reveal a more complete picture of galaxy environments.
△ Less
Submitted 16 March, 2023;
originally announced March 2023.
-
Search for the Epoch of Reionisation with HERA: Upper Limits on the Closure Phase Delay Power Spectrum
Authors:
Pascal M. Keller,
Bojan Nikolic,
Nithyanandan Thyagarajan,
Chris L. Carilli,
Gianni Bernardi,
Ntsikelelo Charles,
Landman Bester,
Oleg M. Smirnov,
Nicholas S. Kern,
Joshua S. Dillon,
Bryna J. Hazelton,
Miguel F. Morales,
Daniel C. Jacobs,
Aaron R. Parsons,
Zara Abdurashidova,
Tyrone Adams,
James E. Aguirre,
Paul Alexander,
Zaki S. Ali,
Rushelle Baartman,
Yanga Balfour,
Adam P. Beardsley,
Tashalee S. Billings,
Judd D. Bowman,
Richard F. Bradley
, et al. (58 additional authors not shown)
Abstract:
Radio interferometers aiming to measure the power spectrum of the redshifted 21 cm line during the Epoch of Reionisation (EoR) need to achieve an unprecedented dynamic range to separate the weak signal from overwhelming foreground emissions. Calibration inaccuracies can compromise the sensitivity of these measurements to the effect that a detection of the EoR is precluded. An alternative to standa…
▽ More
Radio interferometers aiming to measure the power spectrum of the redshifted 21 cm line during the Epoch of Reionisation (EoR) need to achieve an unprecedented dynamic range to separate the weak signal from overwhelming foreground emissions. Calibration inaccuracies can compromise the sensitivity of these measurements to the effect that a detection of the EoR is precluded. An alternative to standard analysis techniques makes use of the closure phase, which allows one to bypass antenna-based direction-independent calibration. Similarly to standard approaches, we use a delay spectrum technique to search for the EoR signal. Using 94 nights of data observed with Phase I of the Hydrogen Epoch of Reionization Array (HERA), we place approximate constraints on the 21 cm power spectrum at $z=7.7$. We find at 95% confidence that the 21 cm EoR brightness temperature is $\le$(372)$^2$ "pseudo" mK$^2$ at 1.14 "pseudo" $h$ Mpc$^{-1}$, where the "pseudo" emphasises that these limits are to be interpreted as approximations to the actual distance scales and brightness temperatures. Using a fiducial EoR model, we demonstrate the feasibility of detecting the EoR with the full array. Compared to standard methods, the closure phase processing is relatively simple, thereby providing an important independent check on results derived using visibility intensities, or related.
△ Less
Submitted 15 February, 2023;
originally announced February 2023.
-
First AI for deep super-resolution wide-field imaging in radio astronomy: unveiling structure in ESO 137--006
Authors:
Arwa Dabbech,
Matthieu Terris,
Adrian Jackson,
Mpati Ramatsoku,
Oleg M. Smirnov,
Yves Wiaux
Abstract:
We introduce the first AI-based framework for deep, super-resolution, wide-field radio-interferometric imaging, and demonstrate it on observations of the ESO~137-006 radio galaxy. The algorithmic framework to solve the inverse problem for image reconstruction builds on a recent ``plug-and-play'' scheme whereby a denoising operator is injected as an image regulariser in an optimisation algorithm, w…
▽ More
We introduce the first AI-based framework for deep, super-resolution, wide-field radio-interferometric imaging, and demonstrate it on observations of the ESO~137-006 radio galaxy. The algorithmic framework to solve the inverse problem for image reconstruction builds on a recent ``plug-and-play'' scheme whereby a denoising operator is injected as an image regulariser in an optimisation algorithm, which alternates until convergence between denoising steps and gradient-descent data-fidelity steps. We investigate handcrafted and learned variants of high-resolution high-dynamic range denoisers. We propose a parallel algorithm implementation relying on automated decompositions of the image into facets and the measurement operator into sparse low-dimensional blocks, enabling scalability to large data and image dimensions. We validate our framework for image formation at a wide field of view containing ESO~137-006, from 19 gigabytes of MeerKAT data at 1053 and 1399 MHz. The recovered maps exhibit significantly more resolution and dynamic range than CLEAN, revealing collimated synchrotron threads close to the galactic core.
△ Less
Submitted 17 January, 2023; v1 submitted 22 July, 2022;
originally announced July 2022.
-
Tricolour: an optimized SumThreshold flagger for MeerKAT
Authors:
Benjamin V. Hugo,
Simon Perkins,
Bruce Merry,
Tom Mauch,
Oleg M. Smirnov
Abstract:
We present Tricolour, a package for Radio Frequency Interference mitigation of wideband finely channelized MeerKAT correlation data. The MeerKAT passband is heavily affected by interference from satellite, mobile, aircraft and terrestrial transponders. Coupled with typical data rates in excess of 100 GiB/hr at 208kHz channelization resolution, mitigation poses a significant processing challenge. O…
▽ More
We present Tricolour, a package for Radio Frequency Interference mitigation of wideband finely channelized MeerKAT correlation data. The MeerKAT passband is heavily affected by interference from satellite, mobile, aircraft and terrestrial transponders. Coupled with typical data rates in excess of 100 GiB/hr at 208kHz channelization resolution, mitigation poses a significant processing challenge. Our flagger is highly configurable, parallel and optimized, employing Dask and Numba technologies to implement the widely used SumThreshold and MAD interference detection algorithms. We find that typical 208kHz channelized datasets can be processed at rates in excess of 400 GiB/hr for a typical L-band flagging strategy on a modern dual-socket Intel Xeon server.
△ Less
Submitted 18 June, 2022;
originally announced June 2022.
-
The extended HI halo of NGC 4945 as seen by MeerKAT
Authors:
Roger Ianjamasimanana,
B. S. Koribalski,
Gyula I. G. Józsa,
Peter Kamphuis,
W. J. G. de Blok,
Dane Kleiner,
Brenda Namumba,
Claude Carignan,
Ralf-Jürgen Dettmar,
Paolo Serra,
Oleg M. Smirnov,
Kshitij Thorat,
Benjamin V. Hugo,
Athanaseus J. T. Ramaila,
Eric Maina,
Filippo M. Maccagni,
Sphesihle Makhathini,
Lexy A. L. Andati,
Dániel Cs. Molnár,
Simon Perkins,
Francesca Loi,
Mpati Ramatsoku,
Marcellin Atemkeng
Abstract:
Observations of the neutral atomic hydrogen (HI) in the nuclear starburst galaxy NGC 4945 with MeerKAT are presented. We find a large amount of halo gas, previously missed by HI observations, accounting for 6.8% of the total HI mass. This is most likely gas blown into the halo by star formation. Our maps go down to a $3σ$ column density level of $5\times10^{18} cm^{-2}$ . We model the HI distribut…
▽ More
Observations of the neutral atomic hydrogen (HI) in the nuclear starburst galaxy NGC 4945 with MeerKAT are presented. We find a large amount of halo gas, previously missed by HI observations, accounting for 6.8% of the total HI mass. This is most likely gas blown into the halo by star formation. Our maps go down to a $3σ$ column density level of $5\times10^{18} cm^{-2}$ . We model the HI distribution using tilted-ring fitting techniques and find a warp on the galaxy's approaching and receding sides. The HI in the northern side of the galaxy appears to be suppressed. This may be the result of ionisation by the starburst activity in the galaxy, as suggested by a previous study. The origin of the warp is unclear but could be due to past interactions or ram pressure stripping. Broad, asymmetric HI absorption lines extending beyond the HI emission velocity channels are present towards the nuclear region of NGC 4945. Such broad lines suggest the existence of a nuclear ring moving at a high circular velocity. This is supported by the clear rotation patterns in the HI absorption velocity field. The asymmetry of the absorption spectra can be caused by outflows or inflows of gas in the nuclear region of NGC 4945. The continuum map shows small extensions on both sides of the galaxy's major axis that might be signs of outflows resulting from the starburst activity.
△ Less
Submitted 1 April, 2022;
originally announced April 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.
-
The detection of a massive chain of dark HI clouds in the GAMA G23 Field
Authors:
Gyula I. G. Jozsa,
T. H. Jarrett,
Michelle Cluver,
O. Ivy Wong,
Okkert Havenga,
H. F. M. Yao,
L. Marchetti,
E. N. Taylor,
Peter Kamphuis,
Filippo M. Maccagni,
Athanaseus J. T. Ramaila,
Paolo Serra,
Oleg M. Smirnov,
Sarah V. White,
Virginia Kilborn,
B. W. Holwerda,
A. M. Hopkins,
S. Brough,
K. A. Pimbblet,
Simon P. Driver,
K. Kuijken
Abstract:
We report on the detection of a large, extended HI cloud complex in the GAMA G23 field, located at a redshift of $z\,\sim\,0.03$, observed as part of the MeerHOGS campaign (a pilot survey to explore the mosaicing capabilities of MeerKAT). The cloud complex, with a total mass of $10^{10.0}\,M_\odot$, lies in proximity to a large galaxy group with $M_\mathrm{dyn}\sim10^{13.5}\,M_\odot$. We identify…
▽ More
We report on the detection of a large, extended HI cloud complex in the GAMA G23 field, located at a redshift of $z\,\sim\,0.03$, observed as part of the MeerHOGS campaign (a pilot survey to explore the mosaicing capabilities of MeerKAT). The cloud complex, with a total mass of $10^{10.0}\,M_\odot$, lies in proximity to a large galaxy group with $M_\mathrm{dyn}\sim10^{13.5}\,M_\odot$. We identify seven HI peak concentrations, interconnected as a tenuous 'chain' structure, extending $\sim 400\,\mathrm{kpc}$ from east-to-west, with the largest (central) concentration containing $10{^{9.7}}\,M_\odot$ in HI gas distributed across $50\,\mathrm{kpc}$. The main source is not detected in ultra-violet, optical or infrared imaging. The implied gas mass-to-light ($M_\mathrm{HI}$/$L_\mathrm{r}$) is extreme ($>$1000) even in comparison to other 'dark clouds'. The complex has very little kinematic structure ($110\,\mathrm{km}\,\mathrm{s}^{-1}$), making it difficult to identify cloud rotation. Assuming pressure support, the total mass of the central concentration is $>10^{10.2}\,M_\odot$, while a lower limit to the dynamical mass in the case of full rotational support is $10^{10.4}\,M_\odot$. If the central concentration is a stable structure, it has to contain some amount of unseen matter, but potentially less than is observed for a typical galaxy. It is, however, not clear whether the structure has any gravitationally stable concentrations. We report a faint UV--optical--infrared source in proximity to one of the smaller concentrations in the gas complex, leading to a possible stellar association. The system nature and origins is enigmatic, potentially being the result of an interaction with or within the galaxy group it appears to be associated with.
△ Less
Submitted 6 December, 2021; v1 submitted 3 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.
-
The Hydrogen Intensity and Real-time Analysis eXperiment: 256-Element Array Status and Overview
Authors:
Devin Crichton,
Moumita Aich,
Adam Amara,
Kevin Bandura,
Bruce A. Bassett,
Carlos Bengaly,
Pascale Berner,
Shruti Bhatporia,
Martin Bucher,
Tzu-Ching Chang,
H. Cynthia Chiang,
Jean-Francois Cliche,
Carolyn Crichton,
Romeel Dave,
Dirk I. L. de Villiers,
Matt A. Dobbs,
Aaron M. Ewall-Wice,
Scott Eyono,
Christopher Finlay,
Sindhu Gaddam,
Ken Ganga,
Kevin G. Gayley,
Kit Gerodias,
Tim Gibbon,
Austin Gumba
, et al. (75 additional authors not shown)
Abstract:
The Hydrogen Intensity and Real-time Analysis eXperiment (HIRAX) is a radio interferometer array currently in development, with an initial 256-element array to be deployed at the South African Radio Astronomy Observatory (SARAO) Square Kilometer Array (SKA) site in South Africa. Each of the 6m, $f/0.23$ dishes will be instrumented with dual-polarisation feeds operating over a frequency range of 40…
▽ More
The Hydrogen Intensity and Real-time Analysis eXperiment (HIRAX) is a radio interferometer array currently in development, with an initial 256-element array to be deployed at the South African Radio Astronomy Observatory (SARAO) Square Kilometer Array (SKA) site in South Africa. Each of the 6m, $f/0.23$ dishes will be instrumented with dual-polarisation feeds operating over a frequency range of 400-800 MHz. Through intensity mapping of the 21 cm emission line of neutral hydrogen, HIRAX will provide a cosmological survey of the distribution of large-scale structure over the redshift range of $0.775 < z < 2.55$ over $\sim$15,000 square degrees of the southern sky. The statistical power of such a survey is sufficient to produce $\sim$7 percent constraints on the dark energy equation of state parameter when combined with measurements from the Planck satellite. Additionally, HIRAX will provide a highly competitive platform for radio transient and HI absorber science while enabling a multitude of cross-correlation studies. In this paper, we describe the science goals of the experiment, overview of the design and status of the sub-components of the telescope system, and describe the expected performance of the initial 256-element array as well as the planned future expansion to the final, 1024-element array.
△ Less
Submitted 17 January, 2022; v1 submitted 28 September, 2021;
originally announced September 2021.
-
A MeerKAT view on galaxy clusters: a radio-optical study of Abell 1300 and MACS J1931.8--2634
Authors:
B. Terni de Gregory,
B. Hugo,
T. Venturi,
G. Bernardi,
D. Dallacasa,
M. Nonino,
S. Makhatini,
V. Parekh,
O. M. Smirnov,
S. Giacintucci,
R. Kale
Abstract:
In this paper we present results from a radio-optical study of the galaxy populations of the galaxy clusters Abell 1300 and MACS J1931.8$-$2634, a merger and a relaxed system respectively both located at $z \sim 0.3$, aimed at finding evidence of merger-induced radio emission. Radio observations are taken at 1.28 GHz with the MeerKAT interferometer during its early-stage commissioning phase, and c…
▽ More
In this paper we present results from a radio-optical study of the galaxy populations of the galaxy clusters Abell 1300 and MACS J1931.8$-$2634, a merger and a relaxed system respectively both located at $z \sim 0.3$, aimed at finding evidence of merger-induced radio emission. Radio observations are taken at 1.28 GHz with the MeerKAT interferometer during its early-stage commissioning phase, and combined with archive optical data. We generated catalogues containing 107 and 162 radio sources in the A$~$1300 and MACS J1931.8--2634 cluster fields respectively, above a 0.2 mJy threshold and within a 30~arcmin radius from the cluster centre (corresponding to 8.1 and 8.8 Mpc respectively). By cross-correlating the radio and optical catalogues, and including spectroscopic information, 9 and 6 sources were found to be cluster members and used to construct the radio luminosity functions respectively for both clusters. The comparison of the radio source catalogues between the two cluster fields leads to a marginal difference, with a $2σ$ statistical significance. We derived the radio luminosity function at 1.28 GHz in both clusters, in the power range $22.81 < \rm {log~P_{1.28~GHz}~(W/Hz)} < 25.95$, and obtained that in A 1300 the radio luminosity function averaged over the full radio power interval is only $3.3 \pm 1.9$ times higher than the MACS J1931.8--2634 one, suggesting no statistical difference in their probability to host nuclear radio emission. We conclude that, at least for the two clusters studied here, the role of cluster mergers in affecting the statistical properties of the radio galaxy population is negligible.
△ Less
Submitted 12 April, 2021;
originally announced April 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.
-
The VLA Frontier Fields Survey: Deep, High-resolution Radio Imaging of the MACS Lensing Clusters at 3 and 6 GHz
Authors:
I. Heywood,
E. J. Murphy,
E. F. Jiménez-Andrade,
L. Armus,
W. D. Cotton,
C. DeCoursey,
M. Dickinson,
T. J. W. Lazio,
E. Momjian,
K. Penner,
I. Smail,
O. M. Smirnov
Abstract:
The Frontier Fields project is an observational campaign targeting six galaxy clusters, with the intention of using the magnification provided by gravitational lensing to study galaxies that are extremely faint or distant. We used the Karl G. Jansky Very Large Array (VLA) at 3 and 6 GHz to observe three Frontier Fields: MACSJ0416.1$-$2403 ($z$ = 0.396), MACSJ0717.5+3745 ($z$ = 0.545), and MACSJ114…
▽ More
The Frontier Fields project is an observational campaign targeting six galaxy clusters, with the intention of using the magnification provided by gravitational lensing to study galaxies that are extremely faint or distant. We used the Karl G. Jansky Very Large Array (VLA) at 3 and 6 GHz to observe three Frontier Fields: MACSJ0416.1$-$2403 ($z$ = 0.396), MACSJ0717.5+3745 ($z$ = 0.545), and MACSJ1149.5+2223 ($z$ = 0.543). The images reach noise levels of $\sim$1 $μ$Jy beam$^{-1}$ with sub-arcsecond resolution ($\sim$2.5 kpc at $z$ = 3), providing a high-resolution view of high-$z$ star-forming galaxies that is unbiased by dust obscuration. We generate dual-frequency continuum images at two different resolutions per band, per cluster, and derive catalogs totalling 1966 compact radio sources. Components within the areas of Hubble Space Telescope and Subaru observations are cross-matched, providing host galaxy identifications for 1296 of them. We detect 13 moderately-lensed (2.1 $<$ $μ$ $<$ 6.5) sources, one of which has a demagnified peak brightness of 0.9 $μ$Jy beam$^{-1}$, making it a candidate for the faintest radio source ever detected. There are 66 radio sources exhibiting complex morphologies, and 58 of these have host galaxy identifications. We reveal that MACSJ1149.5+2223 is not a cluster with a double relic, as the western candidate relic is resolved as a double-lobed radio galaxy associated with a foreground elliptical at $z$ = 0.24. The VLA Frontier Fields project is a public legacy survey. The image and catalog products from this work are freely available.
△ Less
Submitted 13 March, 2021;
originally announced March 2021.
-
Radio multifrequency observations of galaxy clusters. The Abell 399$-$401 pair
Authors:
C. D. Nunhokee,
G. Bernardi,
S. Manti,
F. Govoni,
A. Bonafede,
T. Venturi,
D. Dallacasa,
M. Murgia,
E. Orrú,
R. F. Pizzo,
O. M. Smirnov,
V. Vacca
Abstract:
Galaxy clusters are assembled via merging of smaller structures, in a process that generates shocks and turbulence in the intra cluster medium and produces radio diffuse emission in the form of halos and relics. The cluster pair A399-A401 represents a special case: both clusters host a radio halo. Recent Low Frequency Array (LOFAR) observations at 140 MHz revealed the presence of a radio bridge co…
▽ More
Galaxy clusters are assembled via merging of smaller structures, in a process that generates shocks and turbulence in the intra cluster medium and produces radio diffuse emission in the form of halos and relics. The cluster pair A399-A401 represents a special case: both clusters host a radio halo. Recent Low Frequency Array (LOFAR) observations at 140 MHz revealed the presence of a radio bridge connecting the two clusters along with two relic candidates. These relics include one South of A399 and the other in between the two clusters, in proximity of a shock front detected in X-ray observations. In this paper we present observations of the A399-A401 cluster pair at 1.7, 1.4, 1.2 GHz and 346 MHz from the Westerbork Synthesis Radio Telescope (WSRT). We detect the radio halo in the A399 cluster at 346 MHz, extending up to $\sim 650$ kpc and with a $125 \pm 6$ mJy flux density. Its spectral index between 140 MHz and 346 MHz is $α= 1.75 \pm 0.14$. The two candidate relics are also seen at 346 MHz and we determine their spectral indices to be $α= 1.10 \pm 0.14$ and $α= 1.46 \pm 0.14$. The low surface brightness bridge connecting the two clusters is below the noise level at 346 MHz, therefore we constrain the bridge average spectral index to be steep, i.e. $α> 1.5$ at 95% confidence level. This result favours the scenario where dynamically-induced turbulence is a viable mechanism to reaccelerate a population of mildly relativistic particles and amplify magnetic fields on scales of a few Mpcs.
Key words: galaxies: clusters: general - galaxies: clusters: individual: Abell 399 - radio continuum: general
△ Less
Submitted 15 May, 2023; v1 submitted 4 February, 2021;
originally announced February 2021.
-
Cygnus A jointly calibrated and imaged via non-convex optimisation from VLA data
Authors:
Arwa Dabbech,
Audrey Repetti,
Rick A. Perley,
Oleg M. Smirnov,
Yves Wiaux
Abstract:
Radio interferometric (RI) data are noisy under-sampled spatial Fourier components of the unknown radio sky affected by direction-dependent antenna gains. Failure to model these antenna gains accurately results in a radio sky estimate with limited fidelity and resolution. The RI inverse problem has been recently addressed via a joint calibration and imaging approach which consists in solving a non…
▽ More
Radio interferometric (RI) data are noisy under-sampled spatial Fourier components of the unknown radio sky affected by direction-dependent antenna gains. Failure to model these antenna gains accurately results in a radio sky estimate with limited fidelity and resolution. The RI inverse problem has been recently addressed via a joint calibration and imaging approach which consists in solving a non-convex minimisation task, involving suitable priors for the DDEs, namely temporal and spatial smoothness, and sparsity for the unknown radio map via an $\ell_1$-norm prior, in the context of realistic RI simulations. Building on these developments, we propose to promote sparsity of the radio map via a log-sum prior, enforcing sparsity more strongly than the $\ell_1$-norm. The resulting minimisation task is addressed via a sequence of non-convex minimisation tasks composed of re-weighted $\ell_1$ image priors, which are solved approximately. We demonstrate the efficiency of the approach on RI observations of the celebrated radio galaxy Cygnus~A obtained with the Karl G. Jansky Very Large Array at X, C, and S bands. More precisely, we showcase that the approach enhances data fidelity significantly while achieving high resolution high dynamic range radio maps, confirming the suitability of the priors considered for the unknown DDEs and radio image. As a clear qualitative indication of the high fidelity achieved by the data and the proposed approach, we report the detection of three background sources in the vicinity of Cyg~A, at S band.
△ Less
Submitted 12 May, 2021; v1 submitted 29 January, 2021;
originally announced February 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.
-
Anomalous gas in ESO 149-G003: A MeerKAT-16 View
Authors:
Gyula I. G. Józsa,
Kshitij Thorat,
Peter Kamphuis,
Lerato Sebokolodi,
Eric K. Maina,
Jing Wang,
Daniëlle L. A. Pieterse,
Paul Groot,
Athanaseus J. T. Ramaila,
Paolo Serra,
Lexy A. L. Andati,
W. J. G. de Blok,
Benjamin V. Hugo,
Dane Kleiner,
Filippo M. Maccagni,
Sphesihle Makhathini,
Dániel Cs. Molnár,
Mpati Ramatsoku,
Oleg M. Smirnov,
Steven Bloemen,
Kerry Paterson,
Paul Vreeswijk,
Vanessa McBride,
Marc Klein-Wolt,
Patrick Woudt
, et al. (6 additional authors not shown)
Abstract:
ESO 149-G003 is a close-by, isolated dwarf irregular galaxy. Previous observations with the ATCA indicated the presence of anomalous neutral hydrogen (HI) deviating from the kinematics of a regularly rotating disc. We conducted follow-up observations with the MeerKAT radio telescope during the 16-dish Early Science programme as well as with the MeerLICHT optical telescope. Our more sensitive radio…
▽ More
ESO 149-G003 is a close-by, isolated dwarf irregular galaxy. Previous observations with the ATCA indicated the presence of anomalous neutral hydrogen (HI) deviating from the kinematics of a regularly rotating disc. We conducted follow-up observations with the MeerKAT radio telescope during the 16-dish Early Science programme as well as with the MeerLICHT optical telescope. Our more sensitive radio observations confirm the presence of anomalous gas in ESO 149-G003, and further confirm the formerly tentative detection of an extraplanar HI component in the galaxy. Employing a simple tilted-ring model, in which the kinematics is determined with only four parameters but including morphological asymmetries, we reproduce the galaxy's morphology, which shows a high degree of asymmetry. By comparing our model with the observed HI, we find that in our model we cannot account for a significant (but not dominant) fraction of the gas. From the differences between our model and the observed data cube we estimate that at least 7%-8% of the HI in the galaxy exhibits anomalous kinematics, while we estimate a minimum mass fraction of less than 1% for the morphologically confirmed extraplanar component. We investigate a number of global scaling relations and find that, besides being gas-dominated with a neutral gas-to-stellar mass ratio of 1.7, the galaxy does not show any obvious global peculiarities. Given its isolation, as confirmed by optical observations, we conclude that the galaxy is likely currently acquiring neutral gas. It is either re-accreting gas expelled from the galaxy or accreting pristine intergalactic material.
△ Less
Submitted 6 January, 2021; v1 submitted 3 December, 2020;
originally announced December 2020.
-
MeerKAT-16 HI observation of the dIrr galaxy WLM
Authors:
Roger Ianjamasimanana,
Brenda Namumba,
Athanaseus J. T. Ramaila,
Anna S. Saburova,
Gyula I. G. Jozsa,
Talon Myburgh,
Kshitij Thorat,
Claude Carignan,
Eric Maina,
W. J. G. de Blok,
Lexy A. L. Andati,
Benjamin V. Hugo,
Dane Kleiner,
Peter Kamphuis,
Paolo Serra,
Oleg M. Smirnov,
Filippo M. Maccagni,
Sphesihle Makhathini,
Daniel Cs. Molnar,
Simon Perkins,
Mpati Ramatsoku,
Sarah V. White
Abstract:
We present observations and models of the kinematics and the distribution of the neutral hydrogen (HI) in the isolated dwarf irregular galaxy, Wolf-Lundmark-Melotte (WLM). We observed WLM with the Green Bank Telescope (GBT) and as part of the MeerKAT Early Science Programme, where 16 dishes were available. The HI disc of WLM extends out to a major axis diameter of 30 arcmin (8.5 kpc), and a minor…
▽ More
We present observations and models of the kinematics and the distribution of the neutral hydrogen (HI) in the isolated dwarf irregular galaxy, Wolf-Lundmark-Melotte (WLM). We observed WLM with the Green Bank Telescope (GBT) and as part of the MeerKAT Early Science Programme, where 16 dishes were available. The HI disc of WLM extends out to a major axis diameter of 30 arcmin (8.5 kpc), and a minor axis diameter of 20 arcmin (5.6 kpc) as measured by the GBT. We use the MeerKAT data to model WLM using the TiRiFiC software suite, allowing us to fit different tilted-ring models and select the one that best matches the observation. Our final best-fitting model is a flat disc with a vertical thickness, a constant inclination and dispersion, and a radially-varying surface brightness with harmonic distortions. To simulate bar-like motions, we include second-order harmonic distortions in velocity in the tangential and the vertical directions. We present a model with only circular motions included and a model with non-circular motions. The latter describes the data better. Overall, the models reproduce the global distribution and the kinematics of the gas, except for some faint emission at the 2-sigma level. We model the mass distribution of WLM with a pseudo-isothermal (ISO) and a Navarro-Frenk-White (NFW) dark matter halo models. The NFW and the ISO models fit the derived rotation curves within the formal errors, but with the ISO model giving better reduced chi-square values. The mass distribution in WLM is dominated by dark matter at all radii.
△ Less
Submitted 2 July, 2020;
originally announced July 2020.
-
VLA imaging of the XMM-LSS / VIDEO deep field at 1-2 GHz
Authors:
I. Heywood,
C. L. Hale,
M. J. Jarvis,
S. Makhathini,
J. A. Peters,
M. L. L. Sebokolodi,
O. M. Smirnov
Abstract:
Modern radio telescopes are routinely reaching depths where normal starforming galaxies are the dominant observed population. Realising the potential of radio as a tracer of star formation and black hole activity over cosmic time involves achieving such depths over representative volumes, with radio forming part of a larger multiwavelength campaign. In pursuit of this we used the Karl G. Jansky Ve…
▽ More
Modern radio telescopes are routinely reaching depths where normal starforming galaxies are the dominant observed population. Realising the potential of radio as a tracer of star formation and black hole activity over cosmic time involves achieving such depths over representative volumes, with radio forming part of a larger multiwavelength campaign. In pursuit of this we used the Karl G. Jansky Very Large Array (VLA) to image $\sim$5 deg$^{2}$ of the VIDEO/XMM-LSS extragalactic deep field at 1--2 GHz. We achieve a median depth of 16 $μ$Jy beam$^{-1}$ with an angular resolution of 4.5\arcsec. Comparisons with existing radio observations of XMM-LSS showcase the improved survey speed of the upgraded VLA: we cover 2.5 times the area and increase the depth by $\sim$20\% in 40\% of the time. Direction-dependent calibration and wide-field imaging were required to suppress the error patterns from off-axis sources of even modest brightness. We derive a catalogue containing 5,762 sources from the final mosaic. Sub-band imaging provides in-band spectral indices for 3,458 (60\%) sources, with the average spectrum becoming flatter than the canonical synchrotron slope below 1 mJy. Positional and flux-density accuracy of the observations, and the differential source counts are in excellent agreement with those of existing measurements. A public release of the images and catalogue accompanies this article.
△ Less
Submitted 15 June, 2020;
originally announced June 2020.
-
MeerKATHI -- an end-to-end data reduction pipeline for MeerKAT and other radio telescopes
Authors:
Gyula I. G. Józsa,
Sarah V. White,
Kshitij Thorat,
Oleg M. Smirnov,
Paolo Serra,
Mpati Ramatsoku,
Athanaseus J. T. Ramaila,
Simon J. Perkins,
Dániel Cs. Molnár,
Sphesihle Makhathini,
Filippo M. Maccagni,
Dane Kleiner,
Peter Kamphuis,
Benjamin V. Hugo,
W. J. G. de Blok,
Lexy A. L. Andati
Abstract:
MeerKATHI is the current development name for a radio-interferometric data reduction pipeline, assembled by an international collaboration. We create a publicly available end-to-end continuum- and line imaging pipeline for MeerKAT and other radio telescopes. We implement advanced techniques that are suitable for producing high-dynamic-range continuum images and spectroscopic data cubes. Using cont…
▽ More
MeerKATHI is the current development name for a radio-interferometric data reduction pipeline, assembled by an international collaboration. We create a publicly available end-to-end continuum- and line imaging pipeline for MeerKAT and other radio telescopes. We implement advanced techniques that are suitable for producing high-dynamic-range continuum images and spectroscopic data cubes. Using containerization, our pipeline is platform-independent. Furthermore, we are applying a standardized approach for using a number of different of advanced software suites, partly developed within our group. We aim to use distributed computing approaches throughout our pipeline to enable the user to reduce larger data sets like those provided by radio telescopes such as MeerKAT. The pipeline also delivers a set of imaging quality metrics that give the user the opportunity to efficiently assess the data quality.
△ Less
Submitted 5 June, 2020; v1 submitted 4 June, 2020;
originally announced June 2020.
-
Hydrodynamical Backflow in X-shaped Radio Galaxy PKS 2014-55
Authors:
W. D. Cotton,
K. Thorat,
J. J. Condon,
B . S. Frank,
G . I. G. Józsa,
S. V. White,
R . Deane,
N. Oozeer,
M. Atemkeng,
L. Bester,
B. Fanaroff,
R. S. Kupa,
O. M. Smirnov,
T. Mauch,
V. Krishnan,
F. Camilo
Abstract:
We present MeerKAT 1.28 GHz total-intensity, polarization, and spectral-index images covering the giant (projected length $l \approx 1.57$~Mpc) X-shaped radio source PKS~2014$-$55 with an unprecedented combination of brightness sensitivity and angular resolution. They show the clear "double boomerang" morphology of hydrodynamical backflows from the straight main jets deflected by the large and obl…
▽ More
We present MeerKAT 1.28 GHz total-intensity, polarization, and spectral-index images covering the giant (projected length $l \approx 1.57$~Mpc) X-shaped radio source PKS~2014$-$55 with an unprecedented combination of brightness sensitivity and angular resolution. They show the clear "double boomerang" morphology of hydrodynamical backflows from the straight main jets deflected by the large and oblique hot-gas halo of the host galaxy PGC~064440. The magnetic field orientation in PKS~2014$-$55 follows the flow lines from the jets through the secondary wings. The radio source is embedded in faint ($T_\mathrm{b} \approx 0.5 \mathrm{\,K}$) cocoons having the uniform brightness temperature and sharp outer edges characteristic of subsonic expansion into the ambient intra-group medium. The position angle of the much smaller ($l \sim 25$~kpc) restarted central source is within $5^\circ$ of the main jets, ruling out models that invoke jet re-orientation or two independent jets. Compression and turbulence in the backflows probably produce the irregular and low polarization bright region behind the apex of each boomerang as well as several features in the flow with bright heads and dark tails.
△ Less
Submitted 6 May, 2020;
originally announced May 2020.
-
Foreground modelling via Gaussian process regression: an application to HERA data
Authors:
Abhik Ghosh,
Florent Mertens,
Gianni Bernardi,
Mário G. Santos,
Nicholas S. Kern,
Christopher L. Carilli,
Trienko L. Grobler,
Léon V. E. Koopmans,
Daniel C. Jacobs,
Adrian Liu,
Aaron R. Parsons,
Miguel F. Morales,
James E. Aguirre,
Joshua S. Dillon,
Bryna J. Hazelton,
Oleg M. Smirnov,
Bharat K. Gehlot,
Siyanda Matika,
Paul Alexander,
Zaki S. Ali,
Adam P. Beardsley,
Roshan K. Benefo,
Tashalee S. Billings,
Judd D. Bowman,
Richard F. Bradley
, et al. (48 additional authors not shown)
Abstract:
The key challenge in the observation of the redshifted 21-cm signal from cosmic reionization is its separation from the much brighter foreground emission. Such separation relies on the different spectral properties of the two components, although, in real life, the foreground intrinsic spectrum is often corrupted by the instrumental response, inducing systematic effects that can further jeopardize…
▽ More
The key challenge in the observation of the redshifted 21-cm signal from cosmic reionization is its separation from the much brighter foreground emission. Such separation relies on the different spectral properties of the two components, although, in real life, the foreground intrinsic spectrum is often corrupted by the instrumental response, inducing systematic effects that can further jeopardize the measurement of the 21-cm signal. In this paper, we use Gaussian Process Regression to model both foreground emission and instrumental systematics in $\sim 2$ hours of data from the Hydrogen Epoch of Reionization Array. We find that a simple co-variance model with three components matches the data well, giving a residual power spectrum with white noise properties. These consist of an "intrinsic" and instrumentally corrupted component with a coherence-scale of 20 MHz and 2.4 MHz respectively (dominating the line of sight power spectrum over scales $k_{\parallel} \le 0.2$ h cMpc$^{-1}$) and a baseline dependent periodic signal with a period of $\sim 1$ MHz (dominating over $k_{\parallel} \sim 0.4 - 0.8$h cMpc$^{-1}$) which should be distinguishable from the 21-cm EoR signal whose typical coherence-scales is $\sim 0.8$ MHz.
△ Less
Submitted 12 May, 2020; v1 submitted 13 April, 2020;
originally announced April 2020.
-
A LOFAR Observation of Ionospheric Scintillation from Two Simultaneous Travelling Ionospheric Disturbances
Authors:
Richard A. Fallows,
Biagio Forte,
Ivan Astin,
Tom Allbrook,
Alex Arnold,
Alan Wood,
Gareth Dorrian,
Maaijke Mevius,
Hanna Rothkaehl,
Barbara Matyjasiak,
Andrzej Krankowski,
James M. Anderson,
Ashish Asgekar,
I. Max Avruch,
Mark Bentum,
Mario M. Bisi,
Harvey R. Butcher,
Benedetta Ciardi,
Bartosz Dabrowski,
Sieds Damstra,
Francesco de Gasperin,
Sven Duscha,
Jochen Eislöffel,
Thomas M. O. Franzen,
Michael A. Garrett
, et al. (33 additional authors not shown)
Abstract:
This paper presents the results from one of the first observations of ionospheric scintillation taken using the Low-Frequency Array (LOFAR). The observation was of the strong natural radio source Cas A, taken overnight on 18-19 August 2013, and exhibited moderately strong scattering effects in dynamic spectra of intensity received across an observing bandwidth of 10-80MHz. Delay-Doppler spectra (t…
▽ More
This paper presents the results from one of the first observations of ionospheric scintillation taken using the Low-Frequency Array (LOFAR). The observation was of the strong natural radio source Cas A, taken overnight on 18-19 August 2013, and exhibited moderately strong scattering effects in dynamic spectra of intensity received across an observing bandwidth of 10-80MHz. Delay-Doppler spectra (the 2-D FFT of the dynamic spectrum) from the first hour of observation showed two discrete parabolic arcs, one with a steep curvature and the other shallow, which can be used to provide estimates of the distance to, and velocity of, the scattering plasma. A cross-correlation analysis of data received by the dense array of stations in the LOFAR "core" reveals two different velocities in the scintillation pattern: a primary velocity of ~30m/s with a north-west to south-east direction, associated with the steep parabolic arc and a scattering altitude in the F-region or higher, and a secondary velocity of ~110m/s with a north-east to south-west direction, associated with the shallow arc and a scattering altitude in the D-region. Geomagnetic activity was low in the mid-latitudes at the time, but a weak sub-storm at high latitudes reached its peak at the start of the observation. An analysis of Global Navigation Satellite Systems (GNSS) and ionosonde data from the time reveals a larger-scale travelling ionospheric disturbance (TID), possibly the result of the high-latitude activity, travelling in the north-west to south-east direction, and, simultaneously, a smaller--scale TID travelling in a north-east to south-west direction, which could be associated with atmospheric gravity wave activity. The LOFAR observation shows scattering from both TIDs, at different altitudes and propagating in different directions. To the best of our knowledge this is the first time that such a phenomenon has been reported.
△ Less
Submitted 9 March, 2020;
originally announced March 2020.
-
The 1.28 GHz MeerKAT DEEP2 Image
Authors:
T. Mauch,
W. D. Cotton,
J. J. Condon,
A. M. Matthews,
T. D. Abbott,
R. M. Adam,
M. A. Aldera,
K. M. B. Asad,
E. F. Bauermeister,
T. G. H. Bennett,
H. Bester,
D. H. Botha,
L. R. S. Brederode,
Z. B. Brits,
S. J. Buchner,
J. P. Burger,
F. Camilo,
J. M. Chalmers,
T. Cheetham,
D. de Villiers,
M. S. de Villiers,
M. A. Dikgale-Mahlakoana,
L. J. du Toit,
S. W. P. Esterhuyse,
G. Fadana
, et al. (79 additional authors not shown)
Abstract:
We present the confusion-limited 1.28 GHz MeerKAT DEEP2 image covering one $\approx 68'$ FWHM primary beam area with $7.6''$ FWHM resolution and $0.55 \pm 0.01$ $μ$Jy/beam rms noise. Its J2000 center position $α=04^h 13^m 26.4^s$, $δ=-80^\circ 00' 00''$ was selected to minimize artifacts caused by bright sources. We introduce the new 64-element MeerKAT array and describe commissioning observations…
▽ More
We present the confusion-limited 1.28 GHz MeerKAT DEEP2 image covering one $\approx 68'$ FWHM primary beam area with $7.6''$ FWHM resolution and $0.55 \pm 0.01$ $μ$Jy/beam rms noise. Its J2000 center position $α=04^h 13^m 26.4^s$, $δ=-80^\circ 00' 00''$ was selected to minimize artifacts caused by bright sources. We introduce the new 64-element MeerKAT array and describe commissioning observations to measure the primary beam attenuation pattern, estimate telescope pointing errors, and pinpoint $(u,v)$ coordinate errors caused by offsets in frequency or time. We constructed a 1.4 GHz differential source count by combining a power-law count fit to the DEEP2 confusion $P(D)$ distribution from $0.25$ to $10$ $μ$Jy with counts of individual DEEP2 sources between $10$ $μ$Jy and $2.5$ mJy. Most sources fainter than $S \sim 100$ $μ$Jy are distant star-forming galaxies obeying the FIR/radio correlation, and sources stronger than $0.25$ $μ$Jy account for $\sim93\%$ of the radio background produced by star-forming galaxies. For the first time, the DEEP2 source count has reached the depth needed to reveal the majority of the star formation history of the universe. A pure luminosity evolution of the 1.4 GHz local luminosity function consistent with the Madau & Dickinson (2014) model for the evolution of star-forming galaxies based on UV and infrared data underpredicts our 1.4 GHz source count in the range $-5 \lesssim \log[S(\mathrm{Jy})] \lesssim -4$.
△ Less
Submitted 12 December, 2019;
originally announced December 2019.
-
Inflation of 430-parsec bipolar radio bubbles in the Galactic Centre by an energetic event
Authors:
I. Heywood,
F. Camilo,
W. D. Cotton,
F. Yusef-Zadeh,
T. D. Abbott,
R. M. Adam,
M. A. Aldera,
E. F. Bauermeister,
R. S. Booth,
A. G. Botha,
D. H. Botha,
L. R. S. Brederode,
Z. B. Brits,
S. J. Buchner,
J. P. Burger,
J. M. Chalmers,
T. Cheetham,
D. de Villiers,
M. A. Dikgale-Mahlakoana,
L. J. du Toit,
S. W. P. Esterhuyse,
B. L. Fanaroff,
A. R. Foley,
D. J. Fourie,
R. R. G. Gamatham
, et al. (74 additional authors not shown)
Abstract:
The Galactic Centre contains a supermassive black hole with a mass of 4 million suns within an environment that differs markedly from that of the Galactic disk. While the black hole is essentially quiescent in the broader context of active galactic nuclei, X-ray observations have provided evidence for energetic outbursts from its surroundings. Also, while the levels of star formation in the Galact…
▽ More
The Galactic Centre contains a supermassive black hole with a mass of 4 million suns within an environment that differs markedly from that of the Galactic disk. While the black hole is essentially quiescent in the broader context of active galactic nuclei, X-ray observations have provided evidence for energetic outbursts from its surroundings. Also, while the levels of star formation in the Galactic Centre have been approximately constant over the last few hundred Myr, there is evidence of elevated short-duration bursts, strongly influenced by interaction of the black hole with the enhanced gas density present within the ring-like Central Molecular Zone at Galactic longitude |l| < 0.7 degrees and latitude |b| < 0.2 degrees. The inner 200 pc region is characterized by large amounts of warm molecular gas, a high cosmic ray ionization rate, unusual gas chemistry, enhanced synchrotron emission, and a multitude of radio-emitting magnetised filaments, the origin of which has not been established. Here we report radio imaging that reveals bipolar bubbles spanning 1 degree x 3 degrees (140 parsecs x 430 parsecs), extending above and below the Galactic plane and apparently associated with the Galactic Centre. The structure is edge-brightened and bounded, with symmetry implying creation by an energetic event in the Galactic Centre. We estimate the age of the bubbles to be a few million years, with a total energy of 7 x 10^52 ergs. We postulate that the progenitor event was a major contributor to the increased cosmic-ray density in the Galactic Centre, and is in turn the principal source of the relativistic particles required to power the synchrotron emission of the radio filaments within and in the vicinity of the bubble cavities.
△ Less
Submitted 12 September, 2019;
originally announced September 2019.
-
On the nature of the unusual transient AT 2018cow from HI observations of its host galaxy
Authors:
Michał J. Michałowski,
P. Kamphuis,
J. Hjorth,
D. A. Kann,
A. de Ugarte Postigo,
L. Galbany,
J. P. U. Fynbo,
A. Ghosh,
L. K. Hunt,
H. Kuncarayakti,
E. Le Floc'h,
A. Leśniewska,
K. Misra,
A. Nicuesa Guelbenzu,
E. Palazzi,
J. Rasmussen,
L. Resmi,
A. Rossi,
S. Savaglio,
P. Schady,
S. Schulze,
C. C. Thöne,
D. Watson,
G. I. G. Józsa,
P. Serra
, et al. (1 additional authors not shown)
Abstract:
Unusual stellar explosions represent an opportunity to learn about both stellar and galaxy evolution. Mapping the atomic gas in host galaxies of such transients can lead to an understanding of the conditions triggering them. We provide resolved atomic gas observations of the host galaxy, CGCG137-068, of the unusual, poorly-understood transient AT2018cow searching for clues to understand its nature…
▽ More
Unusual stellar explosions represent an opportunity to learn about both stellar and galaxy evolution. Mapping the atomic gas in host galaxies of such transients can lead to an understanding of the conditions triggering them. We provide resolved atomic gas observations of the host galaxy, CGCG137-068, of the unusual, poorly-understood transient AT2018cow searching for clues to understand its nature. We test whether it is consistent with a recent inflow of atomic gas from the intergalactic medium, as suggested for host galaxies of gamma-ray bursts (GRBs) and some supernovae (SNe). We observed the HI hyperfine structure line of the AT2018cow host with the Giant Metrewave Radio Telescope. There is no unusual atomic gas concentration near the position of AT2018cow. The gas distribution is much more regular than those of GRB/SN hosts. The AT2018cow host has an atomic gas mass lower by 0.24 dex than predicted from its star formation rate (SFR) and is at the lower edge of the galaxy main sequence. In the continuum we detected the emission of AT2018cow and of a star-forming region in the north-eastern part of the bar (away from AT2018cow). This region hosts a third of the galaxy's SFR. The absence of atomic gas concentration close to AT2018cow, along with a normal SFR and regular HI velocity field, sets CGCG137-068 apart from GRB/SN hosts studied in HI. The environment of AT2018cow therefore suggests that its progenitor may not have been a massive star. Our findings are consistent with an origin of the transient that does not require a connection between its progenitor and gas concentration or inflow: an exploding low-mass star, a tidal disruption event, a merger of white dwarfs, or a merger between a neutron star and a giant star. We interpret the recently reported atomic gas ring in CGCG137-068 as a result of internal processes connected with gravitational resonances caused by the bar.
△ Less
Submitted 12 June, 2019; v1 submitted 26 February, 2019;
originally announced February 2019.
-
CubiCal - Fast radio interferometric calibration suite exploiting complex optimisation
Authors:
J. S. Kenyon,
O. M. Smirnov,
T. L. Grobler,
S. J. Perkins
Abstract:
It has recently been shown that radio interferometric gain calibration can be expressed succinctly in the language of complex optimisation. In addition to providing an elegant framework for further development, it exposes properties of the calibration problem which can be exploited to accelerate traditional non-linear least squares solvers such as Gauss-Newton and Levenberg-Marquardt. We extend ex…
▽ More
It has recently been shown that radio interferometric gain calibration can be expressed succinctly in the language of complex optimisation. In addition to providing an elegant framework for further development, it exposes properties of the calibration problem which can be exploited to accelerate traditional non-linear least squares solvers such as Gauss-Newton and Levenberg-Marquardt. We extend existing derivations to chains of Jones terms: products of several gains which model different aberrant effects. In doing so, we find that the useful properties found in the single term case still hold. We also develop several specialised solvers which deal with complex gains parameterised by real values. The newly developed solvers have been implemented in a Python package called CubiCal, which uses a combination of Cython, multiprocessing and shared memory to leverage the power of modern hardware. We apply CubiCal to both simulated and real data, and perform both direction-independent and direction-dependent self-calibration. Finally, we present the results of some rudimentary profiling to show that CubiCal is competitive with respect to existing calibration tools such as MeqTrees.
△ Less
Submitted 15 May, 2018; v1 submitted 9 May, 2018;
originally announced May 2018.
-
Redundant interferometric calibration as a complex optimization problem
Authors:
T. L. Grobler,
G. Bernardi,
J. S. Kenyon,
A. R. Parsons,
O. M. Smirnov
Abstract:
Observations of the redshifted 21-cm line from the epoch of reionization have recently motivated the construction of low frequency radio arrays with highly redundant configurations. These configurations provide an alternative calibration strategy - "redundant calibration" - and boosts sensitivity on specific spatial scales. In this paper, we formulate calibration of redundant interferometric array…
▽ More
Observations of the redshifted 21-cm line from the epoch of reionization have recently motivated the construction of low frequency radio arrays with highly redundant configurations. These configurations provide an alternative calibration strategy - "redundant calibration" - and boosts sensitivity on specific spatial scales. In this paper, we formulate calibration of redundant interferometric arrays as a complex optimization problem. We solve this optimization problem via the Levenberg-Marquardt algorithm. This calibration approach is more robust to initial conditions than current algorithms and, by leveraging an approximate matrix inversion, allows for further optimization and an efficient implementation ("redundant StEfCal"). We also investigated using the preconditioned conjugate gradient method as an alternative to the approximate matrix inverse, but found that its computational performance is not competitive with respect to "redundant StEfCal". The efficient implementation of this new algorithm is made publicly available.
△ Less
Submitted 9 February, 2018; v1 submitted 10 January, 2018;
originally announced January 2018.
-
Cygnus A super-resolved via convex optimisation from VLA data
Authors:
Arwa Dabbech,
Alexandru Onose,
Abdullah Abdulaziz,
Richard A. Perley,
Oleg M. Smirnov,
Yves Wiaux
Abstract:
We leverage the Sparsity Averaging Reweighted Analysis (SARA) approach for interferometric imaging, that is based on convex optimisation, for the super-resolution of Cyg A from observations at the frequencies 8.422GHz and 6.678GHz with the Karl G. Jansky Very Large Array (VLA). The associated average sparsity and positivity priors enable image reconstruction beyond instrumental resolution. An adap…
▽ More
We leverage the Sparsity Averaging Reweighted Analysis (SARA) approach for interferometric imaging, that is based on convex optimisation, for the super-resolution of Cyg A from observations at the frequencies 8.422GHz and 6.678GHz with the Karl G. Jansky Very Large Array (VLA). The associated average sparsity and positivity priors enable image reconstruction beyond instrumental resolution. An adaptive Preconditioned Primal-Dual algorithmic structure is developed for imaging in the presence of unknown noise levels and calibration errors. We demonstrate the superior performance of the algorithm with respect to the conventional CLEAN-based methods, reflected in super-resolved images with high fidelity. The high resolution features of the recovered images are validated by referring to maps of Cyg A at higher frequencies, more precisely 17.324GHz and 14.252GHz. We also confirm the recent discovery of a radio transient in Cyg A, revealed in the recovered images of the investigated data sets. Our matlab code is available online on GitHub.
△ Less
Submitted 14 May, 2018; v1 submitted 24 October, 2017;
originally announced October 2017.
-
A MeerKAT view on galaxy clusters
Authors:
G. Bernardi,
T. Venturi,
R. Cassano,
G. Brunetti,
D. Dallacasa,
B. Fanaroff,
B. Hugo,
S. Makhatini,
N. Oozeer,
O. M. Smirnov,
J. T. L. Zwart
Abstract:
Almost two decades of observations of radio emission in galaxy clusters have proven the existence of relativistic particles and magnetic fields that generate extended synchrotron emission in the form of radio halos. In the current scenario, radio halos are generated through re--acceleration of relativistic electrons by turbulence generated by cluster mergers. Although this theoretical framework ha…
▽ More
Almost two decades of observations of radio emission in galaxy clusters have proven the existence of relativistic particles and magnetic fields that generate extended synchrotron emission in the form of radio halos. In the current scenario, radio halos are generated through re--acceleration of relativistic electrons by turbulence generated by cluster mergers. Although this theoretical framework has received increasingly supporting observational evidence over the last ten years, observations of statistically complete samples are needed in order to fundamentally test model predictions. In this paper we briefly review our 7--element Karoo Radio Telescope observations of a sample of nearby clusters aimed to test the predictions of the turbulent re--acceleration model in small systems ($M_{500} > 4 \times 10^{14}$ M$_{\odot}$). We conclude by presenting two galaxy cluster surveys to be carried out with MeerKAT in order to provide crucial test of models of radio halo formation in nearby ($z < 0.1$) and high redshift ($z > 0.4$) systems respectively.
△ Less
Submitted 25 August, 2017;
originally announced August 2017.
-
Using baseline-dependent window functions for data compression and field-of-interest shaping in radio interferometry
Authors:
M. T. Atemkeng,
O. M. Smirnov,
C. Tasse,
G. Foster,
J. Jonas
Abstract:
In radio interferometry, observed visibilities are intrinsically sampled at some interval in time and frequency. Modern interferometers are capable of producing data at very high time and frequency resolution; practical limits on storage and computation costs require that some form of data compression be imposed. The traditional form of compression is a simple averaging of the visibilities over co…
▽ More
In radio interferometry, observed visibilities are intrinsically sampled at some interval in time and frequency. Modern interferometers are capable of producing data at very high time and frequency resolution; practical limits on storage and computation costs require that some form of data compression be imposed. The traditional form of compression is a simple averaging of the visibilities over coarser time and frequency bins. This has an undesired side effect: the resulting averaged visibilities "decorrelate", and do so differently depending on the baseline length and averaging interval. This translates into a non-trivial signature in the image domain known as "smearing", which manifests itself as an attenuation in amplitude towards off-centre sources. With the increasing fields of view and/or longer baselines employed in modern and future instruments, the trade-off between data rate and smearing becomes increasingly unfavourable. In this work we investigate alternative approaches to low-loss data compression. We show that averaging of the visibility data can be treated as a form of convolution by a boxcar-like window function, and that by employing alternative baseline-dependent window functions a more optimal interferometer smearing response may be induced. In particular, we show improved amplitude response over a chosen field of interest, and better attenuation of sources outside the field of interest. The main cost of this technique is a reduction in nominal sensitivity; we investigate the smearing vs. sensitivity trade-off, and show that in certain regimes a favourable compromise can be achieved. We show the application of this technique to simulated data from the Karl G. Jansky Very Large Array (VLA) and the European Very-long-baseline interferometry Network (EVN).
△ Less
Submitted 14 July, 2016;
originally announced July 2016.
-
Calibration artefacts in radio interferometry. III. Phase-only calibration and primary beam correction
Authors:
T. L. Grobler,
A. J. Stewart,
S. J. Wijnholds,
J. S. Kenyon,
O. M. Smirnov
Abstract:
This is the third installment in a series of papers in which we investigate calibration artefacts. Calibration artefacts (also known as ghosts or spurious sources) are created when we calibrate with an incomplete model. In the first two papers of this series we developed a mathematical framework which enabled us to study the ghosting mechanism itself. An interesting concomitant of the second paper…
▽ More
This is the third installment in a series of papers in which we investigate calibration artefacts. Calibration artefacts (also known as ghosts or spurious sources) are created when we calibrate with an incomplete model. In the first two papers of this series we developed a mathematical framework which enabled us to study the ghosting mechanism itself. An interesting concomitant of the second paper was that ghosts appear in symmetrical pairs. This could possibly account for spurious symmetrization. Spurious symmetrization refers to the appearance of a spurious source (the anti-ghost) symmetrically opposite an unmodelled source around a modelled source. The analysis in the first two papers indicates that the anti-ghost is usually very faint, in particular when a large number of antennas are used. This suggests that spurious symmetrization will mainly occur at an almost undetectable flux level. In this paper, we show that phase-only calibration produces an anti-ghost that is $N$-times (where $N$ denotes the number of antennas in the array) as bright as the one produced by phase and amplitude calibration and that this already bright ghost can be further amplified by the primary beam correction.
△ Less
Submitted 18 July, 2016; v1 submitted 20 June, 2016;
originally announced June 2016.
-
A deep / wide 1-2 GHz snapshot survey of SDSS Stripe 82 using the Karl G. Jansky Very Large Array in a compact hybrid configuration
Authors:
I. Heywood,
M. J. Jarvis,
A. J. Baker,
K. W. Bannister,
C. S. Carvalho,
M. Hardcastle,
M. Hilton,
K. Moodley,
O. M. Smirnov,
D. J. B. Smith,
S. V. White,
E. J. Wollack
Abstract:
We have used the Karl G. Jansky Very Large Array to image ~100 sq. deg. of SDSS Stripe 82 at 1-2 GHz. The survey consists of 1,026 snapshot observations of 2.5 minutes duration, using the hybrid CnB configuration. The survey has good sensitivity to diffuse, low surface brightness structures and extended radio emission, making it highly synergistic with existing 1.4 GHz radio observations of the re…
▽ More
We have used the Karl G. Jansky Very Large Array to image ~100 sq. deg. of SDSS Stripe 82 at 1-2 GHz. The survey consists of 1,026 snapshot observations of 2.5 minutes duration, using the hybrid CnB configuration. The survey has good sensitivity to diffuse, low surface brightness structures and extended radio emission, making it highly synergistic with existing 1.4 GHz radio observations of the region. The principal data products are continuum images, with 16 x 10 arcsecond resolution, and a catalogue containing 11,782 point and Gaussian components resulting from fits to the thresholded Stokes-I brightness distribution, forming approximately 8,948 unique radio sources. The typical effective 1σ noise level is 88 μJy / beam. Spectral index estimates are included, as derived from the 1 GHz of instantaneous bandwidth. Astrometric and photometric accuracy are in excellent agreement with existing narrowband observations. A large-scale simulation is used to investigate clean bias, which we extend into the spectral domain. Clean bias remains an issue for snapshot surveys with the VLA, affecting our total intensity measurements at the ~1σ level. Statistical spectral index measurements are in good agreement with existing measurements derived from matching separate surveys at two frequencies. At flux densities below ~35σ the median in-band spectral index measurements begin to exhibit a bias towards flatness that is dependent on both flux density and the intrinsic spectral index. In-band spectral curvature measurements are likely to be unreliable for all but the very brightest components. Image products and catalogues are publicly available via an FTP server.
△ Less
Submitted 29 May, 2016;
originally announced May 2016.
-
A KAT-7 view of a low-mass sample of galaxy clusters
Authors:
G. Bernardi,
T. Venturi,
R. Cassano,
D. Dallacasa,
G. Brunetti,
V. Cuciti,
M. Johnston-Hollitt,
N. Oozeer,
O. M. Smirnov
Abstract:
Radio observations over the last two decades have provided evidence that diffuse synchrotron emission in the form of megaparsec-scale radio halos in galaxy clusters is likely tracing regions of the intracluster medium where relativistic particles are accelerated during cluster mergers. In this paper we present results of a survey of 14 galaxy clusters carried out with the 7-element Karoo Array Tel…
▽ More
Radio observations over the last two decades have provided evidence that diffuse synchrotron emission in the form of megaparsec-scale radio halos in galaxy clusters is likely tracing regions of the intracluster medium where relativistic particles are accelerated during cluster mergers. In this paper we present results of a survey of 14 galaxy clusters carried out with the 7-element Karoo Array Telescope at 1.86 GHz, aimed to extend the current studies of radio halo occurrence to systems with lower masses (M$_{\rm 500} > 4\times10^{14}$ M${_\odot}$). We found upper limits at the $0.6 - 1.9 \times 10^{24}$ Watt Hz$^{-1}$ level for $\sim 50\%$ of the sample, confirming that bright radio halos in less massive galaxy clusters are statistically rare.
△ Less
Submitted 24 March, 2016;
originally announced March 2016.
-
KAT-7 observations of a mass-selected sample of galaxy clusters
Authors:
G. Bernardi,
T. Venturi,
R. Cassano,
D. Dallacasa,
G. Brunetti,
V. Cuciti,
M. Johnston-Hollitt,
N. Oozeer,
V. Parekh,
O. M. Smirnov
Abstract:
The presence of megaparsec-scale radio halos in galaxy clusters has already been established by many observations over the last two decades. The emerging explanation for the formation of these giant sources of diffuse synchrotron radio emission is that they trace turbulent regions in the intracluster medium, where particles are trapped and accelerated during cluster mergers. Our current observatio…
▽ More
The presence of megaparsec-scale radio halos in galaxy clusters has already been established by many observations over the last two decades. The emerging explanation for the formation of these giant sources of diffuse synchrotron radio emission is that they trace turbulent regions in the intracluster medium, where particles are trapped and accelerated during cluster mergers. Our current observational knowledge is, however, mainly limited to massive systems. Here we present observations of a sample of 14 mass-selected galaxy clusters, i.e. $M_{\rm 500} > 4\times10^{14}$~M${_\odot}$, in the Southern Hemisphere, aimed to study the occurrence of radio halos in low mass clusters and test the correlation between the radio halo power at 1.4 GHz $P_{\rm 1.4}$ and the cluster mass $M_{\rm 500}$. Our observations were performed with the 7-element Karoo Array Telescope at 1.86 GHz. We found three candidates to host diffuse cluster-scale emission and derived upper limits at the level of $0.6 - 1.9 \times 10^{24}$~Watt~Hz$^{-1}$ for $\sim 50\%$ of the clusters in the sample, significantly increasing the number of clusters with radio halo information in the considered mass range. Our results confirm that bright radio halos in less massive galaxy clusters are statistically rare.
△ Less
Submitted 3 December, 2015; v1 submitted 3 November, 2015;
originally announced November 2015.
-
Generalized Formalisms of the Radio Interferometer Measurement Equation
Authors:
Danny C. Price,
Oleg M. Smirnov
Abstract:
The Radio Interferometer Measurement Equation (RIME) is a matrix-based mathematical model that describes the response of a radio interferometer. The Jones calculus it employs is not suitable for describing the analogue components of a telescope. This is because it does not consider the effect of impedance mismatches between components. This paper aims to highlight the limitations of Jones calculus…
▽ More
The Radio Interferometer Measurement Equation (RIME) is a matrix-based mathematical model that describes the response of a radio interferometer. The Jones calculus it employs is not suitable for describing the analogue components of a telescope. This is because it does not consider the effect of impedance mismatches between components. This paper aims to highlight the limitations of Jones calculus, and suggests some alternative methods that are more applicable. We reformulate the RIME with a different basis that includes magnetic and mixed coherency statistics. We present a microwave network inspired 2N-port version of the RIME, and a tensor formalism based upon the electromagnetic tensor from special relativity. We elucidate the limitations of the Jones-matrix-based RIME for describing analogue components. We show how measured scattering parameters of analogue components can be used in a 2N-port version of the RIME. In addition, we show how motion at relativistic speed affects the observed flux. We present reformulations of the RIME that correctly account for magnetic field coherency. These reformulations extend the standard formulation, highlight its limitations, and may have applications in space-based interferometry and precise absolute calibration experiments.
△ Less
Submitted 26 January, 2015;
originally announced January 2015.
-
Calibration artefacts in radio interferometry. I. Ghost sources in WSRT data
Authors:
T. G. Grobler,
C. D. Nunhokee,
O. M. Smirnov,
A. J. van Zyl,
A. G. de Bruyn
Abstract:
This work investigates a particular class of artefacts, or ghost sources, in radio interferometric images. Earlier observations with (and simulations of) the Westerbork Synthesis Radio Telescope (WSRT) suggested that these were due to calibration with incomplete sky models. A theoretical framework is derived that validates this suggestion, and provides predictions of ghost formation in a two-sourc…
▽ More
This work investigates a particular class of artefacts, or ghost sources, in radio interferometric images. Earlier observations with (and simulations of) the Westerbork Synthesis Radio Telescope (WSRT) suggested that these were due to calibration with incomplete sky models. A theoretical framework is derived that validates this suggestion, and provides predictions of ghost formation in a two-source scenario. The predictions are found to accurately match the result of simulations, and qualitatively reproduce the ghosts previously seen in observational data. The theory also provides explanations for many previously puzzling features of these artefacts (regular geometry, PSF-like sidelobes, seeming independence on model flux), and shows that the observed phenomenon of flux suppression affecting unmodelled sources is due to the same mechanism. We demonstrate that this ghost formation mechanism is a fundamental feature of calibration, and exhibits a particularly strong and localized signature due to array redundancy. To some extent this mechanism will affect all observations (including those with non-redundant arrays), though in most cases the ghosts remain hidden below the noise or masked by other instrumental artefacts. The implications of such errors on future deep observations are discussed.
△ Less
Submitted 6 February, 2014;
originally announced February 2014.
-
Sub-millimetre source identifications and the micro-Jansky source population at 8.4 GHz in the William Herschel Deep Field
Authors:
I. Heywood,
R. M. Bielby,
M. D. Hill,
N. Metcalfe,
S. Rawlings,
T. Shanks,
O. M. Smirnov
Abstract:
[Abridged] Sub-mm observations of the William Herschel Deep Field using LABOCA revealed possible counterparts for 2 X-ray absorbed QSOs. The aim here is to exploit EVLA imaging at 8.4 GHz to establish the QSOs as radio/sub-mm sources. The challenge in reducing the EVLA data was the presence of a strong 4C source in the field. A new calibration algorithm was applied to the data to subtract it. The…
▽ More
[Abridged] Sub-mm observations of the William Herschel Deep Field using LABOCA revealed possible counterparts for 2 X-ray absorbed QSOs. The aim here is to exploit EVLA imaging at 8.4 GHz to establish the QSOs as radio/sub-mm sources. The challenge in reducing the EVLA data was the presence of a strong 4C source in the field. A new calibration algorithm was applied to the data to subtract it. The resulting thermal noise limited radio map covers the 16'x16' Extended WHDF. It contains 41 sources above a 4-sigma limit, 17 of which have primary beam corrected flux. The radio observations show that the absorbed AGN with LABOCA detections are coincident with radio sources, confirming the tendency for X-ray absorbed AGN to be sub-mm bright. These sources show strong ultraviolet excess (UVX) suggesting the nuclear sightline is gas- but not dust-absorbed. Of the 3 remaining LABOCA sources within the ~5' half-power beam width, 1 is identified with a faint nuclear X-ray/radio source in a nearby galaxy, 1 with a faint radio source and 1 is unidentified in any other band.
More generally, differential radio source counts are in good agreement with previous observations, showing at S<50 micro-Jy a significant excess over a pure AGN model. In the full area, of 10 sources fainter than this limit, 6 have optical counterparts of which 3 are UVX (i.e. likely QSOs) including the 2 absorbed quasar LABOCA sources. The other faint radio counterparts are not UVX but are only slightly less blue and likely to be star-forming/merging galaxies, predominantly at lower luminosities and redshifts. The 4 faint, optically unidentified radio sources may be either dust obscured QSOs or galaxies. These high-z obscured AGN and lower-z star-forming populations are thus the main candidates to explain the observed excess in faint source counts and hence the excess radio background found previously by the ARCADE2 experiment.
△ Less
Submitted 20 September, 2012;
originally announced September 2012.
-
MeqTrees and direction-dependent effects
Authors:
O. M. Smirnov,
A. G. de Bruyn
Abstract:
Direction-dependent effects (DDEs) represent a major challenge both for calibration of new radio telescopes, and for advancing the state of the art of existing ones. Various approaches to the problem are currently in development. This paper describes several WSRT observations where a pointing error -- commonly thought to be a major contributor to DDEs at higher frequencies -- was deliberately intr…
▽ More
Direction-dependent effects (DDEs) represent a major challenge both for calibration of new radio telescopes, and for advancing the state of the art of existing ones. Various approaches to the problem are currently in development. This paper describes several WSRT observations where a pointing error -- commonly thought to be a major contributor to DDEs at higher frequencies -- was deliberately introduced. These observations were reduced using the MeqTrees package, using two different approaches to correct for DDEs: differential gain solutions, and a direct solution for pointing errors. The results of this are analysed and compared.
△ Less
Submitted 13 October, 2011;
originally announced October 2011.
-
Element gain drifts as an imaging dynamic range limitation in PAF-based interferometers
Authors:
O. M. Smirnov,
M. V. Ivashina
Abstract:
Interferometry with phased-array feeds (PAFs) presents new calibration challenges in comparison with single-pixel feeds. In particular, temporal instability of the compound beam patterns due to element gain drifts (EGDs) can produce calibration artefacts in interferometric images. To translate imaging dynamic range requirements into PAF hardware and calibration requirements, we must learn to relat…
▽ More
Interferometry with phased-array feeds (PAFs) presents new calibration challenges in comparison with single-pixel feeds. In particular, temporal instability of the compound beam patterns due to element gain drifts (EGDs) can produce calibration artefacts in interferometric images. To translate imaging dynamic range requirements into PAF hardware and calibration requirements, we must learn to relate EGD levels to imaging artefact levels. We present a MeqTrees-based simulations framework that addresses this problem, and apply it to the APERTIF prototype currently in development for the WSRT.
△ Less
Submitted 13 October, 2011;
originally announced October 2011.
-
Revisiting the radio interferometer measurement equation. IV. A generalized tensor formalism
Authors:
Oleg M. Smirnov
Abstract:
The radio interferometer measurement equation (RIME), especially in its 2x2 form, has provided a comprehensive matrix-based formalism for describing classical radio interferometry and polarimetry, as shown in the previous three papers of this series. However, recent practical and theoretical developments, such as phased array feeds (PAFs), aperture arrays (AAs) and wide-field polarimetry, are expo…
▽ More
The radio interferometer measurement equation (RIME), especially in its 2x2 form, has provided a comprehensive matrix-based formalism for describing classical radio interferometry and polarimetry, as shown in the previous three papers of this series. However, recent practical and theoretical developments, such as phased array feeds (PAFs), aperture arrays (AAs) and wide-field polarimetry, are exposing limitations of the formalism. This paper aims to develop a more general formalism that can be used to both clearly define the limitations of the matrix RIME, and to describe observational scenarios that lie outside these limitations. Some assumptions underlying the matrix RIME are explicated and analysed in detail. To this purpose, an array correlation matrix (ACM) formalism is explored. This proves of limited use; it is shown that matrix algebra is simply not a sufficiently flexible tool for the job. To overcome these limitations, a more general formalism based on tensors and the Einstein notation is proposed and explored both theoretically, and with a view to practical implementations. The tensor formalism elegantly yields generalized RIMEs describing beamforming, mutual coupling, and wide-field polarimetry in one equation. It is shown that under the explicated assumptions, tensor equations reduce to the 2x2 RIME. From a practical point of view, some methods for implementing tensor equations in an optimal way are proposed and analysed. The tensor RIME is a powerful means of describing observational scenarios not amenable to the matrix RIME. Even in cases where the latter remains applicable, the tensor formalism can be a valuable tool for understanding the limits of such applicability.
△ Less
Submitted 3 June, 2011;
originally announced June 2011.