We present the discovery and subarcsecond localization of a new Fast Radio Burst with the Karl G. Jansky Very Large Array and realfast search system. The FRB was discovered on 2019 June 14 with a dispersion measure of 959 pc/cm3. This is the highest DM of any localized FRB and its measured burst fluence of 0.6 Jy ms is less than nearly all other FRBs. The source is not detected to repeat in 15 hours of VLA observing and 153 hours of CHIME/FRB observing. We describe a suite of statistical and data quality tests we used to verify the significance of the event and its localization precision. Follow-up optical/infrared photometry with Keck and Gemini associate the FRB to a pair of galaxies with $\rm{r}\sim23$ mag. The false-alarm rate for radio transients of this significance that are associated with a host galaxy is roughly $3\times10^{-4}\ \rm{hr}^{-1}$. The two putative host galaxies have similar photometric redshifts of $z_{\rm{phot}}\sim0.6$, but different colors and stellar masses. Comparing the host distance to that implied by the dispersion measure suggests a modest (~ 50 pc/cm3) electron column density associated with the FRB environment or host galaxy/galaxies.

We report on the detection of a burst from FRB180916 by realfast/VLA and present software for interpreting fast radio bursts (FRB) periodicity. We demonstrate a range of periodicity analyses with bursts from FRB180916, FRB121102 and FRB180814. Our results for FRB180916 and FRB121102 are consistent with published results. For FRB180814, we did not detect any significant periodic episodes. The realfast-detected and other high-frequency bursts for FRB180916 tend to lie at the beginning of the activity window, indicating a possible phase-frequency relation. The python package $\texttt{frbpa}$ can be used to reproduce and expand on this analysis to test models for repeating FRBs.

With the upcoming commensal surveys for Fast Radio Bursts (FRBs), and their high candidate rate, usage of machine learning algorithms for candidate classification is a necessity. Such algorithms will also play a pivotal role in sending real-time triggers for prompt follow-ups with other instruments. In this paper, we have used the technique of Transfer Learning to train the state-of-the-art deep neural networks for classification of FRB and Radio Frequency Interference (RFI) candidates. These are convolutional neural networks which work on radio frequency-time and dispersion measure-time images as the inputs. We trained these networks using simulated FRBs and real RFI candidates from telescopes at the Green Bank Observatory. We present 11 deep learning models, each with an accuracy and recall above 99.5% on our test dataset comprising of real RFI and pulsar candidates. As we demonstrate, these algorithms are telescope and frequency agnostic and are able to detect all FRBs with signal-to-noise ratios above 10 in ASKAP and Parkes data. We also provide an open-source python package FETCH (Fast Extragalactic Transient Candidate Hunter) for classification of candidates, using our models. Using FETCH, these models can be deployed along with any commensal search pipeline for real-time candidate classification.

We report results of a novel high-energy follow-up observation of a potential Fast Radio Burst. The radio burst was detected by VLA/realfast and followed-up by the Neil Gehrels Swift Observatory in very low latency utilizing new operational capabilities of Swift (arXiv:2005.01751), with pointed soft X-ray and UV observations beginning at T0+32 minutes, and hard X-ray/gamma-ray event data saved around T0. These observations are $>10$x faster than previous X-ray/UV follow-up of any radio transient to date. No emission is seen coincident with the FRB candidate at T0, with a 0.2s fluence $5\sigma$ upper limit of $1.35\times10^{-8}$ erg cm$^{-2}$ (14-195 keV) for a SGR 1935+2154-like flare, nor at T0+32 minutes down to $3\sigma$ upper limits of 22.18 AB mag in UVOT u band, and $3.33\times10^{-13}$ erg cm$^{-2}$ s$^{-1}$ from 0.3-10 keV for the 2 ks observation. The candidate FRB alone is not significant enough to be considered astrophysical, so this note serves as a technical demonstration. These new Swift operational capabilities will allow future FRB detections to be followed up with Swift at even lower latencies than demonstrated here: 15-20 minutes should be regularly achievable, and 5-10 minutes occasionally achievable. We encourage FRB detecting facilities to release alerts in low latency to enable this science.

We present a new analysis of the profile data from the 47 millisecond pulsars comprising the 12.5-year data set of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), which is presented in a parallel paper (Alam et al.; NG12.5). Our reprocessing is performed using "wideband" timing methods, which use frequency-dependent template profiles, simultaneous time-of-arrival (TOA) and dispersion measure (DM) measurements from broadband observations, and novel analysis techniques. In particular, the wideband DM measurements are used to constrain the DM portion of the timing model. We compare the ensemble timing results to NG12.5 by examining the timing residuals, timing models, and noise model components. There is a remarkable level of agreement across all metrics considered. Our best-timed pulsars produce encouragingly similar results to those from NG12.5. In certain cases, such as high-DM pulsars with profile broadening, or sources that are weak and scintillating, wideband timing techniques prove to be beneficial, leading to more precise timing model parameters by 10-15%. The high-precision multi-band measurements in several pulsars indicate frequency-dependent DMs. The TOA volume is reduced by a factor of 33, which may ultimately facilitate computational speed-ups for complex pulsar timing array analyses. This first wideband pulsar timing data set is a stepping stone, and its consistent results with NG12.5 assure us that such data sets are appropriate for gravitational wave analyses.

We present time-of-arrival measurements and timing models of 47 millisecond pulsars (MSPs) observed from 2004 to 2017 at the Arecibo Observatory and the Green Bank Telescope by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). The observing cadence was three to four weeks for most pulsars over most of this time span, with weekly observations of six sources. These data were collected for use in low-frequency gravitational wave searches and for other astrophysical purposes. We detail our observational methods and present a set of time-of-arrival (TOA) measurements, based on "narrowband" analysis, in which many TOAs are calculated within narrow radio-frequency bands for data collected simultaneously across a wide bandwidth. A separate set of "wideband" TOAs will be presented in a companion paper. We detail a number of methodological changes compared to our previous work which yield a cleaner and more uniformly processed data set. Our timing models include several new astrometric and binary pulsar measurements, including previously unpublished values for the parallaxes of PSRs J1832-0836 and J2322+2057, the secular derivatives of the projected semi-major orbital axes of PSRs J0613-0200 and J2229+2643, and the first detection of the Shapiro delay in PSR J2145-0750. We report detectable levels of red noise in the time series for fourteen pulsars. As a check on timing model reliability, we investigate the stability of astrometric parameters across data sets of different lengths. We report flux density measurements for all pulsars observed. Searches for stochastic and continuous gravitational waves using these data will be subjects of forthcoming publications.

