Local determination of the Hubble constant and the deceleration parameter
David Camarena Valerio Marra






PPGCosmo, Universidade Federal do Espírito Santo, 29075-910, Vitória, ES, Brazil


Núcleo Cosmo-ufes & Departamento de Física, Universidade Federal do Espírito Santo, 29075-910, Vitória, ES, Brazil



The determination of the Hubble constant H_0 from 
the Cosmic Microwave Background by the Planck Collaboration <cit.>  is in tension at 4.2σ with respect to the local determination of H_0 by the SH0ES collaboration <cit.>. 
Here, we improve upon the local determination, which fixes the deceleration parameter  to the standard ΛCDM model value of q_0=-0.55, that is, uses information from observations beyond the local universe.
First, we derive the effective calibration prior on the absolute magnitude M_B of Supernovae Ia, which can be used in cosmological analyses in order to avoid the double counting of low-redshift supernovae.
We find M_B = -19.2334 ± 0.0404 mag.
Then, we use the above M_B prior in order to obtain a determination of the local H_0 which only uses local observations and only assumes the cosmological principle, that is, large-scale homogeneity and isotropy.
This is achieved by adopting an uninformative flat prior for q_0 in the cosmographic expansion of the luminosity distance.
We use the latest Pantheon sample and find H_0= 75.35 ± 1.68  km s^-1 Mpc^-1, which features a 2.2% uncertainty, close to the 1.9% error obtained by the SH0ES Collaboration.
Our determination is at the higher tension of 4.5σ with the latest results from the Planck Collaboration that assume the ΛCDM  model.
Furthermore, we also constrain the deceleration parameter to q_0= -1.08 ± 0.29, which  
disagrees with Planck at the 1.9σ level.
These estimations only use supernovae in the redshift range 0.023≤ z≤ 0.15.








Local determination of the Hubble constant and the deceleration parameter
    Valerio Marra
    December 30, 2023
=========================================================================








§ INTRODUCTION


The ΛCDM model – the standard model of cosmology – has been extremely successful.
Assuming only General Relativity and well-understood linear perturbations about a homogeneous and isotropic background model, it accounts, with just 6 parameters, for basically all cosmological observations on a vast range of scales in space and time.
Its key ingredients are the cosmological constant, a constant of nature, and the dark matter,
a yet-undetected particle which is predicted by, e.g., supersymmetric extensions of the standard model of particle physics.
Although the theoretical basis for both the cosmological constant and dark matter may be rightfully questioned, the standard model has been keeping pragmatically its throne, thanks to its performance and simplicity.

However, since the first release of the Cosmic Microwave Background (CMB) observations by the Planck Collaboration in 2013 <cit.>, the determination of the Hubble  constant H_0 based on the standard model of cosmology started to be in tension with the model-independent determination via calibrated local Supernovae Ia (SN) by the SH0ES Collaboration in 2011 <cit.>.
The initial tension of 2.4σ worsened over the past 6 years. On one side, systematics were better understood and CMB data accumulated <cit.>. On the other side, the sample of local supernovae increased and the anchors used to calibrate them considerably improved <cit.>.
The present situation is that the two determinations of the Hubble   constant –  one by the Planck Collaboration in 2018 <cit.> and the other by the SH0ES Collaboration in late 2019 <cit.> – are now in tension at the considerable 4.2σ level.
See <cit.> for a historical overview 
and <cit.> for the present status.

Much work has been done trying to understand the implications of this tension: it is indeed, by far, the most severe problem the ΛCDM model is facing.
On one hand, the effect of the local structure – the so-called cosmic variance on H_0 – has been thoroughly studied (see <cit.> and references therein), as well as possible re-assessments of the error budget <cit.>.
On the other hand, physics beyond the standard model has been investigated, hoping that this tension could shine a light on possible alternatives to the highly tuned cosmological constant and the yet-undetected dark matter <cit.>.
However, at the moment, it is not clear which kind of physics beyond ΛCDM could solve this crisis <cit.>. 


