Search

The birth mass function of neutron stars encodes rich information about supernova explosions, double star evolution, and properties of matter under extreme conditions. To date, it has remained poorly constrained by observations, however. Applying probabilistic corrections to account for mass accreted by recycled pulsars in binary systems to mass measurements of 90 neutron stars, we find that the birth masses of neutron stars can be described by a unimodal distribution that smoothly turns on at $\mathbf{\unit[1.1]{\mathrm{M}_{\odot}}}$, peaks at $\mathbf{\approx \unit[1.27]{\mathrm{M}_{\odot}}}$, before declining as a steep power law. Such a ``turn-on" power-law distribution is strongly favoured against the widely-adopted empirical double-Gaussian model at the $\mathbf{3\sigma}$ level. The power-law shape may be inherited from the initial mass function of massive stars, but the relative dearth of massive neutron stars implies that single stars with initial masses greater than $\mathbf{\approx \unit[18]{\mathrm{M}_{\odot}}}$ do not form neutron stars, in agreement with the absence of massive red supergiant progenitors to supernovae., Comment: 57 pages, 14 figures, 4 tables

We conduct the first-ever Pulsar Polarization Array (PPA) analysis to detect the ultralight Axion-Like Dark Matter (ALDM) using the polarization data of 22 millisecond pulsars from the third data release of Parkes Pulsar Timing Array. As one of the major dark matter candidates, the ultralight ALDM exhibits a pronounced wave nature on astronomical scales and offers a promising solution to small-scale structure issues within local galaxies. While the linearly polarized pulsar light travels through the ALDM galactic halo, its position angle (PA) can be subject to an oscillation induced by the ALDM Chern-Simons coupling with electromagnetic field. The PPA is thus especially suited for detecting the ultralight ALDM by correlating polarization data across the arrayed pulsars. To accomplish this task, we develop an advanced Bayesian analysis framework that allows us to construct pulsar PA residual time series, model noise contributions properly and search for pulsar cross-correlations. We find that for an ALDM density of $\rho_0=0.4\,\textrm{GeV}/\textrm{cm}^3$, the Parkes PPA offers the best global limits on the ALDM Chern-Simons coupling, namely $\lesssim 10^{-13.5}-10^{-12.2}~{\rm GeV}^{-1}$, for the mass range of $10^{-22} - 10^{-21}~{\rm eV}$. The crucial role of pulsar cross-correlation in recognizing the nature of the derived limits is also highlighted., Comment: 6+15 pages, 10 figures, 2 tables, submitted to the journal

This paper demonstrates the impact of state-of-the-art instrumental calibration techniques on the precision of arrival times obtained from 9.6 years of observations of millisecond pulsars using the Murriyang 64-m CSIRO Parkes Radio Telescope. Our study focuses on 21-cm observations of 25 high-priority pulsars that are regularly observed as part of the Parkes Pulsar Timing Array (PPTA) project, including those predicted to be the most susceptible to calibration errors. We employ Measurement Equation Template Matching (METM) for instrumental calibration and Matrix Template Matching (MTM) for arrival time estimation, resulting in significantly improved timing residuals with up to a sixfold reduction in white noise compared to arrival times estimated using Scalar Template Matching and conventional calibration based on the Ideal Feed Assumption. The median relative reduction in white noise is 33 percent, and the maximum absolute reduction is 4.5 microseconds. For PSR J0437-4715, METM and MTM reduce the best-fit power-law amplitude (2.7 sigma) and spectral index (1.7 sigma) of the red noise in the arrival time residuals, which can can be tentatively interpreted as mitigation of 1/f noise due to otherwise unmodeled steps in polarimetric response. These findings demonstrate the potential to directly enhance the sensitivity of pulsar timing array experiments through more accurate methods of instrumental calibration and arrival time estimation., Comment: 30 pages, 9 figures, accepted for publication in The Astrophysical Journal

The observation of neutron stars enables the otherwise impossible study of fundamental physical processes. The timing of binary radio pulsars is particularly powerful, as it enables precise characterization of their (three-dimensional) positions and orbits. PSR~J0437$-$4715 is an important millisecond pulsar for timing array experiments and is also a primary target for the Neutron Star Interior Composition Explorer (NICER). The main aim of the NICER mission is to constrain the neutron star equation of state by inferring the compactness ($M_p/R$) of the star. Direct measurements of the mass $M_p$ from pulsar timing therefore substantially improve constraints on the radius $R$ and the equation of state. Here we use observations spanning 26 years from Murriyang, the 64-m Parkes radio telescope, to improve the timing model for this pulsar. Among the new precise measurements are the pulsar mass $M_p=1.418\pm 0.044$ $M_{\odot}$, distance $D=156.96 \pm 0.11$ pc, and orbital inclination angle $i=137.506 \pm 0.016^\circ$, which can be used to inform the X-ray pulse profile models inferred from NICER observations. We demonstrate that these results are consistent between multiple data sets from the Parkes Pulsar Timing Array (PPTA), each modeled with different noise assumptions. Using the longest available PPTA data set, we measure an apparent second derivative of the pulsar spin frequency and discuss how this can be explained either by kinematic effects due to the proper motion and radial velocity of the pulsar or excess low-frequency noise such as a gravitational-wave background., Comment: 13 pages, 3 figures, accepted for publication in Astrophysical Journal Letters

Radio emission from magnetars provides a unique probe of the relativistic, magnetized plasma within the near-field environment of these ultra-magnetic neutron stars. The transmitted waves can undergo birefringent and dispersive propagation effects that result in frequency-dependent conversions of linear to circularly polarized radiation and vice-versa, thus necessitating classification when relating the measured polarization to the intrinsic properties of neutron star and fast radio burst (FRB) emission sites. We report the detection of such behavior in 0.7-4 GHz observations of the P = 5.54 s radio magnetar XTE J1810$-$197 following its 2018 outburst. The phenomenon is restricted to a narrow range of pulse phase centered around the magnetic meridian. Its temporal evolution is closely coupled to large-scale variations in magnetic topology that originate from either plastic motion of an active region on the magnetar surface or free precession of the neutron star crust. Our model of the effect deviates from simple theoretical expectations for radio waves propagating through a magnetized plasma. Birefringent self-coupling between the transmitted wave modes, line-of-sight variations in the magnetic field direction and differences in particle charge or energy distributions above the magnetic pole are explored as possible explanations. We discuss potential links between the immediate magneto-ionic environments of magnetars and those of FRB progenitors., Comment: 22 pages (inc. extended data and supplementary materials), 12 figures, 3 tables. Published in Nature Astronomy

