Showing 1–50 of 116 results for author: Karuppusamy, R

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

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

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

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

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)

Variable Scintillation Arcs of Millisecond Pulsars observed with the Large European Array for Pulsars

Authors: R. A. Main, J. Antoniadis, S. Chen, I. Cognard, H. Hu, J. Jang, R. Karuppusamy, M. Kramer, K. Liu, Y. Liu, G. Mall, J. W. McKee, M. B. Mickaliger, D. Perrodin, S. A. Sanidas, B. W. Stappers, T. Sprenger, O. Wucknitz, C. G. Bassa, M. Burgay, R. Concu, M. Gaikwad, G. H. Janssen, K. J. Lee, A. Melis , et al. (4 additional authors not shown)

Abstract: We present the first large sample of scintillation arcs in millisecond pulsars, analysing 12 sources observed with the Large European Array for Pulsars (LEAP), and the Effelsberg 100\,m telescope. We estimate the delays from multipath propagation, measuring significant correlated changes in scattering timescales over a 10-year timespan. Many sources show compact concentrations of power in the seco… ▽ More We present the first large sample of scintillation arcs in millisecond pulsars, analysing 12 sources observed with the Large European Array for Pulsars (LEAP), and the Effelsberg 100\,m telescope. We estimate the delays from multipath propagation, measuring significant correlated changes in scattering timescales over a 10-year timespan. Many sources show compact concentrations of power in the secondary spectrum, which in PSRs J0613$-$0200 and J1600$-$3053 can be tracked between observations, and are consistent with compact scattering at fixed angular positions. Other sources such as PSRs J1643$-$1224 and J0621+1002 show diffuse, asymmetric arcs which are likely related to phase-gradients across the scattering screen. PSR B1937+21 shows at least three distinct screens which dominate at different times and evidence of varying screen axes or multi-screen interactions. We model annual and orbital arc curvature variations in PSR J0613$-$0200, providing a measurement of the longitude of ascending node, resolving the sense of the orbital inclination, where our best fit model is of a screen with variable axis of anisotropy over time, corresponding to changes in the scattering of the source. Unmodeled variations of the screen's axis of anisotropy are likely to be a limiting factor in determining orbital parameters with scintillation, requiring careful consideration of variable screen properties, or independent VLBI measurements. Long-term scintillation studies such as this serve as a complementary tool to pulsar timing, to measure a source of correlated noise for pulsar timing arrays, solve pulsar orbits, and to understand the astrophysical origin of scattering screens. △ Less

Submitted 23 June, 2023; originally announced June 2023.

Comments: Submitted to MNRAS, first revision. 18 pages, 16 figures

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: Ben S. 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 21 March, 2023; originally announced March 2023.

Comments: Accepted for publication in ApJ

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

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

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

Submitted 16 December, 2022; originally announced December 2022.

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

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

Detection of quasi-periodic micro-structure in three millisecond pulsars with the Large European Array for Pulsars

Authors: K. Liu, J. Antoniadis, C. G. Bassa, S. Chen, I. Cognard, M. Gaikwad, H. Hu, J. Jang, G. H. Janssen, R. Karuppusamy, M. Kramer, K. J. Lee, R. A. Main, G. Mall, J. W. McKee, M. B. Mickaliger, D. Perrodin, S. A. Sanidas, B. W. Stappers, L. Wang, W. W. Zhu, M. Burgay, R. Concu, A. Corongiu, A. Melis , et al. (2 additional authors not shown)

Abstract: We report on the detection of quasi-periodic micro-structure in three millisecond pulsars (MSPs), PSRs J1022+1001, J2145-0750 and J1744-1134, using high time resolution data acquired with the Large European Array for Pulsars at a radio frequency of 1.4 GHz. The occurrence rate of quasi-periodic micro-structure is consistent among pulses with different peak flux densities. Using an auto-correlation… ▽ More We report on the detection of quasi-periodic micro-structure in three millisecond pulsars (MSPs), PSRs J1022+1001, J2145-0750 and J1744-1134, using high time resolution data acquired with the Large European Array for Pulsars at a radio frequency of 1.4 GHz. The occurrence rate of quasi-periodic micro-structure is consistent among pulses with different peak flux densities. Using an auto-correlation analysis, we measure the periodicity and width of the micro-structure in these three pulsars. The detected micro-structure from PSRs J1022+1001 and J1744-1134 is often highly linearly polarised. In PSR J1022+1001, the linear polarisation position angles of micro-structure pulses are in general flat with a small degree of variation. Using these results, we further examine the frequency and rotational period dependency of micro-structure properties established in previous work, along with the angular beaming and temporal modulation models that explains the appearance of micro-structure. We also discuss a possible link of micro-structure to the properties of some of the recently discovered fast radio bursts which exhibit a very similar emission morphology. △ Less

Submitted 20 June, 2022; originally announced June 2022.

Comments: 8 pages, 6 figures, 1 table, published in MNRAS

A burst storm from the repeating FRB 20200120E in an M81 globular cluster

Authors: K. Nimmo, J. W. T. Hessels, M. P. Snelders, R. Karuppusamy, D. M. Hewitt, F. Kirsten, B. Marcote, U. Bach, A. Bansod, E. D. Barr, J. Behrend, V. Bezrukovs, S. Buttaccio, R. Feiler, M. P. Gawroński, M. Lindqvist, A. Orbidans, W. Puchalska, N. Wang, T. Winchen, P. Wolak, J. Wu, J. Yuan

