Showing 1–50 of 55 results for author: Perera, B B P

The Anomalous Acceleration of PSR J2043+1711: Long-Period Orbital Companion or Stellar Flyby?

Authors: Thomas Donlon II, Sukanya Chakrabarti, Michael T. Lam, Daniel Huber, Daniel Hey, Enrico Ramirez-Ruiz, Benjamin Shappee, David L. Kaplan, Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Paul R. Brook, H. Thankful Cromartie, Kathryn Crowter, Megan E. DeCesar, Paul B. Demorest, Timothy Dolch, Elizabeth C. Ferrara, William Fiore, Emmanuel Fonseca, Gabriel E. Freedman, Nate Garver-Daniels, Peter A. Gentile , et al. (31 additional authors not shown)

Abstract: Based on the rate of change of its orbital period, PSR J2043+1711 has a substantial peculiar acceleration of 3.5 $\pm$ 0.8 mm/s/yr, which deviates from the acceleration predicted by equilibrium Milky Way models at a $4σ$ level. The magnitude of the peculiar acceleration is too large to be explained by disequilibrium effects of the Milky Way interacting with orbiting dwarf galaxies ($\sim$1 mm/s/yr… ▽ More Based on the rate of change of its orbital period, PSR J2043+1711 has a substantial peculiar acceleration of 3.5 $\pm$ 0.8 mm/s/yr, which deviates from the acceleration predicted by equilibrium Milky Way models at a $4σ$ level. The magnitude of the peculiar acceleration is too large to be explained by disequilibrium effects of the Milky Way interacting with orbiting dwarf galaxies ($\sim$1 mm/s/yr), and too small to be caused by period variations due to the pulsar being a redback. We identify and examine two plausible causes for the anomalous acceleration: a stellar flyby, and a long-period orbital companion. We identify a main-sequence star in \textit{Gaia} DR3 and Pan-STARRS DR2 with the correct mass, distance, and on-sky position to potentially explain the observed peculiar acceleration. However, the star and the pulsar system have substantially different proper motions, indicating that they are not gravitationally bound. However, it is possible that this is an unrelated star that just happens to be located near J2043+1711 along our line of sight (chance probability of 1.6\%). Therefore, we also constrain possible orbital parameters for a circumbinary companion in a hierarchical triple system with J2043+1711; the changes in the spindown rate of the pulsar are consistent with an outer object that has an orbital period of 80 kyr, a companion mass of 0.3 $M_\odot$ (indicative of a white dwarf or low-mass star), and a semi-major axis of 2000 AU. Continued timing and/or future faint optical observations of J2043+1711 may eventually allow us to differentiate between these scenarios. △ Less

Submitted 8 July, 2024; originally announced July 2024.

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

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

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

Submitted 23 May, 2024; originally announced May 2024.

The NANOGrav 15-year data set: Search for Transverse Polarization Modes in the Gravitational-Wave Background

Authors: Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Jeremy Baier, Paul T. Baker, Bence Bécsy, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Rand Burnette, Robin Case, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Megan E. DeCesar, Dallas DeGan, Paul B. Demorest , et al. (74 additional authors not shown)

Abstract: Recently we found compelling evidence for a gravitational wave background with Hellings and Downs (HD) correlations in our 15-year data set. These correlations describe gravitational waves as predicted by general relativity, which has two transverse polarization modes. However, more general metric theories of gravity can have additional polarization modes which produce different interpulsar correl… ▽ More Recently we found compelling evidence for a gravitational wave background with Hellings and Downs (HD) correlations in our 15-year data set. These correlations describe gravitational waves as predicted by general relativity, which has two transverse polarization modes. However, more general metric theories of gravity can have additional polarization modes which produce different interpulsar correlations. In this work we search the NANOGrav 15-year data set for evidence of a gravitational wave background with quadrupolar Hellings and Downs (HD) and Scalar Transverse (ST) correlations. We find that HD correlations are the best fit to the data, and no significant evidence in favor of ST correlations. While Bayes factors show strong evidence for a correlated signal, the data does not strongly prefer either correlation signature, with Bayes factors $\sim 2$ when comparing HD to ST correlations, and $\sim 1$ for HD plus ST correlations to HD correlations alone. However, when modeled alongside HD correlations, the amplitude and spectral index posteriors for ST correlations are uninformative, with the HD process accounting for the vast majority of the total signal. Using the optimal statistic, a frequentist technique that focuses on the pulsar-pair cross-correlations, we find median signal-to-noise-ratios of 5.0 for HD and 4.6 for ST correlations when fit for separately, and median signal-to-noise-ratios of 3.5 for HD and 3.0 for ST correlations when fit for simultaneously. While the signal-to-noise-ratios for each of the correlations are comparable, the estimated amplitude and spectral index for HD are a significantly better fit to the total signal, in agreement with our Bayesian analysis. △ Less

Submitted 18 October, 2023; originally announced October 2023.

Comments: 11 pages, 5 figures

How to Detect an Astrophysical Nanohertz Gravitational-Wave Background

Authors: Bence Bécsy, Neil J. Cornish, Patrick M. Meyers, Luke Zoltan Kelley, Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Katerina Chatziioannou, Tyler Cohen, James M. Cordes, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Megan E. DeCesar, Paul B. Demorest, Timothy Dolch , et al. (71 additional authors not shown)

Abstract: Analysis of pulsar timing data have provided evidence for a stochastic gravitational wave background in the nHz frequency band. The most plausible source of such a background is the superposition of signals from millions of supermassive black hole binaries. The standard statistical techniques used to search for such a background and assess its significance make several simplifying assumptions, nam… ▽ More Analysis of pulsar timing data have provided evidence for a stochastic gravitational wave background in the nHz frequency band. The most plausible source of such a background is the superposition of signals from millions of supermassive black hole binaries. The standard statistical techniques used to search for such a background and assess its significance make several simplifying assumptions, namely: i) Gaussianity; ii) isotropy; and most often iii) a power-law spectrum. However, a stochastic background from a finite collection of binaries does not exactly satisfy any of these assumptions. To understand the effect of these assumptions, we test standard analysis techniques on a large collection of realistic simulated datasets. The dataset length, observing schedule, and noise levels were chosen to emulate the NANOGrav 15-year dataset. Simulated signals from millions of binaries drawn from models based on the Illustris cosmological hydrodynamical simulation were added to the data. We find that the standard statistical methods perform remarkably well on these simulated datasets, despite their fundamental assumptions not being strictly met. They are able to achieve a confident detection of the background. However, even for a fixed set of astrophysical parameters, different realizations of the universe result in a large variance in the significance and recovered parameters of the background. We also find that the presence of loud individual binaries can bias the spectral recovery of the background if we do not account for them. △ Less

Submitted 1 December, 2023; v1 submitted 8 September, 2023; originally announced September 2023.

Comments: 14 pages, 8 figures, version matching published paper

Journal ref: ApJ 959 9 (2023)

Comparing recent PTA results on the nanohertz stochastic gravitational wave background

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

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

Submitted 1 September, 2023; originally announced September 2023.

Comments: 21 pages, 9 figures, submitted to ApJ

The 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 NANOGrav 15-year Gravitational-Wave Background Analysis Pipeline

Authors: Aaron D. Johnson, Patrick M. Meyers, Paul T. Baker, Neil J. Cornish, Jeffrey S. Hazboun, Tyson B. Littenberg, Joseph D. Romano, Stephen R. Taylor, Michele Vallisneri, Sarah J. Vigeland, Ken D. Olum, Xavier Siemens, Justin A. Ellis, Rutger van Haasteren, Sophie Hourihane, Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Bence Bécsy, J. Andrew Casey-Clyde , et al. (71 additional authors not shown)

Abstract: This paper presents rigorous tests of pulsar timing array methods and software, examining their consistency across a wide range of injected parameters and signal strength. We discuss updates to the 15-year isotropic gravitational-wave background analyses and their corresponding code representations. Descriptions of the internal structure of the flagship algorithms \texttt{Enterprise} and \texttt{P… ▽ More This paper presents rigorous tests of pulsar timing array methods and software, examining their consistency across a wide range of injected parameters and signal strength. We discuss updates to the 15-year isotropic gravitational-wave background analyses and their corresponding code representations. Descriptions of the internal structure of the flagship algorithms \texttt{Enterprise} and \texttt{PTMCMCSampler} are given to facilitate understanding of the PTA likelihood structure, how models are built, and what methods are currently used in sampling the high-dimensional PTA parameter space. We introduce a novel version of the PTA likelihood that uses a two-step marginalization procedure that performs much faster when the white noise parameters remain fixed. We perform stringent tests of consistency and correctness of the Bayesian and frequentist analysis software. For the Bayesian analysis, we test prior recovery, injection recovery, and Bayes factors. For the frequentist analysis, we test that the cross-correlation-based optimal statistic, when modified to account for a non-negligible gravitational-wave background, accurately recovers the amplitude of the background. We also summarize recent advances and tests performed on the optimal statistic in the literature from both GWB detection and parameter estimation perspectives. The tests presented here validate current and future analyses of PTA data. △ Less

Submitted 7 July, 2023; v1 submitted 28 June, 2023; originally announced June 2023.

Comments: 30 pages, 10 figures, 1 table; Companion paper to "The NANOGrav 15-year Data Set: Evidence for a Gravitational-Wave Background"; For questions or comments, please email comments@nanograv.org

The NANOGrav 15-year Data Set: Search for Anisotropy in the Gravitational-Wave Background

Authors: Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Bence Bécsy, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Megan E. DeCesar, Paul B. Demorest, Timothy Dolch, Brendan Drachler, Elizabeth C. Ferrara, William Fiore , et al. (68 additional authors not shown)

Abstract: The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) has reported evidence for the presence of an isotropic nanohertz gravitational wave background (GWB) in its 15 yr dataset. However, if the GWB is produced by a population of inspiraling supermassive black hole binary (SMBHB) systems, then the background is predicted to be anisotropic, depending on the distribution of these… ▽ More The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) has reported evidence for the presence of an isotropic nanohertz gravitational wave background (GWB) in its 15 yr dataset. However, if the GWB is produced by a population of inspiraling supermassive black hole binary (SMBHB) systems, then the background is predicted to be anisotropic, depending on the distribution of these systems in the local Universe and the statistical properties of the SMBHB population. In this work, we search for anisotropy in the GWB using multiple methods and bases to describe the distribution of the GWB power on the sky. We do not find significant evidence of anisotropy, and place a Bayesian $95\%$ upper limit on the level of broadband anisotropy such that $(C_{l>0} / C_{l=0}) < 20\%$. We also derive conservative estimates on the anisotropy expected from a random distribution of SMBHB systems using astrophysical simulations conditioned on the isotropic GWB inferred in the 15-yr dataset, and show that this dataset has sufficient sensitivity to probe a large fraction of the predicted level of anisotropy. We end by highlighting the opportunities and challenges in searching for anisotropy in pulsar timing array data. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 19 pages, 11 figures; submitted to Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org