Here, we wish to improve upon the local determination by the SH0ES Collaboration, which adopts a Dirac delta prior on the deceleration parameter, centered at the standard ΛCDM model value of  q_0=-0.55.
In <cit.> it was shown that using the broad (truncated) Gaussian prior  q_0=-0.5± 1 it is indeed possible to obtain a competitive constraint on the Hubble constant.
In what follows, first we derive the calibration prior on the absolute magnitude M_B of Supernovae Ia that was effectively used in the H_0 determination by the SH0ES Collaboration.
Then, using the M_B prior and an uninformative flat prior on q_0, we will obtain a competitive determination of the  H_0.
Indeed, our determination has an uncertainty comparable with the one by SH0ES and it is at the even higher tension of 4.5σ with the latest results from the Planck Collaboration that assume the ΛCDM  model.
We stress that our determination only uses local observations and only assumes the FLRW metric, that is, the cosmological principle, according to which the universe is homogeneous and isotropic at large scales.

This paper is organized as follows.
After introducing low-redshift cosmography in Section <ref>, we will obtain the effective calibration prior in Section <ref> and the determination of the Hubble constant in Section <ref> using the Pantheon dataset.
We will discuss our results in Section <ref> and conclude in Section <ref>.
Finally, Appendix <ref> presents details regarding the derivation of the effective calibration prior, Appendix <ref> reports the results relative to the Supercal supernova dataset, and Appendix <ref> reports the results relative to the early-2019 determination of H_0 by <cit.>.








§ COSMOGRAPHY
 

The apparent magnitude m_B^t of a supernova at redshift z is given by


    m_B^t(z) = 5log_10[d_L(z)/1 Mpc] + 25 + M_B  ,


where d_L(z) is the luminosity distance and M_B the absolute magnitude. Using a cosmographic approach within an FLRW metric – which only assumes large-scale homogeneity and isotropy – one has:


    d_L(z)    =  c z/H_0 [ 1 + (1-q_0)z/2  + O(z^2) ]  ,


where the Hubble   constant and the deceleration parameter are defined, respectively, according to:


    H_0=. ȧ(t)/a(t)|_t_0 ,
        
    q_0=. - ä(t)/H(t)^2 a(t)|_t_0 .


Cosmography is a model-independent approach in the sense that it does not assume a specific model as it is based on the Taylor expansion of the scale factor. However, this does not mean that its parameters do not contain cosmological information.
For example, the deceleration parameter is connected to the parameters of wCDM cosmologies according to:


    q_0    = Ω_ m 0/2 + 1+3 w/2Ω_ de 0
    = 1 + 3 w Ω_ de 0/2
       = 1 - 3  Ω_Λ 0/2≃ -0.55  ,


where the last equality in first line assumes spatial flatness (Ω_ m 0+Ω_ de 0=1),
the first equality in the second line assumes the ΛCDM model (w=-1), and the last one the concordance value of Ω_Λ 0≃ 0.7.



Combining now equations (<ref>-<ref>) one has:


    m_B^t     =  5 log_10[ 1 + (1-q_0)z/2 + O(z^2)  ] 
        + 5 log_10c z H_0^-1/1 Mpc  + 25 + M_B  .


We are neglecting the second-order correction that contains the jerk parameter j_0=.⃛ a/H^3 a|_t_0 because we will consider only low-redshift supernovae.
Indeed, the latest local determinations of H_0 use supernovae in the range 0.023≤ z≤ 0.15, where the minimum redshift is large enough in order to reduce the impact of cosmic variance <cit.> and the maximum redshift is small enough in order to reduce the impact of cosmology in the determination of H_0.
We have to bear in mind that the cosmographic expansion could be problematic if extended to high redshifts (z>1). Indeed, it fails to converge <cit.> and also its accuracy depends on the order adopted <cit.>.
Here, we avoid these possible issues as the maximum redshift is z=0.15 <cit.>.
Figure <ref> shows the percentage difference with respect to the computation that considers the second-order correction, together with the distribution of the Supernovae Ia that are used to determine H_0. The weighted error is 0.2%, negligible compared to the error budget to be discussed below.
Furthermore, the fact that it is safe to neglect the second-order correction implies that this analysis is valid also for spatially curved models.











§ EFFECTIVE CALIBRATION PRIOR
  

The determination of H_0 by the SH0ES Collaboration is a two-step process.
On one hand the anchors, Cepheid hosts and calibrators are combined to produce a constraint on M_B. This step depends on the astrophysical properties of the sources and is independent of cosmology.
On the other hand the Hubble-Flow supernovae of the Supercal sample <cit.> are used to probe the luminosity distance-redshift relation in order to determine H_0.
In this Section, using two different methods, we will determine the calibration prior f(M_B) which was effectively used in the latest 2019 analysis <cit.>.



 §.§ Demarginalization
 

