-
The neutron star mass, distance, and inclination from precision timing of the brilliant millisecond pulsar J0437$-$4715
Authors:
Daniel J. Reardon,
Matthew Bailes,
Ryan M. Shannon,
Chris Flynn,
Jacob Askew,
N. D. Ramesh Bhat,
Zu-Cheng Chen,
Małgorzata Curyło,
Yi Feng,
George B. Hobbs,
Agastya Kapur,
Matthew Kerr,
Xiaojin Liu,
Richard N. Manchester,
Rami Mandow,
Saurav Mishra,
Christopher J. Russell,
Mohsen Shamohammadi,
Lei Zhang,
Andrew Zic
Abstract:
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 Co…
▽ More
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.
△ Less
Submitted 30 July, 2024; v1 submitted 9 July, 2024;
originally announced July 2024.
-
A NICER View of the Nearest and Brightest Millisecond Pulsar: PSR J0437$\unicode{x2013}$4715
Authors:
Devarshi Choudhury,
Tuomo Salmi,
Serena Vinciguerra,
Thomas E. Riley,
Yves Kini,
Anna L. Watts,
Bas Dorsman,
Slavko Bogdanov,
Sebastien Guillot,
Paul S. Ray,
Daniel J. Reardon,
Ronald A. Remillard,
Anna V. Bilous,
Daniela Huppenkothen,
James M. Lattimer,
Nathan Rutherford,
Zaven Arzoumanian,
Keith C. Gendreau,
Sharon M. Morsink,
Wynn C. G. Ho
Abstract:
We report Bayesian inference of the mass, radius and hot X-ray emitting region properties - using data from the Neutron Star Interior Composition ExploreR (NICER) - for the brightest rotation-powered millisecond X-ray pulsar PSR J0437$\unicode{x2013}$4715. Our modeling is conditional on informative tight priors on mass, distance and binary inclination obtained from radio pulsar timing using the Pa…
▽ More
We report Bayesian inference of the mass, radius and hot X-ray emitting region properties - using data from the Neutron Star Interior Composition ExploreR (NICER) - for the brightest rotation-powered millisecond X-ray pulsar PSR J0437$\unicode{x2013}$4715. Our modeling is conditional on informative tight priors on mass, distance and binary inclination obtained from radio pulsar timing using the Parkes Pulsar Timing Array (PPTA) (Reardon et al. 2024), and we use NICER background models to constrain the non-source background, cross-checking with data from XMM-Newton. We assume two distinct hot emitting regions, and various parameterized hot region geometries that are defined in terms of overlapping circles; while simplified, these capture many of the possibilities suggested by detailed modeling of return current heating. For the preferred model identified by our analysis we infer a mass of $M = 1.418 \pm 0.037$ M$_\odot$ (largely informed by the PPTA mass prior) and an equatorial radius of $R = 11.36^{+0.95}_{-0.63}$ km, each reported as the posterior credible interval bounded by the 16% and 84% quantiles. This radius favors softer dense matter equations of state and is highly consistent with constraints derived from gravitational wave measurements of neutron star binary mergers. The hot regions are inferred to be non-antipodal, and hence inconsistent with a pure centered dipole magnetic field.
△ Less
Submitted 9 July, 2024;
originally announced July 2024.
-
Status Report on Global Pulsar-Timing-Array Efforts to Detect Gravitational Waves
Authors:
Joris P. W. Verbiest,
Sarah J. Vigeland,
Nataliya K. Porayko,
Siyuan Chen,
Daniel J. Reardon
Abstract:
The stability of the spin of pulsars and the precision with which these spins can be determined, allows many unique tests of interest to physics and astrophysics. Perhaps the most challenging and revolutionary of these, is the detection of nanohertz gravitational waves. An increasing number of efforts to detect and study long-period gravitational waves by timing an array of pulsars have been ongoi…
▽ More
The stability of the spin of pulsars and the precision with which these spins can be determined, allows many unique tests of interest to physics and astrophysics. Perhaps the most challenging and revolutionary of these, is the detection of nanohertz gravitational waves. An increasing number of efforts to detect and study long-period gravitational waves by timing an array of pulsars have been ongoing for several decades and the field is moving ever closer to actual gravitational-wave science. In this review article, we summarise the state of this field by presenting the current sensitivity to gravitational waves and by reviewing recent progress along the multiple lines of research that are part of the continuous push towards greater sensitivity. We also briefly review some of the most recent efforts at astrophysical interpretation of the most recent GW estimates derived from pulsar timing.
△ Less
Submitted 30 April, 2024;
originally announced April 2024.
-
MeerKAT Pulsar Timing Array parallaxes and proper motions
Authors:
Mohsen Shamohammadi,
Matthew Bailes,
Christopher Flynn,
Daniel J. Reardon,
Ryan M. Shannon,
Sarah Buchner,
Andrew D. Cameron,
Fernando Camilo,
Alessandro Coronigu,
Marisa Geyer,
Michael Kramer,
Matthew Miles,
Renee Spiewak
Abstract:
We have determined positions, proper motions, and parallaxes of $77$ millisecond pulsars (MSPs) from $\sim3$ years of MeerKAT radio telescope observations. Our timing and noise analyses enable us to measure $35$ significant parallaxes ($12$ of them for the first time) and $69$ significant proper motions. Eight pulsars near the ecliptic have an accurate proper motion in ecliptic longitude only. PSR…
▽ More
We have determined positions, proper motions, and parallaxes of $77$ millisecond pulsars (MSPs) from $\sim3$ years of MeerKAT radio telescope observations. Our timing and noise analyses enable us to measure $35$ significant parallaxes ($12$ of them for the first time) and $69$ significant proper motions. Eight pulsars near the ecliptic have an accurate proper motion in ecliptic longitude only. PSR~J0955$-$6150 has a good upper limit on its very small proper motion ($<$0.4 mas yr$^{-1}$). We used pulsars with accurate parallaxes to study the MSP velocities. This yields $39$ MSP transverse velocities, and combined with MSPs in the literature (excluding those in Globular Clusters) we analyse $66$ MSPs in total. We find that MSPs have, on average, much lower velocities than normal pulsars, with a mean transverse velocity of only $78(8)$ km s$^{-1}$ (MSPs) compared with $246(21)$ km s$^{-1}$ (normal pulsars). We found no statistical differences between the velocity distributions of isolated and binary millisecond pulsars. From Galactocentric cylindrical velocities of the MSPs, we derive 3-D velocity dispersions of $σ_ρ$, $σ_φ$, $σ_{z}$ = $63(11)$, $48(8)$, $19(3)$ km s$^{-1}$. We measure a mean asymmetric drift with amplitude $38(11)$ km s$^{-1}$, consistent with expectation for MSPs, given their velocity dispersions and ages. The MSP velocity distribution is consistent with binary evolution models that predict very few MSPs with velocities $>300$ km s$^{-1}$ and a mild anticorrelation of transverse velocity with orbital period.
△ Less
Submitted 12 January, 2024;
originally announced January 2024.
-
An insight into chromatic behaviour of jitter in pulsars and its modelling: A case study of PSR J0437$-$4715
Authors:
A. D. Kulkarni,
R. M. Shannon,
D. J. Reardon,
M. T. Miles,
M. Bailes,
M. Shamohammadi
Abstract:
Pulse-to-pulse profile shape variations introduce correlations in pulsar times of arrival (TOAs) across radio frequency measured at the same observational epoch. This leads to a broadband noise in excess of radiometer noise, which is termed pulse jitter noise. The presence of jitter noise limits the achievable timing precision and decreases the sensitivity of pulsar-timing data sets to signals of…
▽ More
Pulse-to-pulse profile shape variations introduce correlations in pulsar times of arrival (TOAs) across radio frequency measured at the same observational epoch. This leads to a broadband noise in excess of radiometer noise, which is termed pulse jitter noise. The presence of jitter noise limits the achievable timing precision and decreases the sensitivity of pulsar-timing data sets to signals of interest such as nanohertz-frequency gravitational waves. Current white noise models used in pulsar timing analyses attempt to account for this, assuming complete correlation of uncertainties through the arrival times collected in a unique observation and no frequency dependence of jitter (which corresponds to a rank-one covariance matrix). However, previous studies show that the brightest millisecond pulsar at decimetre wavelengths, PSR J0437$-$4715, shows decorrelation and frequency dependence of jitter noise. Here we present a detailed study of the decorrelation of jitter noise in PSR J0437$-$4715 and implement a new technique to model it. We show that the rate of decorrelation due to jitter can be expressed as a power-law in frequency. We analyse the covariance matrix associated with the jitter noise process and find that a higher-rank-approximation is essential to account for the decorrelation and to account for frequency dependence of jitter noise. We show that the use of this novel method significantly improves the estimation of other chromatic noise parameters such as dispersion measure variations. However, we find no significant improvement in errors and estimation of other timing model parameters suggesting that current methods are not biased for other parameters, for this pulsar due to this misspecification. We show that pulse energy variations show a similar decorrelation to the jitter noise, indicating a common origin for both observables.
△ Less
Submitted 8 January, 2024; v1 submitted 7 January, 2024;
originally announced January 2024.
-
Precision constraints on the neutron star equation of state with third-generation gravitational-wave observatories
Authors:
Kris Walker,
Rory Smith,
Eric Thrane,
Daniel J. Reardon
Abstract:
It is currently unknown how matter behaves at the extreme densities found within the cores of neutron stars. Measurements of the neutron star equation of state probe nuclear physics that is otherwise inaccessible in a laboratory setting. Gravitational waves from binary neutron star mergers encode details about this physics, allowing the equation of state to be inferred. Planned third-generation gr…
▽ More
It is currently unknown how matter behaves at the extreme densities found within the cores of neutron stars. Measurements of the neutron star equation of state probe nuclear physics that is otherwise inaccessible in a laboratory setting. Gravitational waves from binary neutron star mergers encode details about this physics, allowing the equation of state to be inferred. Planned third-generation gravitational-wave observatories, having vastly improved sensitivity, are expected to provide tight constraints on the neutron star equation of state. We combine simulated observations of binary neutron star mergers by the third-generation observatories Cosmic Explorer and Einstein Telescope to determine future constraints on the equation of state across a plausible neutron star mass range. In one year of operation, a network consisting of one Cosmic Explorer and the Einstein Telescope is expected to detect $\gtrsim 3\times 10^5$ binary neutron star mergers. By considering only the 75 loudest events, we show that such a network will be able to constrain the neutron star radius to at least $\lesssim 200$ m (90% credibility) in the mass range $1-1.97$ $M_{\odot}$ -- about ten times better than current constraints from LIGO-Virgo-KAGRA and NICER. The constraint is $\lesssim 75$ m (90% credibility) near $1.4-1.6$ $M_{\odot}$ where we assume the the binary neutron star mass distribution is peaked. This constraint is driven primarily from the loudest $\sim 20$ events.
△ Less
Submitted 4 January, 2024;
originally announced January 2024.
-
A MeerKAT view of the double pulsar eclipses -- Geodetic precession of pulsar B and system geometry
Authors:
M. E. Lower,
M. Kramer,
R. M. Shannon,
R. P. Breton,
N. Wex,
S. Johnston,
M. Bailes,
S. Buchner,
H. Hu,
V. Venkatraman Krishnan,
V. A. Blackmon,
F. Camilo,
D. J. Champion,
P. C. C. Freire,
M. Geyer,
A. Karastergiou,
J. van Leeuwen,
M. A. McLaughlin,
D. J. Reardon,
I. H. Stairs
Abstract:
The double pulsar system, PSR J0737$-$3039A/B, consists of two neutron stars bound together in a highly relativistic orbit that is viewed nearly edge-on from the Earth. This alignment results in brief radio eclipses of the fast-rotating pulsar A when it passes behind the toroidal magnetosphere of the slow-rotating pulsar B. The morphology of these eclipses is strongly dependent on the geometric or…
▽ More
The double pulsar system, PSR J0737$-$3039A/B, consists of two neutron stars bound together in a highly relativistic orbit that is viewed nearly edge-on from the Earth. This alignment results in brief radio eclipses of the fast-rotating pulsar A when it passes behind the toroidal magnetosphere of the slow-rotating pulsar B. The morphology of these eclipses is strongly dependent on the geometric orientation and rotation phase of pulsar B, and their time-evolution can be used to constrain the geodetic precession rate of the pulsar. We demonstrate a Bayesian inference framework for modelling eclipse light-curves obtained with MeerKAT between 2019-2023. Using a hierarchical inference approach, we obtained a precession rate of $Ω_{\rm SO}^{\rm B} = {5.16^{\circ}}^{+0.32^{\circ}}_{-0.34^{\circ}}$ yr$^{-1}$ for pulsar B, consistent with predictions from General Relativity to a relative uncertainty of 6.5%. This updated measurement provides a 6.1% test of relativistic spin-orbit coupling in the strong-field regime. We show that a simultaneous fit to all of our observed eclipses can in principle return a $\sim$1.5% test of spin-orbit coupling. However, systematic effects introduced by the current geometric orientation of pulsar B along with inconsistencies between the observed and predicted eclipse light curves result in difficult to quantify uncertainties. Assuming the validity of General Relativity, we definitively show that the spin-axis of pulsar B is misaligned from the total angular momentum vector by $40.6^{\circ} \pm 0.1^{\circ}$ and that the orbit of the system is inclined by approximately $90.5^{\circ}$ from the direction of our line of sight. Our measured geometry for pulsar B suggests the largely empty emission cone contains an elongated horseshoe shaped beam centered on the magnetic axis, and that it may not be re-detected as a radio pulsar until early-2035.
△ Less
Submitted 1 February, 2024; v1 submitted 10 November, 2023;
originally announced November 2023.
-
Pulsar Scintillation through Thick and Thin: Bow Shocks, Bubbles, and the Broader Interstellar Medium
Authors:
S. K. Ocker,
J. M. Cordes,
S. Chatterjee,
D. R. Stinebring,
T. Dolch,
V. Pelgrims,
J. W. McKee,
C. Giannakopoulos,
D. J. Reardon
Abstract:
Observations of pulsar scintillation are among the few astrophysical probes of very small-scale ($\lesssim$ au) phenomena in the interstellar medium (ISM). In particular, characterization of scintillation arcs, including their curvature and intensity distributions, can be related to interstellar turbulence and potentially over-pressurized plasma in local ISM inhomogeneities, such as supernova remn…
▽ More
Observations of pulsar scintillation are among the few astrophysical probes of very small-scale ($\lesssim$ au) phenomena in the interstellar medium (ISM). In particular, characterization of scintillation arcs, including their curvature and intensity distributions, can be related to interstellar turbulence and potentially over-pressurized plasma in local ISM inhomogeneities, such as supernova remnants, HII regions, and bow shocks. Here we present a survey of eight pulsars conducted at the Five-hundred-meter Aperture Spherical Telescope (FAST), revealing a diverse range of scintillation arc characteristics at high sensitivity. These observations reveal more arcs than measured previously for our sample. At least nine arcs are observed toward B1929$+$10 at screen distances spanning $\sim 90\%$ of the pulsar's $361$ pc path-length to the observer. Four arcs are observed toward B0355$+$54, with one arc yielding a screen distance as close as $\sim10^5$ au ($<1$ pc) from either the pulsar or the observer. Several pulsars show highly truncated, low-curvature arcs that may be attributable to scattering near the pulsar. The scattering screen constraints are synthesized with continuum maps of the local ISM and other well-characterized pulsar scintillation arcs, yielding a three-dimensional view of the scattering media in context.
△ Less
Submitted 27 November, 2023; v1 submitted 24 September, 2023;
originally announced September 2023.
-
Flux density monitoring of 89 millisecond pulsars with MeerKAT
Authors:
P. Gitika,
M. Bailes,
R. M. Shannon,
D. J. Reardon,
A. D. Cameron,
M. Shamohammadi,
M. T. Miles,
C. M. L. Flynn,
A. Corongiu,
M. Kramer
Abstract:
We present a flux density study of 89 millisecond pulsars (MSPs) regularly monitored as part of the MeerKAT Pulsar Timing Array (MPTA) using the L-Band receiver with an approximately two week cadence between 2019-2022. For each pulsar, we have determined the mean flux densities at each epoch in eight $\sim$97 MHz sub-bands ranging from 944 to 1625 MHz. From these we have derived their modulation i…
▽ More
We present a flux density study of 89 millisecond pulsars (MSPs) regularly monitored as part of the MeerKAT Pulsar Timing Array (MPTA) using the L-Band receiver with an approximately two week cadence between 2019-2022. For each pulsar, we have determined the mean flux densities at each epoch in eight $\sim$97 MHz sub-bands ranging from 944 to 1625 MHz. From these we have derived their modulation indices, their average and peak-to-median flux densities in each sub-band, as well as their mean spectral indices across the entire frequency range. We find that the vast majority of the MSPs have spectra that are well described by a simple power law, with a mean spectral index of -1.86(6). Using the temporal variation of the flux densities we measured the structure functions and determined the refractive scintillation timescale for seven. The structure functions provide strong evidence that the intrinsic radio luminosities of MSPs are stable. As a population, the average modulation index at 20 cm wavelengths peaks near unity at dispersion measures (DMs) of $\sim$20 pc cm$^{-3}$ and by a DM of 100 pc cm$^{-3}$ are closer to 0.2, due to refractive scintillation. We find that timing arrays can improve their observing efficiency by reacting to scintillation maxima, and that 20 cm FRB surveys should prioritise highly scintillating mid-latitude regions of the Galactic sky where they will find $\sim$30% more events and bursts at greater distances.
△ Less
Submitted 14 September, 2023;
originally announced September 2023.
-
Comparing recent PTA results on the nanohertz stochastic gravitational wave background
Authors:
The International Pulsar Timing Array Collaboration,
G. Agazie,
J. Antoniadis,
A. Anumarlapudi,
A. M. Archibald,
P. Arumugam,
S. Arumugam,
Z. Arzoumanian,
J. Askew,
S. Babak,
M. Bagchi,
M. Bailes,
A. -S. Bak Nielsen,
P. T. Baker,
C. G. Bassa,
A. Bathula,
B. Bécsy,
A. Berthereau,
N. D. R. Bhat,
L. Blecha,
M. Bonetti,
E. Bortolas,
A. Brazier,
P. R. Brook,
M. Burgay
, et al. (220 additional authors not shown)
Abstract:
The Australian, Chinese, European, Indian, and North American pulsar timing array (PTA) collaborations recently reported, at varying levels, evidence for the presence of a nanohertz gravitational wave background (GWB). Given that each PTA made different choices in modeling their data, we perform a comparison of the GWB and individual pulsar noise parameters across the results reported from the PTA…
▽ More
The Australian, Chinese, European, Indian, and North American pulsar timing array (PTA) collaborations recently reported, at varying levels, evidence for the presence of a nanohertz gravitational wave background (GWB). Given that each PTA made different choices in modeling their data, we perform a comparison of the GWB and individual pulsar noise parameters across the results reported from the PTAs that constitute the International Pulsar Timing Array (IPTA). We show that despite making different modeling choices, there is no significant difference in the GWB parameters that are measured by the different PTAs, agreeing within $1σ$. The pulsar noise parameters are also consistent between different PTAs for the majority of the pulsars included in these analyses. We bridge the differences in modeling choices by adopting a standardized noise model for all pulsars and PTAs, finding that under this model there is a reduction in the tension in the pulsar noise parameters. As part of this reanalysis, we "extended" each PTA's data set by adding extra pulsars that were not timed by that PTA. Under these extensions, we find better constraints on the GWB amplitude and a higher signal-to-noise ratio for the Hellings and Downs correlations. These extensions serve as a prelude to the benefits offered by a full combination of data across all pulsars in the IPTA, i.e., the IPTA's Data Release 3, which will involve not just adding in additional pulsars, but also including data from all three PTAs where any given pulsar is timed by more than as single PTA.
△ Less
Submitted 1 September, 2023;
originally announced September 2023.
-
The Parkes Pulsar Timing Array Third Data Release
Authors:
Andrew Zic,
Daniel J. Reardon,
Agastya Kapur,
George Hobbs,
Rami Mandow,
Małgorzata Curyło,
Ryan M. Shannon,
Jacob Askew,
Matthew Bailes,
N. D. Ramesh Bhat,
Andrew Cameron,
Zu-Cheng Chen,
Shi Dai,
Valentina Di Marco,
Yi Feng,
Matthew Kerr,
Atharva Kulkarni,
Marcus E. Lower,
Rui Luo,
Richard N. Manchester,
Matthew T. Miles,
Rowina S. Nathan,
Stefan Osłowski,
Axl F. Rogers,
Christopher J. Russell
, et al. (9 additional authors not shown)
Abstract:
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 ob…
▽ More
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.
△ Less
Submitted 17 October, 2023; v1 submitted 28 June, 2023;
originally announced June 2023.
-
The gravitational-wave background null hypothesis: Characterizing noise in millisecond pulsar arrival times with the Parkes Pulsar Timing Array
Authors:
Daniel J. Reardon,
Andrew Zic,
Ryan M. Shannon,
Valentina Di Marco,
George B. Hobbs,
Agastya Kapur,
Marcus E. Lower,
Rami Mandow,
Hannah Middleton,
Matthew T. Miles,
Axl F. Rogers,
Jacob Askew,
Matthew Bailes,
N. D. Ramesh Bhat,
Andrew Cameron,
Matthew Kerr,
Atharva Kulkarni,
Richard N. Manchester,
Rowina S. Nathan,
Christopher J. Russell,
Stefan Osłowski,
Xing-Jiang Zhu
Abstract:
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 gravita…
▽ More
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 $γ=13/3$, the value predicted for a stochastic GW background from a population of supermassive black hole binaries driven solely by GW emission.
△ Less
Submitted 28 June, 2023;
originally announced June 2023.
-
Search for an isotropic gravitational-wave background with the Parkes Pulsar Timing Array
Authors:
Daniel J. Reardon,
Andrew Zic,
Ryan M. Shannon,
George B. Hobbs,
Matthew Bailes,
Valentina Di Marco,
Agastya Kapur,
Axl F. Rogers,
Eric Thrane,
Jacob Askew,
N. D. Ramesh Bhat,
Andrew Cameron,
Małgorzata Curyło,
William A. Coles,
Shi Dai,
Boris Goncharov,
Matthew Kerr,
Atharva Kulkarni,
Yuri Levin,
Marcus E. Lower,
Richard N. Manchester,
Rami Mandow,
Matthew T. Miles,
Rowina S. Nathan,
Stefan Osłowski
, et al. (4 additional authors not shown)
Abstract:
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 t…
▽ More
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})^α$, we measure $A=3.1^{+1.3}_{-0.9} \times 10^{-15}$ and $α=-0.45 \pm 0.20$ respectively (median and 68% credible interval). For a spectral index of $α=-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 $α=-2/3$, we measure spatial correlations consistent with a GWB, with an estimated false-alarm probability of $p \lesssim 0.02$ (approx. $2σ$). 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.
△ Less
Submitted 28 June, 2023;
originally announced June 2023.
-
Toward robust detections of nanohertz gravitational waves
Authors:
Valentina Di Marco,
Andrew Zic,
Matthew T. Miles,
Daniel J. Reardon,
Eric Thrane,
Ryan M. Shannon
Abstract:
The recent observation of a common red-noise process in pulsar timing arrays (PTAs) suggests that the detection of nanohertz gravitational waves might be around the corner. However, in order to confidently attribute this red process to gravitational waves, one must observe the Hellings-Downs curve -- the telltale angular correlation function associated with a gravitational-wave background. This ef…
▽ More
The recent observation of a common red-noise process in pulsar timing arrays (PTAs) suggests that the detection of nanohertz gravitational waves might be around the corner. However, in order to confidently attribute this red process to gravitational waves, one must observe the Hellings-Downs curve -- the telltale angular correlation function associated with a gravitational-wave background. This effort is complicated by the complex modelling of pulsar noise. Without proper care, mis-specified noise models can lead to false-positive detections. Background estimation using bootstrap methods such as sky scrambles and phase shifts, which use the data to characterize the noise, are therefore important tools for assessing significance. We investigate the ability of current PTA experiments to estimate their background with "quasi-independent" scrambles -- characterized by a statistical "match" below the fiducial value: $|M|<0.1$. We show that sky scrambling is affected by "saturation" after $O(10)$ quasi-independent realizations; subsequent scrambles are no longer quasi-independent. We show phase scrambling saturates after $O(100)$ quasi-independent realizations. With so few independent scrambles, it is difficult to make reliable statements about the $\gtrsim 5 σ$ tail of the null distribution of the detection statistic. We discuss various methods by which one may increase the number of independent scrambles. We also consider an alternative approach wherein one re-frames the background estimation problem so that the significance is calculated using statistically dependent scrambles. The resulting $p$-value is in principle well-defined but may be susceptible to failure if assumptions about the data are incorrect.
△ Less
Submitted 9 August, 2023; v1 submitted 8 May, 2023;
originally announced May 2023.
-
Mass measurements and 3D orbital geometry of PSR J1933$-$6211
Authors:
M. Geyer,
V. Venkatraman Krishnan,
P. C. C. Freire,
M. Kramer,
J. Antoniadis,
M. Bailes,
M. C. i Bernadich,
S. Buchner,
A. D. Cameron,
D. J. Champion,
A. Karastergiou,
M. J. Keith,
M. E. Lower,
S. Osłowski,
A. Possenti,
A. Parthasarathy,
D. J. Reardon,
M. Serylak,
R. M. Shannon,
R. Spiewak,
W. van Straten,
J. P. W. Verbiest
Abstract:
PSR J1933$-$6211 is a 3.5-ms pulsar in a 12.8-d orbit with a white dwarf (WD). Its high proper motion and low dispersion measure result in such significant interstellar scintillation that high signal-to-noise detections require long observing durations or fortuitous timing. We turn to the sensitive MeerKAT telescope and, combined with historic Parkes data, leverage PSR J1933$-$6211's kinematic and…
▽ More
PSR J1933$-$6211 is a 3.5-ms pulsar in a 12.8-d orbit with a white dwarf (WD). Its high proper motion and low dispersion measure result in such significant interstellar scintillation that high signal-to-noise detections require long observing durations or fortuitous timing. We turn to the sensitive MeerKAT telescope and, combined with historic Parkes data, leverage PSR J1933$-$6211's kinematic and relativistic effects to constrain its 3D orbital geometry and the component masses. We obtain precise proper motion and parallax estimates, and measure their effects as secular changes in the Keplerian orbital parameters: a variation in orbital period of $7(1) \times 10^{-13}$ s s$^{-1}$ and a change in projected semi-major axis of $1.60(5) \times 10^{-14}$ s s$^{-1}$. A self-consistent analysis of all kinematic and relativistic effects yields a distance of $1.6^{+0.2}_{-0.3}$ kpc, an orbital inclination, $i = 55(1)$ deg and a longitude of the ascending node, $Ω= 255^{+8}_{-14}$ deg. The probability densities for $Ω$ and $i$ and their symmetric counterparts, ($180-i$, $360-Ω$), are seen to depend on the fiducial orbit used to measure the time of periastron passage. We investigate this unexpected dependence and rule out software-related causes using simulations. Nevertheless, we constrain the pulsar and WD masses to $1.4^{+0.3}_{-0.2}$ M$_\odot$ and $0.43(5)$ M$_\odot$ respectively. These strongly disfavour a helium-dominated WD. The orbital similarities between PSRs J1933$-$6211 and J1614$-$2230 suggest they underwent Case A Roche lobe overflow, an extended evolution while the companion star is still on the Main Sequence. However, with a mass of $\sim 1.4$ M$_\odot$, PSR J1933$-$6211 has not accreted significant matter. This highlights the low accretion efficiency of the spin-up process and suggests that observed neutron star masses are mostly a result of supernova physics.
△ Less
Submitted 18 April, 2023;
originally announced April 2023.
-
Improving pulsar-timing solutions through dynamic pulse fitting
Authors:
Rowina S. Nathan,
Matthew T. Miles,
Gregory Ashton,
Paul D. Lasky,
Eric Thrane,
Daniel J. Reardon,
Ryan M. Shannon,
Andrew D. Cameron
Abstract:
Precision pulsar timing is integral to the detection of the nanohertz stochastic gravitational-wave background as well as understanding the physics of neutron stars. Conventional pulsar timing often uses fixed time and frequency-averaged templates to determine the pulse times of arrival, which can lead to reduced accuracy when the pulse profile evolves over time. We illustrate a dynamic timing met…
▽ More
Precision pulsar timing is integral to the detection of the nanohertz stochastic gravitational-wave background as well as understanding the physics of neutron stars. Conventional pulsar timing often uses fixed time and frequency-averaged templates to determine the pulse times of arrival, which can lead to reduced accuracy when the pulse profile evolves over time. We illustrate a dynamic timing method that fits each observing epoch using basis functions. By fitting each epoch separately, we allow for the evolution of the pulse shape epoch to epoch. We apply our method to PSR J1103$-$5403 and demonstrate that it undergoes mode changing, making it the fourth millisecond pulsar to exhibit such behaviour. Our method, which is able to identify and time a single mode, yields a timing solution with a root-mean-square error of 1.343 $μ\mathrm{s}$, a factor of 1.78 improvement over template fitting on both modes. In addition, the white-noise amplitude is reduced 4.3 times, suggesting that fitting the full data set causes the mode changing to be incorrectly classified as white noise. This reduction in white noise boosts the signal-to-noise ratio of a gravitational-wave background signal for this particular pulsar by 32%. We discuss the possible applications for this method of timing to study pulsar magnetospheres and further improve the sensitivity of searches for nanohertz gravitational waves.
△ Less
Submitted 31 May, 2023; v1 submitted 5 April, 2023;
originally announced April 2023.
-
Determining electron column density fluctuations in a dominant scattering region using pulsar scintillation
Authors:
Daniel J. Reardon,
William A. Coles
Abstract:
Density fluctuations in the ionised interstellar medium have a profound effect on radio pulsar observations, through angular scattering, intensity scintillations, and small changes in time delays from dispersion. Here we show that it is possible to recover the variations in dispersive delays that originate from a dominant scattering region using measurements of the dynamic spectrum of intensity sc…
▽ More
Density fluctuations in the ionised interstellar medium have a profound effect on radio pulsar observations, through angular scattering, intensity scintillations, and small changes in time delays from dispersion. Here we show that it is possible to recover the variations in dispersive delays that originate from a dominant scattering region using measurements of the dynamic spectrum of intensity scintillations, provided that the pulsar velocity and scattering region location are known. We provide a theoretical framework for the technique, which involves estimating the phase gradient from the dynamic spectra and integrating that gradient to obtain phase variations. It can be used to search for "extreme scattering events" (ESEs) in pulsars for which precision dispersion delay measurements are not otherwise possible, or to separate true dispersion variations from apparent variability caused by frequency-dependent pulse shape changes. We demonstrate that it works in practice by recovering an ESE in PSR J1603$-$7202, which is known from precision dispersion delay measurements from pulsar timing. For this pulsar, we find that the phase gradients also track the long-term variations in electron column density observed by pulsar timing, indicating that the column density variations and the scattering are dominated by the same thin scattering screen. We identify a sudden increase in the scintillation strength and magnitude of phase gradients over $\sim$days in 2010, indicating a compact structure. A decrease in the electron density in 2012 was associated with persistent phase gradients and preceded a period of decreased scintillation strength and an absence of scintillation arcs.
△ Less
Submitted 28 March, 2023;
originally announced March 2023.
-
Searching for continuous Gravitational Waves in the second data release of the International Pulsar Timing Array
Authors:
M. Falxa,
S. Babak,
P. T. Baker,
B. Bécsy,
A. Chalumeau,
S. Chen,
Z. Chen,
N. J. Cornish,
L. Guillemot,
J. S. Hazboun,
C. M. F. Mingarelli,
A. Parthasarathy,
A. Petiteau,
N. S. Pol,
A. Sesana,
S. B. Spolaor,
S. R. Taylor,
G. Theureau,
M. Vallisneri,
S. J. Vigeland,
C. A. Witt,
X. Zhu,
J. Antoniadis,
Z. Arzoumanian,
M. Bailes
, et al. (102 additional authors not shown)
Abstract:
The International Pulsar Timing Array 2nd data release is the combination of datasets from worldwide collaborations. In this study, we search for continuous waves: gravitational wave signals produced by individual supermassive black hole binaries in the local universe. We consider binaries on circular orbits and neglect the evolution of orbital frequency over the observational span. We find no evi…
▽ More
The International Pulsar Timing Array 2nd data release is the combination of datasets from worldwide collaborations. In this study, we search for continuous waves: gravitational wave signals produced by individual supermassive black hole binaries in the local universe. We consider binaries on circular orbits and neglect the evolution of orbital frequency over the observational span. We find no evidence for such signals and set sky averaged 95% upper limits on their amplitude h 95 . The most sensitive frequency is 10nHz with h 95 = 9.1 10-15 . We achieved the best upper limit to date at low and high frequencies of the PTA band thanks to improved effective cadence of observations. In our analysis, we have taken into account the recently discovered common red noise process, which has an impact at low frequencies. We also find that the peculiar noise features present in some pulsars data must be taken into account to reduce the false alarm. We show that using custom noise models is essential in searching for continuous gravitational wave signals and setting the upper limit.
△ Less
Submitted 19 March, 2023;
originally announced March 2023.
-
The MSPSR$π$ catalogue: VLBA astrometry of 18 millisecond pulsars
Authors:
H. Ding,
A. T. Deller,
B. W. Stappers,
T. J. W. Lazio,
D. Kaplan,
S. Chatterjee,
W. Brisken,
J. Cordes,
P. C. C. Freire,
E. Fonseca,
I. Stairs,
L. Guillemot,
A. Lyne,
I. Cognard,
D. J. Reardon,
G. Theureau
Abstract:
With unparalleled rotational stability, millisecond pulsars (MSPs) serve as ideal laboratories for numerous astrophysical studies, many of which require precise knowledge of the distance and/or velocity of the MSP. Here, we present the astrometric results for 18 MSPs of the "MSPSR$π$" project focusing exclusively on astrometry of MSPs, which includes the re-analysis of 3 previously published sourc…
▽ More
With unparalleled rotational stability, millisecond pulsars (MSPs) serve as ideal laboratories for numerous astrophysical studies, many of which require precise knowledge of the distance and/or velocity of the MSP. Here, we present the astrometric results for 18 MSPs of the "MSPSR$π$" project focusing exclusively on astrometry of MSPs, which includes the re-analysis of 3 previously published sources. On top of a standardized data reduction protocol, more complex strategies (i.e., normal and inverse-referenced 1D interpolation) were employed where possible to further improve astrometric precision. We derived astrometric parameters using sterne, a new Bayesian astrometry inference package that allows the incorporation of prior information based on pulsar timing where applicable. We measured significant ($>3\,σ$) parallax-based distances for 15 MSPs, including $0.81\pm0.02\,$kpc for PSR J1518+4904 -- the most significant model-independent distance ever measured for a double neutron star system. For each MSP with a well-constrained distance, we estimated its transverse space velocity and radial acceleration. Among the estimated radial accelerations, the updated ones of PSR J1012+5307 and PSR J1738+0333 impose new constraints on dipole gravitational radiation and the time derivative of Newton's gravitational constant. Additionally, significant angular broadening was detected for PSR J1643-1224, which offers an independent check of the postulated association between the HII region Sh 2-27 and the main scattering screen of PSR J1643-1224. Finally, the upper limit of the death line of $γ$-ray-emitting pulsars is refined with the new radial acceleration of the hitherto least energetic $γ$-ray pulsar PSR J1730-2304.
△ Less
Submitted 31 March, 2023; v1 submitted 12 December, 2022;
originally announced December 2022.
-
The MeerKAT Pulsar Timing Array: First Data Release
Authors:
Matthew T. Miles,
Ryan M. Shannon,
Matthew Bailes,
Daniel J. Reardon,
Michael J. Keith,
Andrew D. Cameron,
Aditya Parthasarathy,
Mohsen Shamohammadi,
Renee Spiewak,
Willem van Straten,
Sarah Buchner,
Fernando Camilo,
Marisa Geyer,
Aris Karastergiou,
Michael Kramer,
Maciej Serylak,
Gilles Theureau,
Vivek Venkatraman Krishnan
Abstract:
We present the first 2.5 years of data from the MeerKAT Pulsar Timing Array (MPTA), part of MeerTime, a MeerKAT Large Survey Project. The MPTA aims to precisely measure pulse arrival times from an ensemble of 88 pulsars visible from the Southern Hemisphere, with the goal of contributing to the search, detection and study of nanohertz-frequency gravitational waves as part of the International Pulsa…
▽ More
We present the first 2.5 years of data from the MeerKAT Pulsar Timing Array (MPTA), part of MeerTime, a MeerKAT Large Survey Project. The MPTA aims to precisely measure pulse arrival times from an ensemble of 88 pulsars visible from the Southern Hemisphere, with the goal of contributing to the search, detection and study of nanohertz-frequency gravitational waves as part of the International Pulsar Timing Array. This project makes use of the MeerKAT telescope, and operates with a typical observing cadence of two weeks using the L-band receiver that records data from 856-1712 MHz. We provide a comprehensive description of the observing system, software, and pipelines used and developed for the MeerTime project. The data products made available as part of this data release are from the 78 pulsars that had at least $30$ observations between the start of the MeerTime programme in February 2019 and October 2021. These include both sub-banded and band-averaged arrival times, as well as the initial timing ephemerides, noise models, and the frequency-dependent standard templates (portraits) used to derive pulse arrival times. After accounting for detected noise processes in the data, the frequency-averaged residuals of $67$ of the pulsars achieved a root-mean-square residual precision of $< 1 μ\rm{s}$. We also present a novel recovery of the clock correction waveform solely from pulsar timing residuals, and an exploration into preliminary findings of interest to the international pulsar timing community. The arrival times, standards and full Stokes parameter calibrated pulsar timing archives are publicly available.
△ Less
Submitted 8 December, 2022;
originally announced December 2022.
-
Searches for Shapiro delay in seven binary pulsars using the MeerKAT telescope
Authors:
Mohsen Shamohammadi,
Matthew Bailes,
Paulo C. C. Freire,
Aditya Parthasarathy,
Daniel J. Reardon,
Ryan M. Shannon,
Vivek Venkatraman Krishnan,
Miquel C. i. Bernadich,
Andrew D. Cameron,
David J. Champion,
Alessandro Corongiu,
Christopher Flynn,
Marisa Geyer,
Michael Kramer,
Matthew T. Miles,
Andrea Possenti,
Renee Spiewak
Abstract:
Precision timing of millisecond pulsars in binary systems enables observers to detect the relativistic Shapiro delay induced by space time curvature. When favourably aligned, this enables constraints to be placed on the component masses and system orientation. Here we present the results of timing campaigns on seven binary millisecond pulsars observed with the 64-antenna MeerKAT radio telescope th…
▽ More
Precision timing of millisecond pulsars in binary systems enables observers to detect the relativistic Shapiro delay induced by space time curvature. When favourably aligned, this enables constraints to be placed on the component masses and system orientation. Here we present the results of timing campaigns on seven binary millisecond pulsars observed with the 64-antenna MeerKAT radio telescope that show evidence of Shapiro delay: PSRs~J0101$-$6422, J1101$-$6424, J1125$-$6014, J1514$-$4946, J1614$-$2230, J1732$-$5049, and J1909$-$3744. Evidence for Shapiro delay was found in all of the systems, and for three the orientations and data quality enabled strong constraints on their orbital inclinations and component masses. For PSRs~J1125$-$6014, J1614$-$2230 and J1909$-$3744, we determined pulsar masses to be $M_{\rm p} = 1.68\pm 0.17 \, {\rm M_{\odot}} $, $1.94\pm 0.03 \, {\rm M_{\odot}} $ and $1.45 \pm 0.03 \, {\rm M_{\odot}}$, and companion masses to be $M_{\rm c} = 0.33\pm 0.02 \, {\rm M_{\odot}} $, $0.495\pm 0.005 \, {\rm M_{\odot}} $ and $0.205 \pm 0.003 \, {\rm M_{\odot}}$, respectively. This provides the first independent confirmation of PSR~J1614$-$2230's mass, one of the highest known. The Shapiro delays measured for PSRs~J0101$-$6422, J1101$-$6424, J1514$-$4946, and J1732$-$5049 were only weak, and could not provide interesting component mass limits. Despite a large number of millisecond pulsars being routinely timed, relatively few have accurate masses via Shapiro delays. We use simulations to show that this is expected, and provide a formula for observers to assess how accurately a pulsar mass can be determined. We also discuss the observed correlation between pulsar companion masses and spin period, and the anti-correlation between recycled pulsar mass and their companion masses.
△ Less
Submitted 7 December, 2022;
originally announced December 2022.
-
Wideband timing of the Parkes Pulsar Timing Array UWL data
Authors:
Małgorzata Curyło,
Timothy T. Pennucci,
Matthew Bailes,
N. D. Ramesh Bhat,
Andrew D. Cameron,
Shi Dai,
George Hobbs,
Agastya Kapur,
Richard N. Manchester,
Rami Mandow,
Matthew T. Miles,
Christopher J. Russell,
Daniel J. Reardon,
Ryan M. Shannon,
Renée Spiewak,
Andrew Zic,
Xing-Jiang Zhu
Abstract:
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 compa…
▽ More
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 $μ$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 $μ$s.
△ Less
Submitted 23 November, 2022;
originally announced November 2022.
-
The Thousand-Pulsar-Array programme on MeerKAT -- X. Scintillation arcs of 107 pulsars
Authors:
R. A. Main,
A. Parthasarathy,
S. Johnston,
A. Karastergiou,
A. Basu,
A. D. Cameron,
M. J. Keith,
L. S. Oswald,
B. Posselt,
D. J. Reardon,
X. Song,
P. Weltevrede
Abstract:
We present the detection of 107 pulsars with interstellar scintillation arcs at 856--1712\,MHz, observed with the MeerKAT Thousand Pulsar Array Programme. Scintillation arcs appear to be ubiquitous in clean, high S/N observations, their detection mainly limited by short observing durations and coarse frequency channel resolution. This led the survey to be sensitive to nearby, lightly scattered pul…
▽ More
We present the detection of 107 pulsars with interstellar scintillation arcs at 856--1712\,MHz, observed with the MeerKAT Thousand Pulsar Array Programme. Scintillation arcs appear to be ubiquitous in clean, high S/N observations, their detection mainly limited by short observing durations and coarse frequency channel resolution. This led the survey to be sensitive to nearby, lightly scattered pulsars with high effective velocity -- from a large proper motion, a screen nearby the pulsar, or a screen near the Earth. We measure the arc curvatures in all of our sources, which can be used to give an estimate of screen distances in pulsars with known proper motion, or an estimate of the proper motion. The short scintillation timescale in J1731$-$4744 implies a scattering screen within 12\,pc of the source, strongly suggesting the association between this pulsar and the supernova remnant RCW 114. We measure multiple parabolic arcs of 5 pulsars, all of which are weakly scintillating with high proper motion. Additionally, several sources show hints of inverted arclets suggesting scattering from anisotropic screens. Building on this work, further targeted MeerKAT observations of many of these pulsars will improve understanding of our local scattering environment and the origins of scintillation; annual scintillation curves would lead to robust screen distance measurements, and the evolution of arclets in time and frequency can constrain models of scintillation.
△ Less
Submitted 15 November, 2022;
originally announced November 2022.
-
Consistency of the Parkes Pulsar Timing Array Signal with a Nanohertz Gravitational-wave Background
Authors:
Boris Goncharov,
Eric Thrane,
Ryan M. Shannon,
Jan Harms,
N. D. Ramesh Bhat,
George Hobbs,
Matthew Kerr,
Richard N. Manchester,
Daniel J. Reardon,
Christopher J. Russell,
Xing-Jiang Zhu,
Andrew Zic
Abstract:
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)…
▽ More
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.
△ Less
Submitted 21 June, 2022; v1 submitted 8 June, 2022;
originally announced June 2022.
-
Orbital dynamics and extreme scattering event properties from long-term scintillation observations of PSR J1603-7202
Authors:
Kris Walker,
Daniel J. Reardon,
Eric Thrane,
Rory Smith
Abstract:
We model long-term variations in the scintillation of binary pulsar PSR J1603$-$7202, observed by the 64 m Parkes radio telescope (Murriyang) between 2004 and 2016. We find that the time variation in the scintillation arc curvature is well-modelled by scattering from an anisotropic thin screen of plasma between the Earth and the pulsar. Using our scintillation model, we measure the inclination ang…
▽ More
We model long-term variations in the scintillation of binary pulsar PSR J1603$-$7202, observed by the 64 m Parkes radio telescope (Murriyang) between 2004 and 2016. We find that the time variation in the scintillation arc curvature is well-modelled by scattering from an anisotropic thin screen of plasma between the Earth and the pulsar. Using our scintillation model, we measure the inclination angle and longitude of ascending node of the orbit, yielding a significant improvement over the constraints from pulsar timing. From our measurement of the inclination angle, we place a lower bound on the mass of J1603$-$7202's companion of $\gtrsim 0.5\,\text{M}_\odot$ assuming a pulsar mass of $\gtrsim1.2\,\text{M}_\odot$. We find that the scintillation arcs are most pronounced when the electron column density along the line of sight is increased, and that arcs are present during a known extreme scattering event. We measure the distance to the interstellar plasma and its velocity, and we discuss some structures seen in individual scintillation arcs within the context of our model.
△ Less
Submitted 23 April, 2022;
originally announced April 2022.
-
The MeerTime Pulsar Timing Array -- A Census of Emission Properties and Timing Potential
Authors:
R. Spiewak,
M. Bailes,
M. T. Miles,
A. Parthasarathy,
D. J. Reardon,
M. Shamohammadi,
R. M. Shannon,
N. D. R. Bhat,
S. Buchner,
A. D. Cameron,
F. Camilo,
M. Geyer,
S. Johnston,
A. Karastergiou,
M. Keith,
M. Kramer,
M. Serylak,
W. van Straten,
G. Theureau,
V. Venkatraman Krishnan
Abstract:
MeerTime is a five-year Large Survey Project to time pulsars with MeerKAT, the 64-dish South African precursor to the Square Kilometre Array. The science goals for the programme include timing millisecond pulsars (MSPs) to high precision (< 1 $μ$s) to study the Galactic MSP population and to contribute to global efforts to detect nanohertz gravitational waves with the International Pulsar Timing A…
▽ More
MeerTime is a five-year Large Survey Project to time pulsars with MeerKAT, the 64-dish South African precursor to the Square Kilometre Array. The science goals for the programme include timing millisecond pulsars (MSPs) to high precision (< 1 $μ$s) to study the Galactic MSP population and to contribute to global efforts to detect nanohertz gravitational waves with the International Pulsar Timing Array (IPTA). In order to plan for the remainder of the programme and to use the allocated time most efficiently, we have conducted an initial census with the MeerKAT "L-band" receiver of 189 MSPs visible to MeerKAT and here present their dispersion measures, polarization profiles, polarization fractions, rotation measures, flux density measurements, spectral indices, and timing potential. As all of these observations are taken with the same instrument (which uses coherent dedispersion, interferometric polarization calibration techniques, and a uniform flux scale), they present an excellent resource for population studies. We used wideband pulse portraits as timing standards for each MSP and demonstrated that the MeerTime Pulsar Timing Array (MPTA) can already contribute significantly to the IPTA as it currently achieves better than 1 $μ$s timing accuracy on 89 MSPs (observed with fortnightly cadence). By the conclusion of the initial five-year MeerTime programme in July 2024, the MPTA will be extremely significant in global efforts to detect the gravitational wave background with a contribution to the detection statistic comparable to other long-standing timing programmes.
△ Less
Submitted 8 April, 2022;
originally announced April 2022.
-
The eccentric millisecond pulsar, PSR J0955$-$6150 I: Pulse profile analysis, mass measurements and constraints on binary evolution
Authors:
M. Serylak,
V. Venkatraman Krishnan,
P. C. C. Freire,
T. M. Tauris,
M. Kramer,
M. Geyer,
A. Parthasarathy,
M. Bailes,
M. C. i Bernadich,
S. Buchner,
M. Burgay,
F. Camilo,
A. Karastergiou,
M. E. Lower,
A. Possenti,
D. J. Reardon,
R. M. Shannon,
R. Spiewak,
I. H. Stairs,
W. van Straten
Abstract:
PSR J0955$-$6150 is a member of a class of eccentric MSP+He WD systems (eMSPs), whose binary evolution is poorly understood and believed to be different to that of traditional MSP+He WD systems. Measuring the masses of the stars in this system is important for testing hypotheses for the formation of eMSPs. We have carried out observations of this pulsar with the Parkes and MeerKAT radio telescopes…
▽ More
PSR J0955$-$6150 is a member of a class of eccentric MSP+He WD systems (eMSPs), whose binary evolution is poorly understood and believed to be different to that of traditional MSP+He WD systems. Measuring the masses of the stars in this system is important for testing hypotheses for the formation of eMSPs. We have carried out observations of this pulsar with the Parkes and MeerKAT radio telescopes. Our observations reveal a strong frequency evolution of this pulsar's intensity, with a spectral index ($α$) of $-3.13(2)$. The sensitivity of MeerKAT has resulted in a $>10$-fold improvement in the timing precision compared to older Parkes observations. Combined with the 8-year timing baseline, it has allowed precise measurements of a proper motion and three orbital "post-Keplerian" parameters: the rate of advance of periastron, $\dotω = 0.00152(1) \, {\rm deg} \, yr^{-1}$ and the orthometric Shapiro delay parameters, $h_3 = 0.89(7) \, μ$s and $ς= 0.88(2)$. Assuming general relativity, we obtain $M_{p} = 1.71(2) \, M_{\odot}$ for the mass of the pulsar and $M_{c} = 0.254(2) \, M_{\odot}$ for the mass of the companion; the orbital inclination is 83.2(4) degrees. We find that the spin axis has a misalignment relative to the orbital angular momentum of $> 4.8$ degrees at 99% CI. While the value of $M_{\rm p}$ falls within the wide range observed in eMSPs, $M_{\rm c}$ is significantly smaller than expected, allowing several formation hypotheses being ruled out. $M_{\rm c}$ is also significantly different from the expected value for an ideal low mass X-ray binary evolution scenario. The putative misalignment between the spin axis of the pulsar and the orbital angular momentum suggests that the unknown process that created the orbital eccentricity of the binary was also capable of changing its orbital orientation, an important evidence for understanding the origin of eMSPs.
△ Less
Submitted 1 March, 2022;
originally announced March 2022.
-
Modelling annual scintillation arc variations in PSR J1643-1224 using the Large European Array for Pulsars
Authors:
G. Mall,
R. A. Main,
J. Antoniadis,
C. G. Bassa,
M. Burgay,
S. Chen,
I. Cognard,
R. Concu,
A. Corongiu,
M. Gaikwad,
H. Hu,
G. H. Janssen,
R. Karuppusamy,
K. J. Lee,
K. Liu,
J. W. McKee,
A. Melis,
M. B. Mickaliger,
D. Perrodin,
M. Pilia,
A. Possenti,
D. J. Reardon,
S. A. Sanidas,
T. Sprenger,
B. W. Stappers
, et al. (3 additional authors not shown)
Abstract:
In this work we study variations in the parabolic scintillation arcs of the binary millisecond pulsar PSR J1643-1224 over five years using the Large European Array for Pulsars (LEAP). The 2D power spectrum of scintillation, called the secondary spectrum, often shows a parabolic distribution of power, where the arc curvature encodes the relative velocities and distances of the pulsar, ionised inter…
▽ More
In this work we study variations in the parabolic scintillation arcs of the binary millisecond pulsar PSR J1643-1224 over five years using the Large European Array for Pulsars (LEAP). The 2D power spectrum of scintillation, called the secondary spectrum, often shows a parabolic distribution of power, where the arc curvature encodes the relative velocities and distances of the pulsar, ionised interstellar medium (IISM), and Earth. We observe a clear parabolic scintillation arc which varies in curvature throughout the year. The distribution of power in the secondary spectra are inconsistent with a single scattering screen which is fully 1D, or entirely isotropic. We fit the observed arc curvature variations with two models; an isotropic scattering screen, and a model with two independent 1D screens. We measure the distance to the scattering screen to be in the range 114-223 pc, depending on the model, consistent with the known distance of the foreground large-diameter HII region Sh 2-27 (112+/-17 pc), suggesting that it is the dominant source of scattering. We obtain only weak constraints on the pulsar's orbital inclination and angle of periastron, since the scintillation pattern is not very sensitive to the pulsar's motion, since the screen is much closer to the Earth than the pulsar. More measurements of this kind - where scattering screens can be associated with foreground objects - will help to inform the origins and distribution of scattering screens within our galaxy.
△ Less
Submitted 11 January, 2022;
originally announced January 2022.
-
The International Pulsar Timing Array second data release: Search for an isotropic Gravitational Wave Background
Authors:
J. Antoniadis,
Z. Arzoumanian,
S. Babak,
M. Bailes,
A. -S. Bak Nielsen,
P. T. Baker,
C. G. Bassa,
B. Becsy,
A. Berthereau,
M. Bonetti,
A. Brazier,
P. R. Brook,
M. Burgay,
S. Burke-Spolaor,
R. N. Caballero,
J. A. Casey-Clyde,
A. Chalumeau,
D. J. Champion,
M. Charisi,
S. Chatterjee,
S. Chen,
I. Cognard,
J. M. Cordes,
N. J. Cornish,
F. Crawford
, et al. (101 additional authors not shown)
Abstract:
We searched for an isotropic stochastic gravitational wave background in the second data release of the International Pulsar Timing Array, a global collaboration synthesizing decadal-length pulsar-timing campaigns in North America, Europe, and Australia. In our reference search for a power law strain spectrum of the form $h_c = A(f/1\,\mathrm{yr}^{-1})^α$, we found strong evidence for a spectrally…
▽ More
We searched for an isotropic stochastic gravitational wave background in the second data release of the International Pulsar Timing Array, a global collaboration synthesizing decadal-length pulsar-timing campaigns in North America, Europe, and Australia. In our reference search for a power law strain spectrum of the form $h_c = A(f/1\,\mathrm{yr}^{-1})^α$, we found strong evidence for a spectrally-similar low-frequency stochastic process of amplitude $A = 3.8^{+6.3}_{-2.5}\times10^{-15}$ and spectral index $α= -0.5 \pm 0.5$, where the uncertainties represent 95\% credible regions, using information from the auto- and cross-correlation terms between the pulsars in the array. For a spectral index of $α= -2/3$, as expected from a population of inspiralling supermassive black hole binaries, the recovered amplitude is $A = 2.8^{+1.2}_{-0.8}\times10^{-15}$. Nonetheless, no significant evidence of the Hellings-Downs correlations that would indicate a gravitational-wave origin was found. We also analyzed the constituent data from the individual pulsar timing arrays in a consistent way, and clearly demonstrate that the combined international data set is more sensitive. Furthermore, we demonstrate that this combined data set produces comparable constraints to recent single-array data sets which have more data than the constituent parts of the combination. Future international data releases will deliver increased sensitivity to gravitational wave radiation, and significantly increase the detection probability.
△ Less
Submitted 11 January, 2022;
originally announced January 2022.
-
High-precision search for dark photon dark matter with the Parkes Pulsar Timing Array
Authors:
Xiao Xue,
Zi-Qing Xia,
Xingjiang Zhu,
Yue Zhao,
Jing Shu,
Qiang Yuan,
N. D. Ramesh Bhat,
Andrew D. Cameron,
Shi Dai,
Yi Feng,
Boris Goncharov,
George Hobbs,
Eric Howard,
Richard N. Manchester,
Aditya Parthasarathy,
Daniel J. Reardon,
Christopher J. Russell,
Ryan M. Shannon,
Renée Spiewak,
Nithyanandan Thyagarajan,
Jingbo Wang,
Lei Zhang,
Songbo Zhang
Abstract:
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…
▽ More
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 $ε$ 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.
△ Less
Submitted 26 December, 2021; v1 submitted 14 December, 2021;
originally announced December 2021.
-
Mode changing in J1909$-$3744: the most precisely timed pulsar
Authors:
Matthew T. Miles,
Ryan M. Shannon,
Matthew Bailes,
Daniel J. Reardon,
Sarah Buchner,
Hannah Middleton,
Renee Spiewak
Abstract:
We present baseband radio observations of the millisecond pulsar J1909$-$3744, the most precisely timed pulsar, using the MeerKAT telescope as part of the MeerTime pulsar timing array campaign. During a particularly bright scintillation event the pulsar showed strong evidence of pulse mode changing, among the first millisecond pulsars and the shortest duty cycle millisecond pulsar to do so. Two mo…
▽ More
We present baseband radio observations of the millisecond pulsar J1909$-$3744, the most precisely timed pulsar, using the MeerKAT telescope as part of the MeerTime pulsar timing array campaign. During a particularly bright scintillation event the pulsar showed strong evidence of pulse mode changing, among the first millisecond pulsars and the shortest duty cycle millisecond pulsar to do so. Two modes appear to be present, with the weak (lower signal-to-noise ratio) mode arriving $9.26$ $\pm 3.94$ $μ$s earlier than the strong counterpart. Further, we present a new value of the jitter noise for this pulsar of $8.20 \pm 0.14$ ns in one hour, finding it to be consistent with previous measurements taken with the MeerKAT ($9 \pm 3$ ns) and Parkes ($8.6 \pm 0.8$ ns) telescopes, but inconsistent with the previously most precise measurement taken with the Green Bank telescope ($14 \pm 0.5$ ns). Timing analysis on the individual modes is carried out for this pulsar, and we find an approximate $10\%$ improvement in the timing precision is achievable through timing the strong mode only as opposed to the full sample of pulses. By forming a model of the average pulse from templates of the two modes, we time them simultaneously and demonstrate that this timing improvement can also be achieved in regular timing observations. We discuss the impact an improvement of this degree on this pulsar would have on searches for the stochastic gravitational wave background, as well as the impact of a similar improvement on all MeerTime PTA pulsars.
△ Less
Submitted 1 December, 2021;
originally announced December 2021.
-
Constraining cosmological phase transitions with the Parkes Pulsar Timing Array
Authors:
Xiao Xue,
Ligong Bian,
Jing Shu,
Qiang Yuan,
Xingjiang Zhu,
N. D. Ramesh Bhat,
Shi Dai,
Yi Feng,
Boris Goncharov,
George Hobbs,
Eric Howard,
Richard N. Manchester,
Christopher J. Russell,
Daniel J. Reardon,
Ryan M. Shannon,
Renée Spiewak,
Nithyanandan Thyagarajan,
Jingbo Wang
Abstract:
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 backgroun…
▽ More
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 $(β/H_*)^{-1}\sim 10^{-2}-10^{-1}$ and therefore can be used to constrain the dark and QCD phase transitions.
△ Less
Submitted 25 November, 2022; v1 submitted 6 October, 2021;
originally announced October 2021.
-
The Thousand-Pulsar-Array programme on MeerKAT: -- VI. Pulse widths of a large and diverse sample of radio pulsars
Authors:
B. Posselt,
A. Karastergiou,
S. Johnston,
A. Parthasarathy,
M. J. Keith,
L. S. Oswald,
X. Song,
P. Weltevrede,
E. D. Barr,
S. Buchner,
M. Geyer,
M. Kramer,
D. J. Reardon,
M. Serylak,
R. M. Shannon,
R. Spiewak,
V. Venkatraman Krishnan
Abstract:
We present pulse width measurements for a sample of radio pulsars observed with the MeerKAT telescope as part of the Thousand-Pulsar-Array (TPA) programme in the MeerTime project. For a centre frequency of 1284 MHz, we obtain 762 $W_{10}$ measurements across the total bandwidth of 775 MHz, where $W_{10}$ is the width at the 10% level of the pulse peak. We also measure about 400 $W_{10}$ values in…
▽ More
We present pulse width measurements for a sample of radio pulsars observed with the MeerKAT telescope as part of the Thousand-Pulsar-Array (TPA) programme in the MeerTime project. For a centre frequency of 1284 MHz, we obtain 762 $W_{10}$ measurements across the total bandwidth of 775 MHz, where $W_{10}$ is the width at the 10% level of the pulse peak. We also measure about 400 $W_{10}$ values in each of the four or eight frequency sub-bands. Assuming, the width is a function of the rotation period P, this relationship can be described with a power law with power law index $μ=-0.29\pm 0.03$. However, using orthogonal distance regression, we determine a steeper power law with $μ=-0.63\pm 0.06$. A density plot of the period-width data reveals such a fit to align well with the contours of highest density. Building on a previous population synthesis model, we obtain population-based estimates of the obliquity of the magnetic axis with respect to the rotation axis for our pulsars. Investigating the width changes over frequency, we unambiguously identify a group of pulsars that have width broadening at higher frequencies. The measured width changes show a monotonic behaviour with frequency for the whole TPA pulsar population, whether the pulses are becoming narrower or broader with increasing frequency. We exclude a sensitivity bias, scattering and noticeable differences in the pulse component numbers as explanations for these width changes, and attempt an explanation using a qualitative model of five contributing Gaussian pulse components with flux density spectra that depend on their rotational phase.
△ Less
Submitted 24 September, 2021;
originally announced September 2021.
-
On the evidence for a common-spectrum process in the search for the nanohertz gravitational-wave background with the Parkes Pulsar Timing Array
Authors:
Boris Goncharov,
R. M. Shannon,
D. J. Reardon,
G. Hobbs,
A. Zic,
M. Bailes,
M. Curylo,
S. Dai,
M. Kerr,
M. E. Lower,
R. N. Manchester,
R. Mandow,
H. Middleton,
M. T. Miles,
A. Parthasarathy,
E. Thrane,
N. Thyagarajan,
X. Xue,
X. J. Zhu,
A. D. Cameron,
Y. Feng,
R. Luo,
C. J. Russell,
J. Sarkissian,
R. Spiewak
, et al. (4 additional authors not shown)
Abstract:
A nanohertz-frequency stochastic gravitational-wave background can potentially be detected through the precise timing of an array of millisecond pulsars. This background produces low-frequency noise in the pulse arrival times that would have a characteristic spectrum common to all pulsars and a well-defined spatial correlation. Recently the North American Nanohertz Observatory for Gravitational Wa…
▽ More
A nanohertz-frequency stochastic gravitational-wave background can potentially be detected through the precise timing of an array of millisecond pulsars. This background produces low-frequency noise in the pulse arrival times that would have a characteristic spectrum common to all pulsars and a well-defined spatial correlation. Recently the North American Nanohertz Observatory for Gravitational Waves collaboration (NANOGrav) found evidence for the common-spectrum component in their 12.5-year data set. Here we report on a search for the background using the second data release of the Parkes Pulsar Timing Array. If we are forced to choose between the two NANOGrav models $\unicode{x2014}$ one with a common-spectrum process and one without $\unicode{x2014}$ we find strong support for the common-spectrum process. However, in this paper, we consider the possibility that the analysis suffers from model misspecification. In particular, we present simulated data sets that contain noise with distinctive spectra but show strong evidence for a common-spectrum process under the standard assumptions. The Parkes data show no significant evidence for, or against, the spatially correlated Hellings-Downs signature of the gravitational-wave background. Assuming we did observe the process underlying the spatially uncorrelated component of the background, we infer its amplitude to be $A = 2.2^{+0.4}_{-0.3} \times 10^{-15}$ in units of gravitational-wave strain at a frequency of $1\, \text{yr}^{-1}$. Extensions and combinations of existing and new data sets will improve the prospects of identifying spatial correlations that are necessary to claim a detection of the gravitational-wave background.
△ Less
Submitted 11 August, 2021; v1 submitted 26 July, 2021;
originally announced July 2021.
-
The Parkes pulsar timing array second data release: Timing analysis
Authors:
D. J. Reardon,
R. M. Shannon,
A. D. Cameron,
B. Goncharov,
G. B. Hobbs,
H. Middleton,
M. Shamohammadi,
N. Thyagarajan,
M. Bailes,
N. D. R. Bhat,
S. Dai,
M. Kerr,
R. N. Manchester,
C. J. Russell,
R. Spiewak,
J. B. Wang,
X. J. Zhu
Abstract:
The main goal of pulsar timing array experiments is to detect correlated signals such as nanohertz-frequency gravitational waves. Pulsar timing data collected in dense monitoring campaigns can also be used to study the stars themselves, their binary companions, and the intervening ionised interstellar medium. Timing observations are extraordinarily sensitive to changes in path length between the p…
▽ More
The main goal of pulsar timing array experiments is to detect correlated signals such as nanohertz-frequency gravitational waves. Pulsar timing data collected in dense monitoring campaigns can also be used to study the stars themselves, their binary companions, and the intervening ionised interstellar medium. Timing observations are extraordinarily sensitive to changes in path length between the pulsar and the Earth, enabling precise measurements of the pulsar positions, distances and velocities, and the shapes of their orbits. Here we present a timing analysis of 25 pulsars observed as part of the Parkes Pulsar Timing Array (PPTA) project over time spans of up to 24 yr. The data are from the second data release of the PPTA, which we have extended by including legacy data. We make the first detection of Shapiro delay in four Southern pulsars (PSRs J1017$-$7156, J1125$-$6014, J1545$-$4550, and J1732$-$5049), and of parallax in six pulsars. The prominent Shapiro delay of PSR J1125$-$6014 implies a neutron star mass of $M_p = 1.5 \pm 0.2 M_\odot$ (68% credibility interval). Measurements of both Shapiro delay and relativistic periastron advance in PSR J1600$-$3053 yield a large but uncertain pulsar mass of $M_p = 2.06^{+0.44}_{-0.41}$ M$_\odot$ (68% credibility interval). We measure the distance to PSR J1909$-$3744 to a precision of 10 lyr, indicating that for gravitational wave periods over a decade, the pulsar provides a coherent baseline for pulsar timing array experiments.
△ Less
Submitted 16 July, 2021; v1 submitted 9 July, 2021;
originally announced July 2021.
-
The Thousand-Pulsar-Array programme on MeerKAT -- V. Scattering analysis of single-component pulsars
Authors:
L. S. Oswald,
A. Karastergiou,
B. Posselt,
S. Johnston,
M. Bailes,
S. Buchner,
M. Geyer,
M. J. Keith,
M. Kramer,
A. Parthasarathy,
D. J. Reardon,
M. Serylak,
R. M. Shannon,
R. Spiewak,
W. van Straten,
V. Venkatraman Krishnan
Abstract:
We have measured the scattering timescale, $τ$, and the scattering spectral index, $α$, for 84 single-component pulsars. Observations were carried out with the MeerKAT telescope as part of the Thousand-Pulsar-Array programme in the MeerTime project at frequencies between 0.895 and 1.670 GHz. Our results give a distribution of values for $α$ (defined in terms of $τ$ and frequency $ν$ as…
▽ More
We have measured the scattering timescale, $τ$, and the scattering spectral index, $α$, for 84 single-component pulsars. Observations were carried out with the MeerKAT telescope as part of the Thousand-Pulsar-Array programme in the MeerTime project at frequencies between 0.895 and 1.670 GHz. Our results give a distribution of values for $α$ (defined in terms of $τ$ and frequency $ν$ as $τ\proptoν^{-α}$) for which, upon fitting a Gaussian, we obtain a mean and standard deviation of $\langleα\rangle = 4.0 \pm 0.6$. This is due to our identification of possible causes of inaccurate measurement of $τ$, which, if not filtered out of modelling results, tend to lead to underestimation of $α$. The pulsars in our sample have large dispersion measures and are therefore likely to be distant. We find that a model using an isotropic scatter broadening function is consistent with the data, likely due to the averaging effect of multiple scattering screens along the line of sight. Our sample of scattering parameters provides a strong data set upon which we can build to test more complex and time-dependent scattering phenomena, such as extreme scattering events.
△ Less
Submitted 2 April, 2021;
originally announced April 2021.
-
Eight new millisecond pulsars from the first MeerKAT globular cluster census
Authors:
A. Ridolfi,
T. Gautam,
P. C. C. Freire,
S. M. Ransom,
S. J. Buchner,
A. Possenti,
V. Venkatraman Krishnan,
M. Bailes,
M. Kramer,
B. W. Stappers,
F. Abbate,
E. D. Barr,
M. Burgay,
F. Camilo,
A. Corongiu,
A. Jameson,
P. V. Padmanabh,
L. Vleeschower,
D. J. Champion,
M. Geyer,
A. Karastergiou,
R. Karuppusamy,
A. Parthasarathy,
D. J. Reardon,
M. Serylak
, et al. (2 additional authors not shown)
Abstract:
We have used the central 44 antennas of the new 64-dish MeerKAT radio telescope array to conduct a deep search for new pulsars in the core of nine globular clusters. This has led to the discovery of eight new millisecond pulsars in six different clusters. Two new binaries, 47 Tuc ac and 47 Tuc ad, are eclipsing "spiders", featuring compact orbits ($\lesssim 0.32$ days), very low-mass companions an…
▽ More
We have used the central 44 antennas of the new 64-dish MeerKAT radio telescope array to conduct a deep search for new pulsars in the core of nine globular clusters. This has led to the discovery of eight new millisecond pulsars in six different clusters. Two new binaries, 47 Tuc ac and 47 Tuc ad, are eclipsing "spiders", featuring compact orbits ($\lesssim 0.32$ days), very low-mass companions and regular occultations of their pulsed emission. The other three new binary pulsars (NGC 6624G, M62G, and Ter 5 an) are in wider ($> 0.7$ days) orbits, with companions that are likely to be white dwarfs or neutron stars. NGC 6624G has a large eccentricity of $e\simeq 0.38$, which enabled us to detect the rate of advance of periastron. This suggests that the system is massive, with a total mass of $M{\rm tot} = 2.65 \pm 0.07$ M$_{\odot}$. Likewise, for Ter 5 an, with $e \simeq 0.0066$, we obtain $M{\rm tot}= 2.97 \pm 0.52$ M$_{\odot}$. The other three new discoveries (NGC 6522D, NGC 6624H and NGC 6752F) are faint isolated pulsars. Finally, we have used the whole MeerKAT array and synthesized 288 beams, covering an area of $\sim2$ arcmin in radius around the center of NGC 6624. This has allowed us to localize many of the pulsars in the cluster, demonstrating the beamforming capabilities of the TRAPUM software backend and paving the way for the upcoming MeerKAT globular cluster pulsar survey.
△ Less
Submitted 8 March, 2021;
originally announced March 2021.
-
The Relativistic Binary Programme on MeerKAT: Science objectives and first results
Authors:
M. Kramer,
I. H. Stairs,
V. Venkatraman Krishnan,
P. C. C. Freire,
F. Abbate,
M. Bailes,
M. Burgay,
S. Buchner,
D. J. Champion,
I. Cognard,
T. Gautam,
M. Geyer,
L. Guillemot,
H. Hu,
G. Janssen,
M. E. Lower,
A. Parthasarathy,
A. Possenti,
S. Ransom,
D. J. Reardon,
A. Ridolfi,
M. Serylak,
R. M. Shannon,
R. Spiewak,
G. Theureau
, et al. (13 additional authors not shown)
Abstract:
We describe the ongoing Relativistic Binary programme (RelBin), a part of the MeerTime large survey project with the MeerKAT radio telescope. RelBin is primarily focused on observations of relativistic effects in binary pulsars to enable measurements of neutron star masses and tests of theories of gravity. We selected 25 pulsars as an initial high priority list of targets based on their characteri…
▽ More
We describe the ongoing Relativistic Binary programme (RelBin), a part of the MeerTime large survey project with the MeerKAT radio telescope. RelBin is primarily focused on observations of relativistic effects in binary pulsars to enable measurements of neutron star masses and tests of theories of gravity. We selected 25 pulsars as an initial high priority list of targets based on their characteristics and observational history with other telescopes. In this paper, we provide an outline of the programme, present polarisation calibrated pulse profiles for all selected pulsars as a reference catalogue along with updated dispersion measures. We report Faraday rotation measures for 24 pulsars, twelve of which have been measured for the first time. More than a third of our selected pulsars show a flat position angle swing confirming earlier observations. We demonstrate the ability of the Rotating Vector Model (RVM), fitted here to seven binary pulsars, including the Double Pulsar (PSR J0737$-$3039A), to obtain information about the orbital inclination angle. We present a high time resolution light curve of the eclipse of PSR J0737$-$3039A by the companion's magnetosphere, a high-phase resolution position angle swing for PSR J1141$-$6545, an improved detection of the Shapiro delay of PSR J1811$-$2405, and pulse scattering measurements for PSRs J1227$-$6208, J1757$-$1854, and J1811$-$1736. Finally, we demonstrate that timing observations with MeerKAT improve on existing data sets by a factor of, typically, 2-3, sometimes by an order of magnitude.
△ Less
Submitted 7 May, 2021; v1 submitted 9 February, 2021;
originally announced February 2021.
-
Measurements of pulse jitter and single-pulse variability in millisecond pulsars using MeerKAT
Authors:
A. Parthasarathy,
M. Bailes,
R. M. Shannon,
W. van Straten,
S. Oslowski,
S. Johnston,
R. Spiewak,
D. J. Reardon,
M. Kramer,
V. Venkatraman Krishnan,
T. T. Pennucici,
F. Abbate,
S. Buchner,
F. Camilo,
D. J. Champion,
M. Geyer,
B. Hugo,
A. Jameson,
A. Karastergiou,
M. J. Keith,
M. Serylak
Abstract:
Using the state-of-the-art SKA precursor, the MeerKAT radio telescope, we explore the limits to precision pulsar timing of millisecond pulsars achievable due to pulse stochasticity (jitter). We report new jitter measurements in 15 of the 29 pulsars in our sample and find that the levels of jitter can vary dramatically between them. For some, like the 2.2~ms pulsar PSR J2241--5236, we measure an im…
▽ More
Using the state-of-the-art SKA precursor, the MeerKAT radio telescope, we explore the limits to precision pulsar timing of millisecond pulsars achievable due to pulse stochasticity (jitter). We report new jitter measurements in 15 of the 29 pulsars in our sample and find that the levels of jitter can vary dramatically between them. For some, like the 2.2~ms pulsar PSR J2241--5236, we measure an implied jitter of just $\sim$ 4~ns/hr, while others like the 3.9~ms PSR J0636--3044 are limited to $\sim$ 100 ns/hr. While it is well known that jitter plays a central role to limiting the precision measurements of arrival times for high signal-to-noise ratio observations, its role in the measurement of dispersion measure (DM) has not been reported, particularly in broad-band observations. Using the exceptional sensitivity of MeerKAT, we explored this on the bright millisecond pulsar PSR J0437--4715 by exploring the DM of literally every pulse. We found that the derived single pulse DMs vary by typically 0.0085 cm$^{-3}$ pc from the mean, and that the best DM estimate is limited by the differential pulse jitter across the band. We postulate that all millisecond pulsars will have their own limit on DM precision which can only be overcome with longer integrations. Using high-time resolution filterbank data of 9 $μ$s, we also present a statistical analysis of single pulse phenomenology. Finally, we discuss optimization strategies for the MeerKAT pulsar timing program and its role in the context of the International Pulsar Timing Array (IPTA).
△ Less
Submitted 21 January, 2021;
originally announced January 2021.
-
The Thousand-Pulsar-Array programme on MeerKAT II: observing strategy for pulsar monitoring with subarrays
Authors:
X. Song,
P. Weltevrede,
M. J. Keith,
S. Johnston,
A. Karastergiou,
M. Bailes,
E. D. Barr,
S. Buchner,
M. Geyer,
B. V. Hugo,
A. Jameson,
A. Parthasarathy,
D. J. Reardon,
M. Serylak,
R. M. Shannon,
R. Spiewak,
W. van Straten,
V. Venkatraman Krishnan
Abstract:
The Thousand Pulsar Array (TPA) project currently monitors about 500 pulsars with the sensitive MeerKAT radio telescope by using subarrays to observe multiple sources simultaneously. Here we define the adopted observing strategy, which guarantees that each target is observed long enough to obtain a high fidelity pulse profile, thereby reaching a sufficient precision of a simple pulse shape paramet…
▽ More
The Thousand Pulsar Array (TPA) project currently monitors about 500 pulsars with the sensitive MeerKAT radio telescope by using subarrays to observe multiple sources simultaneously. Here we define the adopted observing strategy, which guarantees that each target is observed long enough to obtain a high fidelity pulse profile, thereby reaching a sufficient precision of a simple pulse shape parameter. This precision is estimated from the contribution of the system noise of the telescope, and the pulse-to-pulse variability of each pulsar, which we quantify under some simplifying assumptions. We test the assumptions and choice of model parameters using data from the MeerKAT 64-dish array, Lovell and Parkes telescopes. We demonstrate that the observing times derived from our method produce high fidelity pulse profiles that meet the needs of the TPA in studying pulse shape variability and pulsar timing. Our method can also be used to compare strategies for observing large numbers of pulsars with telescopes capable of forming multiple subarray configurations. We find that using two 32-dish MeerKAT subarrays is the most efficient strategy for the TPA project. We also find that the ability to observe in different array configurations will become increasingly important for large observing programmes using the Square Kilometre Array telescope.
△ Less
Submitted 7 December, 2020;
originally announced December 2020.
-
Identifying and mitigating noise sources in precision pulsar timing data sets
Authors:
Boris Goncharov,
D. J. Reardon,
R. M. Shannon,
Xing-Jiang Zhu,
Eric Thrane,
M. Bailes,
N. D. R. Bhat,
S. Dai,
G. Hobbs,
M. Kerr,
R. N. Manchester,
S. Osłowski,
A. Parthasarathy,
C. J. Russell,
R. Spiewak,
N. Thyagarajan,
J. B. Wang
Abstract:
Pulsar timing array projects measure the pulse arrival times of millisecond pulsars for the primary purpose of detecting nanohertz-frequency gravitational waves. The measurements include contributions from a number of astrophysical and instrumental processes, which can either be deterministic or stochastic. It is necessary to develop robust statistical and physical models for these noise processes…
▽ More
Pulsar timing array projects measure the pulse arrival times of millisecond pulsars for the primary purpose of detecting nanohertz-frequency gravitational waves. The measurements include contributions from a number of astrophysical and instrumental processes, which can either be deterministic or stochastic. It is necessary to develop robust statistical and physical models for these noise processes because incorrect models diminish sensitivity and may cause a spurious gravitational wave detection. Here we characterise noise processes for the 26 pulsars in the second data release of the Parkes Pulsar Timing Array using Bayesian inference. In addition to well-studied noise sources found previously in pulsar timing array data sets such as achromatic timing noise and dispersion measure variations, we identify new noise sources including time-correlated chromatic noise that we attribute to variations in pulse scattering. We also identify "exponential dip" events in four pulsars, which we attribute to magnetospheric effects as evidenced by pulse profile shape changes observed for three of the pulsars. This includes an event in PSR J1713$+$0747, which had previously been attributed to interstellar propagation. We present noise models to be used in searches for gravitational waves. We outline a robust methodology to evaluate the performance of noise models and identify unknown signals in the data. The detection of variations in pulse profiles highlights the need to develop efficient profile domain timing methods.
△ Less
Submitted 1 November, 2020; v1 submitted 12 October, 2020;
originally announced October 2020.
-
Precision orbital dynamics from interstellar scintillation arcs for PSR J0437-4715
Authors:
Daniel J. Reardon,
William A. Coles,
Matthew Bailes,
N. D. Ramesh Bhat,
Shi Dai,
George B. Hobbs,
Matthew Kerr,
Richard N. Manchester,
Stefan Oslowski,
Aditya Parthasarathy,
Christopher J. Russell,
Ryan M. Shannon,
Renee Spiewak,
Lawrence Toomey,
Artem V. Tuntsov,
Willem van Straten,
Mark A. Walker,
Jingbo Wang,
Lei Zhang,
Xing-Jiang Zhu
Abstract:
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 l…
▽ More
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, $Ω=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.
△ Less
Submitted 27 September, 2020;
originally announced September 2020.
-
The Thousand-Pulsar-Array programme on MeerKAT IV: Polarisation properties of young, energetic pulsars
Authors:
M. Serylak,
S. Johnston,
M. Kramer,
S. Buchner,
A. Karastergiou,
M. J. Keith,
A. Parthasarathy,
P. Weltevrede,
M. Bailes,
E. D. Barr,
F. Camilo,
M. Geyer,
B. V. Hugo,
A. Jameson,
D. J. Reardon,
R. M. Shannon,
R. Spiewak,
W. van Straten,
V. Venkatraman Krishnan
Abstract:
We present observations of 35 high spin-down energy radio pulsars using the MeerKAT telescope. Polarisation profiles and associated parameters are also presented. We derive the geometry for a selection of pulsars which show interpulse emission. We point out that, in several cases, these radio pulsars should also be seen in $γ$-rays but that improved radio timing is required to aid the high-energy…
▽ More
We present observations of 35 high spin-down energy radio pulsars using the MeerKAT telescope. Polarisation profiles and associated parameters are also presented. We derive the geometry for a selection of pulsars which show interpulse emission. We point out that, in several cases, these radio pulsars should also be seen in $γ$-rays but that improved radio timing is required to aid the high-energy detection. We discuss the relationship between the width of the radio profile and its high-energy detectability. Finally, we reflect on the correlation between the spin-down energy and the radio polarisation fraction and the implications this may have for $γ$-ray emission.
△ Less
Submitted 12 September, 2020;
originally announced September 2020.
-
The MeerKAT Telescope as a Pulsar Facility: System verification and early science results from MeerTime
Authors:
M. Bailes,
A. Jameson,
F. Abbate,
E. D. Barr,
N. D. R. Bhat,
L. Bondonneau,
M. Burgay,
S. J. Buchner,
F. Camilo,
D. J. Champion,
I. Cognard,
P. B. Demorest,
P. C. C. Freire,
T. Gautam,
M. Geyer,
J. M. Griessmeier,
L. Guillemot,
H. Hu,
F. Jankowski,
S. Johnston,
A. Karastergiou,
R. Karuppusamy,
D. Kaur,
M. J. Keith,
M. Kramer
, et al. (50 additional authors not shown)
Abstract:
We describe system verification tests and early science results from the pulsar processor (PTUSE) developed for the newly-commissioned 64-dish SARAO MeerKAT radio telescope in South Africa. MeerKAT is a high-gain (~2.8 K/Jy) low-system temperature (~18 K at 20cm) radio array that currently operates from 580-1670 MHz and can produce tied-array beams suitable for pulsar observations. This paper pres…
▽ More
We describe system verification tests and early science results from the pulsar processor (PTUSE) developed for the newly-commissioned 64-dish SARAO MeerKAT radio telescope in South Africa. MeerKAT is a high-gain (~2.8 K/Jy) low-system temperature (~18 K at 20cm) radio array that currently operates from 580-1670 MHz and can produce tied-array beams suitable for pulsar observations. This paper presents results from the MeerTime Large Survey Project and commissioning tests with PTUSE. Highlights include observations of the double pulsar J0737-3039A, pulse profiles from 34 millisecond pulsars from a single 2.5h observation of the Globular cluster Terzan 5, the rotation measure of Ter5O, a 420-sigma giant pulse from the Large Magellanic Cloud pulsar PSR J0540-6919, and nulling identified in the slow pulsar PSR J0633-2015. One of the key design specifications for MeerKAT was absolute timing errors of less than 5 ns using their novel precise time system. Our timing of two bright millisecond pulsars confirm that MeerKAT delivers exceptional timing. PSR J2241-5236 exhibits a jitter limit of <4 ns per hour whilst timing of PSR J1909-3744 over almost 11 months yields an rms residual of 66 ns with only 4 min integrations. Our results confirm that the MeerKAT is an exceptional pulsar telescope. The array can be split into four separate sub-arrays to time over 1000 pulsars per day and the future deployment of S-band (1750-3500 MHz) receivers will further enhance its capabilities.
△ Less
Submitted 28 May, 2020;
originally announced May 2020.
-
The Parkes Pulsar Timing Array Project: Second data release
Authors:
M. Kerr,
D. J. Reardon,
G. Hobbs,
R. M. Shannon,
R. N. Manchester,
S. Dai,
C. J. Russell,
S. -B. Zhang,
W. van Straten,
S. Osłowski,
A. Parthasarathy,
R. Spiewak,
M. Bailes,
N. D. R. Bhat,
A. D. Cameron,
W. A. Coles,
J. Dempsey,
X. Deng,
B. Goncharov,
J. F Kaczmarek,
M. J. Keith,
P. D. Lasky,
M. E. Lower,
B. Preisig,
J. M. Sarkissian
, et al. (5 additional authors not shown)
Abstract:
We describe 14 years of public data from the Parkes Pulsar Timing Array (PPTA), an ongoing project that is producing precise measurements of pulse times of arrival from 26 millisecond pulsars using the 64-m Parkes radio telescope with a cadence of approximately three weeks in three observing bands. A comprehensive description of the pulsar observing systems employed at the telescope since 2004 is…
▽ More
We describe 14 years of public data from the Parkes Pulsar Timing Array (PPTA), an ongoing project that is producing precise measurements of pulse times of arrival from 26 millisecond pulsars using the 64-m Parkes radio telescope with a cadence of approximately three weeks in three observing bands. A comprehensive description of the pulsar observing systems employed at the telescope since 2004 is provided, including the calibration methodology and an analysis of the stability of system components. We attempt to provide full accounting of the reduction from the raw measured Stokes parameters to pulse times of arrival to aid third parties in reproducing our results. This conversion is encapsulated in a processing pipeline designed to track provenance. Our data products include pulse times of arrival for each of the pulsars along with an initial set of pulsar parameters and noise models. The calibrated pulse profiles and timing template profiles are also available. These data represent almost 21,000 hrs of recorded data spanning over 14 years. After accounting for processes that induce time-correlated noise, 22 of the pulsars have weighted root-mean-square timing residuals of < 1 $μ$s in at least one radio band. The data should allow end users to quickly undertake their own gravitational-wave analyses (for example) without having to understand the intricacies of pulsar polarisation calibration or attain a mastery of radio-frequency interference mitigation as is required when analysing raw data files.
△ Less
Submitted 21 March, 2020;
originally announced March 2020.
-
The Thousand-Pulsar-Array programme on MeerKAT I: Science objectives and first results
Authors:
Simon Johnston,
A. Karastergiou,
M. J. Keith,
X. Song,
P. Weltevrede,
F. Abbate,
M. Bailes,
S. Buchner,
F. Camilo,
M. Geyer,
B. Hugo,
A. Jameson. M. Kramer,
A. Parthasarathy,
D. J. Reardon,
A. Ridolfi,
M. Serylak,
R. M. Shannon,
R. Spiewak,
W. van Straten,
V. Venkatraman Krishnan,
F. Jankowski,
B. W. Meyers,
L. Oswald,
B. Posselt,
C. Sobey
, et al. (2 additional authors not shown)
Abstract:
We report here on initial results from the Thousand Pulsar Array (TPA) programme, part of the Large Survey Project "MeerTime" on the MeerKAT telescope. The interferometer is used in tied-array mode in the band from 856 to 1712~MHz, and the wide band coupled with the large collecting area and low receiver temperature make it an excellent telescope for the study of radio pulsars. The TPA is a 5 year…
▽ More
We report here on initial results from the Thousand Pulsar Array (TPA) programme, part of the Large Survey Project "MeerTime" on the MeerKAT telescope. The interferometer is used in tied-array mode in the band from 856 to 1712~MHz, and the wide band coupled with the large collecting area and low receiver temperature make it an excellent telescope for the study of radio pulsars. The TPA is a 5 year project which aims to observe (a) more than 1000 pulsars to obtain high-fidelity pulse profiles, (b) some 500 of these pulsars over multiple epochs, (c) long sequences of single-pulse trains from several hundred pulsars. The scientific outcomes from the programme will include determination of pulsar geometries, the location of the radio emission within the pulsar magnetosphere, the connection between the magnetosphere and the crust and core of the star, tighter constraints on the nature of the radio emission itself as well as interstellar medium studies. First results presented here include updated dispersion measures, 26 pulsars with Faraday rotation measures derived for the first time and a description of interesting emission phenomena observed thus far.
△ Less
Submitted 18 February, 2020;
originally announced February 2020.
-
A pulsar-based timescale from the International Pulsar Timing Array
Authors:
G. Hobbs,
L. Guo,
R. N. Caballero,
W. Coles,
K. J. Lee,
R. N. Manchester,
D. J. Reardon,
D. Matsakis,
M. L. Tong,
Z. Arzoumanian,
M. Bailes,
C. G. Bassa,
N. D. R. Bhat,
A. Brazier,
S. Burke-Spolaor,
D. J. Champion,
S. Chatterjee,
I. Cognard,
S. Dai,
G. Desvignes,
T. Dolch,
R. D. Ferdman,
E. Graikou,
L. Guillemot,
G. H. Janssen
, et al. (34 additional authors not shown)
Abstract:
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…
▽ More
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.
△ Less
Submitted 29 October, 2019;
originally announced October 2019.
-
The International Pulsar Timing Array: Second data release
Authors:
B. B. P. Perera,
M. E. DeCesar,
P. B. Demorest,
M. Kerr,
L. Lentati,
D. J. Nice,
S. Oslowski,
S. M. Ransom,
M. J. Keith,
Z. Arzoumanian,
M. Bailes,
P. T. Baker,
C. G. Bassa,
N. D. R. Bhat,
A. Brazier,
M. Burgay,
S. Burke-Spolaor,
R. N. Caballero,
D. J. Champion,
S. Chatterjee,
S. Chen,
I. Cognard,
J. M. Cordes,
K. Crowter,
S. Dai
, et al. (50 additional authors not shown)
Abstract:
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 a…
▽ More
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.
△ Less
Submitted 10 September, 2019;
originally announced September 2019.
-
Modelling annual and orbital variations in the scintillation of the relativistic binary PSR J1141$-$6545
Authors:
D. J. Reardon,
W. A. Coles,
G. Hobbs,
S. Ord,
M. Kerr,
M. Bailes,
N. D. R. Bhat,
V. Venkatraman Krishnan
Abstract:
We have observed the relativistic binary pulsar PSR J1141$-$6545 over a period of $\sim$6 years using the Parkes 64 m radio telescope, with a focus on modelling the diffractive intensity scintillations to improve the accuracy of the astrometric timing model. The long-term scintillation, which shows orbital and annual variations, allows us to measure parameters that are difficult to measure with pu…
▽ More
We have observed the relativistic binary pulsar PSR J1141$-$6545 over a period of $\sim$6 years using the Parkes 64 m radio telescope, with a focus on modelling the diffractive intensity scintillations to improve the accuracy of the astrometric timing model. The long-term scintillation, which shows orbital and annual variations, allows us to measure parameters that are difficult to measure with pulsar timing alone. These include: the orbital inclination $i$; the longitude of the ascending node $Ω$; and the pulsar system transverse velocity. We use the annual variations to resolve the previous ambiguity in the sense of the inclination angle. Using the correct sense, and a prior probability distribution given by a constraint from pulsar timing ($i=73\pm3^\circ$), we find $Ω=24.8\pm1.8^\circ$ and we estimate the pulsar distance to be $D=10^{+4}_{-3}$ kpc. This then gives us an estimate of this pulsar's proper motion of $μ_α\cosδ=2.9\pm1.0$ mas yr$^{-1}$ in right ascension and $μ_δ=1.8\pm0.6$ mas yr$^{-1}$ in declination. Finally, we obtain measurements of the spatial structure of the interstellar electron density fluctuations, including: the spatial scale and anisotropy of the diffraction pattern; the distribution of scattering material along the line of sight; and spatial variation in the strength of turbulence from epoch to epoch. We find that the scattering is dominated by a thin screen at a distance of $(0.724\pm0.008)D$, with an anisotropy axial ratio $A_{\rm r} = 2.14\pm0.11$.
△ Less
Submitted 5 March, 2019;
originally announced March 2019.
-
Parkes Pulsar Timing Array constraints on ultralight scalar-field dark matter
Authors:
Nataliya K. Porayko,
Xingjiang Zhu,
Yuri Levin,
Lam Hui,
George Hobbs,
Aleksandra Grudskaya,
Konstantin Postnov,
Matthew Bailes,
N. D. Ramesh Bhat,
William Coles,
Shi Dai,
James Dempsey,
Michael J. Keith,
Matthew Kerr,
Michael Kramer,
Paul D. Lasky,
Richard N. Manchester,
Stefan Osłowski,
Aditya Parthasarathy,
Vikram Ravi,
Daniel J. Reardon,
Pablo A. Rosado,
Christopher J. Russell,
Ryan M. Shannon,
Renée Spiewak
, et al. (5 additional authors not shown)
Abstract:
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, pr…
▽ More
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.
△ Less
Submitted 16 December, 2019; v1 submitted 7 October, 2018;
originally announced October 2018.