-
Correcting Bandwidth Depolarization by Extreme Faraday Rotation
Authors:
Maxwell A. Fine,
Cameron L. Van Eck,
Luke Pratley
Abstract:
Measurements of the polarization of radio emission are subject to a number of depolarization effects such as bandwidth depolarization, which is caused by the averaging effect of a finite channel bandwidth combined with the frequency-dependent polarization caused by Faraday rotation. There have been very few mathematical treatments of bandwidth depolarization, especially in the context of the rotat…
▽ More
Measurements of the polarization of radio emission are subject to a number of depolarization effects such as bandwidth depolarization, which is caused by the averaging effect of a finite channel bandwidth combined with the frequency-dependent polarization caused by Faraday rotation. There have been very few mathematical treatments of bandwidth depolarization, especially in the context of the rotation measure (RM) synthesis method for analyzing radio polarization data. We have found a simple equation for predicting if bandwidth depolarization is significant for a given observational configuration. We have derived and tested three methods of modifying RM synthesis to correct for bandwidth depolarization. From these tests we have developed a new algorithm that can detect bandwidth-depolarized signals with higher signal-to-noise than conventional RM synthesis and recover the correct source polarization properties (RM and polarized intensity). We have verified that this algorithm works as expected with real data from the LOFAR Two-metre Sky Survey. To make this algorithm available to the community, we have added it as a new tool in the RM-Tools polarization analysis package.
△ Less
Submitted 6 February, 2023;
originally announced February 2023.
-
The Global Magneto-Ionic Medium Survey (GMIMS): The brightest polarized region in the Southern sky at 75cm and its implications for Radio Loop II
Authors:
Alec J. M. Thomson,
T. L. Landecker,
N. M. McClure-Griffiths,
John M. Dickey,
J. L. Campbell,
Ettore Carretti,
S. E. Clark,
Christoph Federrath,
B. M. Gaensler,
J. L. Han,
Marijke Haverkorn,
Alex. S. Hill,
S. A. Mao,
Anna Ordog,
Luke Pratley,
Wolfgang Reich,
Cameron L. Van Eck,
J. L. West,
M. Wolleben
Abstract:
Using the Global Magneto-Ionic Medium Survey (GMIMS) Low-Band South (LBS) southern sky polarization survey, covering 300 to 480 MHz at 81 arcmin resolution, we reveal the brightest region in the Southern polarized sky at these frequencies. The region, G150-50, covers nearly 20deg$^2$, near (l,b)~(150 deg,-50 deg). Using GMIMS-LBS and complementary data at higher frequencies (~0.6--30 GHz), we appl…
▽ More
Using the Global Magneto-Ionic Medium Survey (GMIMS) Low-Band South (LBS) southern sky polarization survey, covering 300 to 480 MHz at 81 arcmin resolution, we reveal the brightest region in the Southern polarized sky at these frequencies. The region, G150-50, covers nearly 20deg$^2$, near (l,b)~(150 deg,-50 deg). Using GMIMS-LBS and complementary data at higher frequencies (~0.6--30 GHz), we apply Faraday tomography and Stokes QU-fitting techniques. We find that the magnetic field associated with G150-50 is both coherent and primarily in the plane of the sky, and indications that the region is associated with Radio Loop II. The Faraday depth spectra across G150-50 are broad and contain a large-scale spatial gradient. We model the magnetic field in the region as an expanding shell, and we can reproduce both the observed Faraday rotation and the synchrotron emission in the GMIMS-LBS band. Using QU-fitting, we find that the Faraday spectra are produced by several Faraday dispersive sources along the line-of-sight. Alternatively, polarization horizon effects that we cannot model are adding complexity to the high-frequency polarized spectra. The magnetic field structure of Loop II dominates a large fraction of the sky, and studies of the large-scale polarized sky will need to account for this object. Studies of G150-50 with high angular resolution could mitigate polarization horizon effects, and clarify the nature of G150-50.
△ Less
Submitted 23 June, 2021;
originally announced June 2021.
-
Sparse image reconstruction on the sphere: a general approach with uncertainty quantification
Authors:
Matthew A. Price,
Luke Pratley,
Jason D. McEwen
Abstract:
Inverse problems defined naturally on the sphere are becoming increasingly of interest. In this article we provide a general framework for evaluation of inverse problems on the sphere, with a strong emphasis on flexibility and scalability. We consider flexibility with respect to the prior selection (regularization), the problem definition - specifically the problem formulation (constrained/unconst…
▽ More
Inverse problems defined naturally on the sphere are becoming increasingly of interest. In this article we provide a general framework for evaluation of inverse problems on the sphere, with a strong emphasis on flexibility and scalability. We consider flexibility with respect to the prior selection (regularization), the problem definition - specifically the problem formulation (constrained/unconstrained) and problem setting (analysis/synthesis) - and optimization adopted to solve the problem. We discuss and quantify the trade-offs between problem formulation and setting. Crucially, we consider the Bayesian interpretation of the unconstrained problem which, combined with recent developments in probability density theory, permits rapid, statistically principled uncertainty quantification (UQ) in the spherical setting. Linearity is exploited to significantly increase the computational efficiency of such UQ techniques, which in some cases are shown to permit analytic solutions. We showcase this reconstruction framework and UQ techniques on a variety of spherical inverse problems. The code discussed throughout is provided under a GNU general public license, in both C++ and Python.
△ Less
Submitted 11 May, 2021;
originally announced May 2021.
-
Removing non-physical structure in fitted Faraday rotated signals: non-parametric QU-fitting
Authors:
Luke Pratley,
Melanie Johnston-Hollitt,
Bryan M. Gaensler
Abstract:
Next-generation spectro-polarimetric broadband surveys will probe cosmic magnetic fields in unprecedented detail, using the magneto-optical effect known as Faraday rotation. However, non-parametric methods such as RMCLEAN can introduce non-observable linearly polarized flux into a fitted model at negative wavelengths squared. This leads to Faraday rotation structures that are consistent with the o…
▽ More
Next-generation spectro-polarimetric broadband surveys will probe cosmic magnetic fields in unprecedented detail, using the magneto-optical effect known as Faraday rotation. However, non-parametric methods such as RMCLEAN can introduce non-observable linearly polarized flux into a fitted model at negative wavelengths squared. This leads to Faraday rotation structures that are consistent with the observed data, but would be impossible or difficult to measure. We construct a convex non-parametric $QU$-fitting algorithm to constrain the flux at negative wavelengths squared to be zero. This allows the algorithm to recover structures that are limited in complexity to the observable region in wavelength squared. We verify this approach on simulated broadband data sets where we show that it has a lower root mean square error and that it can change the scientific conclusions for real observations. We advise using this prior in next-generation broadband surveys that aim to uncover complex Faraday depth structures. We provide a public Python implementation of the algorithm at \url{https://github.com/Luke-Pratley/Faraday-Dreams}.
△ Less
Submitted 15 October, 2021; v1 submitted 15 October, 2020;
originally announced October 2020.
-
Novel perspectives gained from new reconstruction algorithms
Authors:
Luke Pratley,
Melanie Johnston-Hollitt,
Jason D. McEwen
Abstract:
Since the 1970s, much of traditional interferometric imaging has been built around variations of the CLEAN algorithm, in both terminology, methodology, and algorithm development. Recent developments in applying new algorithms from convex optimization to interferometry has allowed old concepts to be viewed from a new perspective, ranging from image restoration to the development of computationally…
▽ More
Since the 1970s, much of traditional interferometric imaging has been built around variations of the CLEAN algorithm, in both terminology, methodology, and algorithm development. Recent developments in applying new algorithms from convex optimization to interferometry has allowed old concepts to be viewed from a new perspective, ranging from image restoration to the development of computationally distributed algorithms. We present how this has ultimately led the authors to new perspectives in wide-field imaging, allowing for the first full individual non-coplanar corrections applied during imaging over extremely wide-fields of view for the Murchison Widefield Array (MWA) telescope. Furthermore, this same mathematical framework has provided a novel understanding of wide-band polarimetry at low frequencies, where instrumental channel depolarization can be corrected through the new $δλ^2$-projection algorithm. This is a demonstration that new algorithm development outside of traditional radio astronomy is valuable for the new theoretical and practical perspectives gained. These perspectives are timely with the next generation of radio telescopes coming online.
△ Less
Submitted 26 September, 2020;
originally announced September 2020.
-
Magnetism Science with the Square Kilometre Array
Authors:
George Heald,
Sui Ann Mao,
Valentina Vacca,
Takuya Akahori,
Ancor Damas-Segovia,
B. M. Gaensler,
Matthias Hoeft,
Ivan Agudo,
Aritra Basu,
Rainer Beck,
Mark Birkinshaw,
Annalisa Bonafede,
Tyler L. Bourke,
Andrea Bracco,
Ettore Carretti,
Luigina Feretti,
J. M. Girart,
Federica Govoni,
James A. Green,
JinLin Han,
Marijke Haverkorn,
Cathy Horellou,
Melanie Johnston-Hollitt,
Roland Kothes,
Tom Landecker
, et al. (19 additional authors not shown)
Abstract:
The Square Kilometre Array (SKA) will answer fundamental questions about the origin, evolution, properties, and influence of magnetic fields throughout the Universe. Magnetic fields can illuminate and influence phenomena as diverse as star formation, galactic dynamics, fast radio bursts, active galactic nuclei, large-scale structure, and Dark Matter annihilation. Preparations for the SKA are swift…
▽ More
The Square Kilometre Array (SKA) will answer fundamental questions about the origin, evolution, properties, and influence of magnetic fields throughout the Universe. Magnetic fields can illuminate and influence phenomena as diverse as star formation, galactic dynamics, fast radio bursts, active galactic nuclei, large-scale structure, and Dark Matter annihilation. Preparations for the SKA are swiftly continuing worldwide, and the community is making tremendous observational progress in the field of cosmic magnetism using data from a powerful international suite of SKA pathfinder and precursor telescopes. In this contribution, we revisit community plans for magnetism research using the SKA, in the light of these recent rapid developments. We focus in particular on the impact that new radio telescope instrumentation is generating, thus advancing our understanding of key SKA magnetism science areas, as well as the new techniques that are required for processing and interpreting the data. We discuss these recent developments in the context of the ultimate scientific goals for the SKA era.
△ Less
Submitted 4 June, 2020;
originally announced June 2020.
-
Spherical Bayesian mass-mapping with uncertainties: full sky observations on the celestial sphere
Authors:
Matthew A. Price,
Jason D. McEwen,
L. Pratley,
Thomas D. Kitching
Abstract:
To date weak gravitational lensing surveys have typically been restricted to small fields of view, such that the $\textit{flat-sky approximation}$ has been sufficiently satisfied. However, with Stage IV surveys ($\textit{e.g. LSST}$ and $\textit{Euclid}$) imminent, extending mass-mapping techniques to the sphere is a fundamental necessity. As such, we extend the sparse hierarchical Bayesian mass-m…
▽ More
To date weak gravitational lensing surveys have typically been restricted to small fields of view, such that the $\textit{flat-sky approximation}$ has been sufficiently satisfied. However, with Stage IV surveys ($\textit{e.g. LSST}$ and $\textit{Euclid}$) imminent, extending mass-mapping techniques to the sphere is a fundamental necessity. As such, we extend the sparse hierarchical Bayesian mass-mapping formalism presented in previous work to the spherical sky. For the first time, this allows us to construct $\textit{maximum a posteriori}$ spherical weak lensing dark-matter mass-maps, with principled Bayesian uncertainties, without imposing or assuming Gaussianty. We solve the spherical mass-mapping inverse problem in the analysis setting adopting a sparsity promoting Laplace-type wavelet prior, though this theoretical framework supports all log-concave posteriors. Our spherical mass-mapping formalism facilitates principled statistical interpretation of reconstructions. We apply our framework to convergence reconstruction on high resolution N-body simulations with pseudo-Euclid masking, polluted with a variety of realistic noise levels, and show a significant increase in reconstruction fidelity compared to standard approaches. Furthermore we perform the largest joint reconstruction to date of the majority of publicly available shear observational datasets (combining DESY1, KiDS450 and CFHTLens) and find that our formalism recovers a convergence map with significantly enhanced small-scale detail. Within our Bayesian framework we validate, in a statistically rigorous manner, the community's intuition regarding the need to smooth spherical Kaiser-Squires estimates to provide physically meaningful convergence maps. Such approaches cannot reveal the small-scale physical structures that we recover within our framework.
△ Less
Submitted 5 February, 2021; v1 submitted 16 April, 2020;
originally announced April 2020.
-
Offline and online reconstruction for radio interferometric imaging
Authors:
Xiaohao Cai,
Luke Pratley,
Jason D. McEwen
Abstract:
Radio astronomy is transitioning to a big-data era due to the emerging generation of radio interferometric (RI) telescopes, such as the Square Kilometre Array (SKA), which will acquire massive volumes of data. In this article we review methods proposed recently to resolve the ill-posed inverse problem of imaging the raw visibilities acquired by RI telescopes in the big-data scenario. We focus on t…
▽ More
Radio astronomy is transitioning to a big-data era due to the emerging generation of radio interferometric (RI) telescopes, such as the Square Kilometre Array (SKA), which will acquire massive volumes of data. In this article we review methods proposed recently to resolve the ill-posed inverse problem of imaging the raw visibilities acquired by RI telescopes in the big-data scenario. We focus on the recently proposed online reconstruction method [4] and the considerable savings in data storage requirements and computational cost that it yields.
△ Less
Submitted 8 April, 2020;
originally announced April 2020.
-
Wide-band Rotation Measure Synthesis
Authors:
Luke Pratley,
Melanie Johnston-Hollitt
Abstract:
Rotation measure synthesis allows the estimation of Faraday dispersion via a Fourier transform and is the primary tool to probe cosmic magnetic fields. We show this can be considered mathematically equivalent to the one dimensional interferometric intensity measurement equation, albeit in a different Fourier space. As a result, familiar concepts in two dimensional intensity interferometry designed…
▽ More
Rotation measure synthesis allows the estimation of Faraday dispersion via a Fourier transform and is the primary tool to probe cosmic magnetic fields. We show this can be considered mathematically equivalent to the one dimensional interferometric intensity measurement equation, albeit in a different Fourier space. As a result, familiar concepts in two dimensional intensity interferometry designed to correctly account for a range of instrumental conditions can be translated to the analysis of Faraday dispersion. In particular, we show how to model the effect of channel averaging during Faraday reconstruction, which has to date limited the progress of polarimetic science using wide-band measurements. Further, we simulate 1d sparse reconstruction with channel averaging for realistic frequency coverages, and show that it is possible to recover signals with large rotation measure values that were previously excluded from possible detection. This is especially important for low-frequency and wide-band polarimetry. We extended these ideas to introduce mosaicking in Faraday depth into the channel averaging process. This work, thus provides the first framework for correctly undertaking wide-band rotation measure synthesis, including the provision to add data from multiple telescopes, a prospect that should vastly improve the quality and quantity of polarimetric science. This is of particular importance for extreme environments which generate high magnetic fields such as those associated with pulsars and Fast Radio Bursts (FRBs), and will allow such sources to be accurately used as probes of cosmological fields.
△ Less
Submitted 21 January, 2020; v1 submitted 3 June, 2019;
originally announced June 2019.
-
Load balancing for distributed interferometric image reconstruction
Authors:
Luke Pratley,
Jason D. McEwen
Abstract:
We present a new algorithm to perform wide-field radio interferometric image reconstruction, with exact non-coplanar correction, that scales to big-data. This algorithm allows us to image 2 billion visibilities on 50 nodes of a computing cluster for a 25 by 25 degree field of view, in a little over an hour. We build on the recently developed distributed $w$-stacking $w$-projection hybrid algorithm…
▽ More
We present a new algorithm to perform wide-field radio interferometric image reconstruction, with exact non-coplanar correction, that scales to big-data. This algorithm allows us to image 2 billion visibilities on 50 nodes of a computing cluster for a 25 by 25 degree field of view, in a little over an hour. We build on the recently developed distributed $w$-stacking $w$-projection hybrid algorithm, extending it to include a new distributed degridding algorithm that balances the computational load of the $w$-projection gridding kernels. The implementation of our algorithm is made publicly available in the PURIFY software package. Wide-field image reconstruction for data sets of this size cannot be performed effectively using the allocated computational resources without computational load balancing, demonstrating that our algorithms are critical for next-generation wide-field radio interferometers.
△ Less
Submitted 18 March, 2019;
originally announced March 2019.
-
$w$-stacking $w$-projection hybrid algorithm for wide-field interferometric imaging: implementation details and improvements
Authors:
L. Pratley,
M. Johnston-Hollitt,
J. D. McEwen
Abstract:
We present a detailed discussion of the implementation strategies for a recently developed $w$-stacking $w$-projection hybrid algorithm used to reconstruct wide-field interferometric images. In particular, we discuss the methodology used to deploy the algorithm efficiently on a supercomputer via use of a Message Passing Interface (MPI) $k$-means clustering technique to achieve efficient constructi…
▽ More
We present a detailed discussion of the implementation strategies for a recently developed $w$-stacking $w$-projection hybrid algorithm used to reconstruct wide-field interferometric images. In particular, we discuss the methodology used to deploy the algorithm efficiently on a supercomputer via use of a Message Passing Interface (MPI) $k$-means clustering technique to achieve efficient construction and application of non co-planar effects. Additionally, we show that the use of conjugate symmetry increases the algorithms performance by imaging an interferometric observation of Fornax A from the Murchison Widefield Array (MWA). We perform exact non-coplanar wide-field correction for 126.6 million visibilities using 50 nodes of a computing cluster. The $w$-projection kernel construction takes only 15 minutes, demonstrating that the implementation is both fast and efficient.
△ Less
Submitted 13 March, 2019;
originally announced March 2019.
-
Sparse Image Reconstruction for the SPIDER Optical Interferometric Telescope
Authors:
Luke Pratley,
Jason D. McEwen
Abstract:
The concept of a recently proposed small-scale interferometric optical imaging device, an instrument known as the Segmented Planar Imaging Detector for Electro-optical Reconnaissance (SPIDER), is of great interest for its possible applications in astronomy and space science. Due to low weight, low power consumption, and high resolution, the SPIDER telescope could replace the large space telescopes…
▽ More
The concept of a recently proposed small-scale interferometric optical imaging device, an instrument known as the Segmented Planar Imaging Detector for Electro-optical Reconnaissance (SPIDER), is of great interest for its possible applications in astronomy and space science. Due to low weight, low power consumption, and high resolution, the SPIDER telescope could replace the large space telescopes that exist today. Unlike traditional optical interferometry the SPIDER accurately retrieves both phase and amplitude information, making the measurement process analogous to a radio interferometer. State of the art sparse radio interferometric image reconstruction techniques have been gaining traction in radio astronomy and reconstruct accurate images of the radio sky. In this work we describe algorithms from radio interferometric imaging and sparse image reconstruction and demonstrate their application to the SPIDER concept telescope through simulated observation and reconstruction of the optical sky. Such algorithms are important for providing high fidelity images from SPIDER observations, helping to power the SPIDER concept for scientific and astronomical analysis.
△ Less
Submitted 22 March, 2021; v1 submitted 13 March, 2019;
originally announced March 2019.
-
Distributed and parallel sparse convex optimization for radio interferometry with PURIFY
Authors:
Luke Pratley,
Jason D. McEwen,
Mayeul d'Avezac,
Xiaohao Cai,
David Perez-Suarez,
Ilektra Christidi,
Roland Guichard
Abstract:
Next generation radio interferometric telescopes are entering an era of big data with extremely large data sets. While these telescopes can observe the sky in higher sensitivity and resolution than before, computational challenges in image reconstruction need to be overcome to realize the potential of forthcoming telescopes. New methods in sparse image reconstruction and convex optimization techni…
▽ More
Next generation radio interferometric telescopes are entering an era of big data with extremely large data sets. While these telescopes can observe the sky in higher sensitivity and resolution than before, computational challenges in image reconstruction need to be overcome to realize the potential of forthcoming telescopes. New methods in sparse image reconstruction and convex optimization techniques (cf. compressive sensing) have shown to produce higher fidelity reconstructions of simulations and real observations than traditional methods. This article presents distributed and parallel algorithms and implementations to perform sparse image reconstruction, with significant practical considerations that are important for implementing these algorithms for Big Data. We benchmark the algorithms presented, showing that they are considerably faster than their serial equivalents. We then pre-sample gridding kernels to scale the distributed algorithms to larger data sizes, showing application times for 1 Gb to 2.4 Tb data sets over 25 to 100 nodes for up to 50 billion visibilities, and find that the run-times for the distributed algorithms range from 100 milliseconds to 3 minutes per iteration. This work presents an important step in working towards computationally scalable and efficient algorithms and implementations that are needed to image observations of both extended and compact sources from next generation radio interferometers such as the SKA. The algorithms are implemented in the latest versions of the SOPT (https://github.com/astro-informatics/sopt) and PURIFY (https://github.com/astro-informatics/purify) software packages {(Versions 3.1.0)}, which have been released alongside of this article.
△ Less
Submitted 13 December, 2019; v1 submitted 11 March, 2019;
originally announced March 2019.
-
A fast and exact $w$-stacking and $w$-projection hybrid algorithm for wide-field interferometric imaging
Authors:
Luke Pratley,
Melanie Johnston-Hollitt,
Jason D. McEwen
Abstract:
The standard wide-field imaging technique, the $w$-projection, allows correction for wide-fields of view for non-coplanar radio interferometric arrays. However, calculating exact corrections for each measurement has not been possible due to the amount of computation required at high resolution and with the large number of visibilities from current interferometers. The required accuracy and computa…
▽ More
The standard wide-field imaging technique, the $w$-projection, allows correction for wide-fields of view for non-coplanar radio interferometric arrays. However, calculating exact corrections for each measurement has not been possible due to the amount of computation required at high resolution and with the large number of visibilities from current interferometers. The required accuracy and computational cost of these corrections is one of the largest unsolved challenges facing next generation radio interferometers such as the Square Kilometre Array. We show that the same calculation can be performed with a radially symmetric $w$-projection kernel, where we use one dimensional adaptive quadrature to calculate the resulting Hankel transform, decreasing the computation required for kernel generation by several orders of magnitude, whilst preserving the accuracy. We confirm that the radial $w$-projection kernel is accurate to approximately 1% by imaging the zero-spacing with an added $w$-term. We demonstrate the potential of our radially symmetric $w$-projection kernel via sparse image reconstruction, using the software package PURIFY. We develop a distributed $w$-stacking and $w$-projection hybrid algorithm. We apply this algorithm to individually correct for non-coplanar effects in 17.5 million visibilities over a $25$ by $25$ degree field of view MWA observation for image reconstruction. Such a level of accuracy and scalability is not possible with standard $w$-projection kernel generation methods. This demonstrates that we can scale to a large number of measurements with large image sizes whilst still maintaining both speed and accuracy.
△ Less
Submitted 14 February, 2019; v1 submitted 24 July, 2018;
originally announced July 2018.
-
Online radio interferometric imaging: assimilating and discarding visibilities on arrival
Authors:
Xiaohao Cai,
Luke Pratley,
Jason D. McEwen
Abstract:
The emerging generation of radio interferometric (RI) telescopes, such as the Square Kilometre Array (SKA), will acquire massive volumes of data and transition radio astronomy to a big-data era. The ill-posed inverse problem of imaging the raw visibilities acquired by RI telescopes will become significantly more computationally challenging, particularly in terms of data storage and computational c…
▽ More
The emerging generation of radio interferometric (RI) telescopes, such as the Square Kilometre Array (SKA), will acquire massive volumes of data and transition radio astronomy to a big-data era. The ill-posed inverse problem of imaging the raw visibilities acquired by RI telescopes will become significantly more computationally challenging, particularly in terms of data storage and computational cost. Current RI imaging methods, such as CLEAN, its variants, and compressive sensing approaches (sparse regularisation), have yielded excellent reconstruction fidelity. However, scaling these methods to big-data remains difficult if not impossible in some cases. All state-of-the-art methods in RI imaging lack the ability to process data streams as they are acquired during the data observation stage. Such approaches are referred to as online processing methods. We present an online sparse regularisation methodology for RI imaging. Image reconstruction is performed simultaneously with data acquisition, where observed visibilities are assimilated into the reconstructed image as they arrive and then discarded. Since visibilities are processed online, good reconstructions are recovered much faster than standard (offline) methods which cannot start until the data acquisition stage completes. Moreover, the online method provides additional computational savings and, most importantly, dramatically reduces data storage requirements. Theoretically, the reconstructed images are of the same fidelity as those recovered by the equivalent offline approach and, in practice, very similar reconstruction fidelity is achieved. We anticipate online imaging techniques, as proposed here, will be critical in scaling RI imaging to the emerging big-data era of radio astronomy.
△ Less
Submitted 12 December, 2017;
originally announced December 2017.
-
Upper limits on a radio halo in Abell 3667 at 1.4 GHz
Authors:
M. Johnston-Hollitt,
L. Pratley
Abstract:
The presence of a radio halo in the massive, merging cluster Abell 3667 has recently become a focus of debate in the literature following a putative halo detection at 2.4 GHz despite a lack of detection at a range of lower frequencies between 120 MHz and 1.8 GHz. Here we develop a new method to place limits on radio haloes via generation of a realistic synthetic halo based on the brightness distri…
▽ More
The presence of a radio halo in the massive, merging cluster Abell 3667 has recently become a focus of debate in the literature following a putative halo detection at 2.4 GHz despite a lack of detection at a range of lower frequencies between 120 MHz and 1.8 GHz. Here we develop a new method to place limits on radio haloes via generation of a realistic synthetic halo based on the brightness distribution of real haloes. The model generated extends on previous methods in the literature producing a single elliptical halo model, capable of being injected into mosaic as well as single observations. Applying this model to the deepest data available 1.4 GHz data for A3667 we derive an upper limit halo power of P$_{1.4} \leq 5.55 \times 10^{23}$ W Hz$^{-1}$. We discuss the result in the context of current scaling relation between the X-ray and radio properties of galaxy clusters and find that the lack of a halo in A3667 places the cluster on the border of the so-called `off-state' region in which clusters are expected not to host any diffuse emission.
△ Less
Submitted 15 June, 2017;
originally announced June 2017.
-
LOFAR MSSS: The Scaling Relation between AGN Cavity Power and Radio Luminosity at Low Radio Frequencies
Authors:
G. Kokotanekov,
M. Wise,
G. H. Heald,
J. P. McKean,
L. Bîrzan,
D. A. Rafferty,
L. E. H. Godfrey,
M. de Vries,
H. T. Intema,
J. W. Broderick,
M. J. Hardcastle,
A. Bonafede,
A. O. Clarke,
R. J. van Weeren,
H. J. A. Röttgering,
R. Pizzo,
M. Iacobelli,
E. Orrú,
A. Shulevski,
C. J. Riseley,
R. P. Breton,
B. Nikiel-Wroczyński,
S. S. Sridhar,
A. J. Stewart,
A. Rowlinson
, et al. (13 additional authors not shown)
Abstract:
We present a new analysis of the widely used relation between cavity power and radio luminosity in clusters of galaxies with evidence for strong AGN feedback. We study the correlation at low radio frequencies using two new surveys - the First Alternative Data Release of the TIFR GMRT Sky Survey (TGSS ADR1) at 148 MHz and LOFAR's first all-sky survey, the Multifrequency Snapshot Sky Survey (MSSS) a…
▽ More
We present a new analysis of the widely used relation between cavity power and radio luminosity in clusters of galaxies with evidence for strong AGN feedback. We study the correlation at low radio frequencies using two new surveys - the First Alternative Data Release of the TIFR GMRT Sky Survey (TGSS ADR1) at 148 MHz and LOFAR's first all-sky survey, the Multifrequency Snapshot Sky Survey (MSSS) at 140 MHz. We find a scaling relation $P_{\rm cav} \propto L_{148}^β$, with a logarithmic slope of $β= 0.51 \pm 0.14$, which is in good agreement with previous results based on data at 327 MHz. The large scatter present in this correlation confirms the conclusion reached at higher frequencies that the total radio luminosity at a single frequency is a poor predictor of the total jet power. We show that including measurements at 148 MHz alone is insufficient to reliably compute the bolometric radio luminosity and reduce the scatter in the correlation. For a subset of four well-resolved sources, we examine the detected extended structures at low frequencies and compare with the morphology known from higher frequency images and Chandra X-ray maps. In Perseus we discuss details in the structures of the radio mini-halo, while in the 2A 0335+096 cluster we observe new diffuse emission associated with multiple X-ray cavities and likely originating from past activity. For A2199 and MS 0735.6+7421, we confirm that the observed low-frequency radio lobes are confined to the extents known from higher frequencies. This new low-frequency analysis highlights the fact that existing cavity power to radio luminosity relations are based on a relatively narrow range of AGN outburst ages. We discuss how the correlation could be extended using low frequency data from the LOFAR Two-metre Sky Survey (LoTSS) in combination with future, complementary deeper X-ray observations.
△ Less
Submitted 1 June, 2017;
originally announced June 2017.
-
PURIFYing real radio interferometric observations
Authors:
Luke Pratley,
Jason D. McEwen,
Mayeul d'Avezac,
Rafael E. Carrillo,
Alexandru Onose,
Yves Wiaux
Abstract:
Next-generation radio interferometers, such as the Square Kilometre Array (SKA), will revolutionise our understanding of the universe through their unprecedented sensitivity and resolution. However, standard methods in radio interferometry produce reconstructed interferometric images that are limited in quality and they are not scalable for big data. In this work we apply and evaluate alternative…
▽ More
Next-generation radio interferometers, such as the Square Kilometre Array (SKA), will revolutionise our understanding of the universe through their unprecedented sensitivity and resolution. However, standard methods in radio interferometry produce reconstructed interferometric images that are limited in quality and they are not scalable for big data. In this work we apply and evaluate alternative interferometric reconstruction methods that make use of state-of-the-art sparse image reconstruction algorithms motivated by compressive sensing, which have been implemented in the PURIFY software package. In particular, we implement and apply the proximal alternating direction method of multipliers (P-ADMM) algorithm presented in a recent article. We apply PURIFY to real interferometric observations. For all observations PURIFY outperforms the standard CLEAN, where in some cases PURIFY provides an improvement in dynamic range by over an order of magnitude. The latest version of PURIFY, which includes the developments presented in this work, is made publicly available.
△ Less
Submitted 22 February, 2017;
originally announced February 2017.
-
The w-effect in interferometric imaging: from a fast sparse measurement operator to super-resolution
Authors:
Arwa Dabbech,
Laura Wolz,
Luke Pratley,
Jason D. McEwen,
Yves Wiaux
Abstract:
Modern radio telescopes, such as the Square Kilometre Array (SKA), will probe the radio sky over large fields-of-view, which results in large w-modulations of the sky image. This effect complicates the relationship between the measured visibilities and the image under scrutiny. In algorithmic terms, it gives rise to massive memory and computational time requirements. Yet, it can be a blessing in t…
▽ More
Modern radio telescopes, such as the Square Kilometre Array (SKA), will probe the radio sky over large fields-of-view, which results in large w-modulations of the sky image. This effect complicates the relationship between the measured visibilities and the image under scrutiny. In algorithmic terms, it gives rise to massive memory and computational time requirements. Yet, it can be a blessing in terms of reconstruction quality of the sky image. In recent years, several works have shown that large w-modulations promote the spread spectrum effect. Within the compressive sensing framework, this effect increases the incoherence between the sensing basis and the sparsity basis of the signal to be recovered, leading to better estimation of the sky image. In this article, we revisit the w-projection approach using convex optimisation in realistic settings, where the measurement operator couples the w-terms in Fourier and the de-gridding kernels. We provide sparse, thus fast, models of the Fourier part of the measurement operator through adaptive sparsification procedures. Consequently, memory requirements and computational cost are significantly alleviated, at the expense of introducing errors on the radio-interferometric data model. We present a first investigation of the impact of the sparse variants of the measurement operator on the image reconstruction quality. We finally analyse the interesting super-resolution potential associated with the spread spectrum effect of the w-modulation, and showcase it through simulations. Our C++ code is available online on GitHub.
△ Less
Submitted 13 July, 2017; v1 submitted 16 February, 2017;
originally announced February 2017.
-
LOFAR MSSS: Discovery of a 2.56 Mpc giant radio galaxy associated with a disturbed galaxy group
Authors:
A. O. Clarke,
G. Heald,
T. Jarrett,
J. D. Bray,
M. J. Hardcastle,
T. M. Cantwell,
A. M. M. Scaife,
M. Brienza,
A. Bonafede,
R. P. Breton,
J. W. Broderick,
D. Carbone,
J. H. Croston,
J. S. Farnes,
J. J. Harwood,
V. Heesen,
A. Horneffer,
A. J. van der Horst,
M. Iacobelli,
W. Jurusik,
G. Kokotanekov,
J. P. McKean,
L. K. Morabito,
D. D. Mulcahy,
B. S. Nikiel-Wroczynski
, et al. (18 additional authors not shown)
Abstract:
We report on the discovery in the LOFAR Multifrequency Snapshot Sky Survey (MSSS) of a giant radio galaxy (GRG) with a projected size of $2.56 \pm 0.07$ Mpc projected on the sky. It is associated with the galaxy triplet UGC 9555, within which one is identified as a broad-line galaxy in the Sloan Digital Sky Survey (SDSS) at a redshift of $0.05453 \pm 1 \times 10^{-5} $, and with a velocity dispers…
▽ More
We report on the discovery in the LOFAR Multifrequency Snapshot Sky Survey (MSSS) of a giant radio galaxy (GRG) with a projected size of $2.56 \pm 0.07$ Mpc projected on the sky. It is associated with the galaxy triplet UGC 9555, within which one is identified as a broad-line galaxy in the Sloan Digital Sky Survey (SDSS) at a redshift of $0.05453 \pm 1 \times 10^{-5} $, and with a velocity dispersion of $215.86 \pm 6.34$ km/s. From archival radio observations we see that this galaxy hosts a compact flat-spectrum radio source, and we conclude that it is the active galactic nucleus (AGN) responsible for generating the radio lobes. The radio luminosity distribution of the jets, and the broad-line classification of the host AGN, indicate this GRG is orientated well out of the plane of the sky, making its physical size one of the largest known for any GRG. Analysis of the infrared data suggests that the host is a lenticular type galaxy with a large stellar mass ($\log~\mathrm{M}/\mathrm{M}_\odot = 11.56 \pm 0.12$), and a moderate star formation rate ($1.2 \pm 0.3~\mathrm{M}_\odot/\mathrm{year}$). Spatially smoothing the SDSS images shows the system around UGC 9555 to be significantly disturbed, with a prominent extension to the south-east. Overall, the evidence suggests this host galaxy has undergone one or more recent moderate merger events and is also experiencing tidal interactions with surrounding galaxies, which have caused the star formation and provided the supply of gas to trigger and fuel the Mpc-scale radio lobes.
△ Less
Submitted 6 February, 2017;
originally announced February 2017.
-
Robust sparse image reconstruction of radio interferometric observations with purify
Authors:
Luke Pratley,
Jason D. McEwen,
Mayeul d'Avezac,
Rafael E. Carrillo,
Alexandru Onose,
Yves Wiaux
Abstract:
Next-generation radio interferometers, such as the Square Kilometre Array (SKA), will revolutionise our understanding of the universe through their unprecedented sensitivity and resolution. However, to realise these goals significant challenges in image and data processing need to be overcome. The standard methods in radio interferometry for reconstructing images, such as CLEAN, have served the co…
▽ More
Next-generation radio interferometers, such as the Square Kilometre Array (SKA), will revolutionise our understanding of the universe through their unprecedented sensitivity and resolution. However, to realise these goals significant challenges in image and data processing need to be overcome. The standard methods in radio interferometry for reconstructing images, such as CLEAN, have served the community well over the last few decades and have survived largely because they are pragmatic. However, they produce reconstructed inter\-ferometric images that are limited in quality and scalability for big data. In this work we apply and evaluate alternative interferometric reconstruction methods that make use of state-of-the-art sparse image reconstruction algorithms motivated by compressive sensing, which have been implemented in the PURIFY software package. In particular, we implement and apply the proximal alternating direction method of multipliers (P-ADMM) algorithm presented in a recent article. First, we assess the impact of the interpolation kernel used to perform gridding and degridding on sparse image reconstruction. We find that the Kaiser-Bessel interpolation kernel performs as well as prolate spheroidal wave functions, while providing a computational saving and an analytic form. Second, we apply PURIFY to real interferometric observations from the Very Large Array (VLA) and the Australia Telescope Compact Array (ATCA) and find images recovered by PURIFY are higher quality than those recovered by CLEAN. Third, we discuss how PURIFY reconstructions exhibit additional advantages over those recovered by CLEAN. The latest version of PURIFY, with developments presented in this work, is made publicly available.
△ Less
Submitted 16 November, 2017; v1 submitted 7 October, 2016;
originally announced October 2016.
-
ATCA observations of the MACS-Planck Radio Halo Cluster Project - I. New detection of a radio halo in PLCK G285.0-23.7
Authors:
Gerardo Martinez Aviles,
Chiara Ferrari,
Melanie Johnston-Hollitt,
Luke Pratley,
Giulia Macario,
Tiziana Venturi,
Gianfranco Brunetti,
Rossella Cassano,
Daniele Dallacasa,
Huib Intema,
Simona Giacintucci,
Guillaume Hurier,
Nabila Aghanim,
Marian Douspis,
Mathieu Langer
Abstract:
We investigate the possible presence of diffuse radio emission in the intermediate redshift, massive cluster PLCK G285.0-23.7 (z=0.39, M_500 = 8.39 x 10^(14) M_Sun). Our 16cm-band ATCA observations of PLCK G285.0-23.7 allow us to reach a rms noise level of ~11 microJy/beam on the wide-band (1.1-3.1 GHz), full-resolution (~5 arcsec) image of the cluster, making it one of the deepest ATCA images yet…
▽ More
We investigate the possible presence of diffuse radio emission in the intermediate redshift, massive cluster PLCK G285.0-23.7 (z=0.39, M_500 = 8.39 x 10^(14) M_Sun). Our 16cm-band ATCA observations of PLCK G285.0-23.7 allow us to reach a rms noise level of ~11 microJy/beam on the wide-band (1.1-3.1 GHz), full-resolution (~5 arcsec) image of the cluster, making it one of the deepest ATCA images yet published. We also re-image visibilities at lower resolution in order to achieve a better sensitivity to low-surface-brightness extended radio sources. We detect one of the lowest luminosity radio halos known at z>0.35, characterised by a slight offset from the well-studied 1.4 GHz radio power vs. cluster mass correlation. Similarly to most known radio-loud clusters (i.e. those hosting diffuse non-thermal sources), PLCK G285.0-23.7 has a disturbed dynamical state. Our analysis reveals a similarly elongated X-ray and radio morphology. While the size of the radio halo in PLCK G285.0-23.7 is smaller than lower redshift radio-loud clusters in the same mass range, it shows a similar correlation with the cluster virial radius, as expected in the framework of hierarchical structure formation.
△ Less
Submitted 6 September, 2016; v1 submitted 24 August, 2016;
originally announced August 2016.
-
An improved method for polarimetric image restoration in interferometry
Authors:
Luke Pratley,
Melanie Johnston-Hollitt
Abstract:
Interferometric radio astronomy data require the effects of limited coverage in the Fourier plane to be accounted for via a deconvolution process. For the last 40 years this process, known as `cleaning', has been performed almost exclusively on all Stokes parameters individually as if they were independent scalar images. However, here we demonstrate for the case of the linear polarisation…
▽ More
Interferometric radio astronomy data require the effects of limited coverage in the Fourier plane to be accounted for via a deconvolution process. For the last 40 years this process, known as `cleaning', has been performed almost exclusively on all Stokes parameters individually as if they were independent scalar images. However, here we demonstrate for the case of the linear polarisation $\mathcal{P}$, this approach fails to properly account for the complex vector nature resulting in a process which is dependant on the axis under which the deconvolution is performed. We present here an improved method, `Generalised Complex CLEAN', which properly accounts for the complex vector nature of polarised emission and is invariant under rotations of the deconvolution axis. We use two Australia Telescope Compact Array datasets to test standard and complex CLEAN versions of the Högbom and SDI CLEAN algorithms. We show that in general the Complex CLEAN version of each algorithm produces more accurate clean components with fewer spurious detections and lower computation cost due to reduced iterations than the current methods. In particular we find that the Complex SDI CLEAN produces the best results for diffuse polarised sources as compared with standard CLEAN algorithms and other Complex CLEAN algorithms. Given the move to widefield, high resolution polarimetric imaging with future telescopes such as the Square Kilometre Array, we suggest that Generalised Complex CLEAN should be adopted as the deconvolution method for all future polarimetric surveys and in particular that the complex version of a SDI CLEAN should be used.
△ Less
Submitted 5 June, 2016;
originally announced June 2016.
-
LOFAR MSSS: Detection of a low-frequency radio transient in 400 hrs of monitoring of the North Celestial Pole
Authors:
A. J. Stewart,
R. P. Fender,
J. W. Broderick,
T. E. Hassall,
T. Muñoz-Darias,
A. Rowlinson,
J. D. Swinbank,
T. D. Staley,
G. J. Molenaar,
B. Scheers,
T. L. Grobler,
M. Pietka,
G. Heald,
J. P. McKean,
M. E. Bell,
A. Bonafede,
R. P. Breton,
D. Carbone,
Y. Cendes,
A. O. Clarke,
S. Corbel,
F. de Gasperin,
J. Eislöffel,
H. Falcke,
C. Ferrari
, et al. (77 additional authors not shown)
Abstract:
We present the results of a four-month campaign searching for low-frequency radio transients near the North Celestial Pole with the Low-Frequency Array (LOFAR), as part of the Multifrequency Snapshot Sky Survey (MSSS). The data were recorded between 2011 December and 2012 April and comprised 2149 11-minute snapshots, each covering 175 deg^2. We have found one convincing candidate astrophysical tra…
▽ More
We present the results of a four-month campaign searching for low-frequency radio transients near the North Celestial Pole with the Low-Frequency Array (LOFAR), as part of the Multifrequency Snapshot Sky Survey (MSSS). The data were recorded between 2011 December and 2012 April and comprised 2149 11-minute snapshots, each covering 175 deg^2. We have found one convincing candidate astrophysical transient, with a duration of a few minutes and a flux density at 60 MHz of 15-25 Jy. The transient does not repeat and has no obvious optical or high-energy counterpart, as a result of which its nature is unclear. The detection of this event implies a transient rate at 60 MHz of 3.9 (+14.7, -3.7) x 10^-4 day^-1 deg^-2, and a transient surface density of 1.5 x 10^-5 deg^-2, at a 7.9-Jy limiting flux density and ~10-minute time-scale. The campaign data were also searched for transients at a range of other time-scales, from 0.5 to 297 min, which allowed us to place a range of limits on transient rates at 60 MHz as a function of observation duration.
△ Less
Submitted 30 November, 2015;
originally announced December 2015.
-
The LOFAR Multifrequency Snapshot Sky Survey (MSSS) I. Survey description and first results
Authors:
G. H. Heald,
R. F. Pizzo,
E. Orrú,
R. P. Breton,
D. Carbone,
C. Ferrari,
M. J. Hardcastle,
W. Jurusik,
G. Macario,
D. Mulcahy,
D. Rafferty,
A. Asgekar,
M. Brentjens,
R. A. Fallows,
W. Frieswijk,
M. C. Toribio,
B. Adebahr,
M. Arts,
M. R. Bell,
A. Bonafede,
J. Bray,
J. Broderick,
T. Cantwell,
P. Carroll,
Y. Cendes
, et al. (125 additional authors not shown)
Abstract:
We present the Multifrequency Snapshot Sky Survey (MSSS), the first northern-sky LOFAR imaging survey. In this introductory paper, we first describe in detail the motivation and design of the survey. Compared to previous radio surveys, MSSS is exceptional due to its intrinsic multifrequency nature providing information about the spectral properties of the detected sources over more than two octave…
▽ More
We present the Multifrequency Snapshot Sky Survey (MSSS), the first northern-sky LOFAR imaging survey. In this introductory paper, we first describe in detail the motivation and design of the survey. Compared to previous radio surveys, MSSS is exceptional due to its intrinsic multifrequency nature providing information about the spectral properties of the detected sources over more than two octaves (from 30 to 160 MHz). The broadband frequency coverage, together with the fast survey speed generated by LOFAR's multibeaming capabilities, make MSSS the first survey of the sort anticipated to be carried out with the forthcoming Square Kilometre Array (SKA). Two of the sixteen frequency bands included in the survey were chosen to exactly overlap the frequency coverage of large-area Very Large Array (VLA) and Giant Metrewave Radio Telescope (GMRT) surveys at 74 MHz and 151 MHz respectively. The survey performance is illustrated within the "MSSS Verification Field" (MVF), a region of 100 square degrees centered at J2000 (RA,Dec)=(15h,69deg). The MSSS results from the MVF are compared with previous radio survey catalogs. We assess the flux and astrometric uncertainties in the catalog, as well as the completeness and reliability considering our source finding strategy. We determine the 90% completeness levels within the MVF to be 100 mJy at 135 MHz with 108" resolution, and 550 mJy at 50 MHz with 166" resolution. Images and catalogs for the full survey, expected to contain 150,000-200,000 sources, will be released to a public web server. We outline the plans for the ongoing production of the final survey products, and the ultimate public release of images and source catalogs.
△ Less
Submitted 3 September, 2015;
originally announced September 2015.
-
Using SKA Rotation Measures to Reveal the Mysteries of the Magnetised Universe
Authors:
Melanie Johnston-Hollitt,
Federica Govoni,
Rainer Beck,
Siamak Dehghan,
Luke Pratley,
Takuya Akahori,
George Heald,
Ivan Agudo,
Annalisa Bonafede,
Ettore Carretti,
Tracy Clarke,
Sergio Colafrancesco,
Torsten Enßlin,
Luigina Feretti,
Bryan Gaensler,
Marijke Haverkorn,
Sui Ann Mao,
Niels Oppermann,
Lawrence Rudnick,
Anna Scaife,
Dominic Schnitzeler,
Jeroen Stil,
A. Russ Taylor,
Valentina Vacca
Abstract:
We know that magnetic fields are pervasive across all scales in the Universe and over all of cosmic time and yet our understanding of many of the properties of magnetic fields is still limited. We do not yet know when, where or how the first magnetic fields in the Universe were formed, nor do we fully understand their role in fundamental processes such as galaxy formation or cosmic ray acceleratio…
▽ More
We know that magnetic fields are pervasive across all scales in the Universe and over all of cosmic time and yet our understanding of many of the properties of magnetic fields is still limited. We do not yet know when, where or how the first magnetic fields in the Universe were formed, nor do we fully understand their role in fundamental processes such as galaxy formation or cosmic ray acceleration or how they influence the evolution of astrophysical objects. The greatest challenge to addressing these issues has been a lack of deep, broad bandwidth polarimetric data over large areas of the sky. The Square Kilometre Array will radically improve this situation via an all-sky polarisation survey that delivers both high quality polarisation imaging in combination with observations of 7-14 million extragalactic rotation measures. Here we summarise how this survey will improve our understanding of a range of astrophysical phenomena on scales from individual Galactic objects to the cosmic web.
△ Less
Submitted 2 June, 2015;
originally announced June 2015.
-
Using Tailed Radio Galaxies to Probe the Environment and Magnetic Field of Galaxy Clusters in the SKA Era
Authors:
M. Johnston-Hollitt,
S. Dehghan,
L. Pratley
Abstract:
The morphology of tailed radio galaxies is an invaluable source of environmental information, in which a history of the past interactions in the intra-cluster medium, such as complex galaxy motions and cluster merger shocks, are preserved. In recent years, the use of tailed radio galaxies as environmental probes has gained momentum as a method for galaxy cluster detection, examining the dynamics o…
▽ More
The morphology of tailed radio galaxies is an invaluable source of environmental information, in which a history of the past interactions in the intra-cluster medium, such as complex galaxy motions and cluster merger shocks, are preserved. In recent years, the use of tailed radio galaxies as environmental probes has gained momentum as a method for galaxy cluster detection, examining the dynamics of individual clusters, measuring the density and velocity flows in the intra-cluster medium, and for probing cluster magnetic fields. To date instrumental limitations in terms of resolution and sensitivity have confined this research to the local (z < 0.7) Universe. The advent of SKA1 surveys however will allow detection of roughly 1,000,000 tailed radio galaxies and their associated galaxy clusters out to redshifts of 2 or more. This is in fact ten times more than the current number of known clusters in the Universe. Additionally between 50,000 and 100,000 tailed radio galaxies will be sufficiently polarized to allow characterization of the magnetic field of their parent cluster. Such a substantial sample of tailed galaxies will provide an invaluable tool not only for detecting clusters, but also for characterizing their intra-cluster medium, magnetic fields and dynamical state as a function of cosmic time. In this chapter we present an analysis of the usability of tailed radio galaxies as tracers of dense environments extrapolated from existing deep radio surveys.
△ Less
Submitted 4 January, 2015;
originally announced January 2015.
-
Using radio jets of PKS J0334-3900 to probe the intra-cluster medium in A3135
Authors:
L. Pratley,
M. Johnston-Hollitt,
S. Dehghan,
M. Sun
Abstract:
We present a multi-wavelength study of the radio galaxy PKS J0334-3900, which resides at the centre of Abell 3135. Using Australia Telescope Compact Array (ATCA) observations at 1.4, 2.5, 4.6 & 8.6 GHz, we performed a detailed analysis of PKS J0334-3900. The morphology and spectral indices give physical parameters that constrain the dynamical history of the galaxy, which we use to produce a simula…
▽ More
We present a multi-wavelength study of the radio galaxy PKS J0334-3900, which resides at the centre of Abell 3135. Using Australia Telescope Compact Array (ATCA) observations at 1.4, 2.5, 4.6 & 8.6 GHz, we performed a detailed analysis of PKS J0334-3900. The morphology and spectral indices give physical parameters that constrain the dynamical history of the galaxy, which we use to produce a simulation of PKS J0334-3900. This simulation shows that the morphology can be generated by a wind in the intracluster medium (ICM), orbital motion caused by a companion galaxy, and precession of the black hole (BH).
△ Less
Submitted 30 December, 2014;
originally announced January 2015.
-
Using the morphology and magnetic fields of tailed radio galaxies as environmental probes
Authors:
M. Johnston-Hollitt,
S. Dehghan,
L. Pratley
Abstract:
Bent-tailed (BT) radio sources have long been known to trace over densities in the Universe up to z ~ 1 and there is increasing evidence this association persists out to redshifts of 2. The morphology of the jets in BT galaxies is primarily a function of the environment that they have resided in and so BTs provide invaluable clues as to their local conditions. Thus, not only can samples of BT gala…
▽ More
Bent-tailed (BT) radio sources have long been known to trace over densities in the Universe up to z ~ 1 and there is increasing evidence this association persists out to redshifts of 2. The morphology of the jets in BT galaxies is primarily a function of the environment that they have resided in and so BTs provide invaluable clues as to their local conditions. Thus, not only can samples of BT galaxies be used as signposts of large-scale structure, but are also valuable for obtaining a statistical measurement of properties of the intra-cluster medium including the presence of cluster accretion shocks & winds, and as historical anemometers, preserving the dynamical history of their surroundings in their jets. We discuss the use of BTs to unveil large-scale structure and provide an example in which a BT was used to unlock the dynamical history of its host cluster. In addition to their use as density and dynamical indicators, BTs are useful probes of the magnetic field on their environment on scales which are inaccessible to other methods. Here we discuss a novel way in which a particular sub-class of BTs, the so-called `corkscrew' galaxies might further elucidate the coherence lengths of the magnetic fields in their vicinity. Given that BTs are estimated to make up a large population in next generation surveys we posit that the use of jets in this way could provide a unique source of environmental information for clusters and groups up to z = 2.
△ Less
Submitted 30 December, 2014;
originally announced December 2014.
-
Estimating extragalactic Faraday rotation
Authors:
Niels Oppermann,
Henrik Junklewitz,
Maksim Greiner,
Torsten A. Enßlin,
Takuya Akahori,
Ettore Carretti,
Bryan M. Gaensler,
Ariel Goobar,
Lisa Harvey-Smith,
Melanie Johnston-Hollitt,
Luke Pratley,
Dominic H. F. M. Schnitzeler,
Jeroen M. Stil,
Valentina Vacca
Abstract:
(abridged) Observations of Faraday rotation for extragalactic sources probe magnetic fields both inside and outside the Milky Way. Building on our earlier estimate of the Galactic contribution, we set out to estimate the extragalactic contributions. We discuss the problems involved; in particular, we point out that taking the difference between the observed values and the Galactic foreground recon…
▽ More
(abridged) Observations of Faraday rotation for extragalactic sources probe magnetic fields both inside and outside the Milky Way. Building on our earlier estimate of the Galactic contribution, we set out to estimate the extragalactic contributions. We discuss the problems involved; in particular, we point out that taking the difference between the observed values and the Galactic foreground reconstruction is not a good estimate for the extragalactic contributions. We point out a degeneracy between the contributions to the observed values due to extragalactic magnetic fields and observational noise and comment on the dangers of over-interpreting an estimate without taking into account its uncertainty information. To overcome these difficulties, we develop an extended reconstruction algorithm based on the assumption that the observational uncertainties are accurately described for a subset of the data, which can overcome the degeneracy with the extragalactic contributions. We present a probabilistic derivation of the algorithm and demonstrate its performance using a simulation, yielding a high quality reconstruction of the Galactic Faraday rotation foreground, a precise estimate of the typical extragalactic contribution, and a well-defined probabilistic description of the extragalactic contribution for each data point. We then apply this reconstruction technique to a catalog of Faraday rotation observations. We vary our assumptions about the data, showing that the dispersion of extragalactic contributions to observed Faraday depths is most likely lower than 7 rad/m^2, in agreement with earlier results, and that the extragalactic contribution to an individual data point is poorly constrained by the data in most cases.
△ Less
Submitted 16 February, 2015; v1 submitted 14 April, 2014;
originally announced April 2014.
-
Using head-tail galaxies to constrain the intracluster magnetic field: an in-depth study of PKS J0334-3900
Authors:
Luke Pratley,
Melanie Johnston-Hollitt,
Siamak Dehghan,
Ming Sun
Abstract:
We present a multi-wavelength study of the radio galaxy PKS J0334-3900 at the centre of Abell 3135. The spectro-polarimetric radio observations are combined with spectroscopic optical and X-ray data to illustrate the use of Head-Tail radio galaxies to reveal properties of the intracluster medium. ATCA observations at 1.4, 2.5, 4.6 & 8.6 GHz are presented with a detailed analysis of the morphology…
▽ More
We present a multi-wavelength study of the radio galaxy PKS J0334-3900 at the centre of Abell 3135. The spectro-polarimetric radio observations are combined with spectroscopic optical and X-ray data to illustrate the use of Head-Tail radio galaxies to reveal properties of the intracluster medium. ATCA observations at 1.4, 2.5, 4.6 & 8.6 GHz are presented with a detailed analysis of the morphology and spectral indices giving physical parameters to constrain the dynamical history of the galaxy. Using these constraints we produce a simulation of PKS J0334-3900. We find that this Head-Tail morphology can be induced via a combination of orbital motion due to a binary companion and relative motion through the ICM. New Chandra images of A3135 are presented from which we obtain a cluster electron density of n_(e,0) = (1.06 +/- 0.11 x 10^(-3) cm^(-3), a global temperature of 2.4 ^(+0.51)_(-0.38) keV and a lower limit to the radio jet power of PKS J0334-3900 of 1.6 x 10^(44) erg/s. A new redshift analysis of the cluster from available spectroscopic data shows A3135 to be comprised of galaxies with 0.058 < z < 0.066 and gives a new mean cluster redshift of 0.06228 +/- 0.00015. We also uncovered a background subgroup between 0.066 < z < 0.070. Stokes Q and U data of Abell 3135 were used to obtain rotation measure values along the line of sight to PKS J0334-3900. Using our simulation we are able to infer the distance between the jets to be 154 +/- 16 kpc, which when combined with the difference in rotation measure between the jets provides a novel new way to estimate the average magnetic field within a cluster. A lower limit to the cluster B-field was calculated to be 0.09 +/- 0.03 micro Gauss. We show observations of Head-Tail galaxies can be used to infer information on the cluster environment, showing them to be an important class of objects in next generation all sky surveys.
△ Less
Submitted 5 April, 2013; v1 submitted 12 March, 2013;
originally announced March 2013.
-
An improved map of the Galactic Faraday sky
Authors:
N. Oppermann,
H. Junklewitz,
G. Robbers,
M. R. Bell,
T. A. Enßlin,
A. Bonafede,
R. Braun,
J. C. Brown,
T. E. Clarke,
I. J. Feain,
B. M. Gaensler,
A. Hammond,
L. Harvey-Smith,
G. Heald,
M. Johnston-Hollitt,
U. Klein,
P. P. Kronberg,
S. A. Mao,
N. M. McClure-Griffiths,
S. P. O'Sullivan,
L. Pratley,
T. Robishaw,
S. Roy,
D. H. F. M. Schnitzeler,
C. Sotomayor-Beltran
, et al. (6 additional authors not shown)
Abstract:
We aim to summarize the current state of knowledge regarding Galactic Faraday rotation in an all-sky map of the Galactic Faraday depth. For this we have assembled the most extensive catalog of Faraday rotation data of compact extragalactic polarized radio sources to date. In the map making procedure we use a recently developed algorithm that reconstructs the map and the power spectrum of a statist…
▽ More
We aim to summarize the current state of knowledge regarding Galactic Faraday rotation in an all-sky map of the Galactic Faraday depth. For this we have assembled the most extensive catalog of Faraday rotation data of compact extragalactic polarized radio sources to date. In the map making procedure we use a recently developed algorithm that reconstructs the map and the power spectrum of a statistically isotropic and homogeneous field while taking into account uncertainties in the noise statistics. This procedure is able to identify some rotation angles that are offset by an integer multiple of pi. The resulting map can be seen as an improved version of earlier such maps and is made publicly available, along with a map of its uncertainty. For the angular power spectrum we find a power law behavior with a power law index of -2.14 for a Faraday sky where an overall variance profile as a function of Galactic latitude has been removed, in agreement with earlier work. We show that this is in accordance with a 3D Fourier power spectrum P(k) proportional to k^-2.14 of the underlying field n_e times B_r under simplifying geometrical and statistical assumptions.
△ Less
Submitted 5 March, 2012; v1 submitted 26 November, 2011;
originally announced November 2011.