Within a Bayesian framework, the determination of H_0 by the SH0ES Collaboration can be effectively summarized according to the following analysis:


    f(H_0,M_B | SN)    = f(H_0) f(M_B) Ł(SN | H_0, q_0,M_B) / , 
    
    f(H_0 | SN)    = ∫Ṃ_B   f(H_0,M_B | SN)   ,


where the posterior on H_0 was obtained by marginalizing over M_B.
In the equations above, f(H_0) is an improper flat prior on H_0, Ł is the likelihood and  is the evidence. SN stands for Supercal supernovae in the range 0.023≤ z≤ 0.15 <cit.>.
f(M_B) is the informative prior on the supernova absolute magnitude and is the result of the complicated calibration of the local supernovae via the cosmic distance ladder, see <cit.> for details.

The likelihood is given by:


    Ł(SN | H_0, q_0,M_B) = |2 πΣ|^-1/2  e^-1/2χ^2(H_0, q_0,M_B) ,


where the χ^2 function is:


    χ^2= { m_B,i  - m_B^t(z_i) }Σ^-1_ij{ m_B,j  - m_B^t(z_j) } ,


where Σ is the supernova covariance matrix and m_B,i are the observed apparent magnitudes at the redshifts z_i.
Following the methodology of <cit.>, we fix the nuisance parameters that control stretch and color corrections to  α=0.14, β=3.1, correct the apparent supernova magnitudes with hosts above and below log M_ stellar∼ 10 by 0.03 mag fainter and brighter, respectively, and include an intrinsic dispersion of σ_ int=0.1 mag together with a peculiar velocity error of 250 km/s.


The analysis of <cit.> fixes q_0=-0.55 in equation (<ref>) and obtains:


    H^ Re19_0 = 73.5 ± 1.4  km s^-1 Mpc^-1 .


Ref. <cit.> states that the sensitivity of H_0 to knowledge of q_0 is very low as the mean SN redshift is only 0.07.
However, q_0=-0.55 is justified only if one uses information beyond the local universe that can constrain the properties of dark matter and dark energy  (see equation (<ref>)), such as data that include supernovae at z>0.15 <cit.>.
This was necessary as there were not enough local supernovae and/or their calibration was less precise so that it was not possible to constrain well H_0 if q_0 were not fixed. However, as we will show below <cit.>, the latest supernova catalog and anchors allow us to perform an analysis which is only based in the z≤ 0.15 universe.
It is important to stress that the actual analysis of <cit.> is more complicated than the one of equation (<ref>) as it involves many intermediate steps. 


Here, in order to get the calibration prior f(M_B) we solve the integral equation obtained by demanding that equation (<ref>) gives the constraint of equation (<ref>).
Assuming a Gaussian distribution for M_B with mean M̅_B and dispersion σ_M – which is also justified a posteriori by the fact that M_B is tightly constrained by data – it is possible to marginalize analytically the 2D posterior in equation (<ref>).
The result is that H_0 is distributed according to a lognormal distribution with parameters:


    μ_ ln^ dm   = ln 10/5[ M̅_B +ln 10/5( σ_M^2+1/S_0)- S_1/S_0]  ,  
    σ_ ln^ dm   = ln 10/5√(σ_M^2 + 1/S_0) .


One can then match first and second moments of the lognormal distribution with equation (<ref>) so that one has:


    μ_ ln^ Re19≃ 4.2971   ,
        σ_ ln^ Re19≃ 0.019046  ,


where we used equations (<ref>-<ref>).
The calibration prior is then given by:


    M_B^ dm   = 5/ln 10μ_ ln^ Re19+S_1/S_0 -ln 10/5( σ_M^2+1/S_0)   , 
    σ^2_M_B^ dm   =  25/ln^2 10σ_ ln^2  Re19 - 1/S_0 .


The quantities S_0 and S_1 are defined in Appendix <ref>, where the full derivation can be found.
Using the previous equations one obtains the calibration prior  f(M_B):


    M_B^ dm = -19.2334 ± 0.0404  mag .