When galaxies merge, the supermassive black holes in their centers may form binaries and, during the process of merger, emit low-frequency gravitational radiation in the process. In this paper we consider the galaxy 3C66B, which was used as the target of the first multi-messenger search for gravitational waves. Due to the observed periodicities present in the photometric and astrometric data of the source of the source, it has been theorized to contain a supermassive black hole binary. Its apparent 1.05-year orbital period would place the gravitational wave emission directly in the pulsar timing band. Since the first pulsar timing array study of 3C66B, revised models of the source have been published, and timing array sensitivities and techniques have improved dramatically. With these advances, we further constrain the chirp mass of the potential supermassive black hole binary in 3C66B to less than $(1.65\pm0.02) \times 10^9~{M_\odot}$ using data from the NANOGrav 11-year data set. This upper limit provides a factor of 1.6 improvement over previous limits, and a factor of 4.3 over the first search done. Nevertheless, the most recent orbital model for the source is still consistent with our limit from pulsar timing array data. In addition, we are able to quantify the improvement made by the inclusion of source properties gleaned from electromagnetic data to `blind' pulsar timing array searches. With these methods, it is apparent that it is not necessary to obtain exact a priori knowledge of the period of a binary to gain meaningful astrophysical inferences.

Fast radio bursts (FRBs) are brief, bright, extragalactic radio flashes. Their physical origin remains unknown, but dozens of possible models have been postulated. Some FRB sources exhibit repeat bursts. Though over a hundred FRB sources have been discovered to date, only four have been localised and associated with a host galaxy, with just one of the four known to repeat. The properties of the host galaxies, and the local environments of FRBs, provide important clues about their physical origins. However, the first known repeating FRB has been localised to a low-metallicity, irregular dwarf galaxy, and the apparently non-repeating sources to higher-metallicity, massive elliptical or star-forming galaxies, suggesting that perhaps the repeating and apparently non-repeating sources could have distinct physical origins. Here we report the precise localisation of a second repeating FRB source, FRB 180916.J0158+65, to a star-forming region in a nearby (redshift $z = 0.0337 \pm 0.0002$) massive spiral galaxy, whose properties and proximity distinguish it from all known hosts. The lack of both a comparably luminous persistent radio counterpart and a high Faraday rotation measure further distinguish the local environment of FRB 180916.J0158+65 from that of the one previously localised repeating FRB source, FRB 121102. This demonstrates that repeating FRBs have a wide range of luminosities, and originate from diverse host galaxies and local environments.

The regularity of pulsar emissions becomes apparent once we reference the pulses' times of arrivals to the inertial rest frame of the solar system. It follows that errors in the determination of Earth's position with respect to the solar-system barycenter can appear as a time-correlated bias in pulsar-timing residual time series, affecting the searches for low-frequency gravitational waves performed with pulsar timing arrays. Indeed, recent array datasets yield different gravitational-wave background upper limits and detection statistics when analyzed with different solar-system ephemerides. Crucially, the ephemerides do not generally provide usable error representations. In this article we describe the motivation, construction, and application of a physical model of solar-system ephemeris uncertainties, which focuses on the degrees of freedom (Jupiter's orbital elements) most relevant to gravitational-wave searches with pulsar timing arrays. This model, BayesEphem, was used to derive ephemeris-robust results in NANOGrav's 11-yr stochastic-background search, and it provides a foundation for future searches by NANOGrav and other consortia. The analysis and simulations reported here suggest that ephemeris modeling reduces the gravitational-wave sensitivity of the 11-yr dataset; and that this degeneracy will vanish with improved ephemerides and with the longer pulsar timing datasets that will become available in the near future.

We present the results of processing an additional 44% of the High Time Resolution Universe South Low Latitude (HTRU-S LowLat) pulsar survey, the most sensitive blind pulsar survey of the southern Galactic plane to date. Our partially-coherent segmented acceleration search pipeline is designed to enable the discovery of pulsars in short, highly-accelerated orbits, while our 72-min integration lengths will allow us to discover pulsars at the lower end of the pulsar luminosity distribution. We report the discovery of 40 pulsars, including three millisecond pulsar-white dwarf binary systems (PSRs J1537-5312, J1547-5709 and J1618-4624), a black-widow binary system (PSR J1745-23) and a candidate black-widow binary system (PSR J1727-2951), a glitching pulsar (PSR J1706-4434), an eclipsing binary pulsar with a 1.5-yr orbital period (PSR J1653-45), and a pair of long spin-period binary pulsars which display either nulling or intermittent behaviour (PSRs J1812-15 and J1831-04). We show that the total population of 100 pulsars discovered in the HTRU-S LowLat survey to date represents both an older and lower-luminosity population, and indicates that we have yet to reach the bottom of the luminosity distribution function. We present evaluations of the performance of our search technique and of the overall yield of the survey, considering the 94% of the survey which we have processed to date. We show that our pulsar yield falls below earlier predictions by approximately 25% (especially in the case of millisecond pulsars), and discuss explanations for this discrepancy as well as future adaptations in RFI mitigation and searching techniques which may address these shortfalls.

The Very Large Array Sky Survey (VLASS) is a synoptic, all-sky radio sky survey with a unique combination of high angular resolution ($\approx$2.5"), sensitivity (a 1$\sigma$ goal of 70 $\mu$Jy/beam in the coadded data), full linear Stokes polarimetry, time domain coverage, and wide bandwidth (2-4 GHz). The first observations began in September 2017, and observing for the survey will finish in 2024. VLASS will use approximately 5500 hours of time on the Karl G. Jansky Very Large Array (VLA) to cover the whole sky visible to the VLA (Declination $>-40^{\circ}$), a total of 33,885 deg$^2$. The data will be taken in three epochs to allow the discovery of variable and transient radio sources. The survey is designed to engage radio astronomy experts, multi-wavelength astronomers, and citizen scientists alike. By utilizing an "on the fly" interferometry mode, the observing overheads are much reduced compared to a conventional pointed survey. In this paper, we present the science case and observational strategy for the survey, and also results from early survey observations.

The mergers of supermassive black hole binaries (SMBHBs) promise to be incredible sources of gravitational waves (GWs). While the oscillatory part of the merger gravitational waveform will be outside the frequency sensitivity range of pulsar timing arrays (PTAs), the non-oscillatory GW memory effect is detectable. Further, any burst of gravitational waves will produce GW memory, making memory a useful probe of unmodeled exotic sources and new physics. We searched the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 11-year data set for GW memory. This dataset is sensitive to very low frequency GWs of $\sim3$ to $400$ nHz (periods of $\sim11$ yr $-$ $1$ mon). Finding no evidence for GWs, we placed limits on the strain amplitude of GW memory events during the observation period. We then used the strain upper limits to place limits on the rate of GW memory causing events. At a strain of $2.5\times10^{-14}$, corresponding to the median upper limit as a function of source sky position, we set a limit on the rate of GW memory events at $<0.4$ yr$^{-1}$. That strain corresponds to a SMBHB merger with reduced mass of $\eta M \sim 2\times10^{10}M_\odot$ and inclination of $\iota=\pi/3$ at a distance of 1 Gpc. As a test of our analysis, we analyzed the NANOGrav 9-year data set as well. This analysis found an anomolous signal, which does not appear in the 11-year data set. This signal is not a GW, and its origin remains unknown.

We have constructed a new timescale, TT(IPTA16), based on observations of radio pulsars presented in the first data release from the International Pulsar Timing Array (IPTA). We used two analysis techniques with independent estimates of the noise models for the pulsar observations and different algorithms for obtaining the pulsar timescale. The two analyses agree within the estimated uncertainties and both agree with TT(BIPM17), a post-corrected timescale produced by the Bureau International des Poids et Mesures (BIPM). We show that both methods could detect significant errors in TT(BIPM17) if they were present. We estimate the stability of the atomic clocks from which TT(BIPM17) is derived using observations of four rubidium fountain clocks at the US Naval Observatory. Comparing the power spectrum of TT(IPTA16) with that of these fountain clocks suggests that pulsar-based timescales are unlikely to contribute to the stability of the best timescales over the next decade, but they will remain a valuable independent check on atomic timescales. We also find that the stability of the pulsar-based timescale is likely to be limited by our knowledge of solar-system dynamics, and that errors in TT(BIPM17) will not be a limiting factor for the primary goal of the IPTA, which is to search for the signatures of nano-Hertz gravitational waves.

We present results of a search for late-time radio emission and Fast Radio Bursts (FRBs) from a sample of type-I superluminous supernovae (SLSNe-I). We used the Karl G. Jansky Very Large Array to observe ten SLSN-I more than 5 years old at a frequency of 3 GHz. We searched fast-sampled visibilities for FRBs and used the same data to perform a deep imaging search for late-time radio emission expected in models of magnetar-powered supernovae. No FRBs were found. One SLSN-I, PTF10hgi, is detected in deep imaging, corresponding to a luminosity of $1.2\times10^{28}$ erg s$^{-1}$. This luminosity, considered with the recent 6 GHz detection of PTF10hgi in Eftekhari et al (2019), supports the interpretation that it is powered by a young, fast-spinning ($\sim$ ms spin period) magnetar with $\sim$ 15 Msun of partially ionized ejecta. Broadly, our observations are most consistent with SLSNe-I being powered by neutron stars with fast spin periods, although most require more free-free absorption than is inferred for PTF10hgi. We predict that radio observations at higher frequencies or in the near future will detect these systems and begin constraining properties of the young pulsars and their birth environments.

An ensemble of inspiraling supermassive black hole binaries should produce a stochastic background of very low frequency gravitational waves. This stochastic background is predicted to be a power law, with a spectral index of -2/3, and it should be detectable by a network of precisely timed millisecond pulsars, widely distributed on the sky. This paper reports a new "time slicing" analysis of the 11-year data release from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) using 34 millisecond pulsars. Methods to flag potential "false positive" signatures are developed, including techniques to identify responsible pulsars. Mitigation strategies are then presented. We demonstrate how an incorrect noise model can lead to spurious signals, and show how independently modeling noise across 30 Fourier components, spanning NANOGrav's frequency range, effectively diagnoses and absorbs the excess power in gravitational-wave searches. This results in a nominal, and expected, progression of our gravitational-wave statistics. Additionally we show that the first interstellar medium event in PSR J1713+0747 pollutes the common red noise process with low-spectral index noise, and use a tailored noise model to remove these effects.

In this paper, we describe the International Pulsar Timing Array second data release, which includes recent pulsar timing data obtained by three regional consortia: the European Pulsar Timing Array, the North American Nanohertz Observatory for Gravitational Waves, and the Parkes Pulsar Timing Array. We analyse and where possible combine high-precision timing data for 65 millisecond pulsars which are regularly observed by these groups. A basic noise analysis, including the processes which are both correlated and uncorrelated in time, provides noise models and timing ephemerides for the pulsars. We find that the timing precisions of pulsars are generally improved compared to the previous data release, mainly due to the addition of new data in the combination. The main purpose of this work is to create the most up-to-date IPTA data release. These data are publicly available for searches for low-frequency gravitational waves and other pulsar science.

The analogy of the host galaxy of the repeating fast radio burst (FRB) source FRB 121102 and those of long gamma-ray bursts (GRBs) and super-luminous supernovae (SLSNe) has led to the suggestion that young magnetars born in GRBs and SLSNe could be the central engine of repeating FRBs. We test such a hypothesis by performing dedicated observations of the remnants of six GRBs with evidence of having a magnetar central engine using the Arecibo telescope and the Robert C. Byrd Green Bank Telescope (GBT). A total of $\sim 20$ hrs of observations of these sources did not detect any FRB from these remnants. Under the assumptions that all these GRBs left behind a long-lived magnetar and that the bursting rate of FRB 121102 is typical for a magnetar FRB engine, we estimate a non-detection probability of $8.9\times10^{-6}$. Even though these non-detections cannot exclude the young magnetar model of FRBs, we place constraints on the burst rate and luminosity function of FRBs from these GRB targets.

Over the last decade, the continuing decline in the cost of digital computing technology has brought about a dramatic transformation in how digital instrumentation for radio astronomy is developed and operated. In most cases, it is now possible to interface consumer computing hardware, e.g. inexpensive graphics processing units and storage devices, directly to the raw data streams produced by radio telescopes. Such systems bring with them myriad benefits: straightforward upgrade paths, cost savings through leveraging an economy of scale, and a lowered barrier to entry for scientists and engineers seeking to add new instrument capabilities. Additionally, the typical data-interconnect technology used with general-purpose computing hardware -- Ethernet -- naturally permits multiple subscribers to a single raw data stream. This allows multiple science programs to be conducted in parallel. When combined with broad bandwidths and wide primary fields of view, radio telescopes become capable of achieving many science goals simultaneously. Moreover, because many science programs are not strongly dependent on observing cadence and direction (e.g. searches for extraterrestrial intelligence and radio transient surveys), these so-called "commensal" observing programs can dramatically increase the scientific productivity and discovery potential of an observatory. In this whitepaper, we detail a project to add an Ethernet-based commensal observing mode to the Jansky Very Large Array (VLA), and discuss how this mode could be leveraged to conduct a powerful program to constrain the distribution of advanced life in the universe through a search for radio emission indicative of technology. We also discuss other potential science use-cases for the system, and how the system could be used for technology development towards next-generation processing systems for the Next Generation VLA.

Pulsar timing arrays (PTAs) will enable the detection of nanohertz gravitational waves (GWs) from a population of supermassive binary black holes (SMBBHs) in the next $\sim 3-7$ years. In addition, PTAs provide a rare opportunity to probe exotic physics. Potential sources of GWs in the nanohertz band include, cosmic strings and cosmic superstrings, inflation and, phase transitions in the early universe. GW observations will also make possible tests of gravitational theories that, by modifying Einstein's theory of general relativity, attempt to explain the origin of cosmic acceleration and reconcile quantum mechanics and gravity, two of the most profound challenges facing fundamental physics today. Finally, PTAs also provide a new means to probe certain dark matter models.

Gravitational wave memory is theorized to arise from the integrated history of gravitational wave emission, and manifests as a spacetime deformation in the wake of a propagating gravitational wave. We explore the detectability of the memory signals from a population of coalescencing supermassive black hole binaries with pulsar timing arrays and the Laser Interferometer Space Antenna (LISA). We find that current pulsar timing arrays have poor prospects, but it is likely that between 1 and 10 memory events with signal-to-noise ratio in excess of 5 will occur within LISA's planned 4-year mission.

The Parkes Pulsar Timing Array (PPTA) project monitors two dozen millisecond pulsars (MSPs) in order to undertake a variety of fundamental physics experiments using the Parkes 64m radio telescope. Since June 2017 we have been undertaking commensal searches for fast radio bursts (FRBs) during the MSP observations. Here, we report the discovery of four FRBs (171209, 180309, 180311 and 180714). The detected events include an FRB with the highest signal-to-noise ratio ever detected at the Parkes observatory, which exhibits unusual spectral properties. All four FRBs are highly polarized. We discuss the future of commensal searches for FRBs at Parkes.

Observations indicate that nearly all galaxies contain supermassive black holes (SMBHs) at their centers. When galaxies merge, their component black holes form SMBH binaries (SMBHBs), which emit low-frequency gravitational waves (GWs) that can be detected by pulsar timing arrays (PTAs). We have searched the recently-released North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 11-year data set for GWs from individual SMBHBs in circular orbits. As we did not find strong evidence for GWs in our data, we placed 95\% upper limits on the strength of GWs from such sources as a function of GW frequency and sky location. We placed a sky-averaged upper limit on the GW strain of $h_0 < 7.3(3) \times 10^{-15}$ at $f_\mathrm{gw}= 8$ nHz. We also developed a technique to determine the significance of a particular signal in each pulsar using ``dropout' parameters as a way of identifying spurious signals in measurements from individual pulsars. We used our upper limits on the GW strain to place lower limits on the distances to individual SMBHBs. At the most-sensitive sky location, we ruled out SMBHBs emitting GWs with $f_\mathrm{gw}= 8$ nHz within 120 Mpc for $\mathcal{M} = 10^9 \, M_\odot$, and within 5.5 Gpc for $\mathcal{M} = 10^{10} \, M_\odot$. We also determined that there are no SMBHBs with $\mathcal{M} > 1.6 \times 10^9 \, M_\odot$ emitting GWs in the Virgo Cluster. Finally, we estimated the number of potentially detectable sources given our current strain upper limits based on galaxies in Two Micron All-Sky Survey (2MASS) and merger rates from the Illustris cosmological simulation project. Only 34 out of 75,000 realizations of the local Universe contained a detectable source, from which we concluded it was unsurprising that we did not detect any individual sources given our current sensitivity to GWs.

Pulsar timing array (PTA) collaborations in North America, Australia, and Europe, have been exploiting the exquisite timing precision of millisecond pulsars over decades of observations to search for correlated timing deviations induced by gravitational waves (GWs). PTAs are sensitive to the frequency band ranging just below 1 nanohertz to a few tens of microhertz. The discovery space of this band is potentially rich with populations of inspiraling supermassive black-holes binaries, decaying cosmic string networks, relic post-inflation GWs, and even non-GW imprints of axionic dark matter. This article aims to provide an understanding of the exciting open science questions in cosmology, galaxy evolution, and fundamental physics that will be addressed by the detection and study of GWs through PTAs. The focus of the article is on providing an understanding of the mechanisms by which PTAs can address specific questions in these fields, and to outline some of the subtleties and difficulties in each case. The material included is weighted most heavily towards the questions which we expect will be answered in the near-term with PTAs; however, we have made efforts to include most currently anticipated applications of nanohertz GWs.

Gravitational wave (GW) and gravitational slingshot recoil kicks, which are natural products of SMBH evolution in merging galaxies, can produce active galactic "nuclei" that are offset from the centers of their host galaxies. Detections of offset AGN would provide key constraints on SMBH binary mass and spin evolution and on GW event rates. Although numerous offset AGN candidates have been identified, none have been definitively confirmed. Multi-wavelength observations with next-generation telescopes, including systematic large-area surveys, will provide unprecedented opportunities to identify and confirm candidate offset AGN from sub-parsec to kiloparsec scales. We highlight ways in which these observations will open a new avenue for multi-messenger studies in the dawn of low-frequency (~ nHz - mHz) GW astronomy.

Precision timing of large arrays (>50) of millisecond pulsars will detect the nanohertz gravitational-wave emission from supermassive binary black holes within the next ~3-7 years. We review the scientific opportunities of these detections, the requirements for success, and the synergies with electromagnetic instruments operating in the 2020s.

Pulsar timing arrays (PTAs) are on the verge of detecting low-frequency gravitational waves (GWs) from supermassive black hole binaries (SMBHBs). With continued observations of a large sample of millisecond pulsars, PTAs will reach this major milestone within the next decade. Already, SMBHB candidates are being identified by electromagnetic surveys in ever-increasing numbers; upcoming surveys will enhance our ability to detect and verify candidates, and will be instrumental in identifying the host galaxies of GW sources. Multi-messenger (GW and electromagnetic) observations of SMBHBs will revolutionize our understanding of the co-evolution of SMBHs with their host galaxies, the dynamical interactions between binaries and their galactic environments, and the fundamental physics of accretion. Multi-messenger observations can also make SMBHBs 'standard sirens' for cosmological distance measurements out to $z\simeq0.5$. LIGO has already ushered in breakthrough insights in our knowledge of black holes. The multi-messenger detection of SMBHBs with PTAs will be a breakthrough in the years $2020-2030$ and beyond, and prepare us for LISA to help complete our views of black hole demographics and evolution at higher redshifts.

Coalescing, massive black-hole (MBH) binaries are the most powerful sources of gravitational waves (GWs) in the Universe, which makes MBH science a prime focus for ongoing and upcoming GW observatories. The Laser Interferometer Space Antenna (LISA) -- a gigameter scale space-based GW observatory -- will grant us access to an immense cosmological volume, revealing MBHs merging when the first cosmic structures assembled in the Dark Ages. LISA will unveil the yet unknown origin of the first quasars, and detect the teeming population of MBHs of $10^4 - 10^7$ solar masses. forming within protogalactic halos. The Pulsar Timing Array, a galactic-scale GW survey, can access the largest MBHs the Universe, detecting the cosmic GW foreground from inspiraling MBH binaries of about 10^9 solar masses. LISA can measure MBH spins and masses with precision far exceeding that from electromagnetic (EM) probes, and together, both GW observatories will provide the first full census of binary MBHs, and their orbital dynamics, across cosmic time. Detecting the loud gravitational signal of these MBH binaries will also trigger alerts for EM counterpart searches, from decades (PTAs) to hours (LISA) prior to the final merger. By witnessing both the GW and EM signals of MBH mergers, precious information will be gathered about the rich and complex environment in the aftermath of a galaxy collision. The unique GW characterization of MBHs will shed light on the deep link between MBHs of $10^4-10^{10}$ solar masses and the grand design of galaxy assembly, as well as on the complex dynamics that drive MBHs to coalescence.

The discovery of Fast Radio Bursts (FRBs) at cosmological distances has opened a powerful window on otherwise unseen matter in the Universe. In the 2020s, observations of $>10^{4}$ FRBs will assess the baryon contents and physical conditions in the hot/diffuse circumgalactic, intracluster, and intergalactic medium, and test extant compact-object dark matter models.

The field of time-domain astrophysics has entered the era of Multi-messenger Astronomy (MMA). One key science goal for the next decade (and beyond) will be to characterize gravitational wave (GW) and neutrino sources using the next generation of Extremely Large Telescopes (ELTs). These studies will have a broad impact across astrophysics, informing our knowledge of the production and enrichment history of the heaviest chemical elements, constrain the dense matter equation of state, provide independent constraints on cosmology, increase our understanding of particle acceleration in shocks and jets, and study the lives of black holes in the universe. Future GW detectors will greatly improve their sensitivity during the coming decade, as will near-infrared telescopes capable of independently finding kilonovae from neutron star mergers. However, the electromagnetic counterparts to high-frequency (LIGO/Virgo band) GW sources will be distant and faint and thus demand ELT capabilities for characterization. ELTs will be important and necessary contributors to an advanced and complete multi-messenger network.

The nearby galaxy NGC 3115 contains a known radio-emitting, low-luminosity active galactic nucleus (AGN), and was recently claimed to host a candidate AGN displaced 14.3 pc from the galaxy's optical photocenter. Our goal is to understand whether this represents a single offset AGN, an AGN in orbit around a central black hole, or something else. We present a new, sensitive (RMS = 4.4 $\mu$Jy beam$^{-1}$) 10 GHz image, which finds evidence for only one AGN. We place a stringent limit on the radio luminosity of any secondary supermassive black hole of $L_{10~\rm{GHz}}<5.8\times10^{33}$ ergs/s. An analysis of the relative positioning of the radio core, X-ray nucleus, and stellar bulge in this galaxy indicate that the radio source is centrally located, and not offset from the galactic bulge. This provides an argument against a single offset AGN in NGC 3115, however does not provide conclusive evidence against the purported offset AGN as an in-spiralling secondary black hole.

We report on the latest six pulsars discovered through our standard pipeline in the intermediate-latitude region (|b| < 15 deg) of the Parkes High Time Resolution Universe Survey (HTRU). We also present timing solutions for the new discoveries and for 19 further pulsars for which only discovery parameters were previously published. Highlights of the presented sample include the isolated millisecond pulsar J1826-2415, the long-period binary pulsar J1837-0822 in a mildly eccentric 98-day orbit with a > 0.27 M_sun companion, and the nulling pulsar J1638-4233, detected only 10% of the time. Other interesting objects are PSR J1757-1500, exhibiting sporadic mode changes, and PSR J1635-2616 showing one glitch over 6 years. The new discoveries bring the total count of HTRU intermediate-latitude pulsars to 113, 25% of which are recycled pulsars. This is the higest ratio of recycled over ordinary pulsars discoveries of all recent pulsar surveys in this region of the sky. Among HTRU recycled pulsars, four are isolated objects. Comparing the characteristics of Galactic fully-recycled isolated MSPs with those of eclipsing binaries ('spiders'), from which the former are believed to have formed, we highlight a discrepancy in their spatial distribution. This may reflect a difference in the natal kick, hence, possibly, a different formation path. On the other hand, however, isolated fully-recycled MSPs spin periods are, on average, longer than those of spiders, in line with what one would expect, from simple magnetic-dipole spin-down, if the former were indeed evolved from the latter.

We report deep, single-dish 21 cm observations of NGC 1052-DF2, taken with the Green Bank Telescope. NGC 1052-DF2, proposed to be lacking in dark matter, is currently classified as an ultra-diffuse galaxy in the NGC 1052 group. We do not detect the galaxy, and derive an upper limit on the HI, mass. The galaxy is extremely gas-poor, and we find that a $3\sigma \, M_{HI}$ detection at a distance of 19 Mpc and using a line width of 3.2 km $\rm s^{-1}$ would have an upper limit of $M_{HI,lim} < 5.5 \times 10^5$ M$_{\odot}$. At this mass limit, the gas fraction of neutral gas mass to stellar mass is extremely low, at $M_{HI}$/M$_{*}$ $\, < \, 0.0027$. This extremely low gas fraction, comparable to Galactic dwarf spheroidals and gas-poor dwarf ellipticals, implies that either the galaxy is within the virial radius of NGC1052, where its gas has been stripped due to its proximity to the central galaxy, or that NGC 1052-DF2 is at distance large enough to inhibit detection of its gas. We also estimated the upper limit of the HI mass of NGC 1052-DF2 resided at 13 Mpc. This would give an HI mass of $M_{HI,lim} < 2.5 \times 10^5$ M$_{\odot}$, and an HI gas fraction of $M_{HI}$/M$_{*}$r $\, < \, 0.0012$, becoming even more extreme. While the dark matter fraction would be less extreme at this distance, the neutral gas fraction would be unprecedented for an object in a low density environment.

We have performed a new search for radio pulsars in archival data of the intermediate and high Galactic latitude parts of the Southern High Time Resolution Universe pulsar survey. This is the first time the entire dataset has been searched for binary pulsars, an achievement enabled by GPU-accelerated dedispersion and periodicity search codes nearly 50 times faster than the previously used pipeline. Candidate selection was handled entirely by a Machine Learning algorithm, allowing for the assessment of 17.6 million candidates in a few person-days. We have also introduced an outlier detection algorithm for efficient radio-frequency interference (RFI) mitigation on folded data, a new approach that enabled the discovery of pulsars previously masked by RFI. We discuss implications for future searches, particularly the importance of expanding work on RFI mitigation to improve survey completeness. In total we discovered 23 previously unknown sources, including 6 millisecond pulsars and at least 4 pulsars in binary systems. We also found an elusive but credible redback candidate that we have yet to confirm.

We review the potential information that can be gleaned from multi-wavelength and multi-messenger detections of Fast Radio Bursts (FRBs). The primary function of such detections, following FRB localization, will be to reveal the progenitors of FRB events. However, these observations serve a more basic function in host galaxy identification, subsequent redshift measurement, and in disentangling the contributions of various line-of-sight contributions to properties of FRBs measured at radio wavelengths.

Fast radio bursts (FRBs) are millisecond pulses of radio emission of seemingly extragalactic origin. More than 50 FRBs have now been detected, with only one seen to repeat. Here we present a new FRB discovery, FRB 110214, which was detected in the high latitude portion of the High Time Resolution Universe South survey at the Parkes telescope. FRB 110214 has one of the lowest dispersion measures of any known FRB (DM = 168.9$\pm$0.5 pc cm$^{-3}$), and was detected in two beams of the Parkes multi-beam receiver. A triangulation of the burst origin on the sky identified three possible regions in the beam pattern where it may have originated, all in sidelobes of the primary detection beam. Depending on the true location of the burst the intrinsic fluence is estimated to fall in the range of 50 -- 2000 Jy ms, making FRB 110214 one of the highest-fluence FRBs detected with the Parkes telescope. No repeating pulses were seen in almost 100 hours of follow-up observations with the Parkes telescope down to a limiting fluence of 0.3 Jy ms for a 2-ms pulse. Similar low-DM, ultra-bright FRBs may be detected in telescope sidelobes in the future, making careful modeling of multi-beam instrument beam patterns of utmost importance for upcoming FRB surveys.

Gravitational-wave (GW) and gravitational slingshot recoil kicks, which are natural products of SMBH evolution in merging galaxies, can produce active galactic "nuclei" that are offset from the centers of their host galaxies. Detections of offset AGN would provide key constraints on SMBH binary mass and spin evolution and on GW event rates. Although numerous offset AGN candidates have been identified, none have been definitively confirmed. The ngVLA offers unparalleled capabilities to identify and confirm candidate offset AGN from sub-parsec to kiloparsec scales, opening a new avenue for multi-messenger studies in the dawn of low-frequency GW astronomy.

Pulsar-timing analyses are sensitive to errors in the solar-system ephemerides (SSEs) that timing models utilise to estimate the location of the solar-system barycentre, the quasi-inertial reference frame to which all recorded pulse times-of-arrival are referred. Any error in the SSE will affect all pulsars, therefore pulsar timing arrays (PTAs) are a suitable tool to search for such errors and impose independent constraints on relevant physical parameters. We employ the first data release of the International Pulsar Timing Array to constrain the masses of the planet-moons systems and to search for possible unmodelled objects (UMOs) in the solar system. We employ ten SSEs from two independent research groups, derive and compare mass constraints of planetary systems, and derive the first PTA mass constraints on asteroid-belt objects. Constraints on planetary-system masses have been improved by factors of up to 20 from the previous relevant study using the same assumptions, with the mass of the Jovian system measured at 9.5479189(3)$\times10^{-4}$ $M_{\odot}$. The mass of the dwarf planet Ceres is measured at 4.7(4)$\times10^{-10}$ $M_{\odot}$. We also present the first sensitivity curves using real data that place generic limits on the masses of UMOs, which can also be used as upper limits on the mass of putative exotic objects. For example, upper limits on dark-matter clumps are comparable to published limits using independent methods. While the constraints on planetary masses derived with all employed SSEs are consistent, we note and discuss differences in the associated timing residuals and UMO sensitivity curves.

The Next-Generation Very Large Array (ngVLA) has the potential to be a workhorse for the discovery and study of paired supermassive black holes either at large separations (dual) or in tightly bound systems (binary). In this chapter, we outline the science case for the study of these supermassive pairs, and summarize discovery methods that can be used at radio wavelengths to discover them: including morphological, spectral, and time-domain identifications. One critical aspect of this work is that multi-messenger binary black hole studies may be possible with the ngVLA when combined with gravitational-wave searches using pulsar timing array techniques. However, long-baseline interferometery (>>1000 km) will make this possibility more likely by expanding the redshift range at which radio emission arising from two separate black holes may be resolved and studied.

The past year has witnessed discovery of the first identified counterparts to a gravitational wave transient (GW 170817A) and a very high-energy neutrino (IceCube-170922A). These source identifications, and ensuing detailed studies, have realized longstanding dreams of astronomers and physicists to routinely carry out observations of cosmic sources by other than electromagnetic means, and inaugurated the era of "multi-messenger" astronomy. While this new era promises extraordinary physical insights into the universe, it brings with it new challenges, including: highly heterogeneous, high-volume, high-velocity datasets; globe-spanning cross-disciplinary teams of researchers, regularly brought together into transient collaborations; an extraordinary breadth and depth of domain-specific knowledge and computing resources required to anticipate, model, and interpret observations; and the routine need for adaptive, distributed, rapid-response observing campaigns to fully exploit the scientific potential of each source. We argue, therefore, that the time is ripe for the community to conceive and propose an Institute for Multi-Messenger Astrophysics that would coordinate its resources in a sustained and strategic fashion to efficiently address these challenges, while simultaneously serving as a center for education and key supporting activities. In this fashion, we can prepare now to realize the bright future that we see, beyond, through these newly opened windows onto the universe.

We search for an isotropic stochastic gravitational-wave background (GWB) in the newly released $11$-year dataset from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). While we find no significant evidence for a GWB, we place constraints on a GWB from a population of supermassive black-hole binaries, cosmic strings, and a primordial GWB. For the first time, we find that the GWB upper limits and detection statistics are sensitive to the Solar System ephemeris (SSE) model used, and that SSE errors can mimic a GWB signal. We developed an approach that bridges systematic SSE differences, producing the first PTA constraints that are robust against SSE uncertainties. We thus place a $95\%$ upper limit on the GW strain amplitude of $A_\mathrm{GWB}<1.45\times 10^{-15}$ at a frequency of $f=1$ yr$^{-1}$ for a fiducial $f^{-2/3}$ power-law spectrum, and with inter-pulsar correlations modeled. This is a factor of $\sim 2$ improvement over the NANOGrav $9$-year limit, calculated using the same procedure. Previous PTA upper limits on the GWB will need revision in light of SSE systematic uncertainties. We use our constraints to characterize the combined influence on the GWB of the stellar mass-density in galactic cores, the eccentricity of SMBH binaries, and SMBH--galactic-bulge scaling relationships. We constrain cosmic-string tension using recent simulations, yielding an SSE-marginalized $95\%$ upper limit on the cosmic string tension of $G\mu < 5.3\times 10^{-11}$---a factor of $\sim 2$ better than the published NANOGrav $9$-year constraints. Our SSE-marginalized $95\%$ upper limit on the energy density of a primordial GWB (for a radiation-dominated post-inflation Universe) is $\Omega_\mathrm{GWB}(f)h^2<3.4\times10^{-10}$.

We report the first detections of the repeating fast radio burst source FRB 121102 above 5.2 GHz. Observations were performed using the 4$-$8 GHz receiver of the Robert C. Byrd Green Bank Telescope with the Breakthrough Listen digital backend. We present the spectral, temporal and polarization properties of 21 bursts detected within the first 60 minutes of a total 6-hour observations. These observations comprise the highest burst density yet reported in the literature, with 18 bursts being detected in the first 30 minutes. A few bursts clearly show temporal sub-structures with distinct spectral properties. These sub-structures superimpose to provide enhanced peak signal-to-noise ratio at higher trial dispersion measures. Broad features occur in $\sim 1$ GHz wide subbands that typically differ in peak frequency between bursts within the band. Finer-scale structures ($\sim 10-50$ MHz) within these bursts are consistent with that expected from Galactic diffractive interstellar scintillation. The bursts exhibit nearly 100% linear polarization, and a large average rotation measure of 9.359$\pm$0.012 $\times$ 10$^{\rm 4}$ rad m$^{\rm -2}$ (in the observer's frame). No circular polarization was found for any burst. We measure an approximately constant polarization position angle in the 13 brightest bursts. The peak flux densities of the reported bursts have average values (0.2$\pm$0.1 Jy), similar to those seen at lower frequencies ($<3$ GHz), while the average burst widths (0.64$\pm$0.46 ms) are relatively narrower.

We present high-precision timing data over time spans of up to 11 years for 45 millisecond pulsars observed as part of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) project, aimed at detecting and characterizing low-frequency gravitational waves. The pulsars were observed with the Arecibo Observatory and/or the Green Bank Telescope at frequencies ranging from 327 MHz to 2.3 GHz. Most pulsars were observed with approximately monthly cadence, with six high--timing-precision pulsars observed weekly, and all were observed at widely separated frequencies at each observing epoch in order to fit for time-variable dispersion delays. We describe our methods for data processing, time-of-arrival (TOA) calculation, and the implementation of a new, automated method for removing outlier TOAs. We fit a timing model for each pulsar that includes spin, astrometric, and, if necessary, binary parameters, in addition to time-variable dispersion delays and parameters that quantify pulse-profile evolution with frequency. The new timing solutions provide three new parallax measurements, two new Shapiro delay measurements, and two new measurements of large orbital-period variations. We fit models that characterize sources of noise for each pulsar. We find that 11 pulsars show significant red noise, with generally smaller spectral indices than typically measured for non-recycled pulsars, possibly suggesting a different origin. Future papers will use these data to constrain or detect the signatures of gravitational-wave signals.

Pulsars in relativistic binary systems have emerged as fantastic natural laboratories for testing theories of gravity, the most prominent example being the double pulsar, PSR J0737$-$3039. The HTRU-South Low Latitude pulsar survey represents one of the most sensitive blind pulsar surveys taken of the southern Galactic plane to date, and its primary aim has been the discovery of new relativistic binary pulsars. Here we present our binary pulsar searching strategy and report on the survey's flagship discovery, PSR J1757$-$1854. A 21.5-ms pulsar in a relativistic binary with an orbital period of 4.4 hours and an eccentricity of 0.61, this double neutron star (DNS) system is the most accelerated pulsar binary known, and probes a relativistic parameter space not yet explored by previous pulsar binaries.

Radio interferometers have the ability to precisely localize and better characterize the properties of sources. This ability is having a powerful impact on the study of fast radio transients, where a few milliseconds of data is enough to pinpoint a source at cosmological distances. However, recording interferometric data at millisecond cadence produces a terabyte-per-hour data stream that strains networks, computing systems, and archives. This challenge mirrors that of other domains of science, where the science scope is limited by the computational architecture as much as the physical processes at play. Here, we present a solution to this problem in the context of radio transients: realfast, a commensal, fast transient search system at the Jansky Very Large Array. Realfast uses a novel architecture to distribute fast-sampled interferometric data to a 32-node, 64-GPU cluster for real-time imaging and transient detection. By detecting transients in situ, we can trigger the recording of data for those rare, brief instants when the event occurs and reduce the recorded data volume by a factor of 1000. This makes it possible to commensally search a data stream that would otherwise be impossible to record. This system will search for millisecond transients in more than 1000 hours of data per year, potentially localizing several Fast Radio Bursts, pulsars, and other sources of impulsive radio emission. We describe the science scope for realfast, the system design, expected outcomes, and ways real-time analysis can help in other fields of astrophysics.

We report the discovery of PSR J1757$-$1854, a 21.5-ms pulsar in a highly-eccentric, 4.4-h orbit around a neutron star (NS) companion. PSR J1757$-$1854 exhibits some of the most extreme relativistic parameters of any known pulsar, including the strongest relativistic effects due to gravitational-wave (GW) damping, with a merger time of 76 Myr. Following a 1.6-yr timing campaign, we have measured five post-Keplerian (PK) parameters, yielding the two component masses ($m_\text{p}=1.3384(9)\,\text{M}_\odot$ and $m_\text{c}=1.3946(9)\,\text{M}_\odot$) plus three tests of general relativity (GR), which the theory passes. The larger mass of the NS companion provides important clues regarding the binary formation of PSR J1757$-$1854. With simulations suggesting 3-$\sigma$ measurements of both the contribution of Lense-Thirring precession to the rate of change of the semi-major axis and the relativistic deformation of the orbit within $\sim7-9$ years, PSR J1757$-$1854 stands out as a unique laboratory for new tests of gravitational theories.

We report the discovery of four Fast Radio Bursts (FRBs) in the ongoing SUrvey for Pulsars and Extragalactic Radio Bursts (SUPERB) at the Parkes Radio Telescope: FRBs 150610, 151206, 151230 and 160102. Our real-time discoveries have enabled us to conduct extensive, rapid multi-messenger follow-up at 12 major facilities sensitive to radio, optical, X-ray, gamma-ray photons and neutrinos on time scales ranging from an hour to a few months post-burst. No counterparts to the FRBs were found and we provide upper limits on afterglow luminosities. None of the FRBs were seen to repeat. Formal fits to all FRBs show hints of scattering while their intrinsic widths are unresolved in time. FRB 151206 is at low Galactic latitude, FRB 151230 shows a sharp spectral cutoff, and FRB 160102 has the highest dispersion measure (DM = $2596.1\pm0.3$ pc cm$^{-3}$) detected to date. Three of the FRBs have high dispersion measures (DM >$1500$ pc cm$^{-3}$), favouring a scenario where the DM is dominated by contributions from the Intergalactic Medium. The slope of the Parkes FRB source counts distribution with fluences $>2$ Jyms is $\alpha=-2.2^{+0.6}_{-1.2}$ and still consistent with a Euclidean distribution ($\alpha=-3/2$). We also find that the all-sky rate is $1.7^{+1.5}_{-0.9}\times10^3$FRBs/($4\pi$ sr)/day above $\sim2$ Jyms and there is currently no strong evidence for a latitude-dependent FRB sky-rate.

Fast radio bursts are a new class of transient radio phenomena currently detected as millisecond radio pulses with very high dispersion measures. As new radio surveys begin searching for FRBs a large population is expected to be detected in real-time, triggering a range of multi-wavelength and multi-messenger telescopes to search for repeating bursts and/or associated emission. Here we propose a method for disseminating FRB triggers using Virtual Observatory Events (VOEvents). This format was developed and is used successfully for transient alerts across the electromagnetic spectrum and for multi-messenger signals such as gravitational waves. In this paper we outline a proposed VOEvent standard for FRBs that includes the essential parameters of the event and where these parameters should be specified within the structure of the event. An additional advantage to the use of VOEvents for FRBs is that the events can automatically be ingested into the FRB Catalogue (FRBCAT) enabling real-time updates for public use. We welcome feedback from the community on the proposed standard outlined below and encourage those interested to join the nascent working group forming around this topic.

We present VLA images and HST/STIS spectra of sources within the center of the brightest cluster galaxy (BCG) in Abell 2261. These observations were obtained to test the hypothesis that its extremely large, flat core reflects the ejection of its supermassive black hole. Spectra of three of the four most luminous "knots" embedded in the core were taken to test whether one may represent stars bound to a displaced massive black hole. The three knots have radial velocity offsets dV < ~150 km/s from the BCG. Knots 2 and 3 show kinematics, colors, and stellar masses consistent with infalling low-mass galaxies or larger stripped cluster members. Large errors in the stellar velocity dispersion of Knot 1, however, mean that we cannot rule out the hypothesis that it hosts a high-mass black hole. A2261-BCG has a compact, relic radio-source offset by 6.5 kpc (projected) from the optical core's center, but no active radio core that would pinpoint the galaxy's central black hole to a tight 10 GHz flux limit <3.6 uJy. Its spectrum and morphology are suggestive of an active galactic nucleus that switched off >48 Myr ago, with an equipartition condition magnetic field of 15 uG. These observations are still consistent with the hypothesis that the nuclear black hole has been ejected from its core, but the critical task of locating the supermassive black hole or demonstrating that A2261-BCG lacks one remains to be done.

We undertook coordinated campaigns with the Green Bank, Effelsberg, and Arecibo radio telescopes during Chandra X-ray Observatory and XMM-Newton observations of the repeating fast radio burst FRB 121102 to search for simultaneous radio and X-ray bursts. We find 12 radio bursts from FRB 121102 during 70 ks total of X-ray observations. We detect no X-ray photons at the times of radio bursts from FRB 121102 and further detect no X-ray bursts above the measured background at any time. We place a 5$\sigma$ upper limit of $3\times10^{-11}$ erg cm$^{-2}$ on the 0.5--10 keV fluence for X-ray bursts at the time of radio bursts for durations $<700$ ms, which corresponds to a burst energy of $4\times10^{45}$ erg at the measured distance of FRB 121102. We also place limits on the 0.5--10 keV fluence of $5\times10^{-10}$ erg cm$^{-2}$ and $1\times10^{-9}$ erg cm$^{-2}$ for bursts emitted at any time during the XMM-Newton and Chandra observations, respectively, assuming a typical X-ray burst duration of 5 ms. We analyze data from the Fermi Gamma-ray Space Telescope Gamma-ray Burst Monitor and place a 5$\sigma$ upper limit on the 10--100 keV fluence of $4\times10^{-9}$ erg cm$^{-2}$ ($5\times10^{47}$ erg at the distance of FRB 121102) for gamma-ray bursts at the time of radio bursts. We also present a deep search for a persistent X-ray source using all of the X-ray observations taken to date and place a 5$\sigma$ upper limit on the 0.5--10 keV flux of $4\times10^{-15}$ erg s$^{-1}$ cm$^{-2}$ ($3\times10^{41}$ erg~s$^{-1}$ at the distance of FRB 121102). We discuss these non-detections in the context of the host environment of FRB 121102 and of possible sources of fast radio bursts in general.

We describe the Survey for Pulsars and Extragalactic Radio Bursts (SUPERB), an ongoing pulsar and fast transient survey using the Parkes radio telescope. SUPERB involves real-time acceleration searches for pulsars and single-pulse searches for pulsars and fast radio bursts. We report on the observational setup, data analysis, multi-wavelength/messenger connections, survey sensitivities to pulsars and fast radio bursts and the impact of radio frequency interference. We further report on the first 10 pulsars discovered in the project. Among these is PSR~J1306$-$40, a millisecond pulsar in a binary system where it appears to be eclipsed for a large fraction of the orbit. PSR~J1421$-$4407 is another binary millisecond pulsar; its orbital period is $30.7$ days. This orbital period is in a range where only highly eccentric binaries are known, and expected by theory; despite this its orbit has an eccentricity of $10^{-5}$.

Supermassive black hole binaries (SMBHBs) in the 10 million to 10 billion $M_\odot$ range form in galaxy mergers, and live in galactic nuclei with large and poorly constrained concentrations of gas and stars. There are currently no observations of merging SMBHBs--- it is in fact possible that they stall at their final parsec of separation and never merge. While LIGO has detected high frequency GWs, SMBHBs emit GWs in the nanohertz to millihertz band. This is inaccessible to ground-based interferometers, but possible with Pulsar Timing Arrays (PTAs). Using data from local galaxies in the 2 Micron All-Sky Survey, together with galaxy merger rates from Illustris, we find that there are on average $91\pm7$ sources emitting GWs in the PTA band, and $7\pm2$ binaries which will never merge. Local unresolved SMBHBs can contribute to GW background anisotropy at a level of $\sim20\%$, and if the GW background can be successfully isolated, GWs from at least one local SMBHB can be detected in 10 years.

We present optical, near- and mid-infrared imaging of the host galaxy of FRB 121102 with the Gemini North telescope, the Hubble Space Telescope and the Spitzer Space Telescope. The FRB 121102 host galaxy is resolved, revealing a bright star forming region located in the outskirts of the irregular, low-metallicity dwarf galaxy. The star forming region has a half-light radius of 0.68 kpc (0.20 arcsec), encompassing the projected location of the compact (<0.7 pc), persistent radio source that is associated with FRB 121102. The half-light diameter of the dwarf galaxy is 5 to 7 kpc, and broadband spectral energy distribution fitting indicates that it has a total stellar mass of M*~10^8 Msun. The metallicity of the host galaxy is low, 12+log10 ([O/H])=8.0+-0.1. The properties of the host galaxy of FRB 121102 are comparable to those of extreme emission line galaxies, also known to host hydrogen-poor superluminous supernovae and long-duration gamma-ray bursts. The projected location of FRB 121102 within the star forming region supports the proposed connection of FRBs with newly born neutron stars or magnetars.

We present results of the coordinated observing campaign that made the first subarcsecond localization of a Fast Radio Burst, FRB 121102. During this campaign, we made the first simultaneous detection of an FRB burst by multiple telescopes: the VLA at 3 GHz and the Arecibo Observatory at 1.4 GHz. Of the nine bursts detected by the Very Large Array at 3 GHz, four had simultaneous observing coverage at other observatories. We use multi-observatory constraints and modeling of bursts seen only at 3 GHz to confirm earlier results showing that burst spectra are not well modeled by a power law. We find that burst spectra are characterized by a ~500 MHz envelope and apparent radio energy as high as $10^{40}$ erg. We measure significant changes in the apparent dispersion between bursts that can be attributed to frequency-dependent profiles or some other intrinsic burst structure that adds a systematic error to the estimate of DM by up to 1%. We use FRB 121102 as a prototype of the FRB class to estimate a volumetric birth rate of FRB sources $R_{FRB} \approx 5x10^{-5}/N_r$ Mpc$^{-3}$ yr$^{-1}$, where $N_r$ is the number of bursts per source over its lifetime. This rate is broadly consistent with models of FRBs from young pulsars or magnetars born in superluminous supernovae or long gamma-ray bursts, if the typical FRB repeats on the order of thousands of times during its lifetime.

We report on the discovery of a new fast radio burst, FRB 150215, with the Parkes radio telescope on 2015 February 15. The burst was detected in real time with a dispersion measure (DM) of 1105.6$\pm$0.8 pc cm^{-3}, a pulse duration of 2.8^{+1.2}_{-0.5} ms, and a measured peak flux density assuming the burst was at beam center of 0.7^{+0.2}_{-0.1} Jy. The FRB originated at a Galactic longitude and latitude of 24.66^{\circ}, 5.28^{\circ}, 25 degrees away from the Galactic Center. The burst was found to be 43$\pm$5% linearly polarized with a rotation measure (RM) in the range -9 < RM < 12 rad m^{-2} (95% confidence level), consistent with zero. The burst was followed-up with 11 telescopes to search for radio, optical, X-ray, gamma-ray and neutrino emission. Neither transient nor variable emission was found to be associated with the burst and no repeat pulses have been observed in 17.25 hours of observing. The sightline to the burst is close to the Galactic plane and the observed physical properties of FRB 150215 demonstrate the existence of sight lines of anomalously low RM for a given electron column density. The Galactic RM foreground may approach a null value due to magnetic field reversals along the line of sight, a decreased total electron column density from the Milky Way, or some combination of these effects. A lower Galactic DM contribution might explain why this burst was detectable whereas previous searches at low latitude have had lower detection rates than those out of the plane.

The fast radio bursts (FRBs) 110220 and 140514 were detected at telescope pointing locations within 9 arcmin of each other over three years apart, both within the same 14.4 arcmin beam of the Parkes radio telescope. Nevertheless, they generally have not been considered to be from the same source because of a vastly different dispersion measure (DM) for the two bursts by over $380\,{\rm pc\,cm^{-3}}$. Here we consider the hypothesis that these two FRBs are from the same neutron star embedded within a supernova remnant (SNR) that provides an evolving DM as the ejecta expands and becomes more diffuse. Using such a model and the observed DM change, it can be argued that the corresponding SN must have occurred within $\approx10.2$ years of FRB 110220. Furthermore, constraints can be placed on the SN ejecta mass and explosion energy, which appear to require a stripped envelope (Type Ib/c) SN and/or a very energetic explosion. A third FRB from this location would be even more constraining, allowing the component of the DM due to the SNR to be separated from the unchanging DM components due to the host galaxy and intergalactic medium. In the future, if more FRBs are found to repeat, the sort of arguments presented here can be used to test the young neutron star progenitor hypothesis for FRBs.

The precise localization of the repeating fast radio burst (FRB 121102) has provided the first unambiguous association (chance coincidence probability $p\lesssim3\times10^{-4}$) of an FRB with an optical and persistent radio counterpart. We report on optical imaging and spectroscopy of the counterpart and find that it is an extended ($0.6^{\prime\prime}-0.8^{\prime\prime}$) object displaying prominent Balmer and [OIII] emission lines. Based on the spectrum and emission line ratios, we classify the counterpart as a low-metallicity, star-forming, $m_{r^\prime} = 25.1$ AB mag dwarf galaxy at a redshift of $z=0.19273(8)$, corresponding to a luminosity distance of 972 Mpc. From the angular size, the redshift, and luminosity, we estimate the host galaxy to have a diameter $\lesssim4$ kpc and a stellar mass of $M_*\sim4-7\times 10^{7}\,M_\odot$, assuming a mass-to-light ratio between 2 to 3$\,M_\odot\,L_\odot^{-1}$. Based on the H$\alpha$ flux, we estimate the star formation rate of the host to be $0.4\,M_\odot\,\mathrm{yr^{-1}}$ and a substantial host dispersion measure depth $\lesssim 324\,\mathrm{pc\,cm^{-3}}$. The net dispersion measure contribution of the host galaxy to FRB 121102 is likely to be lower than this value depending on geometrical factors. We show that the persistent radio source at FRB 121102's location reported by Marcote et al (2017) is offset from the galaxy's center of light by $\sim$200 mas and the host galaxy does not show optical signatures for AGN activity. If FRB 121102 is typical of the wider FRB population and if future interferometric localizations preferentially find them in dwarf galaxies with low metallicities and prominent emission lines, they would share such a preference with long gamma ray bursts and superluminous supernovae.

The millisecond-duration radio flashes known as Fast Radio Bursts (FRBs) represent an enigmatic astrophysical phenomenon. Recently, the sub-arcsecond localization (~ 100mas precision) of FRB121102 using the VLA has led to its unambiguous association with persistent radio and optical counterparts, and to the identification of its host galaxy. However, an even more precise localization is needed in order to probe the direct physical relationship between the millisecond bursts themselves and the associated persistent emission. Here we report very-long-baseline radio interferometric observations using the European VLBI Network and the 305-m Arecibo telescope, which simultaneously detect both the bursts and the persistent radio emission at milliarcsecond angular scales and show that they are co-located to within a projected linear separation of < 40pc (< 12mas angular separation, at 95% confidence). We detect consistent angular broadening of the bursts and persistent radio source (~ 2-4mas at 1.7GHz), which are both similar to the expected Milky Way scattering contribution. The persistent radio source has a projected size constrained to be < 0.7pc (< 0.2mas angular extent at 5.0GHz) and a lower limit for the brightness temperature of T_b > 5 x 10^7K. Together, these observations provide strong evidence for a direct physical link between FRB121102 and the compact persistent radio source. We argue that a burst source associated with a low-luminosity active galactic nucleus or a young neutron star energizing a supernova remnant are the two scenarios for FRB121102 that best match the observed data.

Fast radio bursts are astronomical radio flashes of unknown physical nature with durations of milliseconds. Their dispersive arrival times suggest an extragalactic origin and imply radio luminosities orders of magnitude larger than any other kind of known short-duration radio transient. Thus far, all FRBs have been detected with large single-dish telescopes with arcminute localizations, and attempts to identify their counterparts (source or host galaxy) have relied on contemporaneous variability of field sources or the presence of peculiar field stars or galaxies. These attempts have not resulted in an unambiguous association with a host or multi-wavelength counterpart. Here we report the sub-arcsecond localization of FRB 121102, the only known repeating burst source, using high-time-resolution radio interferometric observations that directly image the bursts themselves. Our precise localization reveals that FRB 121102 originates within 100 mas of a faint 180 uJy persistent radio source with a continuum spectrum that is consistent with non-thermal emission, and a faint (25th magnitude) optical counterpart. The flux density of the persistent radio source varies by tens of percent on day timescales, and very long baseline radio interferometry yields an angular size less than 1.7 mas. Our observations are inconsistent with the fast radio burst having a Galactic origin or its source being located within a prominent star-forming galaxy. Instead, the source appears to be co-located with a low-luminosity active galactic nucleus or a previously unknown type of extragalactic source. [Truncated] If other fast radio bursts have similarly faint radio and optical counterparts, our findings imply that direct sub-arcsecond localizations of FRBs may be the only way to provide reliable associations.

We extend profile domain pulsar timing to incorporate wide-band effects such as frequency-dependent profile evolution and broadband shape variation in the pulse profile. We also incorporate models for temporal variations in both pulse width and in the separation in phase of the main pulse and interpulse. We perform the analysis with both nested sampling and Hamiltonian Monte Carlo methods. In the latter case we introduce a new parameterisation of the posterior that is extremely efficient in the low signal-to-noise regime and can be readily applied to a wide range of scientific problems. We apply this methodology to a series of simulations, and to between seven and nine yr of observations for PSRs J1713$+$0747, J1744$-$1134, and J1909$-$3744 with frequency coverage that spans 700-3600MHz. We use a smooth model for profile evolution across the full frequency range, and compare smooth and piecewise models for the temporal variations in DM. We find the profile domain framework consistently results in improved timing precision compared to the standard analysis paradigm by as much as 40% for timing parameters. Incorporating smoothness in the DM variations into the model further improves timing precision by as much as 30%. For PSR J1713+0747 we also detect pulse shape variation uncorrelated between epochs, which we attribute to variation intrinsic to the pulsar at a level consistent with previously published analyses. Not accounting for this shape variation biases the measured arrival times at the level of $\sim$30ns, the same order of magnitude as the expected shift due to gravitational-waves in the pulsar timing band.

Statistical interpretation of sparsely sampled event rates has become vital for new transient surveys, particularly those aimed at detecting fast radio bursts (FRBs). We provide an accessible reference for a number of simple, but critical, statistical questions relevant for current transient and FRB research and utilizing the negative binomial model for counts in which the count rate parameter is uncertain or randomly biased from one study to the next. We apply these methods to re-assess and update results from previous FRB surveys, finding as follows. 1) Thirteen FRBs detected across five high-Galactic-latitude (> 30$^\circ$) surveys are highly significant $(p = 5\times 10^{-5})$ evidence of a higher rate relative to the single FRB detected across four low-latitude (< 5$^\circ$) surveys, even after accounting for effects that dampen Galactic plane sensitivity. High- vs. mid-latitude (5 to 15$^\circ$) is marginally significant $(p = 0.03)$. 2) A meta analysis of twelve heterogeneous surveys gives an FRB rate of 2866 sky$^{-1}$ day$^{-1}$ above 1 Jy at high Galactic latitude (95% confidence 1121 to 7328) and 285 sky$^{-1}$ day$^{-1}$ at low/mid latitudes (95% from 48 to 1701). 3) Using the Parkes HTRU high-latitude setup requires 193 observing hours to achieve 50% probability of detecting an FRB and 937 hours to achieve 95% probability, based on the ten detections of (Champion et al. 2016) and appropriately accounting for uncertainty in the unknown Poisson rate. 4) Two quick detections at Parkes from a small number of high-latitude fields (Ravi et al. 2015; Petroff et al. 2015) tentatively favor a look long survey style relative to the scan wide HTRU survey, but only at $p = 0.07$ significance.

Fast radio bursts (FRBs) are millisecond-duration events thought to originate beyond the Milky Way galaxy. Uncertainty surrounding the burst sources, and their propagation through intervening plasma, has limited their use as cosmological probes. We report on a mildly dispersed (dispersion measure 266.5+-0.1 pc cm^-3), exceptionally intense (120+-30 Jy), linearly polarized, scintillating burst (FRB 150807) that we directly localize to 9 arcmin^2. Based on a low Faraday rotation (12.0+-0.7 rad m^-2), we infer negligible magnetization in the circum-burst plasma and constrain the net magnetization of the cosmic web along this sightline to <21 nG, parallel to the line-of-sight. The burst scintillation suggests weak turbulence in the ionized intergalactic medium.

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.

Pulsar timing arrays are now setting increasingly tight limits on the gravitational wave background from binary supermassive black holes. But as upper limits grow more constraining, what can be implied about galaxy evolution? We investigate which astrophysical parameters have the largest impact on strain spectrum predictions and provide a simple framework to directly translate between measured values for the parameters of galaxy evolution and PTA limits on the gravitational wave background of binary supermassive black holes. We find that the most influential observable is the relation between a host galaxy's central bulge and its central black hole, $\mbox{$M_{\bullet}$-$M_{\rm bulge}$}$, which has the largest effect on the mean value of the characteristic strain amplitude. However, the variance of each prediction is dominated by uncertainties in the galaxy stellar mass function. Using this framework with the best published PTA limit, we can set limits on the shape and scatter of the $\mbox{$M_{\bullet}$-$M_{\rm bulge}$}$ relation. We find our limits to be in contention with strain predictions using two leading measurements of this relation. We investigate several possible reasons for this disagreement. If we take the $\mbox{$M_{\bullet}$-$M_{\rm bulge}$}$ relations to be correct within a simple power-law model for the gravitational wave background, then the inconsistency is reconcilable by allowing for an additional "stalling" time between a galaxy merger and evolution of a binary supermassive black hole to sub-parsec scales, with lower limits on this timescale of $\sim 1-2$ Gyr.

Quantum fluctuations of the gravitational field in the early Universe, amplified by inflation, produce a primordial gravitational-wave background across a broad frequency band. We derive constraints on the spectrum of this gravitational radiation, and hence on theories of the early Universe, by combining experiments that cover 29 orders of magnitude in frequency. These include Planck observations of cosmic microwave background temperature and polarization power spectra and lensing, together with baryon acoustic oscillations and big bang nucleosynthesis measurements, as well as new pulsar timing array and ground-based interferometer limits. While individual experiments constrain the gravitational-wave energy density in specific frequency bands, the combination of experiments allows us to constrain cosmological parameters, including the inflationary spectral index, $n_t$, and the tensor-to-scalar ratio, $r$. Results from individual experiments include the most stringent nanohertz limit of the primordial background to date from the Parkes Pulsar Timing Array, $\Omega_{\rm gw}(f)<2.3\times10^{-10}$. Observations of the cosmic microwave background alone limit the gravitational-wave spectral index at 95\% confidence to $n_t\lesssim5$ for a tensor-to-scalar ratio of $r = 0.11$. However, the combination of all the above experiments limits $n_t<0.36$. Future Advanced LIGO observations are expected to further constrain $n_t<0.34$ by 2020. When cosmic microwave background experiments detect a non-zero $r$, our results will imply even more stringent constraints on $n_t$ and hence theories of the early Universe.

In recent years, millisecond duration radio signals originating from distant galaxies appear to have been discovered in the so-called Fast Radio Bursts. These signals are dispersed according to a precise physical law and this dispersion is a key observable quantity which, in tandem with a redshift measurement, can be used for fundamental physical investigations. While every fast radio burst has a dispersion measurement, none before now have had a redshift measurement, due to the difficulty in pinpointing their celestial coordinates. Here we present the discovery of a fast radio burst and the identification of a fading radio transient lasting $\sim 6$ days after the event, which we use to identify the host galaxy; we measure the galaxy's redshift to be $z=0.492\pm0.008$. The dispersion measure and redshift, in combination, provide a direct measurement of the cosmic density of ionised baryons in the intergalactic medium of $\Omega_{\mathrm{IGM}}=4.9 \pm 1.3\%$, in agreement with the expectation from WMAP, and including all of the so-called "missing baryons". The $\sim6$-day transient is largely consistent with a short gamma-ray burst radio afterglow, and its existence and timescale do not support progenitor models such as giant pulses from pulsars, and supernovae. This contrasts with the interpretation of another recently discovered fast radio burst, suggesting there are at least two classes of bursts.

We analyse the stochastic properties of the 49 pulsars that comprise the first International Pulsar Timing Array (IPTA) data release. We use Bayesian methodology, performing model selection to determine the optimal description of the stochastic signals present in each pulsar. In addition to spin-noise and dispersion-measure (DM) variations, these models can include timing noise unique to a single observing system, or frequency band. We show the improved radio-frequency coverage and presence of overlapping data from different observing systems in the IPTA data set enables us to separate both system and band-dependent effects with much greater efficacy than in the individual PTA data sets. For example, we show that PSR J1643$-$1224 has, in addition to DM variations, significant band-dependent noise that is coherent between PTAs which we interpret as coming from time-variable scattering or refraction in the ionised interstellar medium. Failing to model these different contributions appropriately can dramatically alter the astrophysical interpretation of the stochastic signals observed in the residuals. In some cases, the spectral exponent of the spin noise signal can vary from 1.6 to 4 depending upon the model, which has direct implications for the long-term sensitivity of the pulsar to a stochastic gravitational-wave (GW) background. By using a more appropriate model, however, we can greatly improve a pulsar's sensitivity to GWs. For example, including system and band-dependent signals in the PSR J0437$-$4715 data set improves the upper limit on a fiducial GW background by $\sim 60\%$ compared to a model that includes DM variations and spin-noise only.

The highly stable spin of neutron stars can be exploited for a variety of (astro-)physical investigations. In particular arrays of pulsars with rotational periods of the order of milliseconds can be used to detect correlated signals such as those caused by gravitational waves. Three such "Pulsar Timing Arrays" (PTAs) have been set up around the world over the past decades and collectively form the "International" PTA (IPTA). In this paper, we describe the first joint analysis of the data from the three regional PTAs, i.e. of the first IPTA data set. We describe the available PTA data, the approach presently followed for its combination and suggest improvements for future PTA research. Particular attention is paid to subtle details (such as underestimation of measurement uncertainty and long-period noise) that have often been ignored but which become important in this unprecedentedly large and inhomogeneous data set. We identify and describe in detail several factors that complicate IPTA research and provide recommendations for future pulsar timing efforts. The first IPTA data release presented here (and available online) is used to demonstrate the IPTA's potential of improving upon gravitational-wave limits placed by individual PTAs by a factor of ~2 and provides a 2-sigma limit on the dimensionless amplitude of a stochastic GWB of 1.7x10^{-15} at a frequency of 1 yr^{-1}. This is 1.7 times less constraining than the limit placed by (Shannon et al. 2015), due mostly to the more recent, high-quality data they used.

We present high-precision timing observations spanning up to nine years for 37 millisecond pulsars monitored with the Green Bank and Arecibo radio telescopes as part of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) project. We describe the observational and instrumental setups used to collect the data, and methodology applied for calculating pulse times of arrival; these include novel methods for measuring instrumental offsets and characterizing low signal-to-noise ratio timing results. The time of arrival data are fit to a physical timing model for each source, including terms that characterize time-variable dispersion measure and frequency-dependent pulse shape evolution. In conjunction with the timing model fit, we have performed a Bayesian analysis of a parameterized timing noise model for each source, and detect evidence for excess low-frequency, or "red," timing noise in 10 of the pulsars. For 5 of these cases this is likely due to interstellar medium propagation effects rather than intrisic spin variations. Subsequent papers in this series will present further analysis of this data set aimed at detecting or limiting the presence of nanohertz-frequency gravitational wave signals.

We have begun an exciting era for gravitational wave detection, as several world-leading experiments are breaching the threshold of anticipated signal strengths. Pulsar timing arrays (PTAs) are pan-Galactic gravitational wave detectors that are already cutting into the expected strength of gravitational waves from cosmic strings and binary supermassive black holes in the nHz-$\mu$Hz gravitational wave band. These limits are leading to constraints on the evolutionary state of the Universe. Here, we provide a broad review of this field, from how pulsars are used as tools for detection, to astrophysical sources of uncertainty in the signals PTAs aim to see, to the primary current challenge areas for PTA work. This review aims to provide an up-to-date reference point for new parties interested in the field of gravitational wave detection via pulsar timing.

The detection of five new fast radio bursts (FRBs) found in the High Time Resolution Universe high latitude survey is presented. The rate implied is 6$^{+4}_{-3}\times~10^3$ (95%) FRBs sky$^{-1}$ day$^{-1}$ above a fluence of between 0.13 and 5.9 Jy ms for FRBs between 0.128 and 262 ms in duration. One of these FRBs has a clear two-component profile, each component is similar to the known population of single component FRBs and are separated by 2.4(4) ms. All the FRB components appear to be unresolved following deconvolution with a scattering tail and accounting for intra-channel smearing. The two-component FRB also has the highest dispersion measure (1629 pc cm$^{-3}$) of any FRB to-date. Many of the proposed models to explain FRBs use a single high energy event involving compact objects (such as neutron star mergers) and therefore cannot easily explain a two-component FRB. Models that are based on extreme versions of flaring, pulsing or orbital events however could produce multiple component profiles. The compatibility of these models and the FRB rate implied by these detections is discussed.

By regularly monitoring the most stable millisecond pulsars over many years, pulsar timing arrays (PTAs) are positioned to detect and study correlations in the timing behaviour of those pulsars. Gravitational waves (GWs) from supermassive black hole binaries (SMBHBs) are an exciting potentially detectable source of such correlations. We describe a straight-forward technique by which a PTA can be "phased-up" to form time series of the two polarisation modes of GWs coming from a particular direction of the sky. Our technique requires no assumptions regarding the time-domain behaviour of a GW signal. This method has already been used to place stringent bounds on GWs from individual SMBHBs in circular orbits. Here, we describe the methodology and demonstrate the versatility of the technique in searches for a wide variety of GW signals including bursts with unmodeled waveforms. Using the first six years of data from the Parkes Pulsar Timing Array, we conduct an all-sky search for a detectable excess of GW power from any direction. For the lines of sight to several nearby massive galaxy clusters, we carry out a more detailed search for GW bursts with memory, which are distinct signatures of SMBHB mergers. In all cases, we find that the data are consistent with noise.

We report here on key science topics for the Next Generation Very Large Array in the areas of time domain, fundamental physics, and cosmology. Key science cases considered are pulsars in orbit around the Galactic Center massive black hole, Sagittarius A*, electromagnetic counterparts to gravitational waves, and astrometric cosmology. These areas all have the potential for ground-breaking and transformative discovery. Numerous other topics were discussed during the preparation of this report and some of those discussions are summarized here, as well. There is no doubt that further investigation of the science case will reveal rich and compelling opportunities.

We present timing models for 20 millisecond pulsars in the Parkes Pulsar Timing Array. The precision of the parameter measurements in these models has been improved over earlier results by using longer data sets and modelling the non-stationary noise. We describe a new noise modelling procedure and demonstrate its effectiveness using simulated data. Our methodology includes the addition of annual dispersion measure (DM) variations to the timing models of some pulsars. We present the first significant parallax measurements for PSRs J1024-0719, J1045-4509, J1600-3053, J1603-7202, and J1730-2304, as well as the first significant measurements of some post-Keplerian orbital parameters in six binary pulsars, caused by kinematic effects. Improved Shapiro delay measurements have resulted in much improved pulsar mass measurements, particularly for PSRs J0437-4715 and J1909-3744 with $M_p=1.44\pm0.07$ $M_\odot$ and $M_p=1.47\pm0.03$ $M_\odot$ respectively. The improved orbital period-derivative measurement for PSR J0437-4715 results in a derived distance measurement at the 0.16% level of precision, $D=156.79\pm0.25$ pc, one of the most fractionally precise distance measurements of any star to date.

Gravitational waves are expected to be radiated by supermassive black hole binaries formed during galaxy mergers. A stochastic superposition of gravitational waves from all such binary systems will modulate the arrival times of pulses from radio pulsars. Using observations of millisecond pulsars obtained with the Parkes radio telescope, we constrain the characteristic amplitude of this background, $A_{\rm c,yr}$, to be < $1.0\times10^{-15}$ with 95% confidence. This limit excludes predicted ranges for $A_{\rm c,yr}$ from current models with 91-99.7% probability. We conclude that binary evolution is either stalled or dramatically accelerated by galactic-center environments, and that higher-cadence and shorter-wavelength observations would result in an increased sensitivity to gravitational waves.

Several theories exist to explain the source of the bright, millisecond duration pulses known as fast radio bursts (FRBs). If the progenitors of FRBs are non-cataclysmic, such as giant pulses from pulsars, pulsar-planet binaries, or magnetar flares, FRB emission may be seen to repeat. We have undertaken a survey of the fields of eight known FRBs from the High Time Resolution Universe survey to search for repeating pulses. Although no repeat pulses were detected the survey yielded the detection of a new FRB, described in Petroff et al. (2015a). From our observations we rule out periodic repeating sources with periods P $\leq$ 8.6 hours and rule out sources with periods 8.6 < P < 21 hours at the 90% confidence level. At P $\geq$ 21 hours our limits fall off as ~1/P. Dedicated and persistent observations of FRB source fields are needed to rule out repetition on longer timescales, a task well-suited to next generation wide-field transient detectors.

We compute upper limits on the nanohertz-frequency isotropic stochastic gravitational wave background (GWB) using the 9-year data release from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) collaboration. We set upper limits for a GWB from supermassive black hole binaries under power law, broken power law, and free spectral coefficient GW spectrum models. We place a 95\% upper limit on the strain amplitude (at a frequency of yr$^{-1}$) in the power law model of $A_{\rm gw} < 1.5\times 10^{-15}$. For a broken power law model, we place priors on the strain amplitude derived from simulations of Sesana (2013) and McWilliams et al. (2014). We find that the data favor a broken power law to a pure power law with odds ratios of 22 and 2.2 to one for the McWilliams and Sesana prior models, respectively. The McWilliams model is essentially ruled out by the data, and the Sesana model is in tension with the data under the assumption of a pure power law. Using the broken power-law analysis we construct posterior distributions on environmental factors that drive the binary to the GW-driven regime including the stellar mass density for stellar-scattering, mass accretion rate for circumbinary disk interaction, and orbital eccentricity for eccentric binaries, marking the first time that the shape of the GWB spectrum has been used to make astrophysical inferences. We then place the most stringent limits so far on the energy density of relic GWs, $\Omega_\mathrm{gw}(f)\,h^2 < 4.2 \times 10^{-10}$, yielding a limit on the Hubble parameter during inflation of $H_*=1.6\times10^{-2}~m_{Pl}$, where $m_{Pl}$ is the Planck mass. Our limit on the cosmic string GWB, $\Omega_\mathrm{gw}(f)\, h^2 < 2.2 \times 10^{-10}$, translates to a conservative limit of $G\mu<3.3\times 10^{-8}$ - a factor of 4 better than the joint Planck and high-$l$ CMB data from other experiments.

Extreme scattering events (ESEs) in the interstellar medium (ISM) were first observed in regular flux measurements of compact extragalactic sources. They are characterized by a flux variation over a period of weeks, suggesting the passage of a "diverging plasma lens" across the line of sight. Modeling the refraction of such a lens indicates that the structure size must be of order AU and the electron density of order 10s of cm^{-3}. Similar structures have been observed in measurements of pulsar intensity scintillation and group delay. Here we report observations of two ESEs showing increases in both intensity scintillation and dispersion made with the Parkes Pulsar Timing Array (PPTA). These allow us to make more complete models of the ESE, including an estimate of the "outer-scale" of the turbulence in the plasma lens. These observations show clearly that the ESE structure is fully turbulent on an AU scale. They provide some support for the idea that the structures are extended along the line of sight, such as would be the case for a scattering shell. The dispersion measurements also show a variety of AU scale structures which would not be called ESEs, yet involve electron density variations typical of ESEs and likely have the same origin.

We present initial results from the low-latitude Galactic plane region of the High Time Resolution Universe pulsar survey conducted at the Parkes 64-m radio telescope. We discuss the computational challenges arising from the processing of the terabyte-sized survey data. Two new radio interference mitigation techniques are introduced, as well as a partially-coherent segmented acceleration search algorithm which aims to increase our chances of discovering highly-relativistic short-orbit binary systems, covering a parameter space including potential pulsar-black hole binaries. We show that under a constant acceleration approximation, a ratio of data length over orbital period of ~0.1 results in the highest effectiveness for this search algorithm. From the 50 per cent of data processed thus far, we have re-detected 435 previously known pulsars and discovered a further 60 pulsars, two of which are fast-spinning pulsars with periods less than 30ms. PSR J1101-6424 is a millisecond pulsar whose heavy white dwarf (WD) companion and short spin period of 5.1ms indicate a rare example of full-recycling via Case A Roche lobe overflow. PSR J1757-27 appears to be an isolated recycled pulsar with a relatively long spin period of 17ms. In addition, PSR J1244-6359 is a mildly-recycled binary system with a heavy WD companion, PSR J1755-25 has a significant orbital eccentricity of 0.09, and PSR J1759-24 is likely to be a long-orbit eclipsing binary with orbital period of the order of tens of years. Comparison of our newly-discovered pulsar sample to the known population suggests that they belong to an older population. Furthermore, we demonstrate that our current pulsar detection yield is as expected from population synthesis.

"Perytons" are millisecond-duration transients of terrestrial origin, whose frequency-swept emission mimics the dispersion of an astrophysical pulse that has propagated through tenuous cold plasma. In fact, their similarity to FRB 010724 had previously cast a shadow over the interpretation of "fast radio bursts," which otherwise appear to be of extragalactic origin. Until now, the physical origin of the dispersion-mimicking perytons had remained a mystery. We have identified strong out-of-band emission at 2.3--2.5 GHz associated with several peryton events. Subsequent tests revealed that a peryton can be generated at 1.4 GHz when a microwave oven door is opened prematurely and the telescope is at an appropriate relative angle. Radio emission escaping from microwave ovens during the magnetron shut-down phase neatly explain all of the observed properties of the peryton signals. Now that the peryton source has been identified, we furthermore demonstrate that the microwaves on site could not have caused FRB 010724. This and other distinct observational differences show that FRBs are excellent candidates for genuine extragalactic transients.

We present high signal-to-noise ratio, multi-frequency polarization pulse profiles for 24 millisecond pulsars that are being observed as part of the Parkes Pulsar Timing Array (PPTA) project. The pulsars are observed in three bands, centred close to 730, 1400 and 3100 MHz, using a dual-band 10 cm/50 cm receiver and the central beam of the 20 cm multibeam receiver. Observations spanning approximately six years have been carefully calibrated and summed to produce high S/N profiles. This allows us to study the individual profile components and in particular how they evolve with frequency. We also identify previously undetected profile features. For many pulsars we show that pulsed emission extends across almost the entire pulse profile. The pulse component widths and component separations follow a complex evolution with frequency; in some cases these parameters increase and in other cases they decrease with increasing frequency. The evolution with frequency of the polarization properties of the profile is also non-trivial. We provide evidence that the pre- and post-cursors generally have higher fractional linear polarization than the main pulse. We have obtained the spectral index and rotation measure for each pulsar by fitting across all three observing bands. For the majority of pulsars, the spectra follow a single power-law and the position angles follow a $\lambda^2$ relation, as expected. However, clear deviations are seen for some pulsars. We also present phase-resolved measurements of the spectral index, fractional linear polarization and rotation measure. All these properties are shown to vary systematically over the pulse profile.

Among efforts to detect gravitational radiation, pulsar timing arrays are uniquely poised to detect "memory" signatures, permanent perturbations in spacetime from highly energetic astrophysical events such as mergers of supermassive black hole binaries. The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) observes dozens of the most stable millisecond pulsars using the Arecibo and Green Bank radio telescopes in an effort to study, among other things, gravitational wave memory. We herein present the results of a search for gravitational wave bursts with memory (BWMs) using the first five years of NANOGrav observations. We develop original methods for dramatically speeding up searches for BWM signals. In the directions of the sky where our sensitivity to BWMs is best, we would detect mergers of binaries with reduced masses of $10^9$ $M_\odot$ out to distances of 30 Mpc; such massive mergers in the Virgo cluster would be marginally detectable. We find no evidence for BWMs. However, with our non-detection, we set upper limits on the rate at which BWMs of various amplitudes could have occurred during the time spanned by our data--e.g., BWMs with amplitudes greater than $10^{-13}$ must occur at a rate less than 1.5 yr$^{-1}$.

We report on the first millisecond timescale radio interferometric search for the new class of transient known as fast radio bursts (FRBs). We used the Very Large Array (VLA) for a 166-hour, millisecond imaging campaign to detect and precisely localize an FRB. We observed at 1.4 GHz and produced visibilities with 5 ms time resolution over 256 MHz of bandwidth. Dedispersed images were searched for transients with dispersion measures from 0 to 3000 pc/cm3. No transients were detected in observations of high Galactic latitude fields taken from September 2013 though October 2014. Observations of a known pulsar show that images typically had a thermal-noise limited sensitivity of 120 mJy/beam (8 sigma; Stokes I) in 5 ms and could detect and localize transients over a wide field of view. Our nondetection limits the FRB rate to less than 7e4/sky/day (95% confidence) above a fluence limit of 1.2 Jy-ms. Assuming a Euclidean flux distribution, the VLA rate limit is inconsistent with the published rate of Thornton et al. We recalculate previously published rates with a homogeneous consideration of the effects of primary beam attenuation, dispersion, pulse width, and sky brightness. This revises the FRB rate downward and shows that the VLA observations had a roughly 60% chance of detecting a typical FRB and that a 95% confidence constraint would require roughly 500 hours of similar VLA observing. Our survey also limits the repetition rate of an FRB to 2 times less than any known repeating millisecond radio transient.

Fast radio bursts (FRBs) are one of the most tantalizing mysteries of the radio sky; their progenitors and origins remain unknown and until now no rapid multiwavelength follow-up of an FRB has been possible. New instrumentation has decreased the time between observation and discovery from years to seconds, and enables polarimetry to be performed on FRBs for the first time. We have discovered an FRB (FRB 140514) in real-time on 14 May, 2014 at 17:14:11.06 UTC at the Parkes radio telescope and triggered follow-up at other wavelengths within hours of the event. FRB 140514 was found with a dispersion measure (DM) of 562.7(6) cm$^{-3}$ pc, giving an upper limit on source redshift of $z \lesssim 0.5$. FRB 140514 was found to be 21$\pm$7% (3-$\sigma$) circularly polarized on the leading edge with a 1-$\sigma$ upper limit on linear polarization $<10%$. We conclude that this polarization is intrinsic to the FRB. If there was any intrinsic linear polarization, as might be expected from coherent emission, then it may have been depolarized by Faraday rotation caused by passing through strong magnetic fields and/or high density environments. FRB 140514 was discovered during a campaign to re-observe known FRB fields, and lies close to a previous discovery, FRB 110220; based on the difference in DMs of these bursts and time-on-sky arguments, we attribute the proximity to sampling bias and conclude that they are distinct objects. Follow-up conducted by 12 telescopes observing from X-ray to radio wavelengths was unable to identify a variable multiwavelength counterpart, allowing us to rule out models in which FRBs originate from nearby ($z < 0.3$) supernovae and long duration gamma-ray bursts.

We present the discovery of a further five recycled pulsar systems in the mid-Galactic latitude portion of the High Time Resolution Universe (HTRU) Survey. The pulsars have rotational periods ranging from 2 ms to 66 ms, and four are in binary systems with orbital periods between 10.8 hours and 9.0 days. Three of these binary systems are particularly interesting; PSR J1227-6208 has a pulse period of 34.5 ms and the highest mass function of all pulsars with near-circular orbits. The circular orbit suggests that the companion is not another neutron star, so future timing experiments may reveal one of the heaviest white dwarfs ever found ($>$ 1.3 M$_\odot$). Timing observations of PSR J1431$-$4715 indicate that it is eclipsed by its companion which has a mass indicating it belongs to the redback class of eclipsing millisecond pulsars. PSR J1653-2054 has a companion with a minimum mass of only $0.08$ M$_\odot$, placing it among the class of pulsars with low-mass companions. Unlike the majority of such systems, however, no evidence of eclipses is seen at 1.4 GHz.

Anisotropic bursts of gravitational radiation produced by events such as super-massive black hole mergers leave permanent imprints on space. Such gravitational wave "memory" (GWM) signals are, in principle, detectable through pulsar timing as sudden changes in the apparent pulse frequency of a pulsar. If an array of pulsars is monitored as a GWM signal passes over the Earth, the pulsars would simultaneously appear to change pulse frequency by an amount that varies with their sky position in a quadrupolar fashion. Here we describe a search algorithm for such events and apply the algorithm to approximately six years of data from the Parkes Pulsar Timing Array. We find no GWM events and set an upper bound on the rate for events which could have been detected. We show, using simple models of black hole coalescence rates, that this non-detection is not unexpected.

We present results of an all-sky search in the Parkes Pulsar Timing Array (PPTA) Data Release 1 data set for continuous gravitational waves (GWs) in the frequency range from $5\times 10^{-9}$ to $2\times 10^{-7}$ Hz. Such signals could be produced by individual supermassive binary black hole systems in the early stage of coalescence. We phase up the pulsar timing array data set to form, for each position on the sky, two data streams that correspond to the two GW polarizations and then carry out an optimal search for GW signals on these data streams. Since no statistically significant GWs were detected, we place upper limits on the intrinsic GW strain amplitude $h_0$ for a range of GW frequencies. For example, at $10^{-8}$ Hz our analysis has excluded with $95\%$ confidence the presence of signals with $h_0\geqslant 1.7\times 10^{-14}$. Our new limits are about a factor of four more stringent than those of Yardley et al. (2010) based on an earlier PPTA data set and a factor of two better than those reported in the recent Arzoumanian et al. (2014) paper. We also present PPTA directional sensitivity curves and find that for the most sensitive region on the sky, the current data set is sensitive to GWs from circular supermassive binary black holes with chirp masses of $10^{9} M_{\odot}$ out to a luminosity distance of about 100 Mpc. Finally, we set an upper limit of $4 \times 10^{-3} {\rm{Mpc}}^{-3} {\rm{Gyr}}^{-1}$ at $95\%$ confidence on the coalescence rate of nearby ($z \lesssim 0.1$) supermassive binary black holes in circular orbits with chirp masses of $10^{10}M_{\odot}$.

We report the detection of a dispersed Fast Radio Burst (FRB) in archival intermediate-latitude Parkes Radio Telescope data. The burst appears to be of the same physical origin as the four purported extragalactic FRBs reported by Thornton et al. (2013). This burst's arrival time precedes the Thornton et al.~bursts by ten years. We consider that this survey, and many other archival low-latitude (|gb|<30deg) pulsar surveys, have been searched for FRBs but produced fewer detections than the comparatively brief Thornton et al.~search. Such a rate dependence on Galactic position could provide critical supporting evidence for an extragalactic origin for FRBs. To test this, we form an analytic expression to account for Galactic position and survey setup in FRB rate predictions. Employing a sky temperature, scattering, and dispersion model of the Milky Way, we compute the expected number of FRBs if they are isotropically distributed on the sky w.r.t. Galactic position (i.e. local), and if they are of extragalactic origin. We demonstrate that the relative detection rates reject a local origin with a confidence of 99.96% (~3.6 sigma). The extragalactic predictions provide a better agreement, however are still strong discrepancies with the low-latitude detection rate at a confidence of 99.69% (~2.9 sigma). However, for the extragalactic population, the differences in predicted vs.~detected population may be accounted for by a number of factors, which we discuss.

High-sensitivity radio-frequency observations of millisecond pulsars usually show stochastic, broadband, pulse-shape variations intrinsic to the pulsar emission process. These variations induce jitter noise in pulsar timing observations; understanding the properties of this noise is of particular importance for the effort to detect gravitational waves with pulsar timing arrays. We assess the short-term profile and timing stability of 22 millisecond pulsars that are part of the Parkes Pulsar Timing Array sample by examining intra-observation arrival time variability and single-pulse phenomenology. In 7 of the 22 pulsars, in the band centred at approximately 1400MHz, we find that the brightest observations are limited by intrinsic jitter. We find consistent results, either detections or upper limits, for jitter noise in other frequency bands. PSR J1909-3744 shows the lowest levels of jitter noise, which we estimate to contribute $\sim$10 ns root mean square error to the arrival times for hour-duration observations. Larger levels of jitter noise are found in pulsars with wider pulses and distributions of pulse intensities. The jitter noise in PSR J0437-4715 decorrelates over a bandwidth of $\sim$2 GHz. We show that the uncertainties associated with timing pulsar models can be improved by including physically motivated jitter uncertainties. Pulse-shape variations will limit the timing precision at future, more sensitive, telescopes; it is imperative to account for this noise when designing instrumentation and timing campaigns for these facilities.

Fast radio bursts (FRBs) are an emerging class of bright, highly dispersed radio pulses. Recent work by Thornton et al. (2013) has revealed a population of FRBs in the High Time Resolution Universe (HTRU) survey at high Galactic latitudes. A variety of progenitors have been proposed including cataclysmic events at cosmological distances, Galactic flare stars, and terrestrial radio frequency interference. Here we report on a search for FRBs at intermediate Galactic latitudes ($-15^{\circ}$ $< b <$ 15$^{\circ}$) in data taken as part of the HTRU survey. No FRBs were discovered in this region. Several effects such as dispersion, scattering, sky temperature and scintillation decrease the sensitivity by more than 3$\sigma$ in $\sim$20\% of survey pointings. Including all of these effects, we exclude the hypothesis that FRBs are uniformly distributed on the sky with 99\% confidence. This low probability implies that additional factors -- not accounted for by standard Galactic models -- must be included to ease the discrepancy between the detection rates at high and low Galactic latitudes. A revised rate estimate or another strong and heretofore unknown selection effect in Galactic latitude would provide closer agreement between the surveys' detection rates. The dearth of detections at low Galactic latitude disfavors a Galactic origin for these bursts.

The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) project currently observes 43 pulsars using the Green Bank and Arecibo radio telescopes. In this work we use a subset of 17 pulsars timed for a span of roughly five years (2005--2010). We analyze these data using standard pulsar timing models, with the addition of time-variable dispersion measure and frequency-variable pulse shape terms. Within the timing data, we perform a search for continuous gravitational waves from individual supermassive black hole binaries in circular orbits using robust frequentist and Bayesian techniques. We find that there is no evidence for the presence of a detectable continuous gravitational wave; however, we can use these data to place the most constraining upper limits to date on the strength of such gravitational waves. Using the full 17 pulsar dataset we place a 95% upper limit on the sky-averaged strain amplitude of $h_0\lesssim 3.8\times 10^{-14}$ at a frequency of 10 nHz. Furthermore, we place 95% \emph{all sky} lower limits on the luminosity distance to such gravitational wave sources finding that the $d_L \gtrsim 425$ Mpc for sources at a frequency of 10 nHz and chirp mass $10^{10}{\rm M}_{\odot}$. We find that for gravitational wave sources near our best timed pulsars in the sky, the sensitivity of the pulsar timing array is increased by a factor of $\sim$4 over the sky-averaged sensitivity. Finally we place limits on the coalescence rate of the most massive supermassive black hole binaries.

This whitepaper describes how the VLASS could be designed in a manner to allow the identification of candidate dual active galactic nuclei (AGN) at separations <7 kpc. Dual AGN represent a clear marker of two supermassive black holes within an ongoing merger. A dual AGN survey will provide a wealth of studies in structure growth and gravitational-wave science. Radio wavelengths are ideal for identifying close pairs, as disturbed stellar and gaseous material can obscure their presence in optical and shorter wavelengths. With sufficiently high resolution and sensitivity, a large-scale radio imaging survey like the VLASS will uncover many of these systems and provide the means to broadly study the radio properties of candidate dual systems revealed at other wavelengths. We determine that the ideal survey for our purposes will be at as high a resolution as possible, with significantly more science return in A array at L-band or higher, or B array at C-band or higher. We describe a range of potential survey parameters within this document. Based on the analysis outlined in this whitepaper, our ideal survey would create a catalogue of $\gtrsim$100 dual AGN in either: 1) a medium-sensitivity (~1 mJy detection threshold), wide-field (few thousand square degree) survey, or 2) a high-sensitivity (~10 uJy threshold) survey of several hundred square degrees.

We report on the discovery of four millisecond pulsars (MSPs) in the High Time Resolution Universe (HTRU) pulsar survey being conducted at the Parkes 64-m radio telescope. All four MSPs are in binary systems and are likely to have white dwarf companions. In addition, we present updated timing solutions for 12 previously published HTRU MSPs, revealing new observational parameters such as five proper motion measurements and significant temporal dispersion measure variations in PSR J1017-7156. We discuss the case of PSR J1801-3210, which shows no significant period derivative after four years of timing data. Our best-fit solution shows a period derivative of the order of $10^{-23}$, an extremely small number compared to that of a typical MSP. However, it is likely that the pulsar lies beyond the Galactic Centre, and an unremarkable intrinsic period derivative is reduced to close to zero by the Galactic potential acceleration. Furthermore, we highlight the potential to employ PSR J1801-3210 in the strong equivalence principle test due to its wide and circular orbit. In a broader comparison with the known MSP population, we suggest a correlation between higher mass functions and the presence of eclipses in `very low-mass binary pulsars', implying that eclipses are observed in systems with high orbital inclinations. We also suggest that the distribution of the total mass of binary systems is inversely-related to the Galactic height distribution. Finally, we report on the first detection of PSRs J1543-5149 and J1811-2404 as gamma-ray pulsars.

