Showing 1–50 of 128 results for author: Desvignes, G

Galactic Centre Pulsars with the SKAO

Authors: F. Abbate, A. Carleo, S. Chatterjee, J. Cordes, P. B. Demorest, G. Desvignes, R. P. Eatough, E. Hackmann, Z. Hu, M. Kramer, J. Lazio, K. J. Lee, K. Liu, I. Rammala-Zitha, S. M. Ransom, G. Saowanit, L. Shao, P. Torne, R. Wharton, J. Wongphechauxsorn, W. Zhu

Abstract: The detection of a pulsar closely orbiting our Galaxy's supermassive black hole - Sagittarius~A* - is one of the ultimate prizes in pulsar astrophysics. The relativistic effects expected in such a system could far exceed those currently observable in compact binaries such as double neutron stars and pulsar white dwarfs. In addition, pulsars offer the opportunity to study the magneto-ionic properti… ▽ More The detection of a pulsar closely orbiting our Galaxy's supermassive black hole - Sagittarius~A* - is one of the ultimate prizes in pulsar astrophysics. The relativistic effects expected in such a system could far exceed those currently observable in compact binaries such as double neutron stars and pulsar white dwarfs. In addition, pulsars offer the opportunity to study the magneto-ionic properties of Earth's nearest galactic nucleus in unprecedented detail. For these reasons, and more, a multitude of pulsar searches of the Galactic Centre have been undertaken, with the outcome of just seven pulsar detections within a projected distance of 100\,pc from Sagittarius~A*. It is currently understood that a larger underlying population likely exists, but it is not until observations with the SKA have started that this population can be revealed. In this chapter, we look at important updates since the publication of the last SKAO science book and offer a focused view of observing strategies and likely outcomes with the updated SKAO design. △ Less

Submitted 3 July, 2026; originally announced July 2026.

Comments: Published in Advancing Astrophysics with the SKA II (AASKAII), 2026 (arXiv:2606.20366). Report-no:AASKAII/Abbate01. Advancing Astrophysics with the SKA II (AASKAII) outlines the transformative scientific advances that will be enabled by the SKA telescopes. An earlier version was published in The Open Journal of Astrophysics with arXiv ID: arXiv:2512.16155

Report number: AASKAII/Abbate01

Full-polarization millimeter wavelength variability of Sagittarius A* during the 2018 EHT campaign

Authors: Ezequiel Albentosa-Ruiz, Jasmin E. Washington, Nicola Marchili, Iván Martí-Vidal, Ciriaco Goddi, Maciek Wielgus, Alejandro Mus, Angelo Ricarte, Daniel P. Marrone, León D. S. Salas, Yuhei Iwata, Douglas F. Carlos, Alexandra J. Tetarenko, Kotaro Moriyama, Vedant Dhruv, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Uwe Bach, Anne-Kathrin Baczko, David Ball , et al. (250 additional authors not shown)

Abstract: Sagittarius A* (Srg A*), the supermassive black hole at the center of the Milky Way, provides a unique laboratory to study accretion dynamics and plasma processes near the event horizon. We investigated the variability and polarization properties of Srg A* using ALMA observations during the 2018 Event Horizon Telescope campaign. We analyzed high-cadence full-polarization light curves from ALMA at… ▽ More Sagittarius A* (Srg A*), the supermassive black hole at the center of the Milky Way, provides a unique laboratory to study accretion dynamics and plasma processes near the event horizon. We investigated the variability and polarization properties of Srg A* using ALMA observations during the 2018 Event Horizon Telescope campaign. We analyzed high-cadence full-polarization light curves from ALMA at millimeter wavelengths, performed time-series analysis, and investigated the temporal behavior during an X-ray flare observed by Chandra on 2018 April 24. The variability characteristics are compared with expectations from standard accretion flow models. We find low variability in total intensity ($σ/μ< 10\%$), but significantly higher variability in linear and circular polarization (~ 30% and ~ 50%, respectively). A time-series analysis reveals red-noise variability, with power spectral densities between -2 and -3 across all Stokes parameters. Polarized intensity shows stable intra-day timescales, while total intensity exhibits more variable timescales, suggesting distinct emission regions, with polarization likely arising from a coherent structure. On April 24, a statistically significant inter-band delay in polarized intensity coincides with a near-simultaneous X-ray and millimeter peak that deviates from the typical delayed flare scenario. This event also features enhanced millimeter variability and coherent polarization loop evolution. The observed simultaneity challenges standard models of transient synchrotron emission with cooling delays, favoring instead a scenario of continuous energy injection in an optically thin region. Our results offer new constraints on the physical mechanisms driving variability in Srg A*, and provide key observational input for refining theoretical models of accretion and plasma behavior in the vicinity of supermassive black holes. △ Less

Submitted 11 April, 2026; originally announced April 2026.

Journal ref: Astronomy & Astrophysics 708 (2026) A179

Detection of Cyclotron Absorption in the Radio Emission of GPM 1839-10

Authors: Yunpeng Men, Ewan Barr, Yuanhong Qu, Csanad Horvath, Jinchen Jiang, Gregory Desvignes, Natasha Hurley-Walker, Michael Kramer, Rui Luo, Samuel J. McSweeney, Jason Wu

Abstract: GPM 1839-10 is an intriguing long-period radio transient (LPT), distinguished by its activity spanning at least three decades and its highly unusual emission characteristics. These features include orthogonal polarization mode (OPM) switches, down-drifting sub-structures, and distinct linear-to-circular polarization conversion behaviors. In this work, we present follow-up observations utilizing th… ▽ More GPM 1839-10 is an intriguing long-period radio transient (LPT), distinguished by its activity spanning at least three decades and its highly unusual emission characteristics. These features include orthogonal polarization mode (OPM) switches, down-drifting sub-structures, and distinct linear-to-circular polarization conversion behaviors. In this work, we present follow-up observations utilizing the FAST telescope at L-band, yielding a total of seven detected radio pulses. We find a consistent association between OPM switches and a decrease in polarized intensity. This feature strongly supports the hypothesis that the OPM switches are generated by the incoherent summation of OPMs. Our measured Rotation Measures (RMs) are consistent with previous observations, indicating that the magneto-ionic environment is stable. If the source is in a binary system, such stability suggests it may host a weakly magnetized companion. Crucially, we firstly observe clear evidence of a cyclotron absorption feature in one radio pulse, a signature rarely observed in radio sources. This feature allows us to infer that the magnetic field strength at the absorption site has a lower limit of tens of Gauss, which is necessary for the phenomenon to occur. This characteristic can be explained in a scenario where GPM 1839-10 possesses a weakly magnetized companion star. △ Less

Submitted 14 February, 2026; originally announced February 2026.

Comments: 9 pages, 7 figures; Accepted for publication in SCPMA

Long-term timing of the relativistic binary PSR J1906+0746

Authors: L. Vleeschower, B. W. Stappers, M. J. Keith, G. Desvignes, P. C. C. Freire, M. Kramer, J. van Leeuwen, L. Levin, A. G. Lyne, I. H. Stairs, V. Venkatraman Krishnan, Y. Y. Wang

Abstract: We conducted a timing analysis of over 18 years of data on the young ($τ_{\rm c} = 112$ kyr, $P = 114$\,ms) relativistic binary PSR~J1906+0746, using six radio telescopes: Arecibo, FAST, Green Bank, Lovell, MeerKAT, and Nançay. This pulsar is known to orbit a compact high-mass companion with a period of 3.98\,hrs in a mildly eccentric orbit ($e = 0.085$). By combining all data and maintaining a co… ▽ More We conducted a timing analysis of over 18 years of data on the young ($τ_{\rm c} = 112$ kyr, $P = 114$\,ms) relativistic binary PSR~J1906+0746, using six radio telescopes: Arecibo, FAST, Green Bank, Lovell, MeerKAT, and Nançay. This pulsar is known to orbit a compact high-mass companion with a period of 3.98\,hrs in a mildly eccentric orbit ($e = 0.085$). By combining all data and maintaining a coherent timing solution over the full span, we obtained a more precise measurement of the advance of periastron, $\dotω = 7.5841(2)$\,$°$\,yr$^{-1}$, the Einstein delay, $γ= 4.59(2) \times 10^{-4}$\,s, and the secular change in orbital period, $\dot{P}_{\rm b} = -5.65(2) \times 10^{-13}$\,s\,s$^{-1}$. Assuming the validity of general relativity, we obtain a total mass of 2.6133(1)\,$M_{\odot}$ and component masses of 1.316(5) $M_{\odot}$ for the pulsar and 1.297(5) $M_{\odot}$ for the companion, consistent with a double neutron star system. However, when fitting for the secular change in the projected semi-major axis we obtain $\dot{x} = -1.8(6) \times 10^{-13}$\,s\,s$^{-1}$, the component masses are shifted by $\sim 3.5σ$, which is expected from the correlation of $\dot{x}$ and $γ$. The $\dot{x}$ has a similar magnitude to that observed in PSR~J1141$-$6545, which is due to spin-orbit coupling; if confirmed, it would indicate that, as in the latter system, the companion of the PSR~J1906+0746 system is a massive fast-rotating white dwarf formed before the pulsar. Additionally, we report and characterize a large glitch near MJD 56664, with a fractional frequency increase comparable to those observed in the Vela pulsar. △ Less

Submitted 5 February, 2026; originally announced February 2026.

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