The noise in millisecond pulsar (MSP) timing data can include contributions from observing instruments, the interstellar medium, the solar wind, solar system ephemeris errors, and the pulsars themselves. The noise environment must be accurately characterized in order to form the null hypothesis from which signal models can be compared, including the signature induced by nanohertz-frequency gravitational waves (GWs). Here we describe the noise models developed for each of the MSPs in the Parkes Pulsar Timing Array (PPTA) third data release, which have been used as the basis of a search for the isotropic stochastic GW background. We model pulsar spin noise, dispersion measure variations, scattering variations, events in the pulsar magnetospheres, solar wind variability, and instrumental effects. We also search for new timing model parameters and detected Shapiro delays in PSR~J0614$-$3329 and PSR~J1902$-$5105. The noise and timing models are validated by testing the normalized and whitened timing residuals for Gaussianity and residual correlations with time. We demonstrate that the choice of noise models significantly affects the inferred properties of a common-spectrum process. Using our detailed models, the recovered common-spectrum noise in the PPTA is consistent with a power law with a spectral index of $\gamma=13/3$, the value predicted for a stochastic GW background from a population of supermassive black hole binaries driven solely by GW emission., Comment: 18 pages, 10 figures. Accepted for publication in ApJL

We present the third data release from the Parkes Pulsar Timing Array (PPTA) project. The release contains observations of 32 pulsars obtained using the 64-m Parkes "Murriyang" radio telescope. The data span is up to 18 years with a typical cadence of 3 weeks. This data release is formed by combining an updated version of our second data release with $\sim 3$ years of more recent data primarily obtained using an ultra-wide-bandwidth receiver system that operates between 704 and 4032 MHz. We provide calibrated pulse profiles, flux-density dynamic spectra, pulse times of arrival, and initial pulsar timing models. We describe methods for processing such wide-bandwidth observations, and compare this data release with our previous release., Comment: 15 pages, 6 figures. Accepted for publication in PASA

Pulsar timing arrays aim to detect nanohertz-frequency gravitational waves (GWs). A background of GWs modulates pulsar arrival times and manifests as a stochastic process, common to all pulsars, with a signature spatial correlation. Here we describe a search for an isotropic stochastic gravitational-wave background (GWB) using observations of 30 millisecond pulsars from the third data release of the Parkes Pulsar Timing Array (PPTA), which spans 18 years. Using current Bayesian inference techniques we recover and characterize a common-spectrum noise process. Represented as a strain spectrum $h_c = A(f/1 {\rm yr}^{-1})^{\alpha}$, we measure $A=3.1^{+1.3}_{-0.9} \times 10^{-15}$ and $\alpha=-0.45 \pm 0.20$ respectively (median and 68% credible interval). For a spectral index of $\alpha=-2/3$, corresponding to an isotropic background of GWs radiated by inspiraling supermassive black hole binaries, we recover an amplitude of $A=2.04^{+0.25}_{-0.22} \times 10^{-15}$. However, we demonstrate that the apparent signal strength is time-dependent, as the first half of our data set can be used to place an upper limit on $A$ that is in tension with the inferred common-spectrum amplitude using the complete data set. We search for spatial correlations in the observations by hierarchically analyzing individual pulsar pairs, which also allows for significance validation through randomizing pulsar positions on the sky. For a process with $\alpha=-2/3$, we measure spatial correlations consistent with a GWB, with an estimated false-alarm probability of $p \lesssim 0.02$ (approx. $2\sigma$). The long timing baselines of the PPTA and the access to southern pulsars will continue to play an important role in the International Pulsar Timing Array., Comment: 19 pages, 10 figures, Accepted for publication in ApJL

In 2018 an ultra-wide-bandwidth low-frequency (UWL) receiver was installed on the 64-m Parkes Radio Telescope enabling observations with an instantaneous frequency coverage from 704 to 4032 MHz. Here, we present the analysis of a three-year data set of 35 millisecond pulsars observed with the UWL by the Parkes Pulsar Timing Array (PPTA), using wideband timing methods. The two key differences compared to typical narrow-band methods are, firstly, generation of two-dimensional templates accounting for pulse shape evolution with frequency and, secondly, simultaneous measurements of the pulse time-of-arrival (ToA) and dispersion measure (DM). This is the first time that wideband timing has been applied to a uniform data set collected with a single large-fractional bandwidth receiver, for which such techniques were originally developed. As a result of our study, we present a set of profile evolution models and new timing solutions including initial noise analysis. Precision of our ToA and DM measurements is in the range of 0.005 $-$ 2.08 $\mu$s and (0.043$-$14.24)$\times10^{-4}$ cm$^{-3}$ pc, respectively, with 94% of the pulsars achieving a median ToA uncertainty of less than 1 $\mu$s., Comment: 29 pages, 12 figures, Accepted for publication in ApJ