The formation and growth processes of supermassive black holes (SMBHs) are not well constrained. SMBH population models, however, provide specific predictions for the properties of the gravitational-wave background (GWB) from binary SMBHs in merging galaxies throughout the Universe. Using observations from the Parkes Pulsar Timing Array, we constrain the fractional GWB energy density with 95% confidence to be ${\Omega}_{GW}(H_0/73 {\rm km} {\rm s}^{-1} {\rm Mpc}^{-1})^2 < 1.3 \times 10^{-9}$ at a frequency of 2.8 nHz, which is approximately a factor of six more stringent than previous limits. We compare our limit to models of the SMBH population and find inconsistencies at confidence levels between 46% and 91%. For example, the standard galaxy formation model implemented in the Millennium simulations is inconsistent with our limit with 50% probability.

We present a polarimetric analysis of 49 long-period pulsars discovered as part of the High Time Resolution Universe (HTRU) southern survey. The sources exhibit the typical characteristics of "old" pulsars, with low fractional linear and circular polarisation and narrow, multicomponent profiles. Although the position angle swings are generally complex, for two of the analysed pulsars (J1622-3751 and J1710-2616) we obtained an indication of the geometry via the rotating vector model. We were able to determine a value of the rotation measure (RM) for 34 of the sources which, when combined with their dispersion measures (DM), yields an integrated magnetic field strength along the line of sight. With the data presented here, the total number of values of RM associated to pulsars discovered during the HTRU southern survey sums to 51. The RMs are not consistent with the hypothesis of a counter-clockwise direction of the Galactic magnetic field within an annulus included between 4 and 6 kpc from the Galactic centre. A partial agreement with a counter-clockwise sense of the Galactic magnetic field within the spiral arms is however found in the area of the Carina-Sagittarius arm.

