Showing 1–50 of 125 results for author: Keith, M J

Measuring glitch recoveries and braking indices with Bayesian model selection

Authors: Yang Liu, Michael J. Keith, Danai Antonopoulou, Patrick Weltevrede, Benjamin Shaw, Benjamin W. Stappers, Andrew G. Lyne, Mitchell B. Mickaliger, Avishek Basu

Abstract: For a selection of 35 pulsars with large spin-up glitches ($Δν/ν\geq10^{-6}$), which are monitored by the Jodrell Bank Observatory, we analyse 157 glitches and their recoveries. All parameters are measured consistently and we choose the best model to describe the post-glitch recovery based on Bayesian evidence. We present updated glitch epochs, sizes, changes of spin down rate, exponentially recov… ▽ More For a selection of 35 pulsars with large spin-up glitches ($Δν/ν\geq10^{-6}$), which are monitored by the Jodrell Bank Observatory, we analyse 157 glitches and their recoveries. All parameters are measured consistently and we choose the best model to describe the post-glitch recovery based on Bayesian evidence. We present updated glitch epochs, sizes, changes of spin down rate, exponentially recovering components (amplitude and corresponding timescale) when present, as well as pulsars' second frequency derivatives and their glitch associated changes if detected. We discuss the different observed styles of post-glitch recovery as well as some particularly interesting sources. Several correlations are revealed between glitch parameters and pulsar spin parameters, including a very strong correlation between a pulsar's interglitch $|\ddotν|$ and $\dotν$, as well as between the glitch-induced spin-down rate change $Δ\dotν_{\rm p}$ that does not relax exponentially and $\dotν$. We find that the ratio $\left|Δ\dotν_{\mathrm{p}}/\ddotν\right|$ can be used as an estimate of glitch recurrence times, especially for those pulsars for which there are indications of a characteristic glitch size and interglitch waiting time. We calculate the interglitch braking index $n$ and find that pulsars with large glitches typically have $n$ greater than $3$, suggesting that internal torques dominate the rotational evolution between glitches. The external torque, e.g. from electromagnetic dipole radiation, could dominate the observed $\ddotν$ for the youngest pulsars ($\lesssim10^{4}\;\mathrm{yr}$), which may be expected to display $n\sim3$. △ Less

Submitted 13 June, 2024; originally announced June 2024.

Comments: 24 pages, 14 figures, 8 tables. Accepted for publication in MNRAS

The NANOGrav 15 yr Data Set: Chromatic Gaussian Process Noise Models for Six Pulsars

Authors: Bjorn Larsen, Chiara M. F. Mingarelli, Jeffrey S. Hazboun, Aurelien Chalumeau, Deborah C. Good, Joseph Simon, Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Paul R. Brook, H. Thankful Cromartie, Kathryn Crowter, Megan E. DeCesar, Paul B. Demorest, Timothy Dolch, Elizabeth C. Ferrara, William Fiore, Emmanuel Fonseca, Gabriel E. Freedman, Nate Garver-Daniels, Peter A. Gentile, Joseph Glaser, Ross J. Jennings , et al. (39 additional authors not shown)

Abstract: Pulsar timing arrays (PTAs) are designed to detect low-frequency gravitational waves (GWs). GWs induce achromatic signals in PTA data, meaning that the timing delays do not depend on radio-frequency. However, pulse arrival times are also affected by radio-frequency dependent "chromatic" noise from sources such as dispersion measure (DM) and scattering delay variations. Furthermore, the characteriz… ▽ More Pulsar timing arrays (PTAs) are designed to detect low-frequency gravitational waves (GWs). GWs induce achromatic signals in PTA data, meaning that the timing delays do not depend on radio-frequency. However, pulse arrival times are also affected by radio-frequency dependent "chromatic" noise from sources such as dispersion measure (DM) and scattering delay variations. Furthermore, the characterization of GW signals may be influenced by the choice of chromatic noise model for each pulsar. To better understand this effect, we assess if and how different chromatic noise models affect achromatic noise properties in each pulsar. The models we compare include existing DM models used by NANOGrav and noise models used for the European PTA Data Release 2 (EPTA DR2). We perform this comparison using a subsample of six pulsars from the NANOGrav 15 yr data set, selecting the same six pulsars as from the EPTA DR2 six-pulsar dataset. We find that the choice of chromatic noise model noticeably affects the achromatic noise properties of several pulsars. This is most dramatic for PSR J1713+0747, where the amplitude of its achromatic red noise lowers from $\log_{10}A_{\text{RN}} = -14.1^{+0.1}_{-0.1}$ to $-14.7^{+0.3}_{-0.5}$, and the spectral index broadens from $γ_{\text{RN}} = 2.6^{+0.5}_{-0.4}$ to $γ_{\text{RN}} = 3.5^{+1.2}_{-0.9}$. We also compare each pulsar's noise properties with those inferred from the EPTA DR2, using the same models. From the discrepancies, we identify potential areas where the noise models could be improved. These results highlight the potential for custom chromatic noise models to improve PTA sensitivity to GWs. △ Less

Submitted 23 May, 2024; originally announced May 2024.

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

The Thousand-Pulsar-Array programme on MeerKAT XIII: Timing, flux density, rotation measure and dispersion measure timeseries of 597 pulsars

Authors: M. J. Keith, S. Johnston, A. Karastergiou, P. Weltevrede, M. E. Lower, A. Basu, B. Posselt, L. S. Oswald, A. Parthasarathy, A. D. Cameron, M. Serylak, S. Buchner

Abstract: We report here on the timing of 597 pulsars over the last four years with the MeerKAT telescope. We provide Times-of-Arrival, pulsar ephemeris files and per-epoch measurements of the flux density, dispersion measure (DM) and rotation measure (RM) for each pulsar. In addition we use a Gaussian process to model the timing residuals to measure the spin frequency derivative at each epoch. We also repo… ▽ More We report here on the timing of 597 pulsars over the last four years with the MeerKAT telescope. We provide Times-of-Arrival, pulsar ephemeris files and per-epoch measurements of the flux density, dispersion measure (DM) and rotation measure (RM) for each pulsar. In addition we use a Gaussian process to model the timing residuals to measure the spin frequency derivative at each epoch. We also report the detection of 11 glitches in 9 individual pulsars. We find significant DM and RM variations in 87 and 76 pulsars respectively. We find that the DM variations scale approximately linearly with DM, which is broadly in agreement with models of the ionised interstellar medium. The observed RM variations seem largely independent of DM, which may suggest that the RM variations are dominated by variations in the interstellar magnetic field on the line of sight, rather than varying electron density. We also find that normal pulsars have around 5 times greater amplitude of DM variability compared to millisecond pulsars, and surmise that this is due to the known difference in their velocity distributions. △ Less

Submitted 2 April, 2024; originally announced April 2024.

Comments: This is a pre-copyedited, author-produced PDF of an article accepted for publication in MNRAS following peer review. 12 pages 7 figures

The Thousand-Pulsar-Array programme on MeerKAT -- XII. Discovery of long-term pulse profile evolution in 7 young pulsars

Authors: A. Basu, P. Weltevrede, M. J. Keith, S. Johnston, A. Karastergiou, L. S. Oswald, B. Posselt, X. Song, A. D. Cameron

Abstract: A number of pulsars are known to have profile evolution on timescales of months, often correlated with spin-down rate changes. Here, we present the first result from 3 years of monitoring observations from MeerKAT as part of the Thousand Pulsar Array programme. This programme obtains high-fidelity pulse profiles for $\sim$ 500 pulsars, which enabled the detection of subtle changes in seven sources… ▽ More A number of pulsars are known to have profile evolution on timescales of months, often correlated with spin-down rate changes. Here, we present the first result from 3 years of monitoring observations from MeerKAT as part of the Thousand Pulsar Array programme. This programme obtains high-fidelity pulse profiles for $\sim$ 500 pulsars, which enabled the detection of subtle changes in seven sources not previously known to exhibit long-term profile evolution. A 2D Gaussian convolution is used to highlight correlated emission variability in both the pulse phase and observing epoch direction. Simulations show that for one additional source the observed profile variability is likely to originate from stochastic single-pulse shape variability (jitter). We find that it is common for long-term profile variability to be associated with changes in polarization fractions, but not with polarisation position angle (PA) changes. PA changes are expected if emission height changes or precession is responsible for the profile variability. PSR J1741$-$3927 is the only pulsar in our sample that shows correlated PA variability, and this is associated with orthogonal polarization mode activity. For the six %the rest, without correlated PA variability, other pulsars limits on possible emission height changes and impact angle changes are derived. These limits are consistent with the small changes in the total intensity profile shape. None of the sources show detectable spin-down variability correlated with the emission changes, which are thought to be driven by magnetospheric current fluctuations. Therefore the absence of correlated spin-down rate variability allows upper limits to be placed on changes in the magnetospheric charge density. △ Less

Submitted 14 February, 2024; originally announced February 2024.

Comments: 19 pages, 12 figures, 3 Tables. Accepted for publication in MNRAS