The NANOGrav 15-year Data Set: Constraints on Supermassive Black Hole Binaries from the Gravitational Wave Background

Authors: Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Paul T. Baker, Bence Bécsy, Laura Blecha, Alexander Bonilla, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Rand Burnette, Robin Case, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Katerina Chatziioannou, Belinda D. Cheeseboro, Siyuan Chen, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Curt J. Cutler , et al. (89 additional authors not shown)

Abstract: The NANOGrav 15-year data set shows evidence for the presence of a low-frequency gravitational-wave background (GWB). While many physical processes can source such low-frequency gravitational waves, here we analyze the signal as coming from a population of supermassive black hole (SMBH) binaries distributed throughout the Universe. We show that astrophysically motivated models of SMBH binary popul… ▽ More The NANOGrav 15-year data set shows evidence for the presence of a low-frequency gravitational-wave background (GWB). While many physical processes can source such low-frequency gravitational waves, here we analyze the signal as coming from a population of supermassive black hole (SMBH) binaries distributed throughout the Universe. We show that astrophysically motivated models of SMBH binary populations are able to reproduce both the amplitude and shape of the observed low-frequency gravitational-wave spectrum. While multiple model variations are able to reproduce the GWB spectrum at our current measurement precision, our results highlight the importance of accurately modeling binary evolution for producing realistic GWB spectra. Additionally, while reasonable parameters are able to reproduce the 15-year observations, the implied GWB amplitude necessitates either a large number of parameters to be at the edges of expected values, or a small number of parameters to be notably different from standard expectations. While we are not yet able to definitively establish the origin of the inferred GWB signal, the consistency of the signal with astrophysical expectations offers a tantalizing prospect for confirming that SMBH binaries are able to form, reach sub-parsec separations, and eventually coalesce. As the significance grows over time, higher-order features of the GWB spectrum will definitively determine the nature of the GWB and allow for novel constraints on SMBH populations. △ Less

Submitted 18 July, 2023; v1 submitted 28 June, 2023; originally announced June 2023.

Comments: Accepted by Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org. Edited to fix two equation typos (Eq.13 & 21), and minor text typos

The NANOGrav 15-year Data Set: Search for Signals from New Physics

Authors: Adeela Afzal, Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Bence Bécsy, Jose Juan Blanco-Pillado, Laura Blecha, Kimberly K. Boddy, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Rand Burnette, Robin Case, Maria Charisi, Shami Chatterjee, Katerina Chatziioannou, Belinda D. Cheeseboro, Siyuan Chen, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie , et al. (98 additional authors not shown)

Abstract: The 15-year pulsar timing data set collected by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) shows positive evidence for the presence of a low-frequency gravitational-wave (GW) background. In this paper, we investigate potential cosmological interpretations of this signal, specifically cosmic inflation, scalar-induced GWs, first-order phase transitions, cosmic string… ▽ More The 15-year pulsar timing data set collected by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) shows positive evidence for the presence of a low-frequency gravitational-wave (GW) background. In this paper, we investigate potential cosmological interpretations of this signal, specifically cosmic inflation, scalar-induced GWs, first-order phase transitions, cosmic strings, and domain walls. We find that, with the exception of stable cosmic strings of field theory origin, all these models can reproduce the observed signal. When compared to the standard interpretation in terms of inspiraling supermassive black hole binaries (SMBHBs), many cosmological models seem to provide a better fit resulting in Bayes factors in the range from 10 to 100. However, these results strongly depend on modeling assumptions about the cosmic SMBHB population and, at this stage, should not be regarded as evidence for new physics. Furthermore, we identify excluded parameter regions where the predicted GW signal from cosmological sources significantly exceeds the NANOGrav signal. These parameter constraints are independent of the origin of the NANOGrav signal and illustrate how pulsar timing data provide a new way to constrain the parameter space of these models. Finally, we search for deterministic signals produced by models of ultralight dark matter (ULDM) and dark matter substructures in the Milky Way. We find no evidence for either of these signals and thus report updated constraints on these models. In the case of ULDM, these constraints outperform torsion balance and atomic clock constraints for ULDM coupled to electrons, muons, or gluons. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 74 pages, 31 figures, 4 tables; published in Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org

The NANOGrav 15-Year Data Set: Detector Characterization and Noise Budget

Authors: Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Bence Bécsy, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Maria Charisi, Shami Chatterjee, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Megan E. Decesar, Paul B. Demorest, Timothy Dolch, Brendan Drachler, Elizabeth C. Ferrara, William Fiore, Emmanuel Fonseca , et al. (66 additional authors not shown)

Abstract: Pulsar timing arrays (PTAs) are galactic-scale gravitational wave detectors. Each individual arm, composed of a millisecond pulsar, a radio telescope, and a kiloparsecs-long path, differs in its properties but, in aggregate, can be used to extract low-frequency gravitational wave (GW) signals. We present a noise and sensitivity analysis to accompany the NANOGrav 15-year data release and associated… ▽ More Pulsar timing arrays (PTAs) are galactic-scale gravitational wave detectors. Each individual arm, composed of a millisecond pulsar, a radio telescope, and a kiloparsecs-long path, differs in its properties but, in aggregate, can be used to extract low-frequency gravitational wave (GW) signals. We present a noise and sensitivity analysis to accompany the NANOGrav 15-year data release and associated papers, along with an in-depth introduction to PTA noise models. As a first step in our analysis, we characterize each individual pulsar data set with three types of white noise parameters and two red noise parameters. These parameters, along with the timing model and, particularly, a piecewise-constant model for the time-variable dispersion measure, determine the sensitivity curve over the low-frequency GW band we are searching. We tabulate information for all of the pulsars in this data release and present some representative sensitivity curves. We then combine the individual pulsar sensitivities using a signal-to-noise-ratio statistic to calculate the global sensitivity of the PTA to a stochastic background of GWs, obtaining a minimum noise characteristic strain of $7\times 10^{-15}$ at 5 nHz. A power law-integrated analysis shows rough agreement with the amplitudes recovered in NANOGrav's 15-year GW background analysis. While our phenomenological noise model does not model all known physical effects explicitly, it provides an accurate characterization of the noise in the data while preserving sensitivity to multiple classes of GW signals. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 67 pages, 73 figures, 3 tables; published in Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org

The NANOGrav 15-year Data Set: Observations and Timing of 68 Millisecond Pulsars

Authors: Gabriella Agazie, Md Faisal Alam, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Laura Blecha, Victoria Bonidie, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Bence Bécsy, Christopher Chapman, Maria Charisi, Shami Chatterjee, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Megan E. DeCesar, Paul B. Demorest, Timothy Dolch, Brendan Drachler , et al. (75 additional authors not shown)

Abstract: We present observations and timing analyses of 68 millisecond pulsars (MSPs) comprising the 15-year data set of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). NANOGrav is a pulsar timing array (PTA) experiment that is sensitive to low-frequency gravitational waves. This is NANOGrav's fifth public data release, including both "narrowband" and "wideband" time-of-arrival… ▽ More We present observations and timing analyses of 68 millisecond pulsars (MSPs) comprising the 15-year data set of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). NANOGrav is a pulsar timing array (PTA) experiment that is sensitive to low-frequency gravitational waves. This is NANOGrav's fifth public data release, including both "narrowband" and "wideband" time-of-arrival (TOA) measurements and corresponding pulsar timing models. We have added 21 MSPs and extended our timing baselines by three years, now spanning nearly 16 years for some of our sources. The data were collected using the Arecibo Observatory, the Green Bank Telescope, and the Very Large Array between frequencies of 327 MHz and 3 GHz, with most sources observed approximately monthly. A number of notable methodological and procedural changes were made compared to our previous data sets. These improve the overall quality of the TOA data set and are part of the transition to new pulsar timing and PTA analysis software packages. For the first time, our data products are accompanied by a full suite of software to reproduce data reduction, analysis, and results. Our timing models include a variety of newly detected astrometric and binary pulsar parameters, including several significant improvements to pulsar mass constraints. We find that the time series of 23 pulsars contain detectable levels of red noise, 10 of which are new measurements. In this data set, we find evidence for a stochastic gravitational-wave background. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 90 pages, 74 figures, 6 tables; published in Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org

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)