In the computation above we considered the second-order correction in d_L/z and fixed j_0=1 because <cit.> does so, see equation (<ref>). The second-order correction shifts the prior by a small ∼0.1σ.
The lognormal distribution on H_0 is very close to a Gaussian as shown in Figure <ref>, where the lognormal f(H_0 | SN) is compared with a Gaussian, both with mean and dispersion as in equation (<ref>).










 §.§ Deconvolution


One may obtain the calibration prior also via the following approach.
[We thank the referee for suggesting this alternative method.]
The Hubble constant and the supernova absolute magnitude are connected through the following equation <cit.>:


    log_10 H_0(q_0,j_0)   = M_B + 5a_B(q_0,j_0)  +25/5 ,


where the “intercept” a_B is obtained via the Hubble-Flow supernovae:


    a_B(z,m_B,q_0,j_0) =  log_10 [ c z f(z,q_0,j_0) ]- 1/5 m_B  , 
       f =1 + 1-q_0/2z -1-q_0-3q_0^2+j_0/6z^2+ O(z^3)  ,


where z and m_B are observed redshift and apparent magnitude of the supernovae.
In (<ref>) the absolute magnitude M_B and the intercept a_B are independent as the former is obtained via the the astrophysical observations of anchors, Cepheid hosts and calibrators, while the latter via Hubble-Flow supernovae that probe the cosmological luminosity distance-redshift relation.
In particular, fixing q_0=-0.55 and j_0=1 affects a_B but not M_B. This shows clearly that it is possible to derive the constraint on M_B by “subtracting” (or deconvolving) the contribution of a_B from the constraint on H_0.

The constraint on a_B using the Supercal sample <cit.> is <cit.>:


    a_B(q_0=-0.55,j_0=1)= 0.71273 ± 0.00176  ,


which can be obtained via <cit.>:


    a_B =  ∑_i,j(C^-1)_ij  a_B,i/∑_i,j(C^-1)_ij ,
        σ_a_B^2 = 1/∑_i,j(C^-1)_ij ,


where a_B,i=a_B(z_i,m_B,i,q_0,j_0) are the intercept measurements with covariance matrix C.
After propagating the errors on z (mostly peculiar velocity) on m_B, it is C= Σ /25, where Σ is the covariance matrix for m_B.
Regarding m_B and Σ we follow the methodology of <cit.> as described in the previous Section.


One can then proceed in two ways. The simplest approach is just to perform first-order error propagation:


    M_B^ e = 5 ( log_10 H^ Re19_0 - a_B-5)=-19.2322  , 
       σ_M^ e = √(( 5/ln 10σ_H^ Re19_0/H^ Re19_0)^2 - 25 σ_a_B^2 )=0.0404  .


Alternatively, one can make a change of variable so that the higher moments are included.
Adopting Gaussian distributions for M_B and a_B it follows that H_0 is distributed according to the lognormal distribution of equation (<ref>) with parameters:


    μ_ ln^ dc   = ln 10/5   (M̅_B+ 5 a_B +25 )  , 
    
    σ_ ln^ dc   = ln 10/5√(σ_M^2 + 25 σ_a_B^2) ,


so that the calibration prior is given by:


    M_B^ dc   = M_B^ dm +ln 10/5( σ_M^2+1/S_0)   , 
    σ^2_M_B^ dc   = σ^2_M_B^ dm ,


where we used the fact that:


    5 a_B +25 =- S_1/S_0 ,
        
    25 σ_a_B^2 = 1/S_0  .


Using the previous equations one obtains the calibration prior  f(M_B):


    M_B^ dc = -19.2326 ± 0.0404  mag ,


which confirms the validity of the first-order estimate of equation (<ref>).
We see that the demarginalization and deconvolution approaches give very similar results:
the variance of the prior is actually identical while the central value differs by a negligible amount:


    M_B^ dc - M_B^ dm≃5/ln 10( σ_H^ Re19_0/H^ Re19_0)^2 ≃ 8 × 10^-4 .


We will adopt the demarginalization result as the analysis of the next Section generalizes its methodology.




 §.§ Robustness
 

The computation of the calibration prior with the methods of the previous Sections is possible because  its determination is independent of the cosmological Hubble flow that is used to determine H_0.
In particular, M_B is independent of the value adopted for the cosmographic parameters q_0 and j_0 as their effect is to alter the luminosity distance and not M_B.
This would not be the case if SH0ES constrained Hubble flow and M_B together as, in this case, M_B would be correlated with the Hubble-flow supernovae.