Composed of ultralight bosons, fuzzy dark matter provides an intriguing solution to challenges that the standard cold dark matter model encounters on sub-galactic scales. The ultralight dark matter with mass $m\sim10^{-23} \rm{eV}$ will induce a periodic oscillation in gravitational potentials with a frequency in the nanohertz band, leading to observable effects in the arrival times of radio pulses from pulsars. Unlike scalar dark matter, pulsar timing signals induced by the vector dark matter are dependent on the oscillation direction of the vector fields. In this work, we search for ultralight vector dark matter in the mass range of $[2\times 10^{-24}, 2\times 10^{-22}]{\rm{eV}}$ through its gravitational effect in the Parkes Pulsar Timing Array (PPTA) second data release. Since no statistically significant detection is made, we place $95\%$ upper limits on the local dark matter density as $\rho_{\rm{\tiny{VF}}} \lesssim 5{\rm{GeV/cm^{3}}}$ for $m\lesssim 10^{-23}{\rm{eV}}$. As no preferred direction is found for the vector dark matter, these constraints are comparable to those given by the scalar dark matter search with an earlier 12-year data set of PPTA., Comment: 9 pages, 3 figures, 2 tables; accepted for publication as a Letter in Phys. Rev. D

New classes of astronomical objects are often discovered serendipitously. The enormous data volumes produced by recent high-time resolution, radio-telescope surveys imply that efficient algorithms are required for a discovery. Such algorithms are usually tuned to detect specific, known sources. Existing data sets therefore likely contain unknown astronomical sources, which will remain undetected unless algorithms are developed that can detect a more diverse range of signals. We present the Single-dish PARKES data challenge for finding the uneXpected (SPARKESX), a compilation of real and simulated high-time resolution observations. SPARKESX comprises three mock surveys from the Parkes "Murriyang" radio telescope. A broad selection of simulated and injected expected signals (such as pulsars, fast radio bursts), poorly characterised signals (plausible flare star signatures) and unknown unknowns are generated for each survey. The goal of this challenge is to aid in the development of new algorithms that can detect a wide-range of source types. We show how successful a typical pipeline based on the standard pulsar search software, PRESTO, is at finding the injected signals. The dataset is publicly available at https://doi.org/10.25919/fd4f-0g20., Comment: Accepted for publication in MNRAS. 18 pages, 8 figures, 4 table

Pulsar timing array collaborations have recently reported evidence for a noise process with a common spectrum among the millisecond pulsars in the arrays. The spectral properties of this common-noise process are consistent with expectations for an isotropic gravitational-wave background (GWB) from inspiralling supermassive black-hole binaries. However, recent simulation analyses based on Parkes Pulsar Timing Array data indicate that such a detection may arise spuriously. In this paper, we use simulated pulsar timing array datasets to further test the robustness of the inference methods for spectral and spatial correlations from a GWB. Expanding on our previous results, we find strong support (Bayes factors exceeding $10^5$) for the presence of a common-spectrum noise process in datasets where no common process is present, under a wide range of timing noise prescriptions per pulsar. We show that these results are highly sensitive to the choice of Bayesian priors on timing noise parameters, with priors that more closely match the injected distributions of timing noise parameters resulting in diminished support for a common-spectrum noise process. These results emphasize shortcomings in current methods for inferring the presence of a common-spectrum process, and imply that the detection of a common process is not a reliable precursor to detection of the GWB. Future searches for the nanohertz GWB should remain focussed on detecting spatial correlations, and make use of more tailored specifications for a common-spectrum noise process., Comment: 12 pages, 12 figures. Accepted for publication in MNRAS

We report the discovery of a new 5.78 ms-period millisecond pulsar (MSP), PSR J1740-5340B (NGC 6397B), in an eclipsing binary system discovered with the Parkes radio telescope (now also known as Murriyang), Australia, and confirmed with the MeerKAT radio telescope in South Africa. The measured orbital period, 1.97 days, is the longest among all eclipsing binaries in globular clusters (GCs) and consistent with that of the coincident X-ray source U18, previously suggested to be a 'hidden MSP'. Our XMM-Newton observations during NGC 6397B's radio quiescent epochs detected no X-ray flares. NGC 6397B is either a transitional MSP or an eclipsing binary in its initial stage of mass transfer after the companion star left the main sequence. The discovery of NGC 6397B potentially reveals a subgroup of extremely faint and heavily obscured binary pulsars, thus providing a plausible explanation to the apparent dearth of binary neutron stars in core-collapsed GCs as well as a critical constraint on the evolution of GCs.

Pulsar timing array experiments have recently reported strong evidence for a common-spectrum stochastic process with a strain spectral index consistent with that expected of a nanohertz-frequency gravitational-wave background, but with negligible yet non-zero evidence for spatial correlations required for a definitive detection. However, it was pointed out by the Parkes Pulsar Timing Array (PPTA) collaboration that the same models used in recent analyses resulted in strong evidence for a common-spectrum process in simulations where none is present. In this work, we introduce a methodology to distinguish pulsar power spectra with the same amplitude from noise power spectra of similar but distinct amplitudes. The former is the signature of a spatially uncorrelated pulsar term of a nanohertz gravitational-wave background, whereas the latter could represent ensemble pulsar noise properties. We test the methodology on simulated data sets. We find that the reported common process in PPTA pulsars is indeed consistent with the spectral feature of a pulsar term. We recommend this methodology as one of the validity tests that the real astrophysical and cosmological backgrounds should pass, as well as for inferences about the spatially uncorrelated component of the background., Comment: 9 pages, 2 figures

We describe a new software package for simulating channelised, high-time resolution data streams from radio telescopes. The software simulates data from the telescope and observing system taking into account the observation strategy, receiver system and digitisation. The signatures of pulsars, fast radio bursts and flare stars are modelled, including frequency-dependent effects such as scattering and scintillation. We also simulate more generic signals using spline curves and images. Models of radio frequency interference include signals from satellites, terrestrial transmitters and impulsive, broadband signals. The simulated signals can also be injected into real data sets. Uses of this software include the production of machine learning training data sets, development and testing of new algorithms to search for anomalous patterns and to characterise processing pipelines., Comment: 12 pages, 11 figures, MNRAS in press