The NANOGrav 15-year Data Set: Evidence for a Gravitational-Wave Background

Authors: Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Bence Becsy, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Rand Burnette, Robin Case, Maria Charisi, Shami Chatterjee, Katerina Chatziioannou, Belinda D. Cheeseboro, Siyuan Chen, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Curt J. Cutler, Megan E. DeCesar , et al. (89 additional authors not shown)

Abstract: We report multiple lines of evidence for a stochastic signal that is correlated among 67 pulsars from the 15-year pulsar-timing data set collected by the North American Nanohertz Observatory for Gravitational Waves. The correlations follow the Hellings-Downs pattern expected for a stochastic gravitational-wave background. The presence of such a gravitational-wave background with a power-law-spectr… ▽ More We report multiple lines of evidence for a stochastic signal that is correlated among 67 pulsars from the 15-year pulsar-timing data set collected by the North American Nanohertz Observatory for Gravitational Waves. The correlations follow the Hellings-Downs pattern expected for a stochastic gravitational-wave background. The presence of such a gravitational-wave background with a power-law-spectrum is favored over a model with only independent pulsar noises with a Bayes factor in excess of $10^{14}$, and this same model is favored over an uncorrelated common power-law-spectrum model with Bayes factors of 200-1000, depending on spectral modeling choices. We have built a statistical background distribution for these latter Bayes factors using a method that removes inter-pulsar correlations from our data set, finding $p = 10^{-3}$ (approx. $3σ$) for the observed Bayes factors in the null no-correlation scenario. A frequentist test statistic built directly as a weighted sum of inter-pulsar correlations yields $p = 5 \times 10^{-5} - 1.9 \times 10^{-4}$ (approx. $3.5 - 4σ$). Assuming a fiducial $f^{-2/3}$ characteristic-strain spectrum, as appropriate for an ensemble of binary supermassive black-hole inspirals, the strain amplitude is $2.4^{+0.7}_{-0.6} \times 10^{-15}$ (median + 90% credible interval) at a reference frequency of 1/(1 yr). The inferred gravitational-wave background amplitude and spectrum are consistent with astrophysical expectations for a signal from a population of supermassive black-hole binaries, although more exotic cosmological and astrophysical sources cannot be excluded. The observation of Hellings-Downs correlations points to the gravitational-wave origin of this signal. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 30 pages, 18 figures. Published in Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org

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

Discovery and Timing of Millisecond Pulsars with the Arecibo 327 MHz Drift-Scan Survey

Authors: Evan F. Lewis, Timothy E. E. Olszanski, Julia S. Deneva, Paulo C. C. Freire, Maura A. McLaughlin, Kevin Stovall, Manjari Bagchi, Jose G. Martinez, Benetge B. P. Perera

Abstract: We present the discovery and timing solutions of four millisecond pulsars (MSPs) discovered in the Arecibo 327 MHz Drift-Scan Pulsar Survey. Three of these pulsars are in binary systems, consisting of a redback (PSR J2055+1545), a black widow (PSR J1630+3550), and a neutron star-white dwarf binary (PSR J2116+1345). The fourth MSP, PSR J2212+2450, is isolated. We present the multiyear timing soluti… ▽ More We present the discovery and timing solutions of four millisecond pulsars (MSPs) discovered in the Arecibo 327 MHz Drift-Scan Pulsar Survey. Three of these pulsars are in binary systems, consisting of a redback (PSR J2055+1545), a black widow (PSR J1630+3550), and a neutron star-white dwarf binary (PSR J2116+1345). The fourth MSP, PSR J2212+2450, is isolated. We present the multiyear timing solutions as well as polarization properties across a range of radio frequencies for each pulsar. We perform a multiwavelength search for emission from these systems and find an optical counterpart for PSR J2055+1545 in Gaia DR3, as well as a gamma-ray counterpart for PSR J2116+1345 with the Fermi-LAT telescope. Despite the close colocation of PSR J2055+1545 with a Fermi source, we are unable to detect gamma-ray pulsations, likely due to the large orbital variability of the system. This work presents the first two binaries found by this survey with orbital periods shorter than a day; we expect to find more in the 40% of the survey data that have yet to be searched. △ Less

Submitted 19 October, 2023; v1 submitted 16 June, 2023; originally announced June 2023.

Comments: 22 pages, 11 figures, 4 tables. Published in ApJ

Searching for continuous Gravitational Waves in the second data release of the International Pulsar Timing Array

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

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

Submitted 19 March, 2023; originally announced March 2023.

An unusual pulse shape change event in PSR J1713+0747 observed with the Green Bank Telescope and CHIME

Authors: Ross J. Jennings, James M. Cordes, Shami Chatterjee, Maura A. McLaughlin, Paul B. Demorest, Zaven Arzoumanian, Paul T. Baker, Harsha Blumer, Paul R. Brook, Tyler Cohen, Fronefield Crawford, H. Thankful Cromartie, Megan E. DeCesar, Timothy Dolch, Elizabeth C. Ferrara, Emmanuel Fonseca, Deborah C. Good, Jeffrey S. Hazboun, Megan L. Jones, David L. Kaplan, Michael T. Lam, T. Joseph W. Lazio, Duncan R. Lorimer, Jing Luo, Ryan S. Lynch , et al. (19 additional authors not shown)

Abstract: The millisecond pulsar J1713+0747 underwent a sudden and significant pulse shape change between April 16 and 17, 2021 (MJDs 59320 and 59321). Subsequently, the pulse shape gradually recovered over the course of several months. We report the results of continued multi-frequency radio observations of the pulsar made using the Canadian Hydrogen Intensity Mapping Experiment (CHIME) and the 100-meter G… ▽ More The millisecond pulsar J1713+0747 underwent a sudden and significant pulse shape change between April 16 and 17, 2021 (MJDs 59320 and 59321). Subsequently, the pulse shape gradually recovered over the course of several months. We report the results of continued multi-frequency radio observations of the pulsar made using the Canadian Hydrogen Intensity Mapping Experiment (CHIME) and the 100-meter Green Bank Telescope (GBT) in a three-year period encompassing the shape change event, between February 2020 and February 2023. As of February 2023, the pulse shape had returned to a state similar to that seen before the event, but with measurable changes remaining. The amplitude of the shape change and the accompanying TOA residuals display a strong non-monotonic dependence on radio frequency, demonstrating that the event is neither a glitch (the effects of which should be independent of radio frequency, $ν$) nor a change in dispersion measure (DM) alone (which would produce a delay proportional to $ν^{-2}$). However, it does bear some resemblance to the two previous "chromatic timing events" observed in J1713+0747 (Demorest et al. 2013; Lam et al. 2016), as well as to a similar event observed in PSR J1643-1224 in 2015 (Shannon et al. 2016). △ Less

Submitted 31 January, 2024; v1 submitted 21 October, 2022; originally announced October 2022.

Comments: 18 pages, 8 figures. Submitted to ApJ. Data available at https://doi.org/10.5281/zenodo.7236459

A Parkes "Murriyang" Search for Pulsars and Transients in the Large Magellanic Cloud

Authors: Shinnosuke Hisano, Fronefield Crawford, Victoria Bonidie, Md F. Alam, Keitaro Takahashi, Duncan R. Lorimer, Josh P. Ridley, Maura M. McLaughlin, Benetge B. P. Perera

Abstract: We report on the latest results of a Parkes multibeam survey for pulsars and dispersed radio bursts in the Large Magellanic Cloud (LMC). We conducted both periodicity and single-pulse searches at a much larger range of trial dispersion measures (DMs) than previously searched. We detected 229 single pulses with signal-to-noise ratio ($\rm{S/N}) > 7$ that were classified by the deep learning network… ▽ More We report on the latest results of a Parkes multibeam survey for pulsars and dispersed radio bursts in the Large Magellanic Cloud (LMC). We conducted both periodicity and single-pulse searches at a much larger range of trial dispersion measures (DMs) than previously searched. We detected 229 single pulses with signal-to-noise ratio ($\rm{S/N}) > 7$ that were classified by the deep learning network FETCH as being real (with $> 90$% likelihood), of which 9 were from the known giant-pulse-emitting pulsar PSR B0540$-$69. Two possibly repeating sources were detected with DMs suggesting that they lie within the LMC, but these require confirmation. Only 3 of the 220 unknown pulses had S/N greater than 8, and the DM distribution for these detected pulses follows an exponential falloff with increasing DM and does not show any excess of signals at DM values expected for the LMC. These features suggest that the detected pulses are not likely to be real, although they are visually compelling. We also report the discovery of a new pulsar (PSR J0556$-$67) in our periodicity search. This pulsar has a spin period of 791 ms, a DM of 71 cm$^{-3}$ pc, an estimated 1400 MHz flux density of $\sim 0.12$ mJy, and shows no evidence of binary motion. Future observations may be able to confirm whether any of the weak but promising candidates in our single pulse and periodicity searches, including our two possible repeaters, are real or not. △ Less

Submitted 23 February, 2022; v1 submitted 22 February, 2022; originally announced February 2022.

Comments: 17 pages, 11 figures

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

Single-pulse studies of three millisecond pulsars

Authors: N. T. Palliyaguru, B. B. P. Perera, M. A. McLaughlin, S. Oslowski, G. L. Siebert

