Enabling pulsar and fast transient searches using coherent dedispersion

We present an implementation of the coherent dedispersion algorithm capable of dedispersing high-time-resolution radio observations to many different dispersion measures (DMs). This approach allows the removal of the dispersive effects of the interstellar medium and enables searches for pulsed emission from pulsars and other millisecond-duration transients at low observing frequencies and/or high DMs where time broadening of the signal due to dispersive smearing would otherwise severely reduce the sensitivity. The implementation, called 'cdmt', for Coherent Dispersion Measure Trials, exploits the parallel processing capability of general-purpose graphics processing units to accelerate the computations. We describe the coherent dedispersion algorithm and detail how cdmt implements the algorithm to efficiently compute many coherent DM trials. We present the concept of a semi-coherent dedispersion search, where coherently dedispersed trials at coarsely separated DMs are subsequently incoherently dedispersed at finer steps in DM. The software is used in an ongoing LOFAR pilot survey to test the feasibility of performing semi-coherent dedispersion searches for millisecond pulsars at 135MHz. This pilot survey has led to the discovery of a radio millisecond pulsar -- the first at these low frequencies. This is the first time that such a broad and comprehensive search in DM-space has been done using coherent dedispersion, and we argue that future low-frequency pulsar searches using this approach are both scientifically compelling and feasible. Finally, we compare the performance of cdmt with other available alternatives.Comment: 8 pages, 7 figures, submitted to Astronomy and Computin

LOFAR discovery of the fastest-spinning millisecond pulsar in the Galactic field

We report the discovery of PSR J0952$-$0607, a 707-Hz binary millisecond pulsar which is now the fastest-spinning neutron star known in the Galactic field (i.e., outside of a globular cluster). PSR J0952$-$0607 was found using LOFAR at a central observing frequency of 135 MHz, well below the 300 MHz to 3 GHz frequencies typically used in pulsar searches. The discovery is part of an ongoing LOFAR survey targeting unassociated Fermi Large Area Telescope $\gamma$-ray sources. PSR J0952$-$0607 is in a 6.42-hr orbit around a very low-mass companion ($M_\mathrm{c}\gtrsim0.02$ M$_\odot$) and we identify a strongly variable optical source, modulated at the orbital period of the pulsar, as the binary companion. The light curve of the companion varies by 1.6 mag from $r^\prime=22.2$ at maximum to $r^\prime>23.8$, indicating that it is irradiated by the pulsar wind. Swift observations place a 3-$\sigma$ upper limit on the $0.3-10$ keV X-ray luminosity of $L_X < 1.1 \times 10^{31}$ erg s$^{-1}$ (using the 0.97 kpc distance inferred from the dispersion measure). Though no eclipses of the radio pulsar are observed, the properties of the system classify it as a black widow binary. The radio pulsed spectrum of PSR J0952$-$0607, as determined through flux density measurements at 150 and 350 MHz, is extremely steep with $\alpha\sim-3$ (where $S \propto \nu^{\alpha}$). We discuss the growing evidence that the fastest-spinning radio pulsars have exceptionally steep radio spectra, as well as the prospects for finding more sources like PSR J0952$-$0607.Comment: 9 pages, 3 figures, 1 table, published in ApJ letter

A Search for Sub-Millisecond Pulsars

We have conducted a search of 19 southern Galactic globular clusters for sub-millisecond pulsars at 660 MHz with the Parkes 64-m radio telescope. To minimize dispersion smearing we used the CPSR baseband recorder, which samples the 20 MHz observing band at the Nyquist rate. By possessing a complete description of the signal we could synthesize an optimal filterbank in software, and in the case of globular clusters of known dispersion measure, much of the dispersion could be removed using coherent techniques. This allowed for very high time resolution (25.6 us in most cases), making our searches in general sensitive to sub-millisecond pulsars with flux densities greater than about 3 mJy at 50 cm. No new pulsars were discovered, placing important constraints on the proportion of pulsars with very short spin periods in these clusters.Comment: 8 pages, 3 figures, to appear in Ap

Searching for Millisecond Pulsars: Surveys, Techniques and Prospects

Searches for millisecond pulsars (which we here loosely define as those with periods $<$ 20 ms) in the Galactic field have undergone a renaissance in the past five years. New or recently refurbished radio telescopes utilizing cooled receivers and state-of-the art digital data acquisition systems are carrying out surveys of the entire sky at a variety of radio frequencies. Targeted searches for millisecond pulsars in point sources identified by the {\it Fermi} Gamma-ray Space Telescope have proved phenomenally successful, with over 50 discoveries in the past five years. The current sample of millisecond pulsars now numbers almost 200 and, for the first time in 25 years, now outnumbers their counterparts in Galactic globular clusters. While many of these searches are motivated to find pulsars which form part of pulsar timing arrays, a wide variety of interesting systems are now being found. Following a brief overview of the millisecond pulsar phenomenon, we describe these searches and present some of the highlights of the new discoveries in the past decade. We conclude with predictions and prospects for ongoing and future surveys.Comment: 16 pages, 3 figures, accepted for publication in Classical and Quantum gravit

Realfast: Real-Time, Commensal Fast Transient Surveys with the Very Large Array

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.Comment: Accepted to ApJS Special Issue on Data; 11 pages, 4 figure