Fast radio bursts (FRBs) are brief, intense flashes of radio waves from unidentified extragalactic sources. Polarized FRBs originate in highly magnetized environments. We report observations of the repeating FRB 20190520B spanning seventeen months , which show its amount of Faraday rotation is highly variable and twice changes its sign. The FRB also depolarizes below radio frequencies around 1 to 3 GHz. We interpret these properties as due to change in the parallel component of the integrated magnetic field along the line-of-sight, including reversals. This could result from propagation through a turbulent, magnetized screen of plasma located between $10^{-5}$ to 100 parsecs of the FRB source. This is consistent with the bursts passing through the stellar wind of a binary companion of the FRB source., Comment: arXiv:2203.08151 has been merged into this

We demonstrate how radio pulsars can be used as random number generators. Specifically, we focus on publicly verifiable randomness (PVR), in which the same sequence of trusted and verifiable random numbers is obtained by multiple parties. PVR is a critical building block for many processes and algorithms (including cryptography, scientific trials, electoral audits and international treaties). However, current approaches (based on number theory) may soon become vulnerable to quantum computers, motivating a growing demand for PVR based on natural physical phenomena. In this context, we explore pulsars as a potential physical PVR source. We first show that bit sequences extracted from the measured flux densities of a bright millisecond pulsar can pass standardised tests for randomness. We then quantify three illustrative methods of bit-extraction from pulsar flux density sequences, using simultaneous observations of a second pulsar carried out with the Parkes telescope in Australia and the Five-hundred-metre Aperture Spherical radio Telescope (FAST) in China, supported by numerical simulations. We demonstrate that the same bit sequence can indeed be obtained at both observatories, but the ubiquitous presence of radiometer noise needs to be accounted for when determining the expected bit error rate between two independent sequences. We discuss our results in the context of an imaginary use-case in which two mutually distrusting parties wish to obtain the same random bit sequence, exploring potential methods to mitigate against a malicious participant., Comment: 12 pages, 4 figures, accepted to Astronomy & Computing

The nature of dark matter remains obscure in spite of decades of experimental efforts. The mass of dark matter candidates can span a wide range, and its coupling with the Standard Model sector remains uncertain. All these unknowns make the etection of dark matter extremely challenging. Ultralight dark matter, with $m \sim10^{-22}$ eV, is proposed to reconcile the disagreements between observations and predictions from simulations of small-scale structures in the cold dark matter paradigm, while remaining consistent with other observations. Because of its large de Broglie wavelength and large local occupation number within galaxies, ultralight dark matter behaves like a coherently oscillating background field with an oscillating frequency dependent on its mass. If the dark matter particle is a spin-1 dark photon, such as the $U(1)_B$ or $U(1)_{B-L}$ gauge boson, it can induce an external oscillating force and lead to displacements of test masses. Such an effect would be observable in the form of periodic variations in the arrival times of radio pulses from highly stable millisecond pulsars. In this study, we search for evidence of ultralight dark photon dark matter (DPDM) using 14-year high-precision observations of 26 pulsars collected with the Parkes Pulsar Timing Array. While no statistically significant signal is found, we place constraints on coupling constants for the $U(1)_B$ and $U(1)_{B-L}$ DPDM. Compared with other experiments, the limits on the dimensionless coupling constant $\epsilon$ achieved in our study are improved by up to two orders of magnitude when the dark photon mass is smaller than $3\times10^{-22}$~eV ($10^{-22}$~eV) for the $U(1)_{B}$ ($U(1)_{B-L}$) scenario., Comment: 12 pages, 4 figures

A cosmological first-order phase transition is expected to produce a stochastic gravitational wave background. If the phase transition temperature is on the MeV scale, the power spectrum of the induced stochastic gravitational waves peaks around nanohertz frequencies, and can thus be probed with high-precision pulsar timing observations. We search for such a stochastic gravitational wave background with the latest data set of the Parkes Pulsar Timing Array. We find no evidence for a Hellings-Downs spatial correlation as expected for a stochastic gravitational wave background. Therefore, we present constraints on first-order phase transition model parameters. Our analysis shows that pulsar timing is particularly sensitive to the low-temperature ($T \sim 1 - 100$ MeV) phase transition with a duration $(\beta/H_*)^{-1}\sim 10^{-2}-10^{-1}$ and therefore can be used to constrain the dark and QCD phase transitions., Comment: 9 pages, 4 figures, 2 tables

Optical depth variations in the Galactic neutral interstellar medium (ISM) with spatial scales from hundreds to thousands of astronomical units have been observed through HI absorption against pulsars and continuum sources, while extremely small structures with spatial scales of tens of astronomical units remain largely unexplored. The nature and formation of such tiny-scale atomic structures (TSAS) need to be better understood. We report a tentative detection of TSAS with a signal-to-noise ratio of 3.2 toward PSR B1557$-$50 in the second epoch of two Parkes sessions just 0.36 yr apart, which are the closest-spaced spectral observations toward this pulsar. One absorption component showing marginal variations has been identified. Based on the pulsar's proper motion of 14 mas $\rm yr^{-1}$ and the component's kinematic distance of 3.3 kpc, the corresponding characteristic spatial scale is 17 au, which is among the smallest sizes of known TSAS. Assuming a similar line-of-sight (LOS) depth, the tentative TSAS cloud detected here is overdense and overpressured relative to the cold neutral medium (CNM), and can radiatively cool fast enough to be in thermal equilibrium with the ambient environment. We find that turbulence is not sufficient to confine the overpressured TSAS. We explore the LOS elongation that would be required for the tentative TSAS to be at the canonical CNM pressure, and find that it is $\sim5000$ -- much larger than filaments observed in the ISM. We see some evidence of line width and temperature variations in the CNM components observed at the two epochs, as predicted by models of TSAS-like cloud formation colliding warm neutral medium flows., Comment: 11 pages, 2 figures, 2 tables, accepted for publication in the ApJL

