skip to main content


Title: CS-ROMER: a novel compressed sensing framework for Faraday depth reconstruction
ABSTRACT

The reconstruction of Faraday depth structure from incomplete spectral polarization radio measurements using the RM synthesis technique is an underconstrained problem requiring additional regularization. In this paper, we present cs-romer: a novel object-oriented compressed sensing framework to reconstruct Faraday depth signals from spectropolarization radio data. Unlike previous compressed sensing applications, this framework is designed to work directly with data that are irregularly sampled in wavelength-squared space and to incorporate multiple forms of compressed sensing regularization. We demonstrate the framework using simulated data for the VLA telescope under a variety of observing conditions, and we introduce a methodology for identifying the optimal basis function for reconstruction of these data, using an approach that can also be applied to data sets from other telescopes and over different frequency ranges. In this work, we show that the delta basis function provides optimal reconstruction for VLA L-band data and we use this basis with observations of the low-mass galaxy cluster Abell 1314 in order to reconstruct the Faraday depth of its constituent cluster galaxies. We use the cs-romer framework to de-rotate the Galactic Faraday depth contribution directly from the wavelength-squared data and to handle the spectral behaviour of different radio sources in a direction-dependent manner. The results of this analysis show that individual galaxies within Abell 1314 deviate from the behaviour expected for a Faraday-thin screen such as the intra-cluster medium and instead suggest that the Faraday rotation exhibited by these galaxies is dominated by their local environments.

 
more » « less
NSF-PAR ID:
10381795
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
518
Issue:
2
ISSN:
0035-8711
Format(s):
Medium: X Size: p. 1955-1974
Size(s):
["p. 1955-1974"]
Sponsoring Org:
National Science Foundation
More Like this
  1. MeerKAT’s large number (64) of 13.5 m diameter antennas, spanning 8 km with a densely packed 1 km core, create a powerful instrument for wide-area surveys, with high sensitivity over a wide range of angular scales. The MeerKAT Galaxy Cluster Legacy Survey (MGCLS) is a programme of long-track MeerKAT L -band (900−1670 MHz) observations of 115 galaxy clusters, observed for ∼6−10 h each in full polarisation. The first legacy product data release (DR1), made available with this paper, includes the MeerKAT visibilities, basic image cubes at ∼8″ resolution, and enhanced spectral and polarisation image cubes at ∼8″ and 15″ resolutions. Typical sensitivities for the full-resolution MGCLS image products range from ∼3−5 μJy beam −1 . The basic cubes are full-field and span 2° × 2°. The enhanced products consist of the inner 1.2° × 1.2° field of view, corrected for the primary beam. The survey is fully sensitive to structures up to ∼10′ scales, and the wide bandwidth allows spectral and Faraday rotation mapping. Relatively narrow frequency channels (209 kHz) are also used to provide H  I mapping in windows of 0 <  z  < 0.09 and 0.19 <  z  < 0.48. In this paper, we provide an overview of the survey and the DR1 products, including caveats for usage. We present some initial results from the survey, both for their intrinsic scientific value and to highlight the capabilities for further exploration with these data. These include a primary-beam-corrected compact source catalogue of ∼626 000 sources for the full survey and an optical and infrared cross-matched catalogue for compact sources in the primary-beam-corrected areas of Abell 209 and Abell S295. We examine dust unbiased star-formation rates as a function of cluster-centric radius in Abell 209, extending out to 3.5 R 200 . We find no dependence of the star-formation rate on distance from the cluster centre, and we observe a small excess of the radio-to-100 μm flux ratio towards the centre of Abell 209 that may reflect a ram pressure enhancement in the denser environment. We detect diffuse cluster radio emission in 62 of the surveyed systems and present a catalogue of the 99 diffuse cluster emission structures, of which 56 are new. These include mini-halos, halos, relics, and other diffuse structures for which no suitable characterisation currently exists. We highlight some of the radio galaxies that challenge current paradigms, such as trident-shaped structures, jets that remain well collimated far beyond their bending radius, and filamentary features linked to radio galaxies that likely illuminate magnetic flux tubes in the intracluster medium. We also present early results from the H  I analysis of four clusters, which show a wide variety of H  I mass distributions that reflect both sensitivity and intrinsic cluster effects, and the serendipitous discovery of a group in the foreground of Abell 3365. 
    more » « less
  2. null (Ed.)
    ABSTRACT The SuperCLuster Assisted Shear Survey (SuperCLASS) is a legacy programme using the e-MERLIN interferometric array. The aim is to observe the sky at L-band (1.4 GHz) to a r.m.s. of $7\, \mu {\rm Jy}\,$beam−1 over an area of $\sim 1\, {\rm deg}^2$ centred on the Abell 981 supercluster. The main scientific objectives of the project are: (i) to detect the effects of weak lensing in the radio in preparation for similar measurements with the Square Kilometre Array (SKA); (ii) an extinction free census of star formation and AGN activity out to z ∼ 1. In this paper we give an overview of the project including the science goals and multiwavelength coverage before presenting the first data release. We have analysed around 400 h of e-MERLIN data allowing us to create a Data Release 1 (DR1) mosaic of $\sim 0.26\, {\rm deg}^2$ to the full depth. These observations have been supplemented with complementary radio observations from the Karl G. Jansky Very Large Array (VLA) and optical/near infrared observations taken with the Subaru, Canada-France-Hawaii, and Spitzer Telescopes. The main data product is a catalogue of 887 sources detected by the VLA, of which 395 are detected by e-MERLIN and 197 of these are resolved. We have investigated the size, flux, and spectral index properties of these sources finding them compatible with previous studies. Preliminary photometric redshifts, and an assessment of galaxy shapes measured in the radio data, combined with a radio-optical cross-correlation technique probing cosmic shear in a supercluster environment, are presented in companion papers. 
    more » « less
  3. Abstract Background

    Spectral CT material decomposition provides quantitative information but is challenged by the instability of the inversion into basis materials. We have previously proposed the constrained One‐Step Spectral CT Image Reconstruction (cOSSCIR) algorithm to stabilize the material decomposition inversion by directly estimating basis material images from spectral CT data. cOSSCIR was previously investigated on phantom data.

    Purpose

    This study investigates the performance of cOSSCIR using head CT datasets acquired on a clinical photon‐counting CT (PCCT) prototype. This is the first investigation of cOSSCIR for large‐scale, anatomically complex, clinical PCCT data. The cOSSCIR decomposition is preceded by a spectrum estimation and nonlinear counts correction calibration step to address nonideal detector effects.

    Methods

    Head CT data were acquired on an early prototype clinical PCCT system using an edge‐on silicon detector with eight energy bins. Calibration data of a step wedge phantom were also acquired and used to train a spectral model to account for the source spectrum and detector spectral response, and also to train a nonlinear counts correction model to account for pulse pileup effects. The cOSSCIR algorithm optimized the bone and adipose basis images directly from the photon counts data, while placing a grouped total variation (TV) constraint on the basis images. For comparison, basis images were also reconstructed by a two‐step projection‐domain approach of Maximum Likelihood Estimation (MLE) for decomposing basis sinograms, followed by filtered backprojection (MLE + FBP) or a TV minimization algorithm (MLE + TVmin) to reconstruct basis images. We hypothesize that the cOSSCIR approach will provide a more stable inversion into basis images compared to two‐step approaches. To investigate this hypothesis, the noise standard deviation in bone and soft‐tissue regions of interest (ROIs) in the reconstructed images were compared between cOSSCIR and the two‐step methods for a range of regularization constraint settings.

    Results

    cOSSCIR reduced the noise standard deviation in the basis images by a factor of two to six compared to that of MLE + TVmin, when both algorithms were constrained to produce images with the same TV. The cOSSCIR images demonstrated qualitatively improved spatial resolution and depiction of fine anatomical detail. The MLE + TVminalgorithm resulted in lower noise standard deviation than cOSSCIR for the virtual monoenergetic images (VMIs) at higher energy levels and constraint settings, while the cOSSCIR VMIs resulted in lower noise standard deviation at lower energy levels and overall higher qualitative spatial resolution. There were no statistically significant differences in the mean values within the bone region of images reconstructed by the studied algorithms. There were statistically significant differences in the mean values within the soft‐tissue region of the reconstructed images, with cOSSCIR producing mean values closer to the expected values.

    Conclusions

    The cOSSCIR algorithm, combined with our previously proposed spectral model estimation and nonlinear counts correction method, successfully estimated bone and adipose basis images from high resolution, large‐scale patient data from a clinical PCCT prototype. The cOSSCIR basis images were able to depict fine anatomical details with a factor of two to six reduction in noise standard deviation compared to that of the MLE + TVmintwo‐step approach.

     
    more » « less
  4. ABSTRACT

    Hybrid Morphology Radio Sources (HyMoRS) are a very rare subclass of radio galaxies with apparent mixed FR morphology, i.e. these galaxies seem to have an FR-I structure on one side of the core and an FR-II structure on the other side of the core. We systematically searched for HyMoRS using Very Large Array (VLA) Faint Images of the Radio Sky at Twenty-cm (FIRST) survey with 1400 MHz frequency and identified 33 candidate HyMoRS. Our finding significantly increased the known sample size of HyMoRS. HyMoRS may play an essential role in understanding the interaction of jets with the interstellar medium and the much-debated topic of the FR dichotomy. We identified optical counterparts for 29 sources in our catalogue. In our sample of sources, one source (J1106+1355) had quasar-like behaviour. Four sources were BRCLG (Brightest Cluster Galaxies) and six were LRG (Luminous Red Galaxies). We have estimated the spectral index and radio luminosity of HyMoRS in our catalogue, when possible. We found that J1136–0328 was the most luminous source in our sample (log L = 27.01 W Hz−1sr−1). It was also the farthest HyMoRS (with a redshift z = 0.82) in our sample. With the help of a large sample size of discovered sources, various statistical properties of detected galaxies were studied.

     
    more » « less
  5. Abstract We present a detailed analysis of the radio galaxy PKS $2250{-}351$ , a giant of 1.2 Mpc projected size, its host galaxy, and its environment. We use radio data from the Murchison Widefield Array, the upgraded Giant Metre-wavelength Radio Telescope, the Australian Square Kilometre Array Pathfinder, and the Australia Telescope Compact Array to model the jet power and age. Optical and IR data come from the Galaxy And Mass Assembly (GAMA) survey and provide information on the host galaxy and environment. GAMA spectroscopy confirms that PKS $2250{-}351$ lies at $z=0.2115$ in the irregular, and likely unrelaxed, cluster Abell 3936. We find its host is a massive, ‘red and dead’ elliptical galaxy with negligible star formation but with a highly obscured active galactic nucleus dominating the mid-IR emission. Assuming it lies on the local M – $\sigma$ relation, it has an Eddington accretion rate of $\lambda_{\rm EDD}\sim 0.014$ . We find that the lobe-derived jet power (a time-averaged measure) is an order of magnitude greater than the hotspot-derived jet power (an instantaneous measure). We propose that over the lifetime of the observed radio emission ( ${\sim} 300\,$ Myr), the accretion has switched from an inefficient advection-dominated mode to a thin disc efficient mode, consistent with the decrease in jet power. We also suggest that the asymmetric radio morphology is due to its environment, with the host of PKS $2250{-}351$ lying to the west of the densest concentration of galaxies in Abell 3936. 
    more » « less