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 characteristics and observational history with other telescopes. In this paper, we provide an outline of the programme, and present polarization 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, 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.
Context . Scalar-tensor gravity (STG) theories are well-motivated alternatives to general relativity (GR). One class of STG theories, Damour–Esposito–Farèse (DEF) gravity, has a massless scalar field ...with two arbitrary coupling parameters. We are interested in this theory because, despite its simplicity, it predicts a wealth of different phenomena, such as dipolar gravitational wave emission and spontaneous scalarisation of neutron stars (NSs). These phenomena of DEF gravity can be tested by timing binary radio pulsars. In the methods used so far, intermediate phenomenological post-Keplerian (PK) parameters are measured by fitting the corresponding timing model to the timing data whose values are then compared to the predictions from the alternative theory being tested. However, this approach loses information between intermediate steps and does not account for possible correlations between PK parameters. Aims . We aim to develop a new binary pulsar timing model ‘DDSTG’ (called after Damour, Deruelle and STG) to enable more precise tests of STG theories based on a minimal set of binary parameters. The expressions for PK parameters in DEF gravity are self-consistently incorporated into the model. PK parameters depend on two masses which are now directly fitted to the data without intermediate steps. The new technique takes into account all possible correlations between PK parameters naturally. Methods . Grids of physical parameters of NSs were calculated in the framework of DEF gravity for a set of 11 equations of state. Automatic differentiation (AutoDiff) technique was employed, which aids in the calculation of gravitational form factors of NSs with a higher precision than in previous works. The pulsar timing program TEMPO was selected as a framework for the realisation of the DDSTG model. The implemented model is applicable to any type of pulsar companions. We also simulated realistic future radio-timing datasets for a number of large radio observatories for the binary pulsar PSR J2222-0137 and three generic pulsar-black hole (PSR-BH) systems. Results . We applied the DDSTG model to the most recently published observational data for PSR J2222-0137. The obtained limits on DEF gravity parameters for this system confirm and improve previous results. New limits are also the most reliable because DEF gravity is directly fitted to the data. We argue that future observations of PSR J2222-0137 can significantly improve the limits and that PSR-BH systems have the potential to place the tightest limits in certain areas of the DEF gravity parameter space.
The sensitivity of Pulsar Timing Arrays to gravitational waves (GWs) depends on the noise present in the individual pulsar timing data. Noise may be either intrinsic or extrinsic to the pulsar. ...Intrinsic sources of noise will include rotational instabilities, for example. Extrinsic sources of noise include contributions from physical processes which are not sufficiently well modelled, for example, dispersion and scattering effects, analysis errors and instrumental instabilities. We present the results from a noise analysis for 42 millisecond pulsars (MSPs) observed with the European Pulsar Timing Array. For characterizing the low-frequency, stochastic and achromatic noise component, or ‘timing noise’, we employ two methods, based on Bayesian and frequentist statistics. For 25 MSPs, we achieve statistically significant measurements of their timing noise parameters and find that the two methods give consistent results. For the remaining 17 MSPs, we place upper limits on the timing noise amplitude at the 95 per cent confidence level. We additionally place an upper limit on the contribution to the pulsar noise budget from errors in the reference terrestrial time standards (below 1 per cent), and we find evidence for a noise component which is present only in the data of one of the four used telescopes. Finally, we estimate that the timing noise of individual pulsars reduces the sensitivity of this data set to an isotropic, stochastic GW background by a factor of >9.1 and by a factor of >2.3 for continuous GWs from resolvable, inspiralling supermassive black hole binaries with circular orbits.
We report the discovery of 1.97 ms period gamma-ray pulsations from the 75 minute orbital-period binary pulsar now named PSR J1653−0158. The associated Fermi Large Area Telescope gamma-ray source ...4FGL J1653.6−0158 has long been expected to harbor a binary millisecond pulsar. Despite the pulsar-like gamma-ray spectrum and candidate optical/X-ray associations-whose periodic brightness modulations suggested an orbit-no radio pulsations had been found in many searches. The pulsar was discovered by directly searching the gamma-ray data using the GPU-accelerated Einstein@Home distributed volunteer computing system. The multidimensional parameter space was bounded by positional and orbital constraints obtained from the optical counterpart. More sensitive analyses of archival and new radio data using knowledge of the pulsar timing solution yield very stringent upper limits on radio emission. Any radio emission is thus either exceptionally weak, or eclipsed for a large fraction of the time. The pulsar has one of the three lowest inferred surface magnetic-field strengths of any known pulsar with Bsurf 4 × 107 G. The resulting mass function, combined with models of the companion star's optical light curve and spectra, suggests a pulsar mass 2 M . The companion is lightweight with mass ∼0.01 M , and the orbital period is the shortest known for any rotation-powered binary pulsar. This discovery demonstrates the Fermi Large Area Telescope's potential to discover extreme pulsars that would otherwise remain undetected.
Context.
NenuFAR (New extension in Nançay upgrading LOFAR) is a new radio telescope developed and built on the site of the Nançay Radio Observatory. It is designed to observe the largely unexplored ...frequency window from 10 to 85 MHz, offering a high sensitivity across its full bandwidth. NenuFAR has started its “early science” operation in July 2019, with 58% of its final collecting area.
Aims.
Pulsars are one of the major phenomena utilized in the scientific exploitation of this frequency range and represent an important challenge in terms of instrumentation. Designing instrumentation at these frequencies is complicated by the need to compensate for the effects of both the interstellar medium and the ionosphere on the observed signal. We have designed a dedicated backend and developed a complete pulsar observation and data analysis pipeline, which we describe in detail in the present paper, together with first science results illustrating the diversity of the pulsar observing modes.
Methods.
Our real-time pipeline LUPPI (Low frequency Ultimate Pulsar Processing Instrumentation) is able to cope with a high data rate and provide real-time coherent de-dispersion down to the lowest frequencies reached by NenuFAR (10 MHz). The full backend functionality is described, as the available pulsar observing modes (folded, single-pulse, waveform, and dynamic spectrum).
Results.
We also present some of the early science results of NenuFAR on pulsars: the detection of 12 millisecond pulsars (eight of which are detected for the first time below 100 MHz); a high-frequency resolution mapping of the PSR B1919+21 emission profile and a detailed observation of single-pulse substructures from PSR B0809+74 down to 16 MHz; the high rate of giant-pulse emission from the Crab pulsar detected at 68.7 MHz (43 events per minute); and the illustration of the very good timing performance of the instrumentation, which allows us to study dispersion measure variations in great detail.
We report on 22 yr of radio timing observations of the millisecond pulsar J1024−0719 by the telescopes participating in the European Pulsar Timing Array (EPTA). These observations reveal a ...significant second derivative of the pulsar spin frequency and confirm the discrepancy between the parallax and Shklovskii distances that has been reported earlier. We also present optical astrometry, photometry and spectroscopy of 2MASS J10243869−0719190. We find that it is a low-metallicity main-sequence star (K7V spectral type, M/H = −1.0, T
eff = 4050 ± 50 K) and that its position, proper motion and distance are consistent with those of PSR J1024−0719. We conclude that PSR J1024−0719 and 2MASS J10243869−0719190 form a common proper motion pair and are gravitationally bound. The gravitational interaction between the main-sequence star and the pulsar accounts for the spin frequency derivatives, which in turn resolves the distance discrepancy. Our observations suggest that the pulsar and main-sequence star are in an extremely wide (P
b > 200 yr) orbit. Combining the radial velocity of the companion and proper motion of the pulsar, we find that the binary system has a high spatial velocity of 384 ± 45 km s−1 with respect to the local standard of rest and has a Galactic orbit consistent with halo objects. Since the observed main-sequence companion star cannot have recycled the pulsar to millisecond spin periods, an exotic formation scenario is required. We demonstrate that this extremely wide-orbit binary could have evolved from a triple system that underwent an asymmetric supernova explosion, though find that significant fine-tuning during the explosion is required. Finally, we discuss the implications of the long period orbit on the timing stability of PSR J1024−0719 in light of its inclusion in pulsar timing arrays.
Context.
PSR J1528−3146 is a 60.8 ms pulsar orbiting a heavy white dwarf (WD) companion, with an orbital period of 3.18 d. The pulsar was discovered in the early 2000 s in a survey at 1.4 GHz of ...intermediate Galactic latitudes conducted with the Parkes radio telescope. The initial timing analysis of PSR J1528−3146, using data recorded from 2001 and 2004, did not reveal any relativistic perturbations to the orbit of the pulsar or to the propagation of its pulses. However, with an orbital eccentricity of ∼0.0002 and a large companion mass on the order of 1
M
⊙
, this system has been deemed likely to exhibit measurable perturbations.
Aims.
This work is aimed at characterizing the pulsar’s astrometric, spin, and orbital parameters by analyzing timing measurements conducted at the Parkes, MeerKAT, and Nançay radio telescopes over nearly two decades. The measurement of post-Keplerian perturbations to the pulsar’s orbit can be used to constrain the masses of the two component stars of the binary and, in turn, to offer insights into the history of the system.
Methods.
We analyzed timing data from the Parkes, MeerKAT, and Nançay radio telescopes collected over about 16 yr, obtaining a precise rotation ephemeris for PSR J1528−3146. A Bayesian analysis of the timing data was carried out to constrain the masses of the two components and the orientation of the orbit. We further analyzed the polarization properties of the pulsar to constrain the orientation of the magnetic axis and of the line of sight with respect to the spin axis.
Results.
We measured a significant rate of advance of periastron, for the first time, and we set constraints on the Shapiro delay in the system and on the rate of change of the projected semi-major axis of the pulsar’s orbit. The Bayesian analysis yielded measurements for the pulsar and companion masses of
M
p
= 1.61
−0.13
+0.14
M
⊙
and
M
c
= 1.33
−0.07
+0.08
M
⊙
(68% C.L.), respectively, confirming that the companion is indeed massive. This companion mass as well as other characteristics of PSR J1528−3146 indicate that this pulsar is very similar to PSR J2222−0137, a 32.8 ms pulsar orbiting a WD whose heavy mass (∼1.32
M
⊙
) has been considered unique among pulsar-WD systems until now. Our measurements suggest common evolutionary scenarios for PSRs J1528−3146 and J2222−0137.
ABSTRACT
Recently, global pulsar timing arrays have released results from searching for a nano-Hertz gravitational wave background signal. Although there has not been any definite evidence of the ...presence of such a signal in residuals of pulsar timing data yet, with more and improved data in future, a statistically significant detection is expected to be made. Stochastic algorithms are used to sample a very large parameter space to infer results from data. In this paper, we attempt to rule out effects arising from the stochasticity of the sampler in the inference process. We compare different configurations of nested samplers and the more commonly used markov chain monte carlo method to sample the pulsar timing array parameter space and account for times taken by the different samplers on same data. Although we obtain consistent results on parameters from different sampling algorithms, we propose two different samplers for robustness checks on data in the future to account for cross-checks between sampling methods as well as realistic run-times.
Abstract
We report the discovery and the results of follow-up timing observations of PSR J2045+3633 and PSR J2053+4650, two binary pulsars found in the Northern High Time Resolution Universe pulsar ...survey being carried out with the Effelsberg radio telescope. Having spin periods of 31.7 and 12.6 ms, respectively, and both with massive white dwarf companions, M
c > 0.8 M⊙, the pulsars can be classified as mildly recycled. PSR J2045+3633 is remarkable due to its orbital period (32.3 d) and eccentricity e = 0.017 212 44(5), which is amongst the largest ever measured for this class. After almost two years of timing, the large eccentricity has allowed the measurement of the rate of advance of periastron at the 5σ level, 0.0010(2)°yr− 1. Combining this with a detection of the orthometric amplitude of the Shapiro delay, we obtained the following constraints on the component masses (within general relativity):
$M_{\text{p}}\, = \, 1.33^{+0.30}_{-0.28}\, \mathrm{M}_{{\odot }}$
; and
$M_{\text{c}}\, = \, 0.94^{+0.14}_{-0.13}\, \mathrm{M}_{{\odot }}$
. PSR J2053+4650 has a 2.45 d circular orbit inclined to the plane of the sky at an angle
$i\, = \, 85.0^{+0.8}_{-0.9}\,{\rm deg}$
. In this nearly edge-on case the masses can be obtained from the Shapiro delay alone. Our timing observations resulted in a significant detection of this effect giving:
$M_{\text{p}}\, = \, 1.40^{+0.21}_{-0.18}\, \mathrm{M}_{{\odot }}$
; and
$M_{\text{c}}\, = \, 0.86^{+0.07}_{-0.06}\, \mathrm{M}_{{\odot }}$
.