Showing 1–21 of 21 results for author: Champion, D J

Tackling artefacts in the timing of relativistic pulsar binaries: towards the SKA

Authors: Huanchen Hu, Nataliya K. Porayko, Willem van Straten, Michael Kramer, David J. Champion, Michael J. Keith

Abstract: Common signal-processing approximations produce artefacts when timing pulsars in relativistic binary systems, especially edge-on systems with tight orbits, such as the Double Pulsar. In this paper, we use extensive simulations to explore various patterns that arise from the inaccuracies of approximations made when correcting dispersion and Shapiro delay. In a relativistic binary, the velocity of t… ▽ More Common signal-processing approximations produce artefacts when timing pulsars in relativistic binary systems, especially edge-on systems with tight orbits, such as the Double Pulsar. In this paper, we use extensive simulations to explore various patterns that arise from the inaccuracies of approximations made when correcting dispersion and Shapiro delay. In a relativistic binary, the velocity of the pulsar projected onto the line-of-sight varies significantly on short time scales, causing rapid changes in the apparent pulsar spin frequency, which is used to convert dispersive delays to pulsar rotational phase shifts. A well-known example of the consequences of this effect is the artificial variation of dispersion measure (DM) with binary phase, first observed in the Double Pulsar 20 years ago. We show that ignoring the Doppler shift of the spin frequency when computing the dispersive phase shift exactly reproduces the shape and magnitude of the reported DM variations. We also simulate and study two additional effects of much smaller magnitude, which are caused by the assumption that the spin frequency used to correct dispersion is constant over the duration of the sub-integration and over the observed bandwidth. We show that failure to account for these two effects leads to orbital phase-dependent dispersive smearing that leads to apparent orbital DM variations. The functional form of the variation depends on the orbital eccentricity. In addition, we find that a polynomial approximation of the timing model is unable to accurately describe the Shapiro delay of edge-on systems with orbits less than 4 hours, which poses problems for the measurements of timing parameters, most notably the Shapiro delay. This will be a potential issue for sensitive facilities like the FAST and the forthcoming Square Kilometre Array (SKA); therefore, a more accurate phase predictor is indispensable. △ Less

Submitted 27 January, 2025; originally announced January 2025.

Comments: 13 pages, 17 figures. Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 694, A243 (2025)

Modeling non stationary noise in pulsar timing array data analysis

Authors: Mikel Falxa, J. Antoniadis, D. J. Champion, I. Cognard, G. Desvignes, L. Guillemot, H. Hu, G. Janssen, J. Jawor, R. Karuppusamy, M. J. Keith, M. Kramer, K. Lackeos, K. Liu, J. W. McKee, D. Perrodin, S. A. Sanidas, G. M. Shaifullah, G. Theureau

Abstract: Pulsar Timing Array (PTA) collaborations recently reported evidence for the presence of a gravitational wave background (GWB) in their datasets. The main candidate that is expected to produce such a GWB is the population of supermassive black hole binaries (SMBHB). Some analyses showed that the recovered signal may exhibit time-dependent properties, i.e. non-stationarity. In this paper, we propose… ▽ More Pulsar Timing Array (PTA) collaborations recently reported evidence for the presence of a gravitational wave background (GWB) in their datasets. The main candidate that is expected to produce such a GWB is the population of supermassive black hole binaries (SMBHB). Some analyses showed that the recovered signal may exhibit time-dependent properties, i.e. non-stationarity. In this paper, we propose an approximated non-stationary Gaussian process (GP) model obtained from the perturbation of stationary processes. The presented method is applied to the second data release of the European pulsar timing array to search for non-stationary features in the GWB. We analyzed the data in different time slices and showed that the inferred properties of the GWB evolve with time. We find no evidence for such non-stationary behavior and the Bayes factor in favor of the latter is $\mathcal{B}^{NS}_{S} = 1.5$. We argue that the evolution of the GWB properties most likely comes from the \mf{improvement of the observation cadence} with time and \mf{better} characterization of the noise of individual pulsars. Such non-stationary GWB could also be produced by the leakage of non-stationary features in the noise of individual pulsars or by the presence of an eccentric single source. △ Less

Submitted 6 May, 2024; originally announced May 2024.

Comments: 14 pages, 14 figures

Constraints on conformal ultralight dark matter couplings from the European Pulsar Timing Array

Authors: Clemente Smarra, Adrien Kuntz, Enrico Barausse, Boris Goncharov, Diana López Nacir, Diego Blas, Lijing Shao, J. Antoniadis, D. J. Champion, I. Cognard, L. Guillemot, H. Hu, M. Keith, M. Kramer, K. Liu, D. Perrodin, S. A. Sanidas, G. Theureau

