Showing 1–50 of 184 results for author: Champion, D

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.

The TRAPUM Small Magellanic Cloud pulsar survey with MeerKAT: I. Discovery of seven new pulsars and two Pulsar Wind Nebula associations

Authors: E. Carli, L. Levin, B. W. Stappers, E. D. Barr, R. P. Breton, S. Buchner, M. Burgay, M. Geyer, M. Kramer, P. V. Padmanabh, A. Possenti, V. Venkatraman Krishnan, W. Becker, M. D. Filipović, C. Maitra, J. Behrend, D. J. Champion, W. Chen, Y. P. Men, A. Ridolfi

Abstract: The sensitivity of the MeerKAT radio interferometer is an opportunity to probe deeper into the population of rare and faint extragalactic pulsars. The TRAPUM (TRAnsients and PUlsars with MeerKAT) collaboration has conducted a radio-domain search for accelerated pulsars and transients in the Small Magellanic Cloud (SMC). This partially targeted survey, performed at L-band (856-1712 MHz) with the co… ▽ More The sensitivity of the MeerKAT radio interferometer is an opportunity to probe deeper into the population of rare and faint extragalactic pulsars. The TRAPUM (TRAnsients and PUlsars with MeerKAT) collaboration has conducted a radio-domain search for accelerated pulsars and transients in the Small Magellanic Cloud (SMC). This partially targeted survey, performed at L-band (856-1712 MHz) with the core array of the MeerKAT telescope in 2-h integrations, is twice as sensitive as the latest SMC radio pulsar survey. We report the discovery of seven new SMC pulsars, doubling this galaxy's radio pulsar population and increasing the total extragalactic population by nearly a quarter. We also carried out a search for accelerated millisecond pulsars in the SMC Globular Cluster NGC 121 using the full array of MeerKAT. This improved the previous upper limit on pulsed radio emission from this cluster by a factor of six. Our discoveries reveal the first radio pulsar-PWN systems in the SMC, with only one such system previously known outside our galaxy (the "Crab pulsar twin" in the Large Magellanic Cloud, PSR J0540$-$6919). We associate the 59 ms pulsar discovery PSR J0040$-$7337, now the fastest spinning radio pulsar in the SMC, with the bow-shock Pulsar Wind Nebula (PWN) of Supernova Remnant DEM S5. We also present a new young pulsar with a 79 ms period, PSR J0048$-$7317, in a PWN recently discovered in a MeerKAT radio continuum image. Using the multi-beam capability of MeerKAT, we localised our pulsar discoveries, and two previous Murriyang discoveries, to a positional uncertainty of a few arcseconds. △ Less

Submitted 20 May, 2024; originally announced May 2024.

Comments: 32 pages, 14 figures, 10 tables. Accepted for publication in Monthly Notices of the Royal Astronomical Society

TRAPUM search for pulsars in supernova remnants and pulsar wind nebulae -- I. Survey description and initial discoveries

Authors: J. D. Turner, B. W. Stappers, E. Carli, E. D. Barr, W. Becker, J. Behrend, R. P. Breton, S. Buchner, M. Burgay, D. J. Champion, W. Chen, C. J. Clark, D. M. Horn, E. F. Keane, M. Kramer, L. K ünkel, L. Levin, Y. P. Men, P. V. Padmanabh, A. Ridolfi, V. Venkatraman Krishnan

Abstract: We present the description and initial results of the TRAPUM (TRAnsients And PUlsars with MeerKAT) search for pulsars associated with supernova remnants (SNRs), pulsar wind nebulae and unidentified TeV emission. The list of sources to be targeted includes a large number of well-known candidate pulsar locations but also new candidate SNRs identified using a range of criteria. Using the 64-dish Meer… ▽ More We present the description and initial results of the TRAPUM (TRAnsients And PUlsars with MeerKAT) search for pulsars associated with supernova remnants (SNRs), pulsar wind nebulae and unidentified TeV emission. The list of sources to be targeted includes a large number of well-known candidate pulsar locations but also new candidate SNRs identified using a range of criteria. Using the 64-dish MeerKAT radio telescope, we use an interferometric beamforming technique to tile the potential pulsar locations with coherent beams which we search for radio pulsations, above a signal-to-noise of 9, down to an average flux density upper limit of 30 $μ$Jy. This limit is target-dependent due to the contribution of the sky and nebula to the system temperature. Coherent beams are arranged to overlap at their 50 per cent power radius, so the sensitivity to pulsars is not degraded by more than this amount, though realistically averages around 65 per cent if every location in the beam is considered. We report the discovery of two new pulsars; PSR J1831$-$0941 is an adolescent pulsar likely to be the plerionic engine of the candidate PWN G20.0+0.0, and PSR J1818$-$1502 appears to be an old and faint pulsar that we serendipitously discovered near the centre of a SNR already hosting a compact central object. The survey holds importance for better understanding of neutron star birth rates and the energetics of young pulsars. △ Less

Submitted 20 May, 2024; originally announced May 2024.

Comments: 17 pages, 6 figures, 3 tables. Accepted for publication in Monthly Notices of the Royal Astronomical Society

Modeling non stationary noise in pulsar timing array data analysis

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

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

Submitted 6 May, 2024; originally announced May 2024.

Comments: 14 pages, 14 figures

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

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

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

Submitted 2 May, 2024; originally announced May 2024.

Comments: 17 pages, 5 figures

Pipeline Provenance for Analysis, Evaluation, Trust or Reproducibility

Authors: Michael A. C. Johnson, Hans-Rainer Klöckner, Albina Muzafarova, Kristen Lackeos, David J. Champion, Marta Dembska, Sirko Schindler, Marcus Paradies

Abstract: Data volumes and rates of research infrastructures will continue to increase in the upcoming years and impact how we interact with their final data products. Little of the processed data can be directly investigated and most of it will be automatically processed with as little user interaction as possible. Capturing all necessary information of such processing ensures reproducibility of the final… ▽ More Data volumes and rates of research infrastructures will continue to increase in the upcoming years and impact how we interact with their final data products. Little of the processed data can be directly investigated and most of it will be automatically processed with as little user interaction as possible. Capturing all necessary information of such processing ensures reproducibility of the final results and generates trust in the entire process. We present PRAETOR, a software suite that enables automated generation, modelling, and analysis of provenance information of Python pipelines. Furthermore, the evaluation of the pipeline performance, based upon a user defined quality matrix in the provenance, enables the first step of machine learning processes, where such information can be fed into dedicated optimisation procedures. △ Less

Submitted 22 April, 2024; originally announced April 2024.

Comments: 4 pages, 3 figures

Journal ref: Research Notes of the AAS 8.4 (2024): 100

The Relativistic Spin Precession in the Compact Double Neutron Star System PSR~J1946+2052

Authors: Lingqi Meng, Weiwei Zhu, Michael Kramer, Xueli Miao, Gregory Desvignes, Lijing Shao, Huanchen Hu, Paulo C. C. Freire, Yongkun Zhang, Mengyao Xue, Ziyao Fang, David J. Champion, Mao Yuan, Chenchen Miao, Jiarui Niu, Qiuyang Fu, Jumei Yao, Yanjun Guo, Chengmin Zhang

Abstract: We observe systematic profile changes in the visible pulsar of the compact double neutron star system PSR~J1946+2052 using observations with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). The interpulse of PSR~J1946+2052 changed from single-peak to double-peak shape from 2018 to 2021. We attribute this evolution as the result of the relativistic spin precession of the pulsar. Wi… ▽ More We observe systematic profile changes in the visible pulsar of the compact double neutron star system PSR~J1946+2052 using observations with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). The interpulse of PSR~J1946+2052 changed from single-peak to double-peak shape from 2018 to 2021. We attribute this evolution as the result of the relativistic spin precession of the pulsar. With the high sensitivity of FAST, we also measure significant polarization for the first time, allowing us to model this with the precessional rotating vector model. Assuming, to the first order, a circular hollow-cone-like emission beam pattern and taking the validity of general relativity, we derive the binary's orbital inclination angle (${63^\circ}^{+5^\circ}_{-3^\circ}$) and pulsar's spin geometry. Pulsar's spin vector and the orbital angular momentum vector are found to be only slightly misaligned (${0.21^\circ}^{+0.28^\circ}_{-0.10^\circ}$).The quoted uncertainties do not reflect the systematic uncertainties introduced by our model assumptions. By simulating future observations of profile and polarization evolution, we estimate that we could constrain the precession rate within a $43\%$ uncertainty in 9 years. Hence, we suggest that the system's profile evolution could be combined with precise pulsar timing to test general relativity in the future. △ Less

Submitted 26 March, 2024; originally announced March 2024.

