The pulse morphology of fast radio bursts (FRBs) provides key information in both understanding progenitor physics and the plasma medium through which the burst propagates. We present a study of the profiles of 33 bright FRBs detected by the Australian Square Kilometre Array Pathfinder. We identify seven FRBs with measureable intrinsic pulse widths, including two FRBs that have been seen to repeat. In our modest sample we see no evidence for bimodality in the pulse width distribution. We also identify five FRBs with evidence of millisecond timescale pulse broadening caused by scattering in inhomogeneous plasma. We find no evidence for a relationship between pulse broadening and extragalactic dispersion measure. The scattering could be either caused by extreme turbulence in the host galaxy or chance propagation through foreground galaxies. With future high time resolution observations and detailed study of host galaxy properties we may be able to probe line-of-sight turbulence on gigaparsec scales.

Gravitational waves from coalescing neutron stars encode information about nuclear matter at extreme densities, inaccessible by laboratory experiments. The late inspiral is influenced by the presence of tides, which depend on the neutron star equation of state. Neutron star mergers are expected to often produce rapidly-rotating remnant neutron stars that emit gravitational waves. These will provide clues to the extremely hot post-merger environment. This signature of nuclear matter in gravitational waves contains most information in the 2-4 kHz frequency band, which is outside of the most sensitive band of current detectors. We present the design concept and science case for a neutron star extreme matter observatory (NEMO): a gravitational-wave interferometer optimized to study nuclear physics with merging neutron stars. The concept uses high circulating laser power, quantum squeezing and a detector topology specifically designed to achieve the high-frequency sensitivity necessary to probe nuclear matter using gravitational waves. Above one kHz, the proposed strain sensitivity is comparable to full third-generation detectors at a fraction of the cost. Such sensitivity changes expected event rates for detection of post-merger remnants from approximately one per few decades with two A+ detectors to a few per year, and potentially allows for the first gravitational-wave observations of supernovae, isolated neutron stars, and other exotica.

This white paper describes the science case for Very Long Baseline Interferometry (VLBI) and provides suggestions towards upgrade paths for the European VLBI Network (EVN). The EVN is a distributed long-baseline radio interferometric array, that operates at the very forefront of astronomical research. Recent results, together with the new science possibilities outlined in this vision document, demonstrate the EVN's potential to generate new and exciting results that will transform our view of the cosmos. Together with e-MERLIN, the EVN provides a range of baseline lengths that permit unique studies of faint radio sources to be made over a wide range of spatial scales. The science cases are reviewed in six chapters that cover the following broad areas: cosmology, galaxy formation and evolution, innermost regions of active galactic nuclei, explosive phenomena and transients, stars and stellar masers in the Milky Way, celestial reference frames and space applications. The document concludes with identifying the synergies with other radio, as well as multi-band/multi-messenger instruments, and provide the recommendations for future improvements. The appendices briefly describe other radio VLBI arrays, the technological framework for EVN developments, and a selection of spectral lines of astrophysical interest below 100 GHz. The document includes a glossary for non-specialists, and a list of acronyms at the end.

Fast radio bursts (FRBs) are millisecond-scale radio pulses, which originate in distant galaxies and are produced by unknown sources. The mystery remains partially because of the typical difficulty in localising FRBs to host galaxies. Accurate localisations delivered by the Commensal Real-time ASKAP Fast Transients (CRAFT) survey now provide an opportunity to study the host galaxies and potential transient counterparts of FRBs at a large range of wavelengths. In this work, we investigate whether the first three FRBs accurately localised by CRAFT have supernova-like transient counterparts. We obtained two sets of imaging epochs with the Very Large Telescope for three host galaxies, one soon after the burst detection and one several months later. After subtracting these images no optical counterparts were identified in the associated FRB host galaxies, so we instead place limits on the brightness of any potential optical transients. A Monte Carlo approach, in which supernova light curves were modelled and their base properties randomised, was used to estimate the probability of a supernova associated with each FRB going undetected. We conclude that Type Ia and IIn supernovae are unlikely to accompany every apparently non-repeating FRB.

Despite existing constraints, it remains possible that up to $35\%$ of all dark matter is comprised of compact objects, such as the black holes in the 10-100\,M$_\odot$ range whose existence has been confirmed by LIGO. The strong gravitational lensing of transients such as FRBs and GRBs has been suggested as a more sensitive probe for compact dark matter than intensity fluctuations observed in microlensing experiments. Recently ASKAP has reported burst substructure down to $15\mu$s timescales in FRBs in the redshift range $0.3-0.5$. We investigate here the implications of this for the detectability of compact dark matter by FRBs. We find that a sample size of $\sim130$ FRBs would be required to constrain compact dark matter to less than the existing 35$\%$ limit with 95$\%$ confidence, if it were distributed along $\gtrsim 1\,$Gpc-long FRB sightlines through the cosmic web. Conversely, existing constraints on the fraction of compact dark matter permit as many as 1 in $\approx 40$ of all $z \lesssim 0.4$ FRBs to exhibit micro-lensed burst structure. Approximately $170$ FRBs intercepting halos within $\sim 50\,$kpc would be required to exclude the fraction of compact dark matter in each intercepted halo to a similar level. Furthermore, we consider the cumulative effects of lensing of the FRB signal by a macroscopic dark matter distribution. We conclude that lensing from a uniform distribution of compact objects is likely not observable, but suggest that FRBs may set meaningful limits on power-law distributions of dark matter.

We present the full panchromatic afterglow light curve data of GW170817, including new radio data as well as archival optical and X-ray data, between 0.5 and 940 days post-merger. By compiling all archival data, and reprocessing a subset of it, we have ensured that the panchromatic dataset is uniform and therefore immune to the differences in data processing or flux determination methods used by different groups. Simple power-law fits to the uniform afterglow light curve indicate a $t^{0.86\pm0.04}$ rise, a $t^{-1.90\pm0.12}$ decline, and a peak occurring at $155\pm4$ days. The afterglow is optically thin throughout its evolution, consistent with a single spectral index ($-0.569\pm0.002$) across all epochs. This gives a precise and updated estimate of the electron power-law index, $p=2.138\pm0.004$. By studying the diffuse X-ray emission from the host galaxy, we place a conservative upper limit on the hot ionized ISM density, $<$0.01 cm$^{-3}$, consistent with previous afterglow studies. Using the late-time afterglow data we rule out any long-lived neutron star remnant having magnetic field strength between 10$^{10.4}$ G and 10$^{16}$ G. Our fits to the afterglow data using an analytical model that includes VLBI proper motion from Mooley et al (2018), and a structured jet model that ignores the proper motion, indicates that the proper motion measurement needs to be considered while seeking an accurate estimate of the viewing angle.

We present a high-resolution analysis of the host galaxy of fast radio burst FRB 190608, an SBc galaxy at $z=0.11778$ (hereafter HG 190608), to dissect its local environment and its contributions to the FRB properties. Our Hubble Space Telescope WFC3/UVIS image reveals that the sub-arcsecond localization of FRB 190608 is coincident with a knot of star-formation ($\Sigma_{SFR} = 1.2 \times 10^{-2}~ M_{\odot} \, kpc^{-2}$) in one of the prominent spiral arms of HG 190608. This is confirmed by H$\beta$ emission present in our Keck/KCWI integral field spectrum of the galaxy with a surface brightness of $\mu_{H\beta} = (3.35\pm0.18)\times10^{-17}\;erg\;s^{-1}\;cm^{-2}\;arcsec^{-2}$. We infer an extinction-corrected H$\alpha$ surface brightness and compute a dispersion measure from the interstellar medium of HG 190608 of ${DM}_{Host,ISM} = 82 \pm 35~ pc \, cm^{-3}$. The galaxy rotates with a circular velocity $v_{circ} = 141 \pm 8~ km \, s^{-1}$ at an inclination $i_{gas} = 37 \pm 3^\circ$, giving a dynamical mass $M_{halo}^{dyn} \approx 10^{11.96 \pm 0.08}~ M_{\odot}$. A surface photometric analysis of the galaxy using FORS2 imaging suggests a stellar disk inclination of $i_{stellar} = 26 \pm 3^\circ$. The dynamical mass estimate implies a halo contribution to the dispersion measure of ${DM}_{Host,Halo} = 55 \pm 25\; pc \, cm^{-3}$ subject to assumptions on the density profile and fraction of baryons retained. The relatively high temporal broadening ($\tau = 3.3 \pm 0.2 \; ms$ at 1.28 GHz) and rotation measure ($ RM = 353 \pm 2\; rad \; m^{-2}$) (Day et al. 2020) of FRB 190608 may be attributable to both turbulent gas within the spiral arm and gas local to the FRB progenitor. In contrast to previous high-resolution studies of FRB progenitor environments, we find no evidence for disturbed morphology, emission, nor kinematics for FRB 190608.

FRB 190608 was detected by ASKAP and localized to a spiral galaxy at $z_{host}=0.11778$ in the SDSS footprint. The burst has a large dispersion measure ($DM_{FRB}=339.8$ $pc/cm^3$) compared to the expected cosmic average at its redshift. It also has a large rotation measure ($RM_{FRB}=353$ $rad/m^2$) and scattering timescale ($\tau=3.3$ $ms$ at $1.28$ $GHz$). Chittidi et al (2020) perform a detailed analysis of the ultraviolet and optical emission of the host galaxy and estimate the host DM contribution to be $137\pm 43$ $pc/cm^3$. This work complements theirs and reports the analysis of the optical data of galaxies in the foreground of FRB 190608 to explore their contributions to the FRB signal. Together, the two manuscripts delineate an observationally driven, end-to-end study of matter distribution along an FRB sightline; the first study of its kind. Combining KCWI observations and public SDSS data, we estimate the expected cosmic dispersion measure $DM_{cosmic}$ along the sightline to FRB 190608. We first estimate the contribution of hot, ionized gas in intervening virialized halos ($DM_{halos} \approx 7-28$ $pc/cm^3$). Then, using the Monte Carlo Physarum Machine (MCPM) methodology, we produce a 3D map of ionized gas in cosmic web filaments and compute the DM contribution from matter outside halos ($DM_{IGM} \approx 91-126$ $pc/cm^3$). This implies a greater fraction of ionized gas along this sightline is extant outside virialized halos. We also investigate whether the intervening halos can account for the large FRB rotation measure and pulse width and conclude that it is implausible. Both the pulse broadening and the large Faraday rotation likely arise from the progenitor environment or the host galaxy.

More than three quarters of the baryonic content of the Universe resides in a highly diffuse state that is difficult to observe, with only a small fraction directly observed in galaxies and galaxy clusters. Censuses of the nearby Universe have used absorption line spectroscopy to observe these invisible baryons, but these measurements rely on large and uncertain corrections and are insensitive to the majority of the volume, and likely mass. Specifically, quasar spectroscopy is sensitive either to only the very trace amounts of Hydrogen that exists in the atomic state, or highly ionized and enriched gas in denser regions near galaxies. Sunyaev-Zel'dovich analyses provide evidence of some of the gas in filamentary structures and studies of X-ray emission are most sensitive to gas near galaxy clusters. Here we report the direct measurement of the baryon content of the Universe using the dispersion of a sample of localized fast radio bursts (FRBs), thus utilizing an effect that measures the electron column density along each sight line and accounts for every ionised baryon. We augment the sample of published arcsecond-localized FRBs with a further four new localizations to host galaxies which have measured redshifts of 0.291, 0.118, 0.378 and 0.522, completing a sample sufficiently large to account for dispersion variations along the line of sight and in the host galaxy environment to derive a cosmic baryon density of $\Omega_{b} = 0.051_{-0.025}^{+0.021} \, h_{70}^{-1}$ (95% confidence). This independent measurement is consistent with Cosmic Microwave Background and Big Bang Nucleosynthesis values.

Combining high time and frequency resolution full-polarisation spectra of Fast Radio Bursts (FRBs) with knowledge of their host galaxy properties provides an opportunity to study both the emission mechanism generating them and the impact of their propagation through their local environment, host galaxy, and the intergalactic medium. The Australian Square Kilometre Array Pathfinder (ASKAP) telescope has provided the first ensemble of bursts with this information. In this paper, we present the high time and spectral resolution, full polarisation observations of five localised FRBs to complement the results published for the previously studied ASKAP FRB~181112. We find that every FRB is highly polarised, with polarisation fractions ranging from 80 -- 100\%, and that they are generally dominated by linear polarisation. While some FRBs in our sample exhibit properties associated with an emerging archetype (i.e., repeating or apparently non-repeating), others exhibit characteristic features of both, implying the existence of a continuum of FRB properties. When examined at high time resolution, we find that all FRBs in our sample have evidence for multiple sub-components and for scattering at a level greater than expected from the Milky Way. We find no correlation between the diverse range of FRB properties (e.g., scattering time, intrinsic width, and rotation measure) and any global property of their host galaxy. The most heavily scattered bursts reside in the outskirts of their host galaxies, suggesting that the source-local environment rather than the host interstellar medium is likely the dominant origin of the scattering in our sample.

The Australian SKA Pathfinder (ASKAP) telescope has started to localize Fast Radio Bursts (FRBs) to arcsecond accuracy from the detection of a single pulse, allowing their host galaxies to be reliably identified. We discuss the global properties of the host galaxies of the first four FRBs localized by ASKAP, which lie in the redshift range $0.11<z<0.48$. All four are massive galaxies (log( $M_{*}/ M_{\odot}$) $\sim 9.4 -10.4$) with modest star-formation rates of up to $2M_{\odot}$yr$^{-1}$ -- very different to the host galaxy of the first repeating FRB 121102, which is a dwarf galaxy with a high specific star-formation rate. The FRBs localized by ASKAP typically lie in the outskirts of their host galaxies, which appears to rule out FRB progenitor models that invoke active galactic nuclei (AGN) or free-floating cosmic strings. The stellar population seen in these host galaxies also disfavors models in which all FRBs arise from young magnetars produced by superluminous supernovae (SLSNe), as proposed for the progenitor of FRB 121102. A range of other progenitor models (including compact-object mergers and magnetars arising from normal core-collapse supernovae) remain plausible.

PSR J1012+5307, a millisecond pulsar in orbit with a helium white dwarf (WD), has been timed with high precision for about 25 years. One of the main objectives of this long-term timing is to use the large asymmetry in gravitational binding energy between the neutron star and the WD to test gravitational theories. Such tests, however, will be eventually limited by the accuracy of the distance to the pulsar. Here, we present VLBI (very long baseline interferometry) astrometry results spanning approximately 2.5 years for PSR J1012+5307, obtained with the Very Long Baseline Array as part of the MSPSRPI project. These provide the first proper motion and absolute position for PSR J1012+5307 measured in a quasi-inertial reference frame. From the VLBI results, we measure a distance of $0.83^{+0.06}_{-0.02}$kpc (all the estimates presented in the abstract are at 68% confidence) for PSR J1012+5307, which is the most precise obtained to date. Using the new distance, we improve the uncertainty of measurements of the unmodeled contributions to orbital period decay, which, combined with three other pulsars, places new constraints on the coupling constant for dipole gravitational radiation $\kappa_D=(-1.7\pm1.7)\times 10^{-4}$ and the fractional time derivative of Newton's gravitational constant $\dot{G}/G = -1.8^{\,+5.6}_{\,-4.7}\times 10^{-13}\,{\rm yr^{-1}}$ in the local universe. As the uncertainties of the observed decays of orbital period for the four leading pulsar-WD systems become negligible in $\approx10$ years, the uncertainties for $\dot{G}/G$ and $\kappa_D$ will be improved to $\leq1.5\times10^{-13}\,{\rm yr^{-1}}$ and $\leq1.0\times10^{-4}$, respectively, predominantly limited by the distance uncertainties.

The jet opening angle and inclination of GW170817 -- the first detected binary neutron star merger -- were vital to understand its energetics, relation to short gamma-ray bursts, and refinement of the standard siren-based determination of the Hubble constant, $H_0$. These basic quantities were determined through a combination of the radio lightcurve and Very Long Baseline Interferometry (VLBI) measurements of proper motion. In this paper we discuss and quantify the prospects for the use of radio VLBI observations and observations of scintillation-induced variability to measure the source size and proper motion of merger afterglows, and thereby infer properties of the merger including inclination angle, opening angle and energetics. We show that these techniques are complementary as they probe different parts of the circum-merger density/inclination angle parameter space and different periods of the temporal evolution of the afterglow. We also find that while VLBI observations will be limited to the very closest events it will be possible to detect scintillation for a large fraction of events beyond the range of current gravitational wave detectors. Scintillation will also be detectable with next generation telescopes such as the Square Kilometre Array, 2000 antenna Deep Synoptic Array and the next generation Very Large Array, for a large fraction of events detected with third generation gravitational wave detectors. Finally, we discuss prospects for the measurement of the $H_0$ with VLBI observations of neutron star mergers and compare this technique to other standard siren methods.

We have developed a new coherent dedispersion mode to study the emission of Fast Radio Bursts that trigger the voltage capture capability of the Australian SKA Pathfinder (ASKAP) interferometer. In principle the mode can probe emission timescales down to 3 ns with full polarimetric information preserved. Enabled by the new capability, here we present a spectropolarimetric analysis of FRB 181112 detected by ASKAP, localized to a galaxy at redshift 0.47. At microsecond time resolution the burst is resolved into four narrow pulses with a rise time of just $15 \mu$s for the brightest. The pulses have a diversity of morphology, but do not show evidence for temporal broadening by turbulent plasma along the line of sight, nor is there any evidence for periodicity in their arrival times. The pulses are highly polarized (up to 95%), with the polarization position angle varying both between and within pulses. The pulses have apparent rotation measures that vary by $15\pm 2\, {\rm rad \,m^{-2}}$ and apparent dispersion measures that vary by $0.041\pm 0.004\,{\rm pc\,cm^{-3}}$. Conversion between linear and circular polarization is observed across the brightest pulse. We conclude that the FRB 181112 pulses are most consistent with being a direct manifestation of the emission process or the result of propagation through a relativistic plasma close to the source. This demonstrates that our method, which facilitates high-time-resolution polarimetric observations of FRBs, can be used to study not only burst emission processes, but also a diversity of propagation effects present on the gigaparsec paths they traverse.

Using the Very Long Baseline Array and the European Very Long Baseline Interferometry Network we have made a precise measurement of the radio parallax of the black hole X-ray binary MAXI\,J1820+070, providing a model-independent distance to the source. Our parallax measurement of ($0.348\pm0.033$) mas for MAXI J1820+070 translates to a distance of ($2.96\pm0.33$) kpc. This distance implies that the source reached ($15\pm3)\%$ of the Eddington luminosity at the peak of its outburst. Further, we use this distance to refine previous estimates of the jet inclination angle, jet velocity and the mass of the black hole in MAXI J1820+070 to be ($63\pm3)^{\circ}$, ($0.89\pm0.09)c$ and ($9.2\pm1.3) M_{\odot}$, respectively.

We report the first detection of a glitch in the radio pulsar PSR J0908$-$4913 (PSR B0906$-$49) during regular timing observations by the Molonglo Observatory Synthesis Telescope (MOST) as part of the UTMOST project.

IGR J17591$-$2342 is a new accreting millisecond X-ray pulsar (AMXP) that was recently discovered in outburst in 2018. Early observations revealed that the source's radio emission is brighter than that of any other known neutron star low-mass X-ray binary (NS-LMXB) at comparable X-ray luminosity, and assuming its likely $\gtrsim 6$ kpc distance. It is comparably radio bright to black hole LMXBs at similar X-ray luminosities. In this work, we present the results of our extensive radio and X-ray monitoring campaign of the 2018 outburst of IGR J17591$-$2342. In total we collected 10 quasi-simultaneous radio (VLA, ATCA) and X-ray (Swift-XRT) observations, which make IGR J17591$-$2342 one of the best-sampled NS-LMXBs. We use these to fit a power-law correlation index $\beta = 0.37^{+0.42}_{-0.40}$ between observed radio and X-ray luminosities ( $L_\mathrm{R}\propto L_\mathrm{X}^{\beta}$). However, our monitoring revealed a large scatter in IGR J17591$-$2342's radio luminosity (at a similar X-ray luminosity, $L_\mathrm{X} \sim 10^{36}$ erg s$^{-1}$, and spectral state), with $L_\mathrm{R} \sim 4 \times 10^{29}$ erg s$^{-1}$ during the first three reported observations, and up to a factor of 4 lower $L_\mathrm{R}$ during later radio observations. Nonetheless, the average radio luminosity of IGR J17591$-$2342 is still one of the highest among NS-LMXBs, and we discuss possible reasons for the wide range of radio luminosities observed in such systems during outburst. We found no evidence for radio pulsations from IGR J17591$-$2342 in our Green Bank Telescope observations performed shortly after the source returned to quiescence. Nonetheless, we cannot rule out that IGR J17591$-$2342 becomes a radio millisecond pulsar during quiescence.