Abstract: Millisecond pulsars are extremely precise celestial clocks: as they rotate, the beamed radio waves emitted along the axis of their magnetic field can be detected with radio telescopes, which allows for tracking subtle changes in the pulsars' rotation periods. A possible effect on the period of a pulsar is given by a potential coupling to dark matter, in cases where it is modeled with an "ultraligh… ▽ More Millisecond pulsars are extremely precise celestial clocks: as they rotate, the beamed radio waves emitted along the axis of their magnetic field can be detected with radio telescopes, which allows for tracking subtle changes in the pulsars' rotation periods. A possible effect on the period of a pulsar is given by a potential coupling to dark matter, in cases where it is modeled with an "ultralight" scalar field. In this paper, we consider a universal conformal coupling of the dark matter scalar to gravity, which in turn mediates an effective coupling between pulsars and dark matter. If the dark matter scalar field is changing in time, as expected in the Milky Way, this effective coupling produces a periodic modulation of the pulsar rotational frequency. By studying the time series of observed radio pulses collected by the European Pulsar Timing Array experiment, we present constraints on the coupling of dark matter, improving on existing bounds. These bounds can also be regarded as constraints on the parameters of scalar-tensor theories of the Fierz-Jordan-Brans-Dicke and Damour-Esposito-Farèse types in the presence of a (light) mass potential term. △ Less

Submitted 4 October, 2024; v1 submitted 2 May, 2024; originally announced May 2024.

Comments: 17 pages, 5 figures, Updated plots, matching published version

Journal ref: Phys. Rev. D 110 (2024) 043033

A pulsar in a binary with a compact object in the mass gap between neutron stars and black holes

Authors: Ewan D. Barr, Arunima Dutta, Paulo C. C. Freire, Mario Cadelano, Tasha Gautam, Michael Kramer, Cristina Pallanca, Scott M. Ransom, Alessandro Ridolfi, Benjamin W. Stappers, Thomas M. Tauris, Vivek Venkatraman Krishnan, Norbert Wex, Matthew Bailes, Jan Behrend, Sarah Buchner, Marta Burgay, Weiwei Chen, David J. Champion, C. -H. Rosie Chen, Alessandro Corongiu, Marisa Geyer, Y. P. Men, Prajwal V. Padmanabh, Andrea Possenti

Abstract: Among the compact objects observed in gravitational wave merger events a few have masses in the gap between the most massive neutron stars (NSs) and least massive black holes (BHs) known. Their nature and the formation of their merging binaries are not well understood. We report on pulsar timing observations using the Karoo Array Telescope (MeerKAT) of PSR J0514-4002E, an eccentric binary millisec… ▽ More Among the compact objects observed in gravitational wave merger events a few have masses in the gap between the most massive neutron stars (NSs) and least massive black holes (BHs) known. Their nature and the formation of their merging binaries are not well understood. We report on pulsar timing observations using the Karoo Array Telescope (MeerKAT) of PSR J0514-4002E, an eccentric binary millisecond pulsar in the globular cluster NGC 1851 with a total binary mass of $3.887 \pm 0.004$ solar masses. The companion to the pulsar is a compact object and its mass (between $2.09$ and $2.71$ solar masses, 95% confidence interval) is in the mass gap, so it either is a very massive NS or a low-mass BH. We propose the companion was formed by a merger between two earlier NSs. △ Less

Submitted 18 January, 2024; originally announced January 2024.

Comments: 41 pages, 13 figures, 3 tables, to be published in Science

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.

Comments: Abridged abstract. 13 pages, 9 figures and 2 tables

Journal ref: A&A 682, A26 (2024)

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.

Comments: 21 pages, 9 figures, submitted to ApJ

The second data release from the European Pulsar Timing Array: VI. Challenging the ultralight dark matter paradigm

Authors: Clemente Smarra, Boris Goncharov, Enrico Barausse, J. Antoniadis, S. Babak, A. -S. Bak Nielsen, C. G. Bassa, A. Berthereau, M. Bonetti, E. Bortolas, P. R. Brook, M. Burgay, R. N. Caballero, A. Chalumeau, D. J. Champion, S. Chanlaridis, S. Chen, I. Cognard, G. Desvignes, M. Falxa, R. D. Ferdman, A. Franchini, J. R. Gair, E. Graikou, J. -M. Grie , et al. (46 additional authors not shown)

Abstract: Pulsar Timing Array experiments probe the presence of possible scalar or pseudoscalar ultralight dark matter particles through decade-long timing of an ensemble of galactic millisecond radio pulsars. With the second data release of the European Pulsar Timing Array, we focus on the most robust scenario, in which dark matter interacts only gravitationally with ordinary baryonic matter. Our results s… ▽ More Pulsar Timing Array experiments probe the presence of possible scalar or pseudoscalar ultralight dark matter particles through decade-long timing of an ensemble of galactic millisecond radio pulsars. With the second data release of the European Pulsar Timing Array, we focus on the most robust scenario, in which dark matter interacts only gravitationally with ordinary baryonic matter. Our results show that ultralight particles with masses $10^{-24.0}~\text{eV} \lesssim m \lesssim 10^{-23.3}~\text{eV}$ cannot constitute $100\%$ of the measured local dark matter density, but can have at most local density $ρ\lesssim 0.3$ GeV/cm$^3$. △ Less