Even though the methodology is well posed, there could be astrophysical or cosmological effects that could bias the determination of M_B.
A calibration bias could be produced by the coherent flow of the local supernovae, akin to cosmic variance <cit.>.
However, both the determination of equation (<ref>) and the Supercal dataset have already been  corrected for such effects using velocity fields derived from galaxy density fields <cit.>.
Another bias could be due to correlations between supernovae and environmental properties such as star formation rate and metallicity. Ref. <cit.> concluded that these correlations are not significant so that they should not impact equation (<ref>) either <cit.>. 
Therefore, we conclude that our determination of the M_B prior is robust and correctly encodes the supernova calibration via the distance ladder.




 §.§ Calibration prior in cosmological analyses
 


The calibration prior on M_B of equation (<ref>) can be meaningfully used in cosmological analyses instead of the corresponding H_0 determination.
Indeed, there are 175 supernovae in common between the Supercal and Pantheon datasets in the range 0.023≤ z≤ 0.15, and, in the standard analysis, these supernovae are used twice: once for the H_0 determination and once when constraining the cosmological parameters.
This induces a covariance between H_0 and the other parameters which could bias  cosmological inference.
The advantage of adopting  the calibration prior is that these 175 supernovae are used only once, avoiding any potential bias.











§ DETERMINATION OF LOCAL H_0
  




Here, we wish to obtain a determination of H_0 that is based solely on the local universe (z ≤ 0.15) and that is truly independent of any cosmological assumptions besides the cosmological principle and the astrophysics of anchors, Cepheid hosts and calibrators.

We achieve the latter by considering 
an uninformative flat prior on q_0 and marginalizing the posterior also over  q_0, that is:


    f(H_0,q_0, M_B | SN) = f(H_0)f(q_0) f(M_B) Ł/ , 
       f(H_0 | SN) = ∫Ṃ_B q̣_0  f(H_0, q_0,M_B  | SN)  ,


where f(q_0) is the flat ignorance prior on q_0.
The method of <cit.> is recovered if we set f(q_0)= δ(q_0 +0.55).
In this analysis we adopt the more recent Pantheon supernova sample <cit.>, which features 237 supernovae in the redshift range 0.023≤ z ≤ 0.15, see Figure <ref>.
In this case the covariance matrix and the apparent supernova magnitudes do not depend on the nuisance parameters that control stretch and color corrections, and the correction for the apparent magnitudes with hosts above and below log M_ stellar∼ 10 is already applied. Also, the covariance matrix includes all the uncertainties.
The results relative to the Supercal sample are given in Appendix <ref>.


We have used  <cit.> – an open-source  sampler for Markov chain Monte Carlo (MCMC) – to sample the posterior and  <cit.> for analyzing the chains.
The result of this analysis is presented in Figure <ref>. We obtain:


    H_0^ loc   =  75.35 ± 1.68  km s^-1 Mpc^-1 , 
    
    q_0^ loc   = -1.08 ± 0.29  .


The latest CMB-only constraint from the Planck Collaboration <cit.> is:


    H^ P18_0 = 67.36 ± 0.54  km s^-1 Mpc^-1 ,


and the tension with our determination now reaches 4.5 σ:


    |H_0^ loc - H^ P18_0|/√(σ_H_0^2+σ^2_H^ P18_0)≃ 4.52  .







As it is clear from Figure <ref>, low-redshift supernovae are able to constrain not only H_0 but also q_0, whose distribution is peaked at values lower than the standard model one of q_0=-0.55.
However, the deceleration parameter is not yet tightly constrained and the tension between our determination and the one by Planck is at the 1.9σ level.
Note that the ΛCDM model cannot give values of q_0 below -1.

Finally, it is interesting to compare directly the marginalized 2D posterior on H_0 and q_0 from CMB and low-z supernova observations. We obtained the CMB posterior using the MCMC chain from the CMB-only  Planck 2018 analysis available at http://www.esa.int/Planckesa.int/Planck.
The tension in the H_0–q_0 plane is also presented in Figure <ref>.

In order to quantity the tension we adopt the index of inconsistency (IOI) <cit.> which directly generalizes the estimator of equation (<ref>). We find:


    √(2 IOI)   ≡√(δ^T (C_ loc + C_ P18)^-1δ )≃ 4.54  , 
    δ   = {H_0^ loc-H_0^ P18,   q_0^ loc -q_0^ P18} ,