Aql X-1 is one of the best-studied neutron star low-mass X-ray binaries. It was previously targeted using quasi-simultaneous radio and X-ray observations during at least 7 different accretion outbursts. Such observations allow us to probe the interplay between accretion inflow (X-ray) and jet outflow (radio). Thus far, these combined observations have only covered one order of magnitude in radio and X-ray luminosity range; this means that any potential radio - X-ray luminosity correlation, $L_R \propto L_X^{\beta}$, is not well constrained ($\beta \approx$ 0.4-0.9, based on various studies) or understood. Here we present quasi-simultaneous Very Large Array and Swift-XRT observations of Aql X-1's 2016 outburst, with which we probe one order of magnitude fainter in radio and X-ray luminosity compared to previous studies ($6 \times 10^{34} < L_X < 3 \times 10^{35}$ erg s$^{-1}$, i.e., the intermediate to low-luminosity regime between outburst peak and quiescence). The resulting radio non-detections indicate that Aql X-1's radio emission decays more rapidly at low X-ray luminosities than previously assumed - at least during the 2016 outburst. Assuming similar behaviour between outbursts, and combining all available data, this can be modelled as a steep $\beta=1.17^{+0.30}_{-0.21}$ power-law index or as a sharp radio cut-off at $L_X \lesssim 5 \times 10^{35}$ erg s$^{-1}$ (given our deep radio upper limits at X-ray luminosities below this value). We discuss these results in the context of other similar studies.

Present-day galaxies are surrounded by cool and enriched halo gas extending to hundreds of kiloparsecs. This halo gas is thought to be the dominant reservoir of material available to fuel future star formation, but direct constraints on its mass and physical properties have been difficult to obtain. We report the detection of a fast radio burst (FRB 181112) with arcsecond precision, which passes through the halo of a foreground galaxy. Analysis of the burst shows the halo gas has low net magnetization and turbulence. Our results imply predominantly diffuse gas in massive galactic halos, even those hosting active supermassive black holes, contrary to some previous results.

We use Very Long Baseline Interferometry to measure the proper motions of three black hole X-ray binaries (BHXBs). Using these results together with data from the literature and Gaia-DR2 to collate the best available constraints on proper motion, parallax, distance and systemic radial velocity of 16 BHXBs, we determined their three dimensional Galactocentric orbits. We extended this analysis to estimate the probability distribution for the potential kick velocity (PKV) a BHXB system could have received on formation. Constraining the kicks imparted to BHXBs provides insight into the birth mechanism of black holes (BHs). Kicks also have a significant effect on BH-BH merger rates, merger sites, and binary evolution, and can be responsible for spin-orbit misalignment in BH binary systems. $75\%$ of our systems have potential kicks $>70\,\rm{km~s^{-1}}$. This suggests that strong kicks and hence spin-orbit misalignment might be common among BHXBs, in agreement with the observed quasi-periodic X-ray variability in their power density spectra. We used a Bayesian hierarchical methodology to analyse the PKV distribution of the BHXB population, and suggest that a unimodal Gaussian model with a mean of $107\pm16\,\rm{km~s^{-1}}$ is a statistically favourable fit. Such relatively high PKVs would also reduce the number of BHs likely to be retained in globular clusters. We found no significant correlation between the BH mass and PKV, suggesting a lack of correlation between BH mass and the BH birth mechanism. Our Python code allows the estimation of the PKV for any system with sufficient observational constraints.

Fast Radio Bursts (FRBs) are brief radio emissions from distant astronomical sources. Some are known to repeat, but most are single bursts. Non-repeating FRB observations have had insufficient positional accuracy to localize them to an individual host galaxy. We report the interferometric localization of the single pulse FRB 180924 to a position 4 kpc from the center of a luminous galaxy at redshift 0.3214. The burst has not been observed to repeat. The properties of the burst and its host are markedly different from the only other accurately localized FRB source. The integrated electron column density along the line of sight closely matches models of the intergalactic medium, indicating that some FRBs are clean probes of the baryonic component of the cosmic web.