Submitted 25 October, 2023; v1 submitted 28 June, 2023; originally announced June 2023.

Comments: 5 pages + acknowledgements + refs, 2 figures. Text and figures match the version published in PRL

The second data release from the European Pulsar Timing Array: IV. Implications for massive black holes, dark matter and the early Universe

Authors: J. Antoniadis, P. Arumugam, S. Arumugam, P. Auclair, S. Babak, M. Bagchi, A. -S. Bak Nielsen, E. Barausse, C. G. Bassa, A. Bathula, A. Berthereau, M. Bonetti, E. Bortolas, P. R. Brook, M. Burgay, R. N. Caballero, C. Caprini, A. Chalumeau, D. J. Champion, S. Chanlaridis, S. Chen, I. Cognard, M. Crisostomi, S. Dandapat, D. Deb , et al. (89 additional authors not shown)

Abstract: The European Pulsar Timing Array (EPTA) and Indian Pulsar Timing Array (InPTA) collaborations have measured a low-frequency common signal in the combination of their second and first data releases respectively, with the correlation properties of a gravitational wave background (GWB). Such signal may have its origin in a number of physical processes including a cosmic population of inspiralling sup… ▽ More The European Pulsar Timing Array (EPTA) and Indian Pulsar Timing Array (InPTA) collaborations have measured a low-frequency common signal in the combination of their second and first data releases respectively, with the correlation properties of a gravitational wave background (GWB). Such signal may have its origin in a number of physical processes including a cosmic population of inspiralling supermassive black hole binaries (SMBHBs); inflation, phase transitions, cosmic strings and tensor mode generation by non-linear evolution of scalar perturbations in the early Universe; oscillations of the Galactic potential in the presence of ultra-light dark matter (ULDM). At the current stage of emerging evidence, it is impossible to discriminate among the different origins. Therefore, in this paper, we consider each process separately, and investigate the implications of the signal under the hypothesis that it is generated by that specific process. We find that the signal is consistent with a cosmic population of inspiralling SMBHBs, and its relatively high amplitude can be used to place constraints on binary merger timescales and the SMBH-host galaxy scaling relations. If this origin is confirmed, this is the first direct evidence that SMBHBs merge in nature, adding an important observational piece to the puzzle of structure formation and galaxy evolution. As for early Universe processes, the measurement would place tight constraints on the cosmic string tension and on the level of turbulence developed by first-order phase transitions. Other processes would require non-standard scenarios, such as a blue-tilted inflationary spectrum or an excess in the primordial spectrum of scalar perturbations at large wavenumbers. Finally, a ULDM origin of the detected signal is disfavoured, which leads to direct constraints on the abundance of ULDM in our Galaxy. △ Less

Submitted 15 May, 2024; v1 submitted 28 June, 2023; originally announced June 2023.

Comments: 30 pages, 23 figures, replaced to match the version published in Astronomy & Astrophysics, note the change in the numbering order in the series (now paper IV)

The second data release from the European Pulsar Timing Array V. Search for continuous gravitational wave signals

Authors: J. Antoniadis, P. Arumugam, S. Arumugam, S. Babak, M. Bagchi, A. S. Bak Nielsen, C. G. Bassa, A. Bathula, A. Berthereau, M. Bonetti, E. Bortolas, P. R. Brook, M. Burgay, R. N. Caballero, A. Chalumeau, D. J. Champion, S. Chanlaridis, S. Chen, I. Cognard, S. Dandapat, D. Deb, S. Desai, G. Desvignes, N. Dhanda-Batra, C. Dwivedi , et al. (75 additional authors not shown)

Abstract: We present the results of a search for continuous gravitational wave signals (CGWs) in the second data release (DR2) of the European Pulsar Timing Array (EPTA) collaboration. The most significant candidate event from this search has a gravitational wave frequency of 4-5 nHz. Such a signal could be generated by a supermassive black hole binary (SMBHB) in the local Universe. We present the results o… ▽ More We present the results of a search for continuous gravitational wave signals (CGWs) in the second data release (DR2) of the European Pulsar Timing Array (EPTA) collaboration. The most significant candidate event from this search has a gravitational wave frequency of 4-5 nHz. Such a signal could be generated by a supermassive black hole binary (SMBHB) in the local Universe. We present the results of a follow-up analysis of this candidate using both Bayesian and frequentist methods. The Bayesian analysis gives a Bayes factor of 4 in favor of the presence of the CGW over a common uncorrelated noise process, while the frequentist analysis estimates the p-value of the candidate to be 1%, also assuming the presence of common uncorrelated red noise. However, comparing a model that includes both a CGW and a gravitational wave background (GWB) to a GWB only, the Bayes factor in favour of the CGW model is only 0.7. Therefore, we cannot conclusively determine the origin of the observed feature, but we cannot rule it out as a CGW source. We present results of simulations that demonstrate that data containing a weak gravitational wave background can be misinterpreted as data including a CGW and vice versa, providing two plausible explanations of the EPTA DR2 data. Further investigations combining data from all PTA collaborations will be needed to reveal the true origin of this feature. △ Less