Abstract: The repeating fast radio burst (FRB) source FRB 20200120E is exceptional because of its proximity and association with a globular cluster. Here we report $60$ bursts detected with the Effelsberg telescope at 1.4 GHz. We observe large variations in the burst rate, and report the first FRB 20200120E `burst storm', where the source suddenly became active and 53 bursts (fluence $\geq 0.04$ Jy ms) occu… ▽ More The repeating fast radio burst (FRB) source FRB 20200120E is exceptional because of its proximity and association with a globular cluster. Here we report $60$ bursts detected with the Effelsberg telescope at 1.4 GHz. We observe large variations in the burst rate, and report the first FRB 20200120E `burst storm', where the source suddenly became active and 53 bursts (fluence $\geq 0.04$ Jy ms) occurred within only 40 minutes. We find no strict periodicity in the burst arrival times, nor any evidence for periodicity in the source's activity between observations. The burst storm shows a steep energy distribution (power-law index $α= 2.39\pm0.12$) and a bi-modal wait-time distribution, with log-normal means of 0.94$^{+0.07}_{-0.06}$ s and 23.61$^{+3.06}_{-2.71}$ s. We attribute these wait-time distribution peaks to a characteristic event timescale and pseudo-Poisson burst rate, respectively. The secondary wait-time peak at $\sim1$ s is $\sim50\times$ longer than the $\sim24$ ms timescale seen for both FRB 20121102A and FRB 20201124A -- potentially indicating a larger emission region, or slower burst propagation. FRB 20200120E shows order-of-magnitude lower burst durations and luminosities compared with FRB 20121102A and FRB 20201124A. Lastly, in contrast to FRB 20121102A, which has observed dispersion measure (DM) variations of $Δ{\rm DM} >1$ pc cm$^{-3}$ on month-to-year timescales, we determine that FRB 20200120E's DM has remained stable ($Δ{\rm DM} <0.15$ pc cm$^{-3}$) over $>10$ months. Overall, the observational characteristics of FRB 20200120E deviate quantitatively from other active repeaters, but it is unclear whether it is qualitatively a different type of source. △ Less

Submitted 23 January, 2023; v1 submitted 8 June, 2022; originally announced June 2022.

Comments: Accepted for publication in MNRAS

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

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

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

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

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

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

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

Two New Black Widow Millisecond Pulsars In M28

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

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

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

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

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

Modelling annual scintillation arc variations in PSR J1643-1224 using the Large European Array for Pulsars

Authors: G. Mall, R. A. Main, J. Antoniadis, C. G. Bassa, M. Burgay, S. Chen, I. Cognard, R. Concu, A. Corongiu, M. Gaikwad, H. Hu, G. H. Janssen, R. Karuppusamy, K. J. Lee, K. Liu, J. W. McKee, A. Melis, M. B. Mickaliger, D. Perrodin, M. Pilia, A. Possenti, D. J. Reardon, S. A. Sanidas, T. Sprenger, B. W. Stappers , et al. (3 additional authors not shown)

Abstract: In this work we study variations in the parabolic scintillation arcs of the binary millisecond pulsar PSR J1643-1224 over five years using the Large European Array for Pulsars (LEAP). The 2D power spectrum of scintillation, called the secondary spectrum, often shows a parabolic distribution of power, where the arc curvature encodes the relative velocities and distances of the pulsar, ionised inter… ▽ More In this work we study variations in the parabolic scintillation arcs of the binary millisecond pulsar PSR J1643-1224 over five years using the Large European Array for Pulsars (LEAP). The 2D power spectrum of scintillation, called the secondary spectrum, often shows a parabolic distribution of power, where the arc curvature encodes the relative velocities and distances of the pulsar, ionised interstellar medium (IISM), and Earth. We observe a clear parabolic scintillation arc which varies in curvature throughout the year. The distribution of power in the secondary spectra are inconsistent with a single scattering screen which is fully 1D, or entirely isotropic. We fit the observed arc curvature variations with two models; an isotropic scattering screen, and a model with two independent 1D screens. We measure the distance to the scattering screen to be in the range 114-223 pc, depending on the model, consistent with the known distance of the foreground large-diameter HII region Sh 2-27 (112+/-17 pc), suggesting that it is the dominant source of scattering. We obtain only weak constraints on the pulsar's orbital inclination and angle of periastron, since the scintillation pattern is not very sensitive to the pulsar's motion, since the screen is much closer to the Earth than the pulsar. More measurements of this kind - where scattering screens can be associated with foreground objects - will help to inform the origins and distribution of scattering screens within our galaxy. △ Less

Submitted 11 January, 2022; originally announced January 2022.

Comments: 11 pages, 6 figures. Accepted by MNRAS

The International Pulsar Timing Array second data release: Search for an isotropic Gravitational Wave Background

Authors: J. Antoniadis, Z. Arzoumanian, S. Babak, M. Bailes, A. -S. Bak Nielsen, P. T. Baker, C. G. Bassa, B. Becsy, A. Berthereau, M. Bonetti, A. Brazier, P. R. Brook, M. Burgay, S. Burke-Spolaor, R. N. Caballero, J. A. Casey-Clyde, A. Chalumeau, D. J. Champion, M. Charisi, S. Chatterjee, S. Chen, I. Cognard, J. M. Cordes, N. J. Cornish, F. Crawford , et al. (101 additional authors not shown)

Abstract: We searched for an isotropic stochastic gravitational wave background in the second data release of the International Pulsar Timing Array, a global collaboration synthesizing decadal-length pulsar-timing campaigns in North America, Europe, and Australia. In our reference search for a power law strain spectrum of the form $h_c = A(f/1\,\mathrm{yr}^{-1})^α$, we found strong evidence for a spectrally… ▽ More We searched for an isotropic stochastic gravitational wave background in the second data release of the International Pulsar Timing Array, a global collaboration synthesizing decadal-length pulsar-timing campaigns in North America, Europe, and Australia. In our reference search for a power law strain spectrum of the form $h_c = A(f/1\,\mathrm{yr}^{-1})^α$, we found strong evidence for a spectrally-similar low-frequency stochastic process of amplitude $A = 3.8^{+6.3}_{-2.5}\times10^{-15}$ and spectral index $α= -0.5 \pm 0.5$, where the uncertainties represent 95\% credible regions, using information from the auto- and cross-correlation terms between the pulsars in the array. For a spectral index of $α= -2/3$, as expected from a population of inspiralling supermassive black hole binaries, the recovered amplitude is $A = 2.8^{+1.2}_{-0.8}\times10^{-15}$. Nonetheless, no significant evidence of the Hellings-Downs correlations that would indicate a gravitational-wave origin was found. We also analyzed the constituent data from the individual pulsar timing arrays in a consistent way, and clearly demonstrate that the combined international data set is more sensitive. Furthermore, we demonstrate that this combined data set produces comparable constraints to recent single-array data sets which have more data than the constituent parts of the combination. Future international data releases will deliver increased sensitivity to gravitational wave radiation, and significantly increase the detection probability. △ Less