Using the Five-hundred-meter Aperture Spherical radio Telescope (FAST), we have recorded 10^5 single pulses from PSR J1022+1001. We studied the polarization properties, their energy distribution and their times of arrival. This is only possible with the high sensitivity available using FAST. There is no indication that PSR~J1022+1001 exhibits giant pulse, nulling or traditional mode changing phenomena. The energy in the leading and trailing components of the integrated profile is shown to be correlated. The degree of both linear and circular polarization increases with the pulse flux density for individual pulses. Our data indicates that pulse jitter leads to an excess noise in the timing residuals of 67 ns when scaled to one hour, which is consistent with Liu et al. (2015). We have unsuccessfully trialled various methods to improve timing precision through the selection of specific single pulses. Our work demonstrates that FAST can detect individual pulses from pulsars that are observed in order to detect and study gravitational waves. This capability enables detailed studies, and parameterisation, of the noise processes that affect the sensitivity of a pulsar timing array., Comment: Accepted by ApJ

We report the first wideband monitoring observations of PSR J1803-3002A, a relatively bright millisecond pulsar in the globular cluster NGC 6522 with a spin period of 7.1 ms and no known binary companion. These observations were performed using the Parkes 64-m radio telescope with the Ultra-Wideband Low (UWL) receiver system, which covers 704 to 4032 MHz. We confirm that PSR J1803-3002A is an isolated millisecond pulsar located near the cluster center and probe the emission properties of the pulsar over the wide observed band. The mean pulse profile consists of three components, with the outer components becoming more prominent at higher frequencies, and a mean spectral index for the pulsed emission of -1.66+/-0.07 over the observed band. The fractional linear and circular polarization increase with increasing frequency, which is unusual for pulsars. We determine a Faraday rotation measure of -107+/-6 rad m-2 for the pulsar. PSR J1803-3002A is a distant pulsar in the Galactic plane, but there is no evidence of pulse broadening due to interstellar scattering in our observations. These results demonstrate the power of ultra-wideband receiving and signal processing systems.

Fast spinning (e.g., sub-second) neutron star with ultra-strong magnetic fields (or so-called magnetar) is one of the promising origins of repeating fast radio bursts (FRBs). Here we discuss circularly polarised emissions produced by propagation effects in the magnetosphere of fast spinning magnetars. We argue that the polarisation-limiting region is well beyond the light cylinder, suggesting that wave mode coupling effects are unlikely to produce strong circular polarisation for fast spinning magnetars. Cyclotron absorption could be significant if the secondary plasma density is high. However, high degrees of circular polarisation can only be produced with large asymmetries in electrons and positrons. We draw attention to the non-detection of circular polarisation in current observations of known repeating FRBs. We suggest that the circular polarisation of FRBs could provide key information on their origins and help distinguish different radiation mechanisms., Comment: ApJ accepted

Intensity scintillations of radio pulsars are known to originate from interference between waves scattered by the electron density irregularities of interstellar plasma, often leading to parabolic arcs in the two-dimensional power spectrum of the recorded dynamic spectrum. The degree of arc curvature depends on the distance to the scattering plasma and its transverse velocity with respect to the line-of-sight. We report the observation of annual and orbital variations in the curvature of scintillation arcs over a period of 16 years for the bright millisecond pulsar, PSR J0437-4715. These variations are the signature of the relative transverse motions of the Earth, pulsar, and scattering medium, which we model to obtain precise measurements of parameters of the pulsar's binary orbit and the scattering medium itself. We observe two clear scintillation arcs in most of our $>$5000 observations and we show that they originate from scattering by thin screens located at distances $D_1 = 89.8 \pm 0.4$ pc and $D_2 = 124 \pm 3$ pc from Earth. The best-fit scattering model we derive for the brightest arc yields the pulsar's orbital inclination angle $i = 137.1 \pm 0.3^\circ$, and longitude of ascending node, $\Omega=206.3\pm0.4^\circ$. Using scintillation arcs for precise astrometry and orbital dynamics can be superior to modelling variations in the diffractive scintillation timescale, because the arc curvature is independent of variations in the level of turbulence of interstellar plasma. This technique can be used in combination with pulsar timing to determine the full three-dimensional orbital geometries of binary pulsars, and provides parameters essential for testing theories of gravity and constraining neutron star masses., Comment: 19 pages, 10 figures. Accepted for publication in ApJ

We report the discovery of a highly dispersed fast radio burst, FRB~181123, from an analysis of $\sim$1500~hr of drift-scan survey data taken using the Five-hundred-meter Aperture Spherical radio Telescope (FAST). The pulse has three distinct emission components, which vary with frequency across our 1.0--1.5~GHz observing band. We measure the peak flux density to be $>0.065$~Jy and the corresponding fluence $>0.2$~Jy~ms. Based on the observed dispersion measure of 1812~cm$^{-3}$~pc, we infer a redshift of $\sim 1.9$. From this, we estimate the peak luminosity and isotropic energy to be $\lesssim 2\times10^{43}$~erg~s$^{-1}$ and $\lesssim 2\times10^{40}$~erg, respectively. With only one FRB from the survey detected so far, our constraints on the event rate are limited. We derive a 95\% confidence lower limit for the event rate of 900 FRBs per day for FRBs with fluences $>0.025$~Jy~ms. We performed follow-up observations of the source with FAST for four hours and have not found a repeated burst. We discuss the implications of this discovery for our understanding of the physical mechanisms of FRBs., Comment: 7 pages, 3 figures, accepted for publication in ApJL

We report the discovery of an eclipsing binary millisecond pulsar in the globular cluster M92 (NGC6341) with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). PSR J1717+4308A, or M92A, has a pulse frequency of 316.5~Hz (3.16~ms) and a dispersion measure of 35.45 pc cm$^{-3}$. The pulsar is a member of a binary system with an orbital period of 0.20~days around a low-mass companion which has a median mass of $\sim$0.18~\Ms. From observations so far, at least two eclipsing events have been observed in each orbit. The longer one lasted for ~5000~s in the orbital phase range 0.1--0.5. The other lasted for ~500~s and occurred between 1000--2000~s before or after the longer eclipsing event. The lengths of these two eclipsing events also change. These properties suggest that J1717+4308A is a ``red-back'' system with a low-mass main sequence or sub-giant companion. Timing observations of the pulsar and further searches of the data for additional pulsars are ongoing., Comment: 9 pages, 2 figures, submitted to ApJL, comments are welcomed!