We describe the ongoing `Survey for Magnetars, Intermittent pulsars, Rotating radio transients and Fast radio bursts' (SMIRF), performed using the newly refurbished UTMOST telescope. SMIRF repeatedly sweeps the southern Galactic plane performing real-time periodicity and single-pulse searches, and is the first survey of its kind carried out with an interferometer. SMIRF is facilitated by a robotic scheduler which is capable of fully autonomous commensal operations. We report on the SMIRF observational parameters, the data analysis methods, the survey's sensitivities to pulsars, techniques to mitigate radio frequency interference and present some early survey results. UTMOST's wide field of view permits a full sweep of the Galactic plane to be performed every fortnight, two orders of magnitude faster than previous surveys. In the six months of operations from January to June 2018, we have performed $\sim 10$ sweeps of the Galactic plane with SMIRF. Notable blind re-detections include the magnetar PSR J1622$-$4950, the RRAT PSR J0941$-$3942 and the eclipsing pulsar PSR J1748$-$2446A. We also report the discovery of a new pulsar, PSR J1705$-$54. Our follow-up of this pulsar with the UTMOST and Parkes telescopes at an average flux limit of $\leq 20$ mJy and $\leq 0.16$ mJy respectively, categorizes this as an intermittent pulsar with a high nulling fraction of $< 0.002$

We detail a new fast radio burst (FRB) survey with the Molonglo Radio Telescope, in which six FRBs were detected between June 2017 and December 2018. By using a real-time FRB detection system, we captured raw voltages for five of the six events, which allowed for coherent dedispersion and very high time resolution (10.24 $\mu$s) studies of the bursts. Five of the FRBs show temporal broadening consistent with interstellar and/or intergalactic scattering, with scattering timescales ranging from 0.16 to 29.1 ms. One burst, FRB181017, shows remarkable temporal structure, with 3 peaks each separated by 1 ms. We searched for phase-coherence between the leading and trailing peaks and found none, ruling out lensing scenarios. Based on this survey, we calculate an all-sky rate at 843 MHz of $98^{+59}_{-39}$ events sky$^{-1}$ day$^{-1}$ to a fluence limit of 8 Jy-ms: a factor of 7 below the rates estimated from the Parkes and ASKAP telescopes at 1.4 GHz assuming the ASKAP-derived spectral index $\alpha=-1.6$ ($F_{\nu}\propto\nu^{\alpha}$). Our results suggest that FRB spectra may turn over below 1 GHz. Optical, radio and X-ray followup has been made for most of the reported bursts, with no associated transients found. No repeat bursts were found in the survey.

We present the results of a 7-12 GHz phased-array study of the Galactic center magnetar J1745-2900 with the Karl G. Jansky Very Large Array (VLA). Using data from two 6.5 hour observations from September 2014, we find that the average profile is comprised of several distinct components at these epochs and is stable over $\sim$day timescales and $\sim$GHz frequencies. Comparison with additional phased VLA data at 8.7 GHz shows significant profile changes on longer timescales. The average profile at 7-12 GHz is dominated by the jitter of relatively narrow pulses. The pulses in each of the four main profile components seen in September 2014 are uncorrelated in phase and amplitude, though there is a small but significant correlation in the occurrence of pulses in two of the profile components. Using the brightest pulses, we measure the dispersion and scattering parameters of J1745-2900. A joint fit of 38 pulses gives a 10 GHz pulse broadening time of $\tau_{\rm sc, 10} = 0.09 \pm 0.03~\rm ms$ and a dispersion measure of ${\rm DM} = 1760^{+2.4}_{-1.3}~{\rm pc~cm}^{-3}$. Both of these results are consistent with previous measurements, which suggests that the scattering and dispersion measure of J1745-2900 may be stable on timescales of several years.

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.

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.

Gravitational lensing is a powerful tool for quantifying the mass content and distribution in distant galaxies. By using milliarcsecond angular resolution observations of radio-loud gravitationally lensed sources it is also possible to detect and quantify small deviations from a smooth mass density distribution, which can be due to low mass substructures in the lensing galaxy. We present high-resolution global VLBI observations of the gravitationally lensed radio source MG J0751+2716 (at z = 3.2), that shows evidence of both compact and extended structure (core-jet morphology) across several gravitational arcs. These data provide a wealth of observational constraints that are used to determine the inner (baryonic and dark matter) mass profile of a group of galaxies and also investigate the smoothness of the dark matter distribution on mas-scales, which is sensitive to possible structures of $10^{6-7}$ M$_{\odot}$ within the lensing halo or along the line-of-sight. Our lens modelling finds evidence for astrometric anomalies in this system, which suggest presence of extra mass structure in the lens model. To date this kind of detailed studies of gravitational lensing systems like MG J0751+2716 has been limited by the currently small sample of radio-loud gravitational lenses. In this context, we also present a new pilot gravitational lens search in the VLBI survey mJIVE-20, in perspective of future surveys with the next generation of radio interferometers.

We present a novel pilot search for gravitational lenses in the mJIVE-20 survey, which observed $24\,903$ radio sources selected from FIRST with the VLBA at an angular resolution of 5 mas. We have taken the visibility data for an initial $3\,640$ sources that were detected by the mJIVE-20 observations and re-mapped them to make wide-field images, selecting fourteen sources that had multiple components separated by $\geq100$ mas, with a flux-ratio of $\leq15$:$1$ and a surface brightness consistent with gravitational lensing. Two of these candidates are re-discoveries of gravitational lenses found as part of CLASS. The remaining twelve candidates were then re-observed at 1.4 GHz and then simultaneously at 4.1 and 7.1 GHz with the VLBA to measure the spectral index and surface brightness of the individual components as a function of frequency. Ten were rejected as core-jet or core-hotspot(s) systems, with surface brightness distributions and/or spectral indices inconsistent with gravitational lensing, and one was rejected after lens modelling demonstrated that the candidate lensed images failed the parity test. The final lens candidate has an image configuration that is consistent with a simple lens mass model, although further observations are required to confirm the lensing nature. Given the two confirmed gravitational lenses in the mJIVE-20 sample, we find a robust lensing-rate of $1$:($318\pm225$) for a statistical sample of 635 radio sources detected on mas-scales, which is consistent with that found for CLASS.

The compact radio source Sagittarius~A$^*$ (Sgr~A$^*$)in the Galactic Center is the primary supermassive black hole candidate. General relativistic magnetohydrodynamical (GRMHD) simulations of the accretion flow around Sgr\,A$^*$ predict the presence of sub-structure at observing wavelengths of $\sim 3$\,mm and below (frequencies of 86\,GHz and above). For very long baseline interferometry (VLBI) observations of Sgr\,A$^*$ at this frequency the blurring effect of interstellar scattering becomes subdominant, and arrays such as the High Sensitivity Array (HSA) and the global mm-VLBI Array (GMVA) are now capable of resolving potential sub-structure in the source. Such investigations improve our understanding of the emission geometry of the mm-wave emission of Sgr\,A$^*$, which is crucial for constraining theoretical models and for providing a background to interpret 1\,mm VLBI data from the Event Horizon Telescope (EHT). We performed high-sensitivity very long baseline interferometry (VLBI) observations of Sgr\,A$^*$ at 3\,mm using the Very Long Baseline Array (VLBA) and the Large Millimeter Telescope (LMT) in Mexico on two consecutive days in May 2015, with the second epoch including the Green Bank Telescope (GBT). We find an overall source geometry that matches previous findings very closely, showing a deviation in fitted model parameters less than 3\% over a time scale of weeks and suggesting a highly stable global source geometry over time. The reported sub-structure in the 3\,mm emission of Sgr\,A$^*$ is consistent with theoretical expectations of refractive noise on long baselines. However, comparing our findings with recent results from 1\,mm and 7\,mm VLBI observations, which also show evidence for east-west asymmetry, an intrinsic origin cannot be excluded. Confirmation of persistent intrinsic substructure will require further VLBI observations spread out over multiple epochs.

The binary neutron star merger GW170817 was accompanied by radiation across the electromagnetic spectrum and localized to the galaxy NGC 4993 at a distance of 41+/-3 Mpc. The radio and X-ray afterglows of GW170817 exhibited delayed onset, a gradual rise in the emission with time as t^0.8, a peak at about 150 days post-merger, followed by a relatively rapid decline. To date, various models have been proposed to explain the afterglow emission, including a choked-jet cocoon and a successful-jet cocoon (a.k.a. structured jet). However, the observational data have remained inconclusive as to whether GW170817 launched a successful relativistic jet. Here we show, through Very Long Baseline Interferometry, that the compact radio source associated with GW170817 exhibits superluminal motion between two epochs at 75 and 230 days post-merger. This measurement breaks the degeneracy between the models and indicates that, while the early-time radio emission was powered by a wider-angle outflow (cocoon), the late-time emission was most likely dominated by an energetic and narrowly-collimated jet, with an opening angle of <5 degrees, and observed from a viewing angle of about 20 degrees. The imaging of a collimated relativistic outflow emerging from GW170817 adds substantial weight to the growing evidence linking binary neutron star mergers and short gamma-ray bursts.

The second data release from the Gaia mission (Gaia DR2) includes, among its billion entries, as- trometric parameters for binary companions to a number of known pulsars, including white dwarf companions to millisecond pulsars and the non-degenerate components of so-called "black widow" and "redback" systems. We find 22 such counterparts in DR2, of which 12 have statistically sig- nificant measurements of parallax. These DR2 optical proper motions and parallaxes provide new measurements of the distances and transverse velocities of the associated pulsars. For the most part, the results agree with existing radio interferometric and pulsar timing-based astrometry, as well as other distance estimates based on photometry or associations, and for some pulsars they provide the best known distance and velocity estimates. In particular, two of these pulsars have no previous dis- tance measurement: PSR J1227$-$4853, for which Gaia measures a parallax of $0.62\pm 0.16$ mas, and PSR J1431$-$4715, with a Gaia parallax of $0.64\pm 0.16$ mas. Using the Gaia distance measurements, we find that dispersion measure-based distance estimates calculated using the NE2001 and YMW16 Galactic electron density models are on average slightly underestimated, which may be a selection effect due to the over-representation of pulsars at high Galactic latitudes in the present Gaia sample. While the Gaia DR2 results do not quite match the precision that can be achieved by dedicated pulsar timing or radio interferometry, taken together they constitute a small but important improvement to the pulsar distance scale, and the subset of millisecond pulsars with distances measured by Gaia may help improve the sensitivity of pulsar timing arrays to nanohertz gravitational waves.

(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).

We report the detection of a glitch event in the pulsar J1709$-$4429 (also known as B1706$-$44) during regular monitoring observations with the Molonglo Observatory Synthesis Telescope (UTMOST). The glitch was found during timing operations, in which we regularly observe over 400 pulsars with up to daily cadence, while commensally searching for Rotating Radio Transients, pulsars, and FRBs. With a fractional size of $\Delta\nu/\nu \approx 52.4 \times10^{-9}$, the glitch reported here is by far the smallest known for this pulsar, attesting to the efficacy of glitch searches with high cadence using UTMOST.

Einstein's theory of gravity, general relativity, has passed stringent tests in laboratories, elsewhere in the Solar Sytem, and in pulsar binaries. Nevertheless it is known to be incompatible with quantum mechanics and must differ from the true behaviour of matter in strong fields and at small spatial scales. A key aspect of general relativity to test is the strong equivalence principle (SEP), which states that all freely falling objects, regardless of how strong their gravity, experience the same acceleration in the same gravitational field. Essentially all alternatives to general relativity violate this principle at some level. Previous direct tests of the SEP are limited by the weak gravity of the bodies in the Earth-Moon-Sun system or by the weak gravitational pull of the Galaxy on pulsar-white dwarf binaries. PSR~J0337+1715 is a hierarchical stellar triple system, where the inner binary consists of a millisecond radio pulsar in a $1.6$-day orbit with a white dwarf. This inner binary is in a $327$-day orbit with another white dwarf. In this system, the pulsar and the inner companion fall toward the outer companion with an acceleration about $10^8$ times greater than that produced by falling in the Galactic potential, and the pulsar's gravitational binding energy is roughly $10\%$ of its mass. Here we report that in spite of the pulsar's strong gravity, the accelerations experienced by it and the inner white dwarf differ by a fraction of no more than $2.6\times 10^{-6}$ ($95\%$ confidence level). We can roughly compare this to other SEP tests by using the strong-field Nordtvedt parameter $\hat\eta_N$. Our limit on $\hat\eta_N$ is a factor of ten smaller than that obtained from (weak-field) Solar-System SEP tests and a factor of almost a thousand smaller than that obtained from other strong-field SEP tests.

The pulsar/massive star binary system PSR B1259-63 / LS 2883 is one of the best-studied gamma-ray binaries, a class of systems whose bright gamma-ray flaring can provide important insights into high-energy physics. Using the Australian Long Baseline Array we have conducted very long baseline interferometric observations of PSR B1259-63 over 4.4 years, fully sampling the 3.4-year orbital period. From our measured parallax of $0.38\pm0.05$ mas we use a Bayesian approach to infer a distance of $2.6^{+0.4}_{-0.3}$ kpc. We find that the binary orbit is viewed at an angle of $154\pm3$ degrees to the line of sight, implying that the pulsar moves clockwise around its orbit as viewed on the sky. Taking our findings together with previous results from pulsar timing observations, all seven orbital elements for the system are now fully determined. We use our measurement of the inclination angle to constrain the mass of the stellar companion to lie in the range 15-31$M_{\odot}$. Our measured distance and proper motion are consistent with the system having originated in the Cen OB1 association and receiving a modest natal kick, causing it to have moved $\sim$8 pc from its birthplace over the past $\sim3\times10^5$ years. The orientation of the orbit on the plane of the sky matches the direction of motion of the X-ray synchrotron-emitting knot observed by the Chandra X-ray Observatory to be moving away from the system.

The Hubble constant ($H_0$) measures the current expansion rate of the Universe, and plays a fundamental role in cosmology. Tremendous effort has been dedicated over the past decades to measure $H_0$. Notably, Planck cosmic microwave background (CMB) and the local Cepheid-supernovae distance ladder measurements determine $H_0$ with a precision of $\sim 1\%$ and $\sim 2\%$ respectively. A $3$-$\sigma$ level of discrepancy exists between the two measurements, for reasons that have yet to be understood. Gravitational wave (GW) sources accompanied by electromagnetic (EM) counterparts offer a completely independent standard siren (the GW analogue of an astronomical standard candle) measurement of $H_0$, as demonstrated following the discovery of the neutron star merger, GW170817. This measurement does not assume a cosmological model and is independent of a cosmic distance ladder. The first joint analysis of the GW signal from GW170817 and its EM localization led to a measurement of $H_0=74^{+16}_{-8}$ km/s/Mpc (median and symmetric $68\%$ credible interval). In this analysis, the degeneracy in the GW signal between the source distance and the weakly constrained viewing angle dominated the $H_0$ measurement uncertainty. Recently, Mooley et al. (2018) obtained tight constraints on the viewing angle using high angular resolution imaging of the radio counterpart of GW170817. Here we obtain a significantly improved measurement $H_0=68.9^{+4.7}_{-4.6}$ km/s/Mpc by using these new radio observations, combined with the previous GW and EM data. We estimate that 15 more localized GW170817-like events (comparable signal-to-noise ratio, favorable orientation), having radio images and light curve data, will potentially bring resolution to the tension between the Planck and Cepheid-supernova measurements, as compared to 50-100 GW events without such data.

Using the Very Long Baseline Array, we have made astrometric observations of the binary pulsar B1913+16 spanning an 18 month period in 2014 - 2015. From these observations we make the first determination of the annual geometric parallax of B1913+16, measuring $\pi= 0.24^{+0.06}_{-0.08}$ milliarcseconds (68% confidence interval). The inferred parallax probability distribution differs significantly from a Gaussian. Using our parallax measurement and prior information on the spatial and luminosity distributions of the millisecond pulsar population, we infer a distance of $d=4.1^{+2.0}_{-0.7}$ kpc, which is significantly closer than the $9.8 \pm 3.1$ kpc suggested by the pulsar's dispersion measure and analyses of the ionized interstellar medium. While the relatively low significance of the parallax detection ($\sim$3$\sigma$) currently precludes an improved test of general relativity using the orbital decay of PSR B1913+16, ongoing observations with improved control of systematic astrometric errors could reach the 10% distance uncertainty required for this goal. The proper motion measured by our VLBI astrometry differs substantially from that obtained by pulsar timing, a discrepancy that has also been found between the proper motion measurements made by interferometers and pulsar timing for some other pulsars and which we speculate is the result of timing noise or dispersion measure variations in the timing data set. Our parallax and proper motion measurements yield a transverse velocityof $15^{+8}_{-4}$ km s$^{-1}$ in the solar reference frame. Analysis incorporating galactic rotation and solar motion finds the space velocity of the pulsar relative to its standard of rest has a component 22$^{+7}_{-3}$ km s$^{-1}$ perpendicular to the galactic plane and components of order 100 km s$^{-1}$ parallel to the galactic plane.

We present high spatial resolution observations of the continuum emission from the young multiple star system UZ Tau at frequencies from 6 to 340 GHz. To quantify the spatial variation of dust emission in the UZ Tau E circumbinary disk, the observed interferometric visibilities are modeled with a simple parametric prescription for the radial surface brightnesses at each frequency. We find evidence that the spectrum steepens with radius in the disk, manifested as a positive correlation between the observing frequency and the radius that encircles a fixed fraction of the emission ($R_{eff} \propto \nu^{0.34 \pm 0.08}$). The origins of this size--frequency relation are explored in the context of a theoretical framework for the growth and migration of disk solids. While that framework can reproduce a similar size--frequency relation, it predicts a steeper spectrum than is observed. Moreover, it comes closest to matching the data only on timescales much shorter ($\le 1$ Myr) than the putative UZ Tau age (~2-3 Myr). These discrepancies are the direct consequences of the rapid radial drift rates predicted by models of dust evolution in a smooth gas disk. One way to mitigate that efficiency problem is to invoke small-scale gas pressure modulations that locally concentrate drifting solids. If such particle traps reach high continuum optical depths at 30-340 GHz with a ~30-60% filling fraction in the inner disk ($r \lesssim20$ au), they can also explain the observed spatial gradient in the UZ Tau E disk spectrum.

We report a new Fast Radio Burst (FRB) discovered in real-time as part of the UTMOST project at the Molonglo Observatory Synthesis Radio Telescope (MOST). FRB170827 is the first detected with our low-latency ($< 24$ s), machine-learning-based FRB detection system. The FRB discovery was accompanied by the capture of voltage data at the native time and frequency resolution of the observing system, enabling coherent dedispersion and detailed off-line analysis, which have unveiled fine temporal and frequency structure. The dispersion measure (DM) of 176.80 $\pm$ 0.04 pc cm$^{-3}$, is the lowest of the FRB population. The Milky Way contribution along the line of sight is $\sim$ 40 pc cm$^{-3}$, leaving an excess DM of $\sim$ 145 pc cm$^{-3}$. The FRB has a fluence $>$ 20 $\pm$ 7 Jy ms, and is narrow, with a width of $\sim$ 400 $\mu$s at 10$\%$ of its maximum amplitude. However, the burst shows three temporal components, the narrowest of which is $\sim$ 30 $\mu$s, and a scattering timescale of $4.1 \pm 2.7$ $\mu$s. The FRB shows spectral modulations on frequency scales of 1.5 MHz and 0.1 MHz. Both are prominent in the dynamic spectrum, which shows a very bright region of emission between 841 and 843 MHz, and weaker, patchy emission across the entire band. We show the fine spectral structure could arise in the FRB host galaxy, or its immediate vicinity.

We present very long baseline interferometry (VLBI) observations of 179 radio sources in the COSMOS field with extremely high sensitivity using the Green Bank Telescope (GBT) together with the Very Long Baseline Array (VLBA) (VLBA+GBT) at 1.4 GHz, to explore the faint radio population in the flux density regime of tens of $\mu$Jy. Here, the identification of active galactic nuclei (AGN) is based on the VLBI detection of the source, i.e., it is independent of X-ray or infrared properties. The milli-arcsecond resolution provided by the VLBI technique implies that the detected sources must be compact and have large brightness temperatures, and therefore they are most likely AGN (when the host galaxy is located at z$\geq$0.1). On the other hand, this technique allows us to only positively identify when a radio-active AGN is present, i.e., we cannot affirm that there is no AGN when the source is not detected. For this reason, the number of identified AGN using VLBI should be always treated as a lower limit. We present a catalogue containing the 35 radio sources detected with the VLBA+GBT, 10 of which were not previously detected using only the VLBA. We have constructed the radio source counts at 1.4 GHz using the samples of the VLBA and VLBA+GBT detected sources of the COSMOS field to determine a lower limit for the AGN contribution to the faint radio source population. We found an AGN contribution of >40-75% at flux density levels between 150 $\mu$Jy and 1 mJy. This flux density range is characterised by the upturn of the Euclidean-normalised radio source counts, which implies a contribution of a new population. This result supports the idea that the sub-mJy radio population is composed of a significant fraction of radio-emitting AGN, rather than solely by star-forming galaxies, in agreement with previous studies.

Type II radio bursts are evidence of shocks in the solar atmosphere and inner heliosphere that emit radio waves ranging from sub-meter to kilometer lengths. These shocks may be associated with CMEs and reach speeds higher than the local magnetosonic speed. Radio imaging of decameter wavelengths (20-90 MHz) is now possible with LOFAR, opening a new radio window in which to study coronal shocks that leave the inner solar corona and enter the interplanetary medium and to understand their association with CMEs. To this end, we study a coronal shock associated with a CME and type II radio burst to determine the locations at which the radio emission is generated, and we investigate the origin of the band-splitting phenomenon.

Previous optical and radio observations of the binary millisecond pulsar PSR J1640+2224 have come to inconsistent conclusions about the identity of its companion, with some observations suggesting the companion is a low-mass helium-core (He-core) white dwarf (WD), while others indicate it is most likely a high-mass carbon-oxygen (CO) WD. Binary evolution models predict PSR J1640+2224 most likely formed in a low-mass X-ray binary (LMXB) based on the pulsar's short spin period and long-period, low-eccentricity orbit, in which case its companion should be a He-core WD with mass about $0.35 - 0.39 \, M_\odot$, depending on metallicity. If it is instead a CO WD, that would suggest the system has an unusual formation history. In this paper we present the first astrometric parallax measurement for this system from observations made with the Very Long Baseline Array (VLBA), from which we determine the distance to be $1520^{+170}_{-150}\,\mathrm{pc}$. We use this distance and a reanalysis of archival optical observations originally taken in 1995 with the Wide Field Planetary Camera 2 (WFPC2) on the Hubble Space Telescope (HST) in order to measure the WD's mass. We also incorporate improvements in calibration, extinction model, and WD cooling models. We find that the existing observations are not sufficient to tightly constrain the companion mass, but we conclude the WD mass is $>0.4\,M_\odot$ with $>90\%$ confidence. The limiting factor in our analysis is the low signal-to-noise ratio of the original HST observations.

We report on the first sub-arcsecond (0.44 $\times$ 0.41 arcsec$\rm ^2$) angular resolution image at 150 MHz of the A-nucleus in the Luminous Infrared Galaxy Arp$\,$299, from International Low Frequency Array (LOFAR) Telescope observations. The most remarkable finding is that of an intriguing two-sided, filamentary structure emanating from A-nucleus, which we interpret as an outflow that extends up to at least 14 arcseconds from the A-nucleus in the N-S direction ($\approx$ 5 kpc deprojected size) and accounts for almost 40% of the extended emission of the entire galaxy system. We also discuss HST/NICMOS [FeII] 1.64 $\rm \mu m$ and H$\rm_2$ 2.12 $\rm \mu m$ images of Arp$\,$299-A, which show similar features to those unveiled by our 150 MHz LOFAR observations, thus giving string morphological support for the outflow scenario. Finally, we discuss unpublished NaI D spectra that confirm the outflow nature of this structure. From energetic arguments, we rule out the low-luminosity active galactic nucleus in Arp$\,$299-A as a driver for the outflow. On the contrary, the powerful, compact starburst in the central regions of Arp$\,$299-A provides plenty of mechanical energy to sustain an outflow, and we conclude that the intense supernova (SN) activity in the nuclear region of Arp299-A is driving the observed outflow. We estimate that the starburst wind can support a mass-outflow rate in the range (11-63) $\rm M_{\odot} yr^{-1}$ at speeds of up to (370 - 890) $\rm km \, s^{-1}$, and is relatively young, with an estimated kinematic age of (3 - 7) Myr. Those results open an avenue to the use of low-frequency (150 MHz), sub-arcsecond imaging with LOFAR to detect outflows in the central regions of local luminous infrared galaxies.

We have observed a complex and continuous change in the integrated pulse profile of PSR B2217+47, manifested as additional components trailing the main peak. These transient components are detected over 6 years at $150$ MHz using the LOw Frequency ARray (LOFAR), but they are not seen in contemporaneous Lovell observations at $1.5$ GHz. We argue that propagation effects in the ionized interstellar medium (IISM) are the most likely cause. The putative structures in the IISM causing the profile variation are roughly half-way between the pulsar and the Earth and have transverse radii $R \sim 30$ AU. We consider different models for the structures. Under the assumption of spherical symmetry, their implied average electron density is $\overline{n}_e \sim 100$ cm$^{-3}$. Since PSR B2217+47 is more than an order of magnitude brighter than the average pulsar population visible to LOFAR, similar profile variations would not have been identified in most pulsars, suggesting that subtle profile variations in low-frequency profiles might be more common than we have observed to date. Systematic studies of these variations at low frequencies can provide a new tool to investigate the proprieties of the IISM and the limits to the precision of pulsar timing.

Newton's gravitational constant $G$ may vary with time at an extremely low level. The time variability of $G$ will affect the orbital motion of a millisecond pulsar in a binary system and cause a tiny difference between the orbital period-dependent measurement of the kinematic distance and the direct measurement of the annual parallax distance. PSR J0437$-$4715 is the nearest millisecond pulsar and the brightest at radio. To explore the feasibility of achieving a parallax distance accuracy of one light-year, comparable to the recent timing result, with the technique of differential astrometry, we searched for compact radio sources quite close to PSR J0437$-$4715. Using existing data from the Very Large Array and the Australia Telescope Compact Array, we detected two sources with flat spectra, relatively stable flux densities of 0.9 and 1.0 mJy at 8.4 GHz and separations of 13 and 45 arcsec. With a network consisting of the Long Baseline Array and the Kunming 40-m radio telescope, we found that both sources have a point-like structure and a brightness temperature of $\geq$10$^7$ K. According to these radio inputs and the absence of counterparts in the other bands, we argue that they are most likely the compact radio cores of extragalactic active galactic nuclei rather than Galactic radio stars. The finding of these two radio active galactic nuclei will enable us to achieve a sub-pc distance accuracy with the in-beam phase-referencing very-long-baseline interferometric observations and provide one of the most stringent constraints on the time variability of $G$ in the near future.

We present Karl G. Jansky Very Large Array radio frequency observations of the new accreting millisecond X-ray pulsar (AMXP), IGR J16597$-$3704, located in the globular cluster NGC 6256. With these data, we detect a radio counterpart to IGR J16597$-$3704, and determine an improved source position. Pairing our radio observations with quasi-simultaneous Swift/XRT X-ray observations, we place IGR J16597$-$3704 on the radio -- X-ray luminosity plane, where we find that IGR J16597$-$3704 is one of the more radio-quiet neutron star low-mass X-ray binaries known to date. We discuss the mechanisms that may govern radio luminosity (and in turn jet production and evolution) in AMXPs. Further, we use our derived radio position to search for a counterpart in archival Hubble Space Telescope and Chandra X-ray Observatory data, and estimate an upper limit on the X-ray luminosity of IGR J16597$-$3704 during quiescence.

GW170817 is the first gravitational wave detection of a binary neutron star merger. It was accompanied by radiation across the electromagnetic spectrum and localized to the galaxy NGC 4993 at a distance of 40 Mpc. It has been proposed that the observed gamma-ray, X-ray and radio emission is due to an ultra-relativistic jet launched during the merger, directed away from our line of sight. The presence of such a jet is predicted from models positing neutron star mergers as the central engines driving short-hard gamma-ray bursts (SGRBs). Here we show that the radio light curve of GW170817 has no direct signature of an off-axis jet afterglow. While we cannot rule out the existence of a jet pointing elsewhere, the observed gamma-rays could not have originated from such a jet. Instead, the radio data requires a mildly relativistic wide-angle outflow moving towards us. This outflow could be the high velocity tail of the neutron-rich material dynamically ejected during the merger or a cocoon of material that breaks out when a jet transfers its energy to the dynamical ejecta. The cocoon scenario can explain the radio light curve of GW170817 as well as the gamma-rays and X-rays (possibly also ultraviolet and optical emission), and hence is the model most consistent with the observational data. Cocoons may be a ubiquitous phenomenon produced in neutron star mergers, giving rise to a heretofore unidentified population of radio, ultraviolet, X-ray and gamma-ray transients in the local universe.

The discovery of the first electromagnetic counterpart to a gravitational wave signal has generated follow-up observations by over 50 facilities world-wide, ushering in the new era of multi-messenger astronomy. In this paper, we present follow-up observations of the gravitational wave event GW170817 and its electromagnetic counterpart SSS17a/DLT17ck (IAU label AT2017gfo) by 14 Australian telescopes and partner observatories as part of Australian-based and Australian-led research programs. We report early- to late-time multi-wavelength observations, including optical imaging and spectroscopy, mid-infrared imaging, radio imaging, and searches for fast radio bursts. Our optical spectra reveal that the transient source afterglow cooled from approximately 6400K to 2100K over a 7-day period and produced no significant optical emission lines. The spectral profiles, cooling rate, and photometric light curves are consistent with the expected outburst and subsequent processes of a binary neutron star merger. Star formation in the host galaxy probably ceased at least a Gyr ago, although there is evidence for a galaxy merger. Binary pulsars with short (100 Myr) decay times are therefore unlikely progenitors, but pulsars like PSR B1534+12 with its 2.7 Gyr coalescence time could produce such a merger. The displacement (about 2.2 kpc) of the binary star system from the centre of the main galaxy is not unusual for stars in the host galaxy or stars originating in the merging galaxy, and therefore any constraints on the kick velocity imparted to the progenitor are poor.

We present coordinated Chandra X-ray Observatory and Karl G. Jansky Very Large Array observations of the transitional millisecond pulsar PSR J1023+0038 in its low-luminosity accreting state. The unprecedented five hours of strictly simultaneous X-ray and radio continuum coverage for the first time unambiguously show a highly reproducible, anti-correlated variability pattern. The characteristic switches from the X-ray high mode into a low mode are always accompanied by a radio brightening with duration that closely matches the X-ray low mode interval. This behavior cannot be explained by a canonical inflow/outflow accretion model where the radiated emission and the jet luminosity are powered by, and positively correlated with, the available accretion energy. We interpret this phenomenology as alternating episodes of low-level accretion onto the neutron star during the X-ray high mode that are interrupted by rapid ejections of plasma by the active rotation-powered pulsar, possibly initiated by a reconfiguration of the pulsar magnetosphere, that cause a transition to a less luminous X-ray mode. The observed anti-correlation between radio and X-ray luminosity has an additional consequence: transitional MSPs can make excursions into a region of the radio/X-ray luminosity plane previously thought to be occupied solely by black hole X-ray binary sources. This complicates the use of this luminosity relation to identify candidate black holes, suggesting the need for additional discriminants when attempting to establish the true nature of the accretor.

We study the faint radio population using wide-field very long baseline interferometry (VLBI) observations of 2865 known radio sources in the Cosmic Evolution Survey (COSMOS) field. The main objective of the project was to determine where active galactic nuclei (AGN) are present. The combination of number of sources, sensitivity, angular resolution and area covered by this project are unprecedented. We have detected 468 radio sources, expected to be AGNs, with the Very Long Baseline Array (VLBA) at 1.4 GHz. This is, to date, the largest sample assembled of VLBI detected sources in the sub-mJy regime. The input sample was taken from previous observations with the Very Large Array (VLA). We present the catalogue with additional multiwavelength information. We find a detection fraction of 20%, considering only those sources from the input catalogue which were in principle detectable with the VLBA (2361). As a function of redshift, we see no evolution of the detection fraction over the redshift range 0.5<z<3. In addition, we find that faint radio sources typically have a greater fraction of their radio luminosity in a compact core: ~70% of the sub-mJy sources detected with the VLBA have more than half of their total radio luminosity in a VLBI-scale component, whereas this is true for only ~30% of the sources that are brighter than 10 mJy. This suggests that fainter radio sources differ intrinsically from brighter ones. Across our entire sample, we find the predominant morphological classification of the host galaxies of the VLBA detected sources to be early type (57%), although this varies with redshift and at z>1.5 we find that spiral galaxies become the most prevalent (48%). We demonstrate that wide-field VLBI observations, together with new calibration methods such as multi-source self-calibration and mosaicing, result in information which is difficult or impossible to obtain otherwise.

We use VLBA+VLA observations to measure the sizes of the scatter-broadened images of 6 of the most heavily scattered known pulsars: 3 within the Galactic Centre (GC) and 3 elsewhere in the inner Galactic plane. By combining the measured sizes with temporal pulse broadening data from the literature and using the thin-screen approximation, we locate the scattering medium along the line of sight to these 6 pulsars. At least two scattering screens are needed to explain the observations of the GC sample. We show that the screen inferred by previous observations of SGR J1745-2900 and Sgr A*, which must be located far from the GC, falls off in strength on scales < 0.2 degree. A second scattering component closer to (< 2 kpc) or even (tentatively) within (< 700 pc) the GC produces most or all of the temporal broadening observed in the other GC pulsars. Outside the GC, the scattering locations for all three pulsars are ~2 kpc from Earth, consistent with the distance of the Carina-Sagittarius or Scutum spiral arm. For each object the 3D scattering origin coincides with a known HII region (and in one case also a supernova remnant), suggesting that such objects preferentially cause the intense interstellar scattering seen towards the Galactic plane. We show that the HII regions should contribute > 25% of the total dispersion measure (DM) towards these pulsars, and calculate reduced DM distances. Those distances for other pulsars lying behind HII regions may be similarly overestimated.

We present quasi-simultaneous radio (VLA) and X-ray ($Swift$) observations of the neutron star low-mass X-ray binary (NS-LMXB) 1RXS J180408.9$-$342058 (J1804) during its 2015 outburst. We found that the radio jet of J1804 was bright ($232 \pm 4 \mu$Jy at $10$ GHz) during the initial hard X-ray state, before being quenched by more than an order of magnitude during the soft X-ray state ($19 \pm 4 \mu$Jy). The source then was undetected in radio (< $13 \mu$Jy) as it faded to quiescence. In NS-LMXBs, possible jet quenching has been observed in only three sources and the J1804 jet quenching we show here is the deepest and clearest example to date. Radio observations when the source was fading towards quiescence ($L_X = 10^{34-35}$ erg s$^{-1}$) show that J1804 must follow a steep track in the radio/X-ray luminosity plane with $\beta > 0.7$ (where $L_R \propto L_X^{\beta}$). Few other sources have been studied in this faint regime, but a steep track is inconsistent with the suggested behaviour for the recently identified class of transitional millisecond pulsars. J1804 also shows fainter radio emission at $L_X < 10^{35}$ erg s$^{-1}$ than what is typically observed for accreting millisecond pulsars. This suggests that J1804 is likely not an accreting X-ray or transitional millisecond pulsar.

In outburst, neutron star X-ray binaries produce less powerful jets than black holes at a given X-ray luminosity. This has made them more difficult to study as they fade towards quiescence. To explore whether neutron stars power jets at low accretion rates ($L_{\rm X} \lesssim 10^{36}$ erg s$^{-1}$), we investigate the radio and X-ray properties of three accreting millisecond X-ray pulsars (IGR J17511-3057, SAX J1808.4-3658 and IGR J00291+5934) during their outbursts in 2015, and of the non-pulsing neutron star Cen X-4 in quiescence (2015) and in outburst (1979). We did not detect the radio counterpart of IGR J17511-3057 in outburst or of Cen X-4 in quiescence, but did detect IGR J00291+5934 and SAX J1808.4-3658, showing that at least some neutron stars launch jets at low accretion rates. While the radio and X-ray emission in IGR J00291+5934 seem to be tightly correlated, the relationship in SAX J1808.4-3658 is more complicated. We find that SAX J1808.4-3658 produces jets during the reflaring tail, and we explore a toy model to ascertain whether the radio emission could be attributed to the onset of a strong propeller. The lack of a universal radio/X-ray correlation, with different behaviours in different neutron star systems (with various radio/X-ray correlations; some being radio faint and others not), points at distinct disk-jet interactions in individual sources, while always being fainter in the radio band than black holes at the same X-ray luminosity.

Recent studies have shown possible connections between highly magnetized neutron stars ("magnetars"), whose X-ray emission is too bright to be powered by rotational energy, and ordinary radio pulsars. In addition to the magnetar SGR J1745-2900, one of the radio pulsars in the Galactic centre (GC) region, PSR J1746-2850, had timing properties implying a large magnetic field strength and young age, as well as a flat spectrum. All characteristics are similar to those of rare, transient, radio-loud magnetars. Using several deep non-detections from the literature and two new detections, we show that this pulsar is also transient in the radio. Both the flat spectrum and large amplitude variability are inconsistent with the light curves and spectral indices of 3 radio pulsars with high magnetic field strengths. We further use frequent, deep archival imaging observations of the GC in the past 15 years to rule out a possible X-ray outburst with a luminosity exceeding the rotational spin down rate. This source, either a transient magnetar without any detected X-ray counterpart or a young, strongly magnetized radio pulsar producing magnetar-like radio emission, further blurs the line between the two categories. We discuss the implications of this object for the radio emission mechanism in magnetars and for star and compact object formation in the GC.

We have monitored the transient neutron star low-mass X-ray binary 1RXS J180408.9-342058 in quiescence after its ~4.5 month outburst in 2015. The source has been observed using Swift and XMM-Newton. Its X-ray spectra were dominated by a thermal component. The thermal evolution showed a gradual X-ray luminosity decay from ~18 x 10^32 to ~4 x 10^32 (D/5.8 kpc)^2 erg s^{-1} between ~8 and ~379 d in quiescence, and the inferred neutron star surface temperature (for an observer at infinity; using a neutron star atmosphere model) decreased from ~100 to ~71 eV. This can be interpreted as cooling of an accretion-heated neutron star crust. Modelling the observed temperature curve (using NSCOOL) indicated that the source required ~1.9 MeV per accreted nucleon of shallow heating in addition to the standard deep crustal heating to explain its thermal evolution. Alternatively, the decay could also be modelled without the presence of deep crustal heating, only having a shallow heat source (again ~1.9 MeV per accreted nucleon was required). However, the XMM-Newton data statistically required an additional power-law component. This component contributed ~30 per cent of the total unabsorbed flux in 0.5-10 keV energy range. The physical origin of this component is unknown. One possibility is that it arises from low-level accretion. The presence of this component in the spectrum complicates our cooling crust interpretation because it might indicate that the smooth luminosity and temperature decay curves we observed may not be due to crust cooling but due to some other process.

The LOFAR Two-metre Sky Survey (LoTSS) is a deep 120-168 MHz imaging survey that will eventually cover the entire Northern sky. Each of the 3170 pointings will be observed for 8 hrs, which, at most declinations, is sufficient to produce ~5arcsec resolution images with a sensitivity of ~0.1mJy/beam and accomplish the main scientific aims of the survey which are to explore the formation and evolution of massive black holes, galaxies, clusters of galaxies and large-scale structure. Due to the compact core and long baselines of LOFAR, the images provide excellent sensitivity to both highly extended and compact emission. For legacy value, the data are archived at high spectral and time resolution to facilitate subarcsecond imaging and spectral line studies. In this paper we provide an overview of the LoTSS. We outline the survey strategy, the observational status, the current calibration techniques, a preliminary data release, and the anticipated scientific impact. The preliminary images that we have released were created using a fully-automated but direction-independent calibration strategy and are significantly more sensitive than those produced by any existing large-area low-frequency survey. In excess of 44,000 sources are detected in the images that have a resolution of 25arcsec, typical noise levels of less than 0.5 mJy/beam, and cover an area of over 350 square degrees in the region of the HETDEX Spring Field (right ascension 10h45m00s to 15h30m00s and declination 45d00m00s to 57d00m00s).

Gravitational forces are expected to excite spiral density waves in protoplanetary disks, disks of gas and dust orbiting young stars. However, previous observations that showed spiral structure were not able to probe disk midplanes, where most of the mass is concentrated and where planet formation takes place. Using the Atacama Large Millimeter/submillimeter Array we detected a pair of trailing symmetric spiral arms in the protoplanetary disk surrounding the young star Elias 2-27. The arms extend to the disk outer regions and can be traced down to the midplane. These millimeter-wave observations also reveal an emission gap closer to the star than the spiral arms. We argue that the observed spirals trace shocks of spiral density waves in the midplane of this young disk.

Transitional millisecond pulsars (tMSPs) switch, on roughly multi-year timescales, between rotation-powered radio millisecond pulsar (RMSP) and accretion-powered low-mass X-ray binary (LMXB) states. The tMSPs have raised several questions related to the nature of accretion flow in their LMXB state and the mechanism that causes the state switch. The discovery of coherent X-ray pulsations from PSR J1023+0038 (while in the LMXB state) provides us with the first opportunity to perform timing observations and to compare the neutron star's spin variation during this state to the measured spin-down in the RMSP state. Whereas the X-ray pulsations in the LMXB state likely indicate that some material is accreting onto the neutron star's magnetic polar caps, radio continuum observations indicate the presence of an outflow. The fraction of the inflowing material being ejected is not clear, but it may be much larger than that reaching the neutron star's surface. Timing observations can measure the total torque on the neutron star. We have phase-connected nine XMM-Newton observations of PSR J1023+0038 over the last 2.5 years of the LMXB state to establish a precise measurement of spin evolution. We find that the average spin-down rate as an LMXB is 26.8+/-0.4% faster than the rate (-2.39x10^-15 Hz s-1) determined during the RMSP state. This shows that negative angular momentum contributions (dipolar magnetic braking and outflow) exceed positive ones (accreted material), and suggests that the pulsar wind continues to operate at a largely unmodified level. We discuss implications of this tight observational constraint in the context of possible accretion models.

(abridged). We outline LBCS (the LOFAR Long-Baseline Calibrator Survey), whose aim is to identify sources suitable for calibrating the highest-resolution observations made with the International LOFAR Telescope, which include baselines >1000 km. Suitable sources must contain significant correlated flux density (50-100mJy) at frequencies around 110--190~MHz on scales of a few hundred mas. At least for the 200--300-km international baselines, we find around 1 suitable calibrator source per square degree over a large part of the northern sky, in agreement with previous work. This should allow a randomly selected target to be successfully phase calibrated on the international baselines in over 50% of cases. Products of the survey include calibrator source lists and fringe-rate and delay maps of wide areas -- typically a few degrees -- around each source. The density of sources with significant correlated flux declines noticeably with baseline length over the range 200--600~km, with good calibrators on the longest baselines appearing only at the rate of 0.5 per square degree. Coherence times decrease from 1--3 minutes on 200-km baselines to about 1 minute on 600-km baselines, suggesting that ionospheric phase variations contain components with scales of a few hundred kilometres. The longest median coherence time, at just over 3 minutes, is seen on the DE609 baseline, which at 227km is close to being the shortest. We see median coherence times of between 80 and 110 seconds on the four longest baselines (580--600~km), and about 2 minutes for the other baselines. The success of phase transfer from calibrator to target is shown to be influenced by distance, in a manner that suggests a coherence patch at 150-MHz of the order of 1 degree.

We present the results of a closure phase analysis of 3 mm very long baseline interferometry (VLBI) measurements performed on Sagittarius A* (Sgr A*). We have analyzed observations made in May 2015 using the Very Long Baseline Array, the Robert C. Byrd Green Bank Telescope and the Large Millimeter Telescope Alfonso Serrano and obtained non-zero closure phase measurements on several station triangles - indicative of a non-point-symmetric source structure. The data are fitted with an asymmetric source structure model in Sgr A*, represented by a simple two-component model, which favours a fainter component due East of the main source. This result is discussed in light of a scattering screen with substructure or an intrinsically asymmetric source.

We examine the processes triggering star formation and Active Galactic Nucleus (AGN) activity in a sample of 25 low redshift ($z<0.13$) gas-rich galaxy mergers observed at milli-arcsecond resolution with Very Long Baseline Interferometry as part of the mJy Imaging VLBA Exploration at 20cm (mJIVE-20) survey. The high ($>10^7$ K) brightness temperature required for an mJIVE-20 detection allows us to unambiguously identify the radio AGN in our sample. We find three such objects. Our VLBI AGN identifications are classified as Seyferts or LINERs in narrow line optical diagnostic plots; mid-infrared colours of our targets and the comparison of H$\alpha$ star formation rates with integrated radio luminosity are also consistent with the VLBI identifications. We reconstruct star formation histories in our galaxies using optical and UV photometry, and find that these radio AGN are not triggered promptly in the merger process, consistent with previous findings for non-VLBI samples of radio AGN. This delay can significantly limit the efficiency of feedback by radio AGN triggered in galaxy mergers. We find that radio AGN hosts have lower star formation rates than non-AGN radio-selected galaxies at the same starburst age. Conventional and VLBI radio imaging shows these AGN to be compact on arcsecond scales. Our modeling suggests that the actual sizes of AGN-inflated radio lobes may be much larger than this, but these are too faint to be detected in existing observations. Deep radio imaging is required to map out the true extent of the AGN, and to determine whether the low star formation rates in radio AGN hosts are a result of the special conditions required for radio jet triggering, or the effect of AGN feedback.

We analyse new observations with the International Low Frequency Array (LOFAR) telescope, and archival data from the Multi-Element Radio Linked Interferometer Network (MERLIN) and the Karl G. Jansky Very Large Array (VLA). We model the spatially resolved radio spectrum of Arp 220 from 150 MHz to 33 GHz. We present an image of Arp 220 at 150 MHz with resolution $0.65''\times0.35''$, sensitivity 0.15 mJy beam$^{-1}$, and integrated flux density $394\pm59$ mJy. More than 80% of the detected flux comes from extended ($6''\approx$2.2 kpc) steep spectrum ($\alpha=-0.7$) emission, likely from star formation in the molecular disk surrounding the two nuclei. We find elongated features extending $0.3''$ (110 pc) and $0.9''$ (330 pc) from the eastern and western nucleus respectively, which we interpret as evidence for outflows. The extent of radio emission requires acceleration of cosmic rays far outside the nuclei. We find that a simple three component model can explain most of the observed radio spectrum of the galaxy. When accounting for absorption at 1.4 GHz, Arp 220 follows the FIR/radio correlation with $q=2.36$, and we estimate a star formation rate of 220 M$_\odot\text{yr}^{-1}$. We derive thermal fractions at 1 GHz of less than 1% for the nuclei, which indicates that a major part of the UV-photons are absorbed by dust. International LOFAR observations shows great promise to detect steep spectrum outflows and probe regions of thermal absorption. However, in LIRGs the emission detected at 150 MHz does not necessarily come from the main regions of star formation. This implies that high spatial resolution is crucial for accurate estimates of star formation rates for such galaxies at 150 MHz.

1RXS J180408.9-342058 is a transient neutron star low-mass X-ray binary that exhibited a bright accretion outburst in 2015. We present Nustar, Swift, and Chandra observations obtained around the peak of this outburst. The source was in a soft X-ray spectral state and displayed an X-ray luminosity of Lx~(2-3)E37 (D/5.8 kpc)^2 erg cm-2 (0.5-10 keV). The Nustar data reveal a broad Fe-K emission line that we model as relativistically broadened reflection to constrain the accretion geometry. We found that the accretion disk is viewed at an inclination of i~27-35 degrees and extended close to the neutron star, down to Rin~5-7.5 gravitational radii (~11-17 km). This inner disk radius suggests that the neutron star magnetic field strength is B<2E8 G. We find a narrow absorption line in the Chandra/HEG data at an energy of ~7.64 keV with a significance of ~4.8 sigma. This feature could correspond to blue-shifted Fe xxvi and arise from an accretion disk wind, which would imply an outflow velocity of v~0.086c (~25800 km s-1). However, this would be extreme for an X-ray binary and it is unclear if a disk wind should be visible at the low inclination angle that we infer from our reflection analysis. Finally, we discuss how the X-ray and optical properties of 1RXS J180408.9-342058 are consistent with a relatively small (Porb<3 hr) binary orbit.

Time domain radio astronomy observing campaigns frequently generate large volumes of data. Our goal is to develop automated methods that can identify events of interest buried within the larger data stream. The V-FASTR fast transient system was designed to detect rare fast radio bursts (FRBs) within data collected by the Very Long Baseline Array. The resulting event candidates constitute a significant burden in terms of subsequent human reviewing time. We have trained and deployed a machine learning classifier that marks each candidate detection as a pulse from a known pulsar, an artifact due to radio frequency interference, or a potential new discovery. The classifier maintains high reliability by restricting its predictions to those with at least 90% confidence. We have also implemented several efficiency and usability improvements to the V-FASTR web-based candidate review system. Overall, we found that time spent reviewing decreased and the fraction of interesting candidates increased. The classifier now classifies (and therefore filters) 80-90% of the candidates, with an accuracy greater than 98%, leaving only the 10-20% most promising candidates to be reviewed by humans.

The magneto-ionic structures of the interstellar medium of the Milky Way and the intergalactic medium are still poorly understood, especially at distances larger than a few kiloparsecs from the Sun. The three-dimensional (3D) structure of the Galactic magnetic field and electron density distribution may be probed through observations of radio pulsars, primarily owing to their compact nature, high velocities, and highly-polarized short-duration radio pulses. Phase 1 of the SKA, i.e. SKA1, will increase the known pulsar population by an order of magnitude, and the full SKA, i.e. SKA2, will discover pulsars in the most distant regions of our Galaxy. SKA1-VLBI will produce model-independent distances to a large number of pulsars, and wide-band polarization observations by SKA1-LOW and SKA1-MID will yield high precision dispersion measure, scattering measure, and rotation measure estimates along thousands of lines of sight. When combined, these observations will enable detailed tomography of the large-scale magneto-ionic structure of both the Galactic disk and the Galactic halo. Turbulence in the interstellar medium can be studied through the variations of these observables and the dynamic spectra of pulsar flux densities. SKA1-LOW and SKA1-MID will monitor interstellar weather and produce sensitive dynamic and secondary spectra of pulsar scintillation, which can be used to make speckle images of the ISM, study turbulence on scales between ~10^8 and ~10^13 m, and probe pulsar emission regions on scales down to $\sim$10 km. In addition, extragalactic pulsars or fast radio bursts to be discovered by SKA1 and SKA2 can be used to probe the electron density distribution and magnetic fields in the intergalactic medium beyond the Milky Way.

The correlation between radio spectral index and redshift has been exploited to discover high redshift radio galaxies, but its underlying cause is unclear. It is crucial to characterise the particle acceleration and loss mechanisms in high redshift radio galaxies to understand why their radio spectral indices are steeper than their local counterparts. Low frequency information on scales of $\sim$1 arcsec are necessary to determine the internal spectral index variation. In this paper we present the first spatially resolved studies at frequencies below 100 MHz of the $z = 2.4$ radio galaxy 4C 43.15 which was selected based on its ultra-steep spectral index ($\alpha < -1$; $S_{\nu} \sim \nu^{\alpha}$ ) between 365 MHz and 1.4 GHz. Using the International Low Frequency Array (LOFAR) Low Band Antenna we achieve sub-arcsecond imaging resolution at 55 MHz with VLBI techniques. Our study reveals low-frequency radio emission extended along the jet axis, which connects the two lobes. The integrated spectral index for frequencies $<$ 500 MHz is -0.83. The lobes have integrated spectral indices of -1.31$\pm$0.03 and -1.75$\pm$0.01 for frequencies $\geq$1.4 GHz, implying a break frequency between 500 MHz and 1.4 GHz. These spectral properties are similar to those of local radio galaxies. We conclude that the initially measured ultra-steep spectral index is due to a combination of the steepening spectrum at high frequencies with a break at intermediate frequencies.

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.

We present a series of sixteen Very Long Baseline Array (VLBA) high angular resolution observations of a cluster of suspected low-mass young stars in the Monoceros R2 region. Four compact and highly variable radio sources are detected; three of them in only one epoch, the fourth one a total of seven times. This latter source is seen in the direction to the previously known \UCHII\ region VLA~1, and has radio properties that resemble those of magnetically active stars; we shall call it VLA~1$^\star$. We model its displacement on the celestial sphere as a combination of proper motion and trigonometric parallax. The fit obtained using a uniform proper motion yields a parallax $\varpi$ = 1.10 $\pm$ 0.18 mas, but with a fairly high post-fit dispersion. If acceleration terms (probably due to an undetected companion) are included, the quality of the fit improves dramatically, and the best estimate of the parallax becomes $\varpi$ = 1.12 $\pm$ 0.05 mas. The magnitude of the fitted acceleration suggest an orbital period of order a decade. The measured parallax corresponds to a distance $d$ = 893$^{+44}_{-40}$ pc, in very good agreement with previous, indirect, determinations.

The V-FASTR experiment on the Very Long Baseline Array was designed to detect dispersed pulses of milliseconds duration, such as fast radio bursts (FRBs). We use all V-FASTR data through February 2015 to report V-FASTR's upper limits on the rates of FRBs, and compare these with re-derived rates from Parkes FRB detection experiments. V-FASTR's operation at lambda=20 cm allows direct comparison with the 20 cm Parkes rate, and we derive a power-law limit of \gamma<-0.4 (95% confidence limit) on the index of FRB source counts, N(>S)\propto S^\gamma. Using the previously measured FRB rate and the unprecedented amount of survey time spent searching for FRBs at a large range of wavelengths (0.3 cm > \lambda > 90 cm), we also place frequency-dependent limits on the spectral distribution of FRBs. The most constraining frequencies place two-point spectral index limits of \alpha_{20cm}^{4cm} < 5.8 and \alpha_{90cm}^{20cm} > -7.6, where fluence F \propto f^\alpha if we assume true the burst rate reported by Champion et al. (2016) of R(F~0.6 Jy ms) = 7 x 10^3 sky^{-1} day^{-1} (for bursts of ~3 ms duration). This upper limit on \alpha suggests that if FRBs are extragalactic but non-cosmological, that on average they are not experiencing excessive free-free absorption due to a medium with high optical depth (assuming temperature ~8,000 K), which excessively invert their low-frequency spectrum. This in turn implies that the dispersion of FRBs arises in either or both of the intergalactic medium or the host galaxy, rather than from the source itself.

Cosmic rays are the highest energy particles found in nature. Measurements of the mass composition of cosmic rays between 10^{17} eV and 10^{18} eV are essential to understand whether this energy range is dominated by Galactic or extragalactic sources. It has also been proposed that the astrophysical neutrino signal comes from accelerators capable of producing cosmic rays of these energies. Cosmic rays initiate cascades of secondary particles (air showers) in the atmosphere and their masses are inferred from measurements of the atmospheric depth of the shower maximum, Xmax, or the composition of shower particles reaching the ground. Current measurements suffer from either low precision, or a low duty cycle and a high energy threshold. Radio detection of cosmic rays is a rapidly developing technique, suitable for determination of Xmax with a duty cycle of in principle nearly 100%. The radiation is generated by the separation of relativistic charged particles in the geomagnetic field and a negative charge excess in the shower front. Here we report radio measurements of Xmax with a mean precision of 16 g/cm^2 between 10^{17}-10^{17.5} eV. Because of the high resolution in $Xmax we can determine the mass spectrum and find a mixed composition, containing a light mass fraction of ~80%. Unless the extragalactic component becomes significant already below 10^{17.5} eV, our measurements indicate an additional Galactic component dominating at this energy range.

We measure the proper motion of the pulsar PSR J1745-2900 relative to the Galactic Center massive black hole, Sgr A*, using the Very Long Baseline Array (VLBA). The pulsar has a transverse velocity of 236 +/- 11 km s^-1 at position angle 22 +/- 2 deg East of North at a projected separation of 0.097 pc from Sgr A*. Given the unknown radial velocity, this transverse velocity measurement does not conclusively prove that the pulsar is bound to Sgr A*; however, the probability of chance alignment is very small. We do show that the velocity and position is consistent with a bound orbit originating in the clockwise disk of massive stars orbiting Sgr A* and a natal velocity kick of <~ 500 km s^-1. An origin among the isotropic stellar cluster is possible but less probable. If the pulsar remains radio-bright, multi-year astrometry of PSR J1745-2900 can detect its acceleration and determine the full three-dimensional orbit. We also demonstrate that PSR J1745-2900 exhibits the same angular broadening as Sgr A* over a wavelength range of 3.6 cm to 0.7 cm, further confirming that the two sources share the same interstellar scattering properties. Finally, we place the first limits on the presence of a wavelength-dependent shift in the position of Sgr A*, i.e., the core shift, one of the expected properties of optically-thick jet emission. Our results for PSR J1745-2900 support the hypothesis that Galactic Center pulsars will originate from the stellar disk and deepens the mystery regarding the small number of detected Galactic Center pulsars.

To characterize the mechanisms of planet formation it is crucial to investigate the properties and evolution of protoplanetary disks around young stars, where the initial conditions for the growth of planets are set. Our goal is to study grain growth in the disk of the young, intermediate mass star HD163296 where dust processing has already been observed, and to look for evidence of growth by ice condensation across the CO snowline, already identified in this disk with ALMA. Under the hypothesis of optically thin emission we compare images at different wavelengths from ALMA and VLA to measure the opacity spectral index across the disk and thus the maximum grain size. We also use a Bayesian tool based on a two-layer disk model to fit the observations and constrain the dust surface density. The measurements of the opacity spectral index indicate the presence of large grains and pebbles ($\geq$1 cm) in the inner regions of the disk (inside $\sim$50 AU) and smaller grains, consistent with ISM sizes, in the outer disk (beyond 150 AU). Re-analysing ALMA Band 7 Science Verification data we find (radially) unresolved excess continuum emission centered near the location of the CO snowline at $\sim$90 AU. Our analysis suggests a grain size distribution consistent with an enhanced production of large grains at the CO snowline and consequent transport to the inner regions. Our results combined with the excess in infrared scattered light found by Garufi et al. (2014) suggests the presence of a structure at 90~AU involving the whole vertical extent of the disk. This could be evidence for small scale processing of dust at the CO snowline.

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).

Theoretical models of grain growth predict dust properties to change as a function of protoplanetary disk radius, mass, age and other physical conditions. We lay down the methodology for a multi-wavelength analysis of (sub-)mm and cm continuum interferometric observations to constrain self-consistently the disk structure and the radial variation of the dust properties. The computational architecture is massively parallel and highly modular. The analysis is based on the simultaneous fit in the uv-plane of observations at several wavelengths with a model for the disk thermal emission and for the dust opacity. The observed flux density at the different wavelengths is fitted by posing constraints on the disk structure and on the radial variation of the grain size distribution. We apply the analysis to observations of three protoplanetary disks (AS 209, FT Tau, DR Tau) for which a combination of spatially resolved observations in the range ~0.88mm to ~10mm is available (from SMA, CARMA, and VLA), finding evidence of a decreasing maximum dust grain size (a_max) with radius. We derive large a_max values up to 1 cm in the inner disk between 15 and 30 AU and smaller grains with a_max~1 mm in the outer disk (R > 80AU). In this paper we develop a multi-wavelength analysis that will allow this missing quantity to be constrained for statistically relevant samples of disks and to investigate possible correlations with disk or stellar parameters.

We present the results of a four-month campaign searching for low-frequency radio transients near the North Celestial Pole with the Low-Frequency Array (LOFAR), as part of the Multifrequency Snapshot Sky Survey (MSSS). The data were recorded between 2011 December and 2012 April and comprised 2149 11-minute snapshots, each covering 175 deg^2. We have found one convincing candidate astrophysical transient, with a duration of a few minutes and a flux density at 60 MHz of 15-25 Jy. The transient does not repeat and has no obvious optical or high-energy counterpart, as a result of which its nature is unclear. The detection of this event implies a transient rate at 60 MHz of 3.9 (+14.7, -3.7) x 10^-4 day^-1 deg^-2, and a transient surface density of 1.5 x 10^-5 deg^-2, at a 7.9-Jy limiting flux density and ~10-minute time-scale. The campaign data were also searched for transients at a range of other time-scales, from 0.5 to 297 min, which allowed us to place a range of limits on transient rates at 60 MHz as a function of observation duration.

Context. Observing Jupiter's synchrotron emission from the Earth remains today the sole method to scrutinize the distribution and dynamical behavior of the ultra energetic electrons magnetically trapped around the planet (because in-situ particle data are limited in the inner magnetosphere). Aims. We perform the first resolved and low-frequency imaging of the synchrotron emission with LOFAR at 127 MHz. The radiation comes from low energy electrons (~1-30 MeV) which map a broad region of Jupiter's inner magnetosphere. Methods (see article for complete abstract) Results. The first resolved images of Jupiter's radiation belts at 127-172 MHz are obtained along with total integrated flux densities. They are compared with previous observations at higher frequencies and show a larger extent of the synchrotron emission source (>=4 $R_J$). The asymmetry and the dynamic of east-west emission peaks are measured and the presence of a hot spot at lambda_III=230 {\deg} $\pm$ 25 {\deg}. Spectral flux density measurements are on the low side of previous (unresolved) ones, suggesting a low-frequency turnover and/or time variations of the emission spectrum. Conclusions. LOFAR is a powerful and flexible planetary imager. The observations at 127 MHz depict an extended emission up to ~4-5 planetary radii. The similarities with high frequency results reinforce the conclusion that: i) the magnetic field morphology primarily shapes the brightness distribution of the emission and ii) the radiating electrons are likely radially and latitudinally distributed inside about 2 $R_J$. Nonetheless, the larger extent of the brightness combined with the overall lower flux density, yields new information on Jupiter's electron distribution, that may shed light on the origin and mode of transport of these particles.

LOFAR offers the unique capability of observing pulsars across the 10-240 MHz frequency range with a fractional bandwidth of roughly 50%. This spectral range is well-suited for studying the frequency evolution of pulse profile morphology caused by both intrinsic and extrinsic effects: such as changing emission altitude in the pulsar magnetosphere or scatter broadening by the interstellar medium, respectively. The magnitude of most of these effects increases rapidly towards low frequencies. LOFAR can thus address a number of open questions about the nature of radio pulsar emission and its propagation through the interstellar medium. We present the average pulse profiles of 100 pulsars observed in the two LOFAR frequency bands: High Band (120-167 MHz, 100 profiles) and Low Band (15-62 MHz, 26 profiles). We compare them with Westerbork Synthesis Radio Telescope (WSRT) and Lovell Telescope observations at higher frequencies (350 and1400 MHz) in order to study the profile evolution. The profiles are aligned in absolute phase by folding with a new set of timing solutions from the Lovell Telescope, which we present along with precise dispersion measures obtained with LOFAR. We find that the profile evolution with decreasing radio frequency does not follow a specific trend but, depending on the geometry of the pulsar, new components can enter into, or be hidden from, view. Nonetheless, in general our observations confirm the widening of pulsar profiles at low frequencies, as expected from radius-to-frequency mapping or birefringence theories. We offer this catalog of low-frequency pulsar profiles in a user friendly way via the EPN Database of Pulsar Profiles (http://www.epta.eu.org/epndb/).

Energetic feedback from active galactic nuclei (AGN) is an important ingredient for regulating the star-formation history of galaxies in models of galaxy formation, which makes it important to study how AGN feedback actually occurs in practice. In order to catch AGNs in the act of quenching star formation we have used the interstellar NaD absorption lines to look for cold-gas outflows in a sample of 456 nearby galaxies for which we could unambigously ascertain the presence of radio AGN activity, thanks to radio imaging at milli-arcsecond scales. While compact radio emission indicating a radio AGN was found in 103 galaxies (23% of the sample), and 23 objects (5%) exhibited NaD absorption-line kinematics suggestive of cold-gas outflows, not one object showed evidence of a radio AGN and of a cold-gas outflow simultaneously. Radio AGN activity was found predominantly in early-type galaxies, while cold-gas outflows were mainly seen in spiral galaxies with central star-formation or composite star-formation/AGN activity. Optical AGNs also do not seem capable of driving galactic winds in our sample. Our work adds to a picture of the low-redshift Universe where cold-gas outflows in massive galaxies are generally driven by star formation and where radio-AGN activity occurs most often in systems in which the gas reservoir has already been significantly depleted.

We present new results from the Disks@EVLA program for two young stars: CY Tau and DoAr 25. We trace continuum emission arising from their circusmtellar disks from spatially resolved observations, down to tens of AU scales, at {\lambda} = 0.9, 2.8, 8.0, and 9.8 mm for DoAr25 and at {\lambda} = 1.3, 2.8, and 7.1 mm for CY Tau. Additionally, we constrain the amount of emission whose origin is different from thermal dust emission from 5 cm observations. Directly from interferometric data, we find that observations at 7 mm and 1 cm trace emission from a compact disk while millimeter-wave observations trace an extended disk structure. From a physical disk model, where we characterize the disk structure of CY Tau and DoAr 25 at wavelengths shorter than 5 cm, we find that (1) dust continuum emission is optically thin at the observed wavelengths and over the spatial scales studied, (2) a constant value of the dust opacity is not warranted by our observations, and (3) a high-significance radial gradient of the dust opacity spectral index, {\beta}, is consistent with the observed dust emission in both disks, with low-{\beta} in the inner disk and high-{\beta} in the outer disk. Assuming that changes in dust properties arise solely due to changes in the maximum particle size (amax), we constrain radial variations of amax in both disks, from cm-sized particles in the inner disk (R < 40 AU) to millimeter sizes in the outer disk (R > 80 AU). These observational constraints agree with theoretical predictions of the radial-drift barrier, however, fragmentation of dust grains could explain our amax(R) constraints if these disks have lower turbulence and/or if dust can survive high-velocity collisions.

We present the Multifrequency Snapshot Sky Survey (MSSS), the first northern-sky LOFAR imaging survey. In this introductory paper, we first describe in detail the motivation and design of the survey. Compared to previous radio surveys, MSSS is exceptional due to its intrinsic multifrequency nature providing information about the spectral properties of the detected sources over more than two octaves (from 30 to 160 MHz). The broadband frequency coverage, together with the fast survey speed generated by LOFAR's multibeaming capabilities, make MSSS the first survey of the sort anticipated to be carried out with the forthcoming Square Kilometre Array (SKA). Two of the sixteen frequency bands included in the survey were chosen to exactly overlap the frequency coverage of large-area Very Large Array (VLA) and Giant Metrewave Radio Telescope (GMRT) surveys at 74 MHz and 151 MHz respectively. The survey performance is illustrated within the "MSSS Verification Field" (MVF), a region of 100 square degrees centered at J2000 (RA,Dec)=(15h,69deg). The MSSS results from the MVF are compared with previous radio survey catalogs. We assess the flux and astrometric uncertainties in the catalog, as well as the completeness and reliability considering our source finding strategy. We determine the 90% completeness levels within the MVF to be 100 mJy at 135 MHz with 108" resolution, and 550 mJy at 50 MHz with 166" resolution. Images and catalogs for the full survey, expected to contain 150,000-200,000 sources, will be released to a public web server. We outline the plans for the ongoing production of the final survey products, and the ultimate public release of images and source catalogs.

The transitional millisecond pulsar binary system PSR J1023+0038 re-entered an accreting state in 2013 June, in which it bears many similarities to low-mass X-ray binaries (LMXBs) in quiescence or near-quiescence. At a distance of just 1.37 kpc, PSR J1023+0038 offers an unsurpassed ability to study low-level accretion onto a highly-magnetized compact object. We have monitored PSR J1023+0038 intensively using radio imaging with the Karl G. Jansky Very Large Array, the European VLBI Network and LOFAR, seeing rapidly variable, flat spectrum emission that persists over a period of 6 months. The flat spectrum and variability are indicative of synchrotron emission originating in an outflow from the system, most likely in the form of a compact, partially self-absorbed jet, as is seen in LMXBs at higher accretion rates. The radio brightness, however, greatly exceeds extrapolations made from observations of more vigorously accreting neutron star LMXB systems. We postulate that PSR J1023+0038 is undergoing radiatively inefficient 'propeller-mode' accretion, with the jet carrying away a dominant fraction of the liberated accretion luminosity. We confirm that the enhanced gamma-ray emission seen in PSR J1023+0038 since it re-entered an accreting state has been maintained; the increased gamma-ray emission in this state can also potentially be associated with propeller-mode accretion. Similar accretion modes can be invoked to explain the radio and X-ray properties of the other two known transitional millisecond pulsar systems XSS J12270-4859 and PSR J1824-2452I (M28I), suggesting that radiatively inefficient accretion may be an ubiquitous phenomenon amongst (at least one class of) neutron star binaries at low accretion rates.

Our analysis of archival VLBI data of PSR 0834+06 revealed that its scintillation properties can be precisely modelled using the inclined sheet model (Pen & Levin 2014), resulting in two distinct lens planes. These data strongly favour the grazing sheet model over turbulence as the primary source of pulsar scattering. This model can reproduce the parameters of the observed diffractive scintillation with an accuracy at the percent level. Comparison with new VLBI proper motion results in a direct measure of the ionized ISM screen transverse velocity. The results are consistent with ISM velocities local to the PSR 0834+06 sight-line (through the Galaxy). The simple 1D structure of the lenses opens up the possibility of using interstellar lenses as precision probes for pulsar lens mapping, precision transverse motions in the ISM, and new opportunities for removing scattering to improve pulsar timing. We describe the parameters and observables of this double screen system. While relative screen distances can in principle be accurately determined, a global conformal distance degeneracy exists that allows a rescaling of the absolute distance scale. For PSR B0834+06, we present VLBI astrometry results that provide (for the fist time) a direct measurement of the distance of the pulsar. For targets where independent distance measurements are not available, which are the cases for most of the recycled millisecond pulsars that are the targets of precision timing observations, the degeneracy presented in the lens modelling could be broken if the pulsar resides in a binary system.

This chapter provides an overview of the possibilities for transient and variable-source astrophysics with the Square Kilometre Array. While subsequent chapters focus on the astrophysics of individual events, we focus on the broader picture, and how to maximise the science coming from the telescope. The SKA as currently designed will be a fantastic and ground-breaking facility for radio transient studies, but the scientifc yield will be dramatically increased by the addition of (i) near-real-time commensal searches of data streams for events, and (ii) on occasion, rapid robotic response to Target-of-Opprtunity style triggers.

Accreting millisecond X-ray pulsars are an important subset of low-mass X-ray binaries in which coherent X-ray pulsations can be observed during occasional, bright outbursts (X-ray luminosity $L_X\sim 10^{36}$ erg s$^{-1}$). These pulsations show that matter is being channeled onto the neutron star's magnetic poles. However, such sources spend most of their time in a low-luminosity, quiescent state ($L_X\lesssim 10^{34}$ erg s$^{-1}$), where the nature of the accretion flow onto the neutron star (if any) is not well understood. Here we report that the millisecond pulsar/low-mass X-ray binary transition object PSR J1023+0038 intermittently shows coherent X-ray pulsations at luminosities nearly 100 times fainter than observed in any other accreting millisecond X-ray pulsar. We conclude that in spite of its low luminosity PSR J1023+0038 experiences episodes of channeled accretion, a discovery that challenges existing models for accretion onto magnetized neutron stars.

Using observations obtained with the LOw Fequency ARray (LOFAR), the Westerbork Synthesis Radio Telescope (WSRT) and archival Very Large Array (VLA) data, we have traced the radio emission to large scales in the complex source 4C 35.06 located in the core of the galaxy cluster Abell 407. At higher spatial resolution (~4"), the source was known to have two inner radio lobes spanning 31 kpc and a diffuse, low-brightness extension running parallel to them, offset by about 11 kpc (in projection). At 62 MHz, we detect the radio emission of this structure extending out to 210 kpc. At 1.4 GHz and intermediate spatial resolution (~30"), the structure appears to have a helical morphology. We have derived the characteristics of the radio spectral index across the source. We show that the source morphology is most likely the result of at least two episodes of AGN activity separated by a dormant period of around 35 Myr. The AGN is hosted by one of the galaxies located in the cluster core of Abell 407. We propose that it is intermittently active as it moves in the dense environment in the cluster core. Using LOFAR, we can trace the relic plasma from that episode of activity out to greater distances from the core than ever before. Using the the WSRT, we detect HI in absorption against the center of the radio source. The absorption profile is relatively broad (FWHM of 288 km/s), similar to what is found in other clusters. Understanding the duty cycle of the radio emission as well as the triggering mechanism for starting (or restarting) the radio-loud activity can provide important constraints to quantify the impact of AGN feedback on galaxy evolution. The study of these mechanisms at low frequencies using morphological and spectral information promises to bring new important insights in this field.

The PSR J1023+0038 binary system hosts a neutron star and a low-mass, main-sequence-like star. It switches on year timescales between states as an eclipsing radio millisecond pulsar and a low-mass X-ray binary. We present a multi-wavelength observational campaign of PSR J1023+0038 in its most recent low-mass X-ray binary state. Two long XMM-Newton observations reveal that the system spends ~70% of the time in a $\approx$$3\times10^{33}$ erg/s X-ray luminosity mode, which, as shown in Archibald et al. (2014), exhibits coherent X-ray pulsations. This emission is interspersed with frequent lower flux mode intervals with $\approx$$5\times 10^{32}$ erg/s and sporadic flares reaching up to $\approx$$10^{34}$ erg/s, with neither mode showing significant X-ray pulsations. The switches between the three flux modes occur on timescales of order 10 s. In the UV and optical, we observe occasional intense flares coincident with those observed in X-rays. Our radio timing observations reveal no pulsations at the pulsar period during any of the three X-ray modes, presumably due to complete quenching of the radio emission mechanism by the accretion flow. Radio imaging detects highly variable, flat-spectrum continuum emission from PSR J1023+0038, consistent with an origin in a weak jet-like outflow. Our concurrent X-ray and radio continuum data sets do not exhibit any correlated behavior. The observational evidence we present bears qualitative resemblance to the behavior predicted by some existing propeller and trapped disk accretion models although none can fully account for all aspects of the rich phenomenology of this system.

Infrared-faint radio sources (IFRS) form a new class of galaxies characterised by radio flux densities between tenths and tens of mJy and faint or absent infrared counterparts. It has been suggested that these objects are radio-loud active galactic nuclei (AGNs) at significant redshifts (z >~ 2). Whereas the high redshifts of IFRS have been recently confirmed based on spectroscopic data, the evidence for the presence of AGNs in IFRS is mainly indirect. So far, only two AGNs have been unquestionably confirmed in IFRS based on very long baseline interferometry (VLBI) observations. In this work, we test the hypothesis that IFRS contain AGNs in a large sample of sources using VLBI. We observed 57 IFRS with the Very Long Baseline Array (VLBA) down to a detection sensitivity in the sub-mJy regime and detected compact cores in 35 sources. Our VLBA detections increase the number of VLBI-detected IFRS from 2 to 37 and provide strong evidence that most - if not all - IFRS contain AGNs. We find that IFRS have a marginally higher VLBI detection fraction than randomly selected sources with mJy flux densities at arcsec-scales. Moreover, our data provide a positive correlation between compactness - defined as the ratio of milliarcsec- to arcsec-scale flux density - and redshift for IFRS, but suggest a decreasing mean compactness with increasing arcsec-scale radio flux density. Based on these findings, we suggest that IFRS tend to contain young AGNs whose jets have not formed yet or have not expanded, equivalent to very compact objects. We found two IFRS that are resolved into two components. The two components are spatially separated by a few hundred milliarcseconds in both cases. They might be components of one AGN, a binary black hole, or the result of gravitational lensing.

Impulsive radio bursts that are detectable across cosmological distances constitute extremely powerful probes of the ionized Inter-Galactic Medium (IGM), intergalactic magnetic fields, and the properties of space-time itself. Their dispersion measures (DMs) will enable us to detect the "missing" baryons in the low-redshift Universe and make the first measurements of the mean galaxy halo profile, a key parameter in models of galaxy formation and feedback. Impulsive bursts can be used as cosmic rulers at redshifts exceeding 2, and constrain the dark energy equation-of-state parameter, $w(z)$ at redshifts beyond those readily accessible by Type Ia SNe. Both of these goals are realisable with a sample of $\sim 10^4$ fast radio bursts (FRBs) whose positions are localized to within one arcsecond, sufficient to obtain host galaxy redshifts via optical follow-up. It is also hypothesised that gravitational wave events may emit coherent emission at frequencies probed by SKA1-LOW, and the localization of such events at cosmological distances would enable their use as cosmological standard sirens. To perform this science, such bursts must be localized to their specific host galaxies so that their redshifts may be obtained and compared against their dispersion measures, rotation measures, and scattering properties. The SKA can achieve this with a design that has a wide field-of-view, a substantial fraction of its collecting area in a compact configuration (80\% within a 3\,km radius), and a capacity to attach high-time-resolution instrumentation to its signal path.

On a time scale of years to decades, gravitational wave (GW) astronomy will become a reality. Low frequency (nanoHz) GWs are detectable through long-term timing observations of the most stable pulsars. Radio observatories worldwide are currently carrying out observing programmes to detect GWs, with data sets being shared through the International Pulsar Timing Array project. One of the most likely sources of low frequency GWs are supermassive black hole binaries (SMBHBs), detectable as a background due to a large number of binaries, or as continuous or burst emission from individual sources. No GW signal has yet been detected, but stringent constraints are already being placed on galaxy evolution models. The SKA will bring this research to fruition. In this chapter, we describe how timing observations using SKA1 will contribute to detecting GWs, or can confirm a detection if a first signal already has been identified when SKA1 commences observations. We describe how SKA observations will identify the source(s) of a GW signal, search for anisotropies in the background, improve models of galaxy evolution, test theories of gravity, and characterise the early inspiral phase of a SMBHB system. We describe the impact of the large number of millisecond pulsars to be discovered by the SKA; and the observing cadence, observation durations, and instrumentation required to reach the necessary sensitivity. We describe the noise processes that will influence the achievable precision with the SKA. We assume a long-term timing programme using the SKA1-MID array and consider the implications of modifications to the current design. We describe the possible benefits from observations using SKA1-LOW. Finally, we describe GW detection prospects with SKA1 and SKA2, and end with a description of the expectations of GW astronomy.

Since their discovery in the late 1960's the population of known neutron stars (NSs) has grown to ~2500. The last five decades of observations have yielded many surprises and demonstrated that the observational properties of NSs are remarkably diverse. The surveys that will be performed with SKA (the Square Kilometre Array) will produce a further tenfold increase in the number of Galactic NSs known. Moreover, the SKA's broad spectral coverage, sub-arraying and multi-beaming capabilities will allow us to characterise these sources with unprecedented efficiency, in turn enabling a giant leap in the understanding of their properties. Here we review the NS population and outline our strategies for studying each of the growing number of diverse classes that are populating the "NS zoo". Some of the main scientific questions that will be addressed by the much larger statistical samples and vastly improved timing efficiency provided by SKA include: (i) the spin period and spin-down rate distributions (and thus magnetic fields) at birth, and the associated information about the SNe wherein they are formed; (ii) the radio pulsar-magnetar connection; (iii) the link between normal radio pulsars, intermittent pulsars and rotating radio transients; (iv) the slowest possible spin period for a radio pulsar (revealing the conditions at the pulsar death-line); (v) proper motions of pulsars (revealing SN kick physics); (vi) the mass distribution of NSs (vii) the fastest possible spin period for a recycled pulsar (constraining magnetosphere-accretion disc interactions, gravitational wave radiation and the equation-of-state); (viii) the origin of high eccentricity millisecond pulsars (MSPs); (ix) the formation channels for recently identified triple systems; and finally (x) how isolated MSPs are formed. We expect that the SKA will break new ground unveiling exotic systems that will challenge... [abridged]

The Square Kilometre Array (SKA) will use pulsars to enable precise measurements of strong gravity effects in pulsar systems, which yield tests of gravitational theories that cannot be carried out anywhere else. The Galactic census of pulsars will discover dozens of relativistic pulsar systems, possibly including pulsar -- black hole binaries which can be used to test the "cosmic censorship conjecture" and the "no-hair theorem". Also, the SKA's remarkable sensitivity will vastly improve the timing precision of millisecond pulsars, allowing probes of potential deviations from general relativity (GR). Aspects of gravitation to be explored include tests of strong equivalence principles, gravitational dipole radiation, extra field components of gravitation, gravitomagnetism, and spacetime symmetries.

Adding VLBI capability to the SKA arrays will greatly broaden the science of the SKA, and is feasible within the current specifications. SKA-VLBI can be initially implemented by providing phased-array outputs for SKA1-MID and SKA1-SUR and using these extremely sensitive stations with other radio telescopes, and in SKA2 by realising a distributed configuration providing baselines up to thousands of km, merging it with existing VLBI networks. The motivation for and the possible realization of SKA-VLBI is described in this paper.

Radio AGN in the nearby Universe are more likely to be found in galaxies with early-type morphology, the detection rate in spiral or late-type galaxies (LTGs) being around an order of magnitude lower. We combine the mJy Imaging VLBA Exploration at 20cm (mJIVE-20) survey with the Sloan Digital Sky Survey (SDSS), to study the relatively rare population of AGN in LTGs that have nuclear radio luminosities similar to that in their early-type counterparts. The LTG AGN population is preferentially hosted by galaxies that have high stellar masses (M* > 10^10.8 MSun), red colours and low star-formation rates, with little dependence on the detailed morphology or local environment of the host LTG. The merger fraction in the LTG AGN is around 4 times higher than that in the general LTG population, indicating that merging is an important trigger for radio AGN in these systems. The red colours of our systems extend recent work which indicates that merger-triggered AGN in the nearby Universe appear after the peak of the associated starburst, implying that they do not strongly regulate star formation. Finally, we find that in systems where parsec-scale jets are clearly observed in our VLBI images, the jets are perpendicular to the major axis of the galaxy, indicating alignment between the accretion disc and the host galaxy stellar disc.

The nuclear starburst in the nearby galaxy M82 provides an excellent laboratory for understanding the physics of star formation. This galaxy has been extensively observed in the past, revealing tens of radio-bright compact objects embedded in a diffuse free-free absorbing medium. Our understanding of the structure and physics of this medium in M82 can be greatly improved by high-resolution images at low frequencies where the effects of free-free absorption are most prominent. The aims of this study are, firstly, to demonstrate imaging using international baselines of the Low Frequency Array (LOFAR), and secondly, to constrain low-frequency spectra of compact and diffuse emission in the central starburst region of M82 via high-resolution radio imaging at low frequencies. The international LOFAR telescope was used to observe M82 at 110-126MHz and 146-162MHz. Images were obtained using standard techniques from very long baseline interferometry. images were obtained at each frequency range: one only using international baselines, and one only using the longest Dutch (remote) baselines. The 154MHz image obtained using international baselines is a new imaging record in terms of combined image resolution (0.3$"$) and sensitivity ($\sigma$=0.15mJy/beam) at low frequencies ($<327$MHz). We detected 16 objects at 154MHz, six of these also at 118MHz. Four weaker but resolved features are also found: a linear (50pc) filament and three other resolved objects, of which two show a clear shell structure. We do not detect any emission from either supernova 2008iz or from the radio transient source 43.78+59.3. The images obtained using remote baselines show diffuse emission, associated with the outflow in M82, with reduced brightness in the region of the edge-on star-forming disk.

Measuring radio emission from air showers offers a novel way to determine properties of the primary cosmic rays such as their mass and energy. Theory predicts that relativistic time compression effects lead to a ring of amplified emission which starts to dominate the emission pattern for frequencies above ~100 MHz. In this article we present the first detailed measurements of this structure. Ring structures in the radio emission of air showers are measured with the LOFAR radio telescope in the frequency range of 110 - 190 MHz. These data are well described by CoREAS simulations. They clearly confirm the importance of including the index of refraction of air as a function of height. Furthermore, the presence of the Cherenkov ring offers the possibility for a geometrical measurement of the depth of shower maximum, which in turn depends on the mass of the primary particle.

Aims. An efficient means of locating calibrator sources for International LOFAR is developed and used to determine the average density of usable calibrator sources on the sky for subarcsecond observations at 140 MHz. Methods. We used the multi-beaming capability of LOFAR to conduct a fast and computationally inexpensive survey with the full International LOFAR array. Sources were pre-selected on the basis of 325 MHz arcminute-scale flux density using existing catalogues. By observing 30 different sources in each of the 12 sets of pointings per hour, we were able to inspect 630 sources in two hours to determine if they possess a sufficiently bright compact component to be usable as LOFAR delay calibrators. Results. Over 40% of the observed sources are detected on multiple baselines between international stations and 86 are classified as satisfactory calibrators. We show that a flat low-frequency spectrum (from 74 to 325 MHz) is the best predictor of compactness at 140 MHz. We extrapolate from our sample to show that the density of calibrators on the sky that are sufficiently bright to calibrate dispersive and non-dispersive delays for the International LOFAR using existing methods is 1.0 per square degree. Conclusions. The observed density of satisfactory delay calibrator sources means that observations with International LOFAR should be possible at virtually any point in the sky, provided that a fast and efficient search using the methodology described here is conducted prior to the observation to identify the best calibrator.

We explore the processes that trigger local AGN and the role of these AGN in regulating star formation, using ~350 nearby galaxies observed by the mJy Imaging VLBA Exploration at 20cm (mJIVE) survey. The >10^7 K brightness temperature required for an mJIVE detection cannot be achieved via star formation alone, allowing us to unambiguously detect nearby radio AGN and study their role in galaxy evolution. Radio AGN are an order of magnitude more common in early-type galaxies (ETGs) than in their late-type counterparts. The VLBI-detected ETGs in this study have a similar stellar mass distribution to their undetected counterparts, are typically not the central galaxies of clusters and exhibit merger fractions that are significantly higher than in the average ETG. This suggests that these radio AGN (which have VLBI luminosities >10^22 W Hz^-1) are primarily fuelled by mergers, and not by internal stellar mass loss or cooling flows. Our radio AGN are a factor of ~3 times more likely to reside in the UV-optical red sequence than the average ETG. Furthermore, typical AGN lifetimes (a few 10^7 yr) are much shorter than the transit times from blue cloud to red sequence (~1.5 Gyr). This indicates that the AGN are not triggered promptly and appear several dynamical timescales into the associated star formation episode, implying that they typically couple only to residual gas, at a point where star formation has already declined significantly. While evidence for AGN feedback is strong in systems where the black hole is fed by the cooling of hot gas, AGN triggered by mergers appear not to strongly regulate the associated star formation. The inability of the AGN to rapidly quench merger-driven star formation is likely to make merging the dominant mode of star formation in nearby ETGs, in line with evidence for minor mergers being the primary driver of stellar mass growth in these systems.

We present LOFAR Low Band observations of the Bootes and 3C295 fields. Our images made at 34, 46, and 62 MHz reach noise levels of 12, 8, and 5 mJy beam$^{-1}$, making them the deepest images ever obtained in this frequency range. In total, we detect between 300 and 400 sources in each of these images, covering an area of 17 to 52 deg$^{2}$. From the observations we derive Euclidean-normalized differential source counts. The 62 MHz source counts agree with previous GMRT 153 MHz and VLA 74 MHz differential source counts, scaling with a spectral index of $-0.7$. We find that a spectral index scaling of $-0.5$ is required to match up the LOFAR 34 MHz source counts. This result is also in agreement with source counts from the 38 MHz 8C survey, indicating that the average spectral index of radio sources flattens towards lower frequencies. We also find evidence for spectral flattening using the individual flux measurements of sources between 34 and 1400 MHz and by calculating the spectral index averaged over the source population. To select ultra-steep spectrum ($\alpha < -1.1$) radio sources, that could be associated with massive high redshift radio galaxies, we compute spectral indices between 62 MHz, 153 MHz and 1.4 GHz for sources in the Bo\"otes field. We cross-correlate these radio sources with optical and infrared catalogues and fit the spectral energy distribution to obtain photometric redshifts. We find that most of these ultra-steep spectrum sources are located in the $ 0.7 \lesssim z \lesssim 2.5$ range.

We have conducted two pilot surveys for radio pulsars and fast transients with the Low-Frequency Array (LOFAR) around 140 MHz and here report on the first low-frequency fast-radio burst limit and the discovery of two new pulsars. The first survey, the LOFAR Pilot Pulsar Survey (LPPS), observed a large fraction of the northern sky, ~1.4 x 10^4 sq. deg, with 1-hr dwell times. Each observation covered ~75 sq. deg using 7 independent fields formed by incoherently summing the high-band antenna fields. The second pilot survey, the LOFAR Tied-Array Survey (LOTAS), spanned ~600 sq. deg, with roughly a 5-fold increase in sensitivity compared with LPPS. Using a coherent sum of the 6 LOFAR "Superterp" stations, we formed 19 tied-array beams, together covering 4 sq. deg per pointing. From LPPS we derive a limit on the occurrence, at 142 MHz, of dispersed radio bursts of < 150 /day/sky, for bursts brighter than S > 107 Jy for the narrowest searched burst duration of 0.66 ms. In LPPS, we re-detected 65 previously known pulsars. LOTAS discovered two pulsars, the first with LOFAR or any digital aperture array. LOTAS also re-detected 27 previously known pulsars. These pilot studies show that LOFAR can efficiently carry out all-sky surveys for pulsars and fast transients, and they set the stage for further surveying efforts using LOFAR and the planned low-frequency component of the Square Kilometer Array.

We present radio observations of the Moon between $35$ and $80$ MHz to demonstrate a novel technique of interferometrically measuring large-scale diffuse emission extending far beyond the primary beam (global signal) for the first time. In particular, we show that (i) the Moon appears as a negative-flux source at frequencies $35<\nu<80$ MHz since it is `colder' than the diffuse Galactic background it occults, (ii) using the (negative) flux of the lunar disc, we can reconstruct the spectrum of the diffuse Galactic emission with the lunar thermal emission as a reference, and (iii) that reflected RFI (radio-frequency interference) is concentrated at the center of the lunar disc due to specular nature of reflection, and can be independently measured. Our RFI measurements show that (i) Moon-based Cosmic Dawn experiments must design for an Earth-isolation of better than $80$ dB to achieve an RFI temperature $<1$ mK, (ii) Moon-reflected RFI contributes to a dipole temperature less than $20$ mK for Earth-based Cosmic Dawn experiments, (iii) man-made satellite-reflected RFI temperature exceeds $20$ mK if the aggregate cross section of visible satellites exceeds $80$ m$^2$ at $800$ km height, or $5$ m$^2$ at $400$ km height. Currently, our diffuse background spectrum is limited by sidelobe confusion on short baselines (10-15% level). Further refinement of our technique may yield constraints on the redshifted global $21$-cm signal from Cosmic Dawn ($40>z>12$) and the Epoch of Reionization ($12>z>5$).