Abstract: Single-pulse studies are important to understand the pulsar emission mechanism and the noise floor in precision timing. We study total intensity and polarimetry properties of three bright millisecond pulsars - PSRs J1022+1001, J1713+0747, and B1855+09 - that have detectable single pulses at multiple frequencies. We report for the first time the detection of single pulses from PSRs J1022+1001 and J… ▽ More Single-pulse studies are important to understand the pulsar emission mechanism and the noise floor in precision timing. We study total intensity and polarimetry properties of three bright millisecond pulsars - PSRs J1022+1001, J1713+0747, and B1855+09 - that have detectable single pulses at multiple frequencies. We report for the first time the detection of single pulses from PSRs J1022+1001 and J1713+0747 at 4.5 GHz. In addition, for those two pulsars the fraction of linear polarization in the average profile is significantly reduced at 4.5 GHz, compared to 1.38 GHz, which could support the expected deviation from a dipolar field closer to the pulsar surface. There is a hint of orthogonal modes in the single pulses of PSR J1713+0747. More sensitive multi-frequency observations may be useful to confirm these findings. The jitter noise contributions at 1.38 GHz, scaled to one hour, for PSRs J1022+1001, J1713+0747 and B1855+09 are ~135 ns, ~45 ns, and ~60 ns respectively and are consistent with previous studies. We also show that selective bright-pulse timing of PSR J1022+1001 yields improved root-mean-square residuals of ~22 $μ$s, which is a factor of ~3 better than timing using single pulses alone. △ Less

Submitted 16 January, 2023; v1 submitted 24 December, 2021; originally announced January 2022.

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

Strong-field Gravity Tests with the Double Pulsar

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

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

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

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

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

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

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

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

Submitted 9 November, 2021; originally announced November 2021.

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

Search for fast radio transients using Arecibo drift-scan observations at 1.4 GHz

Authors: B. B. P. Perera, A. J. Smith, S. Vaddi, R. Carballo-Rubio, A. McGilvray, A. Venkataraman, D. Anish Roshi, P. K. Manoharan, P. Perillat, E. Lieb, D. R. Lorimer, M. A. McLaughlin, D. Agarwal, K. Aggarwal, S. M. Ransom

Abstract: We conducted a drift-scan observation campaign using the 305-m Arecibo telescope in January and March 2020 when the observatory was temporarily closed during the intense earthquakes and the initial outbreak of the COVID-19 pandemic, respectively. The primary objective of the survey was to search for fast radio transients, including Fast Radio Bursts (FRBs) and Rotating Radio Transients (RRATs). We… ▽ More We conducted a drift-scan observation campaign using the 305-m Arecibo telescope in January and March 2020 when the observatory was temporarily closed during the intense earthquakes and the initial outbreak of the COVID-19 pandemic, respectively. The primary objective of the survey was to search for fast radio transients, including Fast Radio Bursts (FRBs) and Rotating Radio Transients (RRATs). We used the 7-beam ALFA receiver to observe different sections of the sky within the declination region $\sim$(10$-$20) deg on 23 nights and collected 160 hours of data in total. We searched our data for single-pulse transients, covering up to a maximum dispersion measure of 11 000 pc cm$^{-3}$ at which the dispersion delay across the entire bandwidth is equal to the 13 s transit length of our observations. The analysis produced more than 18 million candidates. Machine learning techniques sorted the radio frequency interference and possibly astrophysical candidates, allowing us to visually inspect and confirm the candidate transients. We found no evidence for new astrophysical transients in our data. We also searched for emission from repeated transient signals, but found no evidence for such sources. We detected single pulses from two known pulsars in our observations and their measured flux densities are consistent with the expected values. Based on our observations and sensitivity, we estimated the upper limit for the FRB rate to be $<$2.8$\times10^5$ sky$^{-1}$ day$^{-1}$ above a fluence of 0.16 Jy ms at 1.4 GHz, which is consistent with the rates from other telescopes and surveys. △ Less

Submitted 27 October, 2021; originally announced October 2021.

Comments: 11 pages, 8 figures. Accepted 2021 October 27 for publication in MNRAS

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

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

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

Submitted 25 October, 2021; originally announced October 2021.

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

Study of 72 pulsars discovered in the PALFA survey: Timing analysis, glitch activity, emission variability, and a pulsar in an eccentric binary

Authors: E. Parent, H. Sewalls, P. C. C. Freire, T. Matheny, A. G. Lyne, B. B. P. Perera, F. Cardoso, M. A. McLaughlin, B. Allen, A. Brazier, F. Camilo, S. Chatterjee, J. M. Cordes, F. Crawford, J. S. Deneva, F. A. Dong, R. D. Ferdman, E. Fonseca, J. W. T. Hessels, V. M. Kaspi, B. Knispel, J. van Leeuwen, R. S. Lynch, B. M. Meyers, J. W. McKee , et al. (9 additional authors not shown)

Abstract: We present new discoveries and results from long-term timing of 72 pulsars discovered in the Arecibo PALFA survey, including precise determination of astrometric and spin parameters, and flux density and scatter broadening measurements at 1.4 GHz. Notable discoveries include two young pulsars (characteristic ages $\sim$30 kyr) with no apparent supernova remnant associations, three mode changing, 1… ▽ More We present new discoveries and results from long-term timing of 72 pulsars discovered in the Arecibo PALFA survey, including precise determination of astrometric and spin parameters, and flux density and scatter broadening measurements at 1.4 GHz. Notable discoveries include two young pulsars (characteristic ages $\sim$30 kyr) with no apparent supernova remnant associations, three mode changing, 12 nulling and two intermittent pulsars. We detected eight glitches in five pulsars. Among them is PSR J1939+2609, an apparently old pulsar (characteristic age $\sim$1 Gy), and PSR J1954+2529, which likely belongs to a newly-emerging class of binary pulsars. The latter is the only pulsar among the 72 that is clearly not isolated: a non-recycled neutron star with a 931-ms spin period in an eccentric ($e\,=\,0.114$) wide ($P_b\,=\,82.7\,$d) orbit with a companion of undetermined nature having a minimum mass of $\sim0.6\,M_{\odot}$. Since operations at Arecibo ceased in 2020 August, we give a final tally of PALFA sky coverage, and compare its 207 pulsar discoveries to the known population. On average, they are 50% more distant than other Galactic plane radio pulsars; PALFA millisecond pulsars (MSP) have twice the dispersion measure per unit spin period than the known population of MSP in the Plane. The four intermittent pulsars discovered by PALFA more than double the population of such objects, which should help to improve our understanding of pulsar magnetosphere physics. The statistics for these, RRATS, and nulling pulsars suggest that there are many more of these objects in the Galaxy than was previously thought. △ Less

Submitted 17 November, 2021; v1 submitted 4 August, 2021; originally announced August 2021.

Comments: 39 pages, 18 figures, 9 tables. Accepted for publication in ApJ

Multi-frequency study of the peculiar pulsars PSR B0919+06 and PSR B1859+07

Authors: K. M. Rajwade, B. B. P. Perera, B. W. Stappers, J. Roy, A. Karastergiou, J. M. Rankin