Submitted 25 June, 2024; v1 submitted 28 June, 2023; originally announced June 2023.

Comments: 13 figures, 15 pages, accepted

Journal ref: A&A 690, A118 (2024)

The second data release from the European Pulsar Timing Array I. The dataset and timing analysis

Authors: J. Antoniadis, S. Babak, A. -S. Bak Nielsen, C. G. Bassa, A. Berthereau, M. Bonetti, E. Bortolas, P. R. Brook, M. Burgay, R. N. Caballero, A. Chalumeau, D. J. Champion, S. Chanlaridis, S. Chen, I. Cognard, G. Desvignes, M. Falxa, R. D. Ferdman, A. Franchini, J. R. Gair, B. Goncharov, E. Graikou, J. -M. Grießmeier, L. Guillemot, Y. J. Guo , et al. (44 additional authors not shown)

Abstract: Pulsar timing arrays offer a probe of the low-frequency gravitational wave spectrum (1 - 100 nanohertz), which is intimately connected to a number of markers that can uniquely trace the formation and evolution of the Universe. We present the dataset and the results of the timing analysis from the second data release of the European Pulsar Timing Array (EPTA). The dataset contains high-precision pu… ▽ More Pulsar timing arrays offer a probe of the low-frequency gravitational wave spectrum (1 - 100 nanohertz), which is intimately connected to a number of markers that can uniquely trace the formation and evolution of the Universe. We present the dataset and the results of the timing analysis from the second data release of the European Pulsar Timing Array (EPTA). The dataset contains high-precision pulsar timing data from 25 millisecond pulsars collected with the five largest radio telescopes in Europe, as well as the Large European Array for Pulsars. The dataset forms the foundation for the search for gravitational waves by the EPTA, presented in associated papers. We describe the dataset and present the results of the frequentist and Bayesian pulsar timing analysis for individual millisecond pulsars that have been observed over the last ~25 years. We discuss the improvements to the individual pulsar parameter estimates, as well as new measurements of the physical properties of these pulsars and their companions. This data release extends the dataset from EPTA Data Release 1 up to the beginning of 2021, with individual pulsar datasets with timespans ranging from 14 to 25 years. These lead to improved constraints on annual parallaxes, secular variation of the orbital period, and Shapiro delay for a number of sources. Based on these results, we derived astrophysical parameters that include distances, transverse velocities, binary pulsar masses, and annual orbital parallaxes. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 29 pages, 9 figures, 13 tables, Astronomy & Astrophysics in press

Practical approaches to analyzing PTA data: Cosmic strings with six pulsars

Authors: Hippolyte Quelquejay Leclere, Pierre Auclair, Stanislav Babak, Aurélien Chalumeau, Danièle A. Steer, J. Antoniadis, A. -S. Bak Nielsen, C. G. Bassa, A. Berthereau, M. Bonetti, E. Bortolas, P. R. Brook, M. Burgay, R. N. Caballero, D. J. Champion, S. Chanlaridis, S. Chen, I. Cognard, G. Desvignes, M. Falxa, R. D. Ferdman, A. Franchini, J. R. Gair, B. Goncharov, E. Graikou , et al. (47 additional authors not shown)

Abstract: We search for a stochastic gravitational wave background (SGWB) generated by a network of cosmic strings using six millisecond pulsars from Data Release 2 (DR2) of the European Pulsar Timing Array (EPTA). We perform a Bayesian analysis considering two models for the network of cosmic string loops, and compare it to a simple power-law model which is expected from the population of supermassive blac… ▽ More We search for a stochastic gravitational wave background (SGWB) generated by a network of cosmic strings using six millisecond pulsars from Data Release 2 (DR2) of the European Pulsar Timing Array (EPTA). We perform a Bayesian analysis considering two models for the network of cosmic string loops, and compare it to a simple power-law model which is expected from the population of supermassive black hole binaries. Our main strong assumption is that the previously reported common red noise process is a SGWB. We find that the one-parameter cosmic string model is slightly favored over a power-law model thanks to its simplicity. If we assume a two-component stochastic signal in the data (supermassive black hole binary population and the signal from cosmic strings), we get a $95\%$ upper limit on the string tension of $\log_{10}(Gμ) < -9.9$ ($-10.5$) for the two cosmic string models we consider. In extended two-parameter string models, we were unable to constrain the number of kinks. We test two approximate and fast Bayesian data analysis methods against the most rigorous analysis and find consistent results. These two fast and efficient methods are applicable to all SGWBs, independent of their source, and will be crucial for analysis of extended data sets. △ Less