A Gaussian-processes approach to fitting for time-variable spherical solar wind in pulsar timing data

Authors: Iuliana C. Niţu, Michael J. Keith, Caterina Tiburzi, Marcus Brüggen, David J. Champion, Siyuan Chen, Ismaël Cognard, Gregory Desvignes, Ralf-Jürgen Dettmar, Jean-Mathias Grießmeier, Lucas Guillemot, Yanjun Guo, Matthias Hoeft, Huanchen Hu, Jiwoong Jang, Gemma H. Janssen, Jedrzej Jawor, Ramesh Karuppusamy, Evan F. Keane, Michael Kramer, Jörn Künsemöller, Kristen Lackeos, Kuo Liu, Robert A. Main, James W. McKee , et al. (4 additional authors not shown)

Abstract: Propagation effects are one of the main sources of noise in high-precision pulsar timing. For pulsars below an ecliptic latitude of $5^\circ$, the ionised plasma in the solar wind can introduce dispersive delays of order 100 microseconds around solar conjunction at an observing frequency of 300 MHz. A common approach to mitigate this assumes a spherical solar wind with a time-constant amplitude. H… ▽ More Propagation effects are one of the main sources of noise in high-precision pulsar timing. For pulsars below an ecliptic latitude of $5^\circ$, the ionised plasma in the solar wind can introduce dispersive delays of order 100 microseconds around solar conjunction at an observing frequency of 300 MHz. A common approach to mitigate this assumes a spherical solar wind with a time-constant amplitude. However, this has been shown to be insufficient to describe the solar wind. We present a linear, Gaussian-process piecewise Bayesian approach to fit a spherical solar wind of time-variable amplitude, which has been implemented in the pulsar software run_enterprise. Through simulations, we find that the current EPTA+InPTA data combination is not sensitive to such variations; however, solar wind variations will become important in the near future with the addition of new InPTA data and data collected with the low-frequency LOFAR telescope. We also compare our results for different high-precision timing datasets (EPTA+InPTA, PPTA, and LOFAR) of three millisecond pulsars (J0030$+$0451, J1022$+$1001, J2145$-$0450), and find that the solar-wind amplitudes are generally consistent for any individual pulsar, but they can vary from pulsar to pulsar. Finally, we compare our results with those of an independent method on the same LOFAR data of the three millisecond pulsars. We find that differences between the results of the two methods can be mainly attributed to the modelling of dispersion variations in the interstellar medium, rather than the solar wind modelling. △ Less

Submitted 15 January, 2024; originally announced January 2024.

Comments: Accepted for publication in MNRAS

The High Time Resolution Universe Pulsar Survey -- XVIII. The reprocessing of the HTRU-S Low Lat survey around the Galactic centre using a Fast Folding Algorithm pipeline for accelerated pulsars

Authors: J. Wongphechauxsorn, D. J. Champion, M. Bailes, V. Balakrishnan, E. D. Barr, M. C. i Bernadich, N. D. R. Bhat, M. Burgay, A. D. Cameron, W. Chen, C. M. L. Flynn, A. Jameson, S. Johnston, M. J. Keith, M. Kramer, C. Ng, A. Possenti, R. Sengar, R. M. Shannon, B. Stappers, W. van Straten

Abstract: The HTRU-S Low Latitude survey data within 1$^{\circ}$of the Galactic Centre (GC) were searched for pulsars using the Fast Folding Algorithm (FFA). Unlike traditional Fast Fourier Transform (FFT) pipelines, the FFA optimally folds the data for all possible periods over a given range, which is particularly advantageous for pulsars with low-duty cycle. For the first time, a search over acceleration… ▽ More The HTRU-S Low Latitude survey data within 1$^{\circ}$of the Galactic Centre (GC) were searched for pulsars using the Fast Folding Algorithm (FFA). Unlike traditional Fast Fourier Transform (FFT) pipelines, the FFA optimally folds the data for all possible periods over a given range, which is particularly advantageous for pulsars with low-duty cycle. For the first time, a search over acceleration was included in the FFA to improve its sensitivity to binary pulsars. The steps in dispersion measure (DM) and acceleration were optimised, resulting in a reduction of the number of trials by 86 per cent. This was achieved over a search period range from 0.6-s to 432-s, i.e. 10 per cent of the observation time (4320s), with a maximum DM of 4000 pc cm$^{-3}$ and an acceleration range of $\pm 128$m s$^{-2}$. The search resulted in the re-detections of four known pulsars, including a pulsar which was missed in previous FFT processing of this survey. This result indicates that the FFA pipeline is more sensitive than the FFT pipeline used in the previous processing of the survey within our parameter range. Additionally, we discovered a 1.89-s pulsar, PSR J1746-2829, with a large DM, located~0.5 from the GC. Follow-up observations revealed that this pulsar has a relatively flat spectrum($α=-0.9\pm0.1$) and has a period derivative of $\sim1.3\times10^{-12}$ s s$^{-1}$, implying a surface magnetic field of $\sim5.2\times10^{13}$ G and a characteristic age of $\sim23000$ yr. While the period, spectral index, and surface magnetic field strength are similar to many radio magnetars, other characteristics such as high linear polarization are absent. △ Less

Submitted 21 October, 2023; originally announced October 2023.

Comments: 12 pages, 6 figures, 4 tables, Accepted for publication on Monthly Notices of the Royal Astronomical Society

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 Third Fermi Large Area Telescope Catalog of Gamma-ray Pulsars

Authors: David A. Smith, Philippe Bruel, Colin J. Clark, Lucas Guillemot, Matthew T. Kerr, Paul Ray, Soheila Abdollahi, Marco Ajello, Luca Baldini, Jean Ballet, Matthew Baring, Cees Bassa, Josefa Becerra Gonzalez, Ronaldo Bellazzini, Alessandra Berretta, Bhaswati Bhattacharyya, Elisabetta Bissaldi, Raffaella Bonino, Eugenio Bottacini, Johan Bregeon, Marta Burgay, Toby Burnett, Rob Cameron, Fernando Camilo, Regina Caputo , et al. (134 additional authors not shown)

Abstract: We present 294 pulsars found in GeV data from the Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope. Another 33 millisecond pulsars (MSPs) discovered in deep radio searches of LAT sources will likely reveal pulsations once phase-connected rotation ephemerides are achieved. A further dozen optical and/or X-ray binary systems co-located with LAT sources also likely harbor gamma-ray M… ▽ More We present 294 pulsars found in GeV data from the Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope. Another 33 millisecond pulsars (MSPs) discovered in deep radio searches of LAT sources will likely reveal pulsations once phase-connected rotation ephemerides are achieved. A further dozen optical and/or X-ray binary systems co-located with LAT sources also likely harbor gamma-ray MSPs. This catalog thus reports roughly 340 gamma-ray pulsars and candidates, 10% of all known pulsars, compared to $\leq 11$ known before Fermi. Half of the gamma-ray pulsars are young. Of these, the half that are undetected in radio have a broader Galactic latitude distribution than the young radio-loud pulsars. The others are MSPs, with 6 undetected in radio. Overall, >235 are bright enough above 50 MeV to fit the pulse profile, the energy spectrum, or both. For the common two-peaked profiles, the gamma-ray peak closest to the magnetic pole crossing generally has a softer spectrum. The spectral energy distributions tend to narrow as the spindown power $\dot E$ decreases to its observed minimum near $10^{33}$ erg s$^{-1}$, approaching the shape for synchrotron radiation from monoenergetic electrons. We calculate gamma-ray luminosities when distances are available. Our all-sky gamma-ray sensitivity map is useful for population syntheses. The electronic catalog version provides gamma-ray pulsar ephemerides, properties and fit results to guide and be compared with modeling results. △ Less

Submitted 20 July, 2023; originally announced July 2023.

Comments: 142 pages. Accepted by the Astrophysical Journal Supplement

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

The second data release from the European Pulsar Timing Array II. Customised pulsar noise models for spatially correlated gravitational waves

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. (73 additional authors not shown)

Abstract: The nanohertz gravitational wave background (GWB) is expected to be an aggregate signal of an ensemble of gravitational waves emitted predominantly by a large population of coalescing supermassive black hole binaries in the centres of merging galaxies. Pulsar timing arrays, ensembles of extremely stable pulsars, are the most precise experiments capable of detecting this background. However, the su… ▽ More The nanohertz gravitational wave background (GWB) is expected to be an aggregate signal of an ensemble of gravitational waves emitted predominantly by a large population of coalescing supermassive black hole binaries in the centres of merging galaxies. Pulsar timing arrays, ensembles of extremely stable pulsars, are the most precise experiments capable of detecting this background. However, the subtle imprints that the GWB induces on pulsar timing data are obscured by many sources of noise. These must be carefully characterized to increase the sensitivity to the GWB. In this paper, we present a novel technique to estimate the optimal number of frequency coefficients for modelling achromatic and chromatic noise and perform model selection. We also incorporate a new model to fit for scattering variations in the pulsar timing package temponest and created realistic simulations of the European Pulsar Timing Array (EPTA) datasets that allowed us to test the efficacy of our noise modelling algorithms. We present an in-depth analysis of the noise properties of 25 millisecond pulsars (MSPs) that form the second data release (DR2) of the EPTA and investigate the effect of incorporating low-frequency data from the Indian PTA collaboration. We use enterprise and temponest packages to compare noise models with those reported with the EPTA DR1. We find that, while in some pulsars we can successfully disentangle chromatic from achromatic noise owing to the wider frequency coverage in DR2, in others the noise models evolve in a more complicated way. We also find evidence of long-term scattering variations in PSR J1600$-$3053. Through our simulations, we identify intrinsic biases in our current noise analysis techniques and discuss their effect on GWB searches. The results presented here directly help improve sensitivity to the GWB and are already being used as part of global PTA efforts. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 20 pages, 6 figures, 9 tables