Submitted 11 January, 2022; originally announced January 2022.

Comments: 17 pages, 12 figures, accepted in MNRAS

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)

Removal and replacement of interference in tied-array radio pulsar observations using the spectral kurtosis estimator

Authors: M. Purver, C. G. Bassa, I. Cognard, G. H. Janssen, R. Karuppusamy, M. Kramer, K. J. Lee, K. Liu, J. W. McKee, D. Perrodin, S. Sanidas, R. Smits, B. W. Stappers

Abstract: We describe how to implement the spectral kurtosis method of interference removal (zapping) on a digitized signal of averaged power values. Spectral kurtosis is a hypothesis test, analogous to the t-test, with a null hypothesis that the amplitudes from which power is formed belong to a `good' distribution -- typically Gaussian with zero mean -- where power values are zapped if the hypothesis is re… ▽ More We describe how to implement the spectral kurtosis method of interference removal (zapping) on a digitized signal of averaged power values. Spectral kurtosis is a hypothesis test, analogous to the t-test, with a null hypothesis that the amplitudes from which power is formed belong to a `good' distribution -- typically Gaussian with zero mean -- where power values are zapped if the hypothesis is rejected at a specified confidence level. We derive signal-to-noise ratios (SNRs) as a function of amount of zapping for folded radio pulsar observations consisting of a sum of signals from multiple telescopes in independent radio-frequency interference (RFI) environments, comparing four methods to compensate for lost data with coherent (tied-array) and incoherent summation. For coherently summed amplitudes, scaling amplitudes from non-zapped telescopes achieves a higher SNR than replacing zapped amplitudes with artificial noise. For incoherently summed power values, the highest SNR is given by scaling power from non-zapped telescopes to maintain a constant mean. We use spectral kurtosis to clean a tied-array radio pulsar observation by the Large European Array for Pulsars (LEAP): the signal from one telescope is zapped with time and frequency resolutions of 6.25 $μ$s and 0.16 MHz, removing interference along with 0.27 per cent of `good' data, giving an uncertainty of 0.25 $μ$s in pulse time of arrival (TOA) for PSR J1022+1001. We use a single-telescope observation to demonstrate recovery of the pulse profile shape, with 0.6 per cent of data zapped and a reduction from 1.22 to 0.70 $μ$s in TOA uncertainty. △ Less

Submitted 24 February, 2023; v1 submitted 28 November, 2021; originally announced November 2021.

Comments: 16 pages, 9 figures

Journal ref: MNRAS, 510, 1597 (2022)

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

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

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

Submitted 26 November, 2021; originally announced November 2021.

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

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

Milliarcsecond localisation of the repeating FRB 20201124A

Authors: K. Nimmo, D. M. Hewitt, J. W. T. Hessels, F. Kirsten, B. Marcote, U. Bach, R. Blaauw, M. Burgay, A. Corongiu, R. Feiler, M. P. Gawroński, M. Giroletti, R. Karuppusamy, A. Keimpema, M. A. Kharinov, M. Lindqvist, G. Maccaferri, A. Melnikov, A. Mikhailov, O. S. Ould-Boukattine, Z. Paragi, M. Pilia, A. Possenti, M. P. Snelders, G. Surcis , et al. (6 additional authors not shown)

Abstract: Very long baseline interferometric (VLBI) localisations of repeating fast radio bursts (FRBs) have demonstrated a diversity of local environments: from nearby star-forming regions to globular clusters. Here we report the VLBI localisation of FRB 20201124A using an ad-hoc array of dishes that also participate in the European VLBI Network (EVN). In our campaign, we detected 18 total bursts from FRB… ▽ More Very long baseline interferometric (VLBI) localisations of repeating fast radio bursts (FRBs) have demonstrated a diversity of local environments: from nearby star-forming regions to globular clusters. Here we report the VLBI localisation of FRB 20201124A using an ad-hoc array of dishes that also participate in the European VLBI Network (EVN). In our campaign, we detected 18 total bursts from FRB 20201124A at two separate epochs. By combining the visibilities from both observing epochs, we were able to localise FRB 20201124A with a 1-$σ$ error of 4.5 milliarcseconds (mas). We use the relatively large burst sample to investigate astrometric accuracy, and find that for $\gtrsim20$ baselines ($\gtrsim7$ dishes) that we can robustly reach milliarcsecond precision even using single-burst data sets. Sub-arcsecond precision is still possible for single bursts, even when only $\sim$ six baselines (four dishes) are available. We explore two methods for determining the individual burst positions: the peaks of the dirty maps and a Gaussian fit to the cross fringe pattern on the dirty maps. We found the latter to be more reliable due to the lower mean and standard deviation in the offsets from the FRB position. Our VLBI work places FRB 20201124A 705$\pm$26 mas (1-$σ$ errors) from the optical centre of the host galaxy, and consistent with originating from within the recently-discovered extended radio structure associated with star-formation in the host galaxy. Future high-resolution optical observations, e.g. with Hubble Space Telescope, can determine the proximity of our FRB 20201124A VLBI position to nearby knots of star formation. △ Less

Submitted 2 November, 2021; originally announced November 2021.

Comments: submitted, comments welcome

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

Wide field beamformed observation with MeerKAT

Authors: Weiwei Chen, Ewan Barr, Ramesh Karuppusamy, Michael Kramer, Benjamin Stappers

Abstract: Large-scale beamforming with radio interferometers has the potential to revolutionize the science done with pulsars and fast radio bursts by improving the survey efficiency for these sources. We describe a wide-field beamformer for the MeerKAT radio telescope and outline strategies to optimally design such surveys. A software implementation of these techniques, ${\rm M{\small OSAIC}}$ is introduce… ▽ More Large-scale beamforming with radio interferometers has the potential to revolutionize the science done with pulsars and fast radio bursts by improving the survey efficiency for these sources. We describe a wide-field beamformer for the MeerKAT radio telescope and outline strategies to optimally design such surveys. A software implementation of these techniques, ${\rm M{\small OSAIC}}$ is introduced and its application in the MeerKAT telescope is presented. We show initial results using the beamformer by observing a globular cluster to track several pulsars simultaneously and demonstrate the source localization capability of this observation. △ Less