Submitted 3 May, 2024; v1 submitted 21 June, 2023; originally announced June 2023.

Comments: 14 pages, 6 figures; typo corrected in (5)

Journal ref: Phys. Rev. D 108 (2023), 123527

New constraints on the kinematic, relativistic and evolutionary properties of the PSR J1757$-$1854 double neutron star system

Authors: A. D. Cameron, M. Bailes, D. J. Champion, P. C. C. Freire, M. Kramer, M. A. McLaughlin, C. Ng, A. Possenti, A. Ridolfi, T. M. Tauris, H. M. Wahl, N. Wex

Abstract: PSR J1757$-$1854 is one of the most relativistic double neutron star binary systems known in our Galaxy, with an orbital period of $P_\text{b}=4.4\,\text{hr}$ and an orbital eccentricity of $e=0.61$. As such, it has promised to be an outstanding laboratory for conducting tests of relativistic gravity. We present the results of a 6-yr campaign with the 100-m Green Bank and 64-m Parkes radio telesco… ▽ More PSR J1757$-$1854 is one of the most relativistic double neutron star binary systems known in our Galaxy, with an orbital period of $P_\text{b}=4.4\,\text{hr}$ and an orbital eccentricity of $e=0.61$. As such, it has promised to be an outstanding laboratory for conducting tests of relativistic gravity. We present the results of a 6-yr campaign with the 100-m Green Bank and 64-m Parkes radio telescopes, designed to capitalise on this potential. We identify secular changes in the profile morphology and polarisation of PSR J1757$-$1854, confirming the presence of geodetic precession and allowing the constraint of viewing geometry solutions consistent with General Relativity. We also update PSR J1757$-$1854's timing, including new constraints of the pulsar's proper motion, post-Keplerian parameters and component masses. We conclude that the radiative test of gravity provided by PSR J1757$-$1854 is fundamentally limited to a precision of 0.3 per cent due to the pulsar's unknown distance. A search for pulsations from the companion neutron star is also described, with negative results. We provide an updated evaluation of the system's evolutionary history, finding strong support for a large kick velocity of $w\ge280\,\text{km s}^{-1}$ following the second progenitor supernova. Finally, we reassess PSR J1757$-$1854's potential to provide new relativistic tests of gravity. We conclude that a 3-$σ$ constraint of the change in the projected semi-major axis ($\dot{x}$) associated with Lense-Thirring precession is expected no earlier than 2031. Meanwhile, we anticipate a 3-$σ$ measurement of the relativistic orbital deformation parameter $δ_θ$ as soon as 2026. △ Less

Submitted 24 May, 2023; originally announced May 2023.

Comments: 23 pages, 16 figures, 7 tables

Gravitational signal propagation in the Double Pulsar studied with the MeerKAT telescope

Authors: H. Hu, M. Kramer, D. J. Champion, N. Wex, A. Parthasarathy, T. T. Pennucci, N. K. Porayko, W. van Straten, V. Venkatraman Krishnan, M. Burgay, P. C. C. Freire, R. N. Manchester, A. Possenti, I. H. Stairs, M. Bailes, S. Buchner, A. D. Cameron, F. Camilo, M. Serylak

Abstract: The Double Pulsar, PSR J0737-3039A/B, has offered a wealth of gravitational experiments in the strong-field regime, all of which GR has passed with flying colours. In particular, among current gravity experiments that test photon propagation, the Double Pulsar probes the strongest spacetime curvature. Observations with MeerKAT and, in future, the SKA can greatly improve the accuracy of current tes… ▽ More The Double Pulsar, PSR J0737-3039A/B, has offered a wealth of gravitational experiments in the strong-field regime, all of which GR has passed with flying colours. In particular, among current gravity experiments that test photon propagation, the Double Pulsar probes the strongest spacetime curvature. Observations with MeerKAT and, in future, the SKA can greatly improve the accuracy of current tests and facilitate tests of NLO contributions in both orbital motion and signal propagation. We present our timing analysis of new observations of PSR J0737-3039A, made using the MeerKAT telescope over the last 3 years. The increased timing precision offered by MeerKAT yields a 2 times better measurement of Shapiro delay parameter s and improved mass measurements compared to previous studies. In addition, our results provide an independent confirmation of the NLO signal propagation effects and already surpass the previous measurement from 16-yr data by a factor of 1.65. These effects include the retardation effect due to the movement of B and the deflection of the signal by the gravitational field of B. We also investigate novel effects which are expected. For instance, we search for potential profile variations near superior conjunctions caused by shifts of the line-of-sight due to latitudinal signal deflection and find insignificant evidence with our current data. With simulations, we find that the latitudinal deflection delay is unlikely to be measured with timing because of its correlation with Shapiro delay. Furthermore, although it is currently not possible to detect the expected lensing correction to the Shapiro delay, our simulations suggest that this effect may be measured with the full SKA. Finally, we provide an improved analytical description for the signal propagation in the Double Pulsar system that meets the timing precision expected from future instruments such as the full SKA. △ Less