Journal ref: A&A 678, A49 (2023)

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

The second data release from the European Pulsar Timing Array III. Search for 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. (73 additional authors not shown)

Abstract: We present the results of the search for an isotropic stochastic gravitational wave background (GWB) at nanohertz frequencies using the second data release of the European Pulsar Timing Array (EPTA) for 25 millisecond pulsars and a combination with the first data release of the Indian Pulsar Timing Array (InPTA). We analysed (i) the full 24.7-year EPTA data set, (ii) its 10.3-year subset based on… ▽ More We present the results of the search for an isotropic stochastic gravitational wave background (GWB) at nanohertz frequencies using the second data release of the European Pulsar Timing Array (EPTA) for 25 millisecond pulsars and a combination with the first data release of the Indian Pulsar Timing Array (InPTA). We analysed (i) the full 24.7-year EPTA data set, (ii) its 10.3-year subset based on modern observing systems, (iii) the combination of the full data set with the first data release of the InPTA for ten commonly timed millisecond pulsars, and (iv) the combination of the 10.3-year subset with the InPTA data. These combinations allowed us to probe the contributions of instrumental noise and interstellar propagation effects. With the full data set, we find marginal evidence for a GWB, with a Bayes factor of four and a false alarm probability of $4\%$. With the 10.3-year subset, we report evidence for a GWB, with a Bayes factor of $60$ and a false alarm probability of about $0.1\%$ ($\gtrsim 3σ$ significance). The addition of the InPTA data yields results that are broadly consistent with the EPTA-only data sets, with the benefit of better noise modelling. Analyses were performed with different data processing pipelines to test the consistency of the results from independent software packages. The inferred spectrum from the latest EPTA data from new generation observing systems is rather uncertain and in mild tension with the common signal measured in the full data set. However, if the spectral index is fixed at 13/3, the two data sets give a similar amplitude of ($2.5\pm0.7)\times10^{-15}$ at a reference frequency of $1\,{\rm yr}^{-1}$. By continuing our detection efforts as part of the International Pulsar Timing Array (IPTA), we expect to be able to improve the measurement of spatial correlations and better characterise this signal in the coming years. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 21 pages, 14 figures, 4 appendix figures, accepted for publication in A&A

Journal ref: A&A 678, A50 (2023)

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

Impact of quasi-periodic and steep-spectrum timing noise on the measurement of pulsar timing parameters

Authors: Michael J. Keith, Iuliana C. Niţu