Individual supermassive black hole binaries in non-circular orbits are possible nanohertz gravitational wave sources for the rapidly maturing Pulsar Timing Array experiments. We develop an accurate and efficient approach to compute Pulsar Timing Array signals due to gravitational waves from inspiraling supermassive black hole binaries in relativistic eccentric orbits. Our approach employs a Keplerian-type parametric solution to model third post-Newtonian accurate precessing eccentric orbits while a novel semi-analytic prescription is provided to model the effects of quadrupolar order gravitational wave emission. These inputs lead to a semi-analytic prescription to model such signals, induced by non-spinning black hole binaries inspiralling along arbitrary eccentricity orbits. Additionally, we provide a fully analytic prescription to model Pulsar Timing Array signals from black hole binaries inspiraling along moderately eccentric orbits, influenced by Boetzelet al.[Phys. Rev. D 96,044011(2017)]. These approaches are being incorporated into Enterprise and TEMPO2 for searching the presence of such binaries in Pulsar Timing Array datasets.

The emission from PSR J1107-5907 is erratic. Sometimes the radio pulse is undetectable, at other times the pulsed emission is weak, and for short durations the emission can be very bright. In order to improve our understanding of these state changes, we have identified archival data sets from the Parkes radio telescope in which the bright emission is present, and find that the emission never switches from the bright state to the weak state, but instead always transitions to the off state. Previous work had suggested the identification of the off state as an extreme manifestation of the weak state. However, the connection between the off and bright emission reported here suggests that the emission can be interpreted as undergoing only two emission states: a bursting state consisting of both bright pulses and nulls as well as the weak-emission state., Comment: 7 pages, 6 figures, accepted by ApJ

We report the first wideband observations of pulsars C, D and J in the globular cluster 47Tucanae (NGC 104) using the Ultra-Wideband Low (UWL) receiver system recently installed on the Parkes 64 m radio telescope. The wide frequency range of the UWL receiver (704-4032 MHz), along with the well-calibrated system, allowed us to obtain flux density measurements and polarization pulse profiles. The mean pulse profiles have significant linear and circular polarization, allowing for determination of the Faraday rotation measure for each pulsar. Precise measurements of the dispersion measures show a significant deviation in the value for pulsar D compared to earlier results. Searches for new pulsars in the cluster are on-going and we have determined optimal bands for such searches using the Parkes UWL receiver system.

We identified a highly-polarized, steep-spectrum radio source in a deep image with the Australian Square Kilometre Array Pathfinder (ASKAP) telescope at 888 MHz. After considering and rejecting a stellar origin for this source, we discovered a new millisecond pulsar (MSP) using observations from the Parkes radio telescope. This pulsar has period 2.77 ms and dispersion measure 228.27 pc/cm**3. Although this pulsar does not yet appear to be particularly remarkable, the short spin period, wide profile and high dispersion measure do make it relatively hard to discover through traditional blind periodicity searches. Over the course of several weeks we see changes in the barycentric period of this pulsar that are consistent with orbital motion in a binary system, but the properties of any binary need to be confirmed by further observations. While even a deep ASKAP survey may not identify large numbers of new MSPs compared to the existing population, it would be competitive with existing all-sky surveys and could discover interesting new MSPs at high Galactic latitude without the need for computationally-expensive all-sky periodicity searches., Comment: ApJ, in press

We describe PSR J1926-0652, a pulsar recently discovered with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). Using sensitive single-pulse detections from FAST and long-term timing observations from the Parkes 64-m radio telescope, we probed phenomena on both long and short time scales. The FAST observations covered a wide frequency range from 270 to 800 MHz, enabling individual pulses to be studied in detail. The pulsar exhibits at least four profile components, short-term nulling lasting from 4 to 450 pulses, complex subpulse drifting behaviours and intermittency on scales of tens of minutes. While the average band spacing P3 is relatively constant across different bursts and components, significant variations in the separation of adjacent bands are seen, especially near the beginning and end of a burst. Band shapes and slopes are quite variable, especially for the trailing components and for the shorter bursts. We show that for each burst the last detectable pulse prior to emission ceasing has different properties compared to other pulses. These complexities pose challenges for the classic carousel-type models., Comment: 13pages with 12 figures

To assist with the commissioning (Jiang et al. 2019) of the Five-hundred-meter Aperture Spherical radio Telescope (FAST), we performed a pulsar search, with the primary goal of developing and testing the pulsar data acquisition and processing pipelines. We tested and used three pipelines, two (P1 and P2 hereafter) searched for the periodic signature of pulsars whereas the other one was used to search for bright single pulses (P3 hereafter). A pulsar candidate was discovered in the observation on the 22nd August, 2017, and later confirmed by the Parkes radio telescope on the 10th September, 2017. This pulsar, named PSR J1900-0134, was the first pulsar discovered by FAST. The pulsar has a pulse period of 1.8 s and a dispersion measure (DM) of 188\,pc\,cm$^{-3}$., Comment: Sci. China-Phys. Mech. Astron. 62, 959508 (2019)

Flux density measurements provide fundamental observational parameters that describe a pulsar. In the current pulsar catalogue, 27% of radio pulsars have no flux density measurement in the 20 cm observing band. Here, we present the first measurements of the flux densities in this band for 32 pulsars observed using the Parkes radio telescope and provide updated pulse profiles for these pulsars. We have used both archival and new observations to make these measurements. Various schemes exist for measuring flux densities. We show how the flux densities measured vary between these methods and how the presence of radio-frequency-interference will bias flux density measurements, Comment: Accepted by RAA