We present polarisation data for 180 extragalactic sources extracted from the Australia Telescope 20 GHz (AT20G) survey catalog, and observed with the Australia Telescope Compact Array during a dedicated, high sensitivity run. For the sake of completeness we extracted the polarisation information for 7 extended sources from the 9-yr WMAP coadded maps at 23 GHz. The full sample of 187 sources constitutes a 99% complete sample of extragalactic sources brighter than S(20 GHz)=500mJy at the selection epoch with declination <-30 deg. The sample has a 91.4% detection rate in polarisation at 20GHz (94% if considering the sub-sample of point like sources). We have measurements also at 4.8 and 8.6 GHz within 1 month of the 20GHz observations for 172 sources to reconstruct the spectral properties of the sample in total intensity and in polarisation: 143 of them have a polarisation detection at all three frequencies. We find that there is no statistically significant evidence of a relationship either between the fraction of polarisation and frequency or between the fraction of polarisation and the total intensity flux density. This indicates that Faraday depolarisation is not very important above 4.8 GHz and that the magnetic field is not substantially more ordered in the regions dominating the emission at higher frequencies (up to 20 GHz). We estimate the distribution of the polarisation fraction and the polarised flux density source counts at 20GHz.

Recent years have seen a burgeoning interest in using pulsar timing arrays (PTAs) as gravitational-wave (GW) detectors. To date, that interest has focused mainly on three particularly promising source types: supermassive--black-hole binaries, cosmic strings, and the stochastic background from early-Universe phase transitions. In this paper, by contrast, our aim is to investigate the PTA potential for discovering unanticipated sources. We derive significant constraints on the available discovery space based solely on energetic and statistical considerations: we show that a PTA detection of GWs at frequencies above ~3.e-5 Hz would either be an extraordinary coincidence or violate "cherished beliefs;" we show that for PTAs GW memory can be more detectable than direct GWs, and that, as we consider events at ever higher redshift, the memory effect increasingly dominates an event's total signal-to-noise ratio. The paper includes also a simple analysis of the effects of pulsar red noise in PTA searches, and a demonstration that the effects of periodic GWs in the 10^-8 -- 10^-4.5 Hz band would not be degenerate with small errors in standard pulsar parameters (except in a few narrow bands).