Abstract: Since their discovery more than 50 years ago, broadband radio studies of pulsars have generated a wealth of information about the underlying physics of radio emission. In order to gain some further insights into this elusive emission mechanism, we performed a multi-frequency study of two very well-known pulsars, PSR~B0919+06 and PSR~B1859+07. These pulsars show peculiar radio emission properties w… ▽ More Since their discovery more than 50 years ago, broadband radio studies of pulsars have generated a wealth of information about the underlying physics of radio emission. In order to gain some further insights into this elusive emission mechanism, we performed a multi-frequency study of two very well-known pulsars, PSR~B0919+06 and PSR~B1859+07. These pulsars show peculiar radio emission properties whereby the emission shifts to an earlier rotation phase before returning to the nominal emission phase in a few tens of pulsar rotations (also known as `swooshes'). We confirm the previous claim that the emission during the swoosh is not necessarily absent at low frequencies and the single pulses during a swoosh show varied behaviour at 220~MHz. We also confirm that in PSR~B0919+06, the pulses during the swoosh show a chromatic dependence of the maximum offset from the normal emission phase with the offset following a consistent relationship with observing frequency. We also observe that the flux density spectrum of the radio profile during the swoosh is inverted compared to the normal emission. For PSR~B1859+07, we have discovered a new mode of emission in the pulsar that is potentially quasi-periodic with a different periodicity than is seen in its swooshes. We invoke an emission model previously proposed in the literature and show that this simple model can explain the macroscopic observed characteristics in both pulsars. We also argue that pulsars that exhibit similar variability on short timescales may have the same underlying emission mechanism. △ Less

Submitted 5 July, 2021; originally announced July 2021.

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

Detection of OH 18-cm Emission from Comet C/2020 F3 NEOWISE using the Arecibo Telescope

Authors: Allison J. Smith, D. Anish Roshi, Periasamy Manoharan, Sravani Vaddi, Benetge B. P. Perera, Anna McGilvray

Abstract: We report the detection of emission from the OH 18 cm $Λ$-doublet transitions toward Comet C/2020 F3 NEOWISE using the Arecibo Telescope. The antenna temperatures are 113$\pm$3 mK for the 1667 MHz line and 57$\pm$3 mK for the 1665 MHz line. The beam averaged OH column density (centered on the comet nucleus) derived from the 1667 transition is $N_{OH}$=1.11$\pm0.06\times10^{13}$ cm$^{-2}$. We imple… ▽ More We report the detection of emission from the OH 18 cm $Λ$-doublet transitions toward Comet C/2020 F3 NEOWISE using the Arecibo Telescope. The antenna temperatures are 113$\pm$3 mK for the 1667 MHz line and 57$\pm$3 mK for the 1665 MHz line. The beam averaged OH column density (centered on the comet nucleus) derived from the 1667 transition is $N_{OH}$=1.11$\pm0.06\times10^{13}$ cm$^{-2}$. We implemented the Haser model to derive an OH production rate. The estimated OH production rate using the 1667 MHz transition is Q$_{OH}$=3.6$\pm0.6\times10^{28}$ s$^{-1}$, a factor of 2.4 lower than optically derived values for the same observing day, the difference of which is likely explained by quenching. △ Less

Submitted 7 May, 2021; originally announced May 2021.

Comments: 6 pages, 3 figures

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

The Future Of The Arecibo Observatory: The Next Generation Arecibo Telescope

Authors: D. Anish Roshi, N. Aponte, E. Araya, H. Arce, L. A. Baker, W. Baan, T. M. Becker, J. K. Breakall, R. G. Brown, C. G. M. Brum, M. Busch, D. B. Campbell, T. Cohen, F. Cordova, J. S. Deneva, M. Devogele, T. Dolch, F. O. Fernandez-Rodriguez, T. Ghosh, P. F. Goldsmith, L. I. Gurvits, M. Haynes, C. Heiles, J. W. T. Hessel, D. Hickson , et al. (49 additional authors not shown)

Abstract: The Arecibo Observatory (AO) is a multidisciplinary research and education facility that is recognized worldwide as a leading facility in astronomy, planetary, and atmospheric and space sciences. AO's cornerstone research instrument was the 305-m William E. Gordon telescope. On December 1, 2020, the 305-m telescope collapsed and was irreparably damaged. In the three weeks following the collapse, A… ▽ More The Arecibo Observatory (AO) is a multidisciplinary research and education facility that is recognized worldwide as a leading facility in astronomy, planetary, and atmospheric and space sciences. AO's cornerstone research instrument was the 305-m William E. Gordon telescope. On December 1, 2020, the 305-m telescope collapsed and was irreparably damaged. In the three weeks following the collapse, AO's scientific and engineering staff and the AO users community initiated extensive discussions on the future of the observatory. The community is in overwhelming agreement that there is a need to build an enhanced, next-generation radar-radio telescope at the AO site. From these discussions, we established the set of science requirements the new facility should enable. These requirements can be summarized briefly as: 5 MW of continuous wave transmitter power at 2 - 6 GHz, 10 MW of peak transmitter power at 430 MHz (also at 220MHz under consideration), zenith angle coverage 0 to 48 deg, frequency coverage 0.2 to 30 GHz and increased Field-of-View. These requirements determine the unique specifications of the proposed new instrument. The telescope design concept we suggest consists of a compact array of fixed dishes on a tiltable, plate-like structure with a collecting area equivalent to a 300m dish. This concept, referred to as the Next Generation Arecibo Telescope (NGAT), meets all of the desired specifications and provides significant new science capabilities to all three research groups at AO. This whitepaper presents a sample of the wide variety of the science that can be achieved with the NGAT, the details of the telescope design concept and the need for the new telescope to be located at the AO site. We also discuss other AO science activities that interlock with the NGAT in the white paper. △ Less

Submitted 1 April, 2021; v1 submitted 1 March, 2021; originally announced March 2021.

Comments: 82 pages (executive summary 10 pages), 21 figures, Arecibo observatory white paper (Updated with the complete author list and minor edits)

Understanding and improving the timing of PSR J0737-3039B

Authors: Aristeidis Noutsos, G. Desvignes, M. Kramer, N. Wex, P. C. C. Freire, I. H. Stairs, M. A. McLaughlin, R. N. Manchester, A. Possenti, M. Burgay, A. G. Lyne, R. P. Breton, B. B. P. Perera, R. D. Ferdman

Abstract: The double pulsar (PSR J0737-3039A/B) provides some of the most stringent tests of general relativity (GR) and its alternatives. The success of this system in tests of GR is largely due to the high-precision, long-term timing of its recycled-pulsar member, pulsar A. On the other hand, pulsar B is a young pulsar that exhibits significant short-term and long-term timing variations due to the electro… ▽ More The double pulsar (PSR J0737-3039A/B) provides some of the most stringent tests of general relativity (GR) and its alternatives. The success of this system in tests of GR is largely due to the high-precision, long-term timing of its recycled-pulsar member, pulsar A. On the other hand, pulsar B is a young pulsar that exhibits significant short-term and long-term timing variations due to the electromagnetic-wind interaction with its companion and geodetic precession. Improving pulsar B's timing precision is a key step towards improving the precision in a number of GR tests with PSR J0737-3039A/B. In this paper, red noise signatures in the timing of pulsar B are investigated using roughly a four-year time span, from 2004 to 2008, beyond which time the pulsar's radio beam precessed out of view ... The timing of pulsar B presented in this paper depends on the size of the pulsar's orbit, which was calculated from GR, in order to precisely account for orbital timing delays. Consequently, our timing cannot directly be used to test theories of gravity. However, our modelling of the beam shape and radial wind of pulsar B can indirectly aid future efforts to time this pulsar by constraining part of the additional red noise observed on top of the orbital delays. As such, we conclude that, in the idealised case of zero covariance between our model's parameters and those of the timing model, our model can bring about a factor 2.6 improvement on the measurement precision of the mass ratio, R = mA/mB, between the two pulsars: a theory-independent parameter, which is pivotal in tests of GR. △ Less

Submitted 3 November, 2020; originally announced November 2020.

Comments: 43 pages, 28 figures, Accepted in A&A

Journal ref: A&A 643, A143 (2020)

Asymmetric mass ratios for bright double neutron-star mergers

Authors: R. D. Ferdman, P. C. C. Freire, B. B. P. Perera, N. Pol, F. Camilo, S. Chatterjee, J. M. Cordes, F. Crawford, J. W. T. Hessels, V. M. Kaspi, M. A. McLaughlin, E. Parent, I. H. Stairs, J. van Leeuwen

Abstract: The discovery of a radioactively powered kilonova associated with the binary neutron star merger GW170817 was the first - and still only - confirmed electromagnetic counterpart to a gravitational-wave event. However, observations of late-time electromagnetic emission are in tension with the expectations from standard neutron-star merger models. Although the large measured ejecta mass is potentiall… ▽ More The discovery of a radioactively powered kilonova associated with the binary neutron star merger GW170817 was the first - and still only - confirmed electromagnetic counterpart to a gravitational-wave event. However, observations of late-time electromagnetic emission are in tension with the expectations from standard neutron-star merger models. Although the large measured ejecta mass is potentially explained by a progenitor system that is asymmetric in terms of the stellar component masses, i.e. with a mass ratio $q$ of 0.7-0.8, the known Galactic population of merging double neutron star (DNS) systems (i.e. those that will coalesce within billions of years or less) has, until now, only consisted of nearly equal-mass ($q > 0.9$) binaries. PSR J1913+1102 is a DNS system in a 5-hour, low-eccentricity ($e = 0.09$) orbit, implying an orbital separation of 1.8 solar radii, with the two neutron stars predicted to coalesce in 470 million years due to gravitational-wave emission. Here we report that the masses of the two neutron stars, as measured by a dedicated pulsar timing campaign, are $1.62 \pm 0.03$ and $1.27 \pm 0.03$ solar masses for the pulsar and companion neutron star, respectively; with a measured mass ratio $q = 0.78 \pm 0.03$, it is the most asymmetric DNS among known merging systems. Based on this detection, our population synthesis analysis implies that such asymmetric binaries represent between 2 and 30% (90% confidence) of the total population of merging DNS binaries. The coalescence of a member of this population offers a possible explanation for the anomalous properties of GW170817, including the observed kilonova emission from that event. △ Less

Submitted 8 July, 2020; originally announced July 2020.

Comments: 17 pages, 3 figures, published in Nature on 9 July 2020 (https://www.nature.com/articles/s41586-020-2439-x); this is the authors' final submitted version before review. Includes methods and supplementary information

The NANOGrav 11-year Data Set: Constraints on Planetary Masses Around 45 Millisecond Pulsars

Authors: E. A. Behrens, S. M. Ransom, D. R. Madison, Z. Arzoumanian, K. Crowter, M. E. DeCesar, P. B. Demorest, T. Dolch, J. A. Ellis, R. D. Ferdman, E. C. Ferrara, E. Fonseca, P. A. Gentile, G. Jones, M. L. Jones, M. T. Lam, L. Levin, D. R. Lorimer, R. S. Lynch, M. A. McLaughlin, C. Ng, D. J. Nice, T. T. Pennucci, B. B. P. Perera, P. S. Ray , et al. (5 additional authors not shown)

Abstract: We search for extrasolar planets around millisecond pulsars using pulsar timing data and seek to determine the minimum detectable planetary masses as a function of orbital period. Using the 11-year data set from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), we look for variations from our models of pulse arrival times due to the presence of exoplanets. No planets are… ▽ More We search for extrasolar planets around millisecond pulsars using pulsar timing data and seek to determine the minimum detectable planetary masses as a function of orbital period. Using the 11-year data set from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), we look for variations from our models of pulse arrival times due to the presence of exoplanets. No planets are detected around the millisecond pulsars in the NANOGrav 11-year data set, but taking into consideration the noise levels of each pulsar and the sampling rate of our observations, we develop limits that show we are sensitive to planetary masses as low as that of the moon. We analyzed potential planet periods, P, in the range 7 days < P < 2000 days, with somewhat smaller ranges for some binary pulsars. The planetary mass limit for our median-sensitivity pulsar within this period range is 1 M_moon (P / 100 days)^(-2/3). △ Less

Submitted 24 March, 2020; v1 submitted 1 December, 2019; originally announced December 2019.

Comments: Revised and accepted by ApJ Letters

A pulsar-based timescale from the International Pulsar Timing Array

Authors: G. Hobbs, L. Guo, R. N. Caballero, W. Coles, K. J. Lee, R. N. Manchester, D. J. Reardon, D. Matsakis, M. L. Tong, Z. Arzoumanian, M. Bailes, C. G. Bassa, N. D. R. Bhat, A. Brazier, S. Burke-Spolaor, D. J. Champion, S. Chatterjee, I. Cognard, S. Dai, G. Desvignes, T. Dolch, R. D. Ferdman, E. Graikou, L. Guillemot, G. H. Janssen , et al. (34 additional authors not shown)

Abstract: We have constructed a new timescale, TT(IPTA16), based on observations of radio pulsars presented in the first data release from the International Pulsar Timing Array (IPTA). We used two analysis techniques with independent estimates of the noise models for the pulsar observations and different algorithms for obtaining the pulsar timescale. The two analyses agree within the estimated uncertainties… ▽ More We have constructed a new timescale, TT(IPTA16), based on observations of radio pulsars presented in the first data release from the International Pulsar Timing Array (IPTA). We used two analysis techniques with independent estimates of the noise models for the pulsar observations and different algorithms for obtaining the pulsar timescale. The two analyses agree within the estimated uncertainties and both agree with TT(BIPM17), a post-corrected timescale produced by the Bureau International des Poids et Mesures (BIPM). We show that both methods could detect significant errors in TT(BIPM17) if they were present. We estimate the stability of the atomic clocks from which TT(BIPM17) is derived using observations of four rubidium fountain clocks at the US Naval Observatory. Comparing the power spectrum of TT(IPTA16) with that of these fountain clocks suggests that pulsar-based timescales are unlikely to contribute to the stability of the best timescales over the next decade, but they will remain a valuable independent check on atomic timescales. We also find that the stability of the pulsar-based timescale is likely to be limited by our knowledge of solar-system dynamics, and that errors in TT(BIPM17) will not be a limiting factor for the primary goal of the IPTA, which is to search for the signatures of nano-Hertz gravitational waves. △ Less

Submitted 29 October, 2019; originally announced October 2019.

Comments: accepted by MNRAS

The International Pulsar Timing Array: Second data release

Authors: B. B. P. Perera, M. E. DeCesar, P. B. Demorest, M. Kerr, L. Lentati, D. J. Nice, S. Oslowski, S. M. Ransom, M. J. Keith, Z. Arzoumanian, M. Bailes, P. T. Baker, C. G. Bassa, N. D. R. Bhat, A. Brazier, M. Burgay, S. Burke-Spolaor, R. N. Caballero, D. J. Champion, S. Chatterjee, S. Chen, I. Cognard, J. M. Cordes, K. Crowter, S. Dai , et al. (50 additional authors not shown)

Abstract: In this paper, we describe the International Pulsar Timing Array second data release, which includes recent pulsar timing data obtained by three regional consortia: the European Pulsar Timing Array, the North American Nanohertz Observatory for Gravitational Waves, and the Parkes Pulsar Timing Array. We analyse and where possible combine high-precision timing data for 65 millisecond pulsars which a… ▽ More In this paper, we describe the International Pulsar Timing Array second data release, which includes recent pulsar timing data obtained by three regional consortia: the European Pulsar Timing Array, the North American Nanohertz Observatory for Gravitational Waves, and the Parkes Pulsar Timing Array. We analyse and where possible combine high-precision timing data for 65 millisecond pulsars which are regularly observed by these groups. A basic noise analysis, including the processes which are both correlated and uncorrelated in time, provides noise models and timing ephemerides for the pulsars. We find that the timing precisions of pulsars are generally improved compared to the previous data release, mainly due to the addition of new data in the combination. The main purpose of this work is to create the most up-to-date IPTA data release. These data are publicly available for searches for low-frequency gravitational waves and other pulsar science. △ Less

Submitted 10 September, 2019; originally announced September 2019.

Comments: Submitted to MNRAS and in review, 23 pages, 5 figures

Astro2020 Activities and Projects White Paper: Arecibo Observatory in the Next Decade

Authors: D. Anish Roshi, L. D. Anderson, E. Araya, D. Balser, W. Brisken, C. Brum, D. Campbell, S. Chatterjee, E. Churchwell, J. Condon, J. Cordes, F. Cordova, Y. Fernandez, J. Gago, T. Ghosh, P. F. Goldsmith, C. Heiles, D. Hickson, B. Jeffs, K. M. Jones, J. Lautenbach, B. M. Lewis, R. S. Lynch, P. K. Manoharan, S. Marshall , et al. (20 additional authors not shown)

Abstract: The white paper discusses Arecibo Observatory's plan for facility improvements and activities over the next decade. The facility improvements include: (a) improving the telescope surface, pointing and focusing to achieve superb performance up to ~12.5 GHz; (b) equip the telescope with ultrawide-band feeds; (c) upgrade the instrumentation with a 4 GHz bandwidth high dynamic range digital link and a… ▽ More The white paper discusses Arecibo Observatory's plan for facility improvements and activities over the next decade. The facility improvements include: (a) improving the telescope surface, pointing and focusing to achieve superb performance up to ~12.5 GHz; (b) equip the telescope with ultrawide-band feeds; (c) upgrade the instrumentation with a 4 GHz bandwidth high dynamic range digital link and a universal backend and (d) augment the VLBI facility by integrating the 12m telescope for phase referencing. These upgrades to the Arecibo telescope are critical to keep the national facility in the forefront of research in radio astronomy while maintaining its dominance in radar studies of near-Earth asteroids, planets and satellites. In the next decade, the Arecibo telescope will play a synergistic role with the upcoming facilities such as ngVLA, SKA and the now commissioned FAST telescope. Further, the observatory will be actively engaged in mentoring and training programs for students from a diverse background. △ Less

Submitted 13 July, 2019; originally announced July 2019.

Comments: Submitted as an Activities and Projects White Paper for the Astro2020 decadal survey (13 pages, 3 figures, 1 table)

The dynamics of Galactic centre pulsars: constraining pulsar distances and intrinsic spin-down

Authors: B. B. P. Perera, E. D. Barr, M. B. Mickaliger, A. G. Lyne, D. R. Lorimer, B. W. Stappers, R. P. Eatough, M. Kramer, C. Ng, R. Spiewak, M. Bailes, D. J. Champion, V. Morello, A. Possenti

Abstract: Through high-precision radio timing observations, we show that five recycled pulsars in the direction of the Galactic Centre (GC) have anomalous spin period time derivative ($\dot P$) measurements -- PSRs J1748$-$3009, J1753$-$2819, J1757$-$2745, and J1804$-$2858 show negative values of $\dot P$ and PSR J1801$-$3210 is found to have an exceptionally small value of $\dot P$. We attribute these obse… ▽ More Through high-precision radio timing observations, we show that five recycled pulsars in the direction of the Galactic Centre (GC) have anomalous spin period time derivative ($\dot P$) measurements -- PSRs J1748$-$3009, J1753$-$2819, J1757$-$2745, and J1804$-$2858 show negative values of $\dot P$ and PSR J1801$-$3210 is found to have an exceptionally small value of $\dot P$. We attribute these observed $\dot P$ measurements to acceleration of these pulsars along their lines-of-sight (LOSs) due to the Galactic gravitational field. Using models of the Galactic mass distribution and pulsar velocities, we constrain the distances to these pulsars, placing them on the far-side of the Galaxy, providing the first accurate distance measurements to pulsars located in this region and allowing us to consider the electron density along these LOSs. We find the new electron density model YMW16 to be more consistent with these observations than the previous model NE2001. The LOS dynamics further constrain the model-dependent intrinsic $\dot P$ values for these pulsars and they are consistent with measurements for other known pulsars. In the future, the independent distance measurements to these and other pulsars near the GC would allow us to constrain the Galactic gravitational potential more accurately. △ Less

Submitted 25 April, 2019; originally announced April 2019.

Comments: 17 pages, 11 figures. Accepted 2019 April 25 for publication in MNRAS

Studying the solar system with the International Pulsar Timing Array

Authors: R. N. Caballero, Y. J. Guo, K. J. Lee, P. Lazarus, D. J. Champion, G. Desvignes, M. Kramer, K. Plant, Z. Arzoumanian, M. Bailes, C. G. Bassa, N. D. R. Bhat, A. Brazier, M. Burgay, S. Burke-Spolaor, S. J. Chamberlin, S. Chatterjee, I. Cognard, J. M. Cordes, S. Dai, P. Demorest, T. Dolch, R. D. Ferdman, E. Fonseca, J. R. Gair , et al. (55 additional authors not shown)

Abstract: Pulsar-timing analyses are sensitive to errors in the solar-system ephemerides (SSEs) that timing models utilise to estimate the location of the solar-system barycentre, the quasi-inertial reference frame to which all recorded pulse times-of-arrival are referred. Any error in the SSE will affect all pulsars, therefore pulsar timing arrays (PTAs) are a suitable tool to search for such errors and im… ▽ More Pulsar-timing analyses are sensitive to errors in the solar-system ephemerides (SSEs) that timing models utilise to estimate the location of the solar-system barycentre, the quasi-inertial reference frame to which all recorded pulse times-of-arrival are referred. Any error in the SSE will affect all pulsars, therefore pulsar timing arrays (PTAs) are a suitable tool to search for such errors and impose independent constraints on relevant physical parameters. We employ the first data release of the International Pulsar Timing Array to constrain the masses of the planet-moons systems and to search for possible unmodelled objects (UMOs) in the solar system. We employ ten SSEs from two independent research groups, derive and compare mass constraints of planetary systems, and derive the first PTA mass constraints on asteroid-belt objects. Constraints on planetary-system masses have been improved by factors of up to 20 from the previous relevant study using the same assumptions, with the mass of the Jovian system measured at 9.5479189(3)$\times10^{-4}$ $M_{\odot}$. The mass of the dwarf planet Ceres is measured at 4.7(4)$\times10^{-10}$ $M_{\odot}$. We also present the first sensitivity curves using real data that place generic limits on the masses of UMOs, which can also be used as upper limits on the mass of putative exotic objects. For example, upper limits on dark-matter clumps are comparable to published limits using independent methods. While the constraints on planetary masses derived with all employed SSEs are consistent, we note and discuss differences in the associated timing residuals and UMO sensitivity curves. △ Less

Submitted 27 September, 2018; originally announced September 2018.

Comments: Accepted for publication by the Monthly Notices of the Royal Astronomical Society

Improving timing sensitivity in the microhertz frequency regime: limits from PSR J1713$+$0747 on gravitational waves produced by super-massive black-hole binaries

Authors: B. B. P. Perera, B. W. Stappers, S. Babak, M. J. Keith, J. Antoniadis, C. G. Bassa, R. N. Caballero, D. J. Champion, I. Cognard, G. Desvignes, E. Graikou, L. Guillemot, G. H. Janssen, R. Karuppusamy, M. Kramer, P. Lazarus, L. Lentati, K. Liu, A. G. Lyne, J. W. McKee, S. Oslowski, D. Perrodin, S. A. Sanidas, A. Sesana, G. Shaifullah , et al. (3 additional authors not shown)

Abstract: We search for continuous gravitational waves (CGWs) produced by individual super-massive black-hole binaries (SMBHBs) in circular orbits using high-cadence timing observations of PSR J1713$+$0747. We observe this millisecond pulsar using the telescopes in the European Pulsar Timing Array (EPTA) with an average cadence of approximately 1.6 days over the period between April 2011 and July 2015, incl… ▽ More We search for continuous gravitational waves (CGWs) produced by individual super-massive black-hole binaries (SMBHBs) in circular orbits using high-cadence timing observations of PSR J1713$+$0747. We observe this millisecond pulsar using the telescopes in the European Pulsar Timing Array (EPTA) with an average cadence of approximately 1.6 days over the period between April 2011 and July 2015, including an approximately daily average between February 2013 and April 2014. The high-cadence observations are used to improve the pulsar timing sensitivity across the GW frequency range of $0.008-5$ $μ$Hz. We use two algorithms in the analysis, including a spectral fitting method and a Bayesian approach. For an independent comparison, we also use a previously published Bayesian algorithm. We find that the Bayesian approaches provide optimal results and the timing observations of the pulsar place a 95 per cent upper limit on the sky-averaged strain amplitude of CGWs to be $\lesssim3.5\times10^{-13}$ at a reference frequency of 1 $μ$Hz. We also find a 95 per cent upper limit on the sky-averaged strain amplitude of low-frequency CGWs to be $\lesssim1.4\times10^{-14}$ at a reference frequency of 20~nHz. △ Less

Submitted 27 April, 2018; originally announced April 2018.

Comments: 12 pages, 4 figures. Accepted for publication in MNRAS

A direct measurement of sense of rotation of PSR J0737$-$3039A

Authors: Nihan Pol, Maura McLaughlin, Michael Kramer, Ingrid Stairs, Benetge B. P. Perera, Andrea Possenti

Abstract: We apply the algorithm published by Liang et al. (2014) to describe the Double Pulsar system J0737$-$3039 and extract the sense of rotation of first born recycled pulsar PSR J0737$-$3039A. We find that this pulsar is rotating prograde in its orbit. This is the first direct measurement of the sense of rotation of a pulsar with respect to its orbit and a direct confirmation of the rotating lighthous… ▽ More We apply the algorithm published by Liang et al. (2014) to describe the Double Pulsar system J0737$-$3039 and extract the sense of rotation of first born recycled pulsar PSR J0737$-$3039A. We find that this pulsar is rotating prograde in its orbit. This is the first direct measurement of the sense of rotation of a pulsar with respect to its orbit and a direct confirmation of the rotating lighthouse model for pulsars. This result confirms that the spin angular momentum vector is closely aligned with the orbital angular momentum, suggesting that kick of the supernova producing the second born pulsar J0737$-$3039B was small. △ Less

Submitted 12 December, 2017; originally announced December 2017.

Comments: Accepted for publication in ApJ

Evidence for an intermediate-mass black hole in the globular cluster NGC 6624

Authors: B. B. P. Perera, B. W. Stappers, A. G. Lyne, C. G. Bassa, I. Cognard, L. Guillemot, M. Kramer, G. Theureau, G. Desvignes

Abstract: PSR B1820$-$30A is located in the globular cluster NGC 6624 and is the closest known pulsar to the centre of any globular cluster. We present more than 25 years of high-precision timing observations of this millisecond pulsar and obtain four rotational frequency time derivative measurements. Modelling these higher-order derivatives as being due to orbital motion, we find solutions which indicate t… ▽ More PSR B1820$-$30A is located in the globular cluster NGC 6624 and is the closest known pulsar to the centre of any globular cluster. We present more than 25 years of high-precision timing observations of this millisecond pulsar and obtain four rotational frequency time derivative measurements. Modelling these higher-order derivatives as being due to orbital motion, we find solutions which indicate that the pulsar is in either a low-eccentricity ($0.33\lesssim e\lesssim0.4$) smaller orbit with a low mass companion (such as a main sequence star, white dwarf, neutron star, or stellar mass black hole) or a high-eccentricity ($e\gtrsim0.9$) larger orbit with a massive companion. The cluster mass properties and the observed properties of 4U 1820$-$30 and the other pulsars in the cluster argue against the low-eccentricity possibility. The high-eccentricity solution reveals that the pulsar is most likely orbiting around an intermediate-mass black hole (IMBH) of mass $> 7,500$~M$_\odot$ located at the cluster centre. A gravitational model for the globular cluster, which includes such a central black hole (BH), predicts an acceleration that is commensurate with that measured for the pulsar. It further predicts that the model-dependent minimum mass of the IMBH is $\sim60,000$~M$_\odot$. Accounting for the associated contribution to the observed period derivative indicates that the $γ$-ray efficiency of the pulsar should be between 0.08 and 0.2. Our results suggest that other globular clusters may also contain central black holes and they may be revealed by the study of new pulsars found sufficiently close to their centres. Note that we found an erratum in Section 5 and thus, the $\sim$60,000~M$_\odot$ mass mentioned above has to be replaced by the correct model-dependent mass limit of $\sim$20,000~M$_\odot$. See the erratum appended. △ Less

Submitted 17 May, 2017; v1 submitted 2 May, 2017; originally announced May 2017.

Comments: 15 pages, 10 figures, Accepted by MNRAS on 23 February 2017. Erratum was accepted by MNRAS on 17 May 2017

Journal ref: MNRAS 468, 2114-2127 (2017)

From Spin Noise to Systematics: Stochastic Processes in the First International Pulsar Timing Array Data Release

Authors: L. Lentati, R. M. Shannon, W. A. Coles, J. P. W. Verbiest, R. van Haasteren, J. A. Ellis, R. N. Caballero, R. N. Manchester, Z. Arzoumanian, S. Babak, C. G. Bassa, N. D. R. Bhat, P. Brem, M. Burgay, S. Burke-Spolaor, D. Champion, S. Chatterjee, I. Cognard, J. M. Cordes, S. Dai, P. Demorest, G. Desvignes, T. Dolch, R. D. Ferdman, E. Fonseca , et al. (58 additional authors not shown)

Abstract: We analyse the stochastic properties of the 49 pulsars that comprise the first International Pulsar Timing Array (IPTA) data release. We use Bayesian methodology, performing model selection to determine the optimal description of the stochastic signals present in each pulsar. In addition to spin-noise and dispersion-measure (DM) variations, these models can include timing noise unique to a single… ▽ More We analyse the stochastic properties of the 49 pulsars that comprise the first International Pulsar Timing Array (IPTA) data release. We use Bayesian methodology, performing model selection to determine the optimal description of the stochastic signals present in each pulsar. In addition to spin-noise and dispersion-measure (DM) variations, these models can include timing noise unique to a single observing system, or frequency band. We show the improved radio-frequency coverage and presence of overlapping data from different observing systems in the IPTA data set enables us to separate both system and band-dependent effects with much greater efficacy than in the individual PTA data sets. For example, we show that PSR J1643$-$1224 has, in addition to DM variations, significant band-dependent noise that is coherent between PTAs which we interpret as coming from time-variable scattering or refraction in the ionised interstellar medium. Failing to model these different contributions appropriately can dramatically alter the astrophysical interpretation of the stochastic signals observed in the residuals. In some cases, the spectral exponent of the spin noise signal can vary from 1.6 to 4 depending upon the model, which has direct implications for the long-term sensitivity of the pulsar to a stochastic gravitational-wave (GW) background. By using a more appropriate model, however, we can greatly improve a pulsar's sensitivity to GWs. For example, including system and band-dependent signals in the PSR J0437$-$4715 data set improves the upper limit on a fiducial GW background by $\sim 60\%$ compared to a model that includes DM variations and spin-noise only. △ Less

Submitted 16 February, 2016; originally announced February 2016.

Comments: 29 pages. 16 figures. Accepted for publication in MNRAS

The International Pulsar Timing Array: First Data Release

Authors: J. P. W. Verbiest, L. Lentati, G. Hobbs, R. van Haasteren, P. B. Demorest, G. H. Janssen, J. -B. Wang, G. Desvignes, R. N. Caballero, M. J. Keith, D. J. Champion, Z. Arzoumanian, S. Babak, C. G. Bassa, N. D. R. Bhat, A. Brazier, P. Brem, M. Burgay, S. Burke-Spolaor, S. J. Chamberlin, S. Chatterjee, B. Christy, I. Cognard, J. M. Cordes, S. Dai , et al. (67 additional authors not shown)

Abstract: The highly stable spin of neutron stars can be exploited for a variety of (astro-)physical investigations. In particular arrays of pulsars with rotational periods of the order of milliseconds can be used to detect correlated signals such as those caused by gravitational waves. Three such "Pulsar Timing Arrays" (PTAs) have been set up around the world over the past decades and collectively form the… ▽ More The highly stable spin of neutron stars can be exploited for a variety of (astro-)physical investigations. In particular arrays of pulsars with rotational periods of the order of milliseconds can be used to detect correlated signals such as those caused by gravitational waves. Three such "Pulsar Timing Arrays" (PTAs) have been set up around the world over the past decades and collectively form the "International" PTA (IPTA). In this paper, we describe the first joint analysis of the data from the three regional PTAs, i.e. of the first IPTA data set. We describe the available PTA data, the approach presently followed for its combination and suggest improvements for future PTA research. Particular attention is paid to subtle details (such as underestimation of measurement uncertainty and long-period noise) that have often been ignored but which become important in this unprecedentedly large and inhomogeneous data set. We identify and describe in detail several factors that complicate IPTA research and provide recommendations for future pulsar timing efforts. The first IPTA data release presented here (and available online) is used to demonstrate the IPTA's potential of improving upon gravitational-wave limits placed by individual PTAs by a factor of ~2 and provides a 2-sigma limit on the dimensionless amplitude of a stochastic GWB of 1.7x10^{-15} at a frequency of 1 yr^{-1}. This is 1.7 times less constraining than the limit placed by (Shannon et al. 2015), due mostly to the more recent, high-quality data they used. △ Less

Submitted 11 February, 2016; originally announced February 2016.

Comments: 25 pages, 6 tables, 5 figures. Accepted for publication in MNRAS

Correlated spin-down rates and radio emission in PSR B1859+07

Authors: B. B. P. Perera, B. W. Stappers, P. Weltevrede, A. G. Lyne, J. M. Rankin

Abstract: We study the spin-down changes of PSR B1859$+$07 over a period of more than 28 years of radio observation. We identify that the time derivative of the rotational frequency ($ν$) varies quasi-periodically with a period of $\sim$350 days, switching mainly between two spin-down states. The profile shape of the pulsar is correlated with the $\dot ν$ variation, producing two slightly different profile… ▽ More We study the spin-down changes of PSR B1859$+$07 over a period of more than 28 years of radio observation. We identify that the time derivative of the rotational frequency ($ν$) varies quasi-periodically with a period of $\sim$350 days, switching mainly between two spin-down states. The profile shape of the pulsar is correlated with the $\dot ν$ variation, producing two slightly different profile shapes corresponding to high- and low-$\dot ν$ states. In addition to these two normal emission states, we confirm the occasional flare-state of the pulsar, in which the emission appears early in spin phase compared to that of the common normal emission. The profile of the flare-state is significantly different from that of the two normal emission states. The correlation analysis further shows that the flare-state is not directly linked with the $\dot ν$ changes. With a simple emission beam model, we estimate the emission altitude of the normal emission to be 240~km, and explain the origin of the flare-state as an emission height variation from the leading edge of the beam. We also argue that the emission of these states can be explained with a partially active beam model. In this scenario, the trailing portion of the radio beam is usually active and the normal emission is produced. The flare-state occurs when the leading edge of the beam becomes active while the trailing part is being blocked. This model estimates a fixed emission altitude of 360~km. However, the cause of the flare-state (i.e. the emission height variation, or the time-dependent activity across the radio beam) is not easily explained. △ Less

Submitted 15 October, 2015; originally announced October 2015.

Comments: 9 pages, 8 figures, Accepted by MNRAS on 14 October 2015

Coordinated X-ray, Ultraviolet, Optical, and Radio Observations of the PSR J1023+0038 System in a Low-mass X-ray Binary State

Authors: Slavko Bogdanov, Anne M. Archibald, Cees Bassa, Adam T. Deller, Jules P. Halpern, George Heald, Jason W. T. Hessels, Gemma H. Janssen, Andrew G. Lyne, Javier Moldon, Zsolt Paragi, Alessandro Patruno, Benetge B. P. Perera, Ben W. Stappers, Shriharsh P. Tendulkar, Caroline R. D'Angelo, Rudy Wijnands

Abstract: The PSR J1023+0038 binary system hosts a neutron star and a low-mass, main-sequence-like star. It switches on year timescales between states as an eclipsing radio millisecond pulsar and a low-mass X-ray binary. We present a multi-wavelength observational campaign of PSR J1023+0038 in its most recent low-mass X-ray binary state. Two long XMM-Newton observations reveal that the system spends ~70% of… ▽ More The PSR J1023+0038 binary system hosts a neutron star and a low-mass, main-sequence-like star. It switches on year timescales between states as an eclipsing radio millisecond pulsar and a low-mass X-ray binary. We present a multi-wavelength observational campaign of PSR J1023+0038 in its most recent low-mass X-ray binary state. Two long XMM-Newton observations reveal that the system spends ~70% of the time in a $\approx$$3\times10^{33}$ erg/s X-ray luminosity mode, which, as shown in Archibald et al. (2014), exhibits coherent X-ray pulsations. This emission is interspersed with frequent lower flux mode intervals with $\approx$$5\times 10^{32}$ erg/s and sporadic flares reaching up to $\approx$$10^{34}$ erg/s, with neither mode showing significant X-ray pulsations. The switches between the three flux modes occur on timescales of order 10 s. In the UV and optical, we observe occasional intense flares coincident with those observed in X-rays. Our radio timing observations reveal no pulsations at the pulsar period during any of the three X-ray modes, presumably due to complete quenching of the radio emission mechanism by the accretion flow. Radio imaging detects highly variable, flat-spectrum continuum emission from PSR J1023+0038, consistent with an origin in a weak jet-like outflow. Our concurrent X-ray and radio continuum data sets do not exhibit any correlated behavior. The observational evidence we present bears qualitative resemblance to the behavior predicted by some existing propeller and trapped disk accretion models although none can fully account for all aspects of the rich phenomenology of this system. △ Less

Submitted 17 April, 2015; v1 submitted 16 December, 2014; originally announced December 2014.

Comments: 24 pages, 15 figures; accepted for publication in the Astrophysical Journal

Understanding the spin-down rate changes of PSR B0919+06

Authors: B. B. P. Perera, B. W. Stappers, P. Weltevrede, A. G. Lyne, C. G. Bassa

Abstract: We study the spin-down properties of PSR B0919+06 based on almost 30 years of radio observations. We confirm that the time derivative of the rotational frequency $\dot ν$ is modulated quasi-periodically and show that it exhibits a repeating double-peaked structure throughout the entire observation span. We model the $\dot ν$ variation of the pulsar assuming two spin-down rates with sudden switches… ▽ More We study the spin-down properties of PSR B0919+06 based on almost 30 years of radio observations. We confirm that the time derivative of the rotational frequency $\dot ν$ is modulated quasi-periodically and show that it exhibits a repeating double-peaked structure throughout the entire observation span. We model the $\dot ν$ variation of the pulsar assuming two spin-down rates with sudden switches between them in time. Our results show that the double-peak structure in $\dot ν$ has a repetition time of about 630 days until MJD 52000 (April 2001) and 550 days since then. During this cycle, the pulsar spin varies from the lower spin-down rate to the upper spin-down rate twice with different amounts of time spent in each state, resulting in a further quasi-stable secondary modulation of the two-state switching. This particular spin-down state switching is broadly consistent with free precession of the pulsar, however, a strong evidence linked with this mechanism is not clearly established. We also confirm that the pulsar occasionally emits groups of pulses which appear early in pulse phase, so-called "flares", and these events significantly contribute to the pulse profile shape. We find the $\dot ν$ modulation and the pulse shape variations are correlated throughout the observations. However, the flare-state is not entirely responsible for this correlation. In addition to the flare-state, we detect flare-like events from the pulsar in single pulse observations. During these events, the shift in pulse phase is small compared to that of the main flare-state and clearly visible only in single pulse observations. △ Less

Submitted 20 October, 2014; originally announced October 2014.

Comments: 10 pages, 12 Figures, Accepted by MNRAS on 10 October 2014