-
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 Data Set: Search for Anisotropy in the Gravitational-Wave Background
Authors:
Gabriella Agazie,
Akash Anumarlapudi,
Anne M. Archibald,
Zaven Arzoumanian,
Paul T. Baker,
Bence Bécsy,
Laura Blecha,
Adam Brazier,
Paul R. Brook,
Sarah Burke-Spolaor,
J. Andrew Casey-Clyde,
Maria Charisi,
Shami Chatterjee,
Tyler Cohen,
James M. Cordes,
Neil J. Cornish,
Fronefield Crawford,
H. Thankful Cromartie,
Kathryn Crowter,
Megan E. DeCesar,
Paul B. Demorest,
Timothy Dolch,
Brendan Drachler,
Elizabeth C. Ferrara,
William Fiore
, et al. (68 additional authors not shown)
Abstract:
The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) has reported evidence for the presence of an isotropic nanohertz gravitational wave background (GWB) in its 15 yr dataset. However, if the GWB is produced by a population of inspiraling supermassive black hole binary (SMBHB) systems, then the background is predicted to be anisotropic, depending on the distribution of these…
▽ More
The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) has reported evidence for the presence of an isotropic nanohertz gravitational wave background (GWB) in its 15 yr dataset. However, if the GWB is produced by a population of inspiraling supermassive black hole binary (SMBHB) systems, then the background is predicted to be anisotropic, depending on the distribution of these systems in the local Universe and the statistical properties of the SMBHB population. In this work, we search for anisotropy in the GWB using multiple methods and bases to describe the distribution of the GWB power on the sky. We do not find significant evidence of anisotropy, and place a Bayesian $95\%$ upper limit on the level of broadband anisotropy such that $(C_{l>0} / C_{l=0}) < 20\%$. We also derive conservative estimates on the anisotropy expected from a random distribution of SMBHB systems using astrophysical simulations conditioned on the isotropic GWB inferred in the 15-yr dataset, and show that this dataset has sufficient sensitivity to probe a large fraction of the predicted level of anisotropy. We end by highlighting the opportunities and challenges in searching for anisotropy in pulsar timing array data.
△ Less
Submitted 28 June, 2023;
originally announced June 2023.
-
The NANOGrav 15-year Data Set: Constraints on Supermassive Black Hole Binaries from the Gravitational Wave Background
Authors:
Gabriella Agazie,
Akash Anumarlapudi,
Anne M. Archibald,
Paul T. Baker,
Bence Bécsy,
Laura Blecha,
Alexander Bonilla,
Adam Brazier,
Paul R. Brook,
Sarah Burke-Spolaor,
Rand Burnette,
Robin Case,
J. Andrew Casey-Clyde,
Maria Charisi,
Shami Chatterjee,
Katerina Chatziioannou,
Belinda D. Cheeseboro,
Siyuan Chen,
Tyler Cohen,
James M. Cordes,
Neil J. Cornish,
Fronefield Crawford,
H. Thankful Cromartie,
Kathryn Crowter,
Curt J. Cutler
, et al. (89 additional authors not shown)
Abstract:
The NANOGrav 15-year data set shows evidence for the presence of a low-frequency gravitational-wave background (GWB). While many physical processes can source such low-frequency gravitational waves, here we analyze the signal as coming from a population of supermassive black hole (SMBH) binaries distributed throughout the Universe. We show that astrophysically motivated models of SMBH binary popul…
▽ More
The NANOGrav 15-year data set shows evidence for the presence of a low-frequency gravitational-wave background (GWB). While many physical processes can source such low-frequency gravitational waves, here we analyze the signal as coming from a population of supermassive black hole (SMBH) binaries distributed throughout the Universe. We show that astrophysically motivated models of SMBH binary populations are able to reproduce both the amplitude and shape of the observed low-frequency gravitational-wave spectrum. While multiple model variations are able to reproduce the GWB spectrum at our current measurement precision, our results highlight the importance of accurately modeling binary evolution for producing realistic GWB spectra. Additionally, while reasonable parameters are able to reproduce the 15-year observations, the implied GWB amplitude necessitates either a large number of parameters to be at the edges of expected values, or a small number of parameters to be notably different from standard expectations. While we are not yet able to definitively establish the origin of the inferred GWB signal, the consistency of the signal with astrophysical expectations offers a tantalizing prospect for confirming that SMBH binaries are able to form, reach sub-parsec separations, and eventually coalesce. As the significance grows over time, higher-order features of the GWB spectrum will definitively determine the nature of the GWB and allow for novel constraints on SMBH populations.
△ Less
Submitted 18 July, 2023; v1 submitted 28 June, 2023;
originally announced June 2023.
-
The NANOGrav 15-year Data Set: Search for Signals from New Physics
Authors:
Adeela Afzal,
Gabriella Agazie,
Akash Anumarlapudi,
Anne M. Archibald,
Zaven Arzoumanian,
Paul T. Baker,
Bence Bécsy,
Jose Juan Blanco-Pillado,
Laura Blecha,
Kimberly K. Boddy,
Adam Brazier,
Paul R. Brook,
Sarah Burke-Spolaor,
Rand Burnette,
Robin Case,
Maria Charisi,
Shami Chatterjee,
Katerina Chatziioannou,
Belinda D. Cheeseboro,
Siyuan Chen,
Tyler Cohen,
James M. Cordes,
Neil J. Cornish,
Fronefield Crawford,
H. Thankful Cromartie
, et al. (98 additional authors not shown)
Abstract:
The 15-year pulsar timing data set collected by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) shows positive evidence for the presence of a low-frequency gravitational-wave (GW) background. In this paper, we investigate potential cosmological interpretations of this signal, specifically cosmic inflation, scalar-induced GWs, first-order phase transitions, cosmic string…
▽ More
The 15-year pulsar timing data set collected by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) shows positive evidence for the presence of a low-frequency gravitational-wave (GW) background. In this paper, we investigate potential cosmological interpretations of this signal, specifically cosmic inflation, scalar-induced GWs, first-order phase transitions, cosmic strings, and domain walls. We find that, with the exception of stable cosmic strings of field theory origin, all these models can reproduce the observed signal. When compared to the standard interpretation in terms of inspiraling supermassive black hole binaries (SMBHBs), many cosmological models seem to provide a better fit resulting in Bayes factors in the range from 10 to 100. However, these results strongly depend on modeling assumptions about the cosmic SMBHB population and, at this stage, should not be regarded as evidence for new physics. Furthermore, we identify excluded parameter regions where the predicted GW signal from cosmological sources significantly exceeds the NANOGrav signal. These parameter constraints are independent of the origin of the NANOGrav signal and illustrate how pulsar timing data provide a new way to constrain the parameter space of these models. Finally, we search for deterministic signals produced by models of ultralight dark matter (ULDM) and dark matter substructures in the Milky Way. We find no evidence for either of these signals and thus report updated constraints on these models. In the case of ULDM, these constraints outperform torsion balance and atomic clock constraints for ULDM coupled to electrons, muons, or gluons.
△ Less
Submitted 28 June, 2023;
originally announced June 2023.
-
The NANOGrav 15-year Data Set: Evidence for a Gravitational-Wave Background
Authors:
Gabriella Agazie,
Akash Anumarlapudi,
Anne M. Archibald,
Zaven Arzoumanian,
Paul T. Baker,
Bence Becsy,
Laura Blecha,
Adam Brazier,
Paul R. Brook,
Sarah Burke-Spolaor,
Rand Burnette,
Robin Case,
Maria Charisi,
Shami Chatterjee,
Katerina Chatziioannou,
Belinda D. Cheeseboro,
Siyuan Chen,
Tyler Cohen,
James M. Cordes,
Neil J. Cornish,
Fronefield Crawford,
H. Thankful Cromartie,
Kathryn Crowter,
Curt J. Cutler,
Megan E. DeCesar
, et al. (89 additional authors not shown)
Abstract:
We report multiple lines of evidence for a stochastic signal that is correlated among 67 pulsars from the 15-year pulsar-timing data set collected by the North American Nanohertz Observatory for Gravitational Waves. The correlations follow the Hellings-Downs pattern expected for a stochastic gravitational-wave background. The presence of such a gravitational-wave background with a power-law-spectr…
▽ More
We report multiple lines of evidence for a stochastic signal that is correlated among 67 pulsars from the 15-year pulsar-timing data set collected by the North American Nanohertz Observatory for Gravitational Waves. The correlations follow the Hellings-Downs pattern expected for a stochastic gravitational-wave background. The presence of such a gravitational-wave background with a power-law-spectrum is favored over a model with only independent pulsar noises with a Bayes factor in excess of $10^{14}$, and this same model is favored over an uncorrelated common power-law-spectrum model with Bayes factors of 200-1000, depending on spectral modeling choices. We have built a statistical background distribution for these latter Bayes factors using a method that removes inter-pulsar correlations from our data set, finding $p = 10^{-3}$ (approx. $3σ$) for the observed Bayes factors in the null no-correlation scenario. A frequentist test statistic built directly as a weighted sum of inter-pulsar correlations yields $p = 5 \times 10^{-5} - 1.9 \times 10^{-4}$ (approx. $3.5 - 4σ$). Assuming a fiducial $f^{-2/3}$ characteristic-strain spectrum, as appropriate for an ensemble of binary supermassive black-hole inspirals, the strain amplitude is $2.4^{+0.7}_{-0.6} \times 10^{-15}$ (median + 90% credible interval) at a reference frequency of 1/(1 yr). The inferred gravitational-wave background amplitude and spectrum are consistent with astrophysical expectations for a signal from a population of supermassive black-hole binaries, although more exotic cosmological and astrophysical sources cannot be excluded. The observation of Hellings-Downs correlations points to the gravitational-wave origin of this signal.
△ Less
Submitted 28 June, 2023;
originally announced June 2023.
-
Searching for continuous Gravitational Waves in the second data release of the International Pulsar Timing Array
Authors:
M. Falxa,
S. Babak,
P. T. Baker,
B. Bécsy,
A. Chalumeau,
S. Chen,
Z. Chen,
N. J. Cornish,
L. Guillemot,
J. S. Hazboun,
C. M. F. Mingarelli,
A. Parthasarathy,
A. Petiteau,
N. S. Pol,
A. Sesana,
S. B. Spolaor,
S. R. Taylor,
G. Theureau,
M. Vallisneri,
S. J. Vigeland,
C. A. Witt,
X. Zhu,
J. Antoniadis,
Z. Arzoumanian,
M. Bailes
, et al. (102 additional authors not shown)
Abstract:
The International Pulsar Timing Array 2nd data release is the combination of datasets from worldwide collaborations. In this study, we search for continuous waves: gravitational wave signals produced by individual supermassive black hole binaries in the local universe. We consider binaries on circular orbits and neglect the evolution of orbital frequency over the observational span. We find no evi…
▽ More
The International Pulsar Timing Array 2nd data release is the combination of datasets from worldwide collaborations. In this study, we search for continuous waves: gravitational wave signals produced by individual supermassive black hole binaries in the local universe. We consider binaries on circular orbits and neglect the evolution of orbital frequency over the observational span. We find no evidence for such signals and set sky averaged 95% upper limits on their amplitude h 95 . The most sensitive frequency is 10nHz with h 95 = 9.1 10-15 . We achieved the best upper limit to date at low and high frequencies of the PTA band thanks to improved effective cadence of observations. In our analysis, we have taken into account the recently discovered common red noise process, which has an impact at low frequencies. We also find that the peculiar noise features present in some pulsars data must be taken into account to reduce the false alarm. We show that using custom noise models is essential in searching for continuous gravitational wave signals and setting the upper limit.
△ Less
Submitted 19 March, 2023;
originally announced March 2023.
-
A Parallelized Bayesian Approach To Accelerated Gravitational-Wave Background Characterization
Authors:
Stephen R. Taylor,
Joseph Simon,
Levi Schult,
Nihan Pol,
William G. Lamb
Abstract:
The characterization of nanohertz-frequency gravitational waves (GWs) with pulsar-timing arrays requires a continual expansion of datasets and monitored pulsars. Whereas detection of the stochastic GW background is predicated on measuring a distinctive pattern of inter-pulsar correlations, characterizing the background's spectrum is driven by information encoded in the power spectra of the individ…
▽ More
The characterization of nanohertz-frequency gravitational waves (GWs) with pulsar-timing arrays requires a continual expansion of datasets and monitored pulsars. Whereas detection of the stochastic GW background is predicated on measuring a distinctive pattern of inter-pulsar correlations, characterizing the background's spectrum is driven by information encoded in the power spectra of the individual pulsars' time series. We propose a new technique for rapid Bayesian characterization of the stochastic GW background that is fully parallelized over pulsar datasets. This Factorized Likelihood (FL) technique empowers a modular approach to parameter estimation of the GW background, multi-stage model selection of a spectrally-common stochastic process and quadrupolar inter-pulsar correlations, and statistical cross-validation of measured signals between independent pulsar sub-arrays. We demonstrate the equivalence of this technique's efficacy with the full pulsar-timing array likelihood, yet at a fraction of the required time. Our technique is fast, easily implemented, and trivially allows for new data and pulsars to be combined with legacy datasets without re-analysis of the latter.
△ Less
Submitted 30 March, 2022; v1 submitted 16 February, 2022;
originally announced February 2022.
-
Wide-band timing of GMRT discovered millisecond pulsars
Authors:
Shyam S. Sharma,
Jayanta Roy,
Bhaswati Bhattacharyya,
Lina Levin,
Ben Stappers,
Timothy T. Pennucci,
Levi Schult,
Shubham Singh,
Aswathy Kaninghat
Abstract:
Modeling of frequency-dependent effects, contributed by the turbulence in the free electron density of interstellar plasma, is required to enable the detection of the expected imprints from the stochastic gravitational-wave (GW) background in pulsar timing data. In this work, we present an investigation of temporal variations of interstellar medium for a set of millisecond pulsars (MSPs) with the…
▽ More
Modeling of frequency-dependent effects, contributed by the turbulence in the free electron density of interstellar plasma, is required to enable the detection of the expected imprints from the stochastic gravitational-wave (GW) background in pulsar timing data. In this work, we present an investigation of temporal variations of interstellar medium for a set of millisecond pulsars (MSPs) with the upgraded GMRT aided by large fractional bandwidth at lower observing frequencies. Contrary to the conventional narrow-band analysis using a frequency invariant template profile, we applied $PulsePortraiture$ based wide-band timing analysis while correcting for the evolution of the pulsar profile with frequency. Implementation of $PulsePortraiture$ based wide-band timing method for the GMRT discovered MSPs to probe the DM variations resulted in a DM precision of $10^{-4}\,pc~cm^{-3}$. In general, we achieve similar DM and timing precision from wide-band timing compared to the narrow-band timing with matching temporal variations of DMs. This wide-band timing study of newly discovered MSPs over a wide frequency range highlights the effectiveness of profile-modeling at low frequencies and probes the potential of using them in pulsar timing array.
△ Less
Submitted 24 August, 2022; v1 submitted 12 January, 2022;
originally announced January 2022.
-
The International Pulsar Timing Array second data release: Search for an isotropic Gravitational Wave Background
Authors:
J. Antoniadis,
Z. Arzoumanian,
S. Babak,
M. Bailes,
A. -S. Bak Nielsen,
P. T. Baker,
C. G. Bassa,
B. Becsy,
A. Berthereau,
M. Bonetti,
A. Brazier,
P. R. Brook,
M. Burgay,
S. Burke-Spolaor,
R. N. Caballero,
J. A. Casey-Clyde,
A. Chalumeau,
D. J. Champion,
M. Charisi,
S. Chatterjee,
S. Chen,
I. Cognard,
J. M. Cordes,
N. J. Cornish,
F. Crawford
, et al. (101 additional authors not shown)
Abstract:
We searched for an isotropic stochastic gravitational wave background in the second data release of the International Pulsar Timing Array, a global collaboration synthesizing decadal-length pulsar-timing campaigns in North America, Europe, and Australia. In our reference search for a power law strain spectrum of the form $h_c = A(f/1\,\mathrm{yr}^{-1})^α$, we found strong evidence for a spectrally…
▽ More
We searched for an isotropic stochastic gravitational wave background in the second data release of the International Pulsar Timing Array, a global collaboration synthesizing decadal-length pulsar-timing campaigns in North America, Europe, and Australia. In our reference search for a power law strain spectrum of the form $h_c = A(f/1\,\mathrm{yr}^{-1})^α$, we found strong evidence for a spectrally-similar low-frequency stochastic process of amplitude $A = 3.8^{+6.3}_{-2.5}\times10^{-15}$ and spectral index $α= -0.5 \pm 0.5$, where the uncertainties represent 95\% credible regions, using information from the auto- and cross-correlation terms between the pulsars in the array. For a spectral index of $α= -2/3$, as expected from a population of inspiralling supermassive black hole binaries, the recovered amplitude is $A = 2.8^{+1.2}_{-0.8}\times10^{-15}$. Nonetheless, no significant evidence of the Hellings-Downs correlations that would indicate a gravitational-wave origin was found. We also analyzed the constituent data from the individual pulsar timing arrays in a consistent way, and clearly demonstrate that the combined international data set is more sensitive. Furthermore, we demonstrate that this combined data set produces comparable constraints to recent single-array data sets which have more data than the constituent parts of the combination. Future international data releases will deliver increased sensitivity to gravitational wave radiation, and significantly increase the detection probability.
△ Less
Submitted 11 January, 2022;
originally announced January 2022.