Journal ref: Mon Not R Astron Soc (2026)

Locating the missing large-scale emission in the jet of M87* with short EHT baselines

Authors: Boris Georgiev, Paul Tiede, Sebastiano D. von Fellenberg, Michael Janssen, Iniyan Natarajan, Lindy Blackburn, Jongho Park, Erandi Chavez, Andrew T. West, Kotaro Moriyama, Jun Yi Koay, Hendrik Müller, Dhanya G. Nair, Avery E. Broderick, Maciek Wielgus, Kazunori Akiyama, Ezequiel Albentosa-Ruíz, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Uwe Bach, Anne-Kathrin Baczko , et al. (258 additional authors not shown)

Abstract: In Very-Long Baseline Interferometric arrays, nearly co-located stations probe the largest scales and typically cannot resolve the observed source. In the absence of large-scale structure, closure phases constructed with these stations are zero and, since they are independent of station-based errors, they can be used to probe data issues. Here, we show with an expansion about co-located stations,… ▽ More In Very-Long Baseline Interferometric arrays, nearly co-located stations probe the largest scales and typically cannot resolve the observed source. In the absence of large-scale structure, closure phases constructed with these stations are zero and, since they are independent of station-based errors, they can be used to probe data issues. Here, we show with an expansion about co-located stations, how these trivial closure phases become non-zero with brightness distribution on smaller scales than their short baseline would suggest. When applied to sources that are made up of a bright compact and large-scale diffuse component, the trivial closure phases directly measure the centroid relative to the compact source and higher-order image moments. We present a technique to measure these image moments with minimal model assumptions and validate it on synthetic Event Horizon Telescope (EHT) data. We then apply this technique to 2017 and 2018 EHT observations of M87* and find a weak preference for extended emission in the direction of the large-scale jet. We also apply it to 2021 EHT data and measure the source centroid about 1 mas northwest of the compact ring, consistent with the jet observed at lower frequencies. △ Less

Submitted 19 January, 2026; originally announced January 2026.

Comments: 15 pages, 8 figures, Submitted to Astronomy & Astrophysics

Journal ref: A&A 709, A26 (2026)

Galactic Centre Pulsars with the SKAO

Authors: F. Abbate, A. Carleo, S. Chatterjee, J. Cordes, P. B. Demorest, G. Desvignes, R. P. Eatough, E. Hackmann, Hu Z., M. Kramer, J. Lazio, K. J. Lee, K. Liu, I. Rammala-Zitha, S. M. Ransom, G. Saowanit, L. Shao, P. Torne, R. Wharton, J. Wongphechauxsorn, W. Zhu, The SKAO Pulsar Science Working Group

Abstract: The detection of a pulsar closely orbiting our Galaxy's supermassive black hole - Sagittarius A* - is one of the ultimate prizes in pulsar astrophysics. The relativistic effects expected in such a system could far exceed those currently observable in compact binaries such as double neutron stars and pulsar white dwarfs. In addition, pulsars offer the opportunity to study the magneto-ionic properti… ▽ More The detection of a pulsar closely orbiting our Galaxy's supermassive black hole - Sagittarius A* - is one of the ultimate prizes in pulsar astrophysics. The relativistic effects expected in such a system could far exceed those currently observable in compact binaries such as double neutron stars and pulsar white dwarfs. In addition, pulsars offer the opportunity to study the magneto-ionic properties of Earth's nearest galactic nucleus in unprecedented detail. For these reasons, and more, a multitude of pulsar searches of the Galactic Centre have been undertaken, with the outcome of just seven pulsar detections within a projected distance of 100 pc from Sagittarius A*. It is currently understood that a larger underlying population likely exists, but it is not until observations with the SKA have started that this population can be revealed. In this paper, we look at important updates since the publication of the last SKAO science book and offer a focused view of observing strategies and likely outcomes with the updated SKAO design. △ Less

Submitted 17 December, 2025; originally announced December 2025.

Comments: 12 pages, 6 Figures. Published in the Open Journal of Astrophysics in a special issue on pulsar science with SKAO

Journal ref: Open J. Astrophys. 8 (2025) 154252

Probing jet base emission of M87* with the 2021 Event Horizon Telescope observations

Authors: Saurabh, Hendrik Müller, Sebastiano D. von Fellenberg, Paul Tiede, Michael Janssen, Lindy Blackburn, Avery E. Broderick, Erandi Chavez, Boris Georgiev, Thomas P. Krichbaum, Kotaro Moriyama, Dhanya G. Nair, Iniyan Natarajan, Jongho Park, Andrew Thomas West, Maciek Wielgus, Kazunori Akiyama, Ezequiel Albentosa-Ruíz, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Uwe Bach , et al. (260 additional authors not shown)

Abstract: We investigate the presence and spatial characteristics of the jet base emission in M87* at 230 GHz, enabled by the enhanced uv coverage in the 2021 Event Horizon Telescope (EHT) observations. The addition of the 12-m Kitt Peak Telescope and NOEMA provides two key intermediate-length baselines to SMT and the IRAM 30-m, giving sensitivity to emission structures at scales of $\sim250~μ$as and… ▽ More We investigate the presence and spatial characteristics of the jet base emission in M87* at 230 GHz, enabled by the enhanced uv coverage in the 2021 Event Horizon Telescope (EHT) observations. The addition of the 12-m Kitt Peak Telescope and NOEMA provides two key intermediate-length baselines to SMT and the IRAM 30-m, giving sensitivity to emission structures at scales of $\sim250~μ$as and $\sim2500~μ$as (0.02 pc and 0.2 pc). Without these baselines, earlier EHT observations lacked the capability to constrain emission on large scales, where a "missing flux" of order $\sim1$ Jy is expected. To probe these scales, we analyzed closure phases, robust against station-based gain errors, and modeled the jet base emission using a simple Gaussian offset from the compact ring emission at separations $>100~μ$as. Our analysis reveals a Gaussian feature centered at ($Δ$RA $\approx320~μ$as, $Δ$Dec $\approx60~μ$as), a projected separation of $\approx5500$ AU, with a flux density of only $\sim60$ mJy, implying that most of the missing flux in previous studies must arise from larger scales. Brighter emission at these scales is ruled out, and the data do not favor more complex models. This component aligns with the inferred direction of the large-scale jet and is consistent with emission from the jet base. While our findings indicate detectable jet base emission at 230 GHz, coverage from only two intermediate baselines limits reconstruction of its morphology. We therefore treat the recovered Gaussian as an upper limit on the jet base flux density. Future EHT observations with expanded intermediate-baseline coverage will be essential to constrain the structure and nature of this component. △ Less

Submitted 1 December, 2025; originally announced December 2025.

Comments: 19 pages, 18 figures. Abstract shortened with respect to the manuscript. Accepted in Astronomy & Astrophysics

Journal ref: A&A, 706, A27 (2026)

A survey for radio pulsars and transients in the 10 pc region around Sgr A*

Authors: G. Desvignes, R. P. Eatough, Y. Men, F. Abbate, R. Karuppusamy, M. Kramer, K. Liu, L. Shao, P. Torne, R. S. Wharton

Abstract: Here we report on a new survey for pulsars and transients in the 10 pc region around Sgr A* using the Effelsberg radio telescope at frequencies between 4 to 8 GHz. Our calibrated full-Stokes data were searched for pulsars and transients using PulsarX, TransientX and PRESTO. Polarisation information is used in the scoring of the candidates. Our periodicity acceleration and jerk searches allowed us… ▽ More Here we report on a new survey for pulsars and transients in the 10 pc region around Sgr A* using the Effelsberg radio telescope at frequencies between 4 to 8 GHz. Our calibrated full-Stokes data were searched for pulsars and transients using PulsarX, TransientX and PRESTO. Polarisation information is used in the scoring of the candidates. Our periodicity acceleration and jerk searches allowed us to maintain good sensitivity towards binary pulsars in $\gtrsim$ 10-hr orbits. In addition we performed a dedicated search in linear polarisation for slow transients. While our searches yielded no new discovery beyond the redetection of the magnetar SGR J1745-2900, we report on a faint single pulse candidate in addition to several weak periodicity search candidates. After thoroughly assessing our survey's sensitivity, we determined that it is still not sensitive to a population of millisecond pulsars. Next generation radio interferometers can overcome the limitations of traditional single-dish pulsar searches of the Galactic Centre. △ Less

Submitted 11 July, 2025; originally announced July 2025.

Comments: 10 pages, 9 figures, submitted to A&A, comments welcome

Journal ref: A&A 706, A113 (2026)

Long-Term Astrometric Monitoring of the Galactic Center Magnetar PSR J1745--2900

Authors: Geoffrey C. Bower, Adam T. Deller, Paul B. Demorest, Jason Dexter, Andreas Brunthaler, Gregory Desvignes, Ralph P. Eatough, Heino Falcke, Ciriaco Goddi, Michael Kramer, F. Yusef-Zadeh