Comments: 12 pages, 9 figures, accepted for publication in ApJ

Discovery and timing of ten new millisecond pulsars in the globular cluster Terzan 5

Authors: P. V. Padmanabh, S. M. Ransom, P. C. C. Freire, A. Ridolfi, J. D. Taylor, C. Choza, C. J. Clark, F. Abbate, M. Bailes, E. D. Barr, S. Buchner, M. Burgay, M. E. DeCesar, W. Chen, A. Corongiu, D. J. Champion, A. Dutta, M. Geyer, J. W. T. Hessels, M. Kramer, A. Possenti, I. H. Stairs, B. W. Stappers, V. Venkatraman Krishnan, L. Vleeschower , et al. (1 additional authors not shown)

Abstract: We report the discovery of ten new pulsars in the globular cluster Terzan 5 as part of the Transients and Pulsars with MeerKAT (TRAPUM) Large Survey Project. We observed Terzan 5 at L-band (856--1712 MHz) with the MeerKAT radio telescope for four hours on two epochs, and performed acceleration searches of 45 out of 288 tied-array beams covering the core of the cluster. We obtained phase-connected… ▽ More We report the discovery of ten new pulsars in the globular cluster Terzan 5 as part of the Transients and Pulsars with MeerKAT (TRAPUM) Large Survey Project. We observed Terzan 5 at L-band (856--1712 MHz) with the MeerKAT radio telescope for four hours on two epochs, and performed acceleration searches of 45 out of 288 tied-array beams covering the core of the cluster. We obtained phase-connected timing solutions for nine discoveries, covering nearly two decades of archival observations from the Green Bank Telescope for all but one. Highlights include PSR J1748$-$2446ao which is an eccentric ($e = 0.32$) wide-orbit (orbital period $P_{\rm b} = 57.55$ d) system. We were able to measure the rate of advance of periastron ($\dotω$) for this system allowing us to determine a total mass of $3.17 \pm \, 0.02\, \rm M_{\odot}$. With a minimum companion mass ($M_{\rm c}$) of $\sim 0.8\, \rm M_{\odot}$, PSR J1748$-$2446ao is a candidate double neutron star (DNS) system. If confirmed to be a DNS, it would be the fastest spinning pulsar ($P = 2.27$ ms) and the longest orbital period measured for any known DNS system. PSR J1748$-$2446ap has the second highest eccentricity for any recycled pulsar ($e \sim 0.905$) and for this system we can measure the total mass ($1.997 \pm 0.006\, \rm M_{\odot}$) and also estimate the individual pulsar and companion masses. PSR J1748$-$2446ar is an eclipsing redback (minimum $M_{\rm c} \sim 0.34\, \rm M_{\odot}$) system whose properties confirm it to be the counterpart to a previously published source identified in radio and X-ray imaging. With these discoveries, the total number of confirmed pulsars in Terzan 5 is 49, the highest for any globular cluster so far. These discoveries further enhance the rich set of pulsars known in Terzan 5 and provide scope for a deeper understanding of binary stellar evolution, cluster dynamics and ensemble population studies. △ Less

Submitted 19 June, 2024; v1 submitted 26 March, 2024; originally announced March 2024.

Comments: 23 pages, 11 figures, 5 tables, published in A&A

Journal ref: A&A 686, A166 (2024)

A targeted radio pulsar survey of redback candidates with MeerKAT

Authors: T. Thongmeearkom, C. J. Clark, R. P. Breton, M. Burgay, L. Nieder, P. C. C. Freire, E. D. Barr, B. W. Stappers, S. M. Ransom, S. Buchner, F. Calore, D. J. Champion, I. Cognard, J. -M. Grießmeier, M. Kramer, L. Levin, P. V. Padmanabh, A. Possenti, A. Ridolfi, V. Venkatraman Krishnan, L. Vleeschower

Abstract: Redbacks are millisecond pulsar binaries with low mass, irradiated companions. These systems have a rich phenomenology that can be used to probe binary evolution models, pulsar wind physics, and the neutron star mass distribution. A number of high-confidence redback candidates have been identified through searches for variable optical and X-ray sources within the localisation regions of unidentifi… ▽ More Redbacks are millisecond pulsar binaries with low mass, irradiated companions. These systems have a rich phenomenology that can be used to probe binary evolution models, pulsar wind physics, and the neutron star mass distribution. A number of high-confidence redback candidates have been identified through searches for variable optical and X-ray sources within the localisation regions of unidentified but pulsar-like Fermi-LAT gamma-ray sources. However, these candidates remain unconfirmed until pulsations are detected. As part of the TRAPUM project, we searched for radio pulsations from six of these redback candidates with MeerKAT. We discovered three new radio millisecond pulsars, PSRs J0838$-$2527, J0955$-$3947 and J2333$-$5526, confirming their redback nature. PSR J0838$-$2827 remained undetected for two years after our discovery despite repeated observations, likely due to evaporated material absorbing the radio emission for long periods of time. While, to our knowledge, this system has not undergone a transition to an accreting state, the disappearance, likely caused by extreme eclipses, illustrates the transient nature of spider pulsars and the heavy selection bias in uncovering their radio population. Radio timing enabled the detection of gamma-ray pulsations from all three pulsars, from which we obtained 15-year timing solutions. All of these sources exhibit complex orbital period variations consistent with gravitational quadrupole moment variations in the companion stars. These timing solutions also constrain the binary mass ratios, allowing us to narrow down the pulsar masses. We find that PSR J2333$-$5526 may have a neutron star mass in excess of 2 M$_{\odot}$. △ Less

Submitted 14 March, 2024; originally announced March 2024.

Comments: 19 pages, 7 figures, accepted for publication in MNRAS

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

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

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

Submitted 18 January, 2024; originally announced January 2024.

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

A 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

A MeerKAT view of the double pulsar eclipses -- Geodetic precession of pulsar B and system geometry

Authors: M. E. Lower, M. Kramer, R. M. Shannon, R. P. Breton, N. Wex, S. Johnston, M. Bailes, S. Buchner, H. Hu, V. Venkatraman Krishnan, V. A. Blackmon, F. Camilo, D. J. Champion, P. C. C. Freire, M. Geyer, A. Karastergiou, J. van Leeuwen, M. A. McLaughlin, D. J. Reardon, I. H. Stairs

Abstract: The double pulsar system, PSR J0737$-$3039A/B, consists of two neutron stars bound together in a highly relativistic orbit that is viewed nearly edge-on from the Earth. This alignment results in brief radio eclipses of the fast-rotating pulsar A when it passes behind the toroidal magnetosphere of the slow-rotating pulsar B. The morphology of these eclipses is strongly dependent on the geometric or… ▽ More The double pulsar system, PSR J0737$-$3039A/B, consists of two neutron stars bound together in a highly relativistic orbit that is viewed nearly edge-on from the Earth. This alignment results in brief radio eclipses of the fast-rotating pulsar A when it passes behind the toroidal magnetosphere of the slow-rotating pulsar B. The morphology of these eclipses is strongly dependent on the geometric orientation and rotation phase of pulsar B, and their time-evolution can be used to constrain the geodetic precession rate of the pulsar. We demonstrate a Bayesian inference framework for modelling eclipse light-curves obtained with MeerKAT between 2019-2023. Using a hierarchical inference approach, we obtained a precession rate of $Ω_{\rm SO}^{\rm B} = {5.16^{\circ}}^{+0.32^{\circ}}_{-0.34^{\circ}}$ yr$^{-1}$ for pulsar B, consistent with predictions from General Relativity to a relative uncertainty of 6.5%. This updated measurement provides a 6.1% test of relativistic spin-orbit coupling in the strong-field regime. We show that a simultaneous fit to all of our observed eclipses can in principle return a $\sim$1.5% test of spin-orbit coupling. However, systematic effects introduced by the current geometric orientation of pulsar B along with inconsistencies between the observed and predicted eclipse light curves result in difficult to quantify uncertainties. Assuming the validity of General Relativity, we definitively show that the spin-axis of pulsar B is misaligned from the total angular momentum vector by $40.6^{\circ} \pm 0.1^{\circ}$ and that the orbit of the system is inclined by approximately $90.5^{\circ}$ from the direction of our line of sight. Our measured geometry for pulsar B suggests the largely empty emission cone contains an elongated horseshoe shaped beam centered on the magnetic axis, and that it may not be re-detected as a radio pulsar until early-2035. △ Less

Submitted 1 February, 2024; v1 submitted 10 November, 2023; originally announced November 2023.

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

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

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

A MeerKAT view of the pulsars in the globular cluster NGC 6522