The Sun is an active source of radio emission which is often associated with energetic phenomena such as solar flares and coronal mass ejections (CMEs). At low radio frequencies (<100 MHz), the Sun has not been imaged extensively because of the instrumental limitations of previous radio telescopes. Here, the combined high spatial, spectral and temporal resolution of the Low Frequency Array (LOFAR) was used to study solar Type III radio bursts at 30-90 MHz and their association with CMEs. The Sun was imaged with 126 simultaneous tied-array beams within 5 solar radii of the solar centre. This method offers benefits over standard interferometric imaging since each beam produces high temporal (83 ms) and spectral resolution (12.5 kHz) dynamic spectra at an array of spatial locations centred on the Sun. LOFAR's standard interferometric output is currently limited to one image per second. Over a period of 30 minutes, multiple Type III radio bursts were observed, a number of which were found to be located at high altitudes (4 solar radii from the solar center at 30 MHz) and to have non-radial trajectories. These bursts occurred at altitudes in excess of values predicted by 1D radial electron density models. The non-radial high altitude Type III bursts were found to be associated with the expanding flank of a CME. The CME may have compressed neighbouring streamer plasma producing larger electron densities at high altitudes, while the non-radial burst trajectories can be explained by the deflection of radial magnetic fields as the CME expanded in the low corona.