Pulsar timing arrays are sensitive to gravitational waves from supermassive black hole (SMBH) binaries at orbital separations of << 1pc. There is currently an observational paucity of such systems, although they are central figures in studies of galaxy evolution, merger dynamics, and active nucleus formation. We review the prospects of detecting SMBH binaries through electromagnetic radiative processes thought to be associated with galaxy mergers and late-stage binary evolution. We then discuss the scientific goals of joint pulsar timing and electromagnetic studies of these systems, including the facilitation of binary parameter estimation, identifying galactic hosts of gravitational wave emitters, and relevant studies of merger dynamics and cosmology. The use of upcoming high-precision timing arrays with the International Pulsar Timing Array and the Square Kilometre Array, combined with ongoing electromagnetic observing campaigns to identify active SMBH binaries, provide generous possibilities for multi-messenger astrophysics in the near future.

This paper presents the discovery and timing parameters for five millisecond pulsars (MSPs), four in binary systems with probable white dwarf companions and one isolated, found in ongoing processing of the High Time Resolution Universe Pulsar Survey (HTRU). We also present high quality polarimetric data on four of them. These further discoveries confirm the high potential of our survey in finding pulsars with very short spin periods. At least two of these five MSPs are excellent candidates to be included in the Pulsar Timing Array projects. Thanks to the wealth of MSP discoveries in the HTRU survey, we revisit the question of whether the luminosity distributions of isolated and binary MSPs are different. Using the Cordes and Lazio distance model and our new and catalogue flux density measurements, we find that 41 of the 42 most luminous MSPs in the Galactic disk are in binaries and a statistical analysis suggests that the luminosity functions differ with 99.9% significance. We conclude that the formation process that leads to solitary MSPs affects their luminosities, despite their period and period derivatives being similar to those of pulsars in binary systems.

