Source author record

Jonathan Zwart

Jonathan Zwart appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

4works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

4 published item(s)

preprint2016arXiv

Resolving the blazar CGRaBS J0809+5341 in the presence of telescope systematics

We analyse Very Long Baseline Interferometry (VLBI) observations of the blazar CGRaBS J0809+5341 using Bayesian inference methods. The observation was carried out at 5 GHz using 8 telescopes that form part of the European VLBI Network. Imaging and deconvolution using traditional methods imply that the blazar is unresolved. To search for source structure beyond the diffraction limit, we perform Bayesian model selection between three source models (point, elliptical Gaussian, and circular Gaussian). Our modelling jointly accounts for antenna-dependent gains and system equivalent flux densities. We obtain posterior distributions for the various source and instrumental parameters along with the corresponding uncertainties and correlations between them. We find that there is very strong evidence (>1e9 :1) in favour of elliptical Gaussian structure and using this model derive the apparent brightness temperature distribution of the blazar, accounting for uncertainties in the shape estimates. To test the integrity of our method, we also perform model selection on synthetic observations and use this to develop a Bayesian criterion for the minimum resolvable source size and consequently the maximum measurable brightness temperature for a given interferometer, dependent on the signal-to-noise ratio (SNR) of the data incorporating the aforementioned systematics. We find that calibration errors play an increasingly important role in determining the over-resolution limit for SNR>>100. We show that it is possible to exploit the resolving power of future VLBI arrays down to about 5 per cent of the size of the naturally-weighted restoring beam, if the gain calibration is precise to 1 per cent or less.

preprint2015arXiv

Montblanc: GPU accelerated Radio Interferometer Measurement Equations in support of Bayesian Inference for Radio Observations

We present Montblanc, a GPU implementation of the Radio interferometer measurement equation (RIME) in support of the Bayesian inference for radio observations (BIRO) technique. BIRO uses Bayesian inference to select sky models that best match the visibilities observed by a radio interferometer. To accomplish this, BIRO evaluates the RIME multiple times, varying sky model parameters to produce multiple model visibilities. Chi-squared values computed from the model and observed visibilities are used as likelihood values to drive the Bayesian sampling process and select the best sky model. As most of the elements of the RIME and chi-squared calculation are independent of one another, they are highly amenable to parallel computation. Additionally, Montblanc caters for iterative RIME evaluation to produce multiple chi-squared values. Modified model parameters are transferred to the GPU between each iteration. We implemented Montblanc as a Python package based upon NVIDIA's CUDA architecture. As such, it is easy to extend and implement different pipelines. At present, Montblanc supports point and Gaussian morphologies, but is designed for easy addition of new source profiles. Montblanc's RIME implementation is performant: On an NVIDIA K40, it is approximately 250 times faster than MeqTrees on a dual hexacore Intel E5-2620v2 CPU. Compared to the OSKAR simulator's GPU-implemented RIME components it is 7.7 and 12 times faster on the same K40 for single and double-precision floating point respectively. However, OSKAR's RIME implementation is more general than Montblanc's BIRO-tailored RIME. Theoretical analysis of Montblanc's dominant CUDA kernel suggests that it is memory bound. In practice, profiling shows that is balanced between compute and memory, as much of the data required by the problem is retained in L1 and L2 cache.

preprint2012arXiv

Future Science Prospects for AMI

The Arcminute Microkelvin Imager (AMI) is a telescope specifically designed for high sensitivity measurements of low-surface-brightness features at cm-wavelength and has unique, important capabilities. It consists of two interferometer arrays operating over 13.5-18 GHz that image structures on scales of 0.5-10 arcmin with very low systematics. The Small Array (AMI-SA; ten 3.7-m antennas) couples very well to Sunyaev-Zel'dovich features from galaxy clusters and to many Galactic features. The Large Array (AMI-LA; eight 13-m antennas) has a collecting area ten times that of the AMI-SA and longer baselines, crucially allowing the removal of the effects of confusing radio point sources from regions of low surface-brightness, extended emission. Moreover AMI provides fast, deep object surveying and allows monitoring of large numbers of objects. In this White Paper we review the new science - both Galactic and extragalactic - already achieved with AMI and outline the prospects for much more.

preprint2011arXiv

Non-standard morphological relic patterns in the cosmic microwave background

Statistically anomalous signals in the microwave background have been extensively studied in general in multipole space, and in real space mainly for circular and other simple patterns. In this paper we search for a range of non-trivial patterns in the temperature data from WMAP 7-year observations. We find a very significant detection of a number of such features and discuss their consequences for the essential character of the cosmos.