-
The NANOGrav 15 yr Data Set: Running of the Spectral Index
Authors:
Gabriella Agazie,
Akash Anumarlapudi,
Anne M. Archibald,
Zaven Arzoumanian,
Jeremy George Baier,
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,
Heling Deng,
Lankeswar Dey,
Timothy Dolch
, et al. (80 additional authors not shown)
Abstract:
The NANOGrav 15-year data provides compelling evidence for a stochastic gravitational-wave (GW) background at nanohertz frequencies. The simplest model-independent approach to characterizing the frequency spectrum of this signal consists in a simple power-law fit involving two parameters: an amplitude A and a spectral index γ. In this paper, we consider the next logical step beyond this minimal sp…
▽ More
The NANOGrav 15-year data provides compelling evidence for a stochastic gravitational-wave (GW) background at nanohertz frequencies. The simplest model-independent approach to characterizing the frequency spectrum of this signal consists in a simple power-law fit involving two parameters: an amplitude A and a spectral index γ. In this paper, we consider the next logical step beyond this minimal spectral model, allowing for a running (i.e., logarithmic frequency dependence) of the spectral index, γ_run(f) = γ+ β\ln(f/f_ref). We fit this running-power-law (RPL) model to the NANOGrav 15-year data and perform a Bayesian model comparison with the minimal constant-power-law (CPL) model, which results in a 95% credible interval for the parameter βconsistent with no running, β\in [-0.80,2.96], and an inconclusive Bayes factor, B(RPL vs. CPL) = 0.69 +- 0.01. We thus conclude that, at present, the minimal CPL model still suffices to adequately describe the NANOGrav signal; however, future data sets may well lead to a measurement of nonzero β. Finally, we interpret the RPL model as a description of primordial GWs generated during cosmic inflation, which allows us to combine our results with upper limits from big-bang nucleosynthesis, the cosmic microwave background, and LIGO-Virgo-KAGRA.
△ Less
Submitted 19 August, 2024;
originally announced August 2024.
-
The NANOGrav 15 yr data set: Posterior predictive checks for gravitational-wave detection with pulsar timing arrays
Authors:
Gabriella Agazie,
Akash Anumarlapudi,
Anne M. Archibald,
Zaven Arzoumanian,
Jeremy George Baier,
Paul T. Baker,
Bence Bécsy,
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,
Neil J. Cornish,
Fronefield Crawford,
H. Thankful Cromartie,
Kathryn Crowter,
Megan E. DeCesar,
Paul B. Demorest,
Heling Deng,
Lankeswar Dey
, et al. (77 additional authors not shown)
Abstract:
Pulsar-timing-array experiments have reported evidence for a stochastic background of nanohertz gravitational waves consistent with the signal expected from a population of supermassive--black-hole binaries. Those analyses assume power-law spectra for intrinsic pulsar noise and for the background, as well as a Hellings--Downs cross-correlation pattern among the gravitational-wave--induced residual…
▽ More
Pulsar-timing-array experiments have reported evidence for a stochastic background of nanohertz gravitational waves consistent with the signal expected from a population of supermassive--black-hole binaries. Those analyses assume power-law spectra for intrinsic pulsar noise and for the background, as well as a Hellings--Downs cross-correlation pattern among the gravitational-wave--induced residuals across pulsars. These assumptions are idealizations that may not be realized in actuality. We test them in the NANOGrav 15 yr data set using Bayesian posterior predictive checks: after fitting our fiducial model to real data, we generate a population of simulated data-set replications, and use them to assess whether the optimal-statistic significance, inter-pulsar correlations, and spectral coefficients assume extreme values for the real data when compared to the replications. We confirm that the NANOGrav 15 yr data set is consistent with power-law and Hellings--Downs assumptions. We also evaluate the evidence for the stochastic background using posterior-predictive versions of the frequentist optimal statistic and of Bayesian model comparison, and find comparable significance (3.2\ $σ$ and 3\ $σ$ respectively) to what was previously reported for the standard statistics. We conclude with novel visualizations of the reconstructed gravitational waveforms that enter the residuals for each pulsar. Our analysis strengthens confidence in the identification and characterization of the gravitational-wave background as reported by NANOGrav.
△ Less
Submitted 29 July, 2024;
originally announced July 2024.
-
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 23 August, 2024; v1 submitted 8 July, 2024;
originally announced July 2024.
-
Constraining the selection corrected luminosity function and total pulse count for radio transients
Authors:
Fengqiu Adam Dong,
Antonio Herrera-Martin,
Ingrid Stairs,
Radu V. Craiu,
Kathryn Crowter,
Gwendolyn M. Eadie,
Emmanuel Fonseca,
Deborah Good,
James W. Mckee,
Bradley W. Meyers,
Aaron B. Pearlman,
David C. Stenning
Abstract:
Studying transient phenomena, such as individual pulses from pulsars, has garnered considerable attention in the era of astronomical big data. Of specific interest to this study are Rotating Radio Transients (RRATs), nulling, and intermittent pulsars. This study introduces a new algorithm named LuNfit, tailored to correct the selection biases originating from the telescope and detection pipelines.…
▽ More
Studying transient phenomena, such as individual pulses from pulsars, has garnered considerable attention in the era of astronomical big data. Of specific interest to this study are Rotating Radio Transients (RRATs), nulling, and intermittent pulsars. This study introduces a new algorithm named LuNfit, tailored to correct the selection biases originating from the telescope and detection pipelines. Ultimately LuNfit estimates the intrinsic luminosity distribution and nulling fraction of the single pulses emitted by pulsars. LuNfit relies on Bayesian nested sampling so that the parameter space can be fully explored. Bayesian nested sampling also provides the additional benefit of simplifying model comparisons through the Bayes ratio. The robustness of LuNfit is shown through simulations and applying LuNfit onto pulsars with known nulling fractions. LuNfit is then applied to three RRATs, J0012+5431, J1538+1523, and J2355+1523, extracting their intrinsic luminosity distribution and burst rates. We find that their nulling fraction is 0.4(2), 0.749(5) and 0.995(2) respectively. We further find that a log-normal distribution likely describes the single pulse luminosity distribution of J0012+5431 and J1538+1523, while the Bayes ratio for J2355+1523 slightly favors an exponential distribution. We show the conventional method of correcting selection effects by "scaling up" the missed fraction of radio transients can be unreliable when the mean luminosity of the source is faint relative to the telescope sensitivity. Finally, we discuss the limitations of the current implementation of LuNfit while also delving into potential enhancements that would enable LuNfit to be applied to sources with complex pulse morphologies.
△ Less
Submitted 6 June, 2024;
originally announced June 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.
-
High-cadence Timing of Binary Pulsars with CHIME
Authors:
Chia Min Tan,
Emmanuel Fonseca,
Kathryn Crowter,
Fengqiu Adam Dong,
Victoria M. Kaspi,
Kiyoshi W. Masui,
James W. McKee,
Bradley W. Meyers,
Scott M. Ransom,
Ingrid H. Stairs
Abstract:
We performed near-daily observations on the binary pulsars PSR J0218+4232, PSR J1518+4904 and PSR J2023+2853 with the Canadian Hydrogen Intensity Mapping Experiment (CHIME). For the first time, we detected the Shapiro time delay in all three pulsar-binary systems, using only 2--4 years of CHIME/Pulsar timing data. We measured the pulsar masses to be $1.49^{+0.23}_{-0.20}$ M$_\odot$,…
▽ More
We performed near-daily observations on the binary pulsars PSR J0218+4232, PSR J1518+4904 and PSR J2023+2853 with the Canadian Hydrogen Intensity Mapping Experiment (CHIME). For the first time, we detected the Shapiro time delay in all three pulsar-binary systems, using only 2--4 years of CHIME/Pulsar timing data. We measured the pulsar masses to be $1.49^{+0.23}_{-0.20}$ M$_\odot$, $1.470^{+0.030}_{-0.034}$ M$_\odot$ and $1.50^{+0.49}_{-0.38}$ M$_\odot$ respectively. The companion mass to PSR J0218+4232 was found to be $0.179^{+0.018}_{-0.016}$ M$_\odot$. We constrained the mass of the neutron-star companion of PSR J1518+4904 to be $1.248^{+0.035}_{-0.029}$ M$_\odot$, using the observed apsidal motion as a constraint on mass estimation. The binary companion to PSR J2023+2853 was found to have a mass of $0.93^{+0.17}_{-0.14}$ M$_\odot$; in the context of the near-circular orbit, this mass estimate suggests that the companion to PSR J2023+2853 is likely a high-mass white dwarf. By comparing the timing model obtained for PSR J0218+4232 with previous observations, we found a significant change in the observed orbital period of the system of $\dot{P_{\rm b}} = 0.14(2) \times 10^{-12}$; we determined that this variation arises from ``Shklovskii acceleration" due to relative motion of the binary system, and used this measurement to estimate a distance of $d=(6.7 \pm 1.0)$ kpc to PSR J0218+4232. This work demonstrates the capability of high-cadence observations, enabled by the CHIME/Pulsar system, to detect and refine general-relativistic effects of binary pulsars over short observing timescales.
△ Less
Submitted 12 February, 2024;
originally announced February 2024.
-
The Green Bank North Celestial Cap Survey IX: Timing Follow-up for 128 Pulsars
Authors:
A. E. McEwen,
J. K. Swiggum,
D. L. Kaplan,
C. M. Tan,
B. W. Meyers,
E. Fonseca,
G. Y. Agazie,
P. Chawla,
K. Crowter,
M. E. DeCesar,
T. Dolch,
F. A. Dong,
W. Fiore,
E. Fonseca,
D. C. Good,
A. G. Istrate,
V. M. Kaspi,
V. I. Kondratiev,
J. van Leeuwen,
L. Levin,
E. F. Lewis,
R. S. Lynch,
K. W. Masui,
J. W. McKee,
M. A. McLaughlin
, et al. (6 additional authors not shown)
Abstract:
The Green Bank North Celestial Cap survey is one of the largest and most sensitive searches for pulsars and transient radio objects. Observations for the survey have finished; priorities have shifted toward long-term monitoring of its discoveries. In this study, we have developed a pipeline to handle large datasets of archival observations and connect them to recent, high-cadence observations take…
▽ More
The Green Bank North Celestial Cap survey is one of the largest and most sensitive searches for pulsars and transient radio objects. Observations for the survey have finished; priorities have shifted toward long-term monitoring of its discoveries. In this study, we have developed a pipeline to handle large datasets of archival observations and connect them to recent, high-cadence observations taken using the Canadian Hydrogen Intensity Mapping Experiment (CHIME) telescope. This pipeline handles data for 128 pulsars and has produced measurements of spin, positional, and orbital parameters that connect data over observation gaps as large as 2000 days. We have also measured glitches in the timing residuals for five of the pulsars included and proper motion for 19 sources (13 new). We include updates to orbital parameters for 19 pulsars, including 9 previously unpublished binaries. For two of these binaries, we provide updated measurements of post-Keplerian binary parameters, which result in much more precise estimates of the total masses of both systems. For PSR J0509+3801, the much improved measurement of the Einstein delay yields much improved mass measurements for the pulsar and its companion, 1.399(6)\Msun and 1.412(6)\Msun, respectively. For this system, we have also obtained a measurement of the orbital decay due to the emission of gravitational waves: $\dot{P}_{\rm B} = -1.37(7)\times10^{-12}$, which is in agreement with the rate predicted by general relativity for these masses.
△ Less
Submitted 26 July, 2024; v1 submitted 12 December, 2023;
originally announced December 2023.
-
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.
-
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.
-
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.
-
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.
-
The NANOGrav 15-year Data Set: Bayesian Limits on Gravitational Waves from Individual Supermassive Black Hole Binaries
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,
Robin Case,
J. Andrew Casey-Clyde,
Maria Charisi,
Shami Chatterjee,
Tyler Cohen,
James M. Cordes,
Neil Cornish,
Fronefield Crawford,
H. Thankful Cromartie,
Kathryn Crowter,
Megan DeCesar,
Paul B. Demorest,
Matthew C. Digman,
Timothy Dolch,
Brendan Drachler
, et al. (74 additional authors not shown)
Abstract:
Evidence for a low-frequency stochastic gravitational wave background has recently been reported based on analyses of pulsar timing array data. The most likely source of such a background is a population of supermassive black hole binaries, the loudest of which may be individually detected in these datasets. Here we present the search for individual supermassive black hole binaries in the NANOGrav…
▽ More
Evidence for a low-frequency stochastic gravitational wave background has recently been reported based on analyses of pulsar timing array data. The most likely source of such a background is a population of supermassive black hole binaries, the loudest of which may be individually detected in these datasets. Here we present the search for individual supermassive black hole binaries in the NANOGrav 15-year dataset. We introduce several new techniques, which enhance the efficiency and modeling accuracy of the analysis. The search uncovered weak evidence for two candidate signals, one with a gravitational-wave frequency of $\sim$4 nHz, and another at $\sim$170 nHz. The significance of the low-frequency candidate was greatly diminished when Hellings-Downs correlations were included in the background model. The high-frequency candidate was discounted due to the lack of a plausible host galaxy, the unlikely astrophysical prior odds of finding such a source, and since most of its support comes from a single pulsar with a commensurate binary period. Finding no compelling evidence for signals from individual binary systems, we place upper limits on the strain amplitude of gravitational waves emitted by such systems.
△ Less
Submitted 28 June, 2023;
originally announced June 2023.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
Searching for continuous Gravitational Waves in the second data release of the International Pulsar Timing Array
Authors:
M. Falxa,
S. Babak,
P. T. Baker,
B. Bécsy,
A. Chalumeau,
S. Chen,
Z. Chen,
N. J. Cornish,
L. Guillemot,
J. S. Hazboun,
C. M. F. Mingarelli,
A. Parthasarathy,
A. Petiteau,
N. S. Pol,
A. Sesana,
S. B. Spolaor,
S. R. Taylor,
G. Theureau,
M. Vallisneri,
S. J. Vigeland,
C. A. Witt,
X. Zhu,
J. Antoniadis,
Z. Arzoumanian,
M. Bailes
, et al. (102 additional authors not shown)
Abstract:
The International Pulsar Timing Array 2nd data release is the combination of datasets from worldwide collaborations. In this study, we search for continuous waves: gravitational wave signals produced by individual supermassive black hole binaries in the local universe. We consider binaries on circular orbits and neglect the evolution of orbital frequency over the observational span. We find no evi…
▽ More
The International Pulsar Timing Array 2nd data release is the combination of datasets from worldwide collaborations. In this study, we search for continuous waves: gravitational wave signals produced by individual supermassive black hole binaries in the local universe. We consider binaries on circular orbits and neglect the evolution of orbital frequency over the observational span. We find no evidence for such signals and set sky averaged 95% upper limits on their amplitude h 95 . The most sensitive frequency is 10nHz with h 95 = 9.1 10-15 . We achieved the best upper limit to date at low and high frequencies of the PTA band thanks to improved effective cadence of observations. In our analysis, we have taken into account the recently discovered common red noise process, which has an impact at low frequencies. We also find that the peculiar noise features present in some pulsars data must be taken into account to reduce the false alarm. We show that using custom noise models is essential in searching for continuous gravitational wave signals and setting the upper limit.
△ Less
Submitted 19 March, 2023;
originally announced March 2023.
-
The second set of pulsar discoveries by CHIME/FRB/Pulsar: 14 Rotating Radio Transients and 7 pulsars
Authors:
Fengqiu Adam Dong,
Kathryn Crowter,
Bradley W. Meyers,
Ziggy Pleunis,
Ingrid Stairs,
Chia Min Tan,
Tinyau Timothy Yu,
Patrick J. Boyle,
Amanda M. Cook,
Emmanuel Fonseca,
B. M. Gaensler,
Deborah C. Good,
Victoria Kaspi,
James W. McKee,
Chitrang Patel,
Aaron B. Pearlman
Abstract:
The Canadian Hydrogen Mapping Experiment (CHIME) is a radio telescope located in British Columbia, Canada. The large FOV allows CHIME/FRB to be an exceptional pulsar and Rotating Radio Transient (RRAT) finding machine, despite saving only the metadata of incoming Galactic events. We have developed a pipeline to search for pulsar/RRAT candidates using DBSCAN, a clustering algorithm. Follow-up obser…
▽ More
The Canadian Hydrogen Mapping Experiment (CHIME) is a radio telescope located in British Columbia, Canada. The large FOV allows CHIME/FRB to be an exceptional pulsar and Rotating Radio Transient (RRAT) finding machine, despite saving only the metadata of incoming Galactic events. We have developed a pipeline to search for pulsar/RRAT candidates using DBSCAN, a clustering algorithm. Follow-up observations are then scheduled with the more sensitive CHIME/Pulsar instrument capable of near-daily high time resolution spectra observations. We have developed the CHIME/Pulsar Single Pulse Pipeline to automate the processing of CHIME/Pulsar search-mode data. We report the discovery of 21 new Galactic sources, with 14 RRATs, 6 isolated long-period pulsars and 1 binary system. Owing to CHIME/Pulsar's observations we have obtained timing solutions for 8 of the 14 RRATs along with all the regular pulsars and the binary system. Notably we report that the binary system is in a long orbit of 412 days with a minimum companion mass of 0.1303 solar masses and no evidence of an optical companion within 10" of the pulsar position. This highlights that working synergistically with CHIME/FRB's large survey volume CHIME/Pulsar can obtain arc second localisations for low burst rate RRATs though pulsar timing. We find that the properties of our newly discovered RRATs are consistent with those of the presently known population. They tend to have lower burst rates than those found in previous surveys, which is likely due to survey bias rather than the underlying population.
△ Less
Submitted 27 July, 2023; v1 submitted 17 October, 2022;
originally announced October 2022.
-
CHIME Discovery of a Binary Pulsar with a Massive Non-Degenerate Companion
Authors:
Bridget C. Andersen,
Emmanuel Fonseca,
J. W. McKee,
B. W. Meyers,
Jing Luo,
C. M. Tan,
I. H. Stairs,
Victoria M. Kaspi,
M. H. van Kerkwijk,
Mohit Bhardwaj,
P. J. Boyle,
Kathryn Crowter,
Paul B. Demorest,
Fengqui A. Dong,
Deborah C. Good,
Jane F. Kaczmarek,
Calvin Leung,
Kiyoshi W. Masui,
Arun Naidu,
Cherry Ng,
Chitrang Patel,
Aaron B. Pearlman,
Ziggy Pleunis,
Masoud Rafiei-Ravandi,
Mubdi Rahman
, et al. (3 additional authors not shown)
Abstract:
Of the more than $3{,}000$ radio pulsars currently known, only ${\sim}300$ are in binary systems, and only five of these consist of young pulsars with massive non-degenerate companions. We present the discovery and initial timing, accomplished using the Canadian Hydrogen Intensity Mapping Experiment telescope (CHIME), of the sixth such binary pulsar, PSR J2108+4516, a $0.577$-s radio pulsar in a 2…
▽ More
Of the more than $3{,}000$ radio pulsars currently known, only ${\sim}300$ are in binary systems, and only five of these consist of young pulsars with massive non-degenerate companions. We present the discovery and initial timing, accomplished using the Canadian Hydrogen Intensity Mapping Experiment telescope (CHIME), of the sixth such binary pulsar, PSR J2108+4516, a $0.577$-s radio pulsar in a 269-day orbit of eccentricity 0.09 with a companion of minimum mass $11$ M$_{\odot}$. Notably, the pulsar undergoes periods of substantial eclipse, disappearing from the CHIME $400{-}800$ MHz observing band for a large fraction of its orbit, and displays significant dispersion measure and scattering variations throughout its orbit, pointing to the possibility of a circumstellar disk or very dense stellar wind associated with the companion star. Subarcsecond resolution imaging with the Karl G. Jansky Very Large Array unambiguously demonstrates that the companion is a bright, $V \simeq 11$ OBe star, EM* UHA 138, located at a distance of $3.26(14)$ kpc. Archival optical observations of \companion{} approximately suggest a companion mass ranging from $17.5$ M$_{\odot} < M_{\rm c} < 23$ M$_{\odot}$, in turn constraining the orbital inclination angle to $50.3^{\circ} \lesssim i \lesssim 58.3^{\circ}$. With further multi-wavelength followup, PSR J2108+4516 promises to serve as another rare laboratory for the exploration of companion winds, circumstellar disks, and short-term evolution through extended-body orbital dynamics.
△ Less
Submitted 30 January, 2023; v1 submitted 14 September, 2022;
originally announced September 2022.
-
A Gamma-ray Pulsar Timing Array Constrains the Nanohertz Gravitational Wave Background
Authors:
M. Ajello,
W. B. Atwood,
L. Baldini,
J. Ballet,
G. Barbiellini,
D. Bastieri,
R. Bellazzini,
A. Berretta,
B. Bhattacharyya,
E. Bissaldi,
R. D. Blandford,
E. Bloom,
R. Bonino,
P. Bruel,
R. Buehler,
E. Burns,
S. Buson,
R. A. Cameron,
P. A. Caraveo,
E. Cavazzuti,
N. Cibrario,
S. Ciprini,
C. J. Clark,
I. Cognard,
J. Coronado-Blázquez
, et al. (107 additional authors not shown)
Abstract:
After large galaxies merge, their central supermassive black holes are expected to form binary systems whose orbital motion generates a gravitational wave background (GWB) at nanohertz frequencies. Searches for this background utilize pulsar timing arrays, which perform long-term monitoring of millisecond pulsars (MSPs) at radio wavelengths. We use 12.5 years of Fermi Large Area Telescope data to…
▽ More
After large galaxies merge, their central supermassive black holes are expected to form binary systems whose orbital motion generates a gravitational wave background (GWB) at nanohertz frequencies. Searches for this background utilize pulsar timing arrays, which perform long-term monitoring of millisecond pulsars (MSPs) at radio wavelengths. We use 12.5 years of Fermi Large Area Telescope data to form a gamma-ray pulsar timing array. Results from 35 bright gamma-ray pulsars place a 95\% credible limit on the GWB characteristic strain of $1.0\times10^{-14}$ at 1 yr$^{-1}$, which scales as the observing time span $t_{\mathrm{obs}}^{-13/6}$. This direct measurement provides an independent probe of the GWB while offering a check on radio noise models.
△ Less
Submitted 11 April, 2022;
originally announced April 2022.
-
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.
-
Narrowband searches for continuous and long-duration transient gravitational waves from known pulsars in the LIGO-Virgo third observing run
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
the KAGRA Collaboration,
R. Abbott,
T. D. Abbott,
F. Acernese,
K. Ackley,
C. Adams,
N. Adhikari,
R. X. Adhikari,
V. B. Adya,
C. Affeldt,
D. Agarwal,
M. Agathos,
K. Agatsuma,
N. Aggarwal,
O. D. Aguiar,
L. Aiello,
A. Ain,
P. Ajith,
T. Akutsu,
S. Albanesi,
A. Allocca,
P. A. Altin,
A. Amato
, et al. (1636 additional authors not shown)
Abstract:
Isolated neutron stars that are asymmetric with respect to their spin axis are possible sources of detectable continuous gravitational waves. This paper presents a fully-coherent search for such signals from eighteen pulsars in data from LIGO and Virgo's third observing run (O3). For known pulsars, efficient and sensitive matched-filter searches can be carried out if one assumes the gravitational…
▽ More
Isolated neutron stars that are asymmetric with respect to their spin axis are possible sources of detectable continuous gravitational waves. This paper presents a fully-coherent search for such signals from eighteen pulsars in data from LIGO and Virgo's third observing run (O3). For known pulsars, efficient and sensitive matched-filter searches can be carried out if one assumes the gravitational radiation is phase-locked to the electromagnetic emission. In the search presented here, we relax this assumption and allow the frequency and frequency time-derivative of the gravitational waves to vary in a small range around those inferred from electromagnetic observations. We find no evidence for continuous gravitational waves, and set upper limits on the strain amplitude for each target. These limits are more constraining for seven of the targets than the spin-down limit defined by ascribing all rotational energy loss to gravitational radiation. In an additional search we look in O3 data for long-duration (hours-months) transient gravitational waves in the aftermath of pulsar glitches for six targets with a total of nine glitches. We report two marginal outliers from this search, but find no clear evidence for such emission either. The resulting duration-dependent strain upper limits do not surpass indirect energy constraints for any of these targets.
△ Less
Submitted 27 June, 2022; v1 submitted 21 December, 2021;
originally announced December 2021.
-
Searches for Gravitational Waves from Known Pulsars at Two Harmonics in the Second and Third LIGO-Virgo Observing Runs
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
the KAGRA Collaboration,
R. Abbott,
H. Abe,
F. Acernese,
K. Ackley,
N. Adhikari,
R. X. Adhikari,
V. K. Adkins,
V. B. Adya,
C. Affeldt,
D. Agarwal,
M. Agathos,
K. Agatsuma,
N. Aggarwal,
O. D. Aguiar,
L. Aiello,
A. Ain,
P. Ajith,
T. Akutsu,
S. Albanesi,
R. A. Alfaidi,
A. Allocca,
P. A. Altin
, et al. (1672 additional authors not shown)
Abstract:
We present a targeted search for continuous gravitational waves (GWs) from 236 pulsars using data from the third observing run of LIGO and Virgo (O3) combined with data from the second observing run (O2). Searches were for emission from the $l=m=2$ mass quadrupole mode with a frequency at only twice the pulsar rotation frequency (single harmonic) and the $l=2, m=1,2$ modes with a frequency of both…
▽ More
We present a targeted search for continuous gravitational waves (GWs) from 236 pulsars using data from the third observing run of LIGO and Virgo (O3) combined with data from the second observing run (O2). Searches were for emission from the $l=m=2$ mass quadrupole mode with a frequency at only twice the pulsar rotation frequency (single harmonic) and the $l=2, m=1,2$ modes with a frequency of both once and twice the rotation frequency (dual harmonic). No evidence of GWs was found so we present 95\% credible upper limits on the strain amplitudes $h_0$ for the single harmonic search along with limits on the pulsars' mass quadrupole moments $Q_{22}$ and ellipticities $\varepsilon$. Of the pulsars studied, 23 have strain amplitudes that are lower than the limits calculated from their electromagnetically measured spin-down rates. These pulsars include the millisecond pulsars J0437\textminus4715 and J0711\textminus6830 which have spin-down ratios of 0.87 and 0.57 respectively. For nine pulsars, their spin-down limits have been surpassed for the first time. For the Crab and Vela pulsars our limits are factors of $\sim 100$ and $\sim 20$ more constraining than their spin-down limits, respectively. For the dual harmonic searches, new limits are placed on the strain amplitudes $C_{21}$ and $C_{22}$. For 23 pulsars we also present limits on the emission amplitude assuming dipole radiation as predicted by Brans-Dicke theory.
△ Less
Submitted 20 July, 2022; v1 submitted 25 November, 2021;
originally announced November 2021.
-
Bayesian Solar Wind Modeling with Pulsar Timing Arrays
Authors:
Jeffrey S. Hazboun,
Joseph Simon,
Dustin R. Madison,
Zaven Arzoumanian,
Kathryn Crowter,
Megan E. DeCesar,
Paul B. Demorest,
Timothy Dolch,
Justin A. Ellis,
Robert D. Ferdman,
Elizabeth C. Ferrara,
Emmanuel Fonseca,
Peter A. Gentile,
Glenn Jones,
Megan L. Jones,
Michael T. Lam,
Lina Levin,
Duncan R. Lorimer,
Ryan S. Lynch,
Maura A. McLaughlin,
Cherry Ng,
David J. Nice,
Timothy T. Pennucci,
Scott M. Ransom,
Paul S. Ray
, et al. (5 additional authors not shown)
Abstract:
Using Bayesian analyses we study the solar electron density with the NANOGrav 11-year pulsar timing array (PTA) dataset. Our model of the solar wind is incorporated into a global fit starting from pulse times-of-arrival. We introduce new tools developed for this global fit, including analytic expressions for solar electron column densities and open source models for the solar wind that port into e…
▽ More
Using Bayesian analyses we study the solar electron density with the NANOGrav 11-year pulsar timing array (PTA) dataset. Our model of the solar wind is incorporated into a global fit starting from pulse times-of-arrival. We introduce new tools developed for this global fit, including analytic expressions for solar electron column densities and open source models for the solar wind that port into existing PTA software. We perform an ab initio recovery of various solar wind model parameters. We then demonstrate the richness of information about the solar electron density, $n_E$, that can be gleaned from PTA data, including higher order corrections to the simple $1/r^2$ model associated with a free-streaming wind (which are informative probes of coronal acceleration physics), quarterly binned measurements of $n_E$ and a continuous time-varying model for $n_E$ spanning approximately one solar cycle period. Finally, we discuss the importance of our model for chromatic noise mitigation in gravitational-wave analyses of pulsar timing data and the potential of developing synergies between sophisticated PTA solar electron density models and those developed by the solar physics community.
△ Less
Submitted 17 November, 2021;
originally announced November 2021.
-
The NANOGrav 12.5-Year Data Set: Polarimetry and Faraday Rotation Measures from Observations of Millisecond Pulsars with the Green Bank Telescope
Authors:
Haley M. Wahl,
Maura McLaughlin,
Peter A. Gentile,
Megan L. Jones,
Renée Spiewak,
Zaven Arzoumanian,
Kathryn Crowter,
Paul Demorest,
Megan E. DeCesar,
Timothy Dolch,
Justin A. Ellis,
Robert D. Ferdman,
Elizabeth C. Ferrara,
Emmanuel Fonseca,
Nate Garver-Daniels,
Glenn Jones,
Michael T. Lam,
Lina Levin,
Natalia Lewandowska,
Duncan Lorimer,
Ryan S. Lynch,
Dustin R. Madison,
Cherry Ng,
David J. Nice,
Timothy T. Pennucci
, et al. (6 additional authors not shown)
Abstract:
In this work, we present polarization profiles for 23 millisecond pulsars observed at 820 MHz and 1500 MHz with the Green Bank Telescope as part of the NANOGrav pulsar timing array. We calibrate the data using Mueller matrix solutions calculated from observations of PSRs B1929+10 and J1022+1001. We discuss the polarization profiles, which can be used to constrain pulsar emission geometry, and pres…
▽ More
In this work, we present polarization profiles for 23 millisecond pulsars observed at 820 MHz and 1500 MHz with the Green Bank Telescope as part of the NANOGrav pulsar timing array. We calibrate the data using Mueller matrix solutions calculated from observations of PSRs B1929+10 and J1022+1001. We discuss the polarization profiles, which can be used to constrain pulsar emission geometry, and present both the first published radio polarization profiles for nine pulsars and the discovery of very low intensity average profile components ("microcomponents") in four pulsars. Using the Faraday rotation measures, we measure for each pulsar and use it to calculate the Galactic magnetic field parallel to the line of sight for different lines of sight through the interstellar medium. We fit for linear and sinusoidal trends in time in the dispersion measure and Galactic magnetic field and detect magnetic field variations with a period of one year in some pulsars, but overall find that the variations in these parameters are more consistent with a stochastic origin.
△ Less
Submitted 6 December, 2022; v1 submitted 12 April, 2021;
originally announced April 2021.
-
First discovery of new pulsars and RRATs with CHIME/FRB
Authors:
D. C. Good,
B. C. Andersen,
P. Chawla,
K. Crowter,
F. Q. Dong,
E. Fonseca,
B. W. Meyers,
C. Ng,
Z. Pleunis,
S. M. Ransom,
I. H. Stairs,
C. M. Tan,
M. Bhardwaj,
P. J. Boyle,
M. Dobbs,
B. M. Gaensler,
V. M. Kaspi,
K. W. Masui,
A. Naidu,
M. Rafiei-Ravandi,
P. Scholz,
K. M. Smith,
S. P. Tendulkar
Abstract:
We report the discovery of seven new Galactic pulsars with the Canadian Hydrogen Intensity Mapping Experiment's Fast Radio Burst backend (CHIME/FRB). These sources were first identified via single pulses in CHIME/FRB, then followed up with CHIME/Pulsar. Four sources appear to be rotating radio transients (RRATs), pulsar-like sources with occasional single pulse emission with an underlying periodic…
▽ More
We report the discovery of seven new Galactic pulsars with the Canadian Hydrogen Intensity Mapping Experiment's Fast Radio Burst backend (CHIME/FRB). These sources were first identified via single pulses in CHIME/FRB, then followed up with CHIME/Pulsar. Four sources appear to be rotating radio transients (RRATs), pulsar-like sources with occasional single pulse emission with an underlying periodicity. Of those four sources, three have detected periods ranging from 220 ms to 2.726 s. Three sources have more persistent but still intermittent emission and are likely intermittent or nulling pulsars. We have determined phase-coherent timing solutions for the latter three. These seven sources are the first discovery of previously unknown Galactic sources with CHIME/FRB and highlight the potential of fast radio burst detection instruments to search for intermittent Galactic radio sources.
△ Less
Submitted 19 November, 2021; v1 submitted 3 December, 2020;
originally announced December 2020.
-
The CHIME Pulsar Project: System Overview
Authors:
CHIME/Pulsar Collaboration,
M. Amiri,
K. M. Bandura,
P. J. Boyle,
C. Brar,
J. F. Cliche,
K. Crowter,
D. Cubranic,
P. B. Demorest,
N. T. Denman,
M. Dobbs,
F. Q. Dong,
M. Fandino,
E. Fonseca,
D. C. Good,
M. Halpern,
A. S. Hill,
C. Höfer,
V. M. Kaspi,
T. L. Landecker,
C. Leung,
H. -H. Lin,
J. Luo,
K. W. Masui,
J. W. McKee
, et al. (20 additional authors not shown)
Abstract:
We present the design, implementation and performance of a digital backend constructed for the Canadian Hydrogen Intensity Mapping Experiment (CHIME) that uses accelerated computing to observe radio pulsars and transient radio sources. When operating, the CHIME correlator outputs 10 independent streams of beamformed data for the CHIME/Pulsar backend that digitally track specified celestial positio…
▽ More
We present the design, implementation and performance of a digital backend constructed for the Canadian Hydrogen Intensity Mapping Experiment (CHIME) that uses accelerated computing to observe radio pulsars and transient radio sources. When operating, the CHIME correlator outputs 10 independent streams of beamformed data for the CHIME/Pulsar backend that digitally track specified celestial positions. Each of these independent streams are processed by the CHIME/Pulsar backend system which can coherently dedisperse, in real-time, up to dispersion measure values of 2500 pc/cm$^{-3}$ . The tracking beams and real-time analysis system are autonomously controlled by a priority-based algorithm that schedules both known sources and positions of interest for observation with observing cadences as small as one day. Given the distribution of known pulsars and radio-transient sources, the CHIME/Pulsar system can monitor up to 900 positions once per sidereal day and observe all sources with declinations greater than $-20^\circ$ once every $\sim$2 weeks. We also discuss the science program enabled through the current modes of data acquisition for CHIME/Pulsar that centers on timing and searching experiments.
△ Less
Submitted 10 June, 2021; v1 submitted 13 August, 2020;
originally announced August 2020.
-
Multi-Messenger Gravitational Wave Searches with Pulsar Timing Arrays: Application to 3C66B Using the NANOGrav 11-year Data Set
Authors:
Zaven Arzoumanian,
Paul T. Baker,
Adam Brazier,
Paul R. Brook,
Sarah Burke-Spolaor,
Bence Becsy,
Maria Charisi,
Shami Chatterjee,
James M. Cordes,
Neil J. Cornish,
Fronefield Crawford,
H. Thankful Cromartie,
Kathryn Crowter,
Megan E. DeCesar,
Paul B. Demorest,
Timothy Dolch,
Rodney D. Elliott,
Justin A. Ellis,
Robert D. Ferdman,
Elizabeth C. Ferrara,
Emmanuel Fonseca,
Nathan Garver-Daniels,
Peter A. Gentile,
Deborah C. Good,
Jeffrey S. Hazboun
, et al. (34 additional authors not shown)
Abstract:
When galaxies merge, the supermassive black holes in their centers may form binaries and, during the process of merger, emit low-frequency gravitational radiation in the process. In this paper we consider the galaxy 3C66B, which was used as the target of the first multi-messenger search for gravitational waves. Due to the observed periodicities present in the photometric and astrometric data of th…
▽ More
When galaxies merge, the supermassive black holes in their centers may form binaries and, during the process of merger, emit low-frequency gravitational radiation in the process. In this paper we consider the galaxy 3C66B, which was used as the target of the first multi-messenger search for gravitational waves. Due to the observed periodicities present in the photometric and astrometric data of the source of the source, it has been theorized to contain a supermassive black hole binary. Its apparent 1.05-year orbital period would place the gravitational wave emission directly in the pulsar timing band. Since the first pulsar timing array study of 3C66B, revised models of the source have been published, and timing array sensitivities and techniques have improved dramatically. With these advances, we further constrain the chirp mass of the potential supermassive black hole binary in 3C66B to less than $(1.65\pm0.02) \times 10^9~{M_\odot}$ using data from the NANOGrav 11-year data set. This upper limit provides a factor of 1.6 improvement over previous limits, and a factor of 4.3 over the first search done. Nevertheless, the most recent orbital model for the source is still consistent with our limit from pulsar timing array data. In addition, we are able to quantify the improvement made by the inclusion of source properties gleaned from electromagnetic data to `blind' pulsar timing array searches. With these methods, it is apparent that it is not necessary to obtain exact a priori knowledge of the period of a binary to gain meaningful astrophysical inferences.
△ Less
Submitted 12 August, 2020; v1 submitted 14 May, 2020;
originally announced May 2020.
-
The GBT 350-MHz Drift Scan Pulsar Survey. III. Detection of a magnetic field in the eclipsing material of PSR J2256-1024
Authors:
Kathryn Crowter,
Ingrid H. Stairs,
Christie A. McPhee,
Anne M. Archibald,
Jason Boyles,
Jason Hessels,
Chen Karako-Argaman,
Duncan R. Lorimer,
Ryan S. Lynch,
Maura A. McLaughlin,
Scott M. Ransom,
Mallory S. E. Roberts,
Kevin Stovall,
Joeri van Leeuwen
Abstract:
We present the first measurement of a non-zero magnetic field in the eclipsing material of a black widow pulsar. Black widows are millisecond pulsars which are ablating their companions; therefore they are often proposed as one potential source of isolated millisecond pulsars. PSR J2256-1024 is an eclipsing black widow discovered at radio wavelengths and later also observed in the X-ray and gamma…
▽ More
We present the first measurement of a non-zero magnetic field in the eclipsing material of a black widow pulsar. Black widows are millisecond pulsars which are ablating their companions; therefore they are often proposed as one potential source of isolated millisecond pulsars. PSR J2256-1024 is an eclipsing black widow discovered at radio wavelengths and later also observed in the X-ray and gamma parts of the spectrum. Here we present the radio timing solution for PSR J2256-1024, polarization profiles at 350, 820, and 1500~MHz and an investigation of changes in the polarization profile due to eclipsing material in the system. In the latter we find evidence of Faraday rotation in the linear polarization shortly after eclipse, measuring a rotation measure of 0.44(6) rad per meter squared and a corresponding line-of-sight magnetic field of 3.5(17) mG.
△ Less
Submitted 31 March, 2020;
originally announced April 2020.
-
Modeling the uncertainties of solar-system ephemerides for robust gravitational-wave searches with pulsar timing arrays
Authors:
M. Vallisneri,
S. R. Taylor,
J. Simon,
W. M. Folkner,
R. S. Park,
C. Cutler,
J. A. Ellis,
T. J. W. Lazio,
S. J. Vigeland,
K. Aggarwal,
Z. Arzoumanian,
P. T. Baker,
A. Brazier,
P. R. Brook,
S. Burke-Spolaor,
S. Chatterjee,
J. M. Cordes,
N. J. Cornish,
F. Crawford,
H. T. Cromartie,
K. Crowter,
M. DeCesar,
P. B. Demorest,
T. Dolch,
R. D. Ferdman
, et al. (39 additional authors not shown)
Abstract:
The regularity of pulsar emissions becomes apparent once we reference the pulses' times of arrivals to the inertial rest frame of the solar system. It follows that errors in the determination of Earth's position with respect to the solar-system barycenter can appear as a time-correlated bias in pulsar-timing residual time series, affecting the searches for low-frequency gravitational waves perform…
▽ More
The regularity of pulsar emissions becomes apparent once we reference the pulses' times of arrivals to the inertial rest frame of the solar system. It follows that errors in the determination of Earth's position with respect to the solar-system barycenter can appear as a time-correlated bias in pulsar-timing residual time series, affecting the searches for low-frequency gravitational waves performed with pulsar timing arrays. Indeed, recent array datasets yield different gravitational-wave background upper limits and detection statistics when analyzed with different solar-system ephemerides. Crucially, the ephemerides do not generally provide usable error representations. In this article we describe the motivation, construction, and application of a physical model of solar-system ephemeris uncertainties, which focuses on the degrees of freedom (Jupiter's orbital elements) most relevant to gravitational-wave searches with pulsar timing arrays. This model, BayesEphem, was used to derive ephemeris-robust results in NANOGrav's 11-yr stochastic-background search, and it provides a foundation for future searches by NANOGrav and other consortia. The analysis and simulations reported here suggest that ephemeris modeling reduces the gravitational-wave sensitivity of the 11-yr dataset; and that this degeneracy will vanish with improved ephemerides and with the longer pulsar timing datasets that will become available in the near future.
△ Less
Submitted 6 January, 2020; v1 submitted 2 January, 2020;
originally announced January 2020.
-
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.
-
The NANOGrav 11-Year Data Set: Limits on Gravitational Wave Memory
Authors:
K. Aggarwal,
Z. Arzoumanian,
P. T. Baker,
A. Brazier,
P. R. Brook,
S. Burke-Spolaor,
S. Chatterjee,
J. M. Cordes,
N. J. Cornish,
F. Crawford,
H. T. Cromartie,
K. Crowter,
M. Decesar,
P. B. Demorest,
T. Dolch,
J. A. Ellis,
R. D. Ferdman,
E. C. Ferrara,
E. Fonseca,
N. Garver-Daniels,
P. Gentile,
D. Good,
J. S. Hazboun,
A. M. Holgado,
E. A. Huerta
, et al. (36 additional authors not shown)
Abstract:
The mergers of supermassive black hole binaries (SMBHBs) promise to be incredible sources of gravitational waves (GWs). While the oscillatory part of the merger gravitational waveform will be outside the frequency sensitivity range of pulsar timing arrays (PTAs), the non-oscillatory GW memory effect is detectable. Further, any burst of gravitational waves will produce GW memory, making memory a us…
▽ More
The mergers of supermassive black hole binaries (SMBHBs) promise to be incredible sources of gravitational waves (GWs). While the oscillatory part of the merger gravitational waveform will be outside the frequency sensitivity range of pulsar timing arrays (PTAs), the non-oscillatory GW memory effect is detectable. Further, any burst of gravitational waves will produce GW memory, making memory a useful probe of unmodeled exotic sources and new physics. We searched the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 11-year data set for GW memory. This dataset is sensitive to very low frequency GWs of $\sim3$ to $400$ nHz (periods of $\sim11$ yr $-$ $1$ mon). Finding no evidence for GWs, we placed limits on the strain amplitude of GW memory events during the observation period. We then used the strain upper limits to place limits on the rate of GW memory causing events. At a strain of $2.5\times10^{-14}$, corresponding to the median upper limit as a function of source sky position, we set a limit on the rate of GW memory events at $<0.4$ yr$^{-1}$. That strain corresponds to a SMBHB merger with reduced mass of $ηM \sim 2\times10^{10}M_\odot$ and inclination of $ι=π/3$ at a distance of 1 Gpc.
As a test of our analysis, we analyzed the NANOGrav 9-year data set as well. This analysis found an anomolous signal, which does not appear in the 11-year data set. This signal is not a GW, and its origin remains unknown.
△ Less
Submitted 6 December, 2019; v1 submitted 19 November, 2019;
originally announced November 2019.
-
The NANOGrav 11-Year Data Set: Evolution of Gravitational Wave Background Statistics
Authors:
J. S. Hazboun,
J. Simon,
S. R. Taylor,
M. T. Lam,
S. J. Vigeland,
K. Islo,
J. S. Key,
Z. Arzoumanian,
P. T. Baker,
A. Brazier,
P. R. Brook,
S. Burke-Spolaor,
S. Chatterjee,
J. M. Cordes,
N. J. Cornish,
F. Crawford,
K. Crowter,
H. T. Cromartie,
M. DeCesar,
P. B. Demorest,
T. Dolch,
J. A. Ellis,
R. D. Ferdman,
E. Ferrara,
E. Fonseca
, et al. (38 additional authors not shown)
Abstract:
An ensemble of inspiraling supermassive black hole binaries should produce a stochastic background of very low frequency gravitational waves. This stochastic background is predicted to be a power law, with a spectral index of -2/3, and it should be detectable by a network of precisely timed millisecond pulsars, widely distributed on the sky. This paper reports a new "time slicing" analysis of the…
▽ More
An ensemble of inspiraling supermassive black hole binaries should produce a stochastic background of very low frequency gravitational waves. This stochastic background is predicted to be a power law, with a spectral index of -2/3, and it should be detectable by a network of precisely timed millisecond pulsars, widely distributed on the sky. This paper reports a new "time slicing" analysis of the 11-year data release from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) using 34 millisecond pulsars. Methods to flag potential "false positive" signatures are developed, including techniques to identify responsible pulsars. Mitigation strategies are then presented. We demonstrate how an incorrect noise model can lead to spurious signals, and show how independently modeling noise across 30 Fourier components, spanning NANOGrav's frequency range, effectively diagnoses and absorbs the excess power in gravitational-wave searches. This results in a nominal, and expected, progression of our gravitational-wave statistics. Additionally we show that the first interstellar medium event in PSR J1713+0747 pollutes the common red noise process with low-spectral index noise, and use a tailored noise model to remove these effects.
△ Less
Submitted 20 September, 2019; v1 submitted 18 September, 2019;
originally announced September 2019.
-
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.
-
The NANOGrav 11-Year Data Set: Limits on Gravitational Waves from Individual Supermassive Black Hole Binaries
Authors:
K. Aggarwal,
Z. Arzoumanian,
P. T. Baker,
A. Brazier,
M. R. Brinson,
P. R. Brook,
S. Burke-Spolaor,
S. Chatterjee,
J. M. Cordes,
N. J. Cornish,
F. Crawford,
K. Crowter,
H. T. Cromartie,
M. DeCesar,
P. B. Demorest,
T. Dolch,
J. A. Ellis,
R. D. Ferdman,
E. Ferrara,
E. Fonseca,
N. Garver-Daniels,
P. Gentile,
J. S. Hazboun,
A. M. Holgado,
E. A. Huerta
, et al. (38 additional authors not shown)
Abstract:
Observations indicate that nearly all galaxies contain supermassive black holes (SMBHs) at their centers. When galaxies merge, their component black holes form SMBH binaries (SMBHBs), which emit low-frequency gravitational waves (GWs) that can be detected by pulsar timing arrays (PTAs). We have searched the recently-released North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 1…
▽ More
Observations indicate that nearly all galaxies contain supermassive black holes (SMBHs) at their centers. When galaxies merge, their component black holes form SMBH binaries (SMBHBs), which emit low-frequency gravitational waves (GWs) that can be detected by pulsar timing arrays (PTAs). We have searched the recently-released North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 11-year data set for GWs from individual SMBHBs in circular orbits. As we did not find strong evidence for GWs in our data, we placed 95\% upper limits on the strength of GWs from such sources as a function of GW frequency and sky location. We placed a sky-averaged upper limit on the GW strain of $h_0 < 7.3(3) \times 10^{-15}$ at $f_\mathrm{gw}= 8$ nHz. We also developed a technique to determine the significance of a particular signal in each pulsar using ``dropout' parameters as a way of identifying spurious signals in measurements from individual pulsars. We used our upper limits on the GW strain to place lower limits on the distances to individual SMBHBs. At the most-sensitive sky location, we ruled out SMBHBs emitting GWs with $f_\mathrm{gw}= 8$ nHz within 120 Mpc for $\mathcal{M} = 10^9 \, M_\odot$, and within 5.5 Gpc for $\mathcal{M} = 10^{10} \, M_\odot$. We also determined that there are no SMBHBs with $\mathcal{M} > 1.6 \times 10^9 \, M_\odot$ emitting GWs in the Virgo Cluster. Finally, we estimated the number of potentially detectable sources given our current strain upper limits based on galaxies in Two Micron All-Sky Survey (2MASS) and merger rates from the Illustris cosmological simulation project. Only 34 out of 75,000 realizations of the local Universe contained a detectable source, from which we concluded it was unsurprising that we did not detect any individual sources given our current sensitivity to GWs.
△ Less
Submitted 21 May, 2019; v1 submitted 30 December, 2018;
originally announced December 2018.
-
The NANOGrav 11-year Data Set: Pulse Profile Variability
Authors:
P. R. Brook,
A. Karastergiou,
M. A. McLaughlin,
M. T. Lam,
Z. Arzoumanian,
S. Chatterjee,
J. M. Cordes,
K. Crowter,
M. DeCesar,
P. B. Demorest,
T. Dolch,
J. A. Ellis,
R. D. Ferdman,
E. Ferrara,
E. Fonseca,
P. A. Gentile,
G. Jones,
M. L. Jones,
T. J. W. Lazio,
L. Levin,
D. R. Lorimer,
R. S. Lynch,
C. Ng,
D. J. Nice,
T. T. Pennucci
, et al. (8 additional authors not shown)
Abstract:
Access to 50 years of data has led to the discovery of pulsar emission and rotation variability on timescales of months and years. Most of this long-term variability has been seen in long-period pulsars, with relatively little focus on recycled millisecond pulsars. We have analyzed a 38-pulsar sub-set of the 45 millisecond pulsars in the NANOGrav 11-year data set, in order to review their pulse pr…
▽ More
Access to 50 years of data has led to the discovery of pulsar emission and rotation variability on timescales of months and years. Most of this long-term variability has been seen in long-period pulsars, with relatively little focus on recycled millisecond pulsars. We have analyzed a 38-pulsar sub-set of the 45 millisecond pulsars in the NANOGrav 11-year data set, in order to review their pulse profile stability. The most variability, on any timescale, is seen in PSRs J1713+0747, B1937+21 and J2145-0750. The strongest evidence for long-timescale pulse profile changes is seen in PSRs B1937+21 and J1643-1224. We have focused our analyses on these four pulsars in an attempt to elucidate the causes of their profile variability. Effects of scintillation seem to be responsible for the profile modifications of PSR J2145-0750. We see evidence that imperfect polarization calibration contributes to the profile variability of PSRs J1713+0747 and B1937+21, along with radio frequency interference around 2 GHz, but find that propagation effects also have an influence. The changes seen in PSR J1643-1224 have been reported previously, yet elude explanation beyond their astrophysical nature. Regardless of cause, unmodeled pulse profile changes are detrimental to the accuracy of pulsar timing and must be incorporated into the timing models where possible.
△ Less
Submitted 14 November, 2018; v1 submitted 18 October, 2018;
originally announced October 2018.
-
The NANOGrav 11-year Data Set: Solar Wind Sounding Through Pulsar Timing
Authors:
D. R. Madison,
J. M. Cordes,
Z. Arzoumanian,
S. Chatterjee,
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. M. F. Mingarelli,
C. Ng,
D. J. Nice,
T. T. Pennucci,
S. M. Ransom
, et al. (6 additional authors not shown)
Abstract:
The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) has observed dozens of millisecond pulsars for over a decade. We have accrued a large collection of dispersion measure (DM) measurements sensitive to the total electron content between Earth and the pulsars at each observation. All lines of sight cross through the solar wind which produces correlated DM fluctuations in all…
▽ More
The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) has observed dozens of millisecond pulsars for over a decade. We have accrued a large collection of dispersion measure (DM) measurements sensitive to the total electron content between Earth and the pulsars at each observation. All lines of sight cross through the solar wind which produces correlated DM fluctuations in all pulsars. We develop and apply techniques for extracting the imprint of the solar wind from the full collection of DM measurements in the recently released NANOGrav 11-yr data set. We filter out long time scale DM fluctuations attributable to structure in the interstellar medium and carry out a simultaneous analysis of all pulsars in our sample that can differentiate the correlated signature of the wind from signals unique to individual lines of sight. When treating the solar wind as spherically symmetric and constant in time, we find the electron number density at 1 A.U. to be $7.9\pm0.2$ cm$^{-3}$. Our data shows little evidence of long-term variation in the density of the wind. We argue that our techniques paired with a high cadence, low radio frequency observing campaign of near-ecliptic pulsars would be capable of mapping out large-scale latitudinal structure in the wind.
△ Less
Submitted 21 August, 2018;
originally announced August 2018.
-
The NANOGrav 11-Year Data Set: Arecibo Observatory Polarimetry And Pulse Microcomponents
Authors:
Peter A. Gentile,
Maura A. McLaughlin,
Paul B. Demorest,
Ingrid H. Stairs,
Zaven Arzoumanian,
Kathryn Crowter,
Timothy Dolch,
Megan E. DeCesar,
Justin A. Ellis,
Robert D. Ferdman,
Elizabeth C. Ferrara,
Emmanuel Fonseca,
Marjorie E. Gonzalez,
Glenn Jones,
Megan L. Jones,
Michael T. Lam,
Lina Levin,
Duncan R. Lorimer,
Ryan S. Lynch,
Cherry Ng,
David J. Nice,
Timothy T. Pennucci,
Scott M. Ransom,
Paul S. Ray,
Ren'ee Spiewak
, et al. (3 additional authors not shown)
Abstract:
We present the polarization pulse profiles for 28 pulsars observed with the Arecibo Observatory by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) timing project at 2.1 GHz, 1.4 GHz, and 430 MHz. These profiles represent some of the most sensitive polarimetric millisecond pulsar profiles to date, revealing the existence of microcomponents (that is, pulse components with…
▽ More
We present the polarization pulse profiles for 28 pulsars observed with the Arecibo Observatory by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) timing project at 2.1 GHz, 1.4 GHz, and 430 MHz. These profiles represent some of the most sensitive polarimetric millisecond pulsar profiles to date, revealing the existence of microcomponents (that is, pulse components with peak intensities much lower than the total pulse peak intensity). Although microcomponents have been detected in some pulsars previously, we present microcomponents for PSRs B1937+21, J1713+0747, and J2234+0944 for the first time. These microcomponents can have an impact on pulsar timing, geometry, and flux density determination. We present rotation measures for all 28 pulsars, determined independently at different observation frequencies and epochs, and find the Galactic magnetic fields derived from these rotation measures to be consistent with current models. These polarization profiles were made using measurement equation template matching, which allows us to generate the polarimetric response of the Arecibo Observatory on an epoch-by-epoch basis. We use this method to describe its time variability, and find that the polarimetric responses of the Arecibo Observatory's 1.4 and 2.1 GHz receivers vary significantly with time.
△ Less
Submitted 2 July, 2018;
originally announced July 2018.
-
Tests of Gravitational Symmetries with Pulsar Binary J1713+0747
Authors:
W. W. Zhu,
G. Desvignes,
N. Wex,
R. N. Caballero,
D. J. Champion,
P. B. Demorest,
J. A. Ellis,
G. H. Janssen,
M. Kramer,
A. Krieger,
L. Lentati,
D. J. Nice,
S. M. Ransom,
I. H. Stairs,
B. W. Stappers,
J. P. W. Verbiest,
Z. Arzoumanian,
C. G. Bassa,
M. Burgay,
I. Cognard,
K. Crowter,
T. Dolch,
R. D. Ferdman,
E. Fonseca,
M. E. Gonzalez
, et al. (28 additional authors not shown)
Abstract:
Symmetries play an important role in modern theories of gravity. The strong equivalence principle (SEP) constitutes a collection of gravitational symmetries which are all implemented by general relativity. Alternative theories, however, are generally expected to violate some aspects of SEP. We test three aspects of SEP using observed change rates in the orbital period and eccentricity of binary pu…
▽ More
Symmetries play an important role in modern theories of gravity. The strong equivalence principle (SEP) constitutes a collection of gravitational symmetries which are all implemented by general relativity. Alternative theories, however, are generally expected to violate some aspects of SEP. We test three aspects of SEP using observed change rates in the orbital period and eccentricity of binary pulsar J1713+0747: 1. the gravitational constant's constancy as part of locational invariance of gravitation; 2. the post-Newtonian parameter $\hatα_3$ in gravitational Lorentz invariance; 3. the universality of free fall (UFF) for strongly self-gravitating bodies. Based on the pulsar timing result of the combined dataset from the North American Nanohertz Gravitational Observatory (NANOGrav) and the European Pulsar Timing Array (EPTA), we find $\dot{G}/G = (-0.1 \pm 0.9) \times 10^{-12}\,{\rm yr}^{-1}$, which is weaker than Solar system limits, but applies for strongly self-gravitating objects. Furthermore, we obtain the constraints $|Δ|< 0.002$ for the UFF test and $-3\times10^{-20} < \hatα_3 < 4\times10^{-20}$ at 95% confidence. These are the first direct UFF and $\hatα_3$ tests based on pulsar binaries, and they overcome various limitations of previous tests.
△ Less
Submitted 26 February, 2018;
originally announced February 2018.
-
The NANOGrav 11-year Data Set: Pulsar-timing Constraints On The Stochastic Gravitational-wave Background
Authors:
Z. Arzoumanian,
P. T. Baker,
A. Brazier,
S. Burke-Spolaor,
S. J. Chamberlin,
S. Chatterjee,
B. Christy,
J. M. Cordes,
N. J. Cornish,
F. Crawford,
H. Thankful Cromartie,
K. Crowter,
M. DeCesar,
P. B. Demorest,
T. Dolch,
J. A. Ellis,
R. D. Ferdman,
E. Ferrara,
W. M. Folkner,
E. Fonseca,
N. Garver-Daniels,
P. A. Gentile,
R. Haas,
J. S. Hazboun,
E. A. Huerta
, et al. (35 additional authors not shown)
Abstract:
We search for an isotropic stochastic gravitational-wave background (GWB) in the newly released $11$-year dataset from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). While we find no significant evidence for a GWB, we place constraints on a GWB from a population of supermassive black-hole binaries, cosmic strings, and a primordial GWB. For the first time, we find that…
▽ More
We search for an isotropic stochastic gravitational-wave background (GWB) in the newly released $11$-year dataset from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). While we find no significant evidence for a GWB, we place constraints on a GWB from a population of supermassive black-hole binaries, cosmic strings, and a primordial GWB. For the first time, we find that the GWB upper limits and detection statistics are sensitive to the Solar System ephemeris (SSE) model used, and that SSE errors can mimic a GWB signal. We developed an approach that bridges systematic SSE differences, producing the first PTA constraints that are robust against SSE uncertainties. We thus place a $95\%$ upper limit on the GW strain amplitude of $A_\mathrm{GWB}<1.45\times 10^{-15}$ at a frequency of $f=1$ yr$^{-1}$ for a fiducial $f^{-2/3}$ power-law spectrum, and with inter-pulsar correlations modeled. This is a factor of $\sim 2$ improvement over the NANOGrav $9$-year limit, calculated using the same procedure. Previous PTA upper limits on the GWB will need revision in light of SSE systematic uncertainties. We use our constraints to characterize the combined influence on the GWB of the stellar mass-density in galactic cores, the eccentricity of SMBH binaries, and SMBH--galactic-bulge scaling relationships. We constrain cosmic-string tension using recent simulations, yielding an SSE-marginalized $95\%$ upper limit on the cosmic string tension of $Gμ< 5.3\times 10^{-11}$---a factor of $\sim 2$ better than the published NANOGrav $9$-year constraints. Our SSE-marginalized $95\%$ upper limit on the energy density of a primordial GWB (for a radiation-dominated post-inflation Universe) is $Ω_\mathrm{GWB}(f)h^2<3.4\times10^{-10}$.
△ Less
Submitted 7 June, 2018; v1 submitted 8 January, 2018;
originally announced January 2018.
-
The NANOGrav 11-year Data Set: High-precision timing of 45 Millisecond Pulsars
Authors:
Zaven Arzoumanian,
Adam Brazier,
Sarah Burke-Spolaor,
Sydney Chamberlin,
Shami Chatterjee,
Brian Christy,
James M. Cordes,
Neil J. Cornish,
Fronefield Crawford,
H. Thankful Cromartie,
Kathryn Crowter,
Megan E. DeCesar,
Paul B. Demorest,
Timothy Dolch,
Justin A. Ellis,
Robert D. Ferdman,
Elizabeth C. Ferrara,
Emmanuel Fonseca,
Nathan Garver-Daniels,
Peter A. Gentile,
Daniel Halmrast,
Eliu Huerta,
Fredrick A. Jenet,
Cody Jessup,
Glenn Jones
, et al. (31 additional authors not shown)
Abstract:
We present high-precision timing data over time spans of up to 11 years for 45 millisecond pulsars observed as part of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) project, aimed at detecting and characterizing low-frequency gravitational waves. The pulsars were observed with the Arecibo Observatory and/or the Green Bank Telescope at frequencies ranging from 327 MHz…
▽ More
We present high-precision timing data over time spans of up to 11 years for 45 millisecond pulsars observed as part of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) project, aimed at detecting and characterizing low-frequency gravitational waves. The pulsars were observed with the Arecibo Observatory and/or the Green Bank Telescope at frequencies ranging from 327 MHz to 2.3 GHz. Most pulsars were observed with approximately monthly cadence, with six high--timing-precision pulsars observed weekly, and all were observed at widely separated frequencies at each observing epoch in order to fit for time-variable dispersion delays. We describe our methods for data processing, time-of-arrival (TOA) calculation, and the implementation of a new, automated method for removing outlier TOAs. We fit a timing model for each pulsar that includes spin, astrometric, and, if necessary, binary parameters, in addition to time-variable dispersion delays and parameters that quantify pulse-profile evolution with frequency. The new timing solutions provide three new parallax measurements, two new Shapiro delay measurements, and two new measurements of large orbital-period variations. We fit models that characterize sources of noise for each pulsar. We find that 11 pulsars show significant red noise, with generally smaller spectral indices than typically measured for non-recycled pulsars, possibly suggesting a different origin. Future papers will use these data to constrain or detect the signatures of gravitational-wave signals.
△ Less
Submitted 10 April, 2018; v1 submitted 23 December, 2017;
originally announced January 2018.
-
The NANOGrav Nine-Year Data Set: Measurement and Interpretation of Variations in Dispersion Measures
Authors:
M. L. Jones,
M. A. McLaughlin,
M. T. Lam,
J. M. Cordes,
L. Levin,
S. Chatterjee,
Z. Arzoumanian,
K. Crowter,
P. B. Demorest,
T. Dolch,
J. A. Ellis,
R. D. Ferdman,
E. Fonseca,
M. E. Gonzalez,
G. Jones,
T. J. W. Lazio,
D. J. Nice,
T. T. Pennucci,
S. M. Ransom,
D. R. Stinebring,
I. H. Stairs,
K. Stovall,
J. K. Swiggum,
W. W. Zhu
Abstract:
We analyze dispersion measure (DM) variations of 37 millisecond pulsars in the 9-year NANOGrav data release and constrain the sources of these variations. Variations are significant for nearly all pulsars, with characteristic timescales comparable to or even shorter than the average spacing between observations. Five pulsars have periodic annual variations, 14 pulsars have monotonically increasing…
▽ More
We analyze dispersion measure (DM) variations of 37 millisecond pulsars in the 9-year NANOGrav data release and constrain the sources of these variations. Variations are significant for nearly all pulsars, with characteristic timescales comparable to or even shorter than the average spacing between observations. Five pulsars have periodic annual variations, 14 pulsars have monotonically increasing or decreasing trends, and 13 pulsars show both effects. Several pulsars show correlations between DM excesses and lines of sight that pass close to the Sun. Mapping of the DM variations as a function of the pulsar trajectory can identify localized ISM features and, in one case, an upper limit to the size of the dispersing region of 13.2 AU. Finally, five pulsars show very nearly quadratic structure functions, which could be indicative of an underlying Kolmogorov medium. Four pulsars show roughly Kolmogorov structure functions and another four show structure functions less steep than Kolmogorov. One pulsar has too large an uncertainty to allow comparisons. We discuss explanations for apparent departures from a Kolmogorov-like spectrum, and show that the presence of other trends in the data is the most likely cause.
△ Less
Submitted 9 December, 2016;
originally announced December 2016.
-
The NANOGrav Nine-Year Data Set: Excess Noise in Millisecond Pulsar Arrival Times
Authors:
M. T. Lam,
J. M. Cordes,
S. Chatterjee,
Z. Arzoumanian,
K. Crowter,
P. B. Demorest,
T. Dolch,
J. A. Ellis,
R. D. Ferdman,
E. Fonseca,
M. E. Gonzalez,
G. Jones,
M. L. Jones,
L. Levin,
D. R. Madison,
M. A. McLaughlin,
D. J. Nice,
T. T. Pennucci,
S. M. Ransom,
R. M. Shannon,
X. Siemens,
I. H. Stairs,
K. Stovall,
J. K. Swiggum,
W. W. Zhu
Abstract:
Gravitational wave astronomy using a pulsar timing array requires high-quality millisecond pulsars, correctable interstellar propagation delays, and high-precision measurements of pulse times of arrival. Here we identify noise in timing residuals that exceeds that predicted for arrival time estimation for millisecond pulsars observed by the North American Nanohertz Observatory for Gravitational Wa…
▽ More
Gravitational wave astronomy using a pulsar timing array requires high-quality millisecond pulsars, correctable interstellar propagation delays, and high-precision measurements of pulse times of arrival. Here we identify noise in timing residuals that exceeds that predicted for arrival time estimation for millisecond pulsars observed by the North American Nanohertz Observatory for Gravitational Waves. We characterize the excess noise using variance and structure function analyses. We find that 26 out of 37 pulsars show inconsistencies with a white-noise-only model based on the short timescale analysis of each pulsar and we demonstrate that the excess noise has a red power spectrum for 15 pulsars. We also decompose the excess noise into chromatic (radio-frequency-dependent) and achromatic components. Associating the achromatic red-noise component with spin noise and including additional power-spectrum-based estimates from the literature, we estimate a scaling law in terms of spin parameters (frequency and frequency derivative) and data-span length and compare it to the scaling law of Shannon \& Cordes (2010). We briefly discuss our results in terms of detection of gravitational waves at nanohertz frequencies.
△ Less
Submitted 2 December, 2016; v1 submitted 6 October, 2016;
originally announced October 2016.
-
PSR J1024-0719: A Millisecond Pulsar in an Unusual Long-Period Orbit
Authors:
D. L. Kaplan,
T. Kupfer,
D. J. Nice,
A. Irrgang,
U. Heber,
Z. Arzoumanian,
E. Beklen,
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,
S. Kreuzer,
M. T. Lam,
L. Levin,
D. R. Lorimer,
R. S. Lynch,
M. A. McLaughlin,
A. A. Miller
, et al. (10 additional authors not shown)
Abstract:
PSR J1024$-$0719 is a millisecond pulsar that was long thought to be isolated. However, puzzling results concerning its velocity, distance, and low rotational period derivative have led to reexamination of its properties. We present updated radio timing observations along with new and archival optical data that show PSR J1024$-$0719 is most likely in a long period (2$-$20 kyr) binary system with a…
▽ More
PSR J1024$-$0719 is a millisecond pulsar that was long thought to be isolated. However, puzzling results concerning its velocity, distance, and low rotational period derivative have led to reexamination of its properties. We present updated radio timing observations along with new and archival optical data that show PSR J1024$-$0719 is most likely in a long period (2$-$20 kyr) binary system with a low-mass ($\approx 0.4\,M_\odot$) low-metallicity ($Z \approx -0.9\,$ dex) main sequence star. Such a system can explain most of the anomalous properties of this pulsar. We suggest that this system formed through a dynamical exchange in a globular cluster that ejected it into a halo orbit, consistent with the low observed metallicity for the stellar companion. Further astrometric and radio timing observations such as measurement of the third period derivative could strongly constrain the range of orbital parameters.
△ Less
Submitted 15 July, 2016; v1 submitted 1 April, 2016;
originally announced April 2016.
-
The NANOGrav Nine-year Data Set: Mass and Geometric Measurements of Binary Millisecond Pulsars
Authors:
Emmanuel Fonseca,
Timothy T. Pennucci,
Justin A. Ellis,
Ingrid H. Stairs,
David J. Nice,
Scott M. Ransom,
Paul B. Demorest,
Zaven Arzoumanian,
Kathryn Crowter,
Timothy Dolch,
Robert D. Ferdman,
Marjorie E. Gonzalez,
Glenn Jones,
Megan L. Jones,
Michael T. Lam,
Lina Levin,
Maura A. McLaughlin,
Kevin Stovall,
Joseph K. Swiggum,
Weiwei Zhu
Abstract:
We analyze 24 binary radio pulsars in the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) nine-year data set. We make fourteen significant measurements of Shapiro delay, including new detections in four pulsar-binary systems (PSRs J0613$-$0200, J2017+0603, J2302+4442, and J2317+1439), and derive estimates of the binary-component masses and orbital inclination for these MSP-…
▽ More
We analyze 24 binary radio pulsars in the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) nine-year data set. We make fourteen significant measurements of Shapiro delay, including new detections in four pulsar-binary systems (PSRs J0613$-$0200, J2017+0603, J2302+4442, and J2317+1439), and derive estimates of the binary-component masses and orbital inclination for these MSP-binary systems. We find a wide range of binary pulsar masses, with values as low as $m_{\rm p} = 1.18^{+0.10}_{-0.09}\text{M}_{\odot}$ for PSR J1918$-$0642 and as high as $m_{\rm p} = 1.928^{+0.017}_{-0.017}\text{M}_{\odot}$ for PSR J1614$-$2230 (both 68.3\% credibility). We make an improved measurement of the Shapiro timing delay in the PSR J1918$-$0642 and J2043+1711 systems, measuring the pulsar mass in the latter system to be $m_{\rm p} = 1.41^{+0.21}_{-0.18}\text{M}_{\odot}$ (68.3\% credibility) for the first time. We measure secular variations of one or more orbital elements in many systems, and use these measurements to further constrain our estimates of the pulsar and companion masses whenever possible. In particular, we used the observed Shapiro delay and periastron advance due to relativistic gravity in the PSR J1903+0327 system to derive a pulsar mass of $m_{\rm p} = 1.65^{+0.02}_{-0.02}\text{M}_{\odot}$ (68.3\% credibility). We discuss the implications that our mass measurements have on the overall neutron-star mass distribution, and on the "mass/orbital-period" correlation due to extended mass transfer.
△ Less
Submitted 4 October, 2016; v1 submitted 1 March, 2016;
originally announced March 2016.
-
The NANOGrav Nine-year Data Set: Monitoring Interstellar Scattering Delays
Authors:
Lina Levin,
Maura A. McLaughlin,
Glenn Jones,
James M. Cordes,
Daniel R. Stinebring,
Shami Chatterjee,
Timothy Dolch,
Michael T. Lam,
T. Joseph W. Lazio,
Nipuni Palliyaguru,
Zaven Arzoumanian,
Kathryn Crowter,
Paul B. Demorest,
Justin A. Ellis,
Robert D. Ferdman,
Emmanuel Fonseca,
Marjorie E. Gonzalez,
Megan L. Jones,
David J. Nice,
Timothy T. Pennucci,
Scott M. Ransom,
Ingrid H. Stairs,
Kevin Stovall,
Joseph K. Swiggum,
Weiwei Zhu
Abstract:
We report on an effort to extract and monitor interstellar scintillation parameters in regular timing observations collected for the NANOGrav pulsar timing array. Scattering delays are measured by creating dynamic spectra for each pulsar and observing epoch of wide-band observations centered near 1500 MHz and carried out at the Green Bank Telescope and the Arecibo Observatory. The ~800-MHz wide fr…
▽ More
We report on an effort to extract and monitor interstellar scintillation parameters in regular timing observations collected for the NANOGrav pulsar timing array. Scattering delays are measured by creating dynamic spectra for each pulsar and observing epoch of wide-band observations centered near 1500 MHz and carried out at the Green Bank Telescope and the Arecibo Observatory. The ~800-MHz wide frequency bands imply dramatic changes in scintillation bandwidth across the bandpass, and a stretching routine has been included to account for this scaling. For most of the 10 pulsars for which the scaling has been measured, the bandwidths scale with frequency less steeply than expected for a Kolmogorov medium. We find estimated scattering delay values that vary with time by up to an order of magnitude. The mean measured scattering delays are similar to previously published values and slightly higher than predicted by interstellar medium models. We investigate the possibility of increasing the timing precision by mitigating timing errors introduced by the scattering delays. For most of the pulsars, the uncertainty in the time of arrival of a single timing point is much larger than the maximum variation of the scattering delay, suggesting that diffractive scintillation remains only a negligible part of their noise budget.
△ Less
Submitted 18 January, 2016;
originally announced January 2016.
-
The NANOGrav Nine-Year Data Set: Noise Budget for Pulsar Arrival Times on Intraday Timescales
Authors:
M. T. Lam,
J. M. Cordes,
S. Chatterjee,
Z. Arzoumanian,
K. Crowter,
P. B. Demorest,
T. Dolch,
J. A. Ellis,
R. D. Ferdman,
E. F. Fonseca,
M. E. Gonzalez,
G. Jones,
M. L. Jones,
L. Levin,
D. R. Madison,
M. A. McLaughlin,
D. J. Nice,
T. T. Pennucci,
S. M. Ransom,
X. Siemens,
I. H. Stairs,
K. Stovall,
J. K. Swiggum,
W. W. Zhu
Abstract:
The use of pulsars as astrophysical clocks for gravitational wave experiments demands the highest possible timing precision. Pulse times of arrival (TOAs) are limited by stochastic processes that occur in the pulsar itself, along the line of sight through the interstellar medium, and in the measurement process. On timescales of seconds to hours, the TOA variance exceeds that from template-fitting…
▽ More
The use of pulsars as astrophysical clocks for gravitational wave experiments demands the highest possible timing precision. Pulse times of arrival (TOAs) are limited by stochastic processes that occur in the pulsar itself, along the line of sight through the interstellar medium, and in the measurement process. On timescales of seconds to hours, the TOA variance exceeds that from template-fitting errors due to additive noise. We assess contributions to the total variance from two additional effects: amplitude and phase jitter intrinsic to single pulses and changes in the interstellar impulse response from scattering. The three effects have different dependencies on time, frequency, and pulse signal-to-noise ratio. We use data on 37 pulsars from the North American Nanohertz Observatory for Gravitational Waves to assess the individual contributions to the overall intraday noise budget for each pulsar. We detect jitter in 22 pulsars and estimate the average value of RMS jitter in our pulsars to be $\sim 1\%$ of pulse phase. We examine how jitter evolves as a function of frequency and find evidence for evolution. Finally, we compare our measurements with previous noise parameter estimates and discuss methods to improve gravitational wave detection pipelines.
△ Less
Submitted 28 September, 2016; v1 submitted 28 December, 2015;
originally announced December 2015.
-
The NANOGrav Nine-year Data Set: Astrometric Measurements of 37 Millisecond Pulsars
Authors:
Allison M. Matthews,
David J. Nice,
Emmanuel Fonseca,
Zaven Arzoumanian,
Kathryn Crowter,
Paul B. Demorest,
Timothy Dolch,
Justin A. Ellis,
Robert D. Ferdman,
Marjorie E. Gonzalez,
Glenn Jones,
Megan L. Jones,
Michael T. Lam,
Lina Levin,
Maura A. McLaughlin,
Timothy T. Pennucci,
Scott M. Ransom,
Ingrid H. Stairs,
Kevin Stovall,
Joseph K. Swiggum,
Weiwei Zhu
Abstract:
Using the nine-year radio-pulsar timing data set from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), collected at Arecibo Observatory and the Green Bank Telescope, we have measured the positions, proper motions, and parallaxes for 37 millisecond pulsars. We report twelve significant parallax measurements and distance measurements, and eighteen lower limits on distance…
▽ More
Using the nine-year radio-pulsar timing data set from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), collected at Arecibo Observatory and the Green Bank Telescope, we have measured the positions, proper motions, and parallaxes for 37 millisecond pulsars. We report twelve significant parallax measurements and distance measurements, and eighteen lower limits on distance. We compare these measurements to distances predicted by the NE2001 interstellar electron density model and find them to be in general agreement. We use measured orbital-decay rates and spin-down rates to confirm two of the parallax distances and to place distance upper limits on other sources; these distance limits agree with the parallax distances with one exception, PSR J1024-0719, which we discuss at length. Using the proper motions of the 37 NANOGrav pulsars in combination with other published measurements, we calculate the velocity dispersion of the millisecond pulsar population in Galactocentric coordinates. We find the radial, azimuthal, and perpendicular dispersions to be 46, 40, and 24 km s-1, respectively, in a model that allows for high-velocity outliers; or 81, 58, and 62 km s-1 for the full population. These velocity dispersions are far smaller than those of the canonical pulsar population, and are similar to older Galactic disk populations. This suggests that millisecond pulsar velocities are largely attributable to their being an old population rather than being artifacts of their birth and evolution as neutron star binary systems. The components of these velocity dispersions follow similar proportions to other Galactic populations, suggesting that our results are not biased by selection effects.
△ Less
Submitted 2 January, 2016; v1 submitted 29 September, 2015;
originally announced September 2015.