Authors: F. Abbate, A. Ridolfi, P. C. C. Freire, P. V. Padmanabh, V. Balakrishnan, S. Buchner, L. Zhang, M. Kramer, B. W. Stappers, E. D. Barr, W. Chen, D. Champion, S. Ransom, A. Possenti

Abstract: We present the results of observations aimed at discovering and studying pulsars in the core-collapsed globular cluster (GC) NGC 6522 performed by the MeerTIME and TRAPUM Large Survey Project with the MeerKAT telescope. We have discovered two new isolated pulsars bringing the total number of known pulsars in the cluster to six. PSR J1803$-$3002E is a mildly recycled pulsar with spin period of 17.9… ▽ More We present the results of observations aimed at discovering and studying pulsars in the core-collapsed globular cluster (GC) NGC 6522 performed by the MeerTIME and TRAPUM Large Survey Project with the MeerKAT telescope. We have discovered two new isolated pulsars bringing the total number of known pulsars in the cluster to six. PSR J1803$-$3002E is a mildly recycled pulsar with spin period of 17.9 ms while pulsar PSR J1803$-$3002F is a slow pulsar with spin period of 148.1 ms. The presence of isolated and slow pulsars is expected in NGC 6522 and confirms the predictions of previous theories for clusters at this stage in evolution. We further present a tentative timing solution for the millisecond pulsar (MSP) PSR J1803$-$3002C combining older observations taken with the Parkes 64m radio telescope, Murriyang. This solution implies a relatively small characteristic age of the pulsar in contrast with the old age of the GC. The presence of a slow pulsar and an apparently young MSP, both rare in GCs, suggests that their formation might be linked to the evolutionary stage of the cluster. △ Less

Submitted 5 October, 2023; originally announced October 2023.

Comments: 8 pages, 5 figues, 3 tables. Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 680, A47 (2023)

Comparing recent PTA results on the nanohertz stochastic gravitational wave background

Authors: The International Pulsar Timing Array Collaboration, G. Agazie, J. Antoniadis, A. Anumarlapudi, A. M. Archibald, P. Arumugam, S. Arumugam, Z. Arzoumanian, J. Askew, S. Babak, M. Bagchi, M. Bailes, A. -S. Bak Nielsen, P. T. Baker, C. G. Bassa, A. Bathula, B. Bécsy, A. Berthereau, N. D. R. Bhat, L. Blecha, M. Bonetti, E. Bortolas, A. Brazier, P. R. Brook, M. Burgay , et al. (220 additional authors not shown)

Abstract: The Australian, Chinese, European, Indian, and North American pulsar timing array (PTA) collaborations recently reported, at varying levels, evidence for the presence of a nanohertz gravitational wave background (GWB). Given that each PTA made different choices in modeling their data, we perform a comparison of the GWB and individual pulsar noise parameters across the results reported from the PTA… ▽ More The Australian, Chinese, European, Indian, and North American pulsar timing array (PTA) collaborations recently reported, at varying levels, evidence for the presence of a nanohertz gravitational wave background (GWB). Given that each PTA made different choices in modeling their data, we perform a comparison of the GWB and individual pulsar noise parameters across the results reported from the PTAs that constitute the International Pulsar Timing Array (IPTA). We show that despite making different modeling choices, there is no significant difference in the GWB parameters that are measured by the different PTAs, agreeing within $1σ$. The pulsar noise parameters are also consistent between different PTAs for the majority of the pulsars included in these analyses. We bridge the differences in modeling choices by adopting a standardized noise model for all pulsars and PTAs, finding that under this model there is a reduction in the tension in the pulsar noise parameters. As part of this reanalysis, we "extended" each PTA's data set by adding extra pulsars that were not timed by that PTA. Under these extensions, we find better constraints on the GWB amplitude and a higher signal-to-noise ratio for the Hellings and Downs correlations. These extensions serve as a prelude to the benefits offered by a full combination of data across all pulsars in the IPTA, i.e., the IPTA's Data Release 3, which will involve not just adding in additional pulsars, but also including data from all three PTAs where any given pulsar is timed by more than as single PTA. △ Less

Submitted 1 September, 2023; originally announced September 2023.

Comments: 21 pages, 9 figures, submitted to ApJ

The MPIfR-MeerKAT Galactic Plane Survey II. The eccentric double neutron star system PSR J1208-5936 and a neutron star merger rate update

Authors: M. Colom i Bernadich, V. Balakrishnan, E. Barr, M. Berezina, M. Burgay, S. Buchner, D. J. Champion, W. Chen, G. Desvignes, P. C. C. Freire, K. Grunthal, M. Kramer, Y. Men, P. V. Padmanabh, A. Parthasarathy, D. Pillay, I. Rammala, S. Sengupta, V. Venkatraman Krishnan

