Radio Galaxy Shape Measurement with Hamiltonian Monte Carlo in the Visibility Domain

Abstract

Radio weak lensing, while a highly promising complementary probe to optical weak lensing, will require incredible precision in the measurement of galaxy shape parameters. In this paper, we extend the Bayesian Inference for Radio Observations model fitting approach to measure galaxy shapes directly from visibility data of radio continuum surveys, instead of from image data. We apply a Hamiltonian Monte Carlo (HMC) technique for sampling the posterior, which is more efficient than the standard Monte Carlo Markov Chain method when dealing with a large dimensional parameter space. Adopting the exponential profile for galaxy model fitting allows us to analytically calculate the likelihood gradient required by HMC, allowing a faster and more accurate sampling. The method is tested on SKA1-MID simulated observations at 1.4 GHz of a field containing up to 1000 star-forming galaxies. It is also applied to a simulated observation of the weak lensing precursor survey SuperCLASS. In both cases we obtain reliable measurements of the galaxies' ellipticity and size for all sources with SNR $\ge 10$, and we also find relationships between the convergence properties of the HMC technique and some source parameters. Direct shape measurement in the visibility domain achieves high accuracy at the expected source number densities of the current and next SKA precursor continuum surveys. The proposed method can be easily extended for the fitting of other galaxy and scientific parameters, as well as simultaneously marginalising over systematic and instrumental effects.

arXiv comment15 pages, 10 figures, 3 tables

References in corpus (35)

• A Conceptual Introduction to Hamiltonian Monte Carlo
• Bayesian Galaxy Shape Measurement for Weak Lensing Surveys - III. Application to the Canada-France-Hawaii Telescope Lensing Survey
• The VLA-COSMOS 3 GHz Large Project: Continuum data and source catalog release
• Shapelets: I. A Method for Image Analysis
• Revisiting the radio interferometer measurement equation. I. A full-sky Jones formalism
• The faint radio sky: radio astronomy becomes mainstream
• UVMULTIFIT: A versatile tool for fitting astronomical radio interferometric data
• The DES Science Verification Weak Lensing Shear Catalogues
• Shapelets: II. A Method for Weak Lensing Measurements
• The VLA-COSMOS 3~GHz Large Project: AGN and host-galaxy properties out to z$\lesssim$6
• GREAT3 results I: systematic errors in shear estimation and the impact of real galaxy morphology
• The sub-mJy radio sky in the Extended Chandra Deep Field South: source population
• The Deep SWIRE Field I. 20cm Continuum Radio Observations: A Crowded Sky
• The Tiered Radio Extragalactic Continuum Simulation (T-RECS)
• Revisiting the radio interferometer measurement equation. II. Calibration and direction-dependent effects
• The MeqTrees software system and its use for third-generation calibration of radio interferometers
• Galaxy Evolution in the Radio Band: The Role of Starforming Galaxies and Active Galactic Nuclei
• Fast optimal CMB power spectrum estimation with Hamiltonian sampling
• SKA Weak Lensing I: Cosmological Forecasts and the Power of Radio-Optical Cross-Correlations
• Weak Lensing by Large-Scale Structure with the FIRST Radio Survey
• Limitations for shapelet-based weak-lensing measurements
• Deep JVLA Imaging of GOODS-N at 20cm
• SERPent: Automated reduction and RFI-mitigation software for e-MERLIN
• Sérsic galaxy models in weak lensing shape measurement: model bias, noise bias and their interaction
• The eMERGE Survey I: Very Large Array 5.5 GHz observations of the GOODS-North Field
• Bayesian Inference for Radio Observations
• Radio Weak Gravitational Lensing with VLA and MERLIN
• Montblanc: GPU accelerated Radio Interferometer Measurement Equations in support of Bayesian Inference for Radio Observations
• Shape Reconstruction and Weak Lensing Measurement with Interferometers: A Shapelet Approach
• Analytical Galaxy Profiles for Photometric and Lensing Analysis
• A Submillimeter Perspective on the GOODS Fields (SUPER GOODS) - II. The High Radio Power Population in the GOODS-N
• Radio-Optical Galaxy Shape Correlations in the COSMOS Field
• Radio Weak Lensing Shear Measurement in the Visibility Domain - I. Methodology
• Radio Weak Lensing Shear Measurement in the Visibility Domain - II. Source Extraction
• SKA Engineering Change Proposal: Gridded Visibilities to Enable Precision Cosmology with Radio Weak Lensing

Cited by in corpus (7)

• Prospects for Recovering Galaxy Intrinsic Shapes from Projected Quantities
• SuperCLASS -- III. Weak lensing from radio and optical observations in Data Release 1
• SuperCLASS -- I. The Super CLuster Assisted Shear Survey: Project overview and Data Release 1
• Radio-Optical Galaxy Shape and Shear Correlations in the COSMOS Field using 3 GHz VLA Observations
• Galaxy Image Restoration with Shape Constraint
• Radio Galaxy Detection in the Visibility Domain
• RadioLensfit: an HPC Tool for Accurate Galaxy Shape Measurement with SKA