It is widely accepted that dark matter contributes about a quarter of the critical mass-energy density in our Universe. The nature of dark matter is currently unknown, with the mass of possible constituents spanning nearly one hundred orders of magnitude. The ultralight scalar field dark matter, consisting of extremely light bosons with $m \sim 10^{-22}$ eV and often called "fuzzy" dark matter, provides intriguing solutions to some challenges at sub-Galactic scales for the standard cold dark matter model. As shown by Khmelnitsky and Rubakov, such a scalar field in the Galaxy would produce an oscillating gravitational potential with nanohertz frequencies, resulting in periodic variations in the times of arrival of radio pulses from pulsars. The Parkes Pulsar Timing Array (PPTA) has been monitoring 20 millisecond pulsars at two to three weeks intervals for more than a decade. In addition to the detection of nanohertz gravitational waves, PPTA offers the opportunity for direct searches for fuzzy dark matter in an astrophysically feasible range of masses. We analyze the latest PPTA data set which includes timing observations for 26 pulsars made between 2004 and 2016. We perform a search in this data set for evidence of ultralight dark matter in the Galaxy using Bayesian and Frequentist methods. No statistically significant detection has been made. We therefore place upper limits on the local dark matter density. Our limits, improving on previous searches by a factor of two to five, constrain the dark matter density of ultralight bosons with $m \leq 10^{-23}$ eV to be below $6\,\text{GeV}\,\text{cm}^{-3}$ with 95\% confidence in the Earth neighborhood. Finally, we discuss the prospect of probing the astrophysically favored mass range $m \gtrsim 10^{-22}$ eV with next-generation pulsar timing facilities., Comment: Fixed a typo in Fig. 5; Published in Phys. Rev. D as Editor's Suggestion article, 17 pages, 7 figures, 4 tables

A pulsar timing array (PTA) refers to a program of regular, high-precision timing observations of a widely distributed array of millisecond pulsars. Here we review the status of the three primary PTA projects and the joint International Pulsar Timing Array project. We discuss current results related to ultra-low-frequency gravitational wave searches and highlight opportunities for the near future., Comment: 14 pages, 3 figures; National Science Review

The primary goal of the pulsar timing array projects is to detect ultra-low-frequency gravitational waves. The pulsar data sets are affected by numerous noise processes including varying dispersive delays in the interstellar medium and from the solar wind. The solar wind can lead to rapidly changing variations that, with existing telescopes, can be hard to measure and then remove. In this paper we study the possibility of using a low frequency telescope to aid in such correction for the Parkes Pulsar Timing Array (PPTA) and also discuss whether the ultra-wide-bandwidth receiver for the FAST telescope is sufficient to model the solar wind variations. Our key result is that a single wide-bandwidth receiver can be used to model and remove the effect of the solar wind. However, for pulsars that pass close to the Sun such as PSR J1022+1022, the solar wind is so variable that observations at two telescopes separated by a day are insufficient to correct the solar wind effect., Comment: accepted by RAA

We have observed the Vela pulsar for one year using a Phased Array Feed (PAF) receiver on the 12-metre antenna of the Parkes Test-Bed Facility. These observations have allowed us to investigate the stability of the PAF beam-weights over time, to demonstrate that pulsars can be timed over long periods using PAF technology and to detect and study the most recent glitch event that occurred on 12 December 2016. The beam-weights are shown to be stable to 1% on time scales on the order of three weeks. We discuss the implications of this for monitoring pulsars using PAFs on single dish telescopes., Comment: 6 pages, 4 figures, 2 tables. Accepted for publication in PASA

A pulsar timing array is a Galactic-scale detector of nanohertz gravitational waves (GWs). Its target signals contain two components: the `Earth term' and the `pulsar term' corresponding to GWs incident on the Earth and pulsar respectively. In this work we present a Frequentist method for the detection and localization of continuous waves that takes into account the pulsar term and is significantly faster than existing methods. We investigate the role of pulsar terms by comparing a full-signal search with an Earth-term-only search for non-evolving black hole binaries. By applying the method to synthetic data sets, we find that (i) a full-signal search can slightly improve the detection probability (by about five percent); (ii) sky localization is biased if only Earth terms are searched for and the inclusion of pulsar terms is critical to remove such a bias; (iii) in the case of strong detections (with signal-to-noise ratio $\gtrsim$ 30), it may be possible to improve pulsar distance estimation through GW measurements., Comment: 12 pages, 9 figures, typos corrected. To match the published version. Code implementing this method is available at the PPTA Wiki page

We report the discovery of two new millisecond pulsars (PSRs J0024$-$7204aa and J0024$-$7204ab) in the globular cluster 47\,Tucanae (NGC 104). Our results bring the total number of pulsars in 47\,Tucanae to 25. These pulsars were discovered by reprocessing archival observations from the Parkes radio telescope. We reprocessed the data using a standard search procedure based on the PRESTO software package as well as using a new method in which we incoherently added the power spectra corresponding to $\sim$1100\,hr of observations. The newly discovered PSR~J0024$-$7204aa, has a pulse frequency of $\rm \sim$541\,Hz (corresponding to a $\rm \sim$1.84 ms period), which is higher than any other pulsars currently known in the cluster and ranks 12$^{\rm{th}}$ amongst all the currently known pulsars. The dispersion measure of this pulsar, 24.941(7)\,cm$^{-3}$ pc, is the highest in the cluster. The second discovered pulsar, PSR~J0024$-$7204ab, is an isolated pulsar with a pulse frequency of $\rm \sim$270\,Hz (corresponding to a period of $\rm \sim$3.70 ms)., Comment: 5 pages, 1 table, 5 figures, accepted by MNRAS letter