Abstract: Timing noise in pulsars is often modelled with a Fourier-basis Gaussian process that follows a power law with periodic boundary conditions on the observation time, $T_\mathrm{span}$. However the actual noise processes can extend well below $1/T_\mathrm{span}$, and many pulsars are known to exhibit quasi-periodic timing noise. In this paper we investigate several adaptions that try to account for t… ▽ More Timing noise in pulsars is often modelled with a Fourier-basis Gaussian process that follows a power law with periodic boundary conditions on the observation time, $T_\mathrm{span}$. However the actual noise processes can extend well below $1/T_\mathrm{span}$, and many pulsars are known to exhibit quasi-periodic timing noise. In this paper we investigate several adaptions that try to account for these differences between the observed behaviour and the simple power-law model. Firstly, we propose to include an additional term that models the quasi-periodic spin-down variations known to be present in many pulsars. Secondly, we show that a Fourier basis of $1/2T_\mathrm{span}$ can be more suited for estimating long term timing parameters such as the spin frequency second derivative (F2), and is required when the exponent of the power spectrum is greater than ~4. We also implement a Bayesian version of the generalised least squares `Cholesky' method which has different limitations at low frequency, but find that there is little advantage over Fourier-basis methods. We apply our quasi-periodic spin down model to a sample of pulsars with known spin-down variations and show that this improves parameter estimation of F2 and proper motion for the most pathological cases, but in general the results are consistent with a power-law model. The models are all made available through the run_enterprise software package. △ Less

Submitted 6 June, 2023; originally announced June 2023.

Comments: This article has been accepted for publication in MNRAS Published by Oxford University Press on behalf of the Royal Astronomical Society

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.

Comments: 16 pages, 7 figures. Abstract shortened to adhere to ArXiv limit

Journal ref: A&A 674, A169 (2023)

Searching for Time-Dependent Axion Dark Matter Signals in Pulsars

Authors: R. A. Battye, M. J. Keith, J. I. McDonald, S. Srinivasan, B. W. Stappers, P. Weltevrede

Abstract: Axion dark matter can be converted into photons in the magnetospheres of neutron stars leading to a spectral line centred on the Compton wavelength of the axion. Due to the rotation of the star and the plasma effects in the magnetosphere the signal is predicted to be periodic with significant time variation - a unique smoking gun for axion dark matter. As a proof of principle and to develop the me… ▽ More Axion dark matter can be converted into photons in the magnetospheres of neutron stars leading to a spectral line centred on the Compton wavelength of the axion. Due to the rotation of the star and the plasma effects in the magnetosphere the signal is predicted to be periodic with significant time variation - a unique smoking gun for axion dark matter. As a proof of principle and to develop the methodology, we carry out the first time domain search of the signal using data from PSR J2144$-$3933 taken as part of the MeerTIME project on MeerKAT telescope. We search for specific signal templates using a matched filter technique and discuss when a time-domain analysis (as is typically the case in pulsar observations) gives greater sensitivity to the axion-coupling to photons in comparison to a simple time-averaged total flux study. We do not find any candidate signals and, hence, impose an upper limit on the axion-to-photon coupling of $g_{aγγ}<4\times 10^{-11}\,{\rm GeV}^{-1}$ over the mass range $m_{\rm a}=3.9-4.7\,μ{\rm eV}$ using this data. This limit relies on PSR J2144$-$3933 not being an extremely aligned rotator, as strongly supported by simple arguments based on the observed pulse profile width. We discuss the possibilities of improving this limit using future observations with MeerKAT and also SKA1-mid and the possibility of using other objects. Finally, to evade modelling uncertainties in axion radio signals, we also carry out a generic ``any periodic-signal search" in the data, finding no evidence for an axion signal. △ Less

Submitted 28 November, 2023; v1 submitted 21 March, 2023; originally announced March 2023.

Comments: 21 pages, 14 figures Published in PRD

Journal ref: Physical Review D vol. 108, no. 6, 2023

The Thousand-Pulsar-Array programme on MeerKAT -- VIII. The subpulse modulation of 1198 pulsars

Authors: X. Song, P. Weltevrede, A. Szary, G. Wright, M. J. Keith, A. Basu, S. Johnston, A. Karastergiou, R. A. Main, L. S. Oswald, A. Parthasarathy, B. Posselt, M. Bailes, S. Buchner, B. Hugo, M. Serylak

Abstract: We report on the subpulse modulation properties of 1198 pulsars using the Thousand-Pulsar-Array programme on MeerKAT. About 35% of the analysed pulsars exhibit drifting subpulses which are more pronounced towards the deathline, consistent with previous studies. We estimate that this common phenomenon is detectable in 60% of the overall pulsar population if high quality data were available for all.… ▽ More We report on the subpulse modulation properties of 1198 pulsars using the Thousand-Pulsar-Array programme on MeerKAT. About 35% of the analysed pulsars exhibit drifting subpulses which are more pronounced towards the deathline, consistent with previous studies. We estimate that this common phenomenon is detectable in 60% of the overall pulsar population if high quality data were available for all. This large study reveals the evolution of drifting subpulses across the pulsar population in unprecedented detail. In particular, we find that the modulation period $P_3$ follows a V-shaped evolution with respect to the characteristic age $τ_c$, such that the smallest $P_3$ values, corresponding to the Nyquist period $P_3>\sim2$, are found at $τ_c>\sim10^{7.5}$ yr. The V-shaped evolution can be interpreted and reproduced if young pulsars possess aliased fast intrinsic $P_3$, which monotonically increase, ultimately achieving a slow unaliased $P_3$. Enhancement of irregularities in intrinsic subpulse modulation by aliasing in small $τ_c$ pulsars would explain their observed less well defined $P_3$'s and weaker spectral features. Modelling these results as rotating subbeams, their circulation must slow down as the pulsar evolves. This is the opposite to that expected if circulation is driven by ExB drift. This can be resolved if the observed $P_3$ periodicity is due to a beat between an ExB system and the pulsar period. As a by-product, we identify the correct periods and spin-down rates for 12 pulsars, for which harmonically related values were reported in the literature. △ Less

Submitted 10 January, 2023; originally announced January 2023.

Comments: 18 pages, 13 figures in the main part of the paper. This version includes all online materials from the accepted MNRAS article

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.

Comments: 18 pages, 6 figures

The Thousand-Pulsar-Array program on MeerKAT -- IX. The time-averaged properties of the observed pulsar population

Authors: B. Posselt, A. Karastergiou, S. Johnston, A. Parthasarathy, L. S. Oswald, R. A. Main, A. Basu, M. J. Keith, X. Song, P. Weltevrede, C. Tiburzi, M. Bailes, S. Buchner, M. Geyer, M. Kramer, R. Spiewak, V. Venkatraman Krishnan

Abstract: We present the largest single survey to date of average profiles of radio pulsars, observed and processed using the same telescope and data reduction software. Specifically, we present measurements for 1170 pulsars, observed by the Thousand Pulsar Array (TPA) programme at the 64-dish SARAO MeerKAT radio telescope, in a frequency band from 856 to 1712 MHz. We provide rotation measures (RM), dispers… ▽ More We present the largest single survey to date of average profiles of radio pulsars, observed and processed using the same telescope and data reduction software. Specifically, we present measurements for 1170 pulsars, observed by the Thousand Pulsar Array (TPA) programme at the 64-dish SARAO MeerKAT radio telescope, in a frequency band from 856 to 1712 MHz. We provide rotation measures (RM), dispersion measures, flux densities and polarization properties. The catalogue includes 254 new RMs that substantially increase the total number of known pulsar RMs. Our integration times typically span over 1000 individual rotations per source. We show that the radio (pseudo)luminosity has a strong, shallow dependence on the spin-down energy, proportional to $\dot{E}^{0.15\pm0.04}$, that contradicts some previous proposals of population synthesis studies. In addition, we find a significant correlation between the steepness of the observed flux density spectra and $\dot{E}$, and correlations of the fractional linear polarization with $\dot{E}$, the spectral index, and the pulse width, which we discuss in the context of what is known about pulsar radio emission and how pulsars evolve with time. On the whole, we do not see significant correlations with the estimated surface magnetic field strength, and the correlations with $\dot{E}$ are much stronger than those with the characteristic age. This finding lends support to the suggestion that magnetic dipole braking may not be the dominant factor for the evolution of pulsar rotation over the lifetimes of pulsars. A public data release of the high-fidelity time-averaged pulse profiles in full polarization accompanies our catalogue. △ Less

Submitted 21 November, 2022; originally announced November 2022.

Comments: Accepted for publication in MNRAS. 21 pages, 25 figures, 6 Tables

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.

Comments: 12 pages, 5 figures, 2 tables. Accepted in MNRAS

MeerKAT observations of the reversing drifting subpulses in PSR J1750-3503

Authors: Andrzej Szary, Joeri van Leeuwen, Geoff Wright, Patrick Weltevrede, Crispin H. Agar, Caterina Tiburzi, Yogesh Maan, Michael J. Keith

Abstract: We present an analysis of the subpulse drift in PSR J1750-3503, which is characterized by abrupt transitions of drift direction. As the pulsar does not exhibit other mode changes or clear nulling, it is an ideal candidate system for studying the phenomenon of drift direction change. For $\sim 80\%$ of the time the subpulses are characterized by positive drift - from early to later longitudes - whi… ▽ More We present an analysis of the subpulse drift in PSR J1750-3503, which is characterized by abrupt transitions of drift direction. As the pulsar does not exhibit other mode changes or clear nulling, it is an ideal candidate system for studying the phenomenon of drift direction change. For $\sim 80\%$ of the time the subpulses are characterized by positive drift - from early to later longitudes - while the drift direction is negative in the other $\sim 20\%$. The subpulse separation for single pulses with positive drift, $P_2=(18.8\pm 0.1)^{\circ}$, is higher then for single pulses with negative drift, $P_2=(17.5\pm 0.2)^{\circ}$. When the drift is stable, the measured repetition time of the drift pattern is $P_3^{\rm obs}=(43.5 \pm 0.4) P$, where $P$ is pulsar period. We show that the observed data can be reproduced by a carousel models with subpulse rotation around the magnetic axis using purely dipolar configuration of surface magnetic field. The observed drift characteristics can be modeled assuming that the actual repetition time $P_3<2P$, such that we observe its aliased value. A small variation in $P_3$, of the order of $6\%$ (or less assuming higher alias orders), is enough to reproduce the characteristic drift direction changes we observe. △ Less

Submitted 6 June, 2022; originally announced June 2022.

Comments: 14 pages, 11 figures, accepted for publication in ApJ

Long-term rotational and emission variability of 17 radio pulsars

Authors: B. Shaw, B. W. Stappers, P. Weltevrede, P. R. Brook, A. Karastergiou, C. A. Jordan, M. J. Keith, M. Kramer, A. G. Lyne

Abstract: With the ever-increasing sensitivity and timing baselines of modern radio telescopes, a growing number of pulsars are being shown to exhibit transitions in their rotational and radio emission properties. In many of these cases, the two are correlated with pulsars assuming a unique spin-down rate ($\dotν$) for each of their specific emission states. In this work we revisit 17 radio pulsars previous… ▽ More With the ever-increasing sensitivity and timing baselines of modern radio telescopes, a growing number of pulsars are being shown to exhibit transitions in their rotational and radio emission properties. In many of these cases, the two are correlated with pulsars assuming a unique spin-down rate ($\dotν$) for each of their specific emission states. In this work we revisit 17 radio pulsars previously shown to exhibit spin-down rate variations. Using a Gaussian process regression (GPR) method to model the timing residuals and the evolution of the profile shape, we confirm the transitions already observed and reveal new transitions in 8 years of extended monitoring with greater time resolution and enhanced observing bandwidth. We confirm that 7 of these sources show emission-correlated $\dotν$ transitions ($Δ\dotν$) and we characterise this correlation for one additional pulsar, PSR B1642$-$03. We demonstrate that GPR is able to reveal extremely subtle profile variations given sufficient data quality. We also corroborate the dependence of $Δ\dotν$ amplitude on $\dotν$ and pulsar characteristic age. Linking $Δ\dotν$ to changes in the global magnetospheric charge density $Δρ$, we speculate that $\dotν$ transitions associated with large $Δρ$ values may be exhibiting detectable profile changes with improved data quality, in cases where they have not previously been observed. △ Less

Submitted 25 April, 2022; v1 submitted 22 April, 2022; originally announced April 2022.

Comments: 20 pages, 16 figures, 3 tables

The High Time Resolution Universe Pulsar Survey -- XVII. PSR J1325-6253, a low eccentricity double neutron star system from an ultra-stripped supernova

Authors: R. Sengar, V. Balakrishnan, S. Stevenson, M. Bailes, E. D. Barr, N. D. R. Bhat, M. Burgay, M. C. i Bernadich, A. D. Cameron, D. J. Champion, W. Chen, C. M. L. Flynn, A. Jameson, S. Johnston, M. J. Keith, M. Kramer, V. Morello, C. Ng, A. Possenti, B. Stappers, R. M. Shannon, W. van Straten, J. Wongphechauxsorn

Abstract: The observable population of double neutron star (DNS) systems in the Milky Way allow us to understand the nature of supernovae and binary stellar evolution. Until now, all DNS systems in wide orbits ($ P_{\textrm{orb}}>$ 1~day) have been found to have orbital eccentricities, $e > 0.1$. In this paper, we report the discovery of pulsar PSR J1325$-$6253: a DNS system in a 1.81 day orbit with a surpr… ▽ More The observable population of double neutron star (DNS) systems in the Milky Way allow us to understand the nature of supernovae and binary stellar evolution. Until now, all DNS systems in wide orbits ($ P_{\textrm{orb}}>$ 1~day) have been found to have orbital eccentricities, $e > 0.1$. In this paper, we report the discovery of pulsar PSR J1325$-$6253: a DNS system in a 1.81 day orbit with a surprisingly low eccentricity of just $e = 0.064$. Through 1.4 yr of dedicated timing with the Parkes radio telescope we have been able to measure its rate of advance of periastron, $\dotω=0.138 \pm 0.002$ $\rm deg$ $\rm yr^{-1}$. If this induced $\dotω$ is solely due to general relativity then the total mass of the system is, $M_{\rm sys} = 2.57 \pm 0.06$ M$_{\odot}$. Assuming an edge-on orbit the minimum companion mass is constrained to be $M_\mathrm{c,min}>0.98$ M$_{\odot}$ which implies the pulsar mass is $M_\mathrm{p,max}<1.59 $ M$_{\odot}$. Its location in the $P$-$\dot{P}$ diagram suggests that, like other DNS systems, PSR J1325$-$6253 is a recycled pulsar and if its mass is similar to the known examples ($>1.3$ M$_\odot$), then the companion neutron star is probably less than $\sim1.25$ M$_\odot$ and the system is inclined at about $50^{\circ}$-$60^{\circ}$. The low eccentricity along with the wide orbit of the system strongly favours a formation scenario involving an ultra-stripped supernova explosion. △ Less

Submitted 26 May, 2022; v1 submitted 14 April, 2022; originally announced April 2022.

Comments: 12 pages, 7 figures, 1 table, accepted for publication in MNRAS

A Gamma-ray Pulsar Timing Array Constrains the Nanohertz Gravitational Wave Background

Authors: M. Ajello, W. B. Atwood, L. Baldini, J. Ballet, G. Barbiellini, D. Bastieri, R. Bellazzini, A. Berretta, B. Bhattacharyya, E. Bissaldi, R. D. Blandford, E. Bloom, R. Bonino, P. Bruel, R. Buehler, E. Burns, S. Buson, R. A. Cameron, P. A. Caraveo, E. Cavazzuti, N. Cibrario, S. Ciprini, C. J. Clark, I. Cognard, J. Coronado-Blázquez , et al. (107 additional authors not shown)

Abstract: After large galaxies merge, their central supermassive black holes are expected to form binary systems whose orbital motion generates a gravitational wave background (GWB) at nanohertz frequencies. Searches for this background utilize pulsar timing arrays, which perform long-term monitoring of millisecond pulsars (MSPs) at radio wavelengths. We use 12.5 years of Fermi Large Area Telescope data to… ▽ More After large galaxies merge, their central supermassive black holes are expected to form binary systems whose orbital motion generates a gravitational wave background (GWB) at nanohertz frequencies. Searches for this background utilize pulsar timing arrays, which perform long-term monitoring of millisecond pulsars (MSPs) at radio wavelengths. We use 12.5 years of Fermi Large Area Telescope data to form a gamma-ray pulsar timing array. Results from 35 bright gamma-ray pulsars place a 95\% credible limit on the GWB characteristic strain of $1.0\times10^{-14}$ at 1 yr$^{-1}$, which scales as the observing time span $t_{\mathrm{obs}}^{-13/6}$. This direct measurement provides an independent probe of the GWB while offering a check on radio noise models. △ Less

Submitted 11 April, 2022; originally announced April 2022.

Comments: 3 figures in the main text. 3 figures and 8 tables are in the supplementary material

Discoveries and Timing of Pulsars in NGC 6440

Authors: L. Vleeschower, B. W. Stappers, M. Bailes, E. D. Barr, M. Kramer, S. Ransom, A. Ridolfi, V. Venkatraman Krishnan, A. Possenti, M. J. Keith, M. Burgay, P. C. C. Freire, R. Spiewak, D. J. Champion, M. C. Bezuidenhout, I. C. Niţu, W. Chen, A. Parthasarathy, M. E. DeCesar, S. Buchner, I. H. Stairs, J. W. T. Hessels

Abstract: Using the MeerKAT radio telescope, a series of observations have been conducted to time the known pulsars and search for new pulsars in the globular cluster NGC 6440. As a result, two pulsars have been discovered, NGC 6440G and NGC 6440H, one of which is isolated and the other a non-eclipsing (at frequencies above 962 MHz) "Black Widow", with a very low mass companion (M$_{\rm c}$ > 0.006 M… ▽ More Using the MeerKAT radio telescope, a series of observations have been conducted to time the known pulsars and search for new pulsars in the globular cluster NGC 6440. As a result, two pulsars have been discovered, NGC 6440G and NGC 6440H, one of which is isolated and the other a non-eclipsing (at frequencies above 962 MHz) "Black Widow", with a very low mass companion (M$_{\rm c}$ > 0.006 M$_{\odot}$). It joins the other binary pulsars discovered so far in this cluster which all have low companion masses (M$_{\rm c}$ < 0.30 M$_{\odot}$). We present the results of long-term timing solutions obtained using data from both Green Bank and MeerKAT telescopes for these two new pulsars and an analysis of the pulsars NGC 6440C and NGC 6440D. For the isolated pulsar NGC 6440C, we searched for planets using a Markov Chain Monte Carlo technique. We find evidence for significant unmodelled variations but they cannot be well modelled as planets nor as part of a power-law red-noise process. Studies of the eclipses of the "Redback" pulsar NGC 6440D at two different frequency bands reveal a frequency dependence with longer and asymmetric eclipses at lower frequencies (962-1283 MHz). △ Less

Submitted 1 April, 2022; originally announced April 2022.

Comments: 14 pages, 10 figures, 4 tables. Accepted for publication in MNRAS

Long-term scintillation studies of EPTA pulsars. I. Observations and basic results

Authors: Yulan Liu, Joris P. W. Verbiest, Robert A. Main, Ziwei Wu, Krishnakumar Moochickal Ambalappat, David J. Champion, Ismaël Cognard, Lucas Guillemot, Madhuri Gaikwad, Gemma H. Janssen, Michael Kramer, Michael J. Keith, Ramesh Karuppusamy, Lars Künkel, Kuo Liu, James W. McKee, Mitchell B. Mickaliger, Ben W. Stappers, Golam. M. Shaifullah, Gilles Theureau

Abstract: Interstellar scintillation analysis of pulsars allows us to probe the small-scale distribution and inhomogeneities of the ionized interstellar medium. Our priority is to present the data set and the basic measurements of scintillation parameters of pulsars employing long-term scintillation observations carried out from 2011 January to 2020 August by the European Pulsar Timing Array radio telesco… ▽ More Interstellar scintillation analysis of pulsars allows us to probe the small-scale distribution and inhomogeneities of the ionized interstellar medium. Our priority is to present the data set and the basic measurements of scintillation parameters of pulsars employing long-term scintillation observations carried out from 2011 January to 2020 August by the European Pulsar Timing Array radio telescopes in the 21-cm and 11-cm bands. Additionally, we aim to identify future possible lines of study using this long-term scintillation dataset. We present the long-term time series of $ν_{\rm d}$ and $τ_{\rm d}$ for 13 pulsars. Sanity-checks and comparisons indicate that the scintillation parameters of our work and previously published works are mostly consistent. For two pulsars, PSRs~J1857+0943 and J1939+2134, we were able to obtain measurements of the $ν_{\rm d}$ at both bands, which allows us to derive the time series of frequency scaling indices with a mean and a standard deviation of 2.82$\pm$1.95 and 3.18$\pm$0.60, respectively. We found some interesting features which will be studied in more detail in subsequent papers in this series: (i) in the time series of PSR~J1939+2134, where the scintillation bandwidth sharply increases or decreases associated with a sharp change of dispersion measure; (ii) PSR~J0613$-$0200 and PSR~J0636+5126 show a strong annual variation in the time series of the $τ_{\rm d}$; (iii) PSR~J1939+2134 shows a weak anti-correlation between scintillation timescale and dispersion in WSRT data. △ Less

Submitted 13 April, 2022; v1 submitted 31 March, 2022; originally announced March 2022.

Comments: 18 pages, 6 figures, 4 tables, accepted for publication in A&A

Journal ref: A&A 664, A116 (2022)

A search for planetary companions around 800 pulsars from the Jodrell Bank pulsar timing programme

Authors: Iuliana C. Niţu, Michael J. Keith, Ben W. Stappers, Andrew G. Lyne, Mitchell B. Mickaliger

Abstract: We have searched for planetary companions around 800 pulsars monitored at the Jodrell Bank Observatory, with both circular and eccentric orbits of periods between $20$ days and $17$ years and inclination-dependent planetary masses from $10^{-4}$ to $100\,\mathrm{M}_{\oplus}$. Using a Bayesian framework, we simultaneously model pulsar timing parameters and a stationary noise process with a power-la… ▽ More We have searched for planetary companions around 800 pulsars monitored at the Jodrell Bank Observatory, with both circular and eccentric orbits of periods between $20$ days and $17$ years and inclination-dependent planetary masses from $10^{-4}$ to $100\,\mathrm{M}_{\oplus}$. Using a Bayesian framework, we simultaneously model pulsar timing parameters and a stationary noise process with a power-law power spectral density. We put limits on the projected masses of any planetary companions, which reach as low as 1/100th of the mass of the Moon ($\sim 10^{-4}\,\mathrm{M}_{\oplus}$). We find that two-thirds of our pulsars are highly unlikely to host any companions above $2-8\,\mathrm{M}_{\oplus}$. Our results imply that fewer than $0.5\%$ of pulsars could host terrestrial planets as large as those known to orbit PSR B1257$+$12 ($\sim4\,\mathrm{M}_{\oplus}$); however, the smaller planet in this system ($\sim0.02\,\mathrm{M}_{\oplus}$) would be undetectable in $95\%$ of our sample, hidden by both instrumental and intrinsic noise processes, although it is not clear if such tiny planets could exist in isolation. We detect significant periodicities in 15 pulsars, however we find that intrinsic quasi-periodic magnetospheric effects can mimic the influence of a planet, and for the majority of these cases we believe this to be the origin of the detected periodicity. Notably, we find that the highly periodic oscillations in PSR B0144$+$59 are correlated with changes in the pulse profile and therefore can be attributed to magnetospheric effects. We believe the most plausible candidate for planetary companions in our sample is PSR J2007$+$3120. △ Less

Submitted 2 March, 2022; originally announced March 2022.

Comments: Accepted for publication in MNRAS

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.

Comments: 17 pages, 12 figures, accepted in MNRAS

Narrowband searches for continuous and long-duration transient gravitational waves from known pulsars in the LIGO-Virgo third observing run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, T. D. Abbott, F. Acernese, K. Ackley, C. Adams, N. Adhikari, R. X. Adhikari, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, A. Allocca, P. A. Altin, A. Amato , et al. (1636 additional authors not shown)

Abstract: Isolated neutron stars that are asymmetric with respect to their spin axis are possible sources of detectable continuous gravitational waves. This paper presents a fully-coherent search for such signals from eighteen pulsars in data from LIGO and Virgo's third observing run (O3). For known pulsars, efficient and sensitive matched-filter searches can be carried out if one assumes the gravitational… ▽ More Isolated neutron stars that are asymmetric with respect to their spin axis are possible sources of detectable continuous gravitational waves. This paper presents a fully-coherent search for such signals from eighteen pulsars in data from LIGO and Virgo's third observing run (O3). For known pulsars, efficient and sensitive matched-filter searches can be carried out if one assumes the gravitational radiation is phase-locked to the electromagnetic emission. In the search presented here, we relax this assumption and allow the frequency and frequency time-derivative of the gravitational waves to vary in a small range around those inferred from electromagnetic observations. We find no evidence for continuous gravitational waves, and set upper limits on the strain amplitude for each target. These limits are more constraining for seven of the targets than the spin-down limit defined by ascribing all rotational energy loss to gravitational radiation. In an additional search we look in O3 data for long-duration (hours-months) transient gravitational waves in the aftermath of pulsar glitches for six targets with a total of nine glitches. We report two marginal outliers from this search, but find no clear evidence for such emission either. The resulting duration-dependent strain upper limits do not surpass indirect energy constraints for any of these targets. △ Less

Submitted 27 June, 2022; v1 submitted 21 December, 2021; originally announced December 2021.

Comments: 37 pages, 9 figures, submitted to ApJ

Report number: LIGO-P2100267

Journal ref: ApJ, 932, 133 (2022)

A comparative analysis of pulse time-of-arrival creation methods

Authors: J. Wang, G. M. Shaifullah, J. P. W. Verbiest, C. Tiburzi, D. J. Champion, I. Cognard, M. Gaikwad, E. Graikou, L. Guillemot, H. Hu, R. Karuppusamy, Michael J. Keith, Michael Kramer, Y. Liu, A. G. Lyne, M. B. Mickaliger, B. W. Stappers, G. Theureau

Abstract: Extracting precise pulse times of arrival (TOAs) and their uncertainties is the first and most fundamental step in high-precision pulsar timing. In the classical method, TOAs are derived from total intensity pulse profiles of pulsars via cross-correlation with an idealised `1D' template of that profile. While a number of results have been presented in the literature relying on the ever increasing… ▽ More Extracting precise pulse times of arrival (TOAs) and their uncertainties is the first and most fundamental step in high-precision pulsar timing. In the classical method, TOAs are derived from total intensity pulse profiles of pulsars via cross-correlation with an idealised `1D' template of that profile. While a number of results have been presented in the literature relying on the ever increasing sensitivity of such pulsar timing experiments, there is no consensus on the most reliable methods for TOA creation and, more importantly, the associated TOA uncertainties for each scheme. In this article, we present a comprehensive comparison of TOA determination practices, focusing on the creation of timing templates, TOA determination methods and the most useful TOA bandwidth. The aim is both to present a possible approach towards TOA optimisation as well as the (partial) identification of an optimal TOA-creation scheme and the demonstration of optimisation differences between pulsars and data sets. We compare the values of data-derived template profiles as compared to analytic profiles and evaluate the three most commonly used template-matching methods. Finally, we study the relation between timing precision and TOA bandwidth to identify any potential breaks in that relationship. As a practical demonstration, we apply our selected methods to European Pulsar Timing Array data on the three test pulsars PSRs\ J0218+4232, J1713+0747 and J2145$-$0750. △ Less