Abstract: We present new astrometric observations of the Galactic Center magnetar, PSR J1745-2900, with the Very Long Baseline Array (VLBA). Combined with previously published measurements in 10 epochs that spanned 477 days, the complete data set consists of 25 epochs and 41 independent measurements that span 1984 days. These data constrain the proper motion to an accuracy of $\lesssim 2\%$ and set an upper… ▽ More We present new astrometric observations of the Galactic Center magnetar, PSR J1745-2900, with the Very Long Baseline Array (VLBA). Combined with previously published measurements in 10 epochs that spanned 477 days, the complete data set consists of 25 epochs and 41 independent measurements that span 1984 days. These data constrain the proper motion to an accuracy of $\lesssim 2\%$ and set an upper limit on the absolute value of the magnetar's acceleration of $\lesssim (0.4, 0.2)\, {\rm mas\,y^{-2}}$ in the two celestial coordinates, consistent with the maximum value of $\sim 0.03\,{\rm mas\,y^{-2}}$ expected for an orbit around Sgr A*. Future measurements have the potential to detect the acceleration of PSR J1745-2900 due to Sgr A* should PSR J1745-2900 re-brighten. We consider several potential sources of systematic variations in the astrometric residuals after fitting for standard parameters, including refractive wander, changes in the structure of Sgr A*, and the presence of an unseen binary companion. While a stellar companion model can be fit to the astrometric data, pulse period measurements are inconsistent with that model. No changes in the apparent image size of the magnetar were detected over the duration of these observations, indicating a lack of change in the properties of the line-of-sight scattering during this period. We also show that the upper limit to the mean core shift of Sgr A* is consistent with expectations for a compact jet or symmetric accretion flow. △ Less

Submitted 2 June, 2025; originally announced June 2025.

Comments: Accepted for publication in ApJ

Origin of the ring ellipticity in the black hole images of M87*

Authors: Rohan Dahale, Ilje Cho, Kotaro Moriyama, Kaj Wiik, Paul Tiede, José L. Gómez, Chi-kwan Chan, Roman Gold, Vadim Y. Bernshteyn, Marianna Foschi, Britton Jeter, Hung-Yi Pu, Boris Georgiev, Abhishek V. Joshi, Alejandro Cruz-Osorio, Iniyan Natarajan, Avery E. Broderick, León D. S. Salas, Koushik Chatterjee, Kazunori Akiyama, Ezequiel Albentosa-Ruíz, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua , et al. (251 additional authors not shown)

Abstract: We investigate the origin of the elliptical ring structure observed in the images of the supermassive black hole M87*, aiming to disentangle contributions from gravitational, astrophysical, and imaging effects. Leveraging the enhanced capabilities of the Event Horizon Telescope (EHT) 2018 array, including improved $(u,v)$-coverage from the Greenland Telescope, we measure the ring's ellipticity usi… ▽ More We investigate the origin of the elliptical ring structure observed in the images of the supermassive black hole M87*, aiming to disentangle contributions from gravitational, astrophysical, and imaging effects. Leveraging the enhanced capabilities of the Event Horizon Telescope (EHT) 2018 array, including improved $(u,v)$-coverage from the Greenland Telescope, we measure the ring's ellipticity using five independent imaging methods, obtaining a consistent average value of $τ= 0.08_{-0.02}^{+0.03}$ with a position angle $ξ= 50.1_{-7.6}^{+6.2}$ degrees. To interpret this measurement, we compare against General Relativistic Magnetohydrodynamic (GRMHD) simulations spanning a wide range of physical parameters including thermal or non-thermal electron distribution function, spins, and ion-to-electron temperature ratios in both low and high-density regions. We find no statistically significant correlation between spin and ellipticity in GRMHD images. Instead, we identify a correlation between ellipticity and the fraction of non-ring emission, particularly in non-thermal models and models with higher jet emission. These results indicate that the ellipticity measured from the \m87 emission structure is consistent with that expected from simulations of turbulent accretion flows around black holes, where it is dominated by astrophysical effects rather than gravitational ones. Future high-resolution imaging, including space very long baseline interferometry and long-term monitoring, will be essential to isolate gravitational signatures from astrophysical effects. △ Less

Submitted 15 May, 2025; originally announced May 2025.

Comments: 18 pages, 13 figures

Journal ref: A&A 699, A279 (2025)

The putative center in NGC 1052

Authors: Anne-Kathrin Baczko, Matthias Kadler, Eduardo Ros, Christian M. Fromm, Maciek Wielgus, Manel Perucho, Thomas P. Krichbaum, Mislav Baloković, Lindy Blackburn, Chi-kwan Chan, Sara Issaoun, Michael Janssen, Luca Ricci, Kazunori Akiyama, Ezequiel Albentosa-Ruíz, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Uwe Bach, David Ball, Bidisha Bandyopadhyay, John Barrett , et al. (262 additional authors not shown)

Abstract: Many active galaxies harbor powerful relativistic jets, however, the detailed mechanisms of their formation and acceleration remain poorly understood. To investigate the area of jet acceleration and collimation with the highest available angular resolution, we study the innermost region of the bipolar jet in the nearby low-ionization nuclear emission-line region (LINER) galaxy NGC 1052. We combine… ▽ More Many active galaxies harbor powerful relativistic jets, however, the detailed mechanisms of their formation and acceleration remain poorly understood. To investigate the area of jet acceleration and collimation with the highest available angular resolution, we study the innermost region of the bipolar jet in the nearby low-ionization nuclear emission-line region (LINER) galaxy NGC 1052. We combined observations of NGC 1052 taken with VLBA, GMVA, and EHT over one week in the spring of 2017. For the first time, NGC 1052 was detected with the EHT, providing a size of the central region in-between both jet bases of 250 RS (Schwarzschild radii) perpendicular to the jet axes. This size estimate supports previous studies of the jets expansion profile which suggest two breaks of the profile at around 300 RS and 10000 RS distances to the core. Furthermore, we estimated the magnetic field to be 1.25 Gauss at a distance of 22 μas from the central engine by fitting a synchrotron-self absorption spectrum to the innermost emission feature, which shows a spectral turn-over at about 130 GHz. Assuming a purely poloidal magnetic field, this implies an upper limit on the magnetic field strength at the event horizon of 26000 Gauss, which is consistent with previous measurements. The complex, low-brightness, double-sided jet structure in NGC 1052 makes it a challenge to detect the source at millimeter (mm) wavelengths. However, our first EHT observations have demonstrated that detection is possible up to at least 230 GHz. This study offers a glimpse through the dense surrounding torus and into the innermost central region, where the jets are formed. This has enabled us to finally resolve this region and provide improved constraints on its expansion and magnetic field strength. △ Less

Submitted 15 January, 2025; originally announced January 2025.

Comments: 22 pages, 10 figures, published in A&A

Journal ref: A&A, 692, A205 (2024)

A multi-frequency study of sub-parsec jets with the Event Horizon Telescope

Authors: Jan Röder, Maciek Wielgus, Andrei P. Lobanov, Thomas P. Krichbaum, Dhanya G. Nair, Sang-Sung Lee, Eduardo Ros, Vincent L. Fish, Lindy Blackburn, Chi-kwan Chan, Sara Issaoun, Michael Janssen, Michael D. Johnson, Sheperd S. Doeleman, Geoffrey C. Bower, Geoffrey B. Crew, Remo P. J. Tilanus, Tuomas Savolainen, C. M. Violette Impellizzeri, Antxon Alberdi, Anne-Kathrin Baczko, José L. Gómez, Ru-Sen Lu, Georgios F. Paraschos, Efthalia Traianou , et al. (265 additional authors not shown)

Abstract: The 2017 observing campaign of the Event Horizon Telescope (EHT) delivered the first very long baseline interferometry (VLBI) images at the observing frequency of 230 GHz, leading to a number of unique studies on black holes and relativistic jets from active galactic nuclei (AGN). In total, eighteen sources were observed: the main science targets, Sgr A* and M87 along with various calibrators. We… ▽ More The 2017 observing campaign of the Event Horizon Telescope (EHT) delivered the first very long baseline interferometry (VLBI) images at the observing frequency of 230 GHz, leading to a number of unique studies on black holes and relativistic jets from active galactic nuclei (AGN). In total, eighteen sources were observed: the main science targets, Sgr A* and M87 along with various calibrators. We investigated the morphology of the sixteen AGN in the EHT 2017 data set, focusing on the properties of the VLBI cores: size, flux density, and brightness temperature. We studied their dependence on the observing frequency in order to compare it with the Blandford-Königl (BK) jet model. We modeled the source structure of seven AGN in the EHT 2017 data set using linearly polarized circular Gaussian components and collected results for the other nine AGN from dedicated EHT publications, complemented by lower frequency data in the 2-86 GHz range. Then, we studied the dependences of the VLBI core flux density, size, and brightness temperature on the frequency measured in the AGN host frame. We compared the observations with the BK jet model and estimated the magnetic field strength dependence on the distance from the central black hole. Our results indicate a deviation from the standard BK model, particularly in the decrease of the brightness temperature with the observing frequency. Either bulk acceleration of the jet material, energy transfer from the magnetic field to the particles, or both are required to explain the observations. △ Less

Submitted 9 January, 2025; originally announced January 2025.

Journal ref: A&A 695, A233 (2025)

Searches for signatures of ultra-light axion dark matter in polarimetry data of the European Pulsar Timing Array

Authors: N. K. Porayko, P. Usynina, J. Terol-Calvo, J. Martin Camalich, G. M. Shaifullah, A. Castillo, D. Blas, L. Guillemot, M. Peel, C. Tiburzi, K. Postnov, M. Kramer, J. Antoniadis, S. Babak, A. -S. Bak Nielsen, E. Barausse, C. G. Bassa, C. Blanchard, M. Bonetti, E. Bortolas, P. R. Brook, M. Burgay, R. N. Caballero, A. Chalumeau, D. J. Champion , et al. (52 additional authors not shown)