This study aims to characterise the polarized foreground emission in the ELAIS-N1 field and to address its possible implications for the extraction of the cosmological 21-cm signal from the Low-Frequency Array - Epoch of Reionization (LOFAR-EoR) data. We use the high band antennas of LOFAR to image this region and RM-synthesis to unravel structures of polarized emission at high Galactic latitudes. The brightness temperature of the detected Galactic emission is on average 4 K in polarized intensity and covers the range from -10 to +13rad m^-2 in Faraday depth. The total polarized intensity and polarization angle show a wide range of morphological features. We have also used the Westerbork Synthesis Radio Telescope (WSRT) at 350 MHz to image the same region. The LOFAR and WSRT images show a similar complex morphology, at comparable brightness levels, but their spatial correlation is very low. The fractional polarization at 150 MHz, expressed as a percentage of the total intensity, amounts to 1.5%. There is no indication of diffuse emission in total intensity in the interferometric data, in line with results at higher frequencies. The wide frequency range, good angular resolution and good sensitivity make LOFAR an exquisite instrument for studying Galactic polarized emission at a resolution of 1-2 rad m^-2 in Faraday depth. The different polarised patterns observed at 150 MHz and 350 MHz are consistent with different source distributions along the line of sight wring in a variety of Faraday thin regions of emission. The presence of polarised foregrounds is a serious complication for Epoch of Reionization experiments. To avoid the leakage of polarized emission into total intensity, which can depend on frequency, we need to calibrate the instrumental polarization across the field of view to a small fraction of 1%.