In the centenary year of Einstein's General Theory of Relativity, this paper reviews the current status of gravitational wave astronomy across a spectrum which stretches from attohertz to kilohertz frequencies. Sect. 1 of this paper reviews the historical development of gravitational wave astronomy from Einstein's first prediction to our current understanding the spectrum. It is shown that detection of signals in the audio frequency spectrum can be expected very soon, and that a north-south pair of next generation detectors would provide large scientific benefits. Sect. 2 reviews the theory of gravitational waves and the principles of detection using laser interferometry. The state of the art Advanced LIGO detectors are then described. These detectors have a high chance of detecting the first events in the near future. Sect. 3 reviews the KAGRA detector currently under development in Japan, which will be the first laser interferometer detector to use cryogenic test masses. Sect. 4 of this paper reviews gravitational wave detection in the nanohertz frequency band using the technique of pulsar timing. Sect. 5 reviews the status of gravitational wave detection in the attohertz frequency band, detectable in the polarisation of the cosmic microwave background, and discusses the prospects for detection of primordial waves from the big bang. The techniques described in sects. 1-5 have already placed significant limits on the strength of gravitational wave sources. Sects. 6 and 7 review ambitious plans for future space based gravitational wave detectors in the millihertz frequency band. Sect. 6 presents a roadmap for development of space based gravitational wave detectors by China while sect. 7 discusses a key enabling technology for space interferometry known as time delay interferometry., Comment: published (2015)

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., Comment: accepted for publication in Physical Review X

In a search for periodic variation in the arrival times of pulses from 151 young, energetic pulsars, we have identified seven cases of modulation consistent with one or two harmonics of a single fundamental with time-scale 0.5-1.5 yr. We use simulations to show that these modulations are statistically significant and of high quality (sinusoidal) even when contaminated by the strong stochastic timing noise common to young pulsars. Although planetary companions could induce such modulation, the large implied masses and 2:1 mean motion resonances challenge such an explanation. Instead, the modulation is likely to be intrinsic to the pulsar, arising from quasi-periodic switching between stable magnetospheric states, and we propose that precession of the neutron star may regulate this switching., Comment: 11 pages, 8 figures, 3 tables, accepted in MNRAS

Pulsar timing array experiments search for phenomena that produce angular correlations in the arrival times of signals from millisecond pulsars. The primary goal is to detect an isotropic and stochastic gravitational wave background. We use simulated data to show that this search can be affected by the presence of other spatially correlated noise, such as errors in the reference time standard, errors in the planetary ephemeris, the solar wind and instrumentation issues. All these effects can induce significant false detections of gravitational waves. We test mitigation routines to account for clock errors, ephemeris errors and the solar wind. We demonstrate that it is non-trivial to find an effective mitigation routine for the planetary ephemeris and emphasise that other spatially correlated signals may be present in the data., Comment: 13 pages, 7 figures. Accepted for publication in MNRAS

Complementary to ground-based laser interferometers, pulsar timing array experiments are being carried out to search for nanohertz gravitational waves. Using the world's most powerful radio telescopes, three major international collaborations have collected $\sim$10-year high precision timing data for tens of millisecond pulsars. In this paper we give an overview on pulsar timing experiments, gravitational wave detection in the nanohertz regime, and recent results obtained by various timing array projects., Comment: 11 pages, 5 figures, review article for the special issue of SCIENCE CHINA Physics, Mechanics & Astronomy to mark the Centenary of General Relativity

We have searched a sample of 151 young, energetic pulsars for periodic variation in pulse time-of-arrival arising from the influence of planetary companions. We are sensitive to objects with masses two orders of magnitude lower than those detectable with optical transit timing, but we find no compelling evidence for pulsar planets. For the older pulsars most likely to host planets, we can rule out Mercury analogues in one third of our sample and planets with masses $>0.4 M_{\oplus}$ and periods $P_b<1$ yr in all but 5% of such systems. If pulsar planets form primarily from supernova fallback disks, these limits imply that such disks do not form, are confined to $<0.1$ AU radii, are disrupted, or form planets more slowly ($>2$ Myr) than their protoplanetary counterparts., Comment: 5 pages, 4 figures, accepted to ApJL

The sensitivity of pulsar timing arrays to gravitational waves is, at some level, limited by timing noise. Red timing noise - the stochastic wandering of pulse arrival times with a red spectrum - is prevalent in slow-spinning pulsars and has been identified in many millisecond pulsars. Phenomenological models of timing noise, such as from superfluid turbulence, suggest that the timing noise spectrum plateaus below some critical frequency, $f_c$, potentially aiding the hunt for gravitational waves. We examine this effect for individual pulsars by calculating minimum observation times, $T_{\rm min}(f_c)$, over which the gravitational wave signal becomes larger than the timing noise plateau. We do this in two ways: 1) in a model-independent manner, and 2) by using the superfluid turbulence model for timing noise as an example to illustrate how neutron star parameters can be constrained. We show that the superfluid turbulence model can reproduce the data qualitatively from a number of pulsars observed as part of the Parkes Pulsar Timing Array. We further show how a value of $f_c$, derived either through observations or theory, can be related to $T_{\rm min}$. This provides a diagnostic whereby the usefulness of timing array pulsars for gravitational-wave detection can be quantified., Comment: Accepted for publication in MNRAS

Pulsar timing arrays (PTAs) can be used to search for very low frequency ($10^{-9}$--$10^{-7}$ Hz) gravitational waves (GWs). In this paper we present a general method for the detection and localization of single-source GWs using PTAs. We demonstrate the effectiveness of this new method for three types of signals: monochromatic waves as expected from individual supermassive binary black holes in circular orbits, GWs from eccentric binaries and GW bursts. We also test its implementation in realistic data sets that include effects such as uneven sampling and heterogeneous data spans and measurement precision. It is shown that our method, which works in the frequency domain, performs as well as published time-domain methods. In particular, we find it equivalent to the $\mathcal{F}_{e}$-statistic for monochromatic waves. We also discuss the construction of null streams -- data streams that have null response to GWs, and the prospect of using null streams as a consistency check in the case of detected GW signals. Finally, we present sensitivities to individual supermassive binary black holes in eccentric orbits. We find that a monochromatic search that is designed for circular binaries can efficiently detect eccentric binaries with both high and low eccentricities, while a harmonic summing technique provides greater sensitivities only for binaries with moderate eccentricities., Comment: To match the published version, 17 pages, 13 figures