Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

On the importance of source population models for gravitational-wave cosmology

S. Mastrogiovanni, K. Leyde, C. Karathanasis, E. Chassande-Mottin, D. A. Steer, J. Gair, A. Ghosh, R. Gray, S. Mukherjee, and S. Rinaldi
Phys. Rev. D 104, 062009 – Published 20 September 2021
PDFHTMLExport Citation
Abstract
Authors
Article Text
  • INTRODUCTION
  • HIERARCHICAL BAYESIAN ANALYSIS
  • APPLICATION OF THE INFERENCE SCHEME TO A…
  • CORRELATIONS BETWEEN COSMOLOGICAL AND…
  • IMPACT OF POPULATION MISCALIBRATION ON…
  • IMPACT OF POPULATION ASSUMPTIONS ON THE…
  • IMPACT OF THE POPULATION PARAMETERS ON H…
  • CONCLUSIONS
  • ACKNOWLEDGMENTS
  • APPENDICES
    • References

    Abstract

    Knowledge of the shape of the mass spectrum of compact objects can be used to help break the degeneracy between the mass and redshift of the gravitational wave (GW) sources and thus can be used to infer cosmological parameters in the absence of redshift measurements obtained from electromagnetic observations. In this paper, we study extensively different aspects of this approach, including its computational limits and achievable accuracy. Focusing on ground-based detectors with current and future sensitivities, we first perform the analysis of an extensive set of simulated data using a hierarchical Bayesian scheme that jointly fits the source population and cosmological parameters. We consider a population model (power-law plus Gaussian) which exhibits characteristic scales (extremes of the mass spectrum, presence of an accumulation point modeled by a Gaussian peak) that allow an indirect estimate of the source redshift. Our analysis of this catalog highlights and quantifies the tight interplay between source population and cosmological parameters, as well as the influence of initial assumptions (whether formulated on the source or cosmological parameters). We then validate our results by an “end-to-end” analysis using simulated GW h(t) data and posterior samples generated from Bayesian samplers used for GW parameter estimation, thus mirroring the analysis chain used for observational data for the first time in literature. Our results then lead us to re-examine the estimation of H0 obtained with GWTC-1 in Abbott et al. [LIGO Scientific, Virgo Collaborations, Astrophys. J. 909, 218 (2021)], and we show explicitly how population assumptions impact the final H0 result. Together, our results underline the importance of inferring source population and cosmological parameters simultaneously (and not separately as is often assumed). The only exception, as we discuss, is if an electromagnetic counterpart was to be observed for all the BBH events; then, the population assumptions have less impact on the estimation of cosmological parameters.

    • Figure
    • Figure
    • Figure
    • Figure
    • Figure
    • Figure
    • Figure
    4 More
    • Received 26 March 2021
    • Revised 25 June 2021
    • Accepted 2 September 2021

    DOI:https://doi.org/10.1103/PhysRevD.104.062009

    © 2021 American Physical Society

    Physics Subject Headings (PhySH)

    1. Research Areas
    Classical black holesCosmological parametersCosmologyGravitational wave sourcesGravitational waves
    Gravitation, Cosmology & Astrophysics

    Authors & Affiliations

    S. Mastrogiovanni1, K. Leyde1, C. Karathanasis2, E. Chassande-Mottin1, D. A. Steer1, J. Gair3, A. Ghosh4, R. Gray5, S. Mukherjee6,7,8, and S. Rinaldi9,10

    • 1Université de Paris, CNRS, AstroParticule et Cosmologie (APC), F-75013 Paris, France
    • 2Institut de Física d’Altes Energies (IFAE), Barcelona Institute of Science and Technology, Barcelona, Spain
    • 3Max Planck Institute for Gravitational Physics (Albert Einstein Institute), Am Mühlenberg 1, Potsdam 14476, Germany
    • 4Ghent University, Proeftuinstraat 86, 9000 Gent, Belgium
    • 5SUPA, University of Glasgow, Glasgow G12 8QQ, United Kingdom
    • 6Gravitation Astroparticle Physics Amsterdam (GRAPPA), Anton Pannekoek Institute for Astronomy and Institute for Physics, University of Amsterdam, Science Park 904, 1090 GL Amsterdam, The Netherlands
    • 7Institute Lorentz, Leiden University, P.O. Box 9506, Leiden 2300 RA, The Netherlands
    • 8Delta Institute for Theoretical Physics, Science Park 904, 1090 GL Amsterdam, The Netherlands
    • 9Dipartimento di Fisica “E. Fermi”, Universit di Pisa, I-56127 Pisa, Italy
    • 10INFN, Sezione di Pisa, I-56127 Pisa, Italy

    Article Text (Subscription Required)

    Click to Expand

    References (Subscription Required)

    Click to Expand
    Issue

    Vol. 104, Iss. 6 — 15 September 2021

    Reuse & Permissions
    Access Options
    CHORUS

    Article Available via CHORUS

    Download Accepted Manuscript
    Author publication services for translation and copyediting assistance advertisement

    Authorization Required

    Log In
    Other Options
    ×

    Images

    • Figure 1
      Figure 1

      Simulated population of 1024 observed events, showing the mass distributions (in the detector and source frames) and redshift distribution.

      Reuse & Permissions
    • Figure 2
      Figure 2

      Posterior probability density distributions on the different population parameters (PLG mass model) as more and more GW detections are analyzed (horizontal axis). The horizontal black dashed line indicates the true parameters of the population. The blue posteriors are obtained by fixing Ωm,0=0.308, while the orange posteriors are marginalized over the estimation of Ωm,0.

      Reuse & Permissions
    • Figure 3
      Figure 3

      (Left) Accuracy at the 90% CL on the population parameters (see levels) when we combine more and more GW detections. For this case, we fix Ωm,0=0.308 to the injected value. (Right) Same but varying Ωm,0 in the range ∈[0.1,0.5].

      Reuse & Permissions
    • Figure 4
      Figure 4

      Posterior distribution on H0 and Ωm,0 for 1024 BBH events detected with LIGO and Virgo at current sensitivities. The blue lines show the true parameters.

      Reuse & Permissions
    • Figure 5
      Figure 5

      Cumulative distributions of the population source-frame masses built by stacking the events posterior samples of the detected binaries (top, primary mass; bottom, secondary mass) inferred from the simulated population and assuming different values for H0. The position of the Gaussian peak and the maximum BH mass for BH are indicated by red and green areas, respectively, see Sec. 4b for an interpretation of this plot. The original prior has not been removed.

      Reuse & Permissions
    • Figure 6
      Figure 6

      Posterior distribution on the H0, mmax, and μg for 64 BBH events detected with LIGO and Virgo at current sensitivities. The blue lines show the true parameters. The contours indicate the 1σ and 2σ confidence level intervals.

      Reuse & Permissions
    • Figure 7
      Figure 7

      Posterior probability density distributions for the population parameters and H0 using posterior from 200 events generated with full parameter estimation. The blue lines indicate the population injected values. Levels indicate the 68% and 90% confidence intervals.

      Reuse & Permissions
    • Figure 8
      Figure 8

      Posterior distribution for H0 obtained by fixing mmax and μg in a range around their true values mmax=85  M⊙ and μg=40  M⊙. The black dashed line indicates the true value of H0.

      Reuse & Permissions
    • Figure 9
      Figure 9

      (Top) Posterior predictive check labeled as “expected” for the two source-frame masses and redshift fitting an underlying PLG population with the correct mass model in comparison with 64 observed events, labeled as “detected”. The detected distributions match the posterior predictive checks. (Bottom) Same but fitting a subcomplete PL model to the underlying PLG model. The model struggles to fit the lower end part of the mass distribution and the excess of BHs around 40–50  M⊙.

      Reuse & Permissions
    • Figure 10
      Figure 10

      Posterior distribution on H0 using the six GWTC-1 events with SNR>12 and the GLADE and DES galaxy catalogs. The plot compares the results obtained in [25] with the new results of this paper, see discussion in Sec. 6.

      Reuse & Permissions
    • Figure 11
      Figure 11

      Hubble constant posterior generated from 64 of our synthetic population of BBHs (Sec. 3a) fixing different population models and providing the redshift of the GW source (assumed from an EM counterpart). An incorrect choice of one of the population parameters (see legend) does not affect significantly the H0 estimation. The vertical dashed line indicates the injected value. Note that for this plot the posterior samples are generated taking into account the correlations between masses and luminosity distance [point (ii) above is dropped].

      Reuse & Permissions
    ×

    Sign up to receive regular email alerts from Physical Review D

    Log In

    Cancel
    ×

    Search

    Article Lookup

    Paste a citation or DOI
    Enter a citation
    ×