Submitted 26 November, 2021; originally announced November 2021.

Journal ref: A&A 658, A181 (2022)

Searches for Gravitational Waves from Known Pulsars at Two Harmonics in the Second and Third LIGO-Virgo Observing Runs

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, H. Abe, F. Acernese, K. Ackley, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi, A. Allocca, P. A. Altin , et al. (1672 additional authors not shown)

Abstract: We present a targeted search for continuous gravitational waves (GWs) from 236 pulsars using data from the third observing run of LIGO and Virgo (O3) combined with data from the second observing run (O2). Searches were for emission from the $l=m=2$ mass quadrupole mode with a frequency at only twice the pulsar rotation frequency (single harmonic) and the $l=2, m=1,2$ modes with a frequency of both… ▽ More We present a targeted search for continuous gravitational waves (GWs) from 236 pulsars using data from the third observing run of LIGO and Virgo (O3) combined with data from the second observing run (O2). Searches were for emission from the $l=m=2$ mass quadrupole mode with a frequency at only twice the pulsar rotation frequency (single harmonic) and the $l=2, m=1,2$ modes with a frequency of both once and twice the rotation frequency (dual harmonic). No evidence of GWs was found so we present 95\% credible upper limits on the strain amplitudes $h_0$ for the single harmonic search along with limits on the pulsars' mass quadrupole moments $Q_{22}$ and ellipticities $\varepsilon$. Of the pulsars studied, 23 have strain amplitudes that are lower than the limits calculated from their electromagnetically measured spin-down rates. These pulsars include the millisecond pulsars J0437\textminus4715 and J0711\textminus6830 which have spin-down ratios of 0.87 and 0.57 respectively. For nine pulsars, their spin-down limits have been surpassed for the first time. For the Crab and Vela pulsars our limits are factors of $\sim 100$ and $\sim 20$ more constraining than their spin-down limits, respectively. For the dual harmonic searches, new limits are placed on the strain amplitudes $C_{21}$ and $C_{22}$. For 23 pulsars we also present limits on the emission amplitude assuming dipole radiation as predicted by Brans-Dicke theory. △ Less

