SciELO - Scientific Electronic Library Online

 
vol.68 número4Cosmological-static metric correspondence and Kruskal type solutions from symmetry transformationsInitial value problem for a Caputo space-time fractional Schrödinger equation for the delta potential índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados

Revista

Articulo

Indicadores

Links relacionados

  • No hay artículos similaresSimilares en SciELO

Compartir


Revista mexicana de física

versión impresa ISSN 0035-001X

Rev. mex. fis. vol.68 no.4 México jul./ago. 2022  Epub 19-Mayo-2023

https://doi.org/10.31349/revmexfis.68.040702 

Research

A dynamical dark energy solution to hubble tension in the light of the multimessenger era

C. Escamilla-Rivera1  * 

1Instituto de Ciencias Nucleares, Universidad Nacional Autónoma de México, Circuito Exterior C.U., Apartado Postal 70-543, México 04510, Cuidad de México.


Abstract

We show that the gravitational waves measurements have raised the opportunity to measure H 0 with dark sirens to within 2σ, the accuracy required to resolve the Hubble-Lemaître tension. There are two principal reasons for our results: (1) upgrades to GW LIGO-Virgo transient catalogues GWTC-1 and GWTC-2 enhance their sensitive with only 10% of contamination fraction, and (2) new dark sirens should help to constrain our dynamical EoS. In conjunction, sensitivity upgrades and a new dark energy model will facilitate an accurate inference of the Hubble-Lemaître constant H 0 to better with an ±0.077 error in comparison to the LIGO dark siren with +14.0/-7.0, which would further solidify the role of dark sirens in late dark energy for precision cosmology in the future.

Keywords: Gravitational waves; GWTC-1; GWTC-2

1. Introduction

