Probabilistic methods for radio interferometry data analysis

Probability theory provides a uniquely valid set of rules for plausible reasoning. This enables us to apply this mathematical formalism of probability, also known as Bayesian, with greater flexibility to problems of scientific inference. In this thesis, we are concerned with applying this method to...

Full description

Bibliographic Details
Main Author: Natarajan, Iniyan
Other Authors: Van Der Heyden, Kurt
Format: Doctoral Thesis
Language:English
Published: University of Cape Town 2018
Subjects:
Online Access:http://hdl.handle.net/11427/27243
id ndltd-netd.ac.za-oai-union.ndltd.org-uct-oai-localhost-11427-27243
record_format oai_dc
spelling ndltd-netd.ac.za-oai-union.ndltd.org-uct-oai-localhost-11427-272432020-07-22T05:07:40Z Probabilistic methods for radio interferometry data analysis Natarajan, Iniyan Van Der Heyden, Kurt Smirnov, Oleg M Zwart, Jonathan T L Astronomy Data analysis – methods Interferometric Probability theory provides a uniquely valid set of rules for plausible reasoning. This enables us to apply this mathematical formalism of probability, also known as Bayesian, with greater flexibility to problems of scientific inference. In this thesis, we are concerned with applying this method to the analysis of visibility data from radio interferometers. Any radio interferometry observation can be described using the Radio Interferometry Measurement Equation (RIME). Throughout the thesis, we use the RIME to model the visibilities in performing the probabilistic analysis. We first develop the theory for employing the RIME in performing Bayesian analysis of interferometric data. We then apply this to the problem of super-resolution with radio interferometers by performing model selection successfully between different source structures, all smaller in scale than the size of the point spread function (PSF) of the interferometer, on Westerbork Synthesis Radio Telescope (WSRT) simulations at a frequency of 1.4 GHz. We also quantify the change in the scale of the sources that can be resolved by WSRT at this frequency, with changing signal-to-noise (SNR) of the data, using simulations. Following this, we apply this method to a 5 GHz European VLBI Network (EVN) observation of the flaring blazar CGRaBS J0809+5341, to ascertain the presence of a jet emanating from its core, taking into account the imperfections in the station gain calibration performed on the data, especially on the longest baselines, prior to our analysis. We find that the extended source model is preferred over the point source model with an odds ratio of 109 : 1. Using the flux-density and shape parameter estimates of this model, we also derive the brightness temperature of the blazar (10¹¹-10¹² K), which confirms the presence of a relativistically boosted jet with an intrinsic brightness temperature lower than the apparent brightness temperature, consistent with the literature. We also develop a Bayesian criterion for super-resolution in the presence of baseline-dependent noise and calibration errors and find that these errors play an important role in determining how close one can get to the theoretical super-resolution limit. We then proceed to include fringe-fitting, the process of solving for the time and frequency dependent phase variations introduced by the interstellar medium and the Earth's atmosphere, in our probabilistic approach. Fringe-fitting is one of the first corrections made to Very Long Baseline Interferometry (VLBI) observations, and, by extending our method to include simultaneous fringefitting and source structure estimation, we will be able to perform end-to-end VLBI analysis using our method. To this end, we estimate source amplitude and fringe-fitting phase terms (phase offsets and delays) on 43 GHz Very Long Baseline Array and 230 GHz Event Horizon Telescope (EHT) simulations of point sources. We then perform model selection on a 5 μas extended Gaussian source (one-fourth the size of the PSF) on a synthetic 230 GHz EHT observation. Finally we incorporate turbulent time-varying phase offsets and delays in our model selection and show that the delays can be estimated to within 10-16 per cent error (often better than contemporary software packages) while simultaneously estimating the extended source structure. 2018-02-02T09:29:37Z 2018-02-02T09:29:37Z 2017 Doctoral Thesis Doctoral PhD http://hdl.handle.net/11427/27243 eng application/pdf University of Cape Town Faculty of Science Department of Astronomy
collection NDLTD
language English
format Doctoral Thesis
sources NDLTD
topic Astronomy
Data analysis – methods
Interferometric
spellingShingle Astronomy
Data analysis – methods
Interferometric
Natarajan, Iniyan
Probabilistic methods for radio interferometry data analysis
description Probability theory provides a uniquely valid set of rules for plausible reasoning. This enables us to apply this mathematical formalism of probability, also known as Bayesian, with greater flexibility to problems of scientific inference. In this thesis, we are concerned with applying this method to the analysis of visibility data from radio interferometers. Any radio interferometry observation can be described using the Radio Interferometry Measurement Equation (RIME). Throughout the thesis, we use the RIME to model the visibilities in performing the probabilistic analysis. We first develop the theory for employing the RIME in performing Bayesian analysis of interferometric data. We then apply this to the problem of super-resolution with radio interferometers by performing model selection successfully between different source structures, all smaller in scale than the size of the point spread function (PSF) of the interferometer, on Westerbork Synthesis Radio Telescope (WSRT) simulations at a frequency of 1.4 GHz. We also quantify the change in the scale of the sources that can be resolved by WSRT at this frequency, with changing signal-to-noise (SNR) of the data, using simulations. Following this, we apply this method to a 5 GHz European VLBI Network (EVN) observation of the flaring blazar CGRaBS J0809+5341, to ascertain the presence of a jet emanating from its core, taking into account the imperfections in the station gain calibration performed on the data, especially on the longest baselines, prior to our analysis. We find that the extended source model is preferred over the point source model with an odds ratio of 109 : 1. Using the flux-density and shape parameter estimates of this model, we also derive the brightness temperature of the blazar (10¹¹-10¹² K), which confirms the presence of a relativistically boosted jet with an intrinsic brightness temperature lower than the apparent brightness temperature, consistent with the literature. We also develop a Bayesian criterion for super-resolution in the presence of baseline-dependent noise and calibration errors and find that these errors play an important role in determining how close one can get to the theoretical super-resolution limit. We then proceed to include fringe-fitting, the process of solving for the time and frequency dependent phase variations introduced by the interstellar medium and the Earth's atmosphere, in our probabilistic approach. Fringe-fitting is one of the first corrections made to Very Long Baseline Interferometry (VLBI) observations, and, by extending our method to include simultaneous fringefitting and source structure estimation, we will be able to perform end-to-end VLBI analysis using our method. To this end, we estimate source amplitude and fringe-fitting phase terms (phase offsets and delays) on 43 GHz Very Long Baseline Array and 230 GHz Event Horizon Telescope (EHT) simulations of point sources. We then perform model selection on a 5 μas extended Gaussian source (one-fourth the size of the PSF) on a synthetic 230 GHz EHT observation. Finally we incorporate turbulent time-varying phase offsets and delays in our model selection and show that the delays can be estimated to within 10-16 per cent error (often better than contemporary software packages) while simultaneously estimating the extended source structure.
author2 Van Der Heyden, Kurt
author_facet Van Der Heyden, Kurt
Natarajan, Iniyan
author Natarajan, Iniyan
author_sort Natarajan, Iniyan
title Probabilistic methods for radio interferometry data analysis
title_short Probabilistic methods for radio interferometry data analysis
title_full Probabilistic methods for radio interferometry data analysis
title_fullStr Probabilistic methods for radio interferometry data analysis
title_full_unstemmed Probabilistic methods for radio interferometry data analysis
title_sort probabilistic methods for radio interferometry data analysis
publisher University of Cape Town
publishDate 2018
url http://hdl.handle.net/11427/27243
work_keys_str_mv AT natarajaniniyan probabilisticmethodsforradiointerferometrydataanalysis
_version_ 1719330655027855360