where C are the covariance matrixes on H_0 and q_0 from the analysis of Figure <ref> and δ is the difference vector.
Note that these estimators assume Gaussianity and that the posteriors on H_0 and q_0 are very close to Gaussian. The tension in the H_0–q_0 plane is again very strong: 4.5 σ.






§ DISCUSSION
  

Our determination of the Hubble constant of equation (<ref>) is based on the calibration prior of equation (<ref>) and on the low-redshift cosmographic expansion of equation (<ref>).
While the validity of the calibration prior depends on the standardizable nature of supernovae Ia, a process that includes corrections due to color, stretch and host-galaxy mass, the validity of the cosmographic analysis is solely based on the approximation that the FLRW metric provides a good description of our universe at large scales.
While it is possible to test directly this hypothesis <cit.>, the validity of the FLRW metric is a direct consequence of assuming the Cosmological principle, according to which the universe is homogeneous and isotropic at large scales.
Therefore, as far as cosmological assumptions are concerned, the determination of equation (<ref>) only assumes the Cosmological principle.


While our goal here was to improve the local determination of H_0 by considering a flat prior on q_0, it is also possible to perform a more thorough re-analysis of the cosmic distance ladder and consider the impact of non-Gaussianities <cit.>, hyper-parameters <cit.> and a blinded pipeline <cit.>.
In particular, our results are in agreement with the findings of <cit.>, which, using a Bayesian hierarchical model, also found that data prefer a lower value of q_0≈ -1 and so – owing to their anti-correlation – a higher value of H_0.
In particular, <cit.> adopted a broad (truncated) Gaussian prior q_0=-0.5± 1, similar to our improper flat prior, and obtained a value of H_0 which is ∼1 km s^-1 Mpc^-1 higher than their baseline H_0 value.
This is in good agreement with our results relative to the Supercal sample, the one adopted by <cit.>.
Indeed the determination of equation (<ref>) is ∼1 km s^-1 Mpc^-1 than the H^ Re19_0 of equation (<ref>).


It is also possible to extend the cosmic ladder to BAO and CMB observations.
Ref. <cit.> constrained H_0 via a cosmographic inverse ladder approach,
using BAO measurements to break the degeneracy between M_B and H_0, and finding a  value of the Hubble constant very close to the Planck one: H^ M18_0 = 67.8 ± 1.3  km s^-1 Mpc^-1.
The cosmographic inverse ladder approach propagates the prior on the sound horizon from CMB to the BAO scale. However, the calibration of M_B depends on the luminosity distance-redshift relation of supernovae both at low and high redshift.
This introduces correlations between M_B and the properties of dark energy, which can explain why our cosmology-independent determination of q_0 differs from the constraint q^ M18_0=-0.37 ± 0.15 obtained by <cit.>.



Finally, a low local value of q_0≈ -1 was also obtained in <cit.>, where the cosmic ladder was extended to BAO and CMB observations using the distance-duality relation instead of a high-redshift cosmographic expansion.









§ CONCLUSIONS
  


We obtained the effective  local calibration prior on the absolute magnitude M_B of Supernovae Ia, which we used to obtain a determination of the Hubble constant from the local universe that assumes only large-scale homogeneity and isotropy: H_0= 75.35 ± 1.68  km s^-1 Mpc^-1.
Our determination uses the latest Pantheon sample.
As shown in Figure <ref>, the allowed values are in stronger tension with what the CMB predicts if the standard model of cosmology is valid. The H_0 crisis is approaching the 5σ level.
We stress that our determination relies only on supernovae in the redshift range 0.023≤ z≤ 0.15.
If on one hand this limits the available information as there are many more supernovae at higher redshifts, on the other hand this restricted redshift range makes this analysis local.

We would like to stress that the calibration prior on M_B = -19.2334 ± 0.0404 mag can be meaningfully used in cosmological analyses, instead of the corresponding H_0 determination, in order to avoid the double counting of low-redshift supernovae.



The fact that  both Hubble constant and deceleration parameter differ
seems to suggest that a new phenomenon is affecting both the sound horizon and the  recent expansion rate of the universe.
Early dark energy may be a promising candidate <cit.>.
The increased statistical significance of our determination of  H_0 reinforces this suggestion.