Abstract: Ultra-light axion-like particles (ALPs) can be a viable solution to the dark matter problem. The scalar field associated with ALPs, coupled to the electromagnetic field, acts as an active birefringent medium, altering the polarisation properties of light through which it propagates. In particular, oscillations of the axionic field induce monochromatic variations of the plane of linearly polarised… ▽ More Ultra-light axion-like particles (ALPs) can be a viable solution to the dark matter problem. The scalar field associated with ALPs, coupled to the electromagnetic field, acts as an active birefringent medium, altering the polarisation properties of light through which it propagates. In particular, oscillations of the axionic field induce monochromatic variations of the plane of linearly polarised radiation of astrophysical signals. The radio emission of millisecond pulsars provides an excellent tool to search for such manifestations, given their high fractional linear polarisation and negligible fluctuations of their polarisation properties. We have searched for the evidence of ALPs in the polarimetry measurements of pulsars collected and preprocessed for the European Pulsar Timing Array (EPTA) campaign. Focusing on the twelve brightest sources in linear polarisation, we searched for an astrophysical signal from axions using both frequentist and Bayesian statistical frameworks. For the frequentist analysis, which uses Lomb-Scargle periodograms at its core, no statistically significant signal has been found. The model used for the Bayesian analysis has been adjusted to accommodate multiple deterministic systematics that may be present in the data. A statistically significant signal has been found in the dataset of multiple pulsars with common frequency between $10^{-8}$ Hz and $2\times10^{-8}$ Hz, which can most likely be explained by the residual Faraday rotation in the terrestrial ionosphere. Strong bounds on the coupling constant $g_{aγ}$, in the same ballpark as other searches, have been obtained in the mass range between $6\times10^{-24}$ eV and $5\times10^{-21}$ eV. We conclude by discussing problems that can limit the sensitivity of our search for ultra-light axions in the polarimetry data of pulsars, and possible ways to resolve them. △ Less

Submitted 3 December, 2024; originally announced December 2024.

Comments: 24 pages, 15 figures

First Very Long Baseline Interferometry Detections at 870μm

Authors: Alexander W. Raymond, Sheperd S. Doeleman, Keiichi Asada, Lindy Blackburn, Geoffrey C. Bower, Michael Bremer, Dominique Broguiere, Ming-Tang Chen, Geoffrey B. Crew, Sven Dornbusch, Vincent L. Fish, Roberto García, Olivier Gentaz, Ciriaco Goddi, Chih-Chiang Han, Michael H. Hecht, Yau-De Huang, Michael Janssen, Garrett K. Keating, Jun Yi Koay, Thomas P. Krichbaum, Wen-Ping Lo, Satoki Matsushita, Lynn D. Matthews, James M. Moran , et al. (254 additional authors not shown)

Abstract: The first very long baseline interferometry (VLBI) detections at 870$μ$m wavelength (345$\,$GHz frequency) are reported, achieving the highest diffraction-limited angular resolution yet obtained from the surface of the Earth, and the highest-frequency example of the VLBI technique to date. These include strong detections for multiple sources observed on inter-continental baselines between telescop… ▽ More The first very long baseline interferometry (VLBI) detections at 870$μ$m wavelength (345$\,$GHz frequency) are reported, achieving the highest diffraction-limited angular resolution yet obtained from the surface of the Earth, and the highest-frequency example of the VLBI technique to date. These include strong detections for multiple sources observed on inter-continental baselines between telescopes in Chile, Hawaii, and Spain, obtained during observations in October 2018. The longest-baseline detections approach 11$\,$G$λ$ corresponding to an angular resolution, or fringe spacing, of 19$μ$as. The Allan deviation of the visibility phase at 870$μ$m is comparable to that at 1.3$\,$mm on the relevant integration time scales between 2 and 100$\,$s. The detections confirm that the sensitivity and signal chain stability of stations in the Event Horizon Telescope (EHT) array are suitable for VLBI observations at 870$μ$m. Operation at this short wavelength, combined with anticipated enhancements of the EHT, will lead to a unique high angular resolution instrument for black hole studies, capable of resolving the event horizons of supermassive black holes in both space and time. △ Less

Submitted 9 October, 2024; originally announced October 2024.

Comments: Corresponding author: S. Doeleman

Journal ref: The Astronomical Journal, Volume 168, Issue 3, id.130, 19 pp. 2024

Periodicity search in the timing of the 25 millisecond pulsars from the second data release of the European Pulsar Timing Array

Authors: Iuliana Nitu, Michael Keith, David Champion, Ismael Cognard, Gregory Desvignes, Lucas Guillemot, Yanjun Guo, Huanchen Hu, Jiwoong Jang, Jedrzej Jawor, Ramesh Karuppusamy, Evan Keane, Michael Kramer, Kristen Lackeos, Kuo Liu, Robert Main, Delphine Perrodin, Nataliya Porayko, Golam Shaifullah, Gilles Theureau

Abstract: In this work, we investigated the presence of strictly periodic, as well as quasi-periodic signals, in the timing of the 25 millisecond pulsars from the EPTA DR2 dataset. This is especially interesting in the context of the recent hints of a gravitational wave background in these data, and the necessary further study of red-noise timing processes, which are known to behave quasi-periodically in so… ▽ More In this work, we investigated the presence of strictly periodic, as well as quasi-periodic signals, in the timing of the 25 millisecond pulsars from the EPTA DR2 dataset. This is especially interesting in the context of the recent hints of a gravitational wave background in these data, and the necessary further study of red-noise timing processes, which are known to behave quasi-periodically in some normal pulsars. We used Bayesian timing models developed through the run_enterprise pipeline: a strict periodicity was modelled as the influence of a planetary companion on the pulsar, while a quasi-periodicity was represented as a Fourier-domain Gaussian process. We found that neither model would clearly improve the timing models of the 25 millisecond pulsars in this dataset. This implies that noise and parameter estimates are unlikely to be biased by the presence of a (quasi-)periodicity in the timing data. Nevertheless, the results for PSRs J1744--1134 and J1012+5307 suggest that the standard noise models for these pulsars may not be sufficient. We also measure upper limits for the projected masses of planetary companions around each of the 25 pulsars. The data of PSR J1909--3744 yielded the best mass limits, such that we constrained the 95-percentile to 2*10^{-4} Earth-masses (roughly the mass of the dwarf planet Ceres) for orbital periods between 5 d--17 yr. These are the best pulsar planet mass limits to date. △ Less

Submitted 19 August, 2024; originally announced August 2024.

Comments: Submitted for publication in MNRAS

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

Broadband Multi-wavelength Properties of M87 during the 2018 EHT Campaign including a Very High Energy Flaring Episode

Authors: J. C. Algaba, M. Balokovic, S. Chandra, W. Y. Cheong, Y. Z. Cui, F. D'Ammando, A. D. Falcone, N. M. Ford, M. Giroletti, C. Goddi, M. A. Gurwell, K. Hada, D. Haggard, S. Jorstad, A. Kaur, T. Kawashima, S. Kerby, J. Y. Kim, M. Kino, E. V. Kravchenko, S. S. Lee, R. S. Lu, S. Markoff, J. Michail, J. Neilsen , et al. (721 additional authors not shown)

Abstract: The nearby elliptical galaxy M87 contains one of the only two supermassive black holes whose emission surrounding the event horizon has been imaged by the Event Horizon Telescope (EHT). In 2018, more than two dozen multi-wavelength (MWL) facilities (from radio to gamma-ray energies) took part in the second M87 EHT campaign. The goal of this extensive MWL campaign was to better understand the physi… ▽ More The nearby elliptical galaxy M87 contains one of the only two supermassive black holes whose emission surrounding the event horizon has been imaged by the Event Horizon Telescope (EHT). In 2018, more than two dozen multi-wavelength (MWL) facilities (from radio to gamma-ray energies) took part in the second M87 EHT campaign. The goal of this extensive MWL campaign was to better understand the physics of the accreting black hole M87*, the relationship between the inflow and inner jets, and the high-energy particle acceleration. Understanding the complex astrophysics is also a necessary first step towards performing further tests of general relativity. The MWL campaign took place in April 2018, overlapping with the EHT M87* observations. We present a new, contemporaneous spectral energy distribution (SED) ranging from radio to very high energy (VHE) gamma-rays, as well as details of the individual observations and light curves. We also conduct phenomenological modelling to investigate the basic source properties. We present the first VHE gamma-ray flare from M87 detected since 2010. The flux above 350 GeV has more than doubled within a period of about 36 hours. We find that the X-ray flux is enhanced by about a factor of two compared to 2017, while the radio and millimetre core fluxes are consistent between 2017 and 2018. We detect evidence for a monotonically increasing jet position angle that corresponds to variations in the bright spot of the EHT image. Our results show the value of continued MWL monitoring together with precision imaging for addressing the origins of high-energy particle acceleration. While we cannot currently pinpoint the precise location where such acceleration takes place, the new VHE gamma-ray flare already presents a challenge to simple one-zone leptonic emission model approaches, and emphasises the need for combined image and spectral modelling. △ Less

Submitted 5 December, 2024; v1 submitted 24 April, 2024; originally announced April 2024.

Comments: 46 pages, 23 figures, accepted by Astronomy & Astrophysics on August. 29, 2024

Journal ref: A&A 692, A140 (2024)

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