The recycled pulsar PSR J2222-0137 is one of the closest known neutron stars, with a parallax distance of $267_{-0.9}^{+1.2}\,$pc and an edge-on orbit. We measure the Shapiro delay in the system through pulsar timing with the Green Bank Telescope, deriving a low pulsar mass ($1.20\pm0.14$ $M_\odot$) and a high companion mass ($1.05\pm0.06$ $M_\odot$) consistent with either a low-mass neutron star or a high-mass white dwarf. We can largely reject the neutron star hypothesis on the basis of the system's extremely low eccentricity (3e-4) - too low to have been the product of two supernovae under normal circumstances. However, despite deep optical and near-infrared searches with SOAR and the Keck telescopes we have not discovered the optical counterpart of the system. This is consistent with the white dwarf hypothesis only if the effective temperature is <3000 K, a limit that is robust to distance, mass, and atmosphere uncertainties. This would make the companion to PSR J2222-0137 one of the coolest white dwarfs ever observed. For the implied age to be consistent with the age of the Milky Way requires the white dwarf to have already crystallized and entered the faster Debye-cooling regime.

We present sub-arcsecond resolution observations of continuum emission associated with the GG Tau quadruple star system at wavelengths of 1.3, 2.8, 7.3, and 50 mm. These data confirm that the GG Tau A binary is encircled by a circumbinary ring at a radius of 235 AU with a FWHM width of ~60 AU. We find no clear evidence for a radial gradient in the spectral shape of the ring, suggesting that the particle size distribution is spatially homogeneous on angular scales of ~0.1". A central point source, likely associated with the primary component (GG Tau Aa), exhibits a composite spectrum from dust and free-free emission. Faint emission at 7.3 mm is observed toward the low-mass star GG Tau Ba, although its origin remains uncertain. Using these measurements of the resolved, multifrequency emission structure of the GG Tau A system, models of the far-infrared to radio spectrum are developed to place constraints on the grain size distribution and dust mass in the circumbinary ring. The non-negligible curvature present in the ring spectrum implies a maximum particle size of 1-10 mm, although we are unable to place strong constraints on the distribution shape. The corresponding dust mass is 30-300 M_earth, at a temperature of 20-30 K. We discuss how this significant concentration of relatively large particles in a narrow ring at a large radius might be produced in a local region of higher gas pressures (i.e., a particle "trap") located near the inner edge of the circumbinary disk.

