The detection of non-repeating or irregular events in time-domain radio astronomy has gained importance over the last decade due to the discovery of fast radio bursts. Existing or upcoming radio telescopes are gathering more and more data and consequently the software, which is an important part of these telescopes, must process large data volumes at high data rates. Data has to be searched through to detect new and interesting events, often in real-time. These requirements necessitate new and fast algorithms which must process data quickly and accurately. In this work we present new algorithms for single pulse detection using boxcar filters. We have quantified the signal loss introduced by single pulse detection algorithms which use boxcar filters and based on these results, we have designed two distinct "lossy" algorithms. Our lossy algorithms use an incomplete set of boxcar filters to accelerate detection at the expense of a small reduction in detected signal power. We present formulae for signal loss, descriptions of our algorithms and their parallel implementation on NVIDIA GPUs using CUDA. We also present tests of correctness, tests on artificial data and the performance achieved. Our implementation can process SKA-MID-like data 266$\times$ faster than real-time on a NVIDIA P100 GPU and 500x faster than real-time on a NVIDIA Titan V GPU with a mean signal power loss of 7%. We conclude with prospects for single pulse detection for beyond SKA era, nanosecond time resolution radio astronomy.

The study of binary pulsars enables tests of general relativity. Orbital motion in binary systems causes the apparent pulsar spin frequency to drift, reducing the sensitivity of periodicity searches. Acceleration searches are methods that account for the effect of orbital acceleration. Existing methods are currently computationally expensive, and the vast amount of data that will be produced by next generation instruments such as the Square Kilometre Array (SKA) necessitates real-time acceleration searches, which in turn requires the use of High Performance Computing (HPC) platforms. We present our implementation of the Correlation Technique for the Fourier Domain Acceleration Search (FDAS) algorithm on Graphics Processor Units (GPUs). The correlation technique is applied as a convolution with multiple Finite Impulse Response filters in the Fourier domain. Two approaches are compared: the first uses the NVIDIA cuFFT library for applying Fast Fourier Transforms (FFTs) on the GPU, and the second contains a custom FFT implementation in GPU shared memory. We find that the FFT shared memory implementation performs between 1.5 and 3.2 times faster than our cuFFT-based application for smaller but sufficient filter sizes. It is also 4 to 6 times faster than the existing GPU and OpenMP implementations of FDAS. This work is part of the AstroAccelerate project, a many-core accelerated time-domain signal processing library for radio astronomy.

The innermost stable cicular orbit (ISCO) of an accretion disc orbiting a neutron star (NS) is often assumed a unique prediction of general relativity. However, it has been argued that ISCO also appears around highly elliptic bodies described by Newtonian theory. In this sense, the behaviour of an ISCO around a rotating oblate neutron star is formed by the interplay between relativistic and Newtonian effects. Here we briefly explore the consequences of this interplay using a straightforward analytic approach as well as numerical models that involve modern NS equations of state. We examine the ratio K between the ISCO radius and the radius of the neutron star. We find that, with growing NS spin, the ratio K first decreases, but then starts to increase. This non-monotonic behaviour of K can give rise to a neutron star spin interval in which ISCO appears for two very different ranges of NS mass. This may strongly affect the distribution of neutron stars that have an ISCO (ISCO-NS). When (all) neutron stars are distributed around a high mass M0, the ISCO-NS spin distribution is roughly the same as the spin distribution corresponding to all neutron stars. In contrast, if M0 is low, the ISCO-NS distribution can only have a peak around a high value of spin. Finally, an intermediate value of M0 can imply an ISCO-NS distribution divided into two distinct groups of slow and fast rotators. Our findings have immediate astrophysical applications. They can be used for example to distinguish between different models of high-frequency quasiperiodic oscillations observed in low-mass NS X-ray binaries.