The Hubble-Lemaître constant H 0 is a fundamental cosmological parameter that governs the expansion rate of our Universe. However, since the first release of the Cosmic Microwave Background (CMB) observations by the Planck Collaboration in 20131, the determination of H0 based on the standard model of cosmology ((CDM, which started to be in tensioni with the model-independent determination via calibrated local Supernovae Ia (SNeIa) by the Hubble Space telescope in 20112. An initial tension of around 2.4σ has been the worst over the last years. Attached to this, we are dealing with a more troublesome issue: we are able to perform very precise measurements of H0 in several independent ways, but such methods do not agree among each other at a statistically significant level. This intense debate have led us to two ways to determine the Hubble-Lemaître constant: one established by the Planck Collaboration in 20183 assuming the (CDM model, and the other by the SH0ES Collaboration in 20194. Up to now, H0 inferredii from the fluctuation spectrum of the CMB (H0 = 67.36 ± 0.54 km s-1 Mpc-1) disagrees with the value obtained from the measurement of the luminosity distance and redshift to SNeIa (H0 = 74.03 ± 1.42 km s-1 Mpc-1) at 4.0 -5.8σ significance5. Recently, other probes have been made in order to clarify the issues, but still without success: Gravitationally-lensed time delays6,7 from quasars analysed from the H0LiCOW collaboration9 gives H0=73.3-1.8+1.7 km s-1 Mpc-1a 2.4% with a 3.1( tension with Planck 2018. The Dark Energy Survey (DES) collaboration9 has used SNeIa and Baryon Acoustic Oscillations (BAO) in an inverse distance ladder method to obtain H 0 = 67.8 ± 1.3 km s-1 Mpc-1.10 And finally, upgrades in the analysis of water masers in NGC425811 have give us a value of H 0 = 73.5 ± 1.4 km s-1 Mpc-1, a 4.2( tension with Planck 2018. A full compendium of Hubble-Lemaître constant estimations are detailed in Fig 1.

Figure 1 Hubble Constant estimations according to the publication year. Planck 2018-based, early universe mean value for H0 is 67 km/s/Mpc. The Cepheid-based, late universe mean value is 74 km/s/Mpc. A new alternative to Cepheids, knows as red giant stars that flare with a known intrinsic brightness, only complicated the tension. They indicated a H0 mean value of about 70 km/s/Mpc, a value in the midway between the other two, but in agreement with both the measurements within 2(. Also, we indicated the independent estimation from the multi-messenger detection of a binary neutron star inspiral H0 = 70.0+12.0-8.0 km s-1 Mpc-1. Dark siren detected in the first and second observing runs of LIGO and Virgo estimated H0 = 68.0+14.0-7.0 km s-1 Mpc-1. As we shall discuss in Sec. 5, our proposal allows to determine a value of H0 = 68.0±0.77 km s-1 Mpc-1 using GW transient catalog solely. 

To understand the implications of this discrepancy it is, by far, the most severe problem the standard cosmology is facing now. Several solutions have been proposed to explain, solve or even alleviate such a tension. For example, in12 has been studied independent measurement of the host galaxy’s redshift with the luminosity distance from GW to infer H 0. Also, the KBC void has been considered an interesting proposal to solve naturally the Hubble-Lemaître tension in Milgromian dynamics13. Although there are studies that consider this local feature not important on H 0 measurements14.

Even if some criticisms inherent to the distance ladder method, both for the procedure and for the existence of a local void15 have been raised. Nowadays the main scenario focuses on translating such observational tension in a tension between our description and understanding of both late and early time physics16. From one part, the effect of the local structure on H 0 (so-called cosmic variance) has been thoroughly studied, as well as possible reassessments of its error budget. Furthermore, physics beyond the standard model has been also investigated, in the hope that this discrepancy could reveal possible alternatives to the highly tuned cosmological constant and the yet-undetected dark matter. Confirm or dismiss this discrepancy is of extreme importance as it may point out to new (or missing) physics beyond the recombination era.

Currently, a panorama has been open to understand much deeper this tension issue: the multi-messenger Gravitational Wave (GW) astronomy. The first combined detection of GW and EM waves from the same source has brought light to analyze astrophysical17 and cosmological18,19 level at the same time. As for example, independent estimations can be obtained through the LIGO collaboration20, which value reported is H0=70.0-8.0+12.0 km s-1 Mpc-1 from the detection of a binary neutron star inspiral21. Dark siren detected in the first and second observing runs of LIGO and Virgo estimated H0=68.0-7.0+14.0 km s-1 Mpc-1.22 In comparison to other methods, this result from GW is settled in the middle. Unfortunately, as we can see from the error bars, this measurement cannot help to solve the tension, yet. Also, radio signals23 have been considered to compute a value of H0=68.9-4.6+4.7 km s-1 Mpc-1, which is consistent with both Planck 2018 and local estimators.

Until the next era of surveys is being ready in order to set light on this issue, our main objective should be to find new probes (or new alternative ways to employ current astrophysical data) which might add information to the argument and be competitive with measurements for what concerns the feasible precision. In that regard, new physics via an alternative theory of gravity should describe gravitational phenomena in a very wide range of systems, from the cosmic scales to compact objects. Hence, we can combine several observations to improve bounds on parameters. In the era of multi-messenger astronomy, modelling GW propagation is the new window in alternative theories of gravity and some of its inherent challenges, therefore in this paper we will study the possibility that late dark energy itself through a new EoS proposal can solve the Hubble-Lemaître tension in the late universe. We should mention that possibilities with an extra dark energy in the early universe has been considered24, moreover this involves first order transitions and many parameter assumptions to be fitted with the data. The aim of our study is slightly different: we will consider an exponential-like late dark energy parameterisation that in its first order approximation at present cosmic time recovers the standard dark energy models, e.g. CPL25,26, BA27, etc. and at higher order approximation can be compatible to be constrained using GW data to investigate how such corrections affect the H0 tension at late-time.

This paper is organised as follows: in 2 we present our proposal for a new late dark energy EoS and how it can have cosmological dependencies on the luminosity distance. In 3 we described the method to obtain the Hubble-Lemaître constant using a dark sirens and the requirements needed from GW LIGO-Virgo catalogues. In 4 we introduce the methodology to deal with the GW LIGO-Virgo transient catalogues GWTC-1 and GWTC-2 and study our late dark energy proposal. In 5 we discuss the statistical results for the cosmology derived and finally, in 6 we present our comments.

2. Testing Hubble-Lemaître with dark energy

A possible way to describe the current observational evidences is to introduce dark energy, a hypothetical fluid with negative pressure. Moreover, apart from this negativity property on the pressure, it is still unknown the nature of this dark component. The simplest and standard explanation comes through the introduction of a positive cosmological constant, Λ, which does not evolve with the time, but carries out with itself two major problems: one of which is the vacuum catastrophe and the other is the cosmic coincidence problem. These issues motivate us to consider scenarios beyond the standard Λ-cosmology landscape.

The simple scenario beyond the standard one, is the EoS w(z)-cosmology where w quantified the ratio of pressure p to its density ρ, i.e. w = p/ρ. Of course, a simple recovery of Λ-cosmology can be obtained easily when w = -1. Any deviation from this value can give a different cosmology, moreover, we shall focus on the alternative cosmologies in which the dark energy EoS is evolving with the Hubble flow.

Over the years, in the literature has been presented wide compendiums of w(z) parameterisations as any arbitrary function of the redshift z or the scale factor a of the Friedmann-Lemaître-Robertson-Walker universe. The advantage of such idea gives a complete freedom to select any specific model of interest and test it with the surveys at hand to study whether that model is able to correctly describe the evolution of the universe. In such a way, we can think this as an inverse evolving cosmology: we introduce a dark energy EoS to probe the expansion history of the universe. This is one of the motivation of using late dark energy as a model to solve the Hubble-Lemaître tension. While possibilities with an extra dark energy in the early universe has been considered24, remember that in this era we are restricted by using the only data we have at hand now: CMB data, which assume a (CDM-cosmology. Moreover, a design of late dark energy that can be tested using independent measurements as GW sirens, can shed some light on the tension problem without theoretical dependencies.

The aim of the present study consist in propose an exponential late dark energy parameterisation that in its zero and first order approximation around present time (z = 0) recovers (CDM and the CPL parameterisation, respectively. Furthermore, we allow its higher order approximation to study how such extended corrections can solve the Hubble-Lemaître tension. As it is standard, this kind of proposals emerge from a Taylor series inspiration. Let us start by proposing a late dark energy fluid with the following form:

w(z)=w0+w1ez1+z, (1)

where w(z = 0) = w0 + w1 is the current value of the dark energy EoS and w1 is another free parameter. We can notice some interesting features of this proposal:

  1. The pivoting redshift can be found in order to uncorrelated the parameters w0 and w1.28

  2. By this approach we will get important information on how the non-linear terms in the parametrised EoS affect the viability of the cosmological model in higher-redshifts.

  3. And finally, in comparison to other dynamical late dark energy whose deals with divergencies at low and high z on the Taylor series, our proposal by its construction includes a smoother evolution that can be controlled via the its first, second and third order corrections.

2.1. Cosmological dependencies

At this point, it is interesting to have an insight on how our proposal (1) might depend on the cosmological parameters, and in particular on H0. First of all, one should note that most of the dependence comes from the luminosity distance definition

DL(Θ)=(1+z)l0zcdz'H0E(z',Θ), (2)

and Θ = {w0, w 1} is the vector with the free cosmological parameters to be fitted and E = H/H0. To introduce our proposal (1) we describe the expansion history by the following expression:

H2(z)=H02Ωm(1+z)3+ΩΛf(z), (3)

where

f(z)=exp(30z1+w(z')1+z'dz'). (4)

This expression arise from the energy-momentum tensor conservation law for dark energy. Notice that (CDM is a special case w(z) = -1, so f(z) = 1. (m, is the matter density parameter today; (( = 1 - (m, the dark energy density parameter; and w(z) our proposal, which will be -1, in the case of dark energy as a cosmological constant. Our fiducial model is a flat cosmology, i.e. (k = 0.

H0 with Dark Sirens

To compute the Hubble-Lemaître constant we need to estimate the velocity of the Hubble flow at the position of a galaxy at certain redshift, e.g. at z = 0.01. In order to have a a precise estimation (less than 2%) it is require to correct peculiar velocities of the hosts of GW events, which correspond to a particular case for very low z sources. This is due to the high signal-to-noise ratio which eventually generate a bias in the joint posteriors for set of events. Therefore, for a standard value of a peculiar velocity around 300km/s for a GW host we consider z = 0.01, since the contribution from peculiar velocity is comparable to the term related with the Hubble flow. Using an EM calibration of the cosmic distance ladder, we can derive the distance to this galaxy using the Tully-Fisher relationship. Since in this relationship the method is still model-dependent, we estimate the velocity of the Hubble flow at the position of the galaxy considering corrections induced by local peculiar motions. Finally, we can constrain the value of the Hubble-Lemaître constant by using the distance and velocity distributions derived from the GW and EM data processed via Bayesian statistics. Treating with GW observation is something else entirely. Its method behind allow us to completely jump over the distance ladder and obtain an independent measurement of the distance to cosmological objects. The simultaneous observation of gravitational and EM waves is crucial to this idea since we compare two things: the distance to an object, and the apparent velocity with which it is moving away from the observer. Measure the redshift z of light is easy to do when we have an EM spectrum of a specific object. Moreover, with GW alone, we cannot do it - there is not enough structure in the spectrum to measure a z. Moreover, GW cannot give us standard candles per se, since every one will have a different intrinsic gravitational luminosity, but we can work what that luminosity is if we take into account the way in which the source evolves. The idea of using GW sources as standard sirens was study by B. Schutz in 198629. Since then, the idea been developed substantially since30, given us a step forward to use Bayesian analysis. For example, in Ref.31 it is described an analysis of mock data using binary neutron star mergers to recover an unbiased estimate of H0. In Fig. 2 we illustrate the statistical relationships between the data and the parameters to compute H0 using sirens.

Figure 2 Statistical relationships between the data and the parameters for a GW to compute H0. Blue circles indicate parameters that require a prior. The red circles describe the measured data, which is conditioned on the analysis. Here we assume that we have measurements of the GW data, xGW (set of detectors), a recession velocity (that is, redshift) vr, and the mean peculiar velocity vp in the galaxy’s neighbourhood. Arrows towards a node indicate that the conditional probability density for the node depends on the source parameters e.g. the conditional distribution of the observed GW data depends on the distance and inclination of the source (and additional parameters, here marginalized). 

Moreover, for specific phenomena as binary black hole (BBH) mergers, we do not expected to have EM counterparts but they can determine the luminosity distance to their hosts independently of the cosmic distance ladder method. These are usually referred to as dark sirens. The GW standard dark siren measurement is broadly consistent with other measurements32. If we combine the information from multiple detections, we can improve the accuracy reaching about 1% with O(100) detections in the next decade. Moreover, the BBH sources with no expected EM counterparts can be used with the method outlined in Ref.29: with a database of potential host galaxies identified in a galaxy catalog for each detection, we can build up data by a process of statistical cross-correlation. These results established the path to cosmology using GW observations with and without transient EM counterparts. The first results of the implementation of this method was presented in Ref.33 using a set of simulations with 5% estimate on H0 in a idealised situation where nearby events and complete galaxy catalogs have similar results for a third-generation detector34.

According to this, GW from compact binary coalescences, being standard sirens, can be used to calibrate distances to SNeIa, when both occur in the same galaxy or galaxy cluster35. Furthermore, it can be possible to identify their potential host galaxies either with the help of a galaxy catalog or looking coalescence effects from stellar-mass BBH36. As it is standard in the statistical analysis, we require taking into account the sources of systematic errors in this method, which arises from the incompleteness of the available galaxy catalogs. Nonetheless, LIGO and Virgo at their design sensitivity could constrain H0 with dark sirens to 5% accuracy with (250 detections37. It is in this direction that an independent estimation from the multi-messenger detection of a binary neutron star inspiral21 has obtained H0=70.0-8.0+12.0 km s-1 Mpc-1. Additionally, in our study we include GW searches for coalescing compact binaries from GW transient catalogs GWTC-122 and GWTC-238, which allow us to include larger solar masses and coalescence effects.

4. Methodology

LIGO and Virgo have published two compact binary mergers found during the third observing run: GW19041239 and 40, which are exceptional due to their large mass-asymmetry, with mass ratios around (3 and (9, respectively. This kind of systems have led to the detection of subdominant spherical harmonic modes beyond the quadrupole mode. Therefore, higher modes will be important in the parameter estimation of asymmetric binaries, e.g. of the luminosity distance DL and orbital inclination (. All these measurements can be attached with the new runs that includes upgrades to GW LIGO-Virgo transient catalogs GWTC-1 and GWTC-2, with sensitive of only 10% of contamination fraction.

In this work we are interested in compute local values of H0 with very distant sources. Consider that, if sources are near (e.g. at DL 100Mpc), its measurements could be imprecise due to the bias from peculiar velocities vp of host galaxies. In this regime, any proposal mimicking (CDM can work. Therefore, we will consider our dark energy proposal distributed uniformly in a volume up to a redshift of z = 2.3, including along GW measurements, supernovae and observational Hubble data.

4.1. Distance-ladder measurements

Pantheon SNeIa compilation41: contains 1048 SNeIa distribute in a redshift range of 0.01 < z < 2.26. The constraining power of this supernovae is due that can be used as standarizable candles. Its implementation is via

F(z,Θ)theo=5log10DL(z,Θ)+μ0, (5)

where DL is the luminosity distance given by (2) and μ0 is the nuisance parameters of the distance estimator. We can write (ℱ(() = ℱtheo - ℱobs, using for this purpose the distance modulus ℱobs associated with the observed magnitude and marginalise χSN2 with respect to (0

χSN2(Θ)=(ΔF(Θ))TCSN-1ΔF(Θ)+lnS2π-k2(Θ)S, (6)

where CSN is the total covariance matrix and S is the sum of all entries of CSN-1, weighed by a covariance matrix as k(Θ)=ΔF(Θ)T·CSN-1.[/p]

Observational Hubble data (CC)42: contain l = 51 measurements in the redshift range 0l.07 < z < 2.0. In this sample, 31 data points correspond to passive galaxies and 20 data points are estimated from BAO data under a (CDM prior. Our theoretical setting to construct χH2 is given by

χH2=i=1lHzi,x-Hobs(zi)2σH2(zi), (7)

where Hobs(zi) is the observed value at zi,H(zi) are the observational errors, and H(zi, x) is the value of a theoretical H for the same zi with the parameter vector x.

4.2. Gravitational waves measurements

Gravitational-Wave Transient Catalog GWTC-143. Includes NGWTC-1 = 11 confident detection events in redshift range of z ∈[0.01, 0.49] from GW searches for coalescing compact binaries, all of them with masses M > 1M. These events observed by LIGOiii and Virgoiv come from the first and second observing runs. Eleven events are reported along with their observational siren distances DS in Mpc and redshifts z and other several observational variables, all of them with their corresponding 90% credible variances. Its statistics can be described by

χGWTC-12=j=1NGWTC-1(DS(zi,h,ΩM,w0,w1,n,Ξ0)-DSobs, GWTC-1(zi))2σjobs, GWTC-12, (8)

where the quantities n and ≡0 are both positive and defined as44

Ξ0=limzM*(0)M*(z), (9)

nαM02(Ξ0-1), (10)

with M* = m2pl (G as the effective time-dependent Planck mass, where G is the gravitational strength. Also, this catalog includes variances in z, therefore the squared variances are given by

σjobs, GWTC-12=σjDSobs, GWTC-12+σjzobs, GWTC-12, (11)

where the siren luminosity distance is proposed as

DS(z)=DL(z)exp(-0zδ(z')1+z'dz'), (12)

with

δ1(z)=n(1-Ξ0)1-Ξ0+Ξ0(1+z)n. (13)

The squared variances of the siren distance per redshift are already added.

Gravitational-Wave Transient Catalog GWTC-238. This catalogue includes NGWTC -2 = 39 GW events, with less than 10% of contamination fraction in a range of z ∈[0.03, 0.8] from the first half of the third observable run. As in the latter catalogue, the statistics can be described by

χGWTC-22=j=1NGWTC-2(DS(zi,h,Ωm,w0,w1,n,Ξ0)-DSobs, GWTC-2(zi))2σjobs, GWTC-22, (14)

with variances in z as

σjobs, GWTC-22=σjDSobs, GWTC-22+σjzobs, GWTC-22. (15)

The squared variances of both siren distances and redshift are added up. Due that these events have less than 10% of contamination fraction, the best fit values and their confidence regions might be biased, even though some bias can be controlled when using the distance ladder measurements described above.

5. Results

We constrain our parameterisation (1) with distance ladder measurements, such as Pantheon and CC described above. Also, we consider another dark energy models (see Table I) to compare with our proposal (1). Afterwards, we employed the best fit values obtained (see Table II) to generated constraints on dark sirens in the transient catalog GWTC-1 and GWTC-2. In Fig. 3 we present the comparisons between the constraining results of the data samples after the inclusion of the GW data to distanc/e ladder measurements, where in particular, we show the (1-dimensional marginalised posterior distributions for the parameters of our model (w0, w1) and the ones derived from the GW ((0, n), as well as the 2-dimensional contour plots (C.L.) between the combinations Pantheon+GW, GW and Pantheon+CC+GW. In the following we describe the effects of GW on the ([eq:eos-general]) corresponding to different observational data samples.

Table I Description of the dynamical parameterisations used in this work. In the first column we denote the model and the reference and in the second column their expression for w(z). This table complements the analysis in Table II

Model Parameterisation
Our Model (1) wz=w0+w1ez1+z
Lambda Cold Dark
Matter (ΛCDM) w = -1
Chevallier-Polarski-
Linder (CPL)25,26 wz=w0+z1+zw1.
Barboza-Alcaniz (BA)27 wz=w0+z1+z1+z2w1.
Linar.redshift (LR)46 wz=w0+w1z.
Low Correlation (LC)47 wz=-z+zcw0+z1+zcw01+zzc
Wetterich (WP)48 wz=w01+w1ln1+z2

Table II Best fits cosmological parameters obtained using GW, Pantheon (SN) and CC. We include the ΛCDM model, our proposed model (1) and five dark energy bidimensional parameterisations: CPL [25,26], BA [27], LR [46], LC [47] and WP [48]. A full review of these models with late-universe data is presented in Ref. [45], and references therein. All these values are at 2σ of uncertainties C.L. 

Model Data h m w0 w1 Ξ0 n
(1) GW 0.680 ± 0.077 0.28+0.21-0.14 -1.4+2.2-1.6 -0.1+3.8-4.5 1.06+0.49-0.56 2.17+0.87-2.1
SN+GW 0.7344 ± 0.004 0.299+0.12-0.04 -1.21+0.17-0.15 0.6+2.3-1.0 1.18+0.24-0.28 1.8+1.1-1.8
SN+CC+GW 0.7348 ± 0.004 0.216+0.024-0.018 -1.043+0.074-0.083 0.89 ± 0.74 1.17+0.24-0.29 1.77+0.69-1.8
ΛCDM GW 0.690 ± 0.073 0.26 ± 0.14 - - 1.00 ± 0.40 1.93+0.77-1.9
SN+GW 0.728 ± 0.002 0.286 ± 0.013 - - 1.19+0.24-0.30 1.8+1.2-1.8
SN+CC+GW 0.733 ± 0.002 0.245 ± 0.007 - - 1.17+0.25-0.28 1.7+1.2-1.7
CPL GW 0.698 ± 0.078 0.286+0.21-0.079 -1.4+2.2-1.7 -0.1+2.9-46 1.04+0.48-0.57 2.2+1.5-2.1
SN+GW 0.734 ± 0.004 0.301+0.13-0.046 −1.18 ± 0.17 0.08+1.7-0.62 1.18+0.24-0.30 1.8+1.1-1.8
SN+CC+GW 0.734 ± 0.003 0.169+0.077-0.030 -0.951+0.054-0.074 0.61+0.46-0.18 1.16 ± 0.34 1.75+0.98-1.8
BA GW 0.698 ± 0.078 0.281+0.21-0.093 -1.4+2.3-1.6 -0.2+2.1-4.4 1.06+0.51-0.58 2.2+1.1-2.2
SN+GW 0.7329 ± 0.004 0.332+0.11-0.033 -1.20+0.15-0.19 -0.34+1.3-0.42 1.18+0.25-0.29 1.8+1.6-1.8
SN+CC+GW 0.7331 ± 0.003 0.157+0.088-0.037 -0.913+0.056-0.092 0.324+0.21-0.056 1.17+0.24-0.29 1.8+1.0-1.8
LR GW 0.698+0.085−0.075 0.282+0.21-0.092 -1.4+2.3-1.6 0.2 ± 2.9 1.07 ± 0.49 2.2+1.2-2.1
SN+GW 0.7328 ± 0.004 0.324+0.12-0.034 -1.18+0.15-0.20 0.33+0.48-1.4 1.17+0.25-0.30 1.75+0.69-1.8
SN+CC+GW 0.7314 ± 0.003 0.211+0.043-0.020 -0.937+0.044-0.056 -0.098+0.095-0.24 1.16+0.25-0.28 1.77+0.94-1.8
LC GW 0.691+0.098−0.077 0.28+0.21-0.13 -0.9+1.7-3.3 -0.8+1.8-3.4 1.12 ± 0.48 2.18+0.88-2.1
SN+GW 0.7333 ± 0.004 0.28+0.21 -0.13 −1.18 ± 0.17 -1.21+0.77-0.32 1.17+0.25-0.29 1.8+1.4-1.8
SN+CC+GW 0.7337 ± 0.003 0.172+0.075-0.029 -0.953+0.051-0.069 -0.75+0.18-0.15 1.17+0.26-0.29 1.8+1.4-1.8
WP GW 0.695 ± 0.079 0.27+0.22-0.12 −1.1 ± 2.1 2.3+2.6-1.7 1.14+0.68-0.44 2.2+1.6-2.2
SN+GW 0.7351 ± 0.004 0.22 ± 0.12 -1.14+0.19-0.15 0.67+0.76-0.58 1.17 ± 0.33 1.8+1.0-1.8
SN+CC+GW 0.7353 ± 0.003 0.134+0.08-0.056 −0.958 ± 0.054 0.58+0.36-0.28 1.17+0.25-0.29 1.8+1.1-1.8

Figure 3 For our model (1), we first constrain the cosmological parameters (h, ( m , w 0, w 1, (0 and n) using the datasets Pantheon+GW (including both GWTC-1 and GWTC-2) (black solid lines), GW alone (green C.L.) and Pantheon+CC+GW (blue C.L.), respectively. The best-fit values of the parameters for each dataset are reported in Table II

In Table I, we present five of the most popular dark energy parameterisations in the literature (a study of these models with current data was presented in Ref.45). Also, we include in the analysis of our model (1) and the (CDM model.

In the second column of this table, we described the data combinations employed to perform the fit of the cosmological parameters. One can clearly notice that the treatment of GW solely for (CDM and the five dark energy models shows a Hubble-Lemaîtreconstant value in the periphery with the H 0 inferred from the fluctuation spectrum of the CMB, while in comparison to the other proposals, our model gives a value of H 0 = 68.0(0.077 km s-1 Mpc-1, closer to the dark siren detected in the first and second observing runs of LIGO and Virgo estimated H 0 = 68.0 −7.0 +14.0 km s-1 Mpc-1. The error bars seems to decrease when we consider, along with the GWTC-1 catalogue, the GWTC-2, this consequence is directly connect to the fact that we are dealing with a sensitivity of 10% of contamination fraction in the catalogues and their resolution allow us to measure H0 with dark sirens to within 2(. We also note that the matter density parameter (m0, for GW alone is better constrained compared to the Planck’s estimation and the inclusion of GW to SN and CC improves the parameter space.

6. Conclusions

We have demonstrated a possibility of measuring the Hubble-Lemaître constant to (2(-level precision using dark sirens with the upgrades of the LIGO and Virgo detectors through the transient catalogs GWTC-1 and GWTC-2. Although the standard model (CDM has still some advantage in comparison other dark energy parameterisations (see Table II), GWs have opened the door to update the analysis previously done with supernovae and Hubble data. In fact, from the Hubble-Lemaître constant fit derived, our proposal (1) offers the opportunity to measure H0 with dark sirens GW to within 2(, the accuracy required to resolve the Hubble-Lemaîtretension, to obtain a value of H0 = 68.0(0.077 km s-1 Mpc-1, closer to the dark siren detected in the first and second observing runs of LIGO and Virgo estimated.

Future GW data might give us more clues to discover the actual model of dark energy, whether it continues to be (CDM or another one. The analyses done in this work will be improved with future numerous and accurate data and establish the preference of GR over the siren distance - luminosity distance non equality coming from modified gravity models49. This analysis will be reported elsewhere.

In our opinion, only a complete and coherent theory of dark energy selected by high level of precision and control of systematic in observations can infer the correct cosmological paradigm. We expect that the ideas presented in this work could give further insight into the Hubble-Lemaîtretension issue; but that we shall know from the future work.

Acknowledgments

CE-R acknowledges the Royal Astronomical Society as FRAS 10147, supported by DGAPA-PAPIIT-UNAM Projects IA100220 and TA100122. The simulations were performed in Centro de cómputo Tochtli-ICN-UNAM as part of CosmoNag.

References

1. P. A. R. Ade et al., Planck 2013 results. XVI. Cosmological parameters. Astron. Astrophys., 571 (2014) A16, [ Links ]

2. A. G. Riess et al., A space telescopeand wide field camera 3. The Astrophysical Journal, 730 (2011) 119, [ Links ]

3. N. Aghanim et al., Planck 2018 results. VI. Cosmological parameters. Astron. Astrophys., 641 (2020) A6, [ Links ]

4. A. G. Riess, S. Casertano, W. Yuan, L. M. Macri, and D. Scolnic, Large Magellanic Cloud Cepheid Standards Provide a 1% Foundation for the Determination of the Hubble Constant and Stronger Evidence for Physics beyond ( Astrophys. J. 876 (2019) 85, [ Links ]

5. L. Verde, T. Treu, and A. G. Riess, Tensions between the Early and the Late Universe. (In Nature Astronomy 2019). [ Links ]

6. G. C. F. Chen et al., A SHARP view of H0LiCOW: H0 from three time-delay gravitational lens systems with adaptive optics imaging. (2019). [ Links ]

7. Kenneth C. Wong et al., H0LiCOW XIII. A 2.4% measurement of H0 from lensed quasars: 5:3( tension between early and late-Universe probes. (2019). [ Links ]

8. S. H. Suyu et al., H0LiCOW I. H0 Lenses in COSMOGRAIL’s Wellspring: program overview. Mon. Not. Roy. Astron. Soc., 468 (2017) 2590, [ Links ]

9. T. M. C. Abbott et al., The Dark Energy Survey Data Release 1. Astrophys. J. Suppl., 239 (2018) 18, [ Links ]

10. E. Macaulay et al., First Cosmological Results using Type Ia Supernovae from the Dark Energy Survey: Measurement of the Hubble Constant. Mon. Not. Roy. Astron. Soc., 486 (2019) 2184, [ Links ]

11. M. J. Reid, D. W. Pesce, and A. G. Riess, An Improved Distance to NGC 4258 and its Implications for the Hubble Constant. (2019). [ Links ]

12. S. Borhanian, A. Dhani, A. Gupta, K. G. Arun, and B. S. Sathyaprakash, Dark Sirens to Resolve the 10 Hubble-Lemaître Tension. Astrophys. J. Lett., 905 (2020) L28. [ Links ]

13. M. Haslbauer, I. Banik, and P. Kroupa, The KBC void and Hubble tension contradict 𝛥CDM on a Gpc scale Milgromian dynamics as a possible solution. Mon. Not. Roy. Astron. Soc., 499 (2020) 2845, [ Links ]

14. W. D’Arcy Kenworthy, D. Scolnic, and A. Riess, The Local Perspective on the Hubble Tension: Local Structure Does Not Impact Measurement of the Hubble Constant. Astrophys. J., 875 (2019) 145, [ Links ]

15. T. Shanks, L. Hogarth, and N. Metcalfe, Gaia Cepheid parallaxes and ’Local Hole’ relieve H0 tension. Mon. Not. Roy. Astron. Soc., 484 (2019) L64, [ Links ]

16. E. Di Valentino et al., Cosmology Intertwined II: The Hubble Constant Tension. 8 (2020) [ Links ]

17. P. Mészáros, D. B. Fox, C. Hanna, and K. Murase, Multi-Messenger Astrophysics. Nature Rev. Phys. 1 (2019) 585, [ Links ]

18. J. M. Ezquiaga and M. Zumalacárregui, Dark Energy in light of Multi-Messenger Gravitational-Wave astronomy. Front. Astron. Space Sci., 5 (2018) 44 [ Links ]

19. M. Corman, C. Escamilla-Rivera, and M. A. Hendry, Constraining extra dimensions on cosmological scales with LISA future gravitational wave siren data. 4 (2020). [ Links ]

20. B. P. Abbott et al., A gravitational-wave standard siren measurement of the Hubble constant. Nature, 551 (2017) 85, [ Links ]

21. B. P. Abbott et al., GW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral. Phys. Rev. Lett. 119 (2017) 161101, [ Links ]

22. B. P. Abbott et al., GWTC-1: A Gravitational-Wave Transient Catalog of Compact Binary Mergers Observed by LIGO and Virgo during the First and Second Observing Runs. Phys. Rev. X, 9 (2019) 031040, [ Links ]

23. Kenta Hotokezaka et al., A Hubble constant measurement from superluminal motion of the jet in GW170817. Nature Astron. (2019). [ Links ]

24. F. Niedermann and M. S. Sloth, New Early Dark Energy. 10 (2019) [ Links ]

25. M. Chevallier and D. Polarski, Accelerating universes with scaling dark matter. Int. J. Mod. Phys. D 10 (2001) 213, [ Links ]

26. Eric V. Linder. The Dynamics of Quintessence, The Quintessence of Dynamics. Gen. Rel. Grav., 40 (2008) 329, [ Links ]

27. E. M. Barboza J.r. and J.S. Alcaniz, A parametric model for dark energy. Physics Letters B, 666 (2008) 415, [ Links ]

28. W. Yang, S. Pan, E. Di Valentino, and E. N. Saridakis. Observational constraints on dynamical dark energy with pivoting redshift. Universe 5 (2019) 219, [ Links ]

29. B. F. Schutz, Determining the Hubble Constant from Gravitational Wave Observations. Nature, 323 (1986) 310, [ Links ]

30. D. E. Holz and S. A. Hughes, Using gravitational-wave standard sirens. Astrophys. J. 629 (2005) 15, [ Links ]

31. R. Gray et al., Cosmological inference using gravitational wave standard sirens: A mock data analysis. Phys. Rev. D, 101 (2020) 122001 [ Links ]

32. B. P. Abbott et al., A Gravitational-wave Measurement of the Hubble Constant Following the Second Observing Run of Advanced LIGO and Virgo. Astrophys. J., 909 (2021) 218, [ Links ]

33. W. Del Pozzo, Inference of the cosmological parameters from gravitational waves: application to second generation interferometers. Phys. Rev. D, 86 (2012) 043011, [ Links ]

34. R. Abbott et al., Search of the Early O3 LIGO Data for Continuous Gravitational Waves from the Cassiopeia A and Vela Jr. Supernova Remnants. 11 (2021) [ Links ]

35. A. Gupta, D. Fox, B. S. Sathyaprakash, and B. F. Schutz, Calibrating the cosmic distance ladder using gravitational-wave observations. 10 (2019). [ Links ]

36. A. Nishizawa, Measurement of Hubble constant with stellar-mass binary black holes. Phys. Rev. D, 96 (2017) 101303, [ Links ]

37. R. Nair, S. Bose, and T. Deep Saini, Measuring the Hubble constant: Gravitational wave observations meet galaxy clustering. Phys. Rev. D, 98 (2018) 023502, [ Links ]

38. R. Abbott et al., GWTC-2: Compact Binary Coalescences Observed by LIGO and Virgo During the First Half of the Third Observing Run. 10 (2020) [ Links ]

39. R. Abbott et al., GW190412: Observation of a Binary-Black-Hole Coalescence with Asymmetric Masses. Phys. Rev. D, 102 (2020) 043015, [ Links ]

40. R. Abbott et al., GW190814: Gravitational Waves from the Co-alescence of a 23 Solar Mass Black Hole with a 2.6 Solar Mass Compact Object. Astrophys. J. Lett., 896 (2020) L44, [ Links ]

41. D. M. Scolnic et al., The Complete Light-curve Sample of Spectroscopically Confirmed SNe Ia from Pan-STARRS1 and Cosmological Constraints from the Combined Pantheon Sample. Astrophys. J., 859 (2018) 101, [ Links ]

42. M. Moresco et al., A 6% measurement of the Hubble parameter at z (0.45 direct evidence of the epoch of cosmic re-acceleration. JCAP, 05 (2016) 014, [ Links ]

43. P. Abbott et al., Gwtc-1: a gravitational-wave transient catalog of compact binary mergers observed by ligo and virgo during the first and second observing runs. Phys. Rev. X, 9 (2019) 031040, [ Links ]

44. A. Mitra, J. Mifsud, D. F. Mota, and . Parkinson, Cosmology with the Einstein Telescope: No Slip Gravity Model and Redshift Specifications. 9 (2020). [ Links ]

45. C. Escamilla-Rivera, Status on bidimensional dark energy parameterizations using SNe Ia JLA and BAO datasets. Galaxies, 4 (2016) 8, [ Links ]

46. D. Huterer and M. S. Turner, Probing dark energy: Methods and strategies. Phys. Rev. D, 64 (2001) 123527, [ Links ]

47. Y. Wang, Figure of merit for dark energy constraints from current observational data. Phys. Rev. D, 77 (2008) 123525, [ Links ]

48. Christof Wetterich. Phenomenological parameterization of quintessence. Physics Letters B, 594 (2004) 17. [ Links ]

49. E. Belgacem et al., Testing modified gravity at cosmological distances with lisa standard sirens. Journal of Cosmology and Astroparticle Physics, 2019 (2019) 024 [ Links ]

i Let us consider two different measurements of the same parameter A. Let A1 = a1 ( b1, be the first measurement and A2 = a2 ( b2, the second measurement. With a1, a2 the mean values and b1, b2 the standard deviation. If we assume that these measurements are independent and that the uncertainties have approximately a normal distribution, then the tension between A1 and A2 is given by: T = (a1 - a2(/ (b21 + b22)1/2(. This T quantity indicates the number of standard deviations by which A1 and A2 differ.

ii We say inferred, because the determination of a H0 value directly relies on the assumption of a cosmological model (in the method employed by Planck it is a standard (CDM model), and is strongly connected with early cosmic times physics (prior to recombination).

Received: November 16, 2021; Accepted: December 10, 2021

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License