-
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.
-
Pulsar Timing Arrays require hierarchical models
Authors:
Rutger van Haasteren
Abstract:
Pulsar Timing Array projects have found evidence of a stochastic background of gravitational waves (GWB) using data from an ensemble of pulsars. In the literature, minimal assumptions are made about the signal and noise processes that affect data from these pulsars, such as pulsar spin noise. These assumptions are encoded as uninformative priors in Bayesian searches, though Frequentist approaches…
▽ More
Pulsar Timing Array projects have found evidence of a stochastic background of gravitational waves (GWB) using data from an ensemble of pulsars. In the literature, minimal assumptions are made about the signal and noise processes that affect data from these pulsars, such as pulsar spin noise. These assumptions are encoded as uninformative priors in Bayesian searches, though Frequentist approaches make similar assumptions. Uninformative priors are not suitable for (noise) properties of pulsars in an ensemble, and they bias estimates of model parameters such as gravitational-wave signal parameters. Both Frequentist and Bayesian searches are affected. In this letter, more appropriate priors are proposed in the language of Hierarchical Bayesian Modeling, where the properties of the ensemble of pulsars are jointly described with the properties of the individual components of the ensemble. Results by Pulsar Timing Array projects should be re-evaluated using Hierarchical Models.
△ Less
Submitted 7 June, 2024;
originally announced June 2024.
-
PINT: Maximum-likelihood estimation of pulsar timing noise parameters
Authors:
Abhimanyu Susobhanan,
David Kaplan,
Anne Archibald,
Jing Luo,
Paul Ray,
Timothy Pennucci,
Scott Ransom,
Gabriella Agazie,
William Fiore,
Bjorn Larsen,
Patrick O'Neill,
Rutger van Haasteren,
Akash Anumarlapudi,
Matteo Bachetti,
Deven Bhakta,
Chloe Champagne,
H. Thankful Cromartie,
Paul Demorest,
Ross Jennings,
Matthew Kerr,
Sasha Levina,
Alexander McEwen,
Brent Shapiro-Albert,
Joseph Swiggum
Abstract:
PINT is a pure-Python framework for high-precision pulsar timing developed on top of widely used and well-tested Python libraries, supporting both interactive and programmatic data analysis workflows. We present a new frequentist framework within PINT to characterize the single-pulsar noise processes present in pulsar timing datasets. This framework enables the parameter estimation for both uncorr…
▽ More
PINT is a pure-Python framework for high-precision pulsar timing developed on top of widely used and well-tested Python libraries, supporting both interactive and programmatic data analysis workflows. We present a new frequentist framework within PINT to characterize the single-pulsar noise processes present in pulsar timing datasets. This framework enables the parameter estimation for both uncorrelated and correlated noise processes as well as the model comparison between different timing and noise models in a computationally inexpensive way. We demonstrate the efficacy of the new framework by applying it to simulated datasets as well as a real dataset of PSR B1855+09. We also describe the new features implemented in PINT since it was first described in the literature.
△ Less
Submitted 19 June, 2024; v1 submitted 3 May, 2024;
originally announced May 2024.
-
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: 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.
-
Exploring the Capabilities of Gibbs Sampling in Pulsar Timing Arrays
Authors:
Nima Laal,
William G Lamb,
Joseph D. Romano,
Xavier Siemens,
Stephen R. Taylor,
Rutger van Haasteren
Abstract:
We explore the use of Gibbs sampling in estimating the noise properties of individual pulsars and illustrate its effectiveness using the NANOGrav 11-year data set. We find that Gibbs sampling noise modeling (GM) is more efficient than the current standard Bayesian techniques (SM) for single pulsar analyses by yielding model parameter posteriors with average effective-sample-size ratio (GM/SM) of 6…
▽ More
We explore the use of Gibbs sampling in estimating the noise properties of individual pulsars and illustrate its effectiveness using the NANOGrav 11-year data set. We find that Gibbs sampling noise modeling (GM) is more efficient than the current standard Bayesian techniques (SM) for single pulsar analyses by yielding model parameter posteriors with average effective-sample-size ratio (GM/SM) of 6 across all parameters and pulsars. Furthermore, the output of GM contains posteriors for the Fourier coefficients that can be used to characterize the underlying red noise process of any pulsar's timing residuals, which are absent in current implementations of SM. Through simulations, we demonstrate the potential for such coefficients to measure the spatial cross-correlations between pulsar pairs produced by a gravitational wave background.
△ Less
Submitted 7 September, 2023; v1 submitted 20 May, 2023;
originally announced May 2023.
-
The Need For Speed: Rapid Refitting Techniques for Bayesian Spectral Characterization of the Gravitational Wave Background Using PTAs
Authors:
William G. Lamb,
Stephen R. Taylor,
Rutger van Haasteren
Abstract:
Pulsar timing arrays (PTAs) have recently found evidence for a nanohertz-frequency stochastic gravitational-wave background (SGWB). Constraining its spectral characteristics will reveal its origins. To achieve this, we must understand how data and modeling conditions in each pulsar influence the precision and accuracy of SGWB spectral recovery, typically requiring many Bayesian analyses on real da…
▽ More
Pulsar timing arrays (PTAs) have recently found evidence for a nanohertz-frequency stochastic gravitational-wave background (SGWB). Constraining its spectral characteristics will reveal its origins. To achieve this, we must understand how data and modeling conditions in each pulsar influence the precision and accuracy of SGWB spectral recovery, typically requiring many Bayesian analyses on real data sets and large-scale simulations that are slow and computationally taxing. To combat this, we have developed several new rapid approaches that operate on intermediate SGWB analysis products. These techniques refit SGWB spectral models to previously computed Bayesian posterior estimates of the timing power spectra. We test our new techniques on simulated PTA data sets and the NANOGrav 12.5-year data set, where in the latter our refit posterior achieves a Hellinger distance -- bounded between 0 for identical distributions and 1 for zero overlap -- from the current full production-level pipeline that is < 0.1. Our techniques are ~ $10^2$--$10^4$ times faster than the production-level likelihood and scale much more favorably (sub-linearly) as a PTA is expanded with new pulsars or observations. Our techniques also allow us to demonstrate conclusively that SGWB spectral characterization in PTA data sets is driven by the longest-timed pulsars and the best-measured power spectral densities. Indeed, the common-process spectral properties found in the NANOGrav 12.5-year data set are given by analyzing only the ~14 longest-timed pulsars out of the full 45 pulsar array, and we find that the 'shallowing' of the common-process power-law model occurs when gravitational-wave frequencies higher than ~50 nanohertz are included. The implementation of our techniques is openly available as a software suite to allow fast and flexible PTA SGWB spectral characterization and model selection.
△ Less
Submitted 9 October, 2023; v1 submitted 27 March, 2023;
originally announced March 2023.
-
Efficient Gravitational Wave Searches with Pulsar Timing Arrays using Hamiltonian Monte Carlo
Authors:
Gabriel E. Freedman,
Aaron D. Johnson,
Rutger van Haasteren,
Sarah J. Vigeland
Abstract:
Pulsar timing arrays (PTAs) detect low-frequency gravitational waves (GWs) by looking for correlated deviations in pulse arrival times. Current Bayesian searches use Markov Chain Monte Carlo (MCMC) methods, which struggle to sample the large number of parameters needed to model the PTA and GW signals. As the data span and number of pulsars increase, this problem will only worsen. An alternative Mo…
▽ More
Pulsar timing arrays (PTAs) detect low-frequency gravitational waves (GWs) by looking for correlated deviations in pulse arrival times. Current Bayesian searches use Markov Chain Monte Carlo (MCMC) methods, which struggle to sample the large number of parameters needed to model the PTA and GW signals. As the data span and number of pulsars increase, this problem will only worsen. An alternative Monte Carlo sampling method, Hamiltonian Monte Carlo (HMC), utilizes Hamiltonian dynamics to produce sample proposals informed by first-order gradients of the model likelihood. This in turn allows it to converge faster to high dimensional distributions. We implement HMC as an alternative sampling method in our search for an isotropic stochastic GW background, and show that this method produces equivalent statistical results to similar analyses run with standard MCMC techniques, while requiring 100-200 times fewer samples. We show that the speed of HMC sample generation scales as $\mathcal{O}(N_\mathrm{psr}^{5/4})$ where $N_\mathrm{psr}$ is the number of pulsars, compared to $\mathcal{O}(N_\mathrm{psr}^2)$ for MCMC methods. These factors offset the increased time required to generate a sample using HMC, demonstrating the value of adopting HMC techniques for PTAs.
△ Less
Submitted 7 March, 2023; v1 submitted 2 November, 2022;
originally announced November 2022.
-
PINT: A Modern Software Package for Pulsar Timing
Authors:
Jing Luo,
Scott Ransom,
Paul Demorest,
Paul S. Ray,
Anne Archibald,
Matthew Kerr,
Ross J. Jennings,
Matteo Bachetti,
Rutger van Haasteren,
Chloe A. Champagne,
Jonathan Colen,
Camryn Phillips,
Josef Zimmerman,
Kevin Stovall,
Michael T. Lam,
Fredrick A. Jenet
Abstract:
Over the past few decades, the measurement precision of some pulsar-timing experiments has advanced from ~10 us to ~10 ns, revealing many subtle phenomena. Such high precision demands both careful data handling and sophisticated timing models to avoid systematic error. To achieve these goals, we present PINT (PINT Is Not Tempo3), a high-precision Python pulsar timing data analysis package, which i…
▽ More
Over the past few decades, the measurement precision of some pulsar-timing experiments has advanced from ~10 us to ~10 ns, revealing many subtle phenomena. Such high precision demands both careful data handling and sophisticated timing models to avoid systematic error. To achieve these goals, we present PINT (PINT Is Not Tempo3), a high-precision Python pulsar timing data analysis package, which is hosted on GitHub and available on Python Package Index (PyPI) as pint-pulsar. PINT is well-tested, validated, object-oriented, and modular, enabling interactive data analysis and providing an extensible and flexible development platform for timing applications. It utilizes well-debugged public Python packages (e.g., the NumPy and Astropy libraries) and modern software development schemes (e.g., version control and efficient development with git and GitHub) and a continually expanding test suite for improved reliability, accuracy, and reproducibility. PINT is developed and implemented without referring to, copying, or transcribing the code from other traditional pulsar timing software packages (e.g., TEMPO and TEMPO2) and therefore provides a robust tool for cross-checking timing analyses and simulating pulse arrival times. In this paper, we describe the design, usage, and validation of PINT, and we compare timing results between it and TEMPO and TEMPO2.
△ Less
Submitted 5 January, 2021; v1 submitted 30 November, 2020;
originally announced December 2020.
-
From Bright Binaries To Bumpy Backgrounds: Mapping Realistic Gravitational Wave Skies With Pulsar-Timing Arrays
Authors:
Stephen R. Taylor,
Rutger van Haasteren,
Alberto Sesana
Abstract:
Within the next several years, pulsar-timing array programs will likely usher in the next era of gravitational-wave astronomy through the detection of a stochastic background of nanohertz-frequency gravitational waves, originating from a cosmological population of inspiraling supermassive binary black holes. While the source positions will likely be isotropic to a good approximation, the gravitati…
▽ More
Within the next several years, pulsar-timing array programs will likely usher in the next era of gravitational-wave astronomy through the detection of a stochastic background of nanohertz-frequency gravitational waves, originating from a cosmological population of inspiraling supermassive binary black holes. While the source positions will likely be isotropic to a good approximation, the gravitational-wave angular power distribution will be anisotropic, with the most massive and/or nearby binaries producing signals that may resound above the background. We study such a realistic angular power distribution, developing fast and accurate sky-mapping strategies to localize pixels and extended regions of excess power while simultaneously modeling the background signal from the less massive and more distant ensemble. We find that power anisotropy will be challenging to discriminate from isotropy for realistic gravitational-wave skies, requiring SNR $>10$ in order to favor anisotropy with $10:1$ posterior odds in our case study. Amongst our techniques, modeling the population signal with multiple point sources in addition to an isotropic background provides the most physically-motivated and easily interpreted maps, while spherical-harmonic modeling of the square-root power distribution, $P(\hatΩ)^{1/2}$, performs best in discriminating from overall isotropy. Our techniques are modular and easily incorporated into existing pulsar-timing array analysis pipelines.
△ Less
Submitted 26 August, 2020; v1 submitted 8 June, 2020;
originally announced June 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.
-
A pulsar-based timescale from the International Pulsar Timing Array
Authors:
G. Hobbs,
L. Guo,
R. N. Caballero,
W. Coles,
K. J. Lee,
R. N. Manchester,
D. J. Reardon,
D. Matsakis,
M. L. Tong,
Z. Arzoumanian,
M. Bailes,
C. G. Bassa,
N. D. R. Bhat,
A. Brazier,
S. Burke-Spolaor,
D. J. Champion,
S. Chatterjee,
I. Cognard,
S. Dai,
G. Desvignes,
T. Dolch,
R. D. Ferdman,
E. Graikou,
L. Guillemot,
G. H. Janssen
, et al. (34 additional authors not shown)
Abstract:
We have constructed a new timescale, TT(IPTA16), based on observations of radio pulsars presented in the first data release from the International Pulsar Timing Array (IPTA). We used two analysis techniques with independent estimates of the noise models for the pulsar observations and different algorithms for obtaining the pulsar timescale. The two analyses agree within the estimated uncertainties…
▽ More
We have constructed a new timescale, TT(IPTA16), based on observations of radio pulsars presented in the first data release from the International Pulsar Timing Array (IPTA). We used two analysis techniques with independent estimates of the noise models for the pulsar observations and different algorithms for obtaining the pulsar timescale. The two analyses agree within the estimated uncertainties and both agree with TT(BIPM17), a post-corrected timescale produced by the Bureau International des Poids et Mesures (BIPM). We show that both methods could detect significant errors in TT(BIPM17) if they were present. We estimate the stability of the atomic clocks from which TT(BIPM17) is derived using observations of four rubidium fountain clocks at the US Naval Observatory. Comparing the power spectrum of TT(IPTA16) with that of these fountain clocks suggests that pulsar-based timescales are unlikely to contribute to the stability of the best timescales over the next decade, but they will remain a valuable independent check on atomic timescales. We also find that the stability of the pulsar-based timescale is likely to be limited by our knowledge of solar-system dynamics, and that errors in TT(BIPM17) will not be a limiting factor for the primary goal of the IPTA, which is to search for the signatures of nano-Hertz gravitational waves.
△ Less
Submitted 29 October, 2019;
originally announced October 2019.
-
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 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.
-
Studying the solar system with the International Pulsar Timing Array
Authors:
R. N. Caballero,
Y. J. Guo,
K. J. Lee,
P. Lazarus,
D. J. Champion,
G. Desvignes,
M. Kramer,
K. Plant,
Z. Arzoumanian,
M. Bailes,
C. G. Bassa,
N. D. R. Bhat,
A. Brazier,
M. Burgay,
S. Burke-Spolaor,
S. J. Chamberlin,
S. Chatterjee,
I. Cognard,
J. M. Cordes,
S. Dai,
P. Demorest,
T. Dolch,
R. D. Ferdman,
E. Fonseca,
J. R. Gair
, et al. (55 additional authors not shown)
Abstract:
Pulsar-timing analyses are sensitive to errors in the solar-system ephemerides (SSEs) that timing models utilise to estimate the location of the solar-system barycentre, the quasi-inertial reference frame to which all recorded pulse times-of-arrival are referred. Any error in the SSE will affect all pulsars, therefore pulsar timing arrays (PTAs) are a suitable tool to search for such errors and im…
▽ More
Pulsar-timing analyses are sensitive to errors in the solar-system ephemerides (SSEs) that timing models utilise to estimate the location of the solar-system barycentre, the quasi-inertial reference frame to which all recorded pulse times-of-arrival are referred. Any error in the SSE will affect all pulsars, therefore pulsar timing arrays (PTAs) are a suitable tool to search for such errors and impose independent constraints on relevant physical parameters. We employ the first data release of the International Pulsar Timing Array to constrain the masses of the planet-moons systems and to search for possible unmodelled objects (UMOs) in the solar system. We employ ten SSEs from two independent research groups, derive and compare mass constraints of planetary systems, and derive the first PTA mass constraints on asteroid-belt objects. Constraints on planetary-system masses have been improved by factors of up to 20 from the previous relevant study using the same assumptions, with the mass of the Jovian system measured at 9.5479189(3)$\times10^{-4}$ $M_{\odot}$. The mass of the dwarf planet Ceres is measured at 4.7(4)$\times10^{-10}$ $M_{\odot}$. We also present the first sensitivity curves using real data that place generic limits on the masses of UMOs, which can also be used as upper limits on the mass of putative exotic objects. For example, upper limits on dark-matter clumps are comparable to published limits using independent methods. While the constraints on planetary masses derived with all employed SSEs are consistent, we note and discuss differences in the associated timing residuals and UMO sensitivity curves.
△ Less
Submitted 27 September, 2018;
originally announced September 2018.
-
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.
-
Taming outliers in pulsar-timing datasets with hierarchical likelihoods and Hamiltonian sampling
Authors:
Michele Vallisneri,
Rutger van Haasteren
Abstract:
Pulsar-timing datasets have been analyzed with great success using probabilistic treatments based on Gaussian distributions, with applications ranging from studies of neutron-star structure to tests of general relativity and searches for nanosecond gravitational waves. As for other applications of Gaussian distributions, outliers in timing measurements pose a significant challenge to statistical i…
▽ More
Pulsar-timing datasets have been analyzed with great success using probabilistic treatments based on Gaussian distributions, with applications ranging from studies of neutron-star structure to tests of general relativity and searches for nanosecond gravitational waves. As for other applications of Gaussian distributions, outliers in timing measurements pose a significant challenge to statistical inference, since they can bias the estimation of timing and noise parameters, and affect reported parameter uncertainties. We describe and demonstrate a practical end-to-end approach to perform Bayesian inference of timing and noise parameters robustly in the presence of outliers, and to identify these probabilistically. The method is fully consistent (i.e., outlier-ness probabilities vary in tune with the posterior distributions of the timing and noise parameters), and it relies on the efficient sampling of the hierarchical form of the pulsar-timing likelihood. Such sampling has recently become possible with a "no-U-turn" Hamiltonian sampler coupled to a highly customized reparametrization of the likelihood; this code is described elsewhere, but it is already available online. We recommend our method as a standard step in the preparation of pulsar-timing-array datasets: even if statistical inference is not affected, follow-up studies of outlier candidates can reveal unseen problems in radio observations and timing measurements; furthermore, confidence in the results of gravitational-wave searches will only benefit from stringent statistical evidence that datasets are clean and outlier-free.
△ Less
Submitted 7 September, 2016;
originally announced September 2016.
-
High-precision timing of 42 millisecond pulsars with the European Pulsar Timing Array
Authors:
G. Desvignes,
R. N. Caballero,
L. Lentati,
J. P. W. Verbiest,
D. J. Champion,
B. W. Stappers,
G. H. Janssen,
P. Lazarus,
S. Osłowski,
S. Babak,
C. G. Bassa,
P. Brem,
M. Burgay,
I. Cognard,
J. R. Gair,
E. Graikou,
L. Guillemot,
J. W. T. Hessels,
A. Jessner,
C. Jordan,
R. Karuppusamy,
M. Kramer,
A. Lassus,
K. Lazaridis,
K. J. Lee
, et al. (18 additional authors not shown)
Abstract:
We report on the high-precision timing of 42 radio millisecond pulsars (MSPs) observed by the European Pulsar Timing Array (EPTA). This EPTA Data Release 1.0 extends up to mid-2014 and baselines range from 7-18 years. It forms the basis for the stochastic gravitational-wave background, anisotropic background, and continuous-wave limits recently presented by the EPTA elsewhere. The Bayesian timing…
▽ More
We report on the high-precision timing of 42 radio millisecond pulsars (MSPs) observed by the European Pulsar Timing Array (EPTA). This EPTA Data Release 1.0 extends up to mid-2014 and baselines range from 7-18 years. It forms the basis for the stochastic gravitational-wave background, anisotropic background, and continuous-wave limits recently presented by the EPTA elsewhere. The Bayesian timing analysis performed with TempoNest yields the detection of several new parameters: seven parallaxes, nine proper motions and, in the case of six binary pulsars, an apparent change of the semi-major axis. We find the NE2001 Galactic electron density model to be a better match to our parallax distances (after correction from the Lutz-Kelker bias) than the M2 and M3 models by Schnitzeler (2012). However, we measure an average uncertainty of 80\% (fractional) for NE2001, three times larger than what is typically assumed in the literature. We revisit the transverse velocity distribution for a set of 19 isolated and 57 binary MSPs and find no statistical difference between these two populations. We detect Shapiro delay in the timing residuals of PSRs J1600$-$3053 and J1918$-$0642, implying pulsar and companion masses $m_p=1.22_{-0.35}^{+0.5} \text{M}_{\odot}$, $m_c = 0.21_{-0.04}^{+0.06} \text{M}_{\odot }$ and $m_p=1.25_{-0.4}^{+0.6} \text{M}_{\odot}$, $m_c = 0.23_{-0.05}^{+0.07} \text{M}_{\odot }$, respectively. Finally, we use the measurement of the orbital period derivative to set a stringent constraint on the distance to PSRs J1012$+$5307 and J1909$-$3744, and set limits on the longitude of ascending node through the search of the annual-orbital parallax for PSRs J1600$-$3053 and J1909$-$3744.
△ Less
Submitted 26 February, 2016;
originally announced February 2016.
-
From Spin Noise to Systematics: Stochastic Processes in the First International Pulsar Timing Array Data Release
Authors:
L. Lentati,
R. M. Shannon,
W. A. Coles,
J. P. W. Verbiest,
R. van Haasteren,
J. A. Ellis,
R. N. Caballero,
R. N. Manchester,
Z. Arzoumanian,
S. Babak,
C. G. Bassa,
N. D. R. Bhat,
P. Brem,
M. Burgay,
S. Burke-Spolaor,
D. Champion,
S. Chatterjee,
I. Cognard,
J. M. Cordes,
S. Dai,
P. Demorest,
G. Desvignes,
T. Dolch,
R. D. Ferdman,
E. Fonseca
, et al. (58 additional authors not shown)
Abstract:
We analyse the stochastic properties of the 49 pulsars that comprise the first International Pulsar Timing Array (IPTA) data release. We use Bayesian methodology, performing model selection to determine the optimal description of the stochastic signals present in each pulsar. In addition to spin-noise and dispersion-measure (DM) variations, these models can include timing noise unique to a single…
▽ More
We analyse the stochastic properties of the 49 pulsars that comprise the first International Pulsar Timing Array (IPTA) data release. We use Bayesian methodology, performing model selection to determine the optimal description of the stochastic signals present in each pulsar. In addition to spin-noise and dispersion-measure (DM) variations, these models can include timing noise unique to a single observing system, or frequency band. We show the improved radio-frequency coverage and presence of overlapping data from different observing systems in the IPTA data set enables us to separate both system and band-dependent effects with much greater efficacy than in the individual PTA data sets. For example, we show that PSR J1643$-$1224 has, in addition to DM variations, significant band-dependent noise that is coherent between PTAs which we interpret as coming from time-variable scattering or refraction in the ionised interstellar medium. Failing to model these different contributions appropriately can dramatically alter the astrophysical interpretation of the stochastic signals observed in the residuals. In some cases, the spectral exponent of the spin noise signal can vary from 1.6 to 4 depending upon the model, which has direct implications for the long-term sensitivity of the pulsar to a stochastic gravitational-wave (GW) background. By using a more appropriate model, however, we can greatly improve a pulsar's sensitivity to GWs. For example, including system and band-dependent signals in the PSR J0437$-$4715 data set improves the upper limit on a fiducial GW background by $\sim 60\%$ compared to a model that includes DM variations and spin-noise only.
△ Less
Submitted 16 February, 2016;
originally announced February 2016.
-
The International Pulsar Timing Array: First Data Release
Authors:
J. P. W. Verbiest,
L. Lentati,
G. Hobbs,
R. van Haasteren,
P. B. Demorest,
G. H. Janssen,
J. -B. Wang,
G. Desvignes,
R. N. Caballero,
M. J. Keith,
D. J. Champion,
Z. Arzoumanian,
S. Babak,
C. G. Bassa,
N. D. R. Bhat,
A. Brazier,
P. Brem,
M. Burgay,
S. Burke-Spolaor,
S. J. Chamberlin,
S. Chatterjee,
B. Christy,
I. Cognard,
J. M. Cordes,
S. Dai
, et al. (67 additional authors not shown)
Abstract:
The highly stable spin of neutron stars can be exploited for a variety of (astro-)physical investigations. In particular arrays of pulsars with rotational periods of the order of milliseconds can be used to detect correlated signals such as those caused by gravitational waves. Three such "Pulsar Timing Arrays" (PTAs) have been set up around the world over the past decades and collectively form the…
▽ More
The highly stable spin of neutron stars can be exploited for a variety of (astro-)physical investigations. In particular arrays of pulsars with rotational periods of the order of milliseconds can be used to detect correlated signals such as those caused by gravitational waves. Three such "Pulsar Timing Arrays" (PTAs) have been set up around the world over the past decades and collectively form the "International" PTA (IPTA). In this paper, we describe the first joint analysis of the data from the three regional PTAs, i.e. of the first IPTA data set. We describe the available PTA data, the approach presently followed for its combination and suggest improvements for future PTA research. Particular attention is paid to subtle details (such as underestimation of measurement uncertainty and long-period noise) that have often been ignored but which become important in this unprecedentedly large and inhomogeneous data set. We identify and describe in detail several factors that complicate IPTA research and provide recommendations for future pulsar timing efforts. The first IPTA data release presented here (and available online) is used to demonstrate the IPTA's potential of improving upon gravitational-wave limits placed by individual PTAs by a factor of ~2 and provides a 2-sigma limit on the dimensionless amplitude of a stochastic GWB of 1.7x10^{-15} at a frequency of 1 yr^{-1}. This is 1.7 times less constraining than the limit placed by (Shannon et al. 2015), due mostly to the more recent, high-quality data they used.
△ Less
Submitted 11 February, 2016;
originally announced February 2016.
-
Are we there yet? Time to detection of nanohertz gravitational waves based on pulsar-timing array limits
Authors:
S. R. Taylor,
M. Vallisneri,
J. A. Ellis,
C. M. F. Mingarelli,
T. J. W. Lazio,
R. van Haasteren
Abstract:
Decade-long timing observations of arrays of millisecond pulsars have placed highly constraining upper limits on the amplitude of the nanohertz gravitational-wave stochastic signal from the mergers of supermassive black-hole binaries ($\sim 10^{-15}$ strain at $f = 1/\mathrm{yr}$). These limits suggest that binary merger rates have been overestimated, or that environmental influences from nuclear…
▽ More
Decade-long timing observations of arrays of millisecond pulsars have placed highly constraining upper limits on the amplitude of the nanohertz gravitational-wave stochastic signal from the mergers of supermassive black-hole binaries ($\sim 10^{-15}$ strain at $f = 1/\mathrm{yr}$). These limits suggest that binary merger rates have been overestimated, or that environmental influences from nuclear gas or stars accelerate orbital decay, reducing the gravitational-wave signal at the lowest, most sensitive frequencies. This prompts the question whether nanohertz gravitational waves are likely to be detected in the near future. In this letter, we answer this question quantitatively using simple statistical estimates, deriving the range of true signal amplitudes that are compatible with current upper limits, and computing expected detection probabilities as a function of observation time. We conclude that small arrays consisting of the pulsars with the least timing noise, which yield the tightest upper limits, have discouraging prospects of making a detection in the next two decades. By contrast, we find large arrays are crucial to detection because the quadrupolar spatial correlations induced by gravitational waves can be well sampled by many pulsar pairs. Indeed, timing programs which monitor a large and expanding set of pulsars have an $\sim 80\%$ probability of detecting gravitational waves within the next ten years, under assumptions on merger rates and environmental influences ranging from optimistic to conservative. Even in the extreme case where $90\%$ of binaries stall before merger and environmental coupling effects diminish low-frequency gravitational-wave power, detection is delayed by at most a few years.
△ Less
Submitted 8 March, 2016; v1 submitted 17 November, 2015;
originally announced November 2015.
-
The noise properties of 42 millisecond pulsars from the European Pulsar Timing Array and their impact on gravitational wave searches
Authors:
R. N. Caballero,
K. J. Lee,
L. Lentati,
G. Desvignes,
D. J. Champion,
J. P. W. Verbiest,
G. H. Janssen,
B. W. Stappers,
M. Kramer,
P. Lazarus,
A. Possenti,
C. Tiburzi,
D. Perrodin,
S. Osłowski,
S. Babak,
C. G. Bassa,
P. Brem,
M. Burgay,
I. Cognard,
J. R. Gair,
E. Graikou,
L. Guillemot,
J. W. T. Hessels,
R. Karuppusamy,
A. Lassus
, et al. (14 additional authors not shown)
Abstract:
The sensitivity of Pulsar Timing Arrays to gravitational waves depends on the noise present in the individual pulsar timing data. Noise may be either intrinsic or extrinsic to the pulsar. Intrinsic sources of noise will include rotational instabilities, for example. Extrinsic sources of noise include contributions from physical processes which are not sufficiently well modelled, for example, dispe…
▽ More
The sensitivity of Pulsar Timing Arrays to gravitational waves depends on the noise present in the individual pulsar timing data. Noise may be either intrinsic or extrinsic to the pulsar. Intrinsic sources of noise will include rotational instabilities, for example. Extrinsic sources of noise include contributions from physical processes which are not sufficiently well modelled, for example, dispersion and scattering effects, analysis errors and instrumental instabilities. We present the results from a noise analysis for 42 millisecond pulsars (MSPs) observed with the European Pulsar Timing Array. For characterising the low-frequency, stochastic and achromatic noise component, or "timing noise", we employ two methods, based on Bayesian and frequentist statistics. For 25 MSPs, we achieve statistically significant measurements of their timing noise parameters and find that the two methods give consistent results. For the remaining 17 MSPs, we place upper limits on the timing noise amplitude at the 95% confidence level. We additionally place an upper limit on the contribution to the pulsar noise budget from errors in the reference terrestrial time standards (below 1%), and we find evidence for a noise component which is present only in the data of one of the four used telescopes. Finally, we estimate that the timing noise of individual pulsars reduces the sensitivity of this data set to an isotropic, stochastic GW background by a factor of >9.1 and by a factor of >2.3 for continuous GWs from resolvable, inspiralling supermassive black-hole binaries with circular orbits.
△ Less
Submitted 19 January, 2016; v1 submitted 30 October, 2015;
originally announced October 2015.
-
European Pulsar Timing Array Limits on Continuous Gravitational Waves from Individual Supermassive Black Hole Binaries
Authors:
Stanislav Babak,
Antoine Petiteau,
Alberto Sesana,
Patrick Brem,
Pablo A. Rosado,
Stephen R. Taylor,
Antoine Lassus,
Jason W. T. Hessels,
Cees G. Bassa,
Marta Burgay,
R. Nicolas Caballero,
David J. Champion,
Ismael Cognard,
Gregory Desvignes,
Jonathan R. Gair,
Lucas Guillemot,
Gemma H. Janssen,
Ramesh Karuppusamy,
Michael Kramer,
Patrick Lazarus,
K. J. Lee,
Lindley Lentati,
Kuo Liu,
Chiara M. F. Mingarelli,
Stefan Oslowski
, et al. (11 additional authors not shown)
Abstract:
We have searched for continuous gravitational wave (CGW) signals produced by individually resolvable, circular supermassive black hole binaries (SMBHBs) in the latest EPTA dataset, which consists of ultra-precise timing data on 41 millisecond pulsars. We develop frequentist and Bayesian detection algorithms to search both for monochromatic and frequency-evolving systems. None of the adopted algori…
▽ More
We have searched for continuous gravitational wave (CGW) signals produced by individually resolvable, circular supermassive black hole binaries (SMBHBs) in the latest EPTA dataset, which consists of ultra-precise timing data on 41 millisecond pulsars. We develop frequentist and Bayesian detection algorithms to search both for monochromatic and frequency-evolving systems. None of the adopted algorithms show evidence for the presence of such a CGW signal, indicating that the data are best described by pulsar and radiometer noise only. Depending on the adopted detection algorithm, the 95\% upper limit on the sky-averaged strain amplitude lies in the range $6\times 10^{-15}<A<1.5\times10^{-14}$ at $5{\rm nHz}<f<7{\rm nHz}$. This limit varies by a factor of five, depending on the assumed source position, and the most constraining limit is achieved towards the positions of the most sensitive pulsars in the timing array. The most robust upper limit -- obtained via a full Bayesian analysis searching simultaneously over the signal and pulsar noise on the subset of ours six best pulsars -- is $A\approx10^{-14}$. These limits, the most stringent to date at $f<10{\rm nHz}$, exclude the presence of sub-centiparsec binaries with chirp mass $\cal{M}_c>10^9$M$_\odot$ out to a distance of about 25Mpc, and with $\cal{M}_c>10^{10}$M$_\odot$ out to a distance of about 1Gpc ($z\approx0.2$). We show that state-of-the-art SMBHB population models predict $<1\%$ probability of detecting a CGW with the current EPTA dataset, consistent with the reported non-detection. We stress, however, that PTA limits on individual CGW have improved by almost an order of magnitude in the last five years. The continuing advances in pulsar timing data acquisition and analysis techniques will allow for strong astrophysical constraints on the population of nearby SMBHBs in the coming years.
△ Less
Submitted 11 September, 2015; v1 submitted 7 September, 2015;
originally announced September 2015.
-
The NANOGrav Nine-year Data Set: Limits on the Isotropic Stochastic Gravitational Wave Background
Authors:
Zaven Arzoumanian,
Adam Brazier,
Sarah Burke-Spolaor,
Sydney Chamberlin,
Shami Chatterjee,
Brian Christy,
Jim Cordes,
Neil Cornish,
Paul Demorest,
Xihao Deng,
Tim Dolch,
Justin Ellis,
Rob Ferdman,
Emmanuel Fonseca,
Nate Garver-Daniels,
Fredrick Jenet,
Glenn Jones,
Vicky Kaspi,
Michael Koop,
Michael Lam,
Joseph Lazio,
Lina Levin,
Andrea Lommen,
Duncan Lorimer,
Jin Luo
, et al. (23 additional authors not shown)
Abstract:
We compute upper limits on the nanohertz-frequency isotropic stochastic gravitational wave background (GWB) using the 9-year data release from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) collaboration. We set upper limits for a GWB from supermassive black hole binaries under power law, broken power law, and free spectral coefficient GW spectrum models. We place a 95…
▽ More
We compute upper limits on the nanohertz-frequency isotropic stochastic gravitational wave background (GWB) using the 9-year data release from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) collaboration. We set upper limits for a GWB from supermassive black hole binaries under power law, broken power law, and free spectral coefficient GW spectrum models. We place a 95\% upper limit on the strain amplitude (at a frequency of yr$^{-1}$) in the power law model of $A_{\rm gw} < 1.5\times 10^{-15}$. For a broken power law model, we place priors on the strain amplitude derived from simulations of Sesana (2013) and McWilliams et al. (2014). We find that the data favor a broken power law to a pure power law with odds ratios of 22 and 2.2 to one for the McWilliams and Sesana prior models, respectively. The McWilliams model is essentially ruled out by the data, and the Sesana model is in tension with the data under the assumption of a pure power law. Using the broken power-law analysis we construct posterior distributions on environmental factors that drive the binary to the GW-driven regime including the stellar mass density for stellar-scattering, mass accretion rate for circumbinary disk interaction, and orbital eccentricity for eccentric binaries, marking the first time that the shape of the GWB spectrum has been used to make astrophysical inferences. We then place the most stringent limits so far on the energy density of relic GWs, $Ω_\mathrm{gw}(f)\,h^2 < 4.2 \times 10^{-10}$, yielding a limit on the Hubble parameter during inflation of $H_*=1.6\times10^{-2}~m_{Pl}$, where $m_{Pl}$ is the Planck mass. Our limit on the cosmic string GWB, $Ω_\mathrm{gw}(f)\, h^2 < 2.2 \times 10^{-10}$, translates to a conservative limit of $Gμ<3.3\times 10^{-8}$ - a factor of 4 better than the joint Planck and high-$l$ CMB data from other experiments.
△ Less
Submitted 12 August, 2015;
originally announced August 2015.
-
Limits on anisotropy in the nanohertz stochastic gravitational-wave background
Authors:
S. R. Taylor,
C. M. F. Mingarelli,
J. R. Gair,
A. Sesana,
G. Theureau,
S. Babak,
C. G. Bassa,
P. Brem,
M. Burgay,
R. N. Caballero,
D. J. Champion,
I. Cognard,
G. Desvignes,
L. Guillemot,
J. W. T. Hessels,
G. H. Janssen,
R. Karuppusamy,
M. Kramer,
A. Lassus,
P. Lazarus,
L. Lentati,
K. Liu,
S. Osłowski,
D. Perrodin,
A. Petiteau
, et al. (10 additional authors not shown)
Abstract:
The paucity of observed supermassive black hole binaries (SMBHBs) may imply that the gravitational wave background (GWB) from this population is anisotropic, rendering existing analyses sub-optimal. We present the first constraints on the angular distribution of a nanohertz stochastic GWB from circular, inspiral-driven SMBHBs using the $2015$ European Pulsar Timing Array data [Desvignes et al. (in…
▽ More
The paucity of observed supermassive black hole binaries (SMBHBs) may imply that the gravitational wave background (GWB) from this population is anisotropic, rendering existing analyses sub-optimal. We present the first constraints on the angular distribution of a nanohertz stochastic GWB from circular, inspiral-driven SMBHBs using the $2015$ European Pulsar Timing Array data [Desvignes et al. (in prep.)]. Our analysis of the GWB in the $\sim 2 - 90$ nHz band shows consistency with isotropy, with the strain amplitude in $l>0$ spherical harmonic multipoles $\lesssim 40\%$ of the monopole value. We expect that these more general techniques will become standard tools to probe the angular distribution of source populations.
△ Less
Submitted 29 June, 2015;
originally announced June 2015.
-
The NANOGrav Nine-year Data Set: Observations, Arrival Time Measurements, and Analysis of 37 Millisecond Pulsars
Authors:
Z. Arzoumanian,
A. Brazier,
S. Burke-Spolaor,
S. Chamberlin,
S. Chatterjee,
B. Christy,
J. M. Cordes,
N. Cornish,
K. Crowter,
P. B. Demorest,
T. Dolch,
J. A. Ellis,
R. D. Ferdman,
E. Fonseca,
N. Garver-Daniels,
M. E. Gonzalez,
F. A. Jenet,
G. Jones,
M. Jones,
V. M. Kaspi,
M. Koop,
M. T. Lam,
T. J. W. Lazio,
L. Levin,
A. N. Lommen
, et al. (19 additional authors not shown)
Abstract:
We present high-precision timing observations spanning up to nine years for 37 millisecond pulsars monitored with the Green Bank and Arecibo radio telescopes as part of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) project. We describe the observational and instrumental setups used to collect the data, and methodology applied for calculating pulse times of arrival; th…
▽ More
We present high-precision timing observations spanning up to nine years for 37 millisecond pulsars monitored with the Green Bank and Arecibo radio telescopes as part of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) project. We describe the observational and instrumental setups used to collect the data, and methodology applied for calculating pulse times of arrival; these include novel methods for measuring instrumental offsets and characterizing low signal-to-noise ratio timing results. The time of arrival data are fit to a physical timing model for each source, including terms that characterize time-variable dispersion measure and frequency-dependent pulse shape evolution. In conjunction with the timing model fit, we have performed a Bayesian analysis of a parameterized timing noise model for each source, and detect evidence for excess low-frequency, or "red," timing noise in 10 of the pulsars. For 5 of these cases this is likely due to interstellar medium propagation effects rather than intrisic spin variations. Subsequent papers in this series will present further analysis of this data set aimed at detecting or limiting the presence of nanohertz-frequency gravitational wave signals.
△ Less
Submitted 5 January, 2016; v1 submitted 27 May, 2015;
originally announced May 2015.
-
European Pulsar Timing Array Limits On An Isotropic Stochastic Gravitational-Wave Background
Authors:
Lindley Lentati,
Stephen R. Taylor,
Chiara M. F. Mingarelli,
Alberto Sesana,
Sotiris A. Sanidas,
Alberto Vecchio,
R. Nicolas Caballero,
K. J. Lee,
Rutger van Haasteren,
Stanislav Babak,
Cees G. Bassa,
Patrick Brem,
Marta Burgay,
David J. Champion,
Ismael Cognard,
Gregory Desvignes,
Jonathan R. Gair,
Lucas Guillemot,
Jason W. T. Hessels,
Gemma H. Janssen,
Ramesh Karuppusamy,
Michael Kramer,
Antoine Lassus,
Patrick Lazarus,
Kuo Liu
, et al. (11 additional authors not shown)
Abstract:
We present new limits on an isotropic stochastic gravitational-wave background (GWB) using a six pulsar dataset spanning 18 yr of observations from the 2015 European Pulsar Timing Array data release. Performing a Bayesian analysis, we fit simultaneously for the intrinsic noise parameters for each pulsar, along with common correlated signals including clock, and Solar System ephemeris errors, obtai…
▽ More
We present new limits on an isotropic stochastic gravitational-wave background (GWB) using a six pulsar dataset spanning 18 yr of observations from the 2015 European Pulsar Timing Array data release. Performing a Bayesian analysis, we fit simultaneously for the intrinsic noise parameters for each pulsar, along with common correlated signals including clock, and Solar System ephemeris errors, obtaining a robust 95$\%$ upper limit on the dimensionless strain amplitude $A$ of the background of $A<3.0\times 10^{-15}$ at a reference frequency of $1\mathrm{yr^{-1}}$ and a spectral index of $13/3$, corresponding to a background from inspiralling super-massive black hole binaries, constraining the GW energy density to $Ω_\mathrm{gw}(f)h^2 < 1.1\times10^{-9}$ at 2.8 nHz. We also present limits on the correlated power spectrum at a series of discrete frequencies, and show that our sensitivity to a fiducial isotropic GWB is highest at a frequency of $\sim 5\times10^{-9}$~Hz. Finally we discuss the implications of our analysis for the astrophysics of supermassive black hole binaries, and present 95$\%$ upper limits on the string tension, $Gμ/c^2$, characterising a background produced by a cosmic string network for a set of possible scenarios, and for a stochastic relic GWB. For a Nambu-Goto field theory cosmic string network, we set a limit $Gμ/c^2<1.3\times10^{-7}$, identical to that set by the {\it Planck} Collaboration, when combining {\it Planck} and high-$\ell$ Cosmic Microwave Background data from other experiments. For a stochastic relic background we set a limit of $Ω^\mathrm{relic}_\mathrm{gw}(f)h^2<1.2 \times10^{-9}$, a factor of 9 improvement over the most stringent limits previously set by a pulsar timing array.
△ Less
Submitted 9 September, 2015; v1 submitted 14 April, 2015;
originally announced April 2015.
-
NANOGrav Constraints on Gravitational Wave Bursts with Memory
Authors:
Z. Arzoumanian,
A. Brazier,
S. Burke-Spolaor,
S. J. Chamberlin,
S. Chatterjee,
B. Christy,
J. M. Cordes,
N. J. Cornish,
P. B. Demorest,
X. Deng,
T. Dolch,
J. A. Ellis,
R. D. Ferdman,
E. Fonseca,
N. Garver-Daniels,
F. Jenet,
G. Jones,
V. M. Kaspi,
M. Koop,
M. T. Lam,
T. J. W. Lazio,
L. Levin,
A. N. Lommen,
D. R. Lorimer,
J. Luo
, et al. (17 additional authors not shown)
Abstract:
Among efforts to detect gravitational radiation, pulsar timing arrays are uniquely poised to detect "memory" signatures, permanent perturbations in spacetime from highly energetic astrophysical events such as mergers of supermassive black hole binaries. The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) observes dozens of the most stable millisecond pulsars using the Areci…
▽ More
Among efforts to detect gravitational radiation, pulsar timing arrays are uniquely poised to detect "memory" signatures, permanent perturbations in spacetime from highly energetic astrophysical events such as mergers of supermassive black hole binaries. The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) observes dozens of the most stable millisecond pulsars using the Arecibo and Green Bank radio telescopes in an effort to study, among other things, gravitational wave memory. We herein present the results of a search for gravitational wave bursts with memory (BWMs) using the first five years of NANOGrav observations. We develop original methods for dramatically speeding up searches for BWM signals. In the directions of the sky where our sensitivity to BWMs is best, we would detect mergers of binaries with reduced masses of $10^9$ $M_\odot$ out to distances of 30 Mpc; such massive mergers in the Virgo cluster would be marginally detectable. We find no evidence for BWMs. However, with our non-detection, we set upper limits on the rate at which BWMs of various amplitudes could have occurred during the time spanned by our data--e.g., BWMs with amplitudes greater than $10^{-13}$ must occur at a rate less than 1.5 yr$^{-1}$.
△ Less
Submitted 21 January, 2015;
originally announced January 2015.
-
Low-rank approximations for large stationary covariance matrices, as used in the Bayesian and generalized-least-squares analysis of pulsar-timing data
Authors:
Rutger van Haasteren,
Michele Vallisneri
Abstract:
Many data-analysis problems involve large dense matrices that describe the covariance of stationary noise processes; the computational cost of inverting these matrices, or equivalently of solving linear systems that contain them, is often a practical limit for the analysis. We describe two general, practical, and accurate methods to approximate stationary covariance matrices as low-rank matrix pro…
▽ More
Many data-analysis problems involve large dense matrices that describe the covariance of stationary noise processes; the computational cost of inverting these matrices, or equivalently of solving linear systems that contain them, is often a practical limit for the analysis. We describe two general, practical, and accurate methods to approximate stationary covariance matrices as low-rank matrix products featuring carefully chosen spectral components. These methods can be used to greatly accelerate data-analysis methods in many contexts, such as the Bayesian and generalized-least-squares analysis of pulsar-timing residuals.
△ Less
Submitted 24 July, 2014;
originally announced July 2014.
-
New advances in the Gaussian-process approach to pulsar-timing data analysis
Authors:
Rutger van Haasteren,
Michele Vallisneri
Abstract:
In this work we review the application of the theory of Gaussian processes to the modeling of noise in pulsar-timing data analysis, and we derive various useful and optimized representations for the likelihood expressions that are needed in Bayesian inference on pulsar-timing-array datasets. The resulting viewpoint and formalism lead us to two improved parameter-sampling schemes inspired by Gibbs…
▽ More
In this work we review the application of the theory of Gaussian processes to the modeling of noise in pulsar-timing data analysis, and we derive various useful and optimized representations for the likelihood expressions that are needed in Bayesian inference on pulsar-timing-array datasets. The resulting viewpoint and formalism lead us to two improved parameter-sampling schemes inspired by Gibbs sampling. The new schemes have vastly lower chain autocorrelation lengths than the Markov Chain Monte Carlo methods currently used in pulsar-timing data analysis, potentially speeding up Bayesian inference by orders of magnitude. The new schemes can be used for a full-noise-model analysis of the large datasets assembled by the International Pulsar Timing Array collaboration, which present a serious computational challenge to existing methods.
△ Less
Submitted 7 July, 2014;
originally announced July 2014.
-
Mapping the nano-Hertz gravitational wave sky
Authors:
Neil J. Cornish,
Rutger van Haasteren
Abstract:
We describe a new method for extracting gravitational wave signals from pulsar timing data. We show that any gravitational wave signal can be decomposed into an orthogonal set of sky maps, with the number of maps equal to the number of pulsars in the timing array. These maps may be used as a basis to construct gravitational wave templates for any type of source, including collections of point sour…
▽ More
We describe a new method for extracting gravitational wave signals from pulsar timing data. We show that any gravitational wave signal can be decomposed into an orthogonal set of sky maps, with the number of maps equal to the number of pulsars in the timing array. These maps may be used as a basis to construct gravitational wave templates for any type of source, including collections of point sources. A variant of the standard Hellings-Downs correlation analysis is recovered for statistically isotropic signals. The template based approach allows us to probe potential anisotropies in the signal and produce maps of the gravitational wave sky.
△ Less
Submitted 19 June, 2014; v1 submitted 17 June, 2014;
originally announced June 2014.
-
Model-based asymptotically optimal dispersion measure correction for pulsar timing
Authors:
K. J. Lee,
C. G. Bassa,
G. H. Janssen,
R. Karuppusamy,
M. Kramer,
K. Liu,
D. Perrodin,
R. Smits,
B. W. Stappers,
R. van Haasteren,
L. Lentati
Abstract:
In order to reach the sensitivity required to detect gravitational waves, pulsar timing array experiments need to mitigate as much noise as possible in timing data. A dominant amount of noise is likely due to variations in the dispersion measure. To correct for such variations, we develop a statistical method inspired by the maximum likelihood estimator and optimal filtering. Our method consists o…
▽ More
In order to reach the sensitivity required to detect gravitational waves, pulsar timing array experiments need to mitigate as much noise as possible in timing data. A dominant amount of noise is likely due to variations in the dispersion measure. To correct for such variations, we develop a statistical method inspired by the maximum likelihood estimator and optimal filtering. Our method consists of two major steps. First, the spectral index and amplitude of dispersion measure variations are measured via a time-domain spectral analysis. Second, the linear optimal filter is constructed based on the model parameters found in the first step, and is used to extract the dispersion measure variation waveforms. Compared to current existing methods, this method has better time resolution for the study of short timescale dispersion variations, and generally produces smaller errors in waveform estimations. This method can process irregularly sampled data without any interpolation because of its time-domain nature. Furthermore, it offers the possibility to interpolate or extrapolate the waveform estimation to regions where no data is available. Examples using simulated data sets are included for demonstration.
△ Less
Submitted 8 April, 2014;
originally announced April 2014.
-
NANOGrav Limits on Gravitational Waves from Individual Supermassive Black Hole Binaries in Circular Orbits
Authors:
Z. Arzoumanian,
A. Brazier,
S. Burke-Spolaor,
S. J. Chamberlin,
S. Chatterjee,
J. M. Cordes,
P. B. Demorest,
X. Deng,
T. Dolch,
J. A. Ellis,
R. D. Ferdman,
N. Garver-Daniels,
F. Jenet,
G. Jones,
V. M. Kaspi,
M. Koop,
M. Lam,
T. J. W. Lazio,
A. N. Lommen,
D. R. Lorimer,
J. Luo,
R. S. Lynch,
D. R. Madison,
M. McLaughlin,
S. T. McWilliams
, et al. (14 additional authors not shown)
Abstract:
The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) project currently observes 43 pulsars using the Green Bank and Arecibo radio telescopes. In this work we use a subset of 17 pulsars timed for a span of roughly five years (2005--2010). We analyze these data using standard pulsar timing models, with the addition of time-variable dispersion measure and frequency-variable pul…
▽ More
The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) project currently observes 43 pulsars using the Green Bank and Arecibo radio telescopes. In this work we use a subset of 17 pulsars timed for a span of roughly five years (2005--2010). We analyze these data using standard pulsar timing models, with the addition of time-variable dispersion measure and frequency-variable pulse shape terms. Within the timing data, we perform a search for continuous gravitational waves from individual supermassive black hole binaries in circular orbits using robust frequentist and Bayesian techniques. We find that there is no evidence for the presence of a detectable continuous gravitational wave; however, we can use these data to place the most constraining upper limits to date on the strength of such gravitational waves. Using the full 17 pulsar dataset we place a 95% upper limit on the sky-averaged strain amplitude of $h_0\lesssim 3.8\times 10^{-14}$ at a frequency of 10 nHz. Furthermore, we place 95% \emph{all sky} lower limits on the luminosity distance to such gravitational wave sources finding that the $d_L \gtrsim 425$ Mpc for sources at a frequency of 10 nHz and chirp mass $10^{10}{\rm M}_{\odot}$. We find that for gravitational wave sources near our best timed pulsars in the sky, the sensitivity of the pulsar timing array is increased by a factor of $\sim$4 over the sky-averaged sensitivity. Finally we place limits on the coalescence rate of the most massive supermassive black hole binaries.
△ Less
Submitted 19 May, 2014; v1 submitted 4 April, 2014;
originally announced April 2014.
-
TempoNest: A Bayesian approach to pulsar timing analysis
Authors:
Lindley Lentati,
Paul Alexander,
Michael P. Hobson,
Farhan Feroz,
Rutger van Haasteren,
Kejia Lee,
Ryan M. Shannon
Abstract:
A new Bayesian software package for the analysis of pulsar timing data is presented in the form of TempoNest which allows for the robust determination of the non-linear pulsar timing solution simultaneously with a range of additional stochastic parameters. This includes both red spin noise and dispersion measure variations using either power law descriptions of the noise, or through a model-indepe…
▽ More
A new Bayesian software package for the analysis of pulsar timing data is presented in the form of TempoNest which allows for the robust determination of the non-linear pulsar timing solution simultaneously with a range of additional stochastic parameters. This includes both red spin noise and dispersion measure variations using either power law descriptions of the noise, or through a model-independent method that parameterises the power at individual frequencies in the signal. We use TempoNest to show that at noise levels representative of current datasets in the European Pulsar Timing Array (EPTA) and International Pulsar Timing Array (IPTA) the linear timing model can underestimate the uncertainties of the timing solution by up to an order of magnitude. We also show how to perform Bayesian model selection between different sets of timing model and stochastic parameters, for example, by demonstrating that in the pulsar B1937+21 both the dispersion measure variations and spin noise in the data are optimally modelled by simple power laws. Finally we show that not including the stochastic parameters simultaneously with the timing model can lead to unpredictable variation in the estimated uncertainties, compromising the robustness of the scientific results extracted from such analysis.
△ Less
Submitted 8 October, 2013;
originally announced October 2013.
-
An Efficient Approximation to the Likelihood for Gravitational Wave Stochastic Background Detection Using Pulsar Timing Data
Authors:
Justin Ellis,
Xavier Siemens,
Rutger van Haasteren
Abstract:
Direct detection of gravitational waves by pulsar timing arrays will become feasible over the next few years. In the low frequency regime ($10^{-7}$ Hz -- $10^{-9}$ Hz), we expect that a superposition of gravitational waves from many sources will manifest itself as an isotropic stochastic gravitational wave background. Currently, a number of techniques exist to detect such a signal; however, many…
▽ More
Direct detection of gravitational waves by pulsar timing arrays will become feasible over the next few years. In the low frequency regime ($10^{-7}$ Hz -- $10^{-9}$ Hz), we expect that a superposition of gravitational waves from many sources will manifest itself as an isotropic stochastic gravitational wave background. Currently, a number of techniques exist to detect such a signal; however, many detection methods are computationally challenging. Here we introduce an approximation to the full likelihood function for a pulsar timing array that results in computational savings proportional to the square of the number of pulsars in the array. Through a series of simulations we show that the approximate likelihood function reproduces results obtained from the full likelihood function. We further show, both analytically and through simulations, that, on average, this approximate likelihood function gives unbiased parameter estimates for astrophysically realistic stochastic background amplitudes.
△ Less
Submitted 7 February, 2013;
originally announced February 2013.
-
Analysis of the first IPTA Mock Data Challenge by the EPTA timing data analysis working group
Authors:
Rutger van Haasteren,
Chiara M. F. Mingarelli,
Alberto Vecchio,
Antoine Lassus
Abstract:
This is a summary of the methods we used to analyse the first IPTA Mock Data Challenge (MDC), and the obtained results. We have used a Bayesian analysis in the time domain, accelerated using the recently developed ABC-method which consists of a form of lossy linear data compression. The TOAs were first processed with Tempo2, where the design matrix was extracted for use in a subsequent Bayesian an…
▽ More
This is a summary of the methods we used to analyse the first IPTA Mock Data Challenge (MDC), and the obtained results. We have used a Bayesian analysis in the time domain, accelerated using the recently developed ABC-method which consists of a form of lossy linear data compression. The TOAs were first processed with Tempo2, where the design matrix was extracted for use in a subsequent Bayesian analysis. We used different noise models to analyse the datasets: no red noise, red noise the same for all pulsars, and individual red noise per pulsar. We sampled from the likelihood with four different samplers: "emcee", "t-walk", "Metropolis-Hastings", and "pyMultiNest". All but emcee agreed on the final result, with emcee failing due to artefacts of the high-dimensionality of the problem. An interesting issue we ran into was that the prior of all the 36 (red) noise amplitudes strongly affects the results. A flat prior in the noise amplitude biases the inferred GWB amplitude, whereas a flat prior in log-amplitude seems to work well. This issue is only apparent when using a noise model with individually modelled red noise for all pulsars. Our results for the blind challenges are in good agreement with the injected values. For the GWB amplitudes we found h_c = 1.03 +/- 0.11 [10^{-14}], h_c = 5.70 +/- 0.35 [10^{-14}], and h_c = 6.91 +/- 1.72 [10^{-15}], and for the GWB spectral index we found gamma = 4.28 +/- 0.20, gamma = 4.35 +/- 0.09, and gamma = 3.75 +/- 0.40. We note that for closed challenge 3 there was quite some covariance between the signal and the red noise: if we constrain the GWB spectral index to the usual choice of gamma = 13/3, we obtain the estimates: h_c = 10.0 +/- 0.64 [10^{-15}], h_c = 56.3 +/- 2.42 [10^{-15}], and h_c = 4.83 +/- 0.50 [10^{-15}], with one-sided 2 sigma upper-limits of: h_c <= 10.98 [10^{-15}], h_c <= 60.29 [10^{-15}], and h_c <= 5.65 [10^{-15}].
△ Less
Submitted 28 January, 2013;
originally announced January 2013.
-
Hyper-efficient model-independent Bayesian method for the analysis of pulsar timing data
Authors:
Lindley Lentati,
Paul Alexander,
Michael P. Hobson,
Stephen Taylor,
Jonathon Gair,
Sreekumar T. Balan,
Rutger van Haasteren
Abstract:
A new model independent method is presented for the analysis of pulsar timing data and the estimation of the spectral properties of an isotropic gravitational wave background (GWB). We show that by rephrasing the likelihood we are able to eliminate the most costly aspects of computation normally associated with this type of data analysis. When applied to the International Pulsar Timing Array Mock…
▽ More
A new model independent method is presented for the analysis of pulsar timing data and the estimation of the spectral properties of an isotropic gravitational wave background (GWB). We show that by rephrasing the likelihood we are able to eliminate the most costly aspects of computation normally associated with this type of data analysis. When applied to the International Pulsar Timing Array Mock Data Challenge data sets this results in speedups of approximately 2 to 3 orders of magnitude compared to established methods. We present three applications of the new likelihood. In the low signal to noise regime we sample directly from the power spectrum coefficients of the GWB signal realization. In the high signal to noise regime, where the data can support a large number of coefficients, we sample from the joint probability density of the power spectrum coefficients for the individual pulsars and the GWB signal realization. Critically in both these cases we need make no assumptions about the form of the power spectrum of the GWB, or the individual pulsars. Finally we present a method for characterizing the spatial correlation between pulsars on the sky, making no assumptions about the form of that correlation, and therefore providing the only truly general Bayesian method of confirming a GWB detection from pulsar timing data.
△ Less
Submitted 25 June, 2013; v1 submitted 12 October, 2012;
originally announced October 2012.
-
Accelerating pulsar timing data analysis
Authors:
Rutger van Haasteren
Abstract:
The analysis of pulsar timing data, especially in pulsar timing array (PTA) projects, has encountered practical difficulties: evaluating the likelihood and/or correlation-based statistics can become prohibitively computationally expensive for large datasets. In situations where a stochastic signal of interest has a power spectral density that dominates the noise in a limited bandwidth of the total…
▽ More
The analysis of pulsar timing data, especially in pulsar timing array (PTA) projects, has encountered practical difficulties: evaluating the likelihood and/or correlation-based statistics can become prohibitively computationally expensive for large datasets. In situations where a stochastic signal of interest has a power spectral density that dominates the noise in a limited bandwidth of the total frequency domain (e.g. the isotropic background of gravitational waves), a linear transformation exists that transforms the timing residuals to a basis in which virtually all the information about the stochastic signal of interest is contained in a small fraction of basis vectors. By only considering such a small subset of these "generalised residuals", the dimensionality of the data analysis problem is greatly reduced, which can cause a large speedup in the evaluation of the likelihood: the ABC-method (Acceleration By Compression). The compression fidelity, calculable with crude estimates of the signal and noise, can be used to determine how far a dataset can be compressed without significant loss of information. Both direct tests on the likelihood, and Bayesian analysis of mock data, show that the signal can be recovered as well as with an analysis of uncompressed data. In the analysis of IPTA Mock Data Challenge datasets, speedups of a factor of three orders of magnitude are demonstrated. For realistic PTA datasets the acceleration may become greater than six orders of magnitude due to the low signal to noise ratio.
△ Less
Submitted 15 November, 2012; v1 submitted 1 October, 2012;
originally announced October 2012.
-
Understanding and analysing time-correlated stochastic signals in pulsar timing
Authors:
Rutger van Haasteren,
Yuri Levin
Abstract:
Although it is widely understood that pulsar timing observations generally contain time-correlated stochastic signals (TCSSs; red timing noise is of this type), most data analysis techniques that have been developed make an assumption that the stochastic uncertainties in the data are uncorrelated, i.e. "white". Recent work has pointed out that this can introduce severe bias in determination of tim…
▽ More
Although it is widely understood that pulsar timing observations generally contain time-correlated stochastic signals (TCSSs; red timing noise is of this type), most data analysis techniques that have been developed make an assumption that the stochastic uncertainties in the data are uncorrelated, i.e. "white". Recent work has pointed out that this can introduce severe bias in determination of timing-model parameters, and that better analysis methods should be used. This paper presents a detailed investigation of timing-model fitting in the presence of TCSSs, and gives closed expressions for the post-fit signals in the data. This results in a Bayesian technique to obtain timing-model parameter estimates in the presence of TCSSs, as well as computationally more efficient expressions of their marginalised posterior distribution. A new method to analyse hundreds of mock dataset realisations simultaneously without significant computational overhead is presented, as well as a statistically rigorous method to check the internal consistency of the results. As a by-product of the analysis, closed expressions of the rms introduced by a stochastic background of gravitational-waves in timing-residuals are obtained. Using $T$ as the length of the dataset, and $h_c(1\rm{yr}^{-1})$ as the characteristic strain, this is: $σ_{\rm GWB}^2 = h_{c}(1\rm{yr}^{-1})^2 (9\sqrt[3]{2π^4}Γ(-10/3) / 8008) \rm{yr}^{-4/3} T^{10/3}$.
△ Less
Submitted 2 October, 2012; v1 submitted 27 February, 2012;
originally announced February 2012.
-
Placing limits on the stochastic gravitational-wave background using European Pulsar Timing Array data
Authors:
R. van Haasteren,
Y. Levin,
G. H. Janssen,
K. Lazaridis,
M. Kramer B. W. Stappers,
G. Desvignes,
M. B. Purver,
A. G. Lyne,
R. D. Ferdman,
A. Jessner,
I. Cognard,
G. Theureau,
N. D'Amico,
A. Possenti,
M. Burgay,
A. Corongiu,
J. W. T. Hessels,
R. Smits,
J. P. W. Verbiest
Abstract:
Direct detection of low-frequency gravitational waves ($10^{-9} - 10^{-8}$ Hz) is the main goal of pulsar timing array (PTA) projects. One of the main targets for the PTAs is to measure the stochastic background of gravitational waves (GWB) whose characteristic strain is expected to approximately follow a power-law of the form $h_c(f)=A (f/\hbox{yr}^{-1})^α$, where $f$ is the gravitational-wave fr…
▽ More
Direct detection of low-frequency gravitational waves ($10^{-9} - 10^{-8}$ Hz) is the main goal of pulsar timing array (PTA) projects. One of the main targets for the PTAs is to measure the stochastic background of gravitational waves (GWB) whose characteristic strain is expected to approximately follow a power-law of the form $h_c(f)=A (f/\hbox{yr}^{-1})^α$, where $f$ is the gravitational-wave frequency. In this paper we use the current data from the European PTA to determine an upper limit on the GWB amplitude $A$ as a function of the unknown spectral slope $α$ with a Bayesian algorithm, by modelling the GWB as a random Gaussian process. For the case $α=-2/3$, which is expected if the GWB is produced by supermassive black-hole binaries, we obtain a 95% confidence upper limit on $A$ of $6\times 10^{-15}$, which is 1.8 times lower than the 95% confidence GWB limit obtained by the Parkes PTA in 2006. Our approach to the data analysis incorporates the multi-telescope nature of the European PTA and thus can serve as a useful template for future intercontinental PTA collaborations.
△ Less
Submitted 2 March, 2011;
originally announced March 2011.
-
The European Pulsar Timing Array: current efforts and a LEAP toward the future
Authors:
Robert D. Ferdman,
Rutger van Haasteren,
Cees G. Bassa,
Marta Burgay,
Ismael Cognard,
Alessandro Corongiu,
Nichi D'Amico,
Gregory Desvignes,
Jason W. T. Hessels,
Gemma H. Janssen,
Axel Jessner,
Christine Jordan,
Ramesh Karuppusamy,
Evan F. Keane,
Michael Kramer,
Kosmas Lazaridis,
Yuri Levin,
Andrew G. Lyne,
Maura Pilia,
Andrea Possenti,
Mark Purver,
Ben Stappers,
Sotirios Sanidas,
Roy Smits,
Gilles Theureau
Abstract:
The European Pulsar Timing Array (EPTA) is a multi-institutional, multi-telescope collaboration, with the goal of using high-precision pulsar timing to directly detect gravitational waves. In this article we discuss the EPTA member telescopes, current achieved timing precision, and near-future goals. We report a preliminary upper limit to the amplitude of a gravitational wave background. We also d…
▽ More
The European Pulsar Timing Array (EPTA) is a multi-institutional, multi-telescope collaboration, with the goal of using high-precision pulsar timing to directly detect gravitational waves. In this article we discuss the EPTA member telescopes, current achieved timing precision, and near-future goals. We report a preliminary upper limit to the amplitude of a gravitational wave background. We also discuss the Large European Array for Pulsars, in which the five major European telescopes involved in pulsar timing will be combined to provide a coherent array that will give similar sensitivity to the Arecibo radio telescope, and larger sky coverage.
△ Less
Submitted 17 March, 2010;
originally announced March 2010.
-
The international pulsar timing array project: using pulsars as a gravitational wave detector
Authors:
G. Hobbs,
A. Archibald,
Z. Arzoumanian,
D. Backer,
M. Bailes,
N. D. R. Bhat,
M. Burgay,
S. Burke-Spolaor,
D. Champion,
I. Cognard,
W. Coles,
J. Cordes,
P. Demorest,
G. Desvignes,
R. D. Ferdman,
L. Finn,
P. Freire,
M. Gonzalez,
J. Hessels,
A. Hotan,
G. Janssen,
F. Jenet,
A. Jessner,
C. Jordan,
V. Kaspi
, et al. (33 additional authors not shown)
Abstract:
The International Pulsar Timing Array project combines observations of pulsars from both Northern and Southern hemisphere observatories with the main aim of detecting ultra-low frequency (~10^-9 to 10^-8 Hz) gravitational waves. Here we introduce the project, review the methods used to search for gravitational waves emitted from coalescing supermassive binary black-hole systems in the centres of…
▽ More
The International Pulsar Timing Array project combines observations of pulsars from both Northern and Southern hemisphere observatories with the main aim of detecting ultra-low frequency (~10^-9 to 10^-8 Hz) gravitational waves. Here we introduce the project, review the methods used to search for gravitational waves emitted from coalescing supermassive binary black-hole systems in the centres of merging galaxies and discuss the status of the project.
△ Less
Submitted 27 November, 2009;
originally announced November 2009.
-
Bayesian evidence: can we beat MultiNest using traditional MCMC methods?
Authors:
Rutger van Haasteren
Abstract:
Markov Chain Monte Carlo (MCMC) methods have revolutionised Bayesian data analysis over the years by making the direct computation of posterior probability densities feasible on modern workstations. However, the calculation of the prior predictive, the Bayesian evidence, has proved to be notoriously difficult with standard techniques. In this work a method is presented that lets one calculate th…
▽ More
Markov Chain Monte Carlo (MCMC) methods have revolutionised Bayesian data analysis over the years by making the direct computation of posterior probability densities feasible on modern workstations. However, the calculation of the prior predictive, the Bayesian evidence, has proved to be notoriously difficult with standard techniques. In this work a method is presented that lets one calculate the Bayesian evidence using nothing but the results from standard MCMC algorithms, like Metropolis-Hastings. This new method is compared to other methods like MultiNest, and greatly outperforms the latter in several cases. One of the toy problems considered in this work is the analysis of mock pulsar timing data, as encountered in pulsar timing array projects. This method is expected to be useful as well in other problems in astrophysics, cosmology and particle physics.
△ Less
Submitted 11 November, 2009;
originally announced November 2009.
-
Gravitational-wave memory and pulsar timing arrays
Authors:
Rutger van Haasteren,
Yuri Levin
Abstract:
Pulsar timing arrays (PTAs) are designed to detect gravitational waves with periods from several months to several years, e.g. those produced by by wide supermassive black-hole binaries in the centers of distant galaxies. Here we show that PTAs are also sensitive to mergers of supermassive black holes. While these mergers occur on a timescale too short to be resolvable by a PTA, they generate a…
▽ More
Pulsar timing arrays (PTAs) are designed to detect gravitational waves with periods from several months to several years, e.g. those produced by by wide supermassive black-hole binaries in the centers of distant galaxies. Here we show that PTAs are also sensitive to mergers of supermassive black holes. While these mergers occur on a timescale too short to be resolvable by a PTA, they generate a change of metric due to non-linear gravitational-wave memory which persists for the duration of the experiment and could be detected. We develop the theory of the single-source detection by PTAs, and derive the sensitivity of PTAs to the gravitational-wave memory jumps. We show that mergers of $10^8M_{\odot}$ black holes are $2-σ$-detectable (in a direction, polarization, and time-dependent way) out to co-moving distances of $\sim 1$ billion light years. Modern prediction for black-hole merger rates imply marginal to modest chance of an individual jump detection by currently developed PTAs. The sensitivity is expected to be somewhat higher for futuristic PTA experiments with SKA.
△ Less
Submitted 4 September, 2009;
originally announced September 2009.
-
On measuring the gravitational-wave background using Pulsar Timing Arrays
Authors:
Rutger van Haasteren,
Yuri Levin,
Patrick McDonald,
Tingting Lu
Abstract:
Long-term precise timing of Galactic millisecond pulsars holds great promise for measuring the long-period (months-to-years) astrophysical gravitational waves. Several gravitational-wave observational programs, called Pulsar Timing Arrays (PTA), are being pursued around the world. Here we develop a Bayesian algorithm for measuring the stochastic gravitational-wave background (GWB) from the PTA d…
▽ More
Long-term precise timing of Galactic millisecond pulsars holds great promise for measuring the long-period (months-to-years) astrophysical gravitational waves. Several gravitational-wave observational programs, called Pulsar Timing Arrays (PTA), are being pursued around the world. Here we develop a Bayesian algorithm for measuring the stochastic gravitational-wave background (GWB) from the PTA data. Our algorithm has several strengths: (1) It analyses the data without any loss of information, (2) It trivially removes systematic errors of known functional form, including quadratic pulsar spin-down, annual modulations and jumps due to a change of equipment, (3) It measures simultaneously both the amplitude and the slope of the GWB spectrum, (4) It can deal with unevenly sampled data and coloured pulsar noise spectra. We sample the likelihood function using Markov Chain Monte Carlo (MCMC) simulations. We extensively test our approach on mock PTA datasets, and find that the algorithm has significant benefits over currently proposed counterparts. We show the importance of characterising all red noise components in pulsar timing noise by demonstrating that the presence of a red component would significantly hinder a detection of the GWB. Lastly, we explore the dependence of the signal-to-noise ratio on the duration of the experiment, number of monitored pulsars, and the magnitude of the pulsar timing noise. These parameter studies will help formulate observing strategies for the PTA experiments.
△ Less
Submitted 11 February, 2009; v1 submitted 4 September, 2008;
originally announced September 2008.