For over a decade, the structure of the inner cavity in the transition disk of TW Hydrae has been a subject of debate. Modeling the disk with data obtained at different wavelengths has led to a variety of proposed disk structures. Rather than being inconsistent, the individual models might point to the different faces of physical processes going on in disks, such as dust growth and planet formation. Our aim is to investigate the structure of the transition disk again and to find to what extent we can reconcile apparent model differences. A large set of high-angular-resolution data was collected from near-infrared to centimeter wavelengths. We investigated the existing disk models and established a new self-consistent radiative-transfer model. A genetic fitting algorithm was used to automatize the parameter fitting. Simple disk models with a vertical inner rim and a radially homogeneous dust composition from small to large grains cannot reproduce the combined data set. Two modifications are applied to this simple disk model: (1) the inner rim is smoothed by exponentially decreasing the surface density in the inner ~3 AU, and (2) the largest grains (>100 um) are concentrated towards the inner disk region. Both properties can be linked to fundamental processes that determine the evolution of protoplanetary disks: the shaping by a possible companion and the different regimes of dust-grain growth, respectively. The full interferometric data set from near-infrared to centimeter wavelengths requires a revision of existing models for the TW Hya disk. We present a new model that incorporates the characteristic structures of previous models but deviates in two key aspects: it does not have a sharp edge at 4 AU, and the surface density of large grains differs from that of smaller grains. This is the first successful radiative-transfer-based model for a full set of interferometric data.

RadioAstron space-ground VLBI observations of the pulsar B0950+08, conducted with the 10-m space radio telescope in conjunction with the Arecibo 300-m telescope and Westerbork Synthesis Radio Telescope at a frequency of 324 MHz, were analyzed in order to investigate plasma inhomogeneities in the direction of this nearby pulsar. The observations were conducted at a spacecraft distance of 330,000 km, resulting in a projected baseline of 220,000 km, providing the greatest angular resolution ever achieved at meter wavelengths. Our analysis is based on fundamental behavior of structure and coherence functions. We find that the pulsar shows scintillation on two frequency scales, both much less than the observing frequency; but modulation is less than 100%. We infer that the scattering is weak, but a refracting wedge disperses the scintillation pattern. The refraction angle of this "cosmic prism" is measured as theta_0=1.1 - 4.4 mas, with the refraction direction being approximately perpendicular to the observer velocity. We show that the observed parameters of scintillation effects indicate that two plasma layers lie along the line of sight to the pulsar, at distances of 4.4 - 16.4 pc and 26 - 170 pc, and traveling in different directions relative to the line of sight. Spectra of turbulence for the two layers are found to follow a power law with the indices gamma_1 = gamma_2 = 3.00 +/- 0.08, significantly different from the index expected for a Kolmogorov spectrum of turbulence, gamma=11/3.

A patch of sky in the SDSS Stripe 82 was observed at 1.6 GHz with Very Long Baseline Interferometry (VLBI) using the European VLBI Network (EVN). The data were correlated at the EVN software correlator at JIVE (SFXC). There are fifteen known mJy/sub-mJy radio sources in the target field defined by the primary beam size of a typical 30-m class EVN radio telescope. The source of particular interest is a recently identified high-redshift radio quasar J222843.54+011032.2 (J2228+0110) at redshift z=5.95. Our aim was to investigate the milli-arcsecond (mas) scale properties of all the VLBI-detectable sources within this primary beam area with a diameter of 20 arcmin. The source J2228+0110 was detected with VLBI with a brightness temperature T_b>10^8 K, supporting the active galactic nucleus (AGN) origin of its radio emission, which is conclusive evidence that the source is a radio quasar. In addition, two other target sources were also detected, one of them with no redshift information. Their brightness temperature values (T_b >10^7 K) measured with VLBI suggest a non-thermal synchrotron radiation origin for their radio emission. The detection rate of 20% is broadly consistent with other wide-field VLBI experiments carried out recently. We also derived the accurate equatorial coordinates of the three detected sources using the phase-referencing technique. This experiment is an early attempt of a wide-field science project with SFXC, paving the way for the EVN to conduct a large-scale VLBI survey in the multiple-phase-centre mode.

Gravitationally bound three-body systems have been studied for hundreds of years and are common in our Galaxy. They show complex orbital interactions, which can constrain the compositions, masses, and interior structures of the bodies and test theories of gravity, if sufficiently precise measurements are available. A triple system containing a radio pulsar could provide such measurements, but the only previously known such system, B1620-26 (with a millisecond pulsar, a white dwarf, and a planetary-mass object in an orbit of several decades), shows only weak interactions. Here we report precision timing and multi-wavelength observations of PSR J0337+1715, a millisecond pulsar in a hierarchical triple system with two other stars. Strong gravitational interactions are apparent and provide the masses of the pulsar (1.4378(13) Msun, where Msun is the solar mass and the parentheses contain the uncertainty in the final decimal places) and the two white dwarf companions (0.19751(15) Msun and 0.4101(3) Msun), as well as the inclinations of the orbits (both approximately 39.2 degrees). The unexpectedly coplanar and nearly circular orbits indicate a complex and exotic evolutionary past that differs from those of known stellar systems. The gravitational field of the outer white dwarf strongly accelerates the inner binary containing the neutron star, and the system will thus provide an ideal laboratory in which to test the strong equivalence principle of general relativity.

We present the description and early results of the mJy Imaging VLBA Exploration at 20 cm (mJIVE-20). mJIVE-20 is a large project on the Very Long Baseline Array which uses "filler" time to systematically inspect a large sample of mJy radio sources, pre-selected from the FIRST survey made with the Very Large Array, to identify any compact emission which may be present. The newly available multifield capability of the VLBA makes it possible for us to inspect of order 100 sources per hour of observing time with a 6.75 sigma detection sensitivity of approximately 1 mJy/beam. The results of the mJIVE-20 survey are made publicly available as soon as the data are calibrated. After 18 months of observing, over 20,000 FIRST sources have been inspected, with 4,336 VLBI detections. These initial results suggest that within the range 1 - 200 mJy, fainter sources are somewhat more likely to be dominated by a very compact component than brighter sources. Over half of all arcsecond-scale mJy radio sources contain a compact component, although the fraction of sources which are dominated by milliarcsecond scale structure is smaller at 30-35%, increasing towards lower flux densities. Significant differences are seen depending on the optical classification of the source. Radio sources with a stellar/point-like counterpart in the Sloan Digital Sky Survey (SDSS) are more likely to be detected overall, but this detection likelihood appears to be independent of the arcsecond-scale radio flux density; the trend towards higher radio compactness for fainter sources is confined to SDSS galaxies and non-detections. These results are consistent with a unification model of AGN in which less luminous sources have on average slower radio jets, with lower Doppler suppression of compact core emission over a wider range of viewing angles. (abridged)

The pulsed radio emission from the Galactic Center (GC) magnetar SGR J1745-29 probes the turbulent, magnetized plasma of the GC hyperstrong scattering screen through both angular and temporal broadening. We present the first measurements of the angular size of SGR J1745-29, obtained with the Very Long Baseline Array and the phased Very Large Array at 8.7 and 15.4 GHz. The source sizes are consistent with the scatter--broadened size of Sagittarius A* at each frequency, demonstrating that SGR J1745-29 is also located behind the same hyperstrong scattering medium. Combining the angular broadening with temporal scattering obtained from pulsar observations provides a complete picture of the scattering properties. If the scattering occurs in a thin screen, then it must be at a distance \Delta >~ 5 kpc. A best-fit solution for the distance of a thin screen is \Delta=5.9 +/-0.3 kpc, consistent with being located in the Scutum spiral arm. This is a substantial revision of the previously held model in which the scattering screen is located very close to the GC. As also discussed in Spitler et al., these results suggest that GC searches can detect millisecond pulsars gravitationally bound to Sgr A* with observations at >~ 10 GHz and ordinary pulsars at even lower frequencies.

The centre of our Milky Way harbours the closest candidate for a supermassive black hole. The source is thought to be powered by radiatively inefficient accretion of gas from its environment. This form of accretion is a standard mode of energy supply for most galactic nuclei. X-ray measurements have already resolved a tenuous hot gas component from which it can be fed. However, the magnetization of the gas, a crucial parameter determining the structure of the accretion flow, remains unknown. Strong magnetic fields can influence the dynamics of the accretion, remove angular momentum from the infalling gas, expel matter through relativistic jets and lead to the observed synchrotron emission. Here we report multi-frequency measurements with several radio telescopes of a newly discovered pulsar close to the Galactic Centre and show that its unusually large Faraday rotation indicates a dynamically relevant magnetic field near the black hole. If this field is accreted down to the event horizon it provides enough magnetic flux to explain the observed emission from the black hole, from radio to X-rays.

The binary pulsar J2222-0137 is an enigmatic system containing a partially recycled millisecond pulsar and a companion of unknown nature. Whilst the low eccentricity of the system favors a white dwarf companion, an unusual double neutron star system is also a possibility, and optical observations will be able to distinguish between these possibilities. In order to allow the absolute luminosity (or upper limit) of the companion object to be properly calibrated, we undertook astrometric observations with the Very Long Baseline Array to constrain the system distance via a measurement of annual geometric parallax. With these observations, we measure the parallax of the J2222-0137 system to be 3.742 +0.013 -0.016 milliarcseconds, yielding a distance of 267.3 +1.2 -0.9 pc, and measure the transverse velocity to be 57.1 +0.3 -0.2 km/s. Fixing these parameters in the pulsar timing model made it possible to obtain a measurement of Shapiro delay and hence the system inclination, which shows that the system is nearly edge-on (sin i = 0.9985 +/- 0.0005). Furthermore, we were able to detect the orbital motion of J2222-0137 in our VLBI observations and measure the longitude of ascending node. The VLBI astrometry yields the most accurate distance obtained for a radio pulsar to date, and is furthermore the most accurate parallax for any radio source obtained at "low" radio frequencies (below ~5 GHz, where the ionosphere dominates the error budget). Using the astrometric results, we show the companion to J2222-0137 will be easily detectable in deep optical observations if it is a white dwarf. Finally, we discuss the implications of this measurement for future ultra-high-precision astrometry, in particular in support of pulsar timing arrays.

LOFAR, the LOw-Frequency ARray, is a new-generation radio interferometer constructed in the north of the Netherlands and across europe. Utilizing a novel phased-array design, LOFAR covers the largely unexplored low-frequency range from 10-240 MHz and provides a number of unique observing capabilities. Spreading out from a core located near the village of Exloo in the northeast of the Netherlands, a total of 40 LOFAR stations are nearing completion. A further five stations have been deployed throughout Germany, and one station has been built in each of France, Sweden, and the UK. Digital beam-forming techniques make the LOFAR system agile and allow for rapid repointing of the telescope as well as the potential for multiple simultaneous observations. With its dense core array and long interferometric baselines, LOFAR achieves unparalleled sensitivity and angular resolution in the low-frequency radio regime. The LOFAR facilities are jointly operated by the International LOFAR Telescope (ILT) foundation, as an observatory open to the global astronomical community. LOFAR is one of the first radio observatories to feature automated processing pipelines to deliver fully calibrated science products to its user community. LOFAR's new capabilities, techniques and modus operandi make it an important pathfinder for the Square Kilometre Array (SKA). We give an overview of the LOFAR instrument, its major hardware and software components, and the core science objectives that have driven its design. In addition, we present a selection of new results from the commissioning phase of this new radio observatory.

The VLBA was used at 1.6 GHz to observe a target field 50 arcminutes in diameter including the core of M31. Novel VLBI correlation techniques were used to observe 200 sources simultaneously, of which 16 were detected. We classify all 16 as background AGN based on their X-ray properties and arcsecond- and mas-scale morphology. The detected sources were then analyzed for evidence of scatter-broadening due to the ionized ISM of M31. The detection of a compact background source only 0.25 kpc projected distance from M31 places a constraint on the extent of any extreme scattering region associated with center of M31. However, the two sources closest to the core show evidence of scatter broadening consistent with that which would be observed if a compact source were observed through the inner disk of our Galaxy at the inclination of M31. We interpret this as a detection of the ionized ISM of M31 along two lines of sight. With the increases in bandwidth and sensitivity envisaged for future long-baseline interferometers, this should prove to be a remarkably powerful technique for understanding the ionized ISM in external galaxies.

For Very Long Baseline Interferometry (VLBI), the fringe spacing is extremely narrow compared to the field of view imposed by the primary beam of each element. This means that an extremely large number of resolution units can potentially be imaged from a single observation. We implement and test a technique for efficiently and accurately imaging large VLBI datasets. The DiFX software correlator is used to generate a dataset with extremely high time and frequency resolution. This large dataset is then transformed and averaged multiple times to generate many smaller datasets, each with a phase centre located at a different area of interest. Results of an 8.4 GHz four-station VLBI observation of a field containing multiple sources are presented. Observations of the calibrator 3C345 were used for preliminary tests of accuracy of the shifting algorithm. A high level of accuracy was achieved, making the method suitable even for the most demanding astrometric VLBI observations. One target source (1320+299A) was detected and was used as a phase-reference calibrator in searching for further detections. An image containing 13 billion pixels was constructed by independently imaging 782 visibility datasets covering the entire primary beam of the array. Current implementations of this algorithm and possible future developments in VLBI data analysis are discussed.

Aims: Very Long Baseline Interferometry yields physical constraints on the compact radio emission of young stellar objects. At the same time, such measurements can be used for precise astrometric measurements of parallaxes and proper motions. Here, we aimed to make the first detections of very compact radio emission from class I protostars in the Corona Australis star-forming region. Methods: We have used the Long Baseline Array (LBA) to observe the protostars IRS 5 and IRS 7 in the Corona Australis star-forming region in three separate epochs Results: We report the first firm radio detection of a class I protostar using very long baseline interferometry. We have detected the previously known non-thermal source IRS 5b, part of a binary system. IRS 5a and IRS 5N were undetected, as were all sources in the IRS 7 region. Conclusions: These results underline the unusual nature of IRS 5b as a genuine protostar with confirmed non-thermal radio emission. Also, these observations highlight the potential of the LBA as a tool to provide precision astrometric measurements of individual young stellar objects in southern star-forming regions that are not accessible to the Very Long Baseline Array in the northern hemisphere.

We use VLBI imaging of the interstellar scattering speckle pattern associated with the pulsar PSR 0834+06 to measure the astrometric motion of its emission. The ~ 5AU interstellar baselines, provided by interference between speckles spanning the scattering disk, enable us to detect motions with sub nanoarcsecond accuracy. We measure a small pulse deflection of ~18+/-2 km (not including geometric uncertainties), which is 100 times smaller than the native resolution of this interstellar interferometer. This implies that the emission region is small, and at an altitude of a few hundred km, with the exact value depending on field geometry. This is substantially closer to the star than to the light cylinder. Future VLBI measurements can improve on this finding. This new regime of ultra-precise astrometry may enable precision parallax distance determination of pulsar binary displacements.

A new generation of observational science instruments is dramatically increasing collected data volumes in a range of fields. These instruments include the Square Kilometre Array (SKA), Large Synoptic Survey Telescope (LSST), terrestrial sensor networks, and NASA satellites participating in "decadal survey" missions. Their unprecedented coverage and sensitivity will likely reveal wholly new categories of unexpected and transient events. Commensal methods passively analyze these data streams, recognizing anomalous events of scientific interest and reacting in real time. We report on a case example: V-FASTR, an ongoing commensal experiment at the Very Long Baseline Array (VLBA) that uses online adaptive pattern recognition to search for anomalous fast radio transients. V-FASTR triages a millisecond-resolution stream of data and promotes candidate anomalies for further offline analysis. It tunes detection parameters in real time, injecting synthetic events to continually retrain itself for optimum performance. This self-tuning approach retains sensitivity to weak signals while adapting to changing instrument configurations and noise conditions. The system has operated since July 2011, making it the longest-running real time commensal radio transient experiment to date.

We define a framework for determining constraints on the detection rate of fast transient events from a population of underlying sources, with a view to incorporating beam shape, frequency effects, scattering effects, and detection efficiency into the metric. We then demonstrate a method for combining independent datasets into a single event rate constraint diagram, using a probabilistic approach to the limits on parameter space. We apply this new framework to present the latest results from the V-FASTR experiment, a commensal fast transients search using the Very Long Baseline Array (VLBA). In the 20 cm band, V-FASTR now has the ability to probe the regions of parameter space of importance for the observed Lorimer and Keane fast radio transient candidates, by combining the information from observations with differing bandwidths, and properly accounting for the source dispersion measure, VLBA antenna beam shape, experiment time sampling, and stochastic nature of events. We then apply the framework to combine the results of the V-FASTR and ATA Fly's Eye experiments, demonstrating their complementarity. Expectations for fast transients experiments for the SKA Phase I dish array are then computed, and the impact of large differential bandwidths is discussed.

Active Galactic Nuclei (AGN) play a decisive role in galaxy evolution, particularly so when operating in a radiatively inefficient mode, where they launch powerful jets that reshape their surroundings. However, identifying them is difficult, since radio observations commonly have resolutions of between 1 arcsec and 10 arcsec, which is equally sensitive to radio emission from star-forming activity and from AGN. Very Long Baseline Interferometry (VLBI) observations allow one to filter out all but the most compact non-thermal emission from radio survey data. The observational and computational demands to do this in large surveys have been, until recently, too high to make such undertakings feasible. Only the recent advent of wide-field observing techniques have facilitated such observations, and we here present the results from a survey of 217 radio sources in the Lockman Hole/XMM field. We describe in detail some new aspects of the calibration, including primary beam correction, multi-source self-calibration, and mosaicing. As a result, we detected 65 out of the 217 radio sources and were able to construct, for the first time, the source counts of VLBI-detected AGN. They indicate that at least 15%-25% of the sub-mJy radio sources are AGN-driven, consistent with recent findings using other AGN selection techniques. We have used ancillary data to investigate the AGN hosts. We find that among the sources nearby enough to be resolved in the optical images, 88% (23/26) could be classified as early-type or bulge-dominated galaxies. While 50% of these sources are correctly represented by the SED of an early-type galaxy, for the rest the best fit was obtained with a heavily extinct starburst template, an effect we ascribe to a degeneracy in the fit. Overall, the typical hosts of VLBI-detected sources are in good agreement with being early-type or bulge-dominated galaxies.

We present dust continuum observations of the protoplanetary disk surrounding the pre-main sequence star AS 209, spanning more than an order of magnitude in wavelength from 0.88 to 9.8 mm. The disk was observed with sub-arcsecond angular resolution (0.2"-0.5") to investigate radial variations in its dust properties. At longer wavelengths, the disk emission structure is notably more compact, providing model-independent evidence for changes in the grain properties across the disk. We find that physical models which reproduce the disk emission require a radial dependence of the dust opacity \kappa_{\nu}. Assuming that the observed wavelength-dependent structure can be attributed to radial variations in the dust opacity spectral index (\beta), we find that \beta(R) increases from \beta<0.5 at \sim20 AU to \beta>1.5 for R>80 AU, inconsistent with a constant value of \beta\ across the disk (at the 10\sigma\ level). Furthermore, if radial variations of \kappa_{\nu} are caused by particle growth, we find that the maximum size of the particle-size distribution (a_{max}) increases from sub-millimeter-sized grains in the outer disk (R>70 AU) to millimeter and centimeter-sized grains in the inner disk regions (R< 70 AU). We compare our observational constraint on a_{max}(R) with predictions from physical models of dust evolution in proto-planetary disks. For the dust composition and particle-size distribution investigated here, our observational constraints on a_{max}(R) are consistent with models where the maximum grain size is limited by radial drift.

The recently discovered transitional millisecond pulsar system J1023+0038 exposes a crucial evolutionary phase of recycled neutron stars for multiwavelength study. The system, comprising the neutron star itself, its stellar companion, and the surrounding medium, is visible across the electromagnetic spectrum from the radio to X-ray/gamma-ray regimes and offers insight into the recycling phase of millisecond pulsar evolution. Here, we report on multiple-epoch astrometric observations with the Very Long Baseline Array (VLBA) which give a system parallax of 0.731 +/- 0.022 milliarcseconds (mas) and a proper motion of 17.98 +/- 0.05 mas/yr. By combining our results with previous optical observations, we are able to use the parallax distance of 1368+42-39 pc to estimate the mass of the pulsar as 1.71 +/- 0.16 solar masses, and we are also able to measure the 3D space velocity of the system as 126 +/- 5 km/s. Despite the precise nature of the VLBA measurements, the remaining ~3% distance uncertainty dominates the 0.16 solar mass error on our mass estimate.

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.