Searches for transient astrophysical sources often reveal unexpected classes of objects that are useful physical laboratories. In a recent survey for pulsars and fast transients we have uncovered four millisecond-duration radio transients all more than 40{\deg} from the Galactic plane. The bursts' properties indicate that they are of celestial rather than terrestrial origin. Host galaxy and intergalactic medium models suggest that they have cosmological redshifts of 0.5 to 1, and distances of up to 3 gigaparsecs. No temporally coincident x- or gamma-ray signature was identified in association with the bursts. Characterization of the source population and identification of host galaxies offers an opportunity to determine the baryonic content of the Universe.

We have used millisecond pulsars (MSPs) from the southern High Time Resolution Universe (HTRU) intermediate latitude survey area to simulate the distribution and total population of MSPs in the Galaxy. Our model makes use of the scale factor method, which estimates the ratio of the total number of MSPs in the Galaxy to the known sample. Using our best fit value for the z-height, z=500 pc, we find an underlying population of MSPs of 8.3(\pm 4.2)*10^4 sources down to a limiting luminosity of L_min=0.1 mJy kpc^2 and a luminosity distribution with a steep slope of d\log N/d\log L = -1.45(\pm 0.14). However, at the low end of the luminosity distribution, the uncertainties introduced by small number statistics are large. By omitting very low luminosity pulsars, we find a Galactic population above L_min=0.2 mJy kpc^2 of only 3.0(\pm 0.7)*10^4 MSPs. We have also simulated pulsars with periods shorter than any known MSP, and estimate the maximum number of sub-MSPs in the Galaxy to be 7.8(\pm 5.0)*10^4 pulsars at L=0.1 mJy kpc^2. In addition, we estimate that the high and low latitude parts of the southern HTRU survey will detect 68 and 42 MSPs respectively, including 78 new discoveries. Pulsar luminosity, and hence flux density, is an important input parameter in the model. Some of the published flux densities for the pulsars in our sample do not agree with the observed flux densities from our data set, and we have instead calculated average luminosities from archival data from the Parkes Telescope. We found many luminosities to be very different than their catalogue values, leading to very different population estimates. Large variations in flux density highlight the importance of including scintillation effects in MSP population studies.