Journal ref: ApJ 966 (2024) 46

Ordered magnetic fields around the 3C 84 central black hole

Authors: G. F. Paraschos, J. -Y. Kim, M. Wielgus, J. Röder, T. P. Krichbaum, E. Ros, I. Agudo, I. Myserlis, M. Moscibrodzka, E. Traianou, J. A. Zensus, L. Blackburn, C. -K. Chan, S. Issaoun, M. Janssen, M. D. Johnson, V. L. Fish, K. Akiyama, A. Alberdi, W. Alef, J. C. Algaba, R. Anantua, K. Asada, R. Azulay, U. Bach , et al. (258 additional authors not shown)

Abstract: 3C84 is a nearby radio source with a complex total intensity structure, showing linear polarisation and spectral patterns. A detailed investigation of the central engine region necessitates the use of VLBI above the hitherto available maximum frequency of 86GHz. Using ultrahigh resolution VLBI observations at the highest available frequency of 228GHz, we aim to directly detect compact structures a… ▽ More 3C84 is a nearby radio source with a complex total intensity structure, showing linear polarisation and spectral patterns. A detailed investigation of the central engine region necessitates the use of VLBI above the hitherto available maximum frequency of 86GHz. Using ultrahigh resolution VLBI observations at the highest available frequency of 228GHz, we aim to directly detect compact structures and understand the physical conditions in the compact region of 3C84. We used EHT 228GHz observations and, given the limited (u,v)-coverage, applied geometric model fitting to the data. We also employed quasi-simultaneously observed, multi-frequency VLBI data for the source in order to carry out a comprehensive analysis of the core structure. We report the detection of a highly ordered, strong magnetic field around the central, SMBH of 3C84. The brightness temperature analysis suggests that the system is in equipartition. We determined a turnover frequency of $ν_m=(113\pm4)$GHz, a corresponding synchrotron self-absorbed magnetic field of $B_{SSA}=(2.9\pm1.6)$G, and an equipartition magnetic field of $B_{eq}=(5.2\pm0.6)$G. Three components are resolved with the highest fractional polarisation detected for this object ($m_\textrm{net}=(17.0\pm3.9)$%). The positions of the components are compatible with those seen in low-frequency VLBI observations since 2017-2018. We report a steeply negative slope of the spectrum at 228GHz. We used these findings to test models of jet formation, propagation, and Faraday rotation in 3C84. The findings of our investigation into different flow geometries and black hole spins support an advection-dominated accretion flow in a magnetically arrested state around a rapidly rotating supermassive black hole as a model of the jet-launching system in the core of 3C84. However, systematic uncertainties due to the limited (u,v)-coverage, however, cannot be ignored. △ Less

Submitted 1 February, 2024; originally announced February 2024.

Comments: 15 pages, 6 figures, published in A&A

Journal ref: Issue: A&A Volume 682, February 2024; Article number: L3; Number of pages: 15

Radio Pulse Profile Evolution of Magnetar Swift J1818.0-1607

Authors: Rebecca Fisher, Elliot Butterworth, Kaustubh Rajwade, Ben Stappers, Gregory Desvignes, Ramesh Karuppusamy, Michael Kramer, Kuo Liu, Andrew Lyne, Mitchell Mickaliger, Benjamin Shaw, Patrick Weltevrede

Abstract: The shape and polarisation properties of the radio pulse profiles of radio-loud magnetars provide a unique opportunity to investigate their magnetospheric properties. Gaussian Process Regression analysis was used to investigate the variation in the total intensity shape of the radio pulse profiles of the magnetar Swift J1818.0-1607. The observed profile shape was found to evolve through three mode… ▽ More The shape and polarisation properties of the radio pulse profiles of radio-loud magnetars provide a unique opportunity to investigate their magnetospheric properties. Gaussian Process Regression analysis was used to investigate the variation in the total intensity shape of the radio pulse profiles of the magnetar Swift J1818.0-1607. The observed profile shape was found to evolve through three modes between MJDs 59104 and 59365. The times at which these transitions occurred coincided with changes in the amplitude of modulations in the spin-down rate. The amount of linear and circular polarisation was also found to vary significantly with time. Lomb-Scargle periodogram analysis of the spin-down rate revealed three possibly harmonically related frequencies. This could point to the magnetar experiencing seismic activity. However, no profile features exhibited significant periodicity, suggesting no simple correlations between the profile variability and fluctuations of the spin-down on shorter timescales within the modes. Overall, this implies the mode changes seen are a result of local magnetospheric changes, with other theories, such as precession, less able to explain these observations. △ Less

Submitted 23 January, 2024; originally announced January 2024.

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

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

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

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

Submitted 15 January, 2024; originally announced January 2024.

Comments: Accepted for publication in MNRAS

Quasi-periodic sub-pulse structure as a unifying feature for radio-emitting neutron stars

Authors: Michael Kramer, Kuo Liu, Gregory Desvignes, Ramesh Karuppusamy, Ben W. Stappers

Abstract: Magnetars are highly-magnetised rotating neutron stars that are predominantly observed as high-energy sources. Six of this class of neutron star are known to also emit radio emission, and magnetars are, thus, a favoured model for the origin for at least some of the Fast Radio Bursts (FRBs). If magnetars, or neutron stars in general, are indeed responsible, sharp empirical constraints on the mechan… ▽ More Magnetars are highly-magnetised rotating neutron stars that are predominantly observed as high-energy sources. Six of this class of neutron star are known to also emit radio emission, and magnetars are, thus, a favoured model for the origin for at least some of the Fast Radio Bursts (FRBs). If magnetars, or neutron stars in general, are indeed responsible, sharp empirical constraints on the mechanism producing radio emission are required. Here we report on the detection of polarised quasi-periodic sub-structure in the emission of all well-studied radio-detected magnetars. A correlation previously seen, relating sub-structure in pulsed emission of radio emitting neutron stars to their rotational period, is extended, and shown to now span more than six of orders of magnitude in pulse period. This behaviour is not only seen in magnetars but in members of all classes of radio-emitting rotating neutron stars, regardless of their evolutionary history, their power source or their inferred magnetic field strength. If magnetars are responsible for FRBs, it supports the idea of being able to infer underlying periods from sub-burst timescales in FRBs. △ Less

Submitted 22 November, 2023; originally announced November 2023.

Comments: published in Nature Astronomy, 48 pages (main text, methods, extended data figures & tables and supplementary information), in total 17 figures, 4 tables

PulsarX: a new pulsar searching package -I. A high performance folding program for pulsar surveys

Authors: Yunpeng Men, Ewan Barr, C. J. Clark, Emma Carli, Gregory Desvignes