Submitted 20 July, 2022; v1 submitted 25 November, 2021; originally announced November 2021.

Comments: 37 pages

Report number: LIGO-P2100049

The Thousand-Pulsar-Array programme on MeerKAT VII: Polarisation properties of pulsars in the Magellanic Clouds

Authors: S. Johnston, A. Parthasarathy, R. A. Main, J. P. Ridley, B. S. Koribalski, M. Bailes, S. J. Buchner, M. Geyer, A. Karastergiou, M. J. Keith, M. Kramer, M. Serylak, R. M. Shannon, R. Spiewak, V. Venkatraman Krishnan

Abstract: The Magellanic Clouds are the only external galaxies known to host radio pulsars. The dispersion and rotation measures of pulsars in the Clouds can aid in understanding their structure, and studies of the pulsars themselves can point to potential differences between them and their Galactic counterparts. We use the high sensitivity of the MeerKAT telescope to observe 17 pulsars in the Small and Lar… ▽ More The Magellanic Clouds are the only external galaxies known to host radio pulsars. The dispersion and rotation measures of pulsars in the Clouds can aid in understanding their structure, and studies of the pulsars themselves can point to potential differences between them and their Galactic counterparts. We use the high sensitivity of the MeerKAT telescope to observe 17 pulsars in the Small and Large Magellanic Clouds in addition to five foreground (Galactic) pulsars. We provide polarization profiles for 18 of these pulsars, improved measurements of their dispersion and rotation measures, and derive the mean parallel magnetic field along the lines of sight. The results are broadly in agreement with expectations for the structure and strength of the magnetic field in the Large and Small Magellanic Clouds. The Magellanic Cloud pulsars have profiles which are narrower than expected from the period-width relationship and we show this is due to selection effects in pulsar surveys rather than any intrinsic difference between the population of Galactic and Magellanic objects. △ Less

Submitted 18 November, 2021; originally announced November 2021.

Comments: Accepted for publication in MNRAS

The Jodrell Bank Glitch Catalogue: 106 new rotational glitches in 70 pulsars

Authors: Avishek Basu, Benjamin Shaw, Danai Antonopoulou, Michael J. Keith, Andrew G. Lyne, Mitchell B. Mickaliger, Benjamin W. Stappers, Patrick Weltevrede, Christine A. Jordan