Submitted 23 September, 2022; originally announced September 2022.

Comments: 13 pages, 13 figures, accepted for publication in A&A

Journal ref: A&A 667, A149 (2022)

Robust parameter estimation from pulsar timing data

Authors: A. Samajdar, G. Shaifullah, A. Sesana, J. Antoniadis, M. Burgay, D. J. Champion, S. Chen, M. Kramer, J. W. McKee, M. B. Mickaliger, E. Van der Wateren

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 ver… ▽ More 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. △ Less

Submitted 9 May, 2022; originally announced May 2022.

Comments: 8 pages, 5 figures, 4 tables

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.

Comments: 21 pages 16 figures, published in MNRAS (replaced earlier submission after small changes added in proofs)

Constraining the dense matter equation-of-state with radio pulsars

Authors: Huanchen Hu, Michael Kramer, Norbert Wex, David J. Champion, Marcel S. Kehl

Abstract: Radio pulsars provide some of the most important constraints for our understanding of matter at supranuclear densities. So far, these constraints are mostly given by precision mass measurements of neutron stars (NS). By combining single measurements of the two most massive pulsars, J0348$+$0432 and J0740$+$6620, the resulting lower limit of 1.98 $M_\odot$ (99% confidence) of the maximum NS mass, e… ▽ More Radio pulsars provide some of the most important constraints for our understanding of matter at supranuclear densities. So far, these constraints are mostly given by precision mass measurements of neutron stars (NS). By combining single measurements of the two most massive pulsars, J0348$+$0432 and J0740$+$6620, the resulting lower limit of 1.98 $M_\odot$ (99% confidence) of the maximum NS mass, excludes a large number of equations of state (EOSs). Further EOS constraints, complementary to other methods, are likely to come from the measurement of the moment of inertia (MOI) of binary pulsars in relativistic orbits. The Double Pulsar, PSR J0737$-$3039A/B, is the most promising system for the first measurement of the MOI via pulsar timing. Reviewing this method, based in particular on the first MeerKAT observations of the Double Pulsar, we provide well-founded projections into the future by simulating timing observations with MeerKAT and the SKA. For the first time, we account for the spin-down mass loss in the analysis. Our results suggest that an MOI measurement with 11% accuracy (68% confidence) is possible by 2030. If by 2030 the EOS is sufficiently well known, however, we find that the Double Pulsar will allow for a 7% test of Lense-Thirring precession, or alternatively provide a $\sim3σ$-measurement of the next-to-leading order gravitational wave damping in GR. Finally, we demonstrate that potential new discoveries of double NS systems with orbital periods shorter than that of the Double Pulsar promise significant improvements in these measurements and the constraints on NS matter. △ Less

Submitted 15 July, 2020; originally announced July 2020.

Comments: 13 pages, 8 figures. Accepted by MNRAS

The High Time Resolution Universe Pulsar Survey - XIII. PSR J1757-1854, the most accelerated binary pulsar

Authors: A. D. Cameron, D. J. Champion, M. Kramer, M. Bailes, E. D. Barr, C. G. Bassa, S. Bhandari, N. D. R. Bhat, M. Burgay, S. Burke-Spolaor, R. P. Eatough, C. M. L. Flynn, P. C. C. Freire, A. Jameson, S. Johnston, R. Karuppusamy, M. J. Keith, L. Levin, D. R. Lorimer, A. G. Lyne, M. A. McLaughlin, C. Ng, E. Petroff, A. Possenti, A. Ridolfi , et al. (5 additional authors not shown)

Abstract: We report the discovery of PSR J1757$-$1854, a 21.5-ms pulsar in a highly-eccentric, 4.4-h orbit around a neutron star (NS) companion. PSR J1757$-$1854 exhibits some of the most extreme relativistic parameters of any known pulsar, including the strongest relativistic effects due to gravitational-wave (GW) damping, with a merger time of 76 Myr. Following a 1.6-yr timing campaign, we have measured f… ▽ More We report the discovery of PSR J1757$-$1854, a 21.5-ms pulsar in a highly-eccentric, 4.4-h orbit around a neutron star (NS) companion. PSR J1757$-$1854 exhibits some of the most extreme relativistic parameters of any known pulsar, including the strongest relativistic effects due to gravitational-wave (GW) damping, with a merger time of 76 Myr. Following a 1.6-yr timing campaign, we have measured five post-Keplerian (PK) parameters, yielding the two component masses ($m_\text{p}=1.3384(9)\,\text{M}_\odot$ and $m_\text{c}=1.3946(9)\,\text{M}_\odot$) plus three tests of general relativity (GR), which the theory passes. The larger mass of the NS companion provides important clues regarding the binary formation of PSR J1757$-$1854. With simulations suggesting 3-$σ$ measurements of both the contribution of Lense-Thirring precession to the rate of change of the semi-major axis and the relativistic deformation of the orbit within $\sim7-9$ years, PSR J1757$-$1854 stands out as a unique laboratory for new tests of gravitational theories. △ Less