Abstract: Pulsar surveys with modern radio telescopes are becoming increasingly computationally demanding. This is particularly true for wide field-of-view pulsar surveys with radio interferometers, and those conducted in real or quasi-real time. These demands result in data analysis bottlenecks that can limit the parameter space covered by the surveys and diminish their scientific return. In this paper, we… ▽ More Pulsar surveys with modern radio telescopes are becoming increasingly computationally demanding. This is particularly true for wide field-of-view pulsar surveys with radio interferometers, and those conducted in real or quasi-real time. These demands result in data analysis bottlenecks that can limit the parameter space covered by the surveys and diminish their scientific return. In this paper, we address the computational challenge of `candidate folding' in pulsar searching, presenting a novel, efficient approach designed to optimise the simultaneous folding of large numbers of pulsar candidates. We provide a complete folding pipeline appropriate for large-scale pulsar surveys including radio frequency interference (RFI) mitigation, dedispersion, folding and parameter optimization. By leveraging the Fast Discrete Dispersion Measure Transform (FDMT) algorithm proposed by Zackay et al. (2017), we have developed an optimized, and cache-friendly implementation that we term the pruned FDMT (pFDMT). The pFDMT approach efficiently reuses intermediate processing results and prunes the unused computation paths, resulting in a significant reduction in arithmetic operations. In addition, we propose a novel folding algorithm based on the Tikhonov-regularised least squares method (TLSM) that can improve the time resolution of the pulsar profile. We present the performance of its real-world application as an integral part of two major pulsar search projects conducted with the MeerKAT telescope: the MPIfR-MeerKAT Galactic Plane Survey (MMGPS) and the Transients and Pulsars with MeerKAT (TRAPUM) project. In our processing, for approximately 500 candidates, the theoretical number of dedispersion operations can be reduced by a factor of around 50 when compared to brute-force dedispersion, which scales with the number of candidates. △ Less

Submitted 5 September, 2023; originally announced September 2023.

Comments: Accepted for publication in A&A

Journal ref: A&A 679, A20 (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)

A search for pulsars around Sgr A* in the first Event Horizon Telescope dataset

Authors: Pablo Torne, Kuo Liu, Ralph P. Eatough, Jompoj Wongphechauxsorn, James M. Cordes, Gregory Desvignes, Mariafelicia De Laurentis, Michael Kramer, Scott M. Ransom, Shami Chatterjee, Robert Wharton, Ramesh Karuppusamy, Lindy Blackburn, Michael Janssen, Chi-kwan Chan, Geoffrey B. Crew, Lynn D. Matthews, Ciriaco Goddi, Helge Rottmann, Jan Wagner, Salvador Sanchez, Ignacio Ruiz, Federico Abbate, Geoffrey C. Bower, Juan J. Salamanca , et al. (261 additional authors not shown)

Abstract: The Event Horizon Telescope (EHT) observed in 2017 the supermassive black hole at the center of the Milky Way, Sagittarius A* (Sgr A*), at a frequency of 228.1 GHz ($λ$=1.3 mm). The fundamental physics tests that even a single pulsar orbiting Sgr A* would enable motivate searching for pulsars in EHT datasets. The high observing frequency means that pulsars - which typically exhibit steep emission… ▽ More The Event Horizon Telescope (EHT) observed in 2017 the supermassive black hole at the center of the Milky Way, Sagittarius A* (Sgr A*), at a frequency of 228.1 GHz ($λ$=1.3 mm). The fundamental physics tests that even a single pulsar orbiting Sgr A* would enable motivate searching for pulsars in EHT datasets. The high observing frequency means that pulsars - which typically exhibit steep emission spectra - are expected to be very faint. However, it also negates pulse scattering, an effect that could hinder pulsar detections in the Galactic Center. Additionally, magnetars or a secondary inverse Compton emission could be stronger at millimeter wavelengths than at lower frequencies. We present a search for pulsars close to Sgr A* using the data from the three most-sensitive stations in the EHT 2017 campaign: the Atacama Large Millimeter/submillimeter Array, the Large Millimeter Telescope and the IRAM 30 m Telescope. We apply three detection methods based on Fourier-domain analysis, the Fast-Folding-Algorithm and single pulse search targeting both pulsars and burst-like transient emission; using the simultaneity of the observations to confirm potential candidates. No new pulsars or significant bursts were found. Being the first pulsar search ever carried out at such high radio frequencies, we detail our analysis methods and give a detailed estimation of the sensitivity of the search. We conclude that the EHT 2017 observations are only sensitive to a small fraction ($\lesssim$2.2%) of the pulsars that may exist close to Sgr A*, motivating further searches for fainter pulsars in the region. △ Less

Submitted 29 August, 2023; originally announced August 2023.

Comments: 33 pages, 7 figures, 6 Tables. Accepted for publication in ApJ

Rotation measure variations in Galactic Centre pulsars

Authors: F. Abbate, A. Noutsos, G. Desvignes, R. S. Wharton, P. Torne, M. Kramer, R. P. Eatough, R. Karuppusamy, K. Liu, L. Shao, J. Wongphechauxsorn

Abstract: We report the results of an observational campaign using the Effelsberg 100-m telescope of the pulsars J1746$-$2849, J1746$-$2850, J1746$-$2856 and J1745$-$2912 located in the Central Molecular Zone (CMZ) close to the Galactic centre in order to study rotation measure (RM) variations. We report for the first time the RM value of PSR J1746$-$2850 to be $-12234 \pm 181$ rad m$^{-2}$. This pulsar sho… ▽ More We report the results of an observational campaign using the Effelsberg 100-m telescope of the pulsars J1746$-$2849, J1746$-$2850, J1746$-$2856 and J1745$-$2912 located in the Central Molecular Zone (CMZ) close to the Galactic centre in order to study rotation measure (RM) variations. We report for the first time the RM value of PSR J1746$-$2850 to be $-12234 \pm 181$ rad m$^{-2}$. This pulsar shows significant variations of RM of $300-400$ rad m$^{-2}$ over the course of months to years that suggest a strongly magnetized environment. The structure function analysis of the RM of PSR J1746$-$2850 revealed a steep power-law index of $1.87_{-0.3}^{+0.4}$ comparable to the value expected for isotropic turbulence. This pulsar also showed large dispersion measure (DM) variation of $\sim 50$ pc cm$^{-3}$ in an event lasting a few months where the RM increased by $\sim 200$ rad m$^{-2}$. The large difference in RM between PSR J1746$-$2849 and PSR J1746$-$2850 despite the small angular separation reveals the presence of a magnetic field of at least 70 $μ$G in the CMZ and can explain the lack of polarization in the radio images of the region. These results contribute to our understanding of the magnetic field in the CMZ and show similarities between the RM behaviours of these pulsars and some fast radio bursts (FRBs). △ Less

Submitted 6 July, 2023; originally announced July 2023.

Comments: Accepted for publication on Monthly Notices of the Royal Astronomical Society, 13 pages, 7 figures

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Comments: 13 figures, 15 pages, accepted

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

The second data release from the European Pulsar Timing Array 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

BlackHoleCam -- Testing general relativity with pulsars orbiting Sagittarius A*

Authors: Ralph P. Eatough, Gregory Desvignes, Kuo Liu, Robert S. Wharton, Aristedis Noutsos, Pablo Torne, Ramesh Karuppusamy, Lijing Shao, Michael Kramer, Heino Falcke, Luciano Rezzolla

Abstract: BlackHoleCam is a project funded by a European Research Council Synergy Grant to build a complete astrophysical description of nearby supermassive black holes by using a combination of radio imaging, pulsar observations, stellar astrometry and general relativistic magneto-hydrodynamic models. BlackHoleCam scientists are active partners of the Event Horizon Telescope Consortium. In this talk I will… ▽ More BlackHoleCam is a project funded by a European Research Council Synergy Grant to build a complete astrophysical description of nearby supermassive black holes by using a combination of radio imaging, pulsar observations, stellar astrometry and general relativistic magneto-hydrodynamic models. BlackHoleCam scientists are active partners of the Event Horizon Telescope Consortium. In this talk I will discuss the use of pulsars orbiting Sagittarius A* for tests of General Relativity, the current difficulties in detecting such sources, recent results from the Galactic Centre magnetar PSR J1745-2900 and how BlackHoleCam aims to search for undiscovered pulsars in the Galactic Centre. △ Less

Submitted 2 June, 2023; originally announced June 2023.

Comments: 4 pages, 3 figures, Proceedings of The Fifteenth Marcel Grossmann Meeting on General Relativity, University of Rome La Sapienza, Italy, 1 - 7 July 2018

Comparison of Polarized Radiative Transfer Codes used by the EHT Collaboration

Authors: Cora Prather, Jason Dexter, Monika Moscibrodzka, Hung-Yi Pu, Thomas Bronzwaer, Jordy Davelaar, Ziri Younsi, Charles F. Gammie, Roman Gold, George N. Wong, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Uwe Bach, Anne-Kathrin Baczko, David Ball, Mislav Baloković, John Barrett, Michi Bauböck, Bradford A. Benson, Dan Bintley , et al. (248 additional authors not shown)

Abstract: Interpretation of resolved polarized images of black holes by the Event Horizon Telescope (EHT) requires predictions of the polarized emission observable by an Earth-based instrument for a particular model of the black hole accretion system. Such predictions are generated by general relativistic radiative transfer (GRRT) codes, which integrate the equations of polarized radiative transfer in curve… ▽ More Interpretation of resolved polarized images of black holes by the Event Horizon Telescope (EHT) requires predictions of the polarized emission observable by an Earth-based instrument for a particular model of the black hole accretion system. Such predictions are generated by general relativistic radiative transfer (GRRT) codes, which integrate the equations of polarized radiative transfer in curved spacetime. A selection of ray-tracing GRRT codes used within the EHT collaboration is evaluated for accuracy and consistency in producing a selection of test images, demonstrating that the various methods and implementations of radiative transfer calculations are highly consistent. When imaging an analytic accretion model, we find that all codes produce images similar within a pixel-wise normalized mean squared error (NMSE) of 0.012 in the worst case. When imaging a snapshot from a cell-based magnetohydrodynamic simulation, we find all test images to be similar within NMSEs of 0.02, 0.04, 0.04, and 0.12 in Stokes I, Q, U , and V respectively. We additionally find the values of several image metrics relevant to published EHT results to be in agreement to much better precision than measurement uncertainties. △ Less

Submitted 23 July, 2026; v1 submitted 21 March, 2023; originally announced March 2023.

Comments: Accepted for publication in ApJ

Journal ref: ApJ, 950, 35, 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.

Precession of magnetars: dynamical evolutions and modulations on polarized electromagnetic waves

Authors: Yong Gao, Lijing Shao, Gregory Desvignes, David Ian Jones, Michael Kramer, Garvin Yim

Abstract: Magnetars are conjectured to be highly magnetized neutron stars (NSs). Strong internal magnetic field and elasticity in the crust may deform the stars and lead to free precession. We study the precession dynamics of triaxially-deformed NSs incorporating the near-field and the far-field electromagnetic torques. We obtain timing residuals for different NS geometries and torques. We also investigate… ▽ More Magnetars are conjectured to be highly magnetized neutron stars (NSs). Strong internal magnetic field and elasticity in the crust may deform the stars and lead to free precession. We study the precession dynamics of triaxially-deformed NSs incorporating the near-field and the far-field electromagnetic torques. We obtain timing residuals for different NS geometries and torques. We also investigate the polarized X-ray and radio signals from precessing magnetars. The modulations on the Stokes parameters are obtained for thermal X-rays emitted from the surface of magnetars. For radio signals, we apply the simple rotating vector model (RVM) to give the modulations on the position angle (PA) of the polarization. Our results are comprehensive, ready to be used to search for magnetar precession with timing data and polarizations of X-ray and radio emissions. Future observations of precessing magnetars will give us valuable information on the geometry and the strength of the strong magnetic fields, the emission geometry, as well as the equation of state (EoS) of NSs. △ Less

Submitted 30 November, 2022; originally announced November 2022.

Comments: 18 pages, 18 figures; accepted by MNRAS

Journal ref: MNRAS 519 (2023) 1080

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)

One pulsar, two white dwarfs, and a planet confirming the strong equivalence principle

Authors: Guillaume Voisin, G Luth, I Cognard, P Freire, N Wex, L Guillemot, G Desvignes, M Kramer, G Theureau, M Saillenfest

Abstract: The strong equivalence principle is a cornerstone of general relativity, tested with exquisite accuracy in the Solar system. However, tests in the strong-field regime require a compact object. Currently, PSR J0337+1715 is the unique millisecond pulsar found in a triple stellar system, orbiting two white dwarfs within an area comparable to the orbit of the Earth. This configuration offers the oppor… ▽ More The strong equivalence principle is a cornerstone of general relativity, tested with exquisite accuracy in the Solar system. However, tests in the strong-field regime require a compact object. Currently, PSR J0337+1715 is the unique millisecond pulsar found in a triple stellar system, orbiting two white dwarfs within an area comparable to the orbit of the Earth. This configuration offers the opportunity for a dramatic improvement over previous tests, provided that accurate and regular timing of the pulsar can be achieved. This also requires the development of a new timing model solving numerically the relativistic three-body problem with great accuracy. We report on the analysis of the high-quality dataset gathered on PSR J0337+1715 by the Nan{\c c}ay radiotelescope over the past 8 years. In particular, I will show how we could obtain the most stringent limit to-date on a potential violation of the strong equivalent principle in the strong field regime. I will also introduce preliminary resuts showing that the presence of a small planet in the system may explain a tiny residual signal so far unaccounted for, which if confirmed would make this system exceptionally rich. △ Less

Submitted 19 May, 2022; originally announced May 2022.

Comments: On behalf of: I. Cognard, P.C.C. Freire, N. Wex, L. Guillemot, G. Desvignes, M. Kramer, G. Theureau, M. Saillenfest

Journal ref: 2022 Gravitation session of the 56th Rencontres de Moriond, Jan 2022, La Thuile, Italy

Long term radio and X-ray evolution of the magnetar Swift J1818.0-1607

Authors: K. M. Rajwade, B. W. Stappers, A. G. Lyne, B. Shaw, M. B. Mickaliger, K. Liu, M. Kramer, G. Desvignes, R. Karuppusamy, T. Enoto, T. Güver, Chin-Ping Hu, M. P. Surnis

Abstract: We report on the the long term monitoring campaign of the seemingly youngest magnetar Swift~J1818.0-1607 at radio and X-ray wavelengths over a span of one year. We obtained a coherent timing solution for the magnetar over the same time span. The frequency derivative of the magnetar shows systematic variation with the values oscillating about a mean value of $-$1.37$\times$10$^{-11}$ Hz s$^{-1}$. T… ▽ More We report on the the long term monitoring campaign of the seemingly youngest magnetar Swift~J1818.0-1607 at radio and X-ray wavelengths over a span of one year. We obtained a coherent timing solution for the magnetar over the same time span. The frequency derivative of the magnetar shows systematic variation with the values oscillating about a mean value of $-$1.37$\times$10$^{-11}$ Hz s$^{-1}$. The magnitude of the variation in the frequency derivative reduces with time before converging on the mean value. We were able to identify four states in the spin-frequency derivative that were quantified by the amount of modulation about the mean value and the transition between these states seem to be correlated with the change in the radio emission of the magnetar while no correlation is seen in the average radio profile variability on a shorter timescale (days). The 0.5--12 keV X-ray flux shows a monotonic decrease that can be attributed to thermal emission from a hot-spot on the surface of the neutron star that is reducing in size. Such decrease is consistent with what is seen in other magnetars. The potential correlation between the radio emission mode and the behaviour of the spin-down rate hints to a global change in the magnetopshere of the magnetar akin to the correlation seen in a subset of mode-changing radio pulsars and suggests a physical link between the two sub-populations. △ Less

Submitted 15 February, 2022; originally announced February 2022.

Comments: 9 pages, 6 figures, 2 tables, supplementary material. Accepted for publication in MNRAS

Submillimeter pulsations from the magnetar XTE J1810-197

Authors: Pablo Torne, Graham Bell, Dan Bintley, Gregory Desvignes, David Berry, Jessica T. Dempsey, Paul T. P. Ho, Harriet Parsons, Ralph P. Eatough, Ramesh Karuppusamy, Michael Kramer, Carsten Kramer, Kuo Liu, Gabriel Paubert, Miguel Sanchez-Portal, Karl F. Schuster

Abstract: We present the first detection of pulsations from a neutron star in the submillimeter range. The source is the magnetar XTE J1810-197, observed with the James Clerk Maxwell Telescope (JCMT) on 2020 February 27, 2020 July 9 and 2021 May 15. XTE J1810-197 is detected at 353 GHz ($λ=0.85\,$mm) in the three epochs, but not detected in the simultaneously-observed band at 666 GHz ($λ=0.45\,$mm). We meas… ▽ More We present the first detection of pulsations from a neutron star in the submillimeter range. The source is the magnetar XTE J1810-197, observed with the James Clerk Maxwell Telescope (JCMT) on 2020 February 27, 2020 July 9 and 2021 May 15. XTE J1810-197 is detected at 353 GHz ($λ=0.85\,$mm) in the three epochs, but not detected in the simultaneously-observed band at 666 GHz ($λ=0.45\,$mm). We measure an averaged flux density at 353 GHz of 6.7$\pm$1.0, 4.0$\pm$0.6, and 1.3$\pm$0.3 mJy and set 3$σ$ flux density upper limits at 666 GHz of 11.3, 4.7 and 4.3 mJy, at each of the three observing epochs, respectively. Combining close-in-time observations with the Effelsberg 100m and IRAM 30m telescopes covering non-contiguously from 6 to 225 GHz (5.0 cm$>λ>$1.33 mm), we investigate the spectral shape and frequency range of a potential spectral turn-up predicted by some pulsar radio emission models. The results demonstrate that the beamed radio emission from neutron stars can extend into the submillimeter regime, but are inconclusive on the existence and location of a potential spectral turn-up within the covered frequency range. The observed properties of the submillimeter emission resemble those of the longer wavelengths, and support a coherent mechanism for the production of pulsations at 353 GHz. △ Less

Submitted 19 January, 2022; originally announced January 2022.

Comments: 14 pages, 2 tables, 5 figures. Accepted for publication in ApJ Letters

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

Strong-field Gravity Tests with the Double Pulsar

Authors: M. Kramer, I. H. Stairs, R. N. Manchester, N. Wex, A. T. Deller, W. A. Coles, M. Ali, M. Burgay, F. Camilo, I. Cognard, T. Damour, G. Desvignes, R. D. Ferdman, P. C. C. Freire, S. Grondin, L. Guillemot, G. B. Hobbs, G. Janssen, R. Karuppusamy, D. R. Lorimer, A. G. Lyne, J. W. McKee, M. McLaughlin, L. E. Muench, B. B. P. Perera , et al. (5 additional authors not shown)

Abstract: Continued observations of the Double Pulsar, PSR J0737-3039A/B, consisting of two radio pulsars (A and B) that orbit each other with a period of 2.45hr in a mildly eccentric (e=0.088) binary system, have led to large improvements in the measurement of relativistic effects in this system. With a 16-yr data span, the results enable precision tests of theories of gravity for strongly self-gravitating… ▽ More Continued observations of the Double Pulsar, PSR J0737-3039A/B, consisting of two radio pulsars (A and B) that orbit each other with a period of 2.45hr in a mildly eccentric (e=0.088) binary system, have led to large improvements in the measurement of relativistic effects in this system. With a 16-yr data span, the results enable precision tests of theories of gravity for strongly self-gravitating bodies and also reveal new relativistic effects that have been expected but are now observed for the first time. These include effects of light propagation in strong gravitational fields which are currently not testable by any other method. We observe retardation and aberrational light-bending that allow determination of the pulsar's spin direction. In total, we have detected seven post-Keplerian (PK) parameters, more than for any other binary pulsar. For some of these effects, the measurement precision is so high that for the first time we have to take higher-order contributions into account. These include contributions of A's effective mass loss (due to spin-down) to the observed orbital period decay, a relativistic deformation of the orbit, and effects of the equation of state of super-dense matter on the observed PK parameters via relativistic spin-orbit coupling. We discuss the implications of our findings, including those for the moment of inertia of neutron stars. We present the currently most precise test of general relativity's (GR's) quadrupolar description of gravitational waves, validating GR's prediction at a level of $1.3 \times 10^{-4}$ (95% conf.). We demonstrate the utility of the Double Pulsar for tests of alternative theories by focusing on two specific examples and discuss some implications for studies of the interstellar medium and models for the formation of the Double Pulsar. Finally, we provide context to other types of related experiments and prospects for the future. △ Less

Submitted 14 December, 2021; v1 submitted 13 December, 2021; originally announced December 2021.

Comments: 56 pages, 23 Figures. Published by Physical Review X. Uploaded ArXiv version is authors' (pre-proof) version with abbreviated abstract. Resubmission for updated acknowledgement information. For final published version see https://link.aps.org/doi/10.1103/PhysRevX.11.041050

Journal ref: Phys. Rev. X 11, 041050 (2021)

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

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

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

Submitted 9 November, 2021; originally announced November 2021.

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

Event Horizon Telescope observations of the jet launching and collimation in Centaurus A

Authors: Michael Janssen, Heino Falcke, Matthias Kadler, Eduardo Ros, Maciek Wielgus, Kazunori Akiyama, Mislav Baloković, Lindy Blackburn, Katherine L. Bouman, Andrew Chael, Chi-kwan Chan, Koushik Chatterjee, Jordy Davelaar, Philip G. Edwards, Christian M. Fromm, José L. Gómez, Ciriaco Goddi, Sara Issaoun, Michael D. Johnson, Junhan Kim, Jun Yi Koay, Thomas P. Krichbaum, Jun Liu, Elisabetta Liuzzo, Sera Markoff , et al. (215 additional authors not shown)

Abstract: Very-long-baseline interferometry (VLBI) observations of active galactic nuclei at millimeter wavelengths have the power to reveal the launching and initial collimation region of extragalactic radio jets, down to $10-100$ gravitational radii ($r_g=GM/c^2$) scales in nearby sources. Centaurus A is the closest radio-loud source to Earth. It bridges the gap in mass and accretion rate between the supe… ▽ More Very-long-baseline interferometry (VLBI) observations of active galactic nuclei at millimeter wavelengths have the power to reveal the launching and initial collimation region of extragalactic radio jets, down to $10-100$ gravitational radii ($r_g=GM/c^2$) scales in nearby sources. Centaurus A is the closest radio-loud source to Earth. It bridges the gap in mass and accretion rate between the supermassive black holes (SMBHs) in Messier 87 and our galactic center. A large southern declination of $-43^{\circ}$ has however prevented VLBI imaging of Centaurus A below $λ1$cm thus far. Here, we show the millimeter VLBI image of the source, which we obtained with the Event Horizon Telescope at $228$GHz. Compared to previous observations, we image Centaurus A's jet at a tenfold higher frequency and sixteen times sharper resolution and thereby probe sub-lightday structures. We reveal a highly-collimated, asymmetrically edge-brightened jet as well as the fainter counterjet. We find that Centaurus A's source structure resembles the jet in Messier 87 on ${\sim}500r_g$ scales remarkably well. Furthermore, we identify the location of Centaurus A's SMBH with respect to its resolved jet core at $λ1.3$mm and conclude that the source's event horizon shadow should be visible at THz frequencies. This location further supports the universal scale invariance of black holes over a wide range of masses. △ Less

Submitted 5 November, 2021; originally announced November 2021.

Comments: 27 pages, 9 figures. This is a post-peer-review, pre-copyedit version of an article published in Nature Astronomy. The final authenticated version is available online at: https://doi.org/10.1038/s41550-021-01417-w

Journal ref: Nature Astronomy, July 2021, Volume 5, p. 1017-1028

Radio and X-ray observations of giant pulses from XTE J1810-197

Authors: M. Caleb, K. Rajwade, G. Desvignes, B. W. Stappers, A. G. Lyne, P. Weltevrede, M. Kramer, L. Levin, M. Surnis

Abstract: We present the results of two years of radio and X-ray monitoring of the magnetar XTE J1810$-$197 since the radio re-activation in late 2018. Single pulse analysis of radio observations from the Lovell and MkII telescopes at 1564 MHz and the Effelsberg telescope at 6 GHz has resulted in the detection of a total of 91 giant pulses (GPs) between MJDs 58858 and 59117. These GPs appear to be confined… ▽ More We present the results of two years of radio and X-ray monitoring of the magnetar XTE J1810$-$197 since the radio re-activation in late 2018. Single pulse analysis of radio observations from the Lovell and MkII telescopes at 1564 MHz and the Effelsberg telescope at 6 GHz has resulted in the detection of a total of 91 giant pulses (GPs) between MJDs 58858 and 59117. These GPs appear to be confined to two specific phase ranges (0.473 <= φ<= 0.502$ and 0.541 <= φ<= 0.567). We also observe that the first detection of GP emission corresponds to a minimum in the spin-down rate. Simultaneous radio and X-ray observations were performed on MJDs 59009 and 59096. The 0.5-10 keV X-ray spectrum from NICER is well characterised by a two component blackbody model that can be interpreted as two hot spots on the polar cap of the neutron star. The blackbody temperature decreases with time, consistent with the previous outburst, while the change in the pulsed fraction does not follow the same trend as was seen in the previous outburst. The radio and X-ray flux of XTE J1810-197 are correlated during the initial phase of the outburst (MJD 58450 - MJD 58550) and an increase in the radio flux is observed later that may be correlated to the onset of GPs. We argue that the disparity in the evolution of the current outburst compared to the previous one can be attributed to a change in geometry of the neutron star. △ Less

Submitted 2 November, 2021; originally announced November 2021.

Comments: Accepted for publication in MNRAS

The Variability of the Black-Hole Image in M87 at the Dynamical Time Scale

Authors: Kaushik Satapathy, Dimitrios Psaltis, Feryal Ozel, Lia Medeiros, Sean T. Dougall, Chi-kwan Chan, Maciek Wielgus, Ben S. Prather, George N. Wong, Charles F. Gammie, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Anne-Kathrin Baczko, David R. Ball, Mislav Baloković, John Barrett, Bradford A. Benson, Dan Bintley, Lindy Blackburn, Raymond Blundell , et al. (213 additional authors not shown)

Abstract: The black-hole images obtained with the Event Horizon Telescope (EHT) are expected to be variable at the dynamical timescale near their horizons. For the black hole at the center of the M87 galaxy, this timescale (5-61 days) is comparable to the 6-day extent of the 2017 EHT observations. Closure phases along baseline triangles are robust interferometric observables that are sensitive to the expect… ▽ More The black-hole images obtained with the Event Horizon Telescope (EHT) are expected to be variable at the dynamical timescale near their horizons. For the black hole at the center of the M87 galaxy, this timescale (5-61 days) is comparable to the 6-day extent of the 2017 EHT observations. Closure phases along baseline triangles are robust interferometric observables that are sensitive to the expected structural changes of the images but are free of station-based atmospheric and instrumental errors. We explored the day-to-day variability in closure phase measurements on all six linearly independent non-trivial baseline triangles that can be formed from the 2017 observations. We showed that three triangles exhibit very low day-to-day variability, with a dispersion of $\sim3-5^\circ$. The only triangles that exhibit substantially higher variability ($\sim90-180^\circ$) are the ones with baselines that cross visibility amplitude minima on the $u-v$ plane, as expected from theoretical modeling. We used two sets of General Relativistic magnetohydrodynamic simulations to explore the dependence of the predicted variability on various black-hole and accretion-flow parameters. We found that changing the magnetic field configuration, electron temperature model, or black-hole spin has a marginal effect on the model consistency with the observed level of variability. On the other hand, the most discriminating image characteristic of models is the fractional width of the bright ring of emission. Models that best reproduce the observed small level of variability are characterized by thin ring-like images with structures dominated by gravitational lensing effects and thus least affected by turbulence in the accreting plasmas. △ Less

Submitted 1 November, 2021; originally announced November 2021.

Comments: Accepted for Publication in ApJ

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

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

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

Submitted 25 October, 2021; originally announced October 2021.

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

Multi-epoch searches for relativistic binary pulsars and fast transients in the Galactic Centre

Authors: R. P. Eatough, P. Torne, G. Desvignes, M. Kramer, R. Karuppusamy, B. Klein, L. G. Spitler, K. J. Lee, D. J. Champion, K. Liu, R. S. Wharton, L. Rezzolla, H. Falcke

Abstract: The high stellar density in the central parsecs around the Galactic Centre makes it a seemingly favourable environment for finding relativistic binary pulsars. These include pulsars orbiting other neutron stars, stellar-mass black holes or the central supermassive black hole, Sagittarius A*. Here we present multi-epoch pulsar searches of the Galactic Centre at four observing frequencies, (4.85, 8.… ▽ More The high stellar density in the central parsecs around the Galactic Centre makes it a seemingly favourable environment for finding relativistic binary pulsars. These include pulsars orbiting other neutron stars, stellar-mass black holes or the central supermassive black hole, Sagittarius A*. Here we present multi-epoch pulsar searches of the Galactic Centre at four observing frequencies, (4.85, 8.35, 14.6 18.95) GHz, using the Effelsberg 100-m radio telescope. Observations were conducted one year prior to the discovery of, and during monitoring observations of, the Galactic Centre magnetar PSR J1745-2900. Our data analysis features acceleration searches on progressively shorter time series to maintain sensitivity to relativistic binary pulsars. The multi-epoch observations increase the likelihood of discovering transient or nulling pulsars, or ensure orbital phases are observed at which acceleration search methods work optimally. In ~147 h of separate observations, no previously undiscovered pulsars have been detected. Through calibration observations, we conclude this might be due to insufficient instantaneous sensitivity; caused by the intense continuum emission from the Galactic Centre, its large distance and, at higher frequencies, the aggregate effect of steep pulsar spectral indices and atmospheric contributions to the system temperature. Additionally we find that for millisecond pulsars in wide circular orbits ~<800 d around Sagittarius A*, linear acceleration effects cannot be corrected in deep observations (9 h) with existing software tools. Pulsar searches of the Galactic Centre with the next generation of radio telescopes - such as MeerKat, ngVLA and SKA1-mid - will have improved chances of uncovering this elusive population. △ Less

Submitted 11 August, 2021; originally announced August 2021.

Comments: 17 pages, 12 figures, Accepted for publication in Monthly Notices of the Royal Astronomical Society