Abstract: Pulsar glitches are rapid spin-up events that occur in the rotation of neutron stars, providing a valuable probe into the physics of the interiors of these objects. Long-term monitoring of a large number of pulsars facilitates the detection of glitches and the robust measurements of their parameters. The Jodrell Bank pulsar timing programme regularly monitors more than 800 radio pulsars and has ac… ▽ More Pulsar glitches are rapid spin-up events that occur in the rotation of neutron stars, providing a valuable probe into the physics of the interiors of these objects. Long-term monitoring of a large number of pulsars facilitates the detection of glitches and the robust measurements of their parameters. The Jodrell Bank pulsar timing programme regularly monitors more than 800 radio pulsars and has accrued, in some cases, over 50 years of timing history on individual objects. In this paper we present 106 new glitches in 70 radio pulsars as observed up to the end of 2018. For 70% of these pulsars, the event we report is its only known glitch. For each new glitch we provide measurements of its epoch, amplitude and any detected changes to the spin-down rate of the star. Combining these new glitches with those listed in the Jodrell Bank glitch catalogue we analyse a total sample of 543 glitches in 178 pulsars. We model the distribution of glitch amplitudes and spin-down rate changes using a mixture of two Gaussian components. We corroborate the known dependence of glitch rate and activity on pulsar spin-down rates and characteristic ages, and show that younger pulsars tend to exhibit larger glitches. Pulsars whose spin-down rates between $10^{-14}$ Hz s$^{-1}$ and $10^{-10.5}$ Hz s$^{-1}$ show a mean reversal of 1.8% of their spin-down as a consequence of glitches. Our results are qualitatively consistent with the superfluid vortex unpinning models of pulsar glitches. △ Less

Submitted 12 November, 2021; originally announced November 2021.

Comments: 14 pages, 17 figures, 3 tables. Accepted for publication in MNRAS

Noise analysis in the European Pulsar Timing Array data release 2 and its implications on the gravitational-wave background search

Authors: A. Chalumeau, S. Babak, A. Petiteau, S. Chen, A. Samajdar, R. N. Caballero, G. Theureau, L. Guillemot, G. Desvignes, A. Parthasarathy, K. Liu, G. Shaifullah, H. Hu, E. van der Wateren, J. Antoniadis, A. -S. Bak Nielsen, C. G. Bassa, A. Berthereau, M. Burgay, D. J. Champion, I. Cognard, M. Falxa, R. D. Ferdman, P. C. C. Freire, J. R. Gair , et al. (27 additional authors not shown)

Abstract: The European Pulsar Timing Array (EPTA) collaboration has recently released an extended data set for six pulsars (DR2) and reported evidence for a common red noise signal. Here we present a noise analysis for each of the six pulsars. We consider several types of noise: (i) radio frequency independent, "achromatic", and time-correlated red noise; (ii) variations of dispersion measure and scattering… ▽ More The European Pulsar Timing Array (EPTA) collaboration has recently released an extended data set for six pulsars (DR2) and reported evidence for a common red noise signal. Here we present a noise analysis for each of the six pulsars. We consider several types of noise: (i) radio frequency independent, "achromatic", and time-correlated red noise; (ii) variations of dispersion measure and scattering; (iii) system and band noise; and (iv) deterministic signals (other than gravitational waves) that could be present in the PTA data. We perform Bayesian model selection to find the optimal combination of noise components for each pulsar. Using these custom models we revisit the presence of the common uncorrelated red noise signal previously reported in the EPTA DR2 and show that the data still supports it with a high statistical significance. Next, we confirm that there is no preference for or against the Hellings-Downs spatial correlations expected for the stochastic gravitational-wave background. The main conclusion of the EPTA DR2 paper remains unchanged despite a very significant change in the noise model of each pulsar. However, modelling the noise is essential for the robust detection of gravitational waves and its impact could be significant when analysing the next EPTA data release, which will include a larger number of pulsars and more precise measurements. △ Less

Submitted 9 November, 2021; originally announced November 2021.

Comments: 21 pages, 11 figures, 7 tables, 1 appendix figure and 1 appendix table, accepted for publication to MNRAS

Common-red-signal analysis with 24-yr high-precision timing of the European Pulsar Timing Array: Inferences in the stochastic gravitational-wave background search

Authors: S. Chen, R. N. Caballero, Y. J. Guo, A. Chalumeau, K. Liu, G. Shaifullah, K. J. Lee, S. Babak, G. Desvignes, A. Parthasarathy, H. Hu, E. van der Wateren, J. Antoniadis, A. -S. Bak Nielsen, C. G. Bassa, A. Berthereau, M. Burgay, D. J. Champion, I. Cognard, M. Falxa, R. D. Ferdman, P. C. C. Freire, J. R. Gair, E. Graikou, L. Guillemot , et al. (27 additional authors not shown)

Abstract: We present results from the search for a stochastic gravitational-wave background (GWB) as predicted by the theory of General Relativity using six radio millisecond pulsars from the Data Release 2 (DR2) of the European Pulsar Timing Array (EPTA) covering a timespan up to 24 years. A GWB manifests itself as a long-term low-frequency stochastic signal common to all pulsars, a common red signal (CRS)… ▽ More We present results from the search for a stochastic gravitational-wave background (GWB) as predicted by the theory of General Relativity using six radio millisecond pulsars from the Data Release 2 (DR2) of the European Pulsar Timing Array (EPTA) covering a timespan up to 24 years. A GWB manifests itself as a long-term low-frequency stochastic signal common to all pulsars, a common red signal (CRS), with the characteristic Hellings-Downs (HD) spatial correlation. Our analysis is performed with two independent pipelines, \eprise{} and \tn{}+\ftwo{}, which produce consistent results. A search for a CRS with simultaneous estimation of its spatial correlations yields spectral properties compatible with theoretical GWB predictions, but does not result in the required measurement of the HD correlation, as required for GWB detection. Further Bayesian model comparison between different types of CRSs, including a GWB, finds the most favoured model to be the common uncorrelated red noise described by a power-law with $A = 5.13_{-2.73}^{+4.20} \times 10^{-15}$ and $γ= 3.78_{-0.59}^{+0.69}$ (95\% credible regions). Fixing the spectral index to $γ=13/3$ as expected from the GWB by circular, inspiralling supermassive black-hole binaries results in an amplitude of $A =2.95_{-0.72}^{+0.89} \times 10^{-15}$. We implement three different models, BAYESEPHEM, LINIMOSS and EPHEMGP, to address possible Solar-system ephemeris (SSE) systematics and conclude that our results may only marginally depend on these effects. This work builds on the methods and models from the studies on the EPTA DR1. We show that under the same analysis framework the results remain consistent after the data set extension. △ Less

Submitted 25 October, 2021; originally announced October 2021.

Comments: 25 pages, 10 figures, 7 tables, 2 appendix tables and 1 appendix figure

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.

Comments: Accepted for publication in MNRAS. 20 pages, 19 figures, 7 tables

A broadband radio study of PSR J0250+5854: the slowest-spinning radio pulsar known

Authors: C. H. Agar, P. Weltevrede, L. Bondonneau, J. -M. Grießmeier, J. W. T. Hessels, W. J. Huang, A. Karastergiou, M. J. Keith, V. I. Kondratiev, J. Künsemöller, D. Li, B. Peng, C. Sobey, B. W. Stappers, C. M. Tan, G. Theureau, H. G. Wang, C. M. Zhang, B. Cecconi, J. N. Girard, A. Loh, P. Zarka

Abstract: We present radio observations of the most slowly rotating known radio pulsar PSR J0250+5854. With a 23.5 s period, it is close, or even beyond, the $P$-$\dot{P}$ diagram region thought to be occupied by active pulsars. The simultaneous observations with FAST, the Chilbolton and Effelsberg LOFAR international stations, and NenuFAR represent a five-fold increase in the spectral coverage of this obje… ▽ More We present radio observations of the most slowly rotating known radio pulsar PSR J0250+5854. With a 23.5 s period, it is close, or even beyond, the $P$-$\dot{P}$ diagram region thought to be occupied by active pulsars. The simultaneous observations with FAST, the Chilbolton and Effelsberg LOFAR international stations, and NenuFAR represent a five-fold increase in the spectral coverage of this object, with the detections at 1250 MHz (FAST) and 57 MHz (NenuFAR) being the highest- and lowest-frequency published respectively to date. We measure a flux density of $4\pm2$ $μ$Jy at 1250 MHz and an exceptionally steep spectral index of $-3.5^{+0.2}_{-1.5}$, with a turnover below $\sim$95 MHz. In conjunction with observations of this pulsar with the GBT and the LOFAR Core, we show that the intrinsic profile width increases drastically towards higher frequencies, contrary to the predictions of conventional radius-to-frequency mapping. We examine polarimetric data from FAST and the LOFAR Core and conclude that its polar cap radio emission is produced at an absolute height of several hundreds of kilometres around 1.5 GHz, similar to other rotation-powered pulsars across the population. Its beam is significantly underfilled at lower frequencies, or it narrows because of the disappearance of conal outriders. Finally, the results for PSR J0250+5854 and other slowly spinning rotation-powered pulsars are contrasted with the radio-detected magnetars. We conclude that magnetars have intrinsically wider radio beams than the slow rotation-powered pulsars, and that consequently the latter's lower beaming fraction is what makes objects such as PSR J0250+5854 so scarce. △ Less