Submitted 7 November, 2022; v1 submitted 4 October, 2021; originally announced October 2021.

Journal ref: Journal of Astronomical Instrumentation, 10, 2150013 (2021)

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

Burst timescales and luminosities link young pulsars and fast radio bursts

Authors: K. Nimmo, J. W. T. Hessels, F. Kirsten, A. Keimpema, J. M. Cordes, M. P. Snelders, D. M. Hewitt, R. Karuppusamy, A. M. Archibald, V. Bezukovs, M. Bhardwaj, R. Blaauw, S. T. Buttaccio, T. Cassanelli, J. E. Conway, A. Corongiu, R. Feiler, E. Fonseca, O. Forssen, M. Gawronski, M. Giroletti, M. A. Kharinov, C. Leung, M. Lindqvist, G. Maccaferri , et al. (23 additional authors not shown)

Abstract: Fast radio bursts (FRBs) are extragalactic radio flashes of unknown physical origin. Their high luminosities and short durations require extreme energy densities, like those found in the vicinity of neutron stars and black holes. Studying the burst intensities and polarimetric properties on a wide range of timescales, from milliseconds down to nanoseconds, is key to understanding the emission mech… ▽ More Fast radio bursts (FRBs) are extragalactic radio flashes of unknown physical origin. Their high luminosities and short durations require extreme energy densities, like those found in the vicinity of neutron stars and black holes. Studying the burst intensities and polarimetric properties on a wide range of timescales, from milliseconds down to nanoseconds, is key to understanding the emission mechanism. However, high-time-resolution studies of FRBs are limited by their unpredictable activity levels, available instrumentation and temporal broadening in the intervening ionised medium. Here we show that the repeating FRB 20200120E can produce isolated shots of emission as short as about 60 nanoseconds in duration, with brightness temperatures as high as $3\times 10^{41}$ K (excluding relativistic effects), comparable to `nano-shots' from the Crab pulsar. Comparing both the range of timescales and luminosities, we find that FRB 20200120E observationally bridges the gap between known Galactic young pulsars and magnetars, and the much more distant extragalactic FRBs. This suggests a common magnetically powered emission mechanism spanning many orders of magnitude in timescale and luminosity. In this work, we probe a relatively unexplored region of the short-duration transient phase space; we highlight that there likely exists a population of ultra-fast radio transients at nanosecond to microsecond timescales, which current FRB searches are insensitive to. △ Less

Submitted 29 September, 2021; v1 submitted 24 May, 2021; originally announced May 2021.

Comments: Submitted. Comments welcome

A repeating fast radio burst source in a globular cluster

Authors: F. Kirsten, B. Marcote, K. Nimmo, J. W. T. Hessels, M. Bhardwaj, S. P. Tendulkar, A. Keimpema, J. Yang, M. P. Snelders, P. Scholz, A. B. Pearlman, C. J. Law, W. M. Peters, M. Giroletti, Z. Paragi, C. Bassa, D. M. Hewitt, U. Bach, V. Bezrukovs, M. Burgay, S. T. Buttaccio, J. E. Conway, A. Corongiu, R. Feiler, O. Forssén , et al. (41 additional authors not shown)

Abstract: Fast radio bursts (FRBs) are exceptionally luminous flashes of unknown physical origin, reaching us from other galaxies (Petroff et al. 2019). Most FRBs have only ever been seen once, while others flash repeatedly, though sporadically (Spitler et al. 2016, CHIME/FRB Collaboration et al. 2021). Many models invoke magnetically powered neutron stars (magnetars) as the engines producing FRB emission (… ▽ More Fast radio bursts (FRBs) are exceptionally luminous flashes of unknown physical origin, reaching us from other galaxies (Petroff et al. 2019). Most FRBs have only ever been seen once, while others flash repeatedly, though sporadically (Spitler et al. 2016, CHIME/FRB Collaboration et al. 2021). Many models invoke magnetically powered neutron stars (magnetars) as the engines producing FRB emission (Margalit & Metzger 2018, CHIME/FRB Collaboration et al. 2020). Recently, CHIME/FRB announced the discovery (Bhardwaj et al. 2021) of the repeating FRB 20200120E, coming from the direction of the nearby grand design spiral galaxy M81. Four potential counterparts at other observing wavelengths were identified (Bhardwaj et al. 2021) but no definitive association with these sources, or M81, could be made. Here we report an extremely precise localisation of FRB 20200120E, which allows us to associate it with a globular cluster (GC) in the M81 galactic system and to place it ~2pc offset from the optical center of light of the GC. This confirms (Bhardwaj et al. 2021) that FRB 20200120E is 40 times closer than any other known extragalactic FRB. Because such GCs host old stellar populations, this association strongly challenges FRB models that invoke young magnetars formed in a core-collapse supernova as powering FRB emission. We propose, instead, that FRB 20200120E is a highly magnetised neutron star formed via either accretion-induced collapse of a white dwarf or via merger of compact stars in a binary system (Margalit et al. 2019). Alternative scenarios involving compact binary systems, efficiently formed inside globular clusters, could also be responsible for the observed bursts. △ Less

Submitted 29 September, 2021; v1 submitted 24 May, 2021; originally announced May 2021.

Comments: Submitted. Comments welcome

The Polarized Image of a Synchrotron Emitting Ring of Gas Orbiting a Black Hole

Authors: Ramesh Narayan, Daniel C. M. Palumbo, Michael D. Johnson, Zachary Gelles, Elizabeth Himwich, Dominic O. Chang, Angelo Ricarte, Jason Dexter, Charles F. Gammie, Andrew A. Chael, The Event Horizon Telescope Collaboration, :, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Anne-Kathrin Baczko, David Ball, Mislav Balokovic, John Barrett, Bradford A. Benson, Dan Bintley , et al. (215 additional authors not shown)

Abstract: Synchrotron radiation from hot gas near a black hole results in a polarized image. The image polarization is determined by effects including the orientation of the magnetic field in the emitting region, relativistic motion of the gas, strong gravitational lensing by the black hole, and parallel transport in the curved spacetime. We explore these effects using a simple model of an axisymmetric, equ… ▽ More Synchrotron radiation from hot gas near a black hole results in a polarized image. The image polarization is determined by effects including the orientation of the magnetic field in the emitting region, relativistic motion of the gas, strong gravitational lensing by the black hole, and parallel transport in the curved spacetime. We explore these effects using a simple model of an axisymmetric, equatorial accretion disk around a Schwarzschild black hole. By using an approximate expression for the null geodesics derived by Beloborodov (2002) and conservation of the Walker-Penrose constant, we provide analytic estimates for the image polarization. We test this model using currently favored general relativistic magnetohydrodynamic simulations of M87*, using ring parameters given by the simulations. For a subset of these with modest Faraday effects, we show that the ring model broadly reproduces the polarimetric image morphology. Our model also predicts the polarization evolution for compact flaring regions, such as those observed from Sgr A* with GRAVITY. With suitably chosen parameters, our simple model can reproduce the EVPA pattern and relative polarized intensity in Event Horizon Telescope images of M87*. Under the physically motivated assumption that the magnetic field trails the fluid velocity, this comparison is consistent with the clockwise rotation inferred from total intensity images. △ Less

Submitted 13 May, 2021; v1 submitted 4 May, 2021; originally announced May 2021.

Comments: 29 pages, 11 figures, published in ApJ on May 3

Journal ref: ApJ 912 35 (2021)

An 86-GHz search for Pulsars in the Galactic Center with the Atacama Large Millimeter/submillimeter Array

Authors: Kuo Liu, Gregory Desvignes, Ralph P. Eatough, Ramesh Karuppusamy, Michael Kramer, Pablo Torne, Robert Wharton, Shami Chatterjee, James M. Cordes, Geoffrey B. Crew, Ciriaco Goddi, Scott M. Ransom, Helge Rottmann, Federico Abbate, Geoffrey C. Bower, Christiaan D. Brinkerink, Heino Falcke, Aristeidis Noutsos, Antonio Hernandez-Gomez, Wu Jiang, Michael D. Johnson, Ru-Sen Lu, Yurii Pidopryhora, Luciano Rezzolla, Lijing Shao , et al. (2 additional authors not shown)

Abstract: We report on the first pulsar and transient survey of the Galactic Center (GC) with the Atacama Large Millimeter/submillimeter Array (ALMA). The observations were conducted during the Global Millimeter VLBI Array campaign in 2017 and 2018. We carry out searches using timeseries of both total intensity and other polarization components in the form of Stokes parameters. We incorporate acceleration a… ▽ More We report on the first pulsar and transient survey of the Galactic Center (GC) with the Atacama Large Millimeter/submillimeter Array (ALMA). The observations were conducted during the Global Millimeter VLBI Array campaign in 2017 and 2018. We carry out searches using timeseries of both total intensity and other polarization components in the form of Stokes parameters. We incorporate acceleration and its derivative in the pulsar search, and also search in segments of the entire observation to compensate for potential orbital motion of the pulsar. While no new pulsar is found, our observations yield the polarization profile of the GC magnetar PSR J1745-2900 at mm-wavelength for the first time, which turns out to be nearly 100 % linearly polarized. Additionally, we estimate the survey sensitivity placed by both system and red noise, and evaluate its capability of finding pulsars in orbital motion with either Sgr A* or a binary companion. We show that the survey is sensitive to only the most luminous pulsars in the known population, and future observations with ALMA in Band-1 will deliver significantly deeper survey sensitivity on the GC pulsar population. △ Less

Submitted 18 April, 2021; originally announced April 2021.

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

Searching for pulsars in the Galactic Centre at 3 and 2 mm

Authors: Pablo Torne, Gregory Desvignes, Ralph Eatough, Michael Kramer, Ramesh Karuppusamy, Kuo Liu, Aris Noutsos, Robert Wharton, Carsten Kramer, Santiago Navarro, Gabriel Paubert, Salvador Sanchez, Miguel Sanchez-Portal, Karl Schuster, Heino Falcke, Luciano Rezzolla

Abstract: Pulsars in the Galactic centre promise to enable unparalleled tests of gravity theories and black hole physics and to serve as probes of the stellar formation history and evolution and the interstellar medium in the complex central region of the Milky Way. The community has surveyed the innermost region of the galaxy for decades without detecting a population of pulsars, which is puzzling. A stron… ▽ More Pulsars in the Galactic centre promise to enable unparalleled tests of gravity theories and black hole physics and to serve as probes of the stellar formation history and evolution and the interstellar medium in the complex central region of the Milky Way. The community has surveyed the innermost region of the galaxy for decades without detecting a population of pulsars, which is puzzling. A strong scattering of the pulsed signals in this particular direction has been argued to be a potential reason for the non-detections. Scattering has a strong inverse dependence on observing frequency, therefore an effective way to alleviate its effect is to use higher frequencies in a survey for pulsars in the Galactic centre, in particular, close to the supermassive black hole Sagittarius A*. We present the first pulsar survey at short millimetre wavelengths, using several frequency bands between 84 and 156 GHz (3.57-1.92 mm), targeted to the Galactic centre. The observations were made with the Institut de Radioastronomie Millimetrique (IRAM) 30m Telescope in 28 epochs between 2016 December and 2018 May. This survey is the first that is essentially unaffected by scattering and therefore unbiased in population coverage, including fast-spinning pulsars that might be out of reach of lower-frequency Galactic centre surveys. We discovered no new pulsars and relate this result mainly to the decreased flux density of pulsars at high frequencies, combined with our current sensitivity. However, we demonstrate that surveys at these extremely high radio frequencies are capable of discovering new pulsars, analyse their sensitivity limits with respect to a simulated Galactic centre pulsar population, and discuss the main challenges and possible improvements for similar surveys in the future. △ Less

Submitted 14 April, 2021; v1 submitted 30 March, 2021; originally announced March 2021.

Comments: 14 pages, 6 figures, 4 tables. Accepted for publication in A&A. This version includes language editor corrections

Journal ref: A&A 650, A95 (2021)

Multi-frequency observations of SGR J1935+2154

Authors: M. Bailes, C. G. Bassa, G. Bernardi, S. Buchner, M. Burgay, M. Caleb, A. J. Cooper, G. Desvignes, P. J. Groot, I. Heywood, F. Jankowski, R. Karuppusamy, M. Kramer, M. Malenta, G. Naldi, M. Pilia, G. Pupillo, K. M. Rajwade, L. Spitler, M. Surnis, B. W. Stappers, A. Addis, S. Bloemen, M. C. Bezuidenhout, G. Bianchi , et al. (32 additional authors not shown)

Abstract: Magnetars are a promising candidate for the origin of Fast Radio Bursts (FRBs). The detection of an extremely luminous radio burst from the Galactic magnetar SGR J1935+2154 on 2020 April 28 added credence to this hypothesis. We report on simultaneous and non-simultaneous observing campaigns using the Arecibo, Effelsberg, LOFAR, MeerKAT, MK2 and Northern Cross radio telescopes and the MeerLICHT opt… ▽ More Magnetars are a promising candidate for the origin of Fast Radio Bursts (FRBs). The detection of an extremely luminous radio burst from the Galactic magnetar SGR J1935+2154 on 2020 April 28 added credence to this hypothesis. We report on simultaneous and non-simultaneous observing campaigns using the Arecibo, Effelsberg, LOFAR, MeerKAT, MK2 and Northern Cross radio telescopes and the MeerLICHT optical telescope in the days and months after the April 28 event. We did not detect any significant single radio pulses down to fluence limits between 25 mJy ms and 18 Jy ms. Some observing epochs overlapped with times when X-ray bursts were detected. Radio images made on four days using the MeerKAT telescope revealed no point-like persistent or transient emission at the location of the magnetar. No transient or persistent optical emission was detected over seven days. Using the multi-colour MeerLICHT images combined with relations between DM, NH and reddening we constrain the distance to SGR J1935+2154, to be between 1.5 and 6.5 kpc. The upper limit is consistent with some other distance indicators and suggests that the April 28 burst is closer to two orders of magnitude less energetic than the least energetic FRBs. The lack of single-pulse radio detections shows that the single pulses detected over a range of fluences are either rare, or highly clustered, or both. It may also indicate that the magnetar lies somewhere between being radio-quiet and radio-loud in terms of its ability to produce radio emission efficiently. △ Less

Submitted 10 March, 2021; originally announced March 2021.

Comments: Accepted for publication in MNRAS. Corresponding author B. W. Stappers

Eight new millisecond pulsars from the first MeerKAT globular cluster census

Authors: A. Ridolfi, T. Gautam, P. C. C. Freire, S. M. Ransom, S. J. Buchner, A. Possenti, V. Venkatraman Krishnan, M. Bailes, M. Kramer, B. W. Stappers, F. Abbate, E. D. Barr, M. Burgay, F. Camilo, A. Corongiu, A. Jameson, P. V. Padmanabh, L. Vleeschower, D. J. Champion, M. Geyer, A. Karastergiou, R. Karuppusamy, A. Parthasarathy, D. J. Reardon, M. Serylak , et al. (2 additional authors not shown)

Abstract: We have used the central 44 antennas of the new 64-dish MeerKAT radio telescope array to conduct a deep search for new pulsars in the core of nine globular clusters. This has led to the discovery of eight new millisecond pulsars in six different clusters. Two new binaries, 47 Tuc ac and 47 Tuc ad, are eclipsing "spiders", featuring compact orbits ($\lesssim 0.32$ days), very low-mass companions an… ▽ More We have used the central 44 antennas of the new 64-dish MeerKAT radio telescope array to conduct a deep search for new pulsars in the core of nine globular clusters. This has led to the discovery of eight new millisecond pulsars in six different clusters. Two new binaries, 47 Tuc ac and 47 Tuc ad, are eclipsing "spiders", featuring compact orbits ($\lesssim 0.32$ days), very low-mass companions and regular occultations of their pulsed emission. The other three new binary pulsars (NGC 6624G, M62G, and Ter 5 an) are in wider ($> 0.7$ days) orbits, with companions that are likely to be white dwarfs or neutron stars. NGC 6624G has a large eccentricity of $e\simeq 0.38$, which enabled us to detect the rate of advance of periastron. This suggests that the system is massive, with a total mass of $M{\rm tot} = 2.65 \pm 0.07$ M$_{\odot}$. Likewise, for Ter 5 an, with $e \simeq 0.0066$, we obtain $M{\rm tot}= 2.97 \pm 0.52$ M$_{\odot}$. The other three new discoveries (NGC 6522D, NGC 6624H and NGC 6752F) are faint isolated pulsars. Finally, we have used the whole MeerKAT array and synthesized 288 beams, covering an area of $\sim2$ arcmin in radius around the center of NGC 6624. This has allowed us to localize many of the pulsars in the cluster, demonstrating the beamforming capabilities of the TRAPUM software backend and paving the way for the upcoming MeerKAT globular cluster pulsar survey. △ Less

Submitted 8 March, 2021; originally announced March 2021.

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

LOFAR Detection of 110-188 MHz Emission and Frequency-Dependent Activity from FRB 20180916B

Authors: Z. Pleunis, D. Michilli, C. G. Bassa, J. W. T. Hessels, A. Naidu, B. C. Andersen, P. Chawla, E. Fonseca, A. Gopinath, V. M. Kaspi, V. I. Kondratiev, D. Z. Li, M. Bhardwaj, P. J. Boyle, C. Brar, T. Cassanelli, Y. Gupta, A. Josephy, R. Karuppusamy, A. Keimpema, F. Kirsten, C. Leung, B. Marcote, K. Masui, R. Mckinven , et al. (10 additional authors not shown)

Abstract: FRB 20180916B is a well-studied repeating fast radio burst source. Its proximity (~150 Mpc), along with detailed studies of the bursts, have revealed many clues about its nature -- including a 16.3-day periodicity in its activity. Here we report on the detection of 18 bursts using LOFAR at 110-188 MHz, by far the lowest-frequency detections of any FRB to date. Some bursts are seen down to the lowe… ▽ More FRB 20180916B is a well-studied repeating fast radio burst source. Its proximity (~150 Mpc), along with detailed studies of the bursts, have revealed many clues about its nature -- including a 16.3-day periodicity in its activity. Here we report on the detection of 18 bursts using LOFAR at 110-188 MHz, by far the lowest-frequency detections of any FRB to date. Some bursts are seen down to the lowest-observed frequency of 110 MHz, suggesting that their spectra extend even lower. These observations provide an order-of-magnitude stronger constraint on the optical depth due to free-free absorption in the source's local environment. The absence of circular polarization and nearly flat polarization angle curves are consistent with burst properties seen at 300-1700 MHz. Compared with higher frequencies, the larger burst widths (~40-160 ms at 150 MHz) and lower linear polarization fractions are likely due to scattering. We find ~2-3 rad/m^2 variations in the Faraday rotation measure that may be correlated with the activity cycle of the source. We compare the LOFAR burst arrival times to those of 38 previously published and 22 newly detected bursts from the uGMRT (200-450 MHz) and CHIME/FRB (400-800 MHz). Simultaneous observations show 5 CHIME/FRB bursts when no emission is detected by LOFAR. We find that the burst activity is systematically delayed towards lower frequencies by ~3 days from 600 MHz to 150 MHz. We discuss these results in the context of a model in which FRB 20180916B is an interacting binary system featuring a neutron star and high-mass stellar companion. △ Less

Submitted 4 March, 2021; v1 submitted 15 December, 2020; originally announced December 2020.

Comments: Accepted for publication by ApJL

Highly polarised microstructure from the repeating FRB 20180916B

Authors: K. Nimmo, J. W. T. Hessels, A. Keimpema, A. M. Archibald, J. M. Cordes, R. Karuppusamy, F. Kirsten, D. Z. Li, B. Marcote, Z. Paragi

Abstract: Fast radio bursts (FRBs) are bright, coherent, short-duration radio transients of as-yet unknown extragalactic origin. FRBs exhibit a wide variety of spectral, temporal and polarimetric properties, which can unveil clues into their emission physics and propagation effects in the local medium. Here we present the high-time-resolution (down to 1 $μ$s) polarimetric properties of four 1.7-GHz bursts f… ▽ More Fast radio bursts (FRBs) are bright, coherent, short-duration radio transients of as-yet unknown extragalactic origin. FRBs exhibit a wide variety of spectral, temporal and polarimetric properties, which can unveil clues into their emission physics and propagation effects in the local medium. Here we present the high-time-resolution (down to 1 $μ$s) polarimetric properties of four 1.7-GHz bursts from the repeating FRB 20180916B, which were detected in voltage data during observations with the European VLBI Network (EVN). We observe a range of emission timescales spanning three orders of magnitude, with the shortest component width reaching 3-4 $μ$s (below which we are limited by scattering). This is the shortest timescale measured in any FRB, to date. We demonstrate that all four bursts are highly linearly polarised ($\gtrsim 80\%$), show no evidence for significant circular polarisation ($\lesssim 15\%$), and exhibit a constant polarisation position angle (PPA) during and between bursts. On short timescales ($\lesssim 100$ $μ$s), however, there appear to be subtle (few degree) PPA variations across the burst profiles. These observational results are most naturally explained in an FRB model where the emission is magnetospheric in origin, as opposed to models where the emission originates at larger distances in a relativistic shock. △ Less

Submitted 5 January, 2021; v1 submitted 12 October, 2020; originally announced October 2020.

Comments: Submitted; comments welcome

Revisiting profile instability of J1022+1001

Authors: Prajwal V. Padmanabh, Ewan D. Barr, David J. Champion, Ramesh Karuppusamy, Michael Kramer, Axel Jessner, Patrick Lazarus

Abstract: Millisecond pulsars in timing arrays can act as probes for gravitational wave detection and improving the solar system ephemerides among several other applications. However, the stability of the integrated pulse profiles can limit the precision of the ephemeris parameters and in turn the applications derived from it. It is thus crucial for the pulsars in the array to have stable integrated pulse p… ▽ More Millisecond pulsars in timing arrays can act as probes for gravitational wave detection and improving the solar system ephemerides among several other applications. However, the stability of the integrated pulse profiles can limit the precision of the ephemeris parameters and in turn the applications derived from it. It is thus crucial for the pulsars in the array to have stable integrated pulse profiles. Here we present evidence for long-term profile instability in PSR J1022+1001 which is currently included in the European and Parkes pulsar timing arrays. We apply a new evaluation method to an expanded data set ranging from the Effelsberg Pulsar Observing System back-end used in the 1990s to that of data from the current PSRIX backend at the Effelsberg Radio Telescope. We show that this intrinsic variability in the pulse shape persists over time scales of years. We investigate if systematic instrumental effects like polarisation calibration or signal propagation effects in the interstellar medium causes the observed profile instability. We find that the total variation cannot be fully accounted for by instrumental and propagation effects. This suggests additional intrinsic effects as the origin for the variation. We finally discuss several factors that could lead to the observed behaviour and comment on the consequent implications. △ Less

Submitted 8 October, 2020; originally announced October 2020.

Comments: 10 pages, 8 figures, Accepted for publication in MNRAS

Gravitational Test Beyond the First Post-Newtonian Order with the Shadow of the M87 Black Hole

Authors: Dimitrios Psaltis, Lia Medeiros, Pierre Christian, Feryal Ozel, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Keiichi Asada, Rebecca Azulay, David Ball, Mislav Balokovic, John Barrett, Dan Bintley, Lindy Blackburn, Wilfred Boland, Geoffrey C. Bower, Michael Bremer, Christiaan D. Brinkerink, Roger Brissenden, Silke Britzen, Dominique Broguiere, Thomas Bronzwaer, Do-Young Byun, John E. Carlstrom, Andrew Chael , et al. (163 additional authors not shown)

Abstract: The 2017 Event Horizon Telescope (EHT) observations of the central source in M87 have led to the first measurement of the size of a black-hole shadow. This observation offers a new and clean gravitational test of the black-hole metric in the strong-field regime. We show analytically that spacetimes that deviate from the Kerr metric but satisfy weak-field tests can lead to large deviations in the p… ▽ More The 2017 Event Horizon Telescope (EHT) observations of the central source in M87 have led to the first measurement of the size of a black-hole shadow. This observation offers a new and clean gravitational test of the black-hole metric in the strong-field regime. We show analytically that spacetimes that deviate from the Kerr metric but satisfy weak-field tests can lead to large deviations in the predicted black-hole shadows that are inconsistent with even the current EHT measurements. We use numerical calculations of regular, parametric, non-Kerr metrics to identify the common characteristic among these different parametrizations that control the predicted shadow size. We show that the shadow-size measurements place significant constraints on deviation parameters that control the second post-Newtonian and higher orders of each metric and are, therefore, inaccessible to weak-field tests. The new constraints are complementary to those imposed by observations of gravitational waves from stellar-mass sources. △ Less

Submitted 2 October, 2020; originally announced October 2020.

Comments: Physical Review Letters

Observing superluminous supernovae and long gamma ray bursts as potential birthplaces of repeating fast radio bursts

Authors: G. H. Hilmarsson, L. G. Spitler, E. F. Keane, T. M. Athanasiadis, E. Barr, M. Cruces, X. Deng, S. Heyminck, R. Karuppusamy, M. Kramer, S. P. Sathyanarayanan, V. Ventakraman Krishnan, G. Wieching, J. Wu, O. Wucknitz

Abstract: Superluminous supernovae (SLSNe) and long gamma ray bursts (LGRBs) have been proposed as progenitors of repeating Fast Radio Bursts (FRBs). In this scenario, bursts originate from the interaction between a young magnetar and its surrounding supernova remnant (SNR). Such a model could explain the repeating, apparently non-Poissonian nature of FRB121102, which appears to display quiescent and active… ▽ More Superluminous supernovae (SLSNe) and long gamma ray bursts (LGRBs) have been proposed as progenitors of repeating Fast Radio Bursts (FRBs). In this scenario, bursts originate from the interaction between a young magnetar and its surrounding supernova remnant (SNR). Such a model could explain the repeating, apparently non-Poissonian nature of FRB121102, which appears to display quiescent and active phases. This bursting behaviour is better explained with a Weibull distribution, which includes parametrisation for clustering. We observed 10 SLSNe/LGRBs for 63 hours, looking for repeating FRBs with the Effelsberg-100 m radio telescope, but have not detected any bursts. We scale the burst rate of FRB121102 to an FRB121102-like source inhabiting each of our observed targets, and compare this rate to our upper burst rate limit on a source by source basis. By adopting a fiducial beaming fraction of 0.6, we obtain 99.99\% and 83.4\% probabilities that at least one, and at least half of our observed sources are beamed towards us respectively. One of our SLSN targets, PTF10hgi, is coincident with a persistent radio source, making it a possible analogue to FRB121102. We performed further observations on this source using the Effelsberg-100~m and Parkes-64~m radio telescopes. Assuming that PTF10hgi contains an FRB121102-like source, the probabilities of not detecting any bursts from a Weibull distribution during our observations are 14\% and 16\% for Effelsberg and Parkes respectively. We conclude by showing that a survey of many short observations increases burst detection probability for a source with Weibull distributed bursting activity. △ Less

Submitted 29 September, 2020; originally announced September 2020.

Comments: 11 pages, 5 figures

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 493, Issue 4, p.5170-5180, 2020

Measuring Interstellar Delays of PSR J0613-0200 over 7 years, using the Large European Array for Pulsars

Authors: R. A. Main, S. A. Sanidas, J. Antoniadis, C. Bassa, S. Chen, I. Cognard, M. Gaikwad, H. Hu, G. H. Janssen, R. Karuppusamy, M. Kramer, K. J. Lee, K. Liu, G. Mall, J. W. McKee, M. B. Mickaliger, D. Perrodin, B. W. Stappers, C. Tiburzi, O. Wucknitz, L. Wang, W. W. Zhu

Abstract: Using data from the Large European Array for Pulsars (LEAP), and the Effelsberg telescope, we study the scintillation parameters of the millisecond pulsar J0613-0200 over a 7 year timespan. The "secondary spectrum" -- the 2D power spectrum of scintillation -- presents the scattered power as a function of time delay, and contains the relative velocities of the pulsar, observer, and scattering mater… ▽ More Using data from the Large European Array for Pulsars (LEAP), and the Effelsberg telescope, we study the scintillation parameters of the millisecond pulsar J0613-0200 over a 7 year timespan. The "secondary spectrum" -- the 2D power spectrum of scintillation -- presents the scattered power as a function of time delay, and contains the relative velocities of the pulsar, observer, and scattering material. We detect a persistent parabolic scintillation arc, suggesting scattering is dominated by a thin, anisotropic region. The scattering is poorly described by a simple exponential tail, with excess power at high delays; we measure significant, detectable scattered power at times out to $\sim 5 μs$, and measure the bulk scattering delay to be between 50 to 200\,ns with particularly strong scattering throughout 2013. These delays are too small to detect a change of the pulse profile shape, yet they would change the times-of-arrival as measured through pulsar timing. The arc curvature varies annually, and is well fit by a one-dimensional scattering screen $\sim 40\%$ of the way towards the pulsar, with a changing orientation during the increased scattering in 2013. Effects of uncorrected scattering will introduce time delays correlated over time in individual pulsars, and may need to be considered in gravitational wave analyses. Pulsar timing programs would benefit from simultaneously recording in a way that scintillation can be resolved, in order to monitor the variable time delays caused by multipath propagation. △ Less

Submitted 3 November, 2020; v1 submitted 22 September, 2020; originally announced September 2020.

Comments: 12 pages, 10 figures. Published in MNRAS