We present the first results from the V-FASTR experiment, a commensal search for fast transient radio bursts using the Very Long Baseline Array (VLBA). V-FASTR is unique in that the widely spaced VLBA antennas provide a discriminant against non-astronomical signals and a mechanism for the localization and identification of events that is not possible with single dishes or short baseline interferometers. Thus far V-FASTR has accumulated over 1300 hours of observation time with the VLBA, between 90 cm and 3 mm wavelength (327 MHz - 86 GHz), providing the first limits on fast transient event rates at high radio frequencies (>1.4 GHz). V-FASTR has blindly detected bright individual pulses from seven known pulsars but has not detected any single-pulse events that would indicate high redshift impulsive bursts of radio emission. At 1.4 GHz, V-FASTR puts limits on fast transient event rates comparable with the PALFA survey at the Arecibo telescope, but generally at lower sensitivities, and comparable to the "fly's eye" survey at the Allen Telescope Array, but with less sky coverage. We also illustrate the likely performance of the Phase 1 SKA dish array for an incoherent fast transient search fashioned on V-FASTR.

A candidate accreting massive black hole (BH) with M_BH ~ 10^6 Msun has recently been identified at the center of the dwarf starburst galaxy Henize 2-10 (He 2-10). This discovery offers the first possibility of studying a growing BH in a nearby galaxy resembling those in the earlier universe, and opens up a new class of host galaxies to search for the smallest supermassive BHs. Here we present very long baseline interferometry observations of He 2-10 taken with the Long Baseline Array (LBA) at 1.4 GHz with an angular resolution of ~0.1" x 0.03". A single compact radio source is detected at the precise location of the putative low-luminosity active galactic nucleus. The physical size of the nuclear radio emission is < 3 pc x 1 pc, an order of magnitude smaller than previous constraints from the Very Large Array (VLA), and the brightness temperature of T_B > 3 x 10^5 K confirms a non-thermal origin. These LBA observations indicate that the nuclear radio emission originates from a single object and exclude the possibility of multiple supernova remnants as the origin of the nuclear radio emission previously detected with the VLA at lower resolution. A weaker, more extended, off-nuclear source is also detected with the LBA and a comparison with multi-wavelength ancillary data indicate that, unlike the nuclear source, the off-nuclear source is co-spatial with a super star cluster, lacks a detectable X-ray point-source counterpart, and is almost certainly due to a supernova remnant in the host star cluster.

The formation mechanism of neutron stars with extremely large magnetic field strengths (magnetars) remains unclear. Some formation scenarios predict that magnetars should be born with extremely high space velocities, >1000 km/s. Using the Long Baseline Array in Australia, we have measured the proper motion of the intermittently radio-bright magnetar J1550-5418 (1E 1547.0-5408) to be 9.2 +/- 0.6 mas/yr. For a likely distance of 6 +/- 2 kpc, the implied transverse velocity is 280 + 130 - 120 km/s after correcting for Galactic rotation. Along with the ~200 km/s transverse velocity measured for the magnetar XTE J1810-197, this result suggests that formation pathways producing large magnetic fields do not require very large birth kicks.

Obtaining pulsar parallaxes via relative astrometry (also known as differential astrometry) yields distances and transverse velocities that can be used to probe properties of the pulsar population and the interstellar medium. Large programs are essential to obtain the sample sizes necessary for these population studies, but they must be efficiently conducted to avoid requiring an infeasible amount of observing time. This paper describes the PSRPI astrometric program, including the use of new features in the DiFX software correlator to efficiently locate calibrator sources, selection and observing strategies for a sample of 60 pulsars, initial results, and likely science outcomes. Potential applications of high-precision relative astrometry to measure source structure evolution in defining sources of the International Celestial Referent Frame are also discussed.

PSR B1259-63 is a 48 ms pulsar in a highly eccentric 3.4 year orbit around the young massive star LS 2883. During the periastron passage the system displays transient non-thermal unpulsed emission from radio to very high energy gamma rays. It is one of the three galactic binary systems clearly detected at TeV energies, together with LS 5039 and LS I +61 303. We observed PSR B1259-63 after the 2007 periastron passage with the Australian Long Baseline Array at 2.3 GHz to trace the milliarcsecond (mas) structure of the source at three different epochs. We have discovered extended and variable radio structure. The peak of the radio emission is detected outside the binary system near periastron, at projected distances of 10-20 mas (25-45 AU assuming a distance of 2.3 kpc). The total extent of the emission is ~50 mas (~120 AU). This is the first observational evidence that non-accreting pulsars orbiting massive stars can produce variable extended radio emission at AU scales. Similar structures are also seen in LS 5039 and LS I +61 303, in which the nature of the compact object is unknown. The discovery presented here for the young non-accreting pulsar PSR B1259-63 reinforces the link with these two sources and supports the presence of pulsars in these systems as well. A simple kinematical model considering only a spherical stellar wind can approximately trace the extended structures if the binary system orbit has a longitude of the ascending node of omega ~ -40 deg and a magnetization parameter of sigma~0.005.

Recent discoveries of dispersed, non-periodic impulsive radio signals with single-dish radio telescopes have sparked significant interest in exploring the relatively uncharted space of fast transient radio signals. Here we describe V-FASTR, an experiment to perform a blind search for fast transient radio signals using the Very Long Baseline Array (VLBA). The experiment runs entirely in a commensal mode, alongside normal VLBA observations and operations. It is made possible by the features and flexibility of the DiFX software correlator that is used to process VLBA data. Using the VLBA for this type of experiment offers significant advantages over single-dish experiments, including a larger field of view, the ability to easily distinguish local radio-frequency interference from real signals and the possibility to localize detected events on the sky to milliarcsecond accuracy. We describe our software pipeline, which accepts short integration (~ms) spectrometer data from each antenna in real time during correlation and performs an incoherent dedispersion separately for each antenna, over a range of trial dispersion measures. The dedispersed data are processed by a sophisticated detector and candidate events are recorded. At the end of the correlation, small snippets of the raw data at the time of the events are stored for further analysis. We present the results of our event detection pipeline from some test observations of the pulsars B0329+54 and B0531+21 (the Crab pulsar).

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.

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.

In-correlator techniques offer the possibility of identifying and/or excising radio frequency interference (RFI) from interferometric observations at much higher time and/or frequency resolution than is generally possible with the final visibility dataset. Due to the considerable computational requirements of the correlation procedure, cross-correlators have most commonly been implemented using high-speed digital signal processing boards, which typically require long development times and are difficult to alter once complete. "Software" correlators, on the other hand, make use of commodity server machines and a correlation algorithm coded in a high-level language. They are inherently much more flexible and can be developed - and modified - much more rapidly than purpose-built "hardware" correlators. Software correlators are thus a natural choice for testing new RFI detection and mitigation techniques for interferometers. The ease with which software correlators can be adapted to test RFI detection algorithms is demonstrated by the addition of kurtosis detection and plotting to the widely used DiFX software correlator, which highlights previously unknown short -duration RFI at the Hancock VLBA station.

VLBI observations are a reliable method to identify AGN, since they require high brightness temperatures for a detection to be made. However, because of the tiny fields of view it is unpractical to carry out VLBI observations of many sources using conventional methods. We used an extension of the DiFX software correlator to image with high sensitivity 96 sources in the Chandra Deep Field South, using only 9h of observing time with the VLBA. We detected 20 sources, 8 of which had not been identified as AGN at any other wavelength, despite the comprehensive coverage of this field. The lack of X-ray counterparts to 1/3 of the VLBI-detected sources, despite the sensitivity of co-located X-ray data, demonstrates that X-ray observations cannot be solely relied upon when searching for AGN activity. Surprisingly, we find that sources classified as type 1 QSOs using X-ray data are always detected, in contrast to the 10% radio-loud objects which are found in optically-selected QSOs. We present the continuation of this project with the goal to image 1450 sources in the Lockman Hole/XMM region.

Wide-field surveys are a commonly-used method for studying thousands of objects simultaneously, to investigate, e.g., the joint evolution of star-forming galaxies and active galactic nuclei. VLBI observations can yield valuable input to such studies because they are able to identify AGN. However, VLBI observations of large swaths of the sky are impractical using standard methods, because the fields of view of VLBI observations are of the order of 10" or less. We have embarked on a project to carry out Very Long Baseline Array (VLBA) observations of all 96 known radio sources in one of the best-studied areas in the sky, the Chandra Deep Field South (CDFS). The challenge was to develop methods which could significantly reduce the amount of observing (and post-processing) time. We have developed an extension to the DiFX software correlator which allows one to correlate hundreds of positions within the primary beams. This extension enabled us to target many sources, at full resolution and high sensitivity, using only a small amount of observing time. The combination of wide fields-of-view and high sensitivity across the field in this survey is unprecedented. We have observed with the VLBA a single pointing containing the Chandra Deep Field South, in which 96 radio sources were known from previous observations with the ATCA. From our input sample, 20 were detected with the VLBA. The majority of objects have flux densities in agreement with arcsec-scale observations, implying that their radio emission comes from very small regions. One VLBI-detected object had earlier been classified as a star-forming galaxy. Comparing the VLBI detections to sources found in sensitive, co-located X-ray observations we find that X-ray detections are not a good indicator for VLBI detections. Wide-field VLBI survey science is now coming of age.

We are developing a purely commensal survey experiment for fast (<5s) transient radio sources. Short-timescale transients are associated with the most energetic and brightest single events in the Universe. Our objective is to cover the enormous volume of transients parameter space made available by ASKAP, with an unprecedented combination of sensitivity and field of view. Fast timescale transients open new vistas on the physics of high brightness temperature emission, extreme states of matter and the physics of strong gravitational fields. In addition, the detection of extragalactic objects affords us an entirely new and extremely sensitive probe on the huge reservoir of baryons present in the IGM. We outline here our approach to the considerable challenge involved in detecting fast transients, particularly the development of hardware fast enough to dedisperse and search the ASKAP data stream at or near real-time rates. Through CRAFT, ASKAP will provide the testbed of many of the key technologies and survey modes proposed for high time resolution science with the SKA.

We have invented a novel technique to measure the radio image of a pulsar scattered by the interstellar plasma with 0.1 mas resolution. We extend the "secondary spectrum" analysis of parabolic arcs by Stinebring et al. (2001) to very long baseline interferometry and, when the scattering is anisotropic, we are able to map the scattered brightness astrometrically with much higher resolution than the diffractive limit of the interferometer. We employ this technique to measure an extremely anisotropic scattered image of the pulsar B0834+06 at 327 MHz. We find that the scattering occurs in a compact region about 420 pc from the Earth. This image has two components, both essentially linear and nearly parallel. The primary feature, which is about 16 AU long and less than 0.5 AU in width, is highly inhomogeneous on spatial scales as small as 0.05 AU. The second feature is much fainter and is displaced from the axis of the primary feature by about 9 AU. We find that the velocity of the scattering plasma is 16+-10 km/s approximately parallel to the axis of the linear feature. The origin of the observed anisotropy is unclear and we discuss two very different models. It could be, as has been assumed in earlier work, that the turbulence on spatial scales of ~1000 km is homogeneous but anisotropic. However it may be that the turbulence on these scales is homogeneous and isotropic but the anisotropy is produced by highly elongated (filamentary) inhomogeneities of scale 0.05-16 AU.

Accurate measurement of pulsar distances via astrometry using very long baseline interferometry enables the improvement of Galactic electron density distribution models, improving distance estimates for the vast majority of pulsars for which parallax measurements are unavailable. However, pulsars at southern declinations have been under-represented in previous interferometric astrometry campaigns. In order to redress this imbalance, we have conducted a two-year astrometric campaign targeting eight southern pulsars with the Australian Long Baseline Array. The program summarized in this paper has resulted in the measurement of seven new pulsar parallaxes, with success on objects down to a mean flux density of 0.8 mJy at 1600 MHz. Our results highlight the substantial uncertainties that remain when utilizing free electron density models for individual pulsar distances. Until this study, PSR J0630-2834 was believed to convert 16% of its spin-down energy into x-rays, but our measured parallax distance of 332 (+52 -40) pc has revised this value to <1%. In contrast, PSR J0108-1431 was found to be almost a factor of two more distant than previously thought, making its conversion of spin-down energy to x-rays the most efficient known (>1%). The 8.5 second radio pulsar J2144-3933 was found to be closer than previously predicted, making its apparent 1400 MHz radio luminosity the lowest of any known pulsar (20 microJy kpc^2). We have examined the growing population of neutron stars with accurate parallaxes to determine the effect of distance errors on the underlying neutron star velocity distribution, and find that typical distance errors may be biasing the estimated mean pulsar velocity upwards by 5%, and are likely to exaggerate the distribution's high-velocity tail.

GHz-peaked spectrum (GPS) radio sources are thought to be young objects which later evolve into FR-I and FR-II radio galaxies. We have used the Australia Telescope 20GHz (AT20G) survey catalogue to select a uniform sample of GPS sources with spectral peaks above 5GHz, which should represent the youngest members of this class. In this paper, we present e-VLBI observations of ten such objects which are associated with nearby (z<0.15) galaxies and so represent a new population of local, low--power GPS sources. Our e-VLBI observations were carried out at 4.8GHz with the Australia Telescope Long Baseline Array (LBA) using a real--time software correlator. All ten sources were detected, and were unresolved on scales of ~100mas, implying that they are typically less than 100pc in linear size.

Time-to-science is an important figure of merit for digital instrumentation serving the astronomical community. A digital signal processing (DSP) community is forming that uses shared hardware development, signal processing libraries, and instrument architectures to reduce development time of digital instrumentation and to improve time-to-science for a wide variety of projects. We suggest prioritizing technological development supporting the needs of this nascent DSP community. After outlining several instrument classes that are relying on digital instrumentation development to achieve new science objectives, we identify key areas where technologies pertaining to interoperability and processing flexibility will reduce the time, risk, and cost of developing the digital instrumentation for radio astronomy. These areas represent focus points where support of general-purpose, open-source development for a DSP community should be prioritized in the next decade. Contributors to such technological development may be centers of support for this DSP community, science groups that contribute general-purpose DSP solutions as part of their own instrumentation needs, or engineering groups engaging in research that may be applied to next-generation DSP instrumentation.

The double pulsar J0737-3039A/B is a unique system with which to test gravitational theories in the strong-field regime. However, the accuracy of such tests will be limited by knowledge of the distance and relative motion of the system. Here we present very long baseline interferometry observations which reveal that the distance to PSR J0737-3039A/B is 1150+220-160 pc, more than double previous estimates, and confirm its low transverse velocity (~9 km/s). Combined with a decade of pulsar timing, these results will allow tests of gravitational radiation emission theories at the 0.01% level, putting stringent constraints on theories which predict dipolar gravitational radiation. They also allow insight into the system's formation and the source of its high-energy emission.

(Abridged) This thesis describes the development of DiFX, the first general-purpose software correlator for radio interferometry, and its use with the Australian Long Baseline Array (LBA) to complete the largest VLBI pulsar astrometry program undertaken to date in the Southern Hemisphere. This two year astrometry program has resulted in the measurement of seven new pulsar parallaxes, more than trebling the number of measured VLBI pulsar parallaxes in the Southern Hemisphere. The measurements included a determination of the distance and transverse velocity of PSR J0437-4715 with better than 1% accuracy, enabling improved tests of General Relativity, and the first significant measurement of parallax for the famous double pulsar system PSR J0737-3039A/B, which will allow tests of General Relativity in this system to proceed to the 0.01% level. The DiFX software correlator developed to enable this science has been extensively tested and is now an integral part of the upgraded LBA Major National Research Facility; furthermore, it has been selected to facilitate a substantial sensitivity upgrade for the US Very Long Baseline Array.

[ABRIDGED VERSION] The future of cm and m-wave astronomy lies with the Square Kilometre Array (SKA), a telescope under development by a consortium of 17 countries. The SKA will be 50 times more sensitive than any existing radio facility. A majority of the key science for the SKA will be addressed through large-area imaging of the Universe at frequencies from 300 MHz to a few GHz. The Australian SKA Pathfinder (ASKAP) is aimed squarely in this frequency range, and achieves instantaneous wide-area imaging through the development and deployment of phase-array feed systems on parabolic reflectors. This large field-of-view makes ASKAP an unprecedented synoptic telescope poised to achieve substantial advances in SKA key science. The central core of ASKAP will be located at the Murchison Radio Observatory in inland Western Australia, one of the most radio-quiet locations on the Earth and one of the sites selected by the international community as a potential location for the SKA. Following an introductory description of ASKAP, this document contains 7 chapters describing specific science programmes for ASKAP. The combination of location, technological innovation and scientific program will ensure that ASKAP will be a world-leading radio astronomy facility, closely aligned with the scientific and technical direction of the SKA. A brief summary chapter emphasizes the point, and considers discovery space.

We describe a data reduction pipeline for VLBI astrometric observations of pulsars, implemented using the ParselTongue AIPS interface. The pipeline performs calibration (including ionosphere modeling), phase referencing with proper accounting of reference source structure, amplitude corrections for pulsar scintillation, and position fitting to yield the position, proper motion and parallax. The optimal data weighting scheme to minimize the total error budget of a parallax fit, and how this scheme varies with pulsar parameters such as flux density, is also investigated. The robustness of the techniques employed are demonstrated with the presentation of the first results from a two year astrometry program using the Australian Long Baseline Array (LBA). The parallax of PSR J1559-4438 is determined to be 0.384 +- 0.081 mas (1 sigma), resulting in a distance estimate of 2600 pc which is consistent with earlier DM and HI absorption estimates.

Using the recently upgraded Long Baseline Array, we have measured the trigonometric parallax of PSR J0437-4715 to better than 1% precision, the most precise pulsar distance determination made to date. Comparing this VLBI distance measurement to the kinematic distance obtained from pulsar timing, which is calculated from the pulsar's proper motion and apparent rate of change of orbital period, gives a precise limit on the unmodeled relative acceleration between the Solar System and PSR J0437-4715, which can be used in a variety of applications. Firstly, it shows that Newton's gravitational constant G is stable with time (\dot{G}/G = (-5 +- 26) x 10^{-13} yr^{-1}, 95% confidence). Secondly, if a stochastic gravitational wave background existed at the currently quoted limit, this null result would fail ~50% of the time. Thirdly, it excludes Jupiter-mass planets within 226 AU of the Sun in 50% of the sky (95% confidence). Finally, the ~1% agreement of the parallax and orbital period derivative distances provides a fundamental confirmation of the parallax distance method upon which all astronomical distances are based.

Infrared-Faint Radio Sources represent a new and unexpected class of object which is bright at radio wavelengths but unusually faint at infrared wavelengths. If, like most mJy radio sources, they were either conventional active or star-forming galaxies in the local Universe, we would expect them to be detectable at infrared wavelengths, and so their non-detection by the Spitzer Space Telescope is surprising. Here we report the detection of one of these sources using Very Long Baseline Interferometry, from which we conclude that the sources are driven by Active Galactic Nuclei. We suggest that these sources are either normal radio-loud quasars at high redshift or abnormally obscured radio galaxies.

We describe the development of an FX style correlator for Very Long Baseline Interferometry (VLBI), implemented in software and intended to run in multi-processor computing environments, such as large clusters of commodity machines (Beowulf clusters) or computers specifically designed for high performance computing, such as multi-processor shared-memory machines. We outline the scientific and practical benefits for VLBI correlation, these chiefly being due to the inherent flexibility of software and the fact that the highly parallel and scalable nature of the correlation task is well suited to a multi-processor computing environment. We suggest scientific applications where such an approach to VLBI correlation is most suited and will give the best returns. We report detailed results from the Distributed FX (DiFX) software correlator, running on the Swinburne supercomputer (a Beowulf cluster of approximately 300 commodity processors), including measures of the performance of the system. For example, to correlate all Stokes products for a 10 antenna array, with an aggregate bandwidth of 64 MHz per station and using typical time and frequency resolution presently requires of order 100 desktop-class compute nodes. Due to the effect of Moore's Law on commodity computing performance, the total number and cost of compute nodes required to meet a given correlation task continues to decrease rapidly with time. We show detailed comparisons between DiFX and two existing hardware-based correlators: the Australian Long Baseline Array (LBA) S2 correlator, and the NRAO Very Long Baseline Array (VLBA) correlator. In both cases, excellent agreement was found between the correlators. Finally, we describe plans for the future operation of DiFX on the Swinburne supercomputer, for both astrophysical and geodetic science.

Infrared and submillimetre observations of nearby Vega-like stars have revealed a number of clumpy, asymmetric dust debris disks. Previous studies using semi-analytical and numerical methods have suggested planetary companions of various mass as the likely cause of most examples of disk asymmetry. In this paper, we modify an N-body symplectic gravitational integrator to include radiation terms and conduct medium-resolution parameter searches to identify likely planetary candidates in observed Vega-like systems. We also present high resolution models of Vega and epsilon Eridani, comparing our results to those of previous authors, and a new model for Fomalhaut.