Submitted 14 January, 2018; v1 submitted 21 November, 2017; originally announced November 2017.

Comments: 6 pages, 2 figures, 2 tables

European Pulsar Timing Array Limits on Continuous Gravitational Waves from Individual Supermassive Black Hole Binaries

Authors: Stanislav Babak, Antoine Petiteau, Alberto Sesana, Patrick Brem, Pablo A. Rosado, Stephen R. Taylor, Antoine Lassus, Jason W. T. Hessels, Cees G. Bassa, Marta Burgay, R. Nicolas Caballero, David J. Champion, Ismael Cognard, Gregory Desvignes, Jonathan R. Gair, Lucas Guillemot, Gemma H. Janssen, Ramesh Karuppusamy, Michael Kramer, Patrick Lazarus, K. J. Lee, Lindley Lentati, Kuo Liu, Chiara M. F. Mingarelli, Stefan Oslowski , et al. (11 additional authors not shown)

Abstract: We have searched for continuous gravitational wave (CGW) signals produced by individually resolvable, circular supermassive black hole binaries (SMBHBs) in the latest EPTA dataset, which consists of ultra-precise timing data on 41 millisecond pulsars. We develop frequentist and Bayesian detection algorithms to search both for monochromatic and frequency-evolving systems. None of the adopted algori… ▽ More We have searched for continuous gravitational wave (CGW) signals produced by individually resolvable, circular supermassive black hole binaries (SMBHBs) in the latest EPTA dataset, which consists of ultra-precise timing data on 41 millisecond pulsars. We develop frequentist and Bayesian detection algorithms to search both for monochromatic and frequency-evolving systems. None of the adopted algorithms show evidence for the presence of such a CGW signal, indicating that the data are best described by pulsar and radiometer noise only. Depending on the adopted detection algorithm, the 95\% upper limit on the sky-averaged strain amplitude lies in the range $6\times 10^{-15}<A<1.5\times10^{-14}$ at $5{\rm nHz}<f<7{\rm nHz}$. This limit varies by a factor of five, depending on the assumed source position, and the most constraining limit is achieved towards the positions of the most sensitive pulsars in the timing array. The most robust upper limit -- obtained via a full Bayesian analysis searching simultaneously over the signal and pulsar noise on the subset of ours six best pulsars -- is $A\approx10^{-14}$. These limits, the most stringent to date at $f<10{\rm nHz}$, exclude the presence of sub-centiparsec binaries with chirp mass $\cal{M}_c>10^9$M$_\odot$ out to a distance of about 25Mpc, and with $\cal{M}_c>10^{10}$M$_\odot$ out to a distance of about 1Gpc ($z\approx0.2$). We show that state-of-the-art SMBHB population models predict $<1\%$ probability of detecting a CGW with the current EPTA dataset, consistent with the reported non-detection. We stress, however, that PTA limits on individual CGW have improved by almost an order of magnitude in the last five years. The continuing advances in pulsar timing data acquisition and analysis techniques will allow for strong astrophysical constraints on the population of nearby SMBHBs in the coming years. △ Less

Submitted 11 September, 2015; v1 submitted 7 September, 2015; originally announced September 2015.

Comments: 16 pages, 11 figures, accepted for publication in MNRAS

Limits on anisotropy in the nanohertz stochastic gravitational-wave background

Authors: S. R. Taylor, C. M. F. Mingarelli, J. R. Gair, A. Sesana, G. Theureau, S. Babak, C. G. Bassa, P. Brem, M. Burgay, R. N. Caballero, D. J. Champion, I. Cognard, G. Desvignes, L. Guillemot, J. W. T. Hessels, G. H. Janssen, R. Karuppusamy, M. Kramer, A. Lassus, P. Lazarus, L. Lentati, K. Liu, S. Osłowski, D. Perrodin, A. Petiteau , et al. (10 additional authors not shown)

