A wider audience: Turning VLBI into a survey instrument

Radio observations using the Very Long Baseline Interferometry (VLBI) technique typically have fields of view of only a few arcseconds, due to the computational problems inherent in imaging larger fields. Furthermore, sensitivity limitations restrict observations to very compact and bright objects, which are few and far between on the sky. Thus, while most branches of observational astronomy can carry out sensitive, wide-field surveys, VLBI observations are limited to targeted observations of carefully selected objects. However, recent advances in technology have made it possible to carry out the computations required to target hundreds of sources simultaneously. Furthermore, sensitivity upgrades have dramatically increased the number of objects accessible to VLBI observations. The combination of these two developments have enhanced the survey capabilities of VLBI observations such that it is now possible to observe (almost) any point in the sky with milli-arcsecond resolution. In this talk I review the development of wide-field VLBI, which has made significant progress over the last three years.Comment: Invited review at the General Assembly of the Astronomische Gesellschaf

LOFAR observations of 4C+19.44. On the discovery of low frequency spectral curvature in relativistic jet knots

We present the first LOFAR observations of the radio jet in the quasar 4C+19.44 (a.k.a. PKS 1354+19) obtained with the long baselines. The achieved resolution is very well matched to that of archival Jansky Very Large Array (JVLA) observations at higher radio frequencies as well as the archival X-ray images obtained with {\it Chandra}. We found that, for several knots along the jet, the radio flux densities measured at hundreds of MHz lie well below the values estimated by extrapolating the GHz spectra. This clearly indicates the presence of spectral curvature. Radio spectral curvature has been already observed in different source classes and/or extended radio structures and it has been often interpreted as due to intrinsic processes, as a curved particle energy distribution, rather than absorption mechanisms ({ Razin-Tsytovich} effect, free-free or synchrotron self absorption to name a few). Here we discuss our results according to the scenario where particles undergo stochastic acceleration mechanisms also in quasar jet knots.Comment: 13 pages, 4 tables, 4 figures, pre-proof version, published on the Astrophysical Journal (Harris, et al. 2019 ApJ, 873, 21

Multi-source self-calibration: Unveiling the microJy population of compact radio sources

Context. Very Long Baseline Interferometry (VLBI) data are extremely sensitive to the phase stability of the VLBI array. This is especially important when we reach {\mu}Jy r.m.s. sensitivities. Calibration using standard phase referencing techniques is often used to improve the phase stability of VLBI data but the results are often not optimal. This is evident in blank fields that do not have in-beam calibrators. Aims. We present a calibration algorithm termed Multi-Source Self-Calibration (MSSC) which can be used after standard phase referencing on wide-field VLBI observations. This is tested on a 1.6 GHz wide-field VLBI data set of the Hubble Deep Field-North and the Hubble Flanking Fields. Methods. MSSC uses multiple target sources detected in the field via standard phase referencing techniques and modifies the visibili- ties so that each data set approximates to a point source. These are combined to increase the signal to noise and permit self-calibration. In principle, this should allow residual phase changes caused by the troposphere and ionosphere to be corrected. By means of faceting, the technique can also be used for direction dependent calibration. Results. Phase corrections, derived using MSSC, were applied to a wide-field VLBI data set of the HDF-N comprising of 699 phase centres. MSSC was found to perform considerably better than standard phase referencing and single source self-calibration. All detected sources exhibited dramatic improvements in dynamic range. Using MSSC, one source reached the detection threshold taking the total detected sources to twenty. 60% of these sources can now be imaged with uniform weighting compared to just 45% with standard phase referencing. The Parseltongue code which implements MSSC has been released and made publicly available to the astronomical community (https://github.com/jradcliffe5/multi_self_cal).Comment: 7 pages, 4 figures, accepted to A&

Nowhere to Hide: Radio-faint AGN in the GOODS-N field. I. Initial catalogue and radio properties

(Abridged) Conventional radio surveys of deep fields ordinarily have arc-second scale resolutions often insufficient to reliably separate radio emission in distant galaxies originating from star-formation and AGN-related activity. Very long baseline interferometry (VLBI) can offer a solution by identifying only the most compact radio emitting regions in galaxies at cosmological distances where the high brightness temperatures (in excess of $10^5$ K) can only be reliably attributed to AGN activity. We present the first in a series of papers exploring the faint compact radio population using a new wide-field VLBI survey of the GOODS-N field. The unparalleled sensitivity of the European VLBI Network (EVN) will probe a luminosity range rarely seen in deep wide-field VLBI observations, thus providing insights into the role of AGN to radio luminosities of the order $10^{22}~\mathrm{W\,Hz^{-1}}$ across cosmic time. The newest VLBI techniques are used to completely cover an entire 7'.5 radius area to milliarcsecond resolutions, while bright radio sources ($S > 0.1$ mJy) are targeted up to 25 arcmin from the pointing centre. Multi-source self-calibration, and a primary beam model for the EVN array are used to correct for residual phase errors and primary beam attenuation respectively. This paper presents the largest catalogue of VLBI detected sources in GOODS-N comprising of 31 compact radio sources across a redshift range of 0.11-3.44, almost three times more than previous VLBI surveys in this field. We provide a machine-readable catalogue and introduce the radio properties of the detected sources using complementary data from the e-MERLIN Galaxy Evolution survey (eMERGE).Comment: 15 pages, 8 figures, accepted in A&A. Machine-readable table available upon reques

DiFX2: A more flexible, efficient, robust and powerful software correlator

Software correlation, where a correlation algorithm written in a high-level language such as C++ is run on commodity computer hardware, has become increasingly attractive for small to medium sized and/or bandwidth constrained radio interferometers. In particular, many long baseline arrays (which typically have fewer than 20 elements and are restricted in observing bandwidth by costly recording hardware and media) have utilized software correlators for rapid, cost-effective correlator upgrades to allow compatibility with new, wider bandwidth recording systems and improve correlator flexibility. The DiFX correlator, made publicly available in 2007, has been a popular choice in such upgrades and is now used for production correlation by a number of observatories and research groups worldwide. Here we describe the evolution in the capabilities of the DiFX correlator over the past three years, including a number of new capabilities, substantial performance improvements, and a large amount of supporting infrastructure to ease use of the code. New capabilities include the ability to correlate a large number of phase centers in a single correlation pass, the extraction of phase calibration tones, correlation of disparate but overlapping sub-bands, the production of rapidly sampled filterbank and kurtosis data at minimal cost, and many more. The latest version of the code is at least 15% faster than the original, and in certain situations many times this value. Finally, we also present detailed test results validating the correctness of the new code.Comment: 28 pages, 9 figures, accepted for publication in PAS

Detection of fast radio transients with multiple stations: a case study using the Very Long Baseline Array

Recent investigations reveal an important new class of transient radio phenomena that occur on sub-millisecond timescales. Often transient surveys' data volumes are too large to archive exhaustively. Instead, an on-line automatic system must excise impulsive interference and detect candidate events in real-time. This work presents a case study using data from multiple geographically distributed stations to perform simultaneous interference excision and transient detection. We present several algorithms that incorporate dedispersed data from multiple sites, and report experiments with a commensal real-time transient detection system on the Very Long Baseline Array (VLBA). We test the system using observations of pulsar B0329+54. The multiple-station algorithms enhanced sensitivity for detection of individual pulses. These strategies could improve detection performance for a future generation of geographically distributed arrays such as the Australian Square Kilometre Array Pathfinder and the Square Kilometre Array.Comment: 12 pages, 14 figures. Accepted for Ap

A new approach to modelling impacts on rubble pile asteroid simulants

Many asteroids with low bulk densities must have a rubble pile structure and internal voids. Although little is known about their internal structure, numerical simulations of impact events on these asteroids rely on assumptions on how the voids are distributed. We present a new approach to model impacts on rubble pile asteroids that explicitly takes into account their internal structure. The formation of the asteroid is modelled as a rubble pile aggregate of spherical pebbles of different sizes. This aggregate is then converted into a high-resolution smoothed particle hydrodynamics (SPH) model, accounting for macroporosity inside the pebbles. We compare impact-event outcomes for a large set of internal configurations to explore the parameter space of our model-building process. The analysis of the fragment size distribution and the disruption threshold quantifies the specific influence of each input parameter. The size distribution of the pebbles used in our model is a simple power law, containing three free parameters: the slope α, the lower cut-off radius rmin and the upper cut-off radius rmax. The influence of all three parameters on the outcome is assessed in this paper. The existence of void space in our model increases the resistance against collisional disruption, a behaviour previously reported based on numerical simulations using a continuum description of porous material (Holsapple 2009). We show, for a set of asteroid collisions typical for small asteroids in the main belt, that no a priori knowledge of the exact size distribution of the pebbles inside the asteroid is needed, as the choice of the corresponding parameters does not directly correlate with the impact outcome

Microarcsecond VLBI pulsar astrometry with PSRπ\pi II. parallax distances for 57 pulsars

We present the results of PSR$\pi$, a large astrometric project targeting radio pulsars using the Very Long Baseline Array (VLBA). From our astrometric database of 60 pulsars, we have obtained parallax-based distance measurements for all but 3, with a parallax precision of typically 40 $\mu$as and approaching 10 $\mu$as in the best cases. Our full sample doubles the number of radio pulsars with a reliable ($\gtrsim$5$\sigma$) model-independent distance constraint. Importantly, many of the newly measured pulsars are well outside the solar neighbourhood, and so PSR$\pi$ brings a near-tenfold increase in the number of pulsars with a reliable model-independent distance at $d>2$ kpc. Using our sample along with previously published results, we show that even the most recent models of the Galactic electron density distribution model contain significant shortcomings, particularly at high Galactic latitudes. When comparing our results to pulsar timing, two of the four millisecond pulsars in our sample exhibit significant discrepancies in the estimates of proper motion obtained by at least one pulsar timing array. With additional VLBI observations to improve the absolute positional accuracy of our reference sources and an expansion of the number of millisecond pulsars, we will be able to extend the comparison of proper motion discrepancies to a larger sample of pulsar reference positions, which will provide a much more sensitive test of the applicability of the solar system ephemerides used for pulsar timing. Finally, we use our large sample to estimate the typical accuracy attainable for differential astrometry with the VLBA when observing pulsars, showing that for sufficiently bright targets observed 8 times over 18 months, a parallax uncertainty of 4 $\mu$as per arcminute of separation between the pulsar and calibrator can be expected.Comment: updated to version accepted by ApJ: 30 pages, 20 figures, 9 table

Microarcsecond VLBI pulsar astrometry with PSRPI I. Two binary millisecond pulsars with white dwarf companions

Model-independent distance constraints to binary millisecond pulsars (MSPs) are of great value to both the timing observations of the radio pulsars, and multiwavelength observations of their companion stars. Very Long Baseline Interferometry (VLBI) astrometry can be employed to provide these model-independent distances with very high precision via the detection of annual geometric parallax. Using the Very Long Baseline Array, we have observed two binary millisecond pulsars, PSR J1022+1001 and J2145-0750, over a two-year period and measured their distances to be 700 +14 -10 pc and 613 +16 -14 pc respectively. We use the well-calibrated distance in conjunction with revised analysis of optical photometry to tightly constrain the nature of their massive (M ~ 0.85 Msun) white dwarf companions. Finally, we show that several measurements of their parallax and proper motion of PSR J1022+1001 and PSR J2145-0750 obtained by pulsar timing array projects are incorrect, differing from the more precise VLBI values by up to 5 sigma. We investigate possible causes for the discrepancy, and find that imperfect modeling of the solar wind is a likely candidate for the timing model errors given the low ecliptic latitude of these two pulsars.Comment: 14 pages, 9 figures, 6 tables; minor revisions in response to referee comments to match version accepted by Ap