Abstract: The MMGPS-L is the most sensitive pulsar survey in the Southern Hemisphere. We present a follow-up study of one of these new discoveries, PSR J1208-5936, a 28.71-ms recycled pulsar in a double neutron star system with an orbital period of Pb=0.632 days and an eccentricity of e=0.348. Through timing of almost one year of observations, we detected the relativistic advance of periastron (0.918(1) deg… ▽ More The MMGPS-L is the most sensitive pulsar survey in the Southern Hemisphere. We present a follow-up study of one of these new discoveries, PSR J1208-5936, a 28.71-ms recycled pulsar in a double neutron star system with an orbital period of Pb=0.632 days and an eccentricity of e=0.348. Through timing of almost one year of observations, we detected the relativistic advance of periastron (0.918(1) deg/yr), resulting in a total system mass of Mt=2.586(5) Mo. We also achieved low-significance constraints on the amplitude of the Einstein delay and Shapiro delay, in turn yielding constraints on the pulsar mass (Mp=1.26(+0.13/-0.25) Mo), the companion mass (Mc=1.32(+0.25/-0.13) Mo, and the inclination angle (i=57(12) degrees). This system is highly eccentric compared to other Galactic field double neutron stars with similar periods, possibly hinting at a larger-than-usual supernova kick during the formation of the second-born neutron star. The binary will merge within 7.2(2) Gyr due to the emission of gravitational waves. With the improved sensitivity of the MMGPS-L, we updated the Milky Way neutron star merger rate to be 25(+19/-9) Myr$^{-1}$ within 90% credible intervals, which is lower than previous studies based on known Galactic binaries owing to the lack of further detections despite the highly sensitive nature of the survey. This implies a local cosmic neutron star merger rate of 293(+222/-103} Gpc/yr, consistent with LIGO and Virgo O3 observations. With this, we predict the observation of 10(+8/-4) neutron star merger events during the LIGO-Virgo-KAGRA O4 run. We predict the uncertainties on the component masses and the inclination angle will be reduced to 5x10$^{-3}$ Mo and 0.4 degrees after two decades of timing, and that in at least a decade from now the detection of the shift in Pb and the sky proper motion will serve to make an independent constraint of the distance to the system. △ Less

Submitted 8 September, 2023; v1 submitted 31 August, 2023; originally announced August 2023.

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

Variability, polarimetry, and timing properties of single pulses from PSR J2222-0137 using FAST

Authors: X. L. Miao, W. W. Zhu, M. Kramer, P. C. C. Freire, L. Shao, M. Yuan, L. Q. Meng, Z. W. Wu, C. C. Miao, Y. J. Guo, D. J. Champion, E. Fonseca, J. M. Yao, M. Y. Xue, J. R. Niu, H. Hu, C. M. Zhang

Abstract: In our work, we analyse $5\times10^{4}$ single pulses from the recycled pulsar PSR J2222$-$0137 in one of its scintillation maxima observed by the Five-hundred-meter Aperture Spherical radio Telescope (FAST). PSR J2222$-$0137 is one of the nearest and best studies of binary pulsars and a unique laboratory for testing gravitational theories. We report single pulses' energy distribution and polariza… ▽ More In our work, we analyse $5\times10^{4}$ single pulses from the recycled pulsar PSR J2222$-$0137 in one of its scintillation maxima observed by the Five-hundred-meter Aperture Spherical radio Telescope (FAST). PSR J2222$-$0137 is one of the nearest and best studies of binary pulsars and a unique laboratory for testing gravitational theories. We report single pulses' energy distribution and polarization from the pulsar's main-pulse region. The single pulse energy follows the log-normal distribution. We resolve a steep polarization swing, but at the current time resolution ($64\,μ{\rm s}$), we find no evidence for the orthogonal jump in the main-pulse region, as has been suspected. We find a potential sub-pulse drifting period of $P_{3} \sim 3.5\,P$. We analyse the jitter noise from different integrated numbers of pulses and find that its $σ_{j}$ is $270\pm{9}\,{\rm ns}$ for 1-hr integration at 1.25 GHz. This result is useful for optimizing future timing campaigns with FAST or other radio telescopes. △ Less

Submitted 21 August, 2023; originally announced August 2023.

Comments: 11 pages, 14 figures, accepted by Monthly Notices of the Royal Astronomical Society

Journal ref: MNRAS 526 (2023) 2156

Periodic interstellar scintillation variations of PSRs~J0613$-$0200 and J0636+5128 associated with the Local Bubble shell

Authors: Yulan Liu, Robert A. Main, Joris P. W. Verbiest, Ziwei Wu, Krishnakumar M. Ambalappat, Jiguang Lu, David J. Champion, Ismaël Cognard, Lucas Guillemot, Kuo Liu, James W. McKee, Nataliya Porayko, Golam. M. Shaifullah, Gilles Theureau

Abstract: Annual variations of interstellar scintillation can be modelled to constrain parameters of the ionized interstellar medium. If a pulsar is in a binary system, then investigating the orbital parameters is possible through analysis of the orbital variation of scintillation. In observations carried out from 2011 January to 2020 August by the European Pulsar Timing Array radio telescopes, PSRs~J0613… ▽ More Annual variations of interstellar scintillation can be modelled to constrain parameters of the ionized interstellar medium. If a pulsar is in a binary system, then investigating the orbital parameters is possible through analysis of the orbital variation of scintillation. In observations carried out from 2011 January to 2020 August by the European Pulsar Timing Array radio telescopes, PSRs~J0613$-$0200 and J0636+5128 show strong annual variations in their scintillation velocity, while the former additionally exhibits an orbital fluctuation. Bayesian theory and Markov-chain-Monte-Carlo methods are used to interpret these periodic variations. We assume a thin and anisotropic scattering screen model, and discuss the mildly and extremely anisotropic scattering cases. PSR~J0613$-$0200 is best described by mildly anisotropic scattering, while PSR~J0636+5128 exhibits extremely anisotropic scattering. We measure the distance, velocity and degree of anisotropy of the scattering screen for our two pulsars, finding that scattering screen distances from Earth for PSRs~J0613$-$0200 and J0636+5128 are 316$^{+28}_{-20}$\,pc and 262$^{+96}_{-38}$\,pc, respectively. The positions of these scattering screens are coincident with the shell of the Local Bubble towards both pulsars. These associations add to the growing evidence of the Local Bubble shell as a dominant region of scattering along many sightlines. △ Less

Submitted 18 August, 2023; v1 submitted 19 July, 2023; originally announced July 2023.

Comments: Accepted by SCIENCE CHINA Physics, Mechanics & Astronomy ( SCPMA)

Reciprocating Magnetic Fields in the Pulsar Wind Observed from the Black Widow Pulsar J1720-0534

Authors: Chen-Chen Miao, Victoria Blackmon, Wei-Wei Zhu, Dong-Zi Li, Mingyu Ge, Xiao-Peng You, Maura McLaughlin, Di Li, Na Wang, Pei Wang, Jia-Rui Niu, M. Cruces, Jian-Ping Yuan, Jun-Tao Bai, D. J. Champion, Yu-Tong Chen, Ming-Min Chi, P. C. C. Freire, Yi Feng, Zhen-Ye Gan, M. Kramer, Fei-Fei Kou, Yu-Xi Li, Xue-Li Miao, Ling-Qi Meng , et al. (19 additional authors not shown)

Abstract: We report the radio observations of the eclipsing black widow pulsar J1720-0534, a 3.26 ms pulsar in orbit with a low mass companion of mass 0.029 to 0.034 M$_{\odot}$. We obtain the phase-connected timing ephemeris and polarization profile of this millisecond pulsar (MSP) using the Five-hundred-meter Aperture Spherical Radio Telescope (FAST), the Green Bank Telescope (GBT), and the Parkes Telesco… ▽ More We report the radio observations of the eclipsing black widow pulsar J1720-0534, a 3.26 ms pulsar in orbit with a low mass companion of mass 0.029 to 0.034 M$_{\odot}$. We obtain the phase-connected timing ephemeris and polarization profile of this millisecond pulsar (MSP) using the Five-hundred-meter Aperture Spherical Radio Telescope (FAST), the Green Bank Telescope (GBT), and the Parkes Telescope. For the first time from such a system, an oscillatory polarisation angle change was observed from a particular eclipse egress with partial depolarization, indicating 10-milliGauss-level reciprocating magnetic fields oscillating in a length scale of 5000 km (assuming an orbital inclination angle of 90 degrees) outside the companion's magnetosphere. The dispersion measure variation observed during the ingresses and egresses shows the rapid raising of the electron density in the shock boundary between the companion's magnetosphere and the surrounding pulsar wind. We suggest that the observed oscillatory magnetic fields originate from the pulsar wind outside the companion's magnetosphere. △ Less

Submitted 28 August, 2023; v1 submitted 2 July, 2023; originally announced July 2023.

Comments: 15 pages, 8 figures, 1 table, accepted by RAA

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

High-altitude Magnetospheric Emissions from Two Pulsars

Authors: Mao Yuan, Weiwei Zhu, Michael Kramer, Bo Peng, Jiguang Lu, Renxin Xu, Lijing Shao, Hong-guang Wang, Lingqi Meng, Jiarui Niu, Rushuang Zhao, Chenchen Miao, Xueli Miao, Mengyao Xue, Yi Feng, Pei Wang, Di Li, Chengmin Zhang, David J. Champion, Emmanuel Fonseca, Huanchen Hu, Jumei Yao, Paulo C. C. Freire, Yanjun Guo

Abstract: We discover three new weak pulse components in two known pulsars, one in PSR J0304+1932 and two in PSR J1518+4904. These components are emitted about half way between the main emission beam and the interpulse beam (beam from the opposite pole). They are separated from their main pulse peak by $99^{\circ}\pm{3}^{\circ}$ for J0304+1932, $123^{\circ}.6\pm{0^{\circ}.7}$ (leading) and… ▽ More We discover three new weak pulse components in two known pulsars, one in PSR J0304+1932 and two in PSR J1518+4904. These components are emitted about half way between the main emission beam and the interpulse beam (beam from the opposite pole). They are separated from their main pulse peak by $99^{\circ}\pm{3}^{\circ}$ for J0304+1932, $123^{\circ}.6\pm{0^{\circ}.7}$ (leading) and $93^{^{\circ}}\pm 0^{\circ}.4$ (trailing) for J1518+4904, respectively. Their peak-intensity ratios to main pulses are: $\sim$ 0.06% for J0304+1932, $\sim$ 0.17% and $\sim$ 0.83% for J1518+4904. We also analyzed flux fluctuation and profile variation of the emissions for two pulsars. The results show correlations between the weak pulses and their main pulses, indicating that these emissions come from the same pole. We estimated the emission altitude of these weak pulses and derived a height of about half of the pulsar's light-cylinder radius. These pulse components are a unique sample of high-altitude emissions from pulsars, and challenge the current pulsar emission models. △ Less

Submitted 4 July, 2023; v1 submitted 8 June, 2023; originally announced June 2023.

Comments: 10 pages, 13 figures, Accepted for publication in ApJ

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

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

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

Submitted 24 May, 2023; originally announced May 2023.

Comments: 23 pages, 16 figures, 7 tables

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 continuous Gravitational Waves in the second data release of the International Pulsar Timing Array

Authors: M. Falxa, S. Babak, P. T. Baker, B. Bécsy, A. Chalumeau, S. Chen, Z. Chen, N. J. Cornish, L. Guillemot, J. S. Hazboun, C. M. F. Mingarelli, A. Parthasarathy, A. Petiteau, N. S. Pol, A. Sesana, S. B. Spolaor, S. R. Taylor, G. Theureau, M. Vallisneri, S. J. Vigeland, C. A. Witt, X. Zhu, J. Antoniadis, Z. Arzoumanian, M. Bailes , et al. (102 additional authors not shown)

Abstract: The International Pulsar Timing Array 2nd data release is the combination of datasets from worldwide collaborations. In this study, we search for continuous waves: gravitational wave signals produced by individual supermassive black hole binaries in the local universe. We consider binaries on circular orbits and neglect the evolution of orbital frequency over the observational span. We find no evi… ▽ More The International Pulsar Timing Array 2nd data release is the combination of datasets from worldwide collaborations. In this study, we search for continuous waves: gravitational wave signals produced by individual supermassive black hole binaries in the local universe. We consider binaries on circular orbits and neglect the evolution of orbital frequency over the observational span. We find no evidence for such signals and set sky averaged 95% upper limits on their amplitude h 95 . The most sensitive frequency is 10nHz with h 95 = 9.1 10-15 . We achieved the best upper limit to date at low and high frequencies of the PTA band thanks to improved effective cadence of observations. In our analysis, we have taken into account the recently discovered common red noise process, which has an impact at low frequencies. We also find that the peculiar noise features present in some pulsars data must be taken into account to reduce the false alarm. We show that using custom noise models is essential in searching for continuous gravitational wave signals and setting the upper limit. △ Less

Submitted 19 March, 2023; originally announced March 2023.

The MPIfR-MeerKAT Galactic Plane survey I -- System setup and early results

Authors: P. V. Padmanabh, E. D. Barr, S. S. Sridhar, M. R. Rugel, A. Damas-Segovia, A. M. Jacob, V. Balakrishnan, M. Berezina, M. C. i Bernadich, A. Brunthaler, D. J. Champion, P. C. C. Freire, S. Khan, H. -R. Klöckner, M. Kramer, Y. K. Ma, S. A. Mao, Y. P. Men, K. M. Menten, S. Sengupta, V. Venkatraman Krishnan, O. Wucknitz, F. Wyrowski, M. C. Bezuidenhout, S. Buchner , et al. (8 additional authors not shown)

Abstract: Galactic plane radio surveys play a key role in improving our understanding of a wide range of astrophysical phenomena. Performing such a survey using the latest interferometric telescopes produces large data rates necessitating a shift towards fully or quasi-real-time data analysis with data being stored for only the time required to process them. We present here the overview and setup for the 30… ▽ More Galactic plane radio surveys play a key role in improving our understanding of a wide range of astrophysical phenomena. Performing such a survey using the latest interferometric telescopes produces large data rates necessitating a shift towards fully or quasi-real-time data analysis with data being stored for only the time required to process them. We present here the overview and setup for the 3000 hour Max-Planck-Institut fuer Radioastronomie (MPIfR) MeerKAT Galactic Plane survey (MMGPS). The survey is unique by operating in a commensal mode, addressing key science objectives of the survey including the discovery of new pulsars and transients as well as studies of Galactic magnetism, the interstellar medium and star formation rates. We explain the strategy coupled with the necessary hardware and software infrastructure needed for data reduction in the imaging, spectral and time domains. We have so far discovered 78 new pulsars including 17 confirmed binary systems of which two are potential double neutron star systems. We have also developed an imaging pipeline sensitive to the order of a few tens of micro-Jansky with a spatial resolution of a few arcseconds. Further science operations with an in-house built S-Band receiver operating between 1.7-3.5 GHz are about to commence. Early spectral line commissioning observations conducted at S-Band, targeting transitions of the key molecular gas tracer CH at 3.3 GHz already illustrate the spectroscopic capabilities of this instrument. These results lay a strong foundation for future surveys with telescopes like the Square Kilometre Array (SKA). △ Less

Submitted 21 June, 2023; v1 submitted 16 March, 2023; originally announced March 2023.

Comments: 25 pages, 10 figures, Accepted in MNRAS

Missing for 20 years: MeerKAT re-detects the elusive binary pulsar M30B

Authors: Vishnu Balakrishnan, Paulo Freire, Scott Ransom, Alessandro Ridolfi, Ewan Barr, Weiwei Chen, Vivek Venkatraman Krishnan, David J. Champion, Michael Kramer, Tasha Gautam, Prajwal Padmanabh, Yunpeng Men, Federico Abbate, Benjamin Stappers, Ingrid Stairs, Evan Keane, Andrea Possenti

Abstract: PSR J2140$-$2311B is a 13-ms pulsar discovered in 2001 in a 7.8-hour Green Bank Telescope (GBT) observation of the core-collapsed globular cluster M30 and predicted to be in a highly eccentric binary orbit. This pulsar has eluded detection since then, therefore its precise orbital parameters have remained a mystery until now. In this work, we present the confirmation of this pulsar using observati… ▽ More PSR J2140$-$2311B is a 13-ms pulsar discovered in 2001 in a 7.8-hour Green Bank Telescope (GBT) observation of the core-collapsed globular cluster M30 and predicted to be in a highly eccentric binary orbit. This pulsar has eluded detection since then, therefore its precise orbital parameters have remained a mystery until now. In this work, we present the confirmation of this pulsar using observations taken with the UHF receivers of the MeerKAT telescope as part of the TRAPUM Large Survey Project. Taking advantage of the beamforming capability of our backends, we have localized it, placing it $1.2(1)^\prime$ from the cluster centre. Our observations have enabled the determination of its orbit: it is highly eccentric ($e = 0.879$) with an orbital period of $6.2$ days. We also measured the rate of periastron advance, $\dotω = 0.078 \pm 0.002\, \rm deg \, yr^{-1}$. Assuming that this effect is fully relativistic, general relativity provides an estimate of the total mass of the system, $M_{\rm TOT} = 2.53 \pm 0.08$ M$_{\odot}$, consistent with the lightest double neutron star systems known. Combining this with the mass function of the system gives the pulsar and companion masses of $m_p < 1.43 \, \rm M_{\odot}$ and $m_c > 1.10 \, \rm M_{\odot}$ respectively. The massive, undetected companion could either be a massive WD or a NS. M30B likely formed as a result of a secondary exchange encounter. Future timing observations will allow the determination of a phase-coherent timing solution, vastly improving our uncertainty in $\dotω$ and likely enabling the detection of additional relativistic effects which will determine $m_p$ and $m_c$. △ Less

Submitted 12 January, 2023; originally announced January 2023.

Comments: Accepted for publication in the Astrophysical Journal Letters (ApJL)

PSR~J1910$-$5959A: A rare gravitational laboratory for testing white dwarf models

Authors: A. Corongiu, V. Venkatraman Krishnan, P. C. C. Freire, M. Kramer, A. Possenti, M. Geyer, A. Ridolfi, F. Abbate, M. Bailes, E. D. Barr, V. Balakrishnan, S. Buchner, D. J. Champion, W. Chen, B. V. Hugo, A. Karastergiou, A. G. Lyne, R. N. Manchester, P. V. Padmanabh, A. Parthasarathy, S. M. Ransom, J. M. Sarkissian, M. Serylak, W. van Straten

Abstract: PSRJ1910-5959A (J1910A) is a binary millisecond pulsar in a 0.837 day circular orbit around a helium white dwarf (HeWD) companion. This pulsar is located 6.3 arcmin away from the centre of the globular cluster NGC6752. Given the large offset, the association of the pulsar to NGC6752 has been debated. We have made use of two decades of archival Parkes 64-m "Murriyang" telescope data and recently ca… ▽ More PSRJ1910-5959A (J1910A) is a binary millisecond pulsar in a 0.837 day circular orbit around a helium white dwarf (HeWD) companion. This pulsar is located 6.3 arcmin away from the centre of the globular cluster NGC6752. Given the large offset, the association of the pulsar to NGC6752 has been debated. We have made use of two decades of archival Parkes 64-m "Murriyang" telescope data and recently carried out observations with the MeerKAT telescope. We obtained Pulse times of arrival using standard data reduction techniques and analysed using Bayesian pulsar timing techniques. We analysed the pulsar's total intensity and polarisation profile, to study the interstellar scattering along the line of sight, and the pulsar's geometry by applying the rotating vector model. We obtain precise measurements of several post-Keplerian parameters: the range $r=0.202(6)T_\odot$ and shape s=0.999823(4) of the Shapiro delay, from which we infer the orbital inclination to be $88.9^{+0.15}_{-0.14}°$ and the masses of both the pulsar and the companion to be $1.55(7)M_{\odot}$ and $0.202(6)M_{\odot}$ respectively; a secular change in the orbital period $\dot{P}_{\rm b}=-53^{+7.4}_{-6.0}\times 10^{-15}$\,s\,s$^{-1}$ that proves the association to NGGC6752 and a secular change in the projected semi-major axis of the pulsar $\dot{x}= -40.7^{+7.3}_{-8.2}\times10^{-16}$\,s\,s$^{-1}$ that is likely caused by the spin-orbit interaction from a misaligned HeWD spin, at odds with the likely isolated binary evolution of the system. We also discuss some theoretical models for the structure and evolution of WDs in NS-WD binaries by using J1910A's companion as a test bed. J1910A is a rare system for which several parameters of both the pulsar and the HeWD companion can be accurately measured. As such, it is a test bed to discriminate between alternative models for HeWD structure and cooling. △ Less

Submitted 10 February, 2023; v1 submitted 10 January, 2023; originally announced January 2023.

Comments: 18 pages, 12 figures. Abstract shortened and rephrased because of arXiv characters limits. Accepted for publication on A&A on January 9th 2023 Updated version after linguistic revision

Journal ref: A&A 671, A72 (2023)

MeerKAT discovery of 13 new pulsars in Omega Centauri

Authors: W. Chen, P. C. C. Freire, A. Ridolfi, E. D. Barr, B. Stappers, M. Kramer, A. Possenti, S. M. Ransom, L. Levin, R. P. Breton, M. Burgay, F. Camilo, S. Buchner, D. J. Champion, F. Abbate, V. Venkatraman Krishnan, P. V. Padmanabh, T. Gautam, L. Vleeschower, M. Geyer, J-M. Grießmeier, Y. P. Men, V. Balakrishnan, M. C. Bezuidenhout

Abstract: The most massive globular cluster in our Galaxy, Omega Centauri, is an interesting target for pulsar searches, because of its multiple stellar populations and the intriguing possibility that it was once the nucleus of a galaxy that was absorbed into the Milky Way. The recent discoveries of pulsars in this globular cluster and their association with known X-ray sources was a hint that, given the la… ▽ More The most massive globular cluster in our Galaxy, Omega Centauri, is an interesting target for pulsar searches, because of its multiple stellar populations and the intriguing possibility that it was once the nucleus of a galaxy that was absorbed into the Milky Way. The recent discoveries of pulsars in this globular cluster and their association with known X-ray sources was a hint that, given the large number of known X-ray sources, there is a much larger undiscovered pulsar population. We used the superior sensitivity of the MeerKAT radio telescope to search for pulsars in Omega Centauri. In this paper, we present some of the first results of this survey, including the discovery of 13 new pulsars; the total number of known pulsars in this cluster currently stands at 18. At least half of them are in binary systems and preliminary orbital constraints suggest that most of the binaries have light companions. We also discuss the ratio between isolated and binaries pulsars and how they were formed in this cluster. △ Less

Submitted 10 January, 2023; originally announced January 2023.

The TRAPUM L-band survey for pulsars in Fermi-LAT gamma-ray sources

Authors: C. J. Clark, R. P. Breton, E. D. Barr, M. Burgay, T. Thongmeearkom, L. Nieder, S. Buchner, B. Stappers, M. Kramer, W. Becker, M. Mayer, A. Phosrisom, A. Ashok, M. C. Bezuidenhout, F. Calore, I. Cognard, P. C. C. Freire, M. Geyer, J. -M. Grießmeier, R. Karuppusamy, L. Levin, P. V. Padmanabh, A. Possenti, S. Ransom, M. Serylak , et al. (13 additional authors not shown)

Abstract: More than 100 millisecond pulsars (MSPs) have been discovered in radio observations of gamma-ray sources detected by the Fermi Large Area Telescope (LAT), but hundreds of pulsar-like sources remain unidentified. Here we present the first results from the targeted survey of Fermi-LAT sources being performed by the Transients and Pulsars with MeerKAT (TRAPUM) Large Survey Project. We observed 79 sou… ▽ More More than 100 millisecond pulsars (MSPs) have been discovered in radio observations of gamma-ray sources detected by the Fermi Large Area Telescope (LAT), but hundreds of pulsar-like sources remain unidentified. Here we present the first results from the targeted survey of Fermi-LAT sources being performed by the Transients and Pulsars with MeerKAT (TRAPUM) Large Survey Project. We observed 79 sources identified as possible gamma-ray pulsar candidates by a Random Forest classification of unassociated sources from the 4FGL catalogue. Each source was observed for 10 minutes on two separate epochs using MeerKAT's L-band receiver (856-1712 MHz), with typical pulsed flux density sensitivities of $\sim$100$\,μ$Jy. Nine new MSPs were discovered, eight of which are in binary systems, including two eclipsing redbacks and one system, PSR J1526$-$2744, that appears to have a white dwarf companion in an unusually compact 5 hr orbit. We obtained phase-connected timing solutions for two of these MSPs, enabling the detection of gamma-ray pulsations in the Fermi-LAT data. A follow-up search for continuous gravitational waves from PSR J1526$-$2744 in Advanced LIGO data using the resulting Fermi-LAT timing ephemeris yielded no detection, but sets an upper limit on the neutron star ellipticity of $2.45\times10^{-8}$. We also detected X-ray emission from the redback PSR J1803$-$6707 in data from the first eROSITA all-sky survey, likely due to emission from an intra-binary shock. △ Less

Submitted 16 December, 2022; originally announced December 2022.

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

Journal ref: MNRAS, 519, 5590-5606 (2023)

Searches for Shapiro delay in seven binary pulsars using the MeerKAT telescope

Authors: Mohsen Shamohammadi, Matthew Bailes, Paulo C. C. Freire, Aditya Parthasarathy, Daniel J. Reardon, Ryan M. Shannon, Vivek Venkatraman Krishnan, Miquel C. i. Bernadich, Andrew D. Cameron, David J. Champion, Alessandro Corongiu, Christopher Flynn, Marisa Geyer, Michael Kramer, Matthew T. Miles, Andrea Possenti, Renee Spiewak

Abstract: Precision timing of millisecond pulsars in binary systems enables observers to detect the relativistic Shapiro delay induced by space time curvature. When favourably aligned, this enables constraints to be placed on the component masses and system orientation. Here we present the results of timing campaigns on seven binary millisecond pulsars observed with the 64-antenna MeerKAT radio telescope th… ▽ More Precision timing of millisecond pulsars in binary systems enables observers to detect the relativistic Shapiro delay induced by space time curvature. When favourably aligned, this enables constraints to be placed on the component masses and system orientation. Here we present the results of timing campaigns on seven binary millisecond pulsars observed with the 64-antenna MeerKAT radio telescope that show evidence of Shapiro delay: PSRs~J0101$-$6422, J1101$-$6424, J1125$-$6014, J1514$-$4946, J1614$-$2230, J1732$-$5049, and J1909$-$3744. Evidence for Shapiro delay was found in all of the systems, and for three the orientations and data quality enabled strong constraints on their orbital inclinations and component masses. For PSRs~J1125$-$6014, J1614$-$2230 and J1909$-$3744, we determined pulsar masses to be $M_{\rm p} = 1.68\pm 0.17 \, {\rm M_{\odot}} $, $1.94\pm 0.03 \, {\rm M_{\odot}} $ and $1.45 \pm 0.03 \, {\rm M_{\odot}}$, and companion masses to be $M_{\rm c} = 0.33\pm 0.02 \, {\rm M_{\odot}} $, $0.495\pm 0.005 \, {\rm M_{\odot}} $ and $0.205 \pm 0.003 \, {\rm M_{\odot}}$, respectively. This provides the first independent confirmation of PSR~J1614$-$2230's mass, one of the highest known. The Shapiro delays measured for PSRs~J0101$-$6422, J1101$-$6424, J1514$-$4946, and J1732$-$5049 were only weak, and could not provide interesting component mass limits. Despite a large number of millisecond pulsars being routinely timed, relatively few have accurate masses via Shapiro delays. We use simulations to show that this is expected, and provide a formula for observers to assess how accurately a pulsar mass can be determined. We also discuss the observed correlation between pulsar companion masses and spin period, and the anti-correlation between recycled pulsar mass and their companion masses. △ Less

Submitted 7 December, 2022; originally announced December 2022.

Comments: 19 pages, 8 figures

TRAPUM upper limits on pulsed radio emission for SMC X-ray pulsar J0058-7218

Authors: E. Carli, L. Levin, B. W. Stappers, E. D. Barr, R. P. Breton, S. Buchner, M. Burgay, M. Kramer, P. V. Padmanabh, A. Possenti, V. Venkatraman Krishnan, J. Behrend, D. J. Champion, W. Chen, Y. P. Men

Abstract: The TRAPUM collaboration has used the MeerKAT telescope to conduct a search for pulsed radio emission from the young Small Magellanic Cloud pulsar J0058-7218 located in the supernova remnant IKT 16, following its discovery in X-rays with XMM-Newton. We report no significant detection of dispersed, pulsed radio emission from this source in three 2-hour L-band observations using the core dishes of M… ▽ More The TRAPUM collaboration has used the MeerKAT telescope to conduct a search for pulsed radio emission from the young Small Magellanic Cloud pulsar J0058-7218 located in the supernova remnant IKT 16, following its discovery in X-rays with XMM-Newton. We report no significant detection of dispersed, pulsed radio emission from this source in three 2-hour L-band observations using the core dishes of MeerKAT, setting an upper limit of 7.0 μJy on its mean flux density at 1284 MHz. This is nearly 7 times deeper than previous radio searches for this pulsar in Parkes L-band observations. This suggests that the radio emission of PSR J0058-7218 is not beamed towards Earth or that PSR J0058-7218 is similar to a handful of Pulsar Wind Nebulae systems that have a very low radio efficiency, such as PSR B0540-6919, the Large Magellanic Cloud Crab pulsar analogue. We have also searched for bright, dispersed, single radio pulses and found no candidates above a fluence of 93 mJy ms at 1284 MHz. △ Less

Submitted 10 October, 2022; originally announced October 2022.

Comments: 7 pages, 2 figures, 1 table. Accepted for publication in Monthly Notices of the Royal Astronomical Society

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

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

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

Submitted 23 September, 2022; originally announced September 2022.

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

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

The SPAN512 mid-latitude pulsar survey at the Nançay Radio Telescope

Authors: Gregory Desvignes, Ismael Cognard, David A. Smith, David Champion, Lucas Guillemot, Michael Kramer, Patrice Lespagnol, Franck Octau, Gilles Theureau

Abstract: The large number of ongoing surveys for pulsars and transients at various radio observatories is motivated by the science obtained from these sources. Timing and polarisation analysis of relativistic binaries can place strong constraints on theories of gravity. The observation of a growing number of millisecond pulsars (MSPs) spread over the celestial sphere may allow the detection of a stochastic… ▽ More The large number of ongoing surveys for pulsars and transients at various radio observatories is motivated by the science obtained from these sources. Timing and polarisation analysis of relativistic binaries can place strong constraints on theories of gravity. The observation of a growing number of millisecond pulsars (MSPs) spread over the celestial sphere may allow the detection of a stochastic gravitational wave background arising from supermassive black hole binaries. A more complete sample of young pulsars improves our knowledge of neutron star birth and evolution. Transients such as fast radio bursts can serve to probe the intergalactic medium. The SPAN512 pulsar survey covers intermediate Galactic latitudes using the L-band receiver of the Nançay Radio Telescope (NRT). The survey covers 224 sq. deg. of the sky for a total exposure time of 2200 h. Population syntheses predict the discovery of 3 to 19 new normal pulsars and a few MSPs. We present detailed modelling of the NRT beam with its L-band receiver and its sensitivity which we used to precisely assess the expected survey yield. We used the flexible Pulsar Arecibo L-band Feed Array data processing pipeline to search the 47 TB of SPAN512 data for pulsars and transients. The SPAN512 survey discovered two new MSPs and one new middle-aged pulsar. We focus on the analysis of the 2.4-ms spin period pulsar J2205+6012 for which we also report the detection of gamma-ray pulsations. Its narrow pulse width (35 $μ$s at an observing frequency of 2.55 GHz) allows for sub-microsecond timing precision over 8 years, with exciting prospects for pulsar timing array programs. △ Less

Submitted 5 September, 2022; originally announced September 2022.

Comments: 13 pages, 11 figures, 5 tables. Accepted in A&A

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

Radio detection of an elusive millisecond pulsar in the Globular Cluster NGC 6397

Authors: Lei Zhang, Alessandro Ridolfi, Harsha Blumer, Paulo Freire, Richard N. Manchester, Maura McLaughlin, Kyle Kremer, Andrew D. Cameron, Zhiyu Zhang, Jan Behrend, Marta Burgay, Sarah Buchner, David J. Champion, Weiwei Chen, Shi Dai, Yi Feng, Xiaoting Fu, Meng Guo, George Hobbs, Evan F. Keane, Michael Kramer, Lina Levin, Xiangdong Li, Mengmeng Ni, Jingshan Pan , et al. (10 additional authors not shown)

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

Submitted 16 July, 2022; originally announced July 2022.

Arecibo and FAST Timing Follow-up of twelve Millisecond Pulsars Discovered in Commensal Radio Astronomy FAST Survey

Authors: C. C. Miao, W. W. Zhu, D. Li, P. C. C. Freire, J. R. Niu, P. Wang, J. P. Yuan, M. Y. Xue, A. D. Cameron, D. J. Champion, M. Cruces, Y. T. Chen, M. M. Chi, X. F. Cheng, S. J. Dang, M. F. Ding, Y. Feng, Z. Y. Gan, G. Hobbs, M. Kramer, Z. J. Liu, Y. X. Li, Z. K. Luo, X. L. Miao, L. Q. Meng , et al. (24 additional authors not shown)

Abstract: We report the phase-connected timing ephemeris, polarization pulse profiles, Faraday rotation measurements, and Rotating-Vector-Model (RVM) fitting results of twelve millisecond pulsars (MSPs) discovered with the Five-hundred-meter Aperture Spherical radio Telescope (FAST) in the Commensal radio Astronomy FAST survey (CRAFTS). The timing campaigns were carried out with FAST and Arecibo over three… ▽ More We report the phase-connected timing ephemeris, polarization pulse profiles, Faraday rotation measurements, and Rotating-Vector-Model (RVM) fitting results of twelve millisecond pulsars (MSPs) discovered with the Five-hundred-meter Aperture Spherical radio Telescope (FAST) in the Commensal radio Astronomy FAST survey (CRAFTS). The timing campaigns were carried out with FAST and Arecibo over three years. Eleven of the twelve pulsars are in neutron star - white dwarf binary systems, with orbital periods between 2.4 and 100 d. Ten of them have spin periods, companion masses, and orbital eccentricities that are consistent with the theoretical expectations for MSP - Helium white dwarf (He WD) systems. The last binary pulsar (PSR J1912$-$0952) has a significantly smaller spin frequency and a smaller companion mass, the latter could be caused by a low orbital inclination for the system. Its orbital period of 29 days is well within the range of orbital periods where some MSP - He WD systems have shown anomalous eccentricities, however, the eccentricity of PSR J1912$-$0952 is typical of what one finds for the remaining MSP - He WD systems. △ Less

Submitted 9 May, 2022; originally announced May 2022.

Comments: 11 pages, 5 figures, MNRAS accepted

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

Four pulsar discoveries in NGC 6624 by TRAPUM using MeerKAT

Authors: F. Abbate, A. Ridolfi, E. D. Barr, S. Buchner, M. Burgay, D. J. Champion, W. Chen, P. C. C. Freire, T. Gautam, J. M. Grießmeier, L. Künkel, M. Kramer, P. V. Padmanabh, A. Possenti, S. Ransom, M. Serylak, B. W. Stappers, V. Venkatraman Krishnan, J. Behrend, R. P. Breton, L. Levin, Y. Men

Abstract: We report 4 new pulsars discovered in the core-collapsed globular cluster (GC) NGC 6624 by the TRAPUM Large Survey Project with the MeerKAT telescope. All of the new pulsars found are isolated. PSR J1823$-$3021I and PSR J1823$-$3021K are millisecond pulsars with period of respectively 4.319 ms and 2.768 ms. PSR J1823$-$3021J is mildly recycled with a period of 20.899 ms, and PSR J1823$-$3022 is a… ▽ More We report 4 new pulsars discovered in the core-collapsed globular cluster (GC) NGC 6624 by the TRAPUM Large Survey Project with the MeerKAT telescope. All of the new pulsars found are isolated. PSR J1823$-$3021I and PSR J1823$-$3021K are millisecond pulsars with period of respectively 4.319 ms and 2.768 ms. PSR J1823$-$3021J is mildly recycled with a period of 20.899 ms, and PSR J1823$-$3022 is a long period pulsar with a period of 2.497 s. The pulsars J1823$-$3021I, J1823$-$3021J, and J1823$-$3021K have position and dispersion measure (DM) compatible with being members of the GC and are therefore associated with NGC 6624. Pulsar J1823$-$3022 is the only pulsar bright enough to be re-detected in archival observations of the cluster. This allowed the determination of a timing solution that spans over two decades. It is not possible at the moment to claim the association of pulsar J1823$-$3022 with the GC given the long period and large offset in position ($\sim 3$ arcminutes) and DM (with a fractional difference of 11 percent compared the average of the pulsars in NGC 6624). The discoveries made use of the beamforming capability of the TRAPUM backend to generate multiple beams in the same field of view which allows sensitive searches to be performed over a few half-light radii from the cluster center and can simultaneously localise the discoveries. The discoveries reflect the properties expected for pulsars in core-collapsed GCs. △ Less

Submitted 11 April, 2022; originally announced April 2022.

Comments: Accepted for publication on Monthly Notices of the Royal Astronomical Society. 11 pages, 6 figures

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)

News and views regarding PSR J1757-1854, a highly-relativistic binary pulsar

Authors: A. D. Cameron, M. Bailes, V. Balakrishnan, D. J. Champion, P. C. C. Freire, M. Kramer, N. Wex, S. Johnston, A. G. Lyne, B. W. Stappers, M. A. McLaughlin, N. Pol, H. Wahl, C. Ng, A. Possenti, A. Ridolfi

Abstract: We provide an update on the ongoing monitoring and study of the highly-relativistic double neutron star binary system PSR J1757-1854, a 21.5-ms pulsar in a highly eccentric, 4.4-hour orbit. The extreme nature of this pulsar's orbit allows it to probe a parameter space largely unexplored by other relativistic binary pulsars. For example, it displays one of the highest gravitational wave (GW) lumino… ▽ More We provide an update on the ongoing monitoring and study of the highly-relativistic double neutron star binary system PSR J1757-1854, a 21.5-ms pulsar in a highly eccentric, 4.4-hour orbit. The extreme nature of this pulsar's orbit allows it to probe a parameter space largely unexplored by other relativistic binary pulsars. For example, it displays one of the highest gravitational wave (GW) luminosities of any known binary pulsar, as well as the highest rate or orbital decay due to GW damping. PSR J1757-1854 is also notable in that it is an excellent candidate for exploring new tests of General Relativity and other gravitational theories, with possible measurements of both Lense-Thirring precession and relativistic orbital deformation (through the post-Keplerian parameter $δ_θ$) anticipated within the next 3-5 years. Here we present a summary of the latest interim results from the ongoing monitoring of this pulsar as part of an international, multi-telescope campaign. This includes an update of the pulsar's long-term timing and post-Keplerian parameters, new constraints on the pulsar's proper motion and corresponding Shklovskii kinematic correction, and new limits on the pulsar's geodetic precession as determined by monitoring for secular changes in the pulse profile. We also highlight prospects for future work, including an updated timeline on new relativistic tests following the introduction of MeerKAT observations. △ Less

Submitted 29 March, 2022; originally announced March 2022.

Comments: 11 pages, 3 figures, 1 table. Accepted for publication in the Proceedings of the 16th Marcel Grossmann meeting, July 5-10 2021, World Scientific

TRAPUM discovery of thirteen new pulsars in NGC 1851 using MeerKAT

Authors: A. Ridolfi, P. C. C. Freire, T. Gautam, S. M. Ransom, E. D. Barr, S. Buchner, M. Burgay, F. Abbate, V. Venkatraman Krishnan, L. Vleeschower, A. Possenti, B. W. Stappers, M. Kramer, W. Chen, P. V. Padmanabh, D. J. Champion, M. Bailes, L. Levin, E. F. Keane, R. P. Breton, M. Bezuidenhout, J. -M. Grießmeier, L. Künkel, Y. Men, F. Camilo , et al. (5 additional authors not shown)

Abstract: We report the discovery of 13 new pulsars in the globular cluster NGC 1851 by the TRAPUM Large Survey Project using the MeerKAT radio telescope. The discoveries consist of six isolated millisecond pulsars (MSPs) and seven binary pulsars, of which six are MSPs and one is mildly recycled. For all the pulsars, we present the basic kinematic, astrometric, and orbital parameters, where applicable, as w… ▽ More We report the discovery of 13 new pulsars in the globular cluster NGC 1851 by the TRAPUM Large Survey Project using the MeerKAT radio telescope. The discoveries consist of six isolated millisecond pulsars (MSPs) and seven binary pulsars, of which six are MSPs and one is mildly recycled. For all the pulsars, we present the basic kinematic, astrometric, and orbital parameters, where applicable, as well as their polarimetric properties, when these are measurable. Two of the binary MSPs (PSR J0514-4002D and PSR J0514-4002E) are in wide and extremely eccentric (e > 0.7) orbits with a heavy white dwarf and a neutron star as their companion, respectively. With these discoveries, NGC 1851 is now tied with M28 as the cluster with the third largest number of known pulsars (14). Its pulsar population shows remarkable similarities with that of M28, Terzan 5 and other clusters with comparable structural parameters. The newly-found pulsars are all located in the innermost regions of NGC 1851 and will likely enable, among other things, detailed studies of the cluster structure and dynamics. △ Less

Submitted 23 March, 2022; originally announced March 2022.

Comments: 12 pages, 6 figures, 3 tables. Accepted for publication in Astronomy & Astrophysics

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

Two New Black Widow Millisecond Pulsars In M28

Authors: Andrew Douglas, Prajwal Padmanabh, Scott Ransom, Alessandro Ridolfi, Paulo Freire, Vivek Venkatraman Krishnan, Ewan Barr, Cristina Pallanca, Mario Cadelano, Andrea Possenti, Ingrid Stairs, Jason Hessels, Megan DeCesar, Ryan Lynch, Matthew Bailes, Marta Burgay, David Champion, Ramesh Karuppusamy, Michael Kramer, Benjamin Stappers, Laila Vleeschower

Abstract: We report the discovery of two Black Widow millisecond pulsars in the globular cluster M28 with the MeerKAT telescope. PSR J1824$-$2452M (M28M) is a 4.78-ms pulsar in a $5.82\,$hour orbit and PSR J1824$-$2452N (M28N) is a 3.35-ms pulsar in a $4.76\,$hour orbit. Both pulsars have dispersion measures near $119.30\,$pc$\,$cm$^{-3}$ and have low mass companion stars ($\sim$$0.01-0.03\,$M$_\odot$), whi… ▽ More We report the discovery of two Black Widow millisecond pulsars in the globular cluster M28 with the MeerKAT telescope. PSR J1824$-$2452M (M28M) is a 4.78-ms pulsar in a $5.82\,$hour orbit and PSR J1824$-$2452N (M28N) is a 3.35-ms pulsar in a $4.76\,$hour orbit. Both pulsars have dispersion measures near $119.30\,$pc$\,$cm$^{-3}$ and have low mass companion stars ($\sim$$0.01-0.03\,$M$_\odot$), which do not cause strong radio eclipses or orbital variations. Including these systems, there are now five known black widow pulsars in M28. The pulsar searches were conducted as a part of an initial phase of MeerKAT's globular cluster census (within the TRAPUM Large Survey Project). These faint discoveries demonstrate the advantages of MeerKAT's survey sensitivity over previous searches and we expect to find additional pulsars in continued searches of this cluster. △ Less

Submitted 17 February, 2022; v1 submitted 26 January, 2022; originally announced January 2022.

Comments: Accepted to ApJ on 12/28/2021 2/17/2022: Added Acknowledgements, Software, and Facilities sections. Fixed missing 'Total Proper Motion' units in table 2

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

Coherent Search for Binary Pulsars across all Five Keplerian Parameters in Radio Observations using the template-bank algorithm

Authors: Vishnu Balakrishnan, David Champion, Ewan Barr, Michael Kramer, V. Venkatraman Krishnan, Ralph P. Eatough, Rahul Sengar, Matthew Bailes

Abstract: Relativistic binary pulsars orbiting white dwarfs and neutron stars have already provided excellent tests of gravity. However, despite observational efforts, a pulsar orbiting a black hole has remained elusive. One possible explanation is the extreme Doppler smearing caused by the pulsar's orbital motion which changes its apparent spin frequency during an observation. The classical solution to thi… ▽ More Relativistic binary pulsars orbiting white dwarfs and neutron stars have already provided excellent tests of gravity. However, despite observational efforts, a pulsar orbiting a black hole has remained elusive. One possible explanation is the extreme Doppler smearing caused by the pulsar's orbital motion which changes its apparent spin frequency during an observation. The classical solution to this problem has been to assume a constant acceleration or jerk for the entire observation. However, this assumption breaks down when the observation samples a large fraction of the orbit. This limits the length of search observations, and hence their sensitivity. This provides a strong motivation to develop techniques that can find compact binaries in longer observations. Here we present a GPU-based radio pulsar search pipeline that can perform a coherent search for binary pulsars by directly searching over three or five Keplerian parameters using the template-bank algorithm. We compare the sensitivity obtained from our pipeline with acceleration and jerk search pipelines for simulated pulsar-stellar-mass black hole binaries and observations of PSR J0737-3039A. We also discuss the computational feasibility of our pipeline for untargeted pulsar surveys and targeted searches. Our benchmarks indicate that circular orbit searches for P-BH binaries with spin-period P$_{\rm spin} \geq 20 \rm ms$ covering the 3-10 T$\mathrm{_{obs}}$ regime are feasible for the High Time Resolution Universe pulsar survey. Additionally, an elliptical orbit search in Globular clusters for P$_{\rm spin} \geq 20 \rm ms$ pulsars orbiting intermediate-mass black holes in the 5-10 T$\mathrm{_{obs}}$ regime is feasible for observations shorter than 2 hours with an eccentricity limit of 0.1. △ Less

Submitted 22 December, 2021; originally announced December 2021.

Comments: Accepted for publication by MNRAS