Next-generation radio arrays, including the Square Kilometre Array (SKA) and its pathfinders, will open up new avenues for exciting transient science at radio wavelengths. Their innovative designs, comprising a large number of small elements, pose several challenges in digital processing and optimal observing strategies. The Giant Metrewave Radio Telescope (GMRT) presents an excellent test-bed for developing and validating suitable observing modes and strategies for transient experiments with future arrays. Here we describe the first phase of the ongoing development of a transient detection system for GMRT that is planned to eventually function in a commensal mode with other observing programs. It capitalizes on the GMRT's interferometric and sub-array capabilities, and the versatility of a new software backend. We outline considerations in the plan and design of transient exploration programs with interferometric arrays, and describe a pilot survey that was undertaken to aid in the development of algorithms and associated analysis software. This survey was conducted at 325 and 610 MHz, and covered 360 deg$^2$ of the sky with short dwell times. It provides large volumes of real data that can be used to test the efficacies of various algorithms and observing strategies applicable for transient detection. We present examples that illustrate the methodologies of detecting short-duration transients, including the use of sub-arrays for higher resilience to spurious events of terrestrial origin, localization of candidate events via imaging and the use of a phased array for improved signal detection and confirmation. In addition to demonstrating applications of interferometric arrays for fast transient exploration, our efforts mark important steps in the roadmap toward SKA-era science.

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.

Six years ago, the discovery of Rotating Radio Transients (RRATs) marked what appeared to be a new type of sparsely-emitting pulsar. Since 2006, more than 70 of these objects have been discovered in single-pulse searches of archival and new surveys. With a continual inflow of new information about the RRAT population in the form of new discoveries, multi-frequency follow-ups, coherent timing solutions, and pulse rate statistics, a view is beginning to form of the place in the pulsar population RRATs hold. Here we review the properties of neutron stars discovered through single pulse searches. We first seek to clarify the definition of the term RRAT, emphasising that "the RRAT population" encompasses several phenomenologies. A large subset of RRATs appear to represent the tail of an extended distribution of pulsar nulling fractions and activity cycles; these objects present several key open questions remaining in this field.

Signals from radio pulsars show a wavelength-dependent delay due to dispersion in the interstellar plasma. At a typical observing wavelength, this delay can vary by tens of microseconds on five-year time scales, far in excess of signals of interest to pulsar timing arrays, such as that induced by a gravitational-wave background. Measurement of these delay variations is not only crucial for the detection of such signals, but also provides an unparallelled measurement of the turbulent interstellar plasma at au scales. In this paper we demonstrate that without consideration of wavelength- independent red-noise, 'simple' algorithms to correct for interstellar dispersion can attenuate signals of interest to pulsar timing arrays. We present a robust method for this correction, which we validate through simulations, and apply it to observations from the Parkes Pulsar Timing Array. Correction for dispersion variations comes at a cost of increased band-limited white noise. We discuss scheduling to minimise this additional noise, and factors, such as scintillation, that can exacerbate the problem. Comparison with scintillation measurements confirms previous results that the spectral exponent of electron density variations in the interstellar medium often appears steeper than expected. We also find a discrete change in dispersion measure of PSR J1603-7202 of ~2x10^{-3} cm^{-3}pc for about 250 days. We speculate that this has a similar origin to the 'extreme scattering events' seen in other sources. In addition, we find that four pulsars show a wavelength-dependent annual variation, indicating a persistent gradient of electron density on an au spatial scale, which has not been reported previously.

A "pulsar timing array" (PTA), in which observations of a large sample of pulsars spread across the celestial sphere are combined, allows investigation of "global" phenomena such as a background of gravitational waves or instabilities in atomic timescales that produce correlated timing residuals in the pulsars of the array. The Parkes Pulsar Timing Array (PPTA) is an implementation of the PTA concept based on observations with the Parkes 64-m radio telescope. A sample of 20 millisecond pulsars is being observed at three radio-frequency bands, 50cm (~700 MHz), 20cm (~1400 MHz) and 10cm (~3100 MHz), with observations at intervals of 2 - 3 weeks. Regular observations commenced in early 2005. This paper describes the systems used for the PPTA observations and data processing, including calibration and timing analysis. The strategy behind the choice of pulsars, observing parameters and analysis methods is discussed. Results are presented for PPTA data in the three bands taken between 2005 March and 2011 March. For ten of the 20 pulsars, rms timing residuals are less than 1 microsec for the best band after fitting for pulse frequency and its first time derivative. Significant "red" timing noise is detected in about half of the sample. We discuss the implications of these results on future projects including the International Pulsar Timing Array (IPTA) and a PTA based on the Square Kilometre Array. We also present an "extended PPTA" data set that combines PPTA data with earlier Parkes timing data for these pulsars.

We present 18 GHz Australia Telescope Compact Array imaging of the Mpc-scale quasar jet PKS 0637-752 with angular resolution ~0.58 arcseconds. We draw attention to a spectacular train of quasi-periodic knots along the inner 11 arcseconds of the jet, with average separation ~1.1 arcsec (7.6 kpc projected). We consider two classes of model to explain the periodic knots: those that involve a static pattern through which the jet plasma travels (e.g. stationary shocks); and those that involve modulation of the jet engine. Interpreting the knots as re-confinement shocks implies the jet kinetic power Q ~ 10^{46} erg/s, but the constant knot separation along the jet is not expected in a realistic external density profile. For models involving modulation of the jet engine, we find that the required modulation period is 2 x 10^3 yr < \tau < 3 x 10^5 yr. The lower end of this range is applicable if the jet remains highly relativistic on kpc-scales, as implied by the IC/CMB model of jet X-ray emission. We suggest that the quasi-periodic jet structure in PKS 0637-752 may be analogous to the quasi-periodic jet modulation seen in the microquasar GRS 1915+105, believed to result from limit cycle behaviour in an unstable accretion disk. If variations in the accretion rate are driven by a binary black hole, the predicted orbital radius is 0.7 < a < 30 pc, which corresponds to a maximum angular separation of ~0.1 - 5 mas.

We present 75 pulsars discovered in the mid-latitude portion of the High Time Resolution Universe survey, 54 of which have full timing solutions. All the pulsars have spin periods greater than 100 ms, and none of those with timing solutions are in binaries. Two display particularly interesting behaviour; PSR J1054-5944 is found to be an intermittent pulsar, and PSR J1809-0119 has glitched twice since its discovery. In the second half of the paper we discuss the development and application of an artificial neural network in the data-processing pipeline for the survey. We discuss the tests that were used to generate scores and find that our neural network was able to reject over 99% of the candidates produced in the data processing, and able to blindly detect 85% of pulsars. We suggest that improvements to the accuracy should be possible if further care is taken when training an artificial neural network; for example ensuring that a representative sample of the pulsar population is used during the training process, or the use of different artificial neural networks for the detection of different types of pulsars.

Using observations of pulsars from the Parkes Pulsar Timing Array (PPTA) project we develop the first pulsar-based timescale that has a precision comparable to the uncertainties in international atomic timescales. Our ensemble of pulsars provides an Ensemble Pulsar Scale (EPS) analogous to the free atomic timescale Echelle Atomique Libre (EAL). The EPS can be used to detect fluctuations in atomic timescales and therefore can lead to a new realisation of Terrestrial Time, TT(PPTA11). We successfully follow features known to affect the frequency of the International Atomic Timescale (TAI) and we find marginally significant differences between TT(PPTA11) and TT(BIPM11). We discuss the various phenomena that lead to a correlated signal in the pulsar timing residuals and therefore limit the stability of the pulsar timescale.

Here we report on observations of the radio magnetar PSR J1622-4950 at frequencies from 1.4 to 17 GHz. We show that although its flux density is varying up to a factor of ~10 within a few days, it has on average decreased by a factor of 2 over the last 700 days. At the same time, timing analysis indicates a trend of decreasing spin-down rate over our entire data set, again of about a factor of 2 over 700 days, but also an erratic variability in the spin-down rate within this time span. Integrated pulse profiles are often close to 100 per cent linearly polarized, but large variations in both the profile shape and fractional polarization are regularly observed. Furthermore, the behaviour of the position angle of the linear polarization is very complex - offsets in both the absolute position angle and the phase of the position angle sweep are often seen and the occasional presence of orthogonal mode jumps further complicates the picture. However, model fitting indicates that the magnetic and rotation axes are close to aligned. Finally, a single pulse analysis has been carried out at four observing frequencies, demonstrating that the wide pulse profile is built up of narrow spikes of emission, with widths that scale inversely with observing frequency. All three of the known radio magnetars seem to have similar characteristics, with highly polarized emission, time-variable flux density and pulse profiles, and with spectral indices close to zero.

We report on the pulse-to-pulse energy distributions and phase-resolved modulation properties for catalogued pulsars in the southern High Time Resolution Universe intermediate-latitude survey. We selected the 315 pulsars detected in a single-pulse search of this survey, allowing a large sample unbiased regarding any rotational parameters of neutron stars. We found that the energy distribution of many pulsars is well-described by a log-normal distribution, with few deviating from a small range in log-normal scale and location parameters. Some pulsars exhibited multiple energy states corresponding to mode changes, and implying that some observed "nulling" may actually be a mode-change effect. PSRJ1900-2600 was found to emit weakly in its previously-identified "null" state. We found evidence for another state-change effect in two pulsars, which show bimodality in their nulling time scales; that is, they switch between a continuous-emission state and a single-pulse-emitting state. Large modulation occurs in many pulsars across the full integrated profile, with increased sporadic bursts at leading and trailing sub-beam edges. Some of these high-energy outbursts may indicate the presence of "giant pulse" phenomena. We found no correlation with modulation and pulsar period, age, or other parameters. Finally, the deviation of integrated pulse energy from its average value was generally quite small, despite the significant phase-resolved modulation in some pulsars; we interpret this as tenuous evidence of energy regulation between distinct pulsar sub-beams.

We present multi-wavelength observations of the radio magnetar PSR J1622-4950 and its environment. Observations of PSR J1622-4950 with Chandra (in 2007 and 2009) and XMM (in 2011) show that the X-ray flux of PSR J1622-4950 has decreased by a factor of ~50 over 3.7 years, decaying exponentially with a characteristic time of 360 +/- 11 days. This behavior identifies PSR J1622-4950 as a possible addition to the small class of transient magnetars. The X-ray decay likely indicates that PSR J1622-4950 is recovering from an X-ray outburst that occurred earlier in 2007, before the 2007 Chandra observations. Observations with the Australia Telescope Compact Array show strong radio variability, including a possible radio flaring event at least one and a half years after the 2007 X-ray outburst that may be a direct result of this X-ray event. Radio observations with the Molonglo Observatory Synthesis Telescope reveal that PSR J1622-4950 is 8' southeast of a diffuse radio arc, G333.9+0.0, which appears non-thermal in nature and which could possibly be a previously undiscovered supernova remnant. If G333.9+0.0 is a supernova remnant then the estimates of its size and age, combined with the close proximity and reasonable implied velocity of PSR J1622-4950, suggests that these two objects could be physically associated.

We present an analysis of high-precision pulsar timing data taken as part of the North American Nanohertz Observatory for Gravitational waves (NANOGrav) project. We have observed 17 pulsars for a span of roughly five years using the Green Bank and Arecibo radio telescopes. We analyze these data using standard pulsar timing models, with the addition of time-variable dispersion measure and frequency-variable pulse shape terms. Sub-microsecond timing residuals are obtained in nearly all cases, and the best root-mean-square timing residuals in this set are ~30-50 ns. We present methods for analyzing post-fit timing residuals for the presence of a gravitational wave signal with a specified spectral shape. These optimally take into account the timing fluctuation power removed by the model fit, and can be applied to either data from a single pulsar, or to a set of pulsars to detect a correlated signal. We apply these methods to our dataset to set an upper limit on the strength of the nHz-frequency stochastic supermassive black hole gravitational wave background of h_c (1 yr^-1) < 7x10^-15 (95%). This result is dominated by the timing of the two best pulsars in the set, PSRs J1713+0747 and J1909-3744.

Preliminary specifications for the Square Kilometre Array (SKA) call for 25% of the total collecting area of the dish array to be located at distances greater than 180 km from the core, with a maximum baseline of at least 3000 km. The array will provide angular resolution ~ 40 - 2 mas at 0.5 - 10 GHz with image sensitivity reaching < 50 nJy/beam in an 8 hour integration with 500 MHz bandwidth. Given these specifications, the high angular resolution component of the SKA will be capable of detecting brightness temperatures < 200 K with milliarcsecond-scale angular resolution. The aim of this article is to bring together in one place a discussion of the broad range of new and important high angular resolution science that will be enabled by the SKA, and in doing so, address the merits of long baselines as part of the SKA. We highlight the fact that high angular resolution requiring baselines greater than 1000 km provides a rich science case with projects from many areas of astrophysics, including important contributions to key SKA science.

We present the discovery of six millisecond pulsars (MSPs) in the High Time Resolution Universe (HTRU) survey for pulsars and fast transients carried out with the Parkes radio telescope. All six are in binary systems with approximately circular orbits and are likely to have white dwarf companions. PSR J1017-7156 has a high flux density and a narrow pulse width, making it ideal for precision timing experiments. PSRs J1446-4701 and J1125-5825 are coincident with gamma-ray sources, and folding the high-energy photons with the radio timing ephemeris shows evidence of pulsed gamma-ray emission. PSR J1502-6752 has a spin period of 26.7 ms, and its low period derivative implies that it is a recycled pulsar. The orbital parameters indicate it has a very low mass function, and therefore a companion mass much lower than usually expected for such a mildly recycled pulsar. In addition we present polarisation profiles for all 12 MSPs discovered in the HTRU survey to date. Similar to previous observations of MSPs, we find that many have large widths and a wide range of linear and circular polarisation fractions. Their polarisation profiles can be highly complex, and although the observed position angles often do not obey the rotating vector model, we present several examples of those that do. We speculate that the emission heights of MSPs are a substantial fraction of the light cylinder radius in order to explain broad emission profiles, which then naturally leads to a large number of cases where emission from both poles is observed.

Millisecond pulsars are thought to be neutron stars that have been spun-up by accretion of matter from a binary companion. Although most are in binary systems, some 30% are solitary, and their origin is therefore mysterious. PSR J1719-1438, a 5.7 ms pulsar, was detected in a recent survey with the Parkes 64m radio telescope. We show that it is in a binary system with an orbital period of 2.2 h. Its companion's mass is near that of Jupiter, but its minimum density of 23 g cm$^{-3}$ suggests that it may be an ultra-low mass carbon white dwarf. This system may thus have once been an Ultra Compact Low-Mass X-ray Binary, where the companion narrowly avoided complete destruction.

We report on variations in the mean position angle of the 20 millisecond pulsars being observed as part of the Parkes Pulsar Timing Array (PPTA) project. It is found that the observed variations are dominated by changes in the Faraday rotation occurring in the Earth's ionosphere. Two ionospheric models are used to correct for the ionospheric contribution and it is found that one based on the International Reference Ionosphere gave the best results. Little or no significant long-term variation in interstellar RM was found with limits typically about 0.1 rad m$^{-2}$ yr$^{-1}$ in absolute value. In a few cases, apparently significant RM variations over timescales of a few 100 days or more were seen. These are unlikely to be due to localised magnetised regions crossing the line of sight since the implied magnetic fields are too high. Most probably they are statistical fluctuations due to random spatial and temporal variations in the interstellar electron density and magnetic field along the line of sight.