Abstract: The paucity of observed supermassive black hole binaries (SMBHBs) may imply that the gravitational wave background (GWB) from this population is anisotropic, rendering existing analyses sub-optimal. We present the first constraints on the angular distribution of a nanohertz stochastic GWB from circular, inspiral-driven SMBHBs using the $2015$ European Pulsar Timing Array data [Desvignes et al. (in… ▽ More The paucity of observed supermassive black hole binaries (SMBHBs) may imply that the gravitational wave background (GWB) from this population is anisotropic, rendering existing analyses sub-optimal. We present the first constraints on the angular distribution of a nanohertz stochastic GWB from circular, inspiral-driven SMBHBs using the $2015$ European Pulsar Timing Array data [Desvignes et al. (in prep.)]. Our analysis of the GWB in the $\sim 2 - 90$ nHz band shows consistency with isotropy, with the strain amplitude in $l>0$ spherical harmonic multipoles $\lesssim 40\%$ of the monopole value. We expect that these more general techniques will become standard tools to probe the angular distribution of source populations. △ Less

Submitted 29 June, 2015; originally announced June 2015.

Comments: 6 pages, 2 figures, 1 table. Accepted for publication in Physical Review Letters

A new limit on local Lorentz invariance violation of gravity from solitary pulsars

Authors: Lijing Shao, R. Nicolas Caballero, Michael Kramer, Norbert Wex, David J. Champion, Axel Jessner

Abstract: Gravitational preferred frame effects are generally predicted by alternative theories that exhibit an isotropic violation of local Lorentz invariance of gravity. They are described by three parameters in the parametrized post-Newtonian formalism. One of their strong-field generalizations, $\hat α_2$, induces a precession of a pulsar's spin around its movement direction with respect to the preferre… ▽ More Gravitational preferred frame effects are generally predicted by alternative theories that exhibit an isotropic violation of local Lorentz invariance of gravity. They are described by three parameters in the parametrized post-Newtonian formalism. One of their strong-field generalizations, $\hat α_2$, induces a precession of a pulsar's spin around its movement direction with respect to the preferred frame. We constrain $\hat α_2$ by using the non-detection of such a precession using the characteristics of the pulse profile. In our analysis we use a large number of observations from the 100-m Effelsberg radio telescope, which cover a time span of approximately 15 years. By combining data from two solitary millisecond pulsars, PSRs B1937+21 and J1744-1134, we get a limit of $|\hat α_2| < 1.6 \times 10^{-9}$ at 95% confidence level, which is more than two orders of magnitude better than its best weak-field counterpart from the Solar system. △ Less

Submitted 9 July, 2013; originally announced July 2013.

Comments: 25 pages, 9 figures; accepted by Classical and Quantum Gravity

Journal ref: Class. Quantum Grav. 30 (2013) 165019

Discovery of Nine Gamma-Ray Pulsars in Fermi-LAT Data Using a New Blind Search Method

Authors: H. J. Pletsch, L. Guillemot, B. Allen, M. Kramer, C. Aulbert, H. Fehrmann, P. S. Ray, E. D. Barr, A. Belfiore, F. Camilo, P. A. Caraveo, O. Celik, D. J. Champion, M. Dormody, R. P. Eatough, E. C. Ferrara, P. C. C. Freire, J. W. T. Hessels, M. Keith, M. Kerr, A. de Luca, A. G. Lyne, M. Marelli, M. A. McLaughlin, D. Parent , et al. (6 additional authors not shown)

Abstract: We report the discovery of nine previously unknown gamma-ray pulsars in a blind search of data from the Fermi Large Area Telescope (LAT). The pulsars were found with a novel hierarchical search method originally developed for detecting continuous gravitational waves from rapidly rotating neutron stars. Designed to find isolated pulsars spinning at up to kHz frequencies, the new method is computati… ▽ More We report the discovery of nine previously unknown gamma-ray pulsars in a blind search of data from the Fermi Large Area Telescope (LAT). The pulsars were found with a novel hierarchical search method originally developed for detecting continuous gravitational waves from rapidly rotating neutron stars. Designed to find isolated pulsars spinning at up to kHz frequencies, the new method is computationally efficient, and incorporates several advances, including a metric-based gridding of the search parameter space (frequency, frequency derivative and sky location) and the use of photon probability weights. The nine pulsars have spin frequencies between 3 and 12 Hz, and characteristic ages ranging from 17 kyr to 3 Myr. Two of them, PSRs J1803-2149 and J2111+4606, are young and energetic Galactic-plane pulsars (spin-down power above 6e35 erg/s and ages below 100 kyr). The seven remaining pulsars, PSRs J0106+4855, J0622+3749, J1620-4927, J1746-3239, J2028+3332, J2030+4415, J2139+4716, are older and less energetic; two of them are located at higher Galactic latitudes (|b| > 10 deg). PSR J0106+4855 has the largest characteristic age (3 Myr) and the smallest surface magnetic field (2e11 G) of all LAT blind-search pulsars. PSR J2139+4716 has the lowest spin-down power (3e33 erg/s) among all non-recycled gamma-ray pulsars ever found. Despite extensive multi-frequency observations, only PSR J0106+4855 has detectable pulsations in the radio band. The other eight pulsars belong to the increasing population of radio-quiet gamma-ray pulsars. △ Less

Submitted 2 November, 2011; originally announced November 2011.

Comments: 20 pages, 13 figures. Accepted for publication in the Astrophysical Journal

Report number: AEI-2011-080

Journal ref: The Astrophysical Journal, 744, 105 (2012)