DC thanks CAPES for financial support.
VM thanks CNPq and FAPES for partial financial support.









§ DERIVATION OF THE DEMARGINALIZATION
  

Here, we explain how we obtained equations (<ref>-<ref>).
First, we define:


    m̃_B^t ≡ m_B^t -M_B + logh  ,


where we used the shorthand notation logh =5 log_10 H_0.
We can then rewrite equation (<ref>) as:


    χ^2   = { W_i+ logh- M_B  }Σ^-1_ij{ W_i +logh- M_B  }
       = S_0(M_B - logh - S_1/S_0)^2  + const  ,


where const consists of all the terms that do not depend on H_0 and M_B, and where we defined:


    W_i    =m_B,i - m̃_B^t(z_i)  , 
    
    V_i    = 1  , 
    
    S_0   =  V ·Σ^-1· V^T , 
    
    S_1   =  W ·Σ^-1· V^T .


The unnormalized posterior of equation (<ref>) is then:


    f(H_0,M_B | SN) ∝ e^-(M-M̅_B)^2/2 σ_M^2  e^-1/2χ^2(H_0,M_B) ,


where an improper flat prior on H_0 was adopted.
Equation (<ref>) can then be integrated over M_B so that the normalized posterior on H_0 is distributed according to the lognormal distribution:


    f(H_0 | SN)    = ∫Ṃ_B   f(H_0,M_B | SN) 
       =  1/H_0 √(2π)σ_ lnexp - (ln H_0-μ_ ln)^2/2σ_ ln ^2,


with parameters:


    μ_ ln^ dm   = ln 10/5[ M̅_B +ln 10/5( σ_M^2+1/S_0)- S_1/S_0]  ,  
    σ_ ln^ dm   = ln 10/5√(σ_M^2 + 1/S_0) ,


Finally, mean and variance of the lognormal distribution are:


    H̅_0     = e^μ_ ln + 1/2σ_ ln^2 , 
    σ^2_H_0   =(e^σ_ ln^2 - 1) e^2μ_ ln + σ_ ln^2 .









§ ANALYSIS WITH THE SUPERCAL SAMPLE
  



Here, we present the results using the Supercal supernova catalog <cit.> that was used by <cit.>. It features 217 supernovae in the redshift range 0.023≤ z ≤ 0.15, see Figure <ref>.
Figure <ref> shows the results of this analysis, which are compared with the ones of Section <ref> that use the more recent Pantheon sample.
We obtain:


    H_0^ loc   =  74.62 ± 1.54  km s^-1 Mpc^-1 ,  
    
    q_0^ loc   = -0.90 ± 0.22  ,


which are in tension with the determinations of H_0 and q_0 by Planck at the  4.5σ and 1.7σ levels, respectively. The tension in the H_0–q_0 plane is 4.5σ.
There is good agreement between the results that use the Supercal and Pantheon datasets.



The previous results were obtained using the methodology of <cit.>, that is, fixing the nuisance parameters that control stretch and color corrections to  α=0.14, β=3.1, correcting the apparent supernova magnitudes with hosts above and below log M_ stellar∼ 10 by 0.03 mag fainter and brighter, respectively, and including an intrinsic dispersion of σ_ int=0.1 mag together with a peculiar velocity error of 250 km/s.
If instead we adopt a flat prior on α and β we obtain basically the same results:  


    H_0^ loc   =  74.85 ± 1.57  km s^-1 Mpc^-1 ,  
    
    q_0^ loc   = -0.98 ± 0.24  .










§ ANALYSIS RELATIVE TO RIESS ET AL. 2019
  

Here, we report the results relative to the early-2019 measurement by Riess et al. (2019) <cit.>:


    H^ R19_0 = 74.03 ± 1.42  km s^-1 Mpc^-1 .


The calibration prior is:


    M_B^ dm = -19.2178 ± 0.0407  mag ,


which gives the determinations:


    H_0^ loc   =  75.89 ± 1.70  km s^-1 Mpc^-1 , 
    
    q_0^ loc   = -1.08 ± 0.29  ,


which are in tension with the corresponding CMB-only constraints from the Planck Collaboration at the 4.8 σ and 1.9σ levels, respectively. The tension in the H_0–q_0 plane is at the 4.8 σ level.



utphys