Submitted 1 September, 2021; originally announced September 2021.

Comments: 14 pages, 6 figures

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.

Comments: Accepted for publication in MNRAS. 14 pages, 14 figures, 2 tables

The slow rise and recovery of the 2019 Crab pulsar glitch

Authors: B. Shaw, M. J. Keith, A. G. Lyne, M. B. Mickaliger, B. W. Stappers, J. D. Turner, P. Weltevrede

Abstract: We present updated measurements of the Crab pulsar glitch of 2019 July 23 using a dataset of pulse arrival times spanning $\sim$5 months. On MJD 58687, the pulsar underwent its seventh largest glitch observed to date, characterised by an instantaneous spin-up of $\sim$1 $μ$Hz. Following the glitch the pulsar's rotation frequency relaxed exponentially towards pre-glitch values over a timescale of a… ▽ More We present updated measurements of the Crab pulsar glitch of 2019 July 23 using a dataset of pulse arrival times spanning $\sim$5 months. On MJD 58687, the pulsar underwent its seventh largest glitch observed to date, characterised by an instantaneous spin-up of $\sim$1 $μ$Hz. Following the glitch the pulsar's rotation frequency relaxed exponentially towards pre-glitch values over a timescale of approximately one week, resulting in a permanent frequency increment of $\sim$0.5 $μ$Hz. Due to our semi-continuous monitoring of the Crab pulsar, we were able to partially resolve a fraction of the total spin-up. This delayed spin-up occurred exponentially over a timescale of $\sim$18 hours. This is the sixth Crab pulsar glitch for which part of the initial rise was resolved in time and this phenomenon has not been observed in any other glitching pulsars, offering a unique opportunity to study the microphysical processes governing interactions between the neutron star interior and the crust. △ Less

Submitted 26 March, 2021; v1 submitted 24 March, 2021; originally announced March 2021.

Comments: 5 pages, 3 figures

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.

Comments: 16 pages, 9 figures

Journal ref: https://ui.adsabs.harvard.edu/abs/2021MNRAS.tmp...69P/abstract

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.

Comments: 13 pages, 8 figures. Accepted for publication in MNRAS

A precise mass measurement of PSR J2045+3633

Authors: J. W. McKee, P. C. C. Freire, M. Berezina, D. J. Champion, I. Cognard, E. Graikou, L. Guillemot, M. J. Keith, M. Kramer, A. G. Lyne, B. W. Stappers, T. M. Tauris, G. Theureau

Abstract: We present the results of a timing analysis undertaken with the goal of obtaining an improved mass measurement of the recycled pulsar J2045+3633. Using regular high-cadence observations with the Effelsberg, Nançay, and Lovell radio telescopes, together with targeted campaigns with the Arecibo Telescope and Effelsberg, we have assembled a 6-yr timing data set for this pulsar. We measure highly sign… ▽ More We present the results of a timing analysis undertaken with the goal of obtaining an improved mass measurement of the recycled pulsar J2045+3633. Using regular high-cadence observations with the Effelsberg, Nançay, and Lovell radio telescopes, together with targeted campaigns with the Arecibo Telescope and Effelsberg, we have assembled a 6-yr timing data set for this pulsar. We measure highly significant values for the proper motion and the related rate of change of orbital semi-major axis ($\dot{x}$), and have obtained high precision values of the rate of advance of periastron time ($\dotω$), and two of the Shapiro delay parameters ($h_{3}$ and $ς$). This has allowed us to improve the measurements of the pulsar and companion masses by an order of magnitude, yielding (with $1σ$ uncertainties) $1.251^{+0.021}_{-0.021}\,\text{M}_{\odot}$ for PSR J2045+3633, and $0.873^{+0.016}_{-0.014}\,\text{M}_{\odot}$ for its white dwarf companion, and has allowed us to place improved constraints on the geometrical orientation of the binary system. Using our measurements of the binary component masses and the orbital size, we consider possible evolutionary scenarios for the system. △ Less

Submitted 27 September, 2020; v1 submitted 25 September, 2020; originally announced September 2020.

Comments: 17 pages, 6 figures, accepted for publication in MNRAS

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.

Comments: 11 pages, 4 figures, accepted for publication by the MNRAS

High-cadence observations and variable spin behaviour of magnetar Swift J1818.0-1607 after its outburst

Authors: David Champion, Ismael Cognard, Marilyn Cruces, Gregory Desvignes, Fabian Jankowski, Ramesh Karuppusamy, Michael J. Keith, Chryssa Kouveliotou, Michael Kramer, Kuo Liu, Andrew G. Lyne, Mitchell B. Mickaliger, Brendan O'Connor, Aditya Parthasarathy, Nataliya Porayko, Kaustubh Rajwade, Ben W. Stappers, Pablo Torne, Alexander J. van der Horst, Patrick Weltevrede

Abstract: We report on multi-frequency radio observations of the new magnetar Swift J1818.0-1607, following it for more than one month with high cadence. The observations commenced less than 35 hours after its registered first outburst. We obtained timing, polarisation and spectral information. Swift J1818.0-1607 has an unusually steep spectrum for a radio emitting magnetar and also has a relatively narrow… ▽ More We report on multi-frequency radio observations of the new magnetar Swift J1818.0-1607, following it for more than one month with high cadence. The observations commenced less than 35 hours after its registered first outburst. We obtained timing, polarisation and spectral information. Swift J1818.0-1607 has an unusually steep spectrum for a radio emitting magnetar and also has a relatively narrow and simple pulse profile. The position angle swing of the polarisation is flat over the pulse profile, possibly suggesting that our line-of-sight grazes the edge of the emission beam. This may also explain the steep spectrum. The spin evolution shows large variation in the spin-down rate, associated with four distinct timing events over the course of our observations. Those events may be related to the appearance and disappearance of a second pulse component. The first timing event coincides with our actual observations, while we did not detect significant changes in the emission properties which could reveal further magnetospheric changes. Characteristic ages inferred from the timing measurements over the course of months vary by nearly an order of magnitude. A longer-term spin-down measurement over approximately 100 days suggests an characteristic age of about 500 years, larger than previously reported. Though Swift J1818.0-1607 could still be one of the youngest neutron stars (and magnetars) detected so far, we caution using the characteristic age as a true-age indicator given the caveats behind its calculation. △ Less

Submitted 8 September, 2020; originally announced September 2020.

Comments: 14 pages, 10 figures, accepted by MNRAS

Gravitational-wave constraints on the equatorial ellipticity of millisecond pulsars

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, R. Abbott, T. D. Abbott, S. Abraham, F. Acernese, K. Ackley, A. Adams, C. Adams, R. X. Adhikari, V. B. Adya, C. Affeldt, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, G. Allen, A. Allocca, P. A. Altin, A. Amato, S. Anand, A. Ananyeva , et al. (1311 additional authors not shown)

Abstract: We present a search for continuous gravitational waves from five radio pulsars, comprising three recycled pulsars (PSR J0437-4715, PSR J0711-6830, and PSR J0737-3039A) and two young pulsars: the Crab pulsar (J0534+2200) and the Vela pulsar (J0835-4510). We use data from the third observing run of Advanced LIGO and Virgo combined with data from their first and second observing runs. For the first t… ▽ More We present a search for continuous gravitational waves from five radio pulsars, comprising three recycled pulsars (PSR J0437-4715, PSR J0711-6830, and PSR J0737-3039A) and two young pulsars: the Crab pulsar (J0534+2200) and the Vela pulsar (J0835-4510). We use data from the third observing run of Advanced LIGO and Virgo combined with data from their first and second observing runs. For the first time we are able to match (for PSR J0437-4715) or surpass (for PSR J0711-6830) the indirect limits on gravitational-wave emission from recycled pulsars inferred from their observed spin-downs, and constrain their equatorial ellipticities to be less than $10^{-8}$. For each of the five pulsars, we perform targeted searches that assume a tight coupling between the gravitational-wave and electromagnetic signal phase evolution. We also present constraints on PSR J0711-6830, the Crab pulsar and the Vela pulsar from a search that relaxes this assumption, allowing the gravitational-wave signal to vary from the electromagnetic expectation within a narrow band of frequencies and frequency derivatives. △ Less

Submitted 13 October, 2020; v1 submitted 28 July, 2020; originally announced July 2020.

Comments: 23 pages, 6 figures, 3 tables

Report number: LIGO-P2000029

Journal ref: 2020 ApJL 902 L21

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.

Comments: 20 pages, 16 Figures, 4 Tables, accepted for publication in PASA

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.

Comments: 26 pages, 6 figures. Accepted for publication in PASA. The data release can be downloaded from https://doi.org/10.25919/5db90a8bdeb59