We present the search methods and initial results for transient radio signals in the High Time Resolution Universe (HTRU) Survey. The HTRU survey's single-pulse search, the software designed to perform the search, and a determination of the HTRU survey's sensitivity to single pulses are described. Initial processing of a small fraction of the survey has produced 11 discoveries, all of which are sparsely-emitting neutron stars, as well as provided confirmation of two previously unconfirmed neutron stars. Most of the newly discovered objects lie in regions surveyed previously, indicating both the improved sensitivity of the HTRU survey observing system and the dynamic nature of the radio sky. The cycles of active and null states in nulling pulsars, rotating radio transients (RRATs), and long-term intermittent pulsars are explored in the context of determining the relationship between these populations, and of the sensitivity of a search to the various radio-intermittent neutron star populations. This analysis supports the case that many RRATs are in fact high-null-fraction pulsars (i. e. with null fraction >~0.95), and indicates that intermittent pulsars appear distinct from nulling pulsars in their activity cycle timescales. We find that in the measured population, there is a deficit of pulsars with typical emission timescales greater than ~300 s that is not readily explained by selection effects. The HTRU low-latitude survey will be capable of addressing whether this deficit is physical. We predict that the HTRU survey will explore pulsars with a broad range of nulling fractions (up to and beyond 0.999), and at its completion is likely to increase the currently known RRATs by a factor of more than two.

We search for the signature of an isotropic stochastic gravitational-wave background in pulsar timing observations using a frequency-domain correlation technique. These observations, which span roughly 12 yr, were obtained with the 64-m Parkes radio telescope augmented by public domain observations from the Arecibo Observatory. A wide range of signal processing issues unique to pulsar timing and not previously presented in the literature are discussed. These include the effects of quadratic removal, irregular sampling, and variable errors which exacerbate the spectral leakage inherent in estimating the steep red spectrum of the gravitational-wave background. These observations are found to be consistent with the null hypothesis, that no gravitational-wave background is present, with 76 percent confidence. We show that the detection statistic is dominated by the contributions of only a few pulsars because of the inhomogeneity of this data set. The issues of detecting the signature of a gravitational-wave background with future observations are discussed.

Polarization profiles are presented for 20 millisecond pulsars that are being observed as part of the Parkes Pulsar Timing Array project. The observations used the Parkes multibeam receiver with a central frequency of 1369 MHz and the Parkes digital filterbank pulsar signal-processing system PDFB2. Because of the large total observing time, the summed polarization profiles have very high signal/noise ratios and show many previously undetected profile features. Thirteen of the 20 pulsars show emission over more than half of the pulse period. Polarization variations across the profiles are complex and the observed position angle variations are generally not in accord with the rotating-vector model for pulsar polarization. Never-the-less, the polarization properties are broadly similar to those of normal (non-millisecond) pulsars, suggesting that the basic radio emission mechanism is the same in both classes of pulsar. The results support the idea that radio emission from millisecond pulsars originates high in the pulsar magnetosphere, probably close to the emission regions for high-energy X-ray and gamma-ray emission. Rotation measures were obtained for all 20 pulsars, eight of which had no previously published measurements.

We present the discovery of 5 millisecond pulsars found in the mid-Galactic latitude portion of the High Time Resolution Universe (HTRU) Survey. The pulsars have rotational periods from ~2.3 to ~7.5 ms, and all are in binary systems with orbital periods ranging from ~0.3 to ~150 d. In four of these systems, the most likely companion is a white dwarf, with minimum masses of ~0.2 Solar Masses. The other pulsar, J1731-1847, has a very low mass companion and exhibits eclipses, and is thus a member of the "black widow" class of pulsar binaries. These eclipses have been observed in bands centred near frequencies of 700, 1400 and 3000 MHz, from which measurements have been made of the electron density in the eclipse region. These measurements have been used to examine some possible eclipse mechanisms. The eclipse and other properties of this source are used to perform a comparison with the other known eclipsing and "black widow" pulsars. These new discoveries occupy a short-period and high-dispersion measure (DM) region of parameter space, which we demonstrate is a direct consequence of the high time and frequency resolution of the HTRU survey. The large implied distances to our new discoveries makes observation of their companions unlikely with both current optical telescopes and the Fermi Gamma-ray Space Telescope. The extremely circular orbits make any advance of periastron measurements highly unlikely. No relativistic Shapiro delays are obvious in any of the systems, although the low flux densities would make their detection difficult unless the orbits were fortuitously edge-on.

The Australia Telescope 20 GHz (AT20G) survey is a blind survey of the whole Southern sky at 20 GHz with follow-up observations at 4.8, 8.6, and 20 GHz carried out with the Australia Telescope Compact Array (ATCA). In this paper we present an analysis of radio spectral properties in total intensity and polarisation, sizes, optical identifications, and redshifts of the sample of the 5808 extragalactic sources in the survey catalogue of confirmed sources over the whole Southern sky excluding the strip at Galactic latitude |b|<1.5deg. The sample has a flux density limit of 40 mJy. Completeness has been measured as a function of scan region and flux density. Averaging over the whole survey area the follow-up survey is 78% complete above 50mJy and 93% complete above 100mJy. 3332 sources with declination <-15deg have good quality almost simultaneous observations at 4.8, 8.6, and 20GHz. The spectral analysis shows that the sample is dominated by flat-spectrum sources. The fraction of flat-spectrum sources decreases from 81% for 20GHz flux densities S>500mJy, to 60% for S<100mJy. There is also a clear spectral steepening at higher frequencies with the median spectral index decreasing from -0.16 between 4.8 and 8.6GHz to -0.28 between 8.6 and 20GHz. Simultaneous observations in polarisation are available for all the sources at all the frequencies. 768 sources have a good quality detection of polarised flux density at 20GHz; 467 of them were also detected in polarisation at 4.8 and/or at 8.6GHz so that it has been possible to compare the spectral behaviour in total intensity and polarisation. We have found that the polarised fraction increases slightly with frequency and decreases with flux density. Cross matches and comparisons have been made with other catalogues at lower radio frequencies, and in the optical, X-ray and gamma-ray bands. Redshift estimates are available for 825 sources.

We have embarked on a survey for pulsars and fast transients using the 13-beam Multibeam receiver on the Parkes radio telescope. Installation of a digital backend allows us to record 400 MHz of bandwidth for each beam, split into 1024 channels and sampled every 64 us. Limits of the receiver package restrict us to a 340 MHz observing band centred at 1352 MHz. The factor of eight improvement in frequency resolution over previous multibeam surveys allows us to probe deeper into the Galactic plane for short duration signals such as the pulses from millisecond pulsars. We plan to survey the entire southern sky in 42641 pointings, split into low, mid and high Galactic latitude regions, with integration times of 4200, 540 and 270 s respectively. Simulations suggest that we will discover 400 pulsars, of which 75 will be millisecond pulsars. With ~30% of the mid-latitude survey complete, we have re-detected 223 previously known pulsars and discovered 27 pulsars, 5 of which are millisecond pulsars. The newly discovered millisecond pulsars tend to have larger dispersion measures than those discovered in previous surveys, as expected from the improved time and frequency resolution of our instrument.

Three years ago, the report of a solitary radio burst was thought to be the first discovery of a rare, impulsive event of unknown extragalactic origin (Lorimer et al. 2007). The extragalactic interpretation was based on the swept-frequency nature of the event, which followed the dispersive delay expected from an extragalactic pulse. We report here on the detection of 16 pulses, the bulk of which exhibit a frequency sweep with a shape and magnitude resembling the Lorimer Burst. These new events were detected in a sidelobe of the Parkes Telescope and are of clearly terrestrial origin, with properties unlike any known sources of terrestrial broad-band radio emission. The new detections cast doubt on the extragalactic interpretation of the original burst, and call for further sophistication in radio-pulse survey techniques to identify the origin of the anomalous terrestrial signals and definitively distinguish future extragalactic pulse detections from local signals. The ambiguous origin of these seemingly dispersed, swept-frequency signals suggest that radio-pulse searches using multiple detectors will be the only experiments able to provide definitive information about the origin of new swept-frequency radio burst detections.

We present the sensitivity of the Parkes Pulsar Timing Array to gravitational waves emitted by individual super-massive black-hole binary systems in the early phases of coalescing at the cores of merged galaxies. Our analysis includes a detailed study of the effects of fitting a pulsar timing model to non-white timing residuals. Pulsar timing is sensitive at nanoHertz frequencies and hence complementary to LIGO and LISA. We place a sky-averaged constraint on the merger rate of nearby ($z < 0.6$) black-hole binaries in the early phases of coalescence with a chirp mass of $10^{10}\,\rmn{M}_\odot$ of less than one merger every seven years. The prospects for future gravitational-wave astronomy of this type with the proposed Square Kilometre Array telescope are discussed.

Using archival VLBI data for 3114 radio-luminous active galactic nuclei, we searched for binary supermassive black holes using a radio spectral index mapping technique which targets spatially resolved, double radio-emitting nuclei. Only one source was detected as a double nucleus. This result is compared with a cosmological merger rate model and interpreted in terms of (1) implications for post-merger timescales for centralisation of the two black holes, (2) implications for the possibility of "stalled" systems, and (3) the relationship of radio activity in nuclei to mergers. Our analysis suggests that the binary evolution of paired supermassive black holes (both of masses >= 1e8 Msun) spends less than 500 Myr in progression from the merging of galactic stellar cores to within the purported stalling radius for supermassive black hole pairs. The data show no evidence for an excess of stalled binary systems at small separations. We see circumstantial evidence that the relative state of radio emission between paired supermassive black holes is correlated within orbital separations of 2.5 kpc.

High-precision pulsar timing relies on a solar-system ephemeris in order to convert times of arrival (TOAs) of pulses measured at an observatory to the solar system barycenter. Any error in the conversion to the barycentric TOAs leads to a systematic variation in the observed timing residuals; specifically, an incorrect planetary mass leads to a predominantly sinusoidal variation having a period and phase associated with the planet's orbital motion about the Sun. By using an array of pulsars (PSRs J0437-4715, J1744-1134, J1857+0943, J1909-3744), the masses of the planetary systems from Mercury to Saturn have been determined. These masses are consistent with the best-known masses determined by spacecraft observations, with the mass of the Jovian system, 9.547921(2)E-4 Msun, being significantly more accurate than the mass determined from the Pioneer and Voyager spacecraft, and consistent with but less accurate than the value from the Galileo spacecraft. While spacecraft are likely to produce the most accurate measurements for individual solar system bodies, the pulsar technique is sensitive to planetary system masses and has the potential to provide the most accurate values of these masses for some planets.

As part of a survey for radio pulsars with the Parkes 64-m telescope we have discovered PSR J1622-4950, a pulsar with a 4.3-s rotation period. Follow-up observations show that the pulsar has the highest inferred surface magnetic field of the known radio pulsars (B ~ 3e14 G), exhibits significant timing noise and appears to have an inverted spectrum. Unlike the vast majority of the known pulsar population, PSR J1622-4950 appears to switch off for many hundreds of days and even in its on-state exhibits extreme variability in its flux density. Furthermore, the integrated pulse profile changes shape with epoch. All of these properties are remarkably similar to the only two magnetars previously known to emit radio pulsations. The position of PSR J1622-4950 is coincident with an X-ray source that, unlike the other radio pulsating magnetars, was found to be in quiescence. We conclude that our newly discovered pulsar is a magnetar - the first to be discovered via its radio emission.

The Parkes Pulsar Timing Array project aims to make a direct detection of a gravitational-wave background through timing of millisecond pulsars. In this article, the main requirements for that endeavour are described and recent and ongoing progress is outlined. We demonstrate that the timing properties of millisecond pulsars are adequate and that technological progress is timely to expect a successful detection of gravitational waves within a decade, or alternatively to rule out all current predictions for gravitational wave backgrounds formed by supermassive black-hole mergers.

The International Pulsar Timing Array project combines observations of pulsars from both Northern and Southern hemisphere observatories with the main aim of detecting ultra-low frequency (~10^-9 to 10^-8 Hz) gravitational waves. Here we introduce the project, review the methods used to search for gravitational waves emitted from coalescing supermassive binary black-hole systems in the centres of merging galaxies and discuss the status of the project.

We present the results of a search for transient radio bursts of between 0.125 and 32 millisecond duration in two archival pulsar surveys of intermediate galactic latitudes with the Parkes multibeam receiver. Fourteen new neutron stars have been discovered, seven of which belong to the recently identified "rotating radio transients" (RRATs) class. Here we describe our search methodology, and discuss the new detections in terms of how the RRAT population relates to the general population of pulsars. The new detections indicate (1) that the galactic z-distribution of RRATs in the surveys closely resembles the distribution of pulsars, with objects up to 0.86 kpc from the galactic plane; (2) where measurable, the RRAT pulse widths are similar to that of individual pulses from pulsars of similar period, implying a similar beaming fraction; and (3) our new detections span a variety of nulling fractions, and thus we postulate that the RRATs may simply be nulling pulsars that are only "on" for less than a pulse period. Finally, the newly discovered object PSR J0941-39 may represent a link between pulsars and RRATs. This bizarre object was discovered as an RRAT, but in follow-up observations often appeared as a bright (~10 mJy) pulsar with a low nulling fraction. It is obvious therefore that a neutron star can oscillate between being an RRAT and a pulsar. Crucially, the sites of the RRAT pulses are coincident with the pulsar's emission, implying that the two emission mechanisms are linked, and that RRATs are not just pulsars observed from different orientations.

We present the full source catalogue from the Australia Telescope 20 GHz (AT20G) Survey. The AT20G is a blind radio survey carried out at 20 GHz with the Australia Telescope Compact Array (ATCA) from 2004 to 2008, and covers the whole sky south of declination 0 deg. The AT20G source catalogue presented here is an order of magnitude larger than any previous catalogue of high-frequency radio sources, and includes 5890 sources above a 20 GHz flux-density limit of 40 mJy. All AT20G sources have total intensity and polarisation measured at 20 GHz, and most sources south of declination -15 deg also have near-simultaneous flux-density measurements at 5 and 8 GHz. A total of 1559 sources were detected in polarised total intensity at one or more of the three frequencies. We detect a small but significant population of non-thermal sources that are either undetected or have only weak detections in low-frequency catalogues. We introduce the term Ultra-Inverted Spectrum (UIS) to describe these radio sources, which have a spectral index alpha(5, 20) > +0.7 and which constitute roughly 1.2 per cent of the AT20G sample. The 20 GHz flux densities measured for the strongest AT20G sources are in excellent agreement with the WMAP 5-year source catalogue of Wright et al. (2009), and we find that the WMAP source catalogue is close to complete for sources stronger than 1.5 Jy at 23 GHz.

Analysis of high-precision timing observations of an array of approx. 20 millisecond pulsars (a so-called "timing array") may ultimately result in the detection of a stochastic gravitational-wave background. The feasibility of such a detection and the required duration of this type of experiment are determined by the achievable rms of the timing residuals and the timing stability of the pulsars involved. We present results of the first long-term, high-precision timing campaign on a large sample of millisecond pulsars used in gravitational-wave detection projects. We show that the timing residuals of most pulsars in our sample do not contain significant low-frequency noise that could limit the use of these pulsars for decade-long gravitational-wave detection efforts. For our most precisely timed pulsars, intrinsic instabilities of the pulsars or the observing system are shown to contribute to timing irregularities on a five-year timescale below the 100 ns level. Based on those results, realistic sensitivity curves for planned and ongoing timing array efforts are determined. We conclude that prospects for detection of a gravitational-wave background through pulsar timing array efforts within five years to a decade are good.

The PULSE@Parkes project has been designed to monitor the rotation of radio pulsars over time spans of days to years. The observations are obtained using the Parkes 64-m and 12-m radio telescopes by Australian and international high school students. These students learn the basis of radio astronomy and undertake small projects with their observations. The data are fully calibrated and obtained with the state-of-the-art pulsar hardware available at Parkes. The final data sets are archived and are currently being used to carry out studies of 1) pulsar glitches, 2) timing noise, 3) pulse profile stability over long time scales and 4) the extreme nulling phenomenon. The data are also included in other projects such as gamma-ray observatory support and for the Parkes Pulsar Timing Array project. In this paper we describe the current status of the project and present the first scientific results from the Parkes 12-m radio telescope. We emphasise that this project offers a straightforward means to enthuse high school students and the general public about radio astronomy while obtaining scientifically valuable data sets.

We present wide-field imaging and polarimetry at 20GHz of seven of the most extended, bright (Stot >= 0.50 Jy), high-frequency selected radio sources in the southern sky with declinations < -30 deg. Accompanying the data are brief reviews of the literature for each source, The results presented here aid in the statistical completeness of the Australia Telescope 20GHz Survey's bright source sample. The data are of crucial interest for future cosmic microwave background missions as a collection of information about candidate calibrator sources. We are able to obtain data for seven of the nine sources identified by our selection criteria. We report that Pictor A is thus far the best extragalactic calibrator candidate for the Low Frequency Instrument of the Planck European Space Agency mission due to its high level of integrated polarized flux density (0.50+/-0.06 Jy) on a scale of 10 arcmin. Six of the seven sources have a clearly detected compact radio core, with either a null or less than two percent detection of polarized emission from the nucleus. Most sources with detected jets have magnetic field alignments running in a longitudinal configuration, however PKS1333-33 exhibits transverse fields and an orthogonal change in field geometry from nucleus to jets.

The first direct detection of gravitational waves may be made through observations of pulsars. The principal aim of pulsar timing array projects being carried out worldwide is to detect ultra-low frequency gravitational waves (f ~ 10^-9 to 10^-8 Hz). Such waves are expected to be caused by coalescing supermassive binary black holes in the cores of merged galaxies. It is also possible that a detectable signal could have been produced in the inflationary era or by cosmic strings. In this paper we review the current status of the Parkes Pulsar Timing Array project (the only such project in the Southern hemisphere) and compare the pulsar timing technique with other forms of gravitational-wave detection such as ground- and space-based interferometer systems.