Free Access
Issue
A&A
Volume 592, August 2016
Article Number A152
Number of page(s) 11
Section Cosmology (including clusters of galaxies)
DOI https://doi.org/10.1051/0004-6361/201628572
Published online 19 August 2016

© ESO, 2016

1. Introduction

Supernova (SN) data have provided evidence that distant sources (z> 0.3) appear dimmer than predicted in a universe with matter only, in comparison with nearby sources. This dimming led to the interpretation that the Universe is expanding in an accelerating fashion. This acceleration has been attributed to a dark energy component, whose simplest solution is a vacuum energy or equivalently a cosmological constant. However, this dimming could be caused by a variation of any other component that affects the luminosity distance, namely by inhomogeneities in the energy densities or in the cosmic expansion (Amendola et al. 2013).

In order to explain the underlying mechanism of the late-time accelerated expansion of the universe, one possibility is to give up of the cosmological principle and allow for an anisotropic expansion of the universe. The idea that we live in a locally underdense region, hence creating a “Hubble bubble”, can explain the cosmic acceleration at late times, subject to specific conditions (Zehavi et al. 1998; Caldwell & Stebbins 2008). From the vast literature on the cosmological principle, it has been found that most cosmological observations can accommodate violations of the cosmological principle (see for example Zibin et al. 2008; Tomita & Inoue 2009; Marra & Pääkkönen 2010; Marra & Notari 2011; Moss et al. 2011). From the theoretical point of view, violations of the cosmological principle can be explained in the context of Lemaître-Tolman-Bondi (LTB) void models. These void models are spherically symmetric and radially inhomogeneous, and can mimic the cosmic expansion of the concordance ΛCDM (Lan et al. 2010; Liu & Zhang 2014).

In Carvalho & Marques (2016), we introduced a method to probe the inhomogeneity of the large-scale structure by measuring the angular distribution in the cosmological parameters that affect the luminosity distance, using SN data. Variation in the cosmological parameters across pixels in the sky implies inhomogeneity in the cosmic expansion. This inhomogeneity was then used to measure the extra component of cosmic acceleration predicted by backreaction, which derives from averaging over an inhomogeneous ensemble of homogeneous pixels, each pixel expanding at a different rate. However, this measurement presupposed an a priori dark energy component in each pixel. In order to investigate alternative interpretations of the cosmic acceleration, it is conceptually more consistent to use a parametrisation that does not assume a specific component as the cause of acceleration (Turner & Riess 2002).

In this manuscript, we use the luminosity distance expressed in terms of time derivatives of the scale factor a(z), in particular the Hubble parameter H0, the deceleration parameter q0 and the jerk parameter j0. Instead of assuming an equation of state and inferring the evolution of the scale factor via the Friedmann equation, the reasoning is to take the data on the scale factor and infer a cosmological equation of state via the Friedmann equation. This parametrisation records the cosmic expansion without regard to its cause, thus being independent of the cosmological equation of state adequate to test interpretations of the cosmic acceleration alternative to the cosmological constant. These parameters can be related to the Taylor expansion of the cosmological equation of state about the present values { ρ0,p0 } expressed up to linear order as p=p0+κ0(ρρ0)+O[(ρρ0)2],\begin{eqnarray} p=p_{0}+\kappa_{0}(\rho-\rho_{0})+O\left[(\rho-\rho_{0})^2\right], \end{eqnarray}(1)hence yielding information about the present values of w0 and κ0 defined as (Visser 2004) w0=pρ|0w0(H0,q0),κ0=dpdρ|0κ0(H0,q0,j0).\begin{eqnarray} w_{0}={p\over \rho}\Big \vert_{0}\equiv w_{0}(H_{0},q_{0}), \quad \kappa_{0}={{\rm d}p\over {\rm d}\rho}\Big \vert_{0}\equiv \kappa_{0}(H_{0},q_{0},j_{0}). \end{eqnarray}(2)Whereas w0 (and hence q0) contains information about the present value of p/ρ,κ0 (and hence j0) contains information about how w0 can evolve.

The purpose of this manuscript is to reapply the method first presented in Carvalho & Marques (2016) to estimate the parameters { H0,q0,j0 } , instead of { H0MΛ } , by fitting the luminosity distance to SN data, and to compare the results from the two estimations in view of an interpretation of the cosmic acceleration independent of a particular energy content of the Universe. (Although the comparisons are made for the case that Ωκ = 1−ΩM−ΩΛ, for completion we also include the results for the case where Ωκ = 0.) Some studies have used the kinematic parametrization to measure inhomogeneities from SN data (see e.g. Schwarz & Weinhorst 2007; Kalus et al. 2013). Most studies, however, have aimed at finding hemispherical anisotropies assuming a ΛCDM energy content (e.g. Blomqvist et al. 2010; Mariano & Perivolaropoulos 2012; Heneka et al. 2014; Jiménez et al. 2015; Bengaly et al. 2015a; Javanmardi et al. 2015; Migkas & Plionis 2016 and references therein).

This paper is organised as follows. In Sect. 2 we describe the data and estimate the observables by performing both a global and a local parameter estimation, obtaining fiducial values and maps respectively for the estimated parameters. We introduce an inhomogeneity test by rotating the supernova subsampling per pixel. In Sect. 3 we compute the power spectrum of the maps of the parameters using two methods for the noise bias removal. In Sect. 4, for the same toy model of backreaction used in Carvalho & Marques (2016), we compute the average values of { q0,j0 } and discuss possible cosmological implications. In Sect. 5 we draw conclusions.

2. Parameter estimation

We use the type Ia supernova sample compiled by the Joint Light-curve Analysis (JLA) collaborative effort (Betoule et al. 2014) from different supernova surveys, totalling NSNe = 740 supernovae (SNe) with redshift z ∈ [0.010,1.30] and distributed on the sky according to Fig. 1 in Carvalho & Marques (2016)1.

We use the luminosity distance dL with a(z) expressed as a Taylor expansion about the present value a(z0) (Visser 2004) dL(z;H0,q0,j0)=czH0×[1+12(1q0)z16(1q03q02+j0+kc2H02a02)z2\begin{eqnarray} &&d_{\rm L}(z;H_0,q_0,j_0)={cz\over H_{0}} \nonumber\\ &&\times \,\, \Biggl[1+{1\over 2}(1-q_0)z -{1\over 6}\left(1-q_0-3q_0^2+j_0+{kc^2\over{H_0^2 a_0^2}}\right)z^2\nonumber\\ &&+O(z^3)\Biggr], \label{eqn:dl} \end{eqnarray}(3)to estimate the cosmological parameters that minimise the chi-square of the fit of the theoretical distance modulus μtheo(z;H0,q0,j0)=5log[dL(z;H0,q0,j0)]+25,\begin{eqnarray} \mu_{\rm theo}(z;H_{0},q_0,j_0) =5\log[d_{\rm L}(z;H_{0},q_0,j_0)]+25 \label{eqn:mu_theo} , \end{eqnarray}(4)(dL in units of Mpc) computed for the trial values of { H0,q0,j0 } , to the measured distance modulus μdata(z)=mB(z)(MBαx1+βc),\begin{eqnarray} \mu_{\rm data}(z)=m_{B}(z)-(M_{B}-\alpha x_{1}+\beta c), \label{eqn:mu_data} \end{eqnarray}(5)computed for the light-curve parameters { mB,x1,c } estimated from each SN’s observed magnitude, and for MB = −19.05 ± 0.02,α = 0.141 ± 0.006 and β = 3.101 ± 0.075 estimated for all SNe (Betoule et al. 2014).

We recall that the { H0,q0 } expansion in Eq. (3) (i.e. up to the second term in the right-hand side) does not allow for an estimation of q0 that is both accurate and precise (Neben & Turner 2013). This is because for low redshifts where the Taylor expansion is most accurate, there is poor leverage for a precise estimation of q0; conversely, for high redshifts where the leverage is larger and allows for a more precise estimation, the Taylor expansion is less accurate. Considering the {H0,q0,j0} expansion in Eq. (3) (i.e. all the terms in the right-hand side), then we expect that the uncertainty will be pushed to the estimation of j0. An expansion in higher-order time derivatives of the expansion factor was considered in Aviles et al. (2012). An exhaustive comparison of luminosity distances can be found in Cattoën & Visser (2008).

Before proceeding further, we check the validity range of Eq. (3) using a spatially flat ΛCDM model as reference model. In this case, we have q0 = (ΩM/ 2)−ΩΛ and j0 = ΩM + ΩΛ = 1. We compute the relative difference between the luminosity distance dL(z) in Eq. (3) and the theoretical luminosity distance in the ΛCDM model dL(z)=c(1+z)H00zdzE(z),\begin{eqnarray} d_{\rm L}(z)=\frac{c(1+z)}{H_{0}} \int_{0}^{z}\frac{{\rm d}z}{E(z)}, \label{eqn:dl_z} \end{eqnarray}(6)where E(z) = [ΩM(1 + z)3 + ΩΛ] 1 / 2 for ΩM = 0.308 and H0 = 67.8 km s-1 Mpc-1 (Planck Collaboration XIII 2016). (See solid line in Fig. 1.) We find that, for z ∈ [0.1,1.1] , the relative difference lies in the interval [−0.02%,−3.9%] ; conversely, for z ∈ ] 1.1,1.30] , the relative difference can reach about −6.4%. We also compute the relative difference between the luminosity distance dL(y), with yz/ (1 + z), in Cattoën & Visser (2008) and the theoretical luminosity distance in the ΛCDM model. (See dashed line in Fig. 1.) In this case we find that, for z ∈ [0.1,1.1] , the relative difference lies in the interval [−0.02%,−2.1%] ; conversely, for z ∈ ] 1.1,1.30] , the relative difference is about −2.6%. We also observe that, in the intermediate redshift range z ∈ [0.1,0.75] where most SNe were detected, Eq. (3) performs slightly better than dL(y); conversely, dL(y) clearly performs better for z> 1. We also compute the corresponding relative differences in μ, which is the quantity that we use in the subsequent estimation and which depends on dL logarithmically according to Eq. (4). (See inset plot of Fig. 1.) We observe that the relative differences in μ are less than −0.4%, which implies that the differences in the luminosity density translate in a negligible effect on our results.

thumbnail Fig. 1

Relative difference between expressions of the luminosity distance. The solid line shows the relative difference between the luminosity distance dL(z) in Eq. (3) and the theoretical luminosity distance in the ΛCDM model in Eq. (6) for ΩM = 0.308 and H0 = 67.8 km s-1 Mpc-1 (Planck Collaboration XIII 2016). The dashed line shows the relative difference between the luminosity distance dL(y), with yz/ (1 + z), in Cattoën & Visser (2008) and the theoretical luminosity distance in the ΛCDM model. The inset plot shows the corresponding differences in μ, given by Eq. (4), which are indicated by the corresponding line styles.

We now proceed to estimate the parameters { H0,q0,j0 } using the method described in Appendix A of Carvalho & Marques (2016). We assume that at the present epoch c/H0a0 ≪ 1, regardless of the value of k, which means that by setting the prior k = 0 and estimating j0 we are actually estimating j0 + kc2/ (H0a0)2 (Neben & Turner 2013). For a model consisting of an incoherent mixture of matter, where each component is described by an equation of state pi = wiρi, this assumption implies that j0q0(1 + 2q0). We also assume that the measurements of the magnitude (and consequently of μ) at different redshifts are independent, which is equivalent to assuming that the parameters’ distribution is Gaussian as a first approximation. We compute the error σμdata by assuming uncorrelated errors for the observables and by error propagating according to Eq. (5).

2.1. Global parameter estimation

We first perform a parameter estimation using the complete SN sample, called the global estimation and corresponding to npixel= 1 pixel, from which we estimate the maximum likelihood values for the parameters xi = { H0,q0,j0 } . (Whenever unspecified, H0 is measured in units of km s-1 Mpc-1.)

For the global estimation, we ran nj = 10 realizations of Markov chains of length 104. The starting point of each realization is randomly generated. By removing the first 20% of entries in the chains to keep the burnt–in phase only, thinning down the remaining 80% to a half by removing one of each consecutive entry in order to remove correlations within the chain, and finally averaging over the various chains, we obtain the following results: xifid={H0fid,q0fid,j0fid}={71.06,0.540,0.533}±{0.46,0.094,0.503}\hbox{$x_{i}^{\text{fid}}=\{H_0^{\text{fid}},q_0^{\text{fid}},j_0^{\text{fid}}\}=\{71.06,-0.540,0.533\}\pm \{0.46,0.094,0.503\}$} (see Table 1). We will use these results as the fiducial values in the subsequent calculations. The errors contain the dispersion in each chain and the dispersion among the averages of the different chains added in quadrature.

The measurements of { H0,q0,j0 } are consistent with the measurements of { H0,q0 = ΩM/ 2−ΩΛ,j0 = ΩM + ΩΛ } obtained in Carvalho & Marques (2016), which we include in Table 1 for convenience. The measurements of { H0,q0 } are equally precise, favouring q0< 0 at over 5σ. The measurement of j0 is comparatively less precise, as also observed in Neben & Turner (2013) in the context of the kinematic parametrisation, but nonetheless it favours j0> 0 at 1σ. These results imply that, in the interval z ∈ [0.01,1.30] covered by the SN sample, the cosmic expansion accelerated and that previously it had decelerated, meaning that there was a time when the acceleration changed sign, which supports the evidence found in Riess et al. (2004).

For illustration, in the left panel of Fig. 2 we plot μtheo computed for the estimated values xifid\hbox{$x_{i}^{\text{fid}}$} as the solid black line. In the right panel we plot the Markov chain of one realization.

We recall that Riess et al. (2004) found q0 ∈ [−1.3,−0.2] (using the gold sample) which implies j0 ∈ [−0.3,5.9] , and q0 ∈ [−1.4,−0.3] (using the combined gold+silver sample) which implies j0 ∈ [−0.1,6.4] . Moreover, Caldwell & Kamionkowski (2004) found q0 ∈ [−1.1,−0.2] and j0 ∈ [−0.5,3.9] . More recent results include Riess et al.’s H0 = 73.8 ± 2.4 for q0 = −0.55 and j0 = 1 using the Hubble Space Telescope set (Riess et al. 2011), and Neben & Turner’s q0 = −0.64 ± 0.14 and j0 = 1.4 ± 0.8 using the Constitution set (Neben & Turner 2013). (See Table 1.)

Table 1

Values for the parameters estimated from the JLA type Ia SN sample.

thumbnail Fig. 2

Left panel: distance modulus as a function of redshift. The dots are coloured according to the original SN survey, as described in Betoule et al. (2014), and indicate the distance modulus μdata computed from each SN in the sample. Blue: several low-z. Green: SDSS-II. Yellow: SNLS. Red: HST. The black line indicates the distance modulus computed theoretically for the values { H0,q0,j0 } estimated using the complete sample. The black dashed and black dotted lines indicate the distance modulus for two extreme cases. Right panel: Markov chain of one realization. The dots indicate the entries in the Markov chain. The position in the graph indicates the values of (q0,j0) of each entry according to the axes, whereas the colour indicates the value of H0 as detailed in the side colour bar. The dashed line represents j0 = q0(1 + 2q0), which separates between flat and curved space–time.

2.2. Local parameter estimation

We then divide the SNe over a pixelated map of the sky with pixels of equal surface area according to the HEALPix pixelation (Gorski et al. 2005), and perform a parameter estimation using the SN subsample that falls into each pixel, called the local estimation. The number of pixels that guarantees non–empty SN subsamples in all pixels is npixel= 12 pixels. The number of SNe in each pixel is indicated in Fig. 3 in Carvalho & Marques (2016). Each pixel is assumed to be described by a Friedmann–Lemaître–Roberston–Walker metric so that the full sky is an inhomogeneous ensemble of disjoint, locally homogeneous regions.

For the local estimation, we distinguish two cases as introduced in Carvalho & Marques (2016): the “Cosmic variance” estimation where in each subsample we use the original redshifts and positions, and the “Shuffle SNe” estimation where in each subsample we randomly shuffle the SNe in redshift while keeping the original positions in the sky. The “Shuffle SNe” estimation is hypothesised to be a measure of the noise bias due to the inhomogeneous coverage of the sky by the SN surveys, from which there results inhomogeneity in the SN subsampling. We model this noise bias by the local estimation obtained from randomizing the dependence of redshift with position, while keeping the original SN positions in the sky, and then averaging over the various randomizations.

In each pixel k and for each case of local estimation, we ran nj = 300 realizations of Markov chains of length 104 and repeated the procedure described above for the global estimation. From the local estimation, there result maps for each estimated parameter xi with the same pixelation as the SN subsamples, denoted by \hbox{$\bar x_{ik}=\big<x_{ijk}\big>_j$} in the “Cosmic variance” estimation and by biasik=xijkbiasj\hbox{$\bar x_{ik}^{\text{bias}}=\big<x_{ijk}^{\text{bias}}\big>_j$} in the “Shuffle SNe” estimation, where the brackets denote averaging over the j realizations. Subtracting the noise bias biasik\hbox{$\bar x_{ik}^{\text{bias}}$} off \hbox{$\bar x_{ik},$} we obtain unbiased maps unbiasikikbiasik+xifid.\hbox{$\bar x_{ik}^{\text{unbias}}\equiv \bar x_{ik}-\bar x_{ik}^{\text{bias}}+x_i^{\text{fid}}.$} For convenience, we also compare the local estimation with the fiducial values by defining the difference maps Δik=ikxifid,\hbox{$\Delta \bar x_{ik}=\bar x_{ik}-x_{i}^{\text{fid}},$} and the unbiased difference maps as Δunbiasikunbiasikxifid=ikbiasik.\hbox{$\Delta \bar x_{ik}^{\text{unbias}}\equiv \bar x_{ik}^{\text{unbias}}-x_{i}^{\text{fid}}=\bar x_{ik}-\bar x_{ik}^{\text{bias}}.$}

The difference maps are shown in Fig. 3 (top panel sets), before (left panel set) and after (right panel set) the noise bias subtraction. Comparing the difference maps with the fiducial values, we measure fluctuations of order 0.1–5% for H0, 1–187% for q0 and 1–184% for j0 before the noise bias subtraction; after the noise bias subtraction, we measure fluctuations of order 0.1–7% for H0, 1–136% for q0, and 1–221% for j0. We recall that Wiegand & Schwarz (2012) found δH0/H0<~5% using galaxy surveys, which is consistent with our results.

For comparison, in the bottom panels of Fig. 3, we reproduce the difference maps of { H0,q0 = ΩM/ 2−ΩΛ,j0 = ΩM + ΩΛ } from the { H0MΛ } estimation in Carvalho & Marques (2016). We recall that the { H0MΛ } estimation yielded fluctuations about the fiducial values of order 0.1–3% for H0, 0.1–63% for q0, and 0.001–34% for j0 before the noise bias subtraction, and fluctuations about the fiducial values of order 0.1–5% for H0, 0.1–32% for q0, and 1–27% for j0 after the noise bias subtraction. This amounts to an increase in the fluctuations by a factor of four and two orders of magnitude respectively before and after the noise bias subtraction, between the { H0MΛ } and the { H0,q0,j0 } estimations. The increase in the fluctuations supports the reasoning that accuracy in the estimation of { H0,q0 } is compromised by the inclusion of j0 as an estimated parameter. However, this is also a consequence of using a parametrisation that is independent of the cosmological equation of state, hence of assuming less about the cause of acceleration.

For further comparison, we compute the largest fluctuation in the maps, defined as \hbox{$\text{max}(\Delta \bar x_{ik})\equiv \text{max}(\bar x_{ik})-\text{min}(\bar x_{ik}),$} which yields max( { ΔH0q0j0 } ) = { 5.58,1.27,1.90 } and max( { ΔH0q0j0 } ) = { 6.40,1.31,1.97 } , before and after the noise bias subtraction respectively. We recall that the { H0MΛ } estimation yielded max( { ΔH0q0j0 } ) = { 3.84,0.452,0.440 } and max( { ΔH0q0j0 } ) = { 4.40,0.416,0.450 } , before and after the noise bias subtraction respectively. In comparison with the results from other SN studies, namely Kalus et al.’s 2.0 < max(ΔH0) < 3.4 km s-1 Mpc-1 from the Union 2 data for fixed q0 = −0.601 (Kalus et al. 2013), and Bengaly et al.’s max(ΔH0) = { 4.6,5.6 } km s-1 Mpc-1 and max(Δq0) = { 2.56,1.62 } respectively from the Union 2.1 data and the JLA data (Bengaly et al. 2015a), we observe that: a) our results for H0 remain intermediate between the previous results; b) our results for q0 are intermediate between those of Bengaly et al. (2015a) for the two data sets; and c) our results for j0 are to the best of our knowledge the first measurements of the kind.

A visual inspection of the pixels through which the Galactic plane crosses does not reveal consistently larger/smaller values of the fluctuations, as also noted in Carvalho & Marques (2016), thus suggesting negligible correlation with the Galactic plane in comparison with the measurement errors (Neben & Turner 2013).

thumbnail Fig. 3

Fluctuations of the estimated parameters as a function of the number of SNe per pixel, before (left panel sets) and after (right panel sets) noise bias removal. Each panel set consists of five panels. At each pixel we plot: a) in the first panel, the number of SNe; b) in the second panel, the value of the parameters { q0,j0 } respectively as the solid red and solid blue lines, with the dashed red and dashed blue lines indicating the corresponding fiducial values; c) in the third panel, the standard deviation { σq0,σj0 } respectively as solid red and solid blue lines; d) in the fourth plot, the value of H0 as the solid black line, with the dashed black line marking the fiducial value; e) in the fifth panel, the standard deviation σH0 as the solid black line. Top panel sets: results from the { H0,q0,j0 } estimation for one pixelation. Centre panel sets: results from the { H0,q0,j0 } estimation for the average over the different pixelations. Bottom panel sets: results from the { H0MΛ } estimation for Ωκ = 1−ΩM−ΩΛ.

In order to verify the dependence of the results on the subsampling of the SN surveys per pixel, we perform a test of the effect of rotating the HEALPix pixelation. In particular, starting from the current pixelation and keeping the pixel size constant (i.e. keeping the same number of pixels on the sky), we consider n divisions of each pixel and rotate the pixelation (np−1) times in the same direction, the npth rotation corresponding to the initial pixelation. At each such rotation by a fraction of 1 /np of the pixel, we obtain a different pixelation of the sky and hence a different subsampling of the SN surveys per pixel, totalling np different pixelations. We choose np = 9 to guarantee a sufficient number of rotations and simultaneously little degeneracy among the rotations.

For each such pixelation p, we perform a local parameter estimation as described above, consisting of a “Cosmic covariance” estimation and a “Shuffle SNe” estimation, and hence obtaining an unbiased parameter estimation in each pixel k,unbiasikp.\hbox{$\bar x^{\text{unbias}}_{ikp}.$} We then compute the mean unbiased maps averaged over the different pixelations as unbiasik=unbiasikp,\hbox{$\bar x_{ik}^{\text{unbias}}=\big<\bar x_{ikp}^{\text{unbias}}\big>,$} weighted by the inverse of the variance in each pixel, over the different pixelations. From the average over the different pixelations there results a mean pixel derived by averaging the (n−1) rotations of the same pixel in the initial pixelation. For convenience, we define the unbiased difference maps as Δunbiasikunbiasikxifid.\hbox{$\Delta \bar x_{ik}^{\text{unbias}}\equiv \bar x_{ik}^{\text{unbias}}-x_{i}^{\text{fid}}.$} The resulting difference maps are shown in Fig. 3 centre panel sets, before (left panel set) and after (right panel set) the noise bias subtraction. Comparing the difference maps with the fiducial values, we measure fluctuations of order 0.1–7% for H0, 0.1–180% for q0, and 1–116% for j0 before the noise bias subtraction; after the noise bias subtraction we measure fluctuations of order 0.1–7% for H0, 0.1–112% for q0, and 1–135% for j0. The noise bias removal brings the pixel values closer to the corresponding values estimated from the complete sample, hence decreasing the fluctuations across the sky, albeit increasing slightly the fluctuations in j0. Simultaneously, it increases the error by up to {10,65,5}% respectively. These results also seem to indicate the validity of biasik\hbox{$\bar x_{ik}^{\text{bias}}$} as a measure of the noise bias due to the inhomogeneous SN sampling.

Similarly, the largest fluctuation in the maps yields max({ΔH0q0j0 } ) = {6.01,1.19,1.22} and max({ΔH0q0j0}) = {5.40,1.05,1.27} , before and after the noise bias subtraction.

Hence by averaging over different pixelations of the sample, we obtain a decrease in the fluctuations of the parameters across the sky. However, the mean fluctuations are still larger (by a factor of four and two orders of magnitude respectively before and after the noise bias subtraction) than those obtained in the { H0MΛ } estimation. In the subsequent calculations, we will use the parameters’ maps obtained from all np pixelations.

3. Power spectra

In order to probe the distribution of the fluctuations by scale, we compute the angular power spectrum of the difference maps normalised to the corresponding fiducial parameter value. The equations below follow closely those in Carvalho & Marques (2016) except for the extra degree of complexity introduced by the different pixelations. In particular, for each parameter xi and for each realization j of the “Cosmic variance” local estimation, we define the map with value δxijkp(xijkpxifid)/xifid\hbox{$\delta x_{ijkp}\equiv (x_{ijkp}-x^{\text{fid}}_{i})/x^{\text{fid}}_{i}$} at each pixel k, for each pixelation p. An estimator of the power spectrum is ℓ,ijp=1(2+1)m=|ℓm,ijp|2,\begin{eqnarray} \hat C_{\ell,ijp}={1\over{(2\ell+1)}}\sum_{m=-\ell}^{\ell}\vert \hat a_{\ell m,ijp}\vert^2, \end{eqnarray}(7)where the harmonic coefficients âℓm,ijp for a pixelated map of Npix pixels are given by ℓm,ijp=4πNpixk=1NpixδxijkpYℓm(nk).\begin{eqnarray} \hat a_{\ell m,ijp}={4\pi\over N_{\text{pix}}}\sum_{k=1}^{N_{\text{pix}}}\delta x_{ijkp}Y^{\ast}_{\ell m}(\boldsymbol{n}_k). \end{eqnarray}(8)We compute the mean power spectrum Ĉℓ,ip for the parameters { H0,q0,j0 } by averaging the power spectra Ĉℓ,ijp of the maps δxijkp = { δH0q0j0 } jkp,cosmic_var over the j realizations of the “Cosmic variance” local estimation ℓ,ip=ℓ,ijpj,\begin{eqnarray} \hat C_{\ell,ip}=\big<\hat C_{\ell,ijp}\big>_j, \end{eqnarray}(9)with variance Var[ℓ,ip]=Var[ℓ,ip]sample+Var[ℓ,ip]estimator,\begin{eqnarray} \text{Var}[\hat C_{\ell,ip}]=\text{Var}[\hat C_{\ell,ip}]^{\text{sample}}+\text{Var}[\hat C_{\ell,ip}]^{\text{estimator}}, \end{eqnarray}(10)where Var[ℓ,ip]estimator=22+12ℓ,ip·\begin{eqnarray} \text{Var}[\hat C_{\ell,ip}]^{\text{estimator}}={2\over {2\ell+1}}\hat C_{\ell,ip}^2\cdot \end{eqnarray}(11)We then compute the mean power spectrum Ĉℓ,i over the p pixelations ℓ,i=ℓ,ipp,\begin{eqnarray} \hat C_{\ell,i}=\big<\hat C_{\ell,ip}\big>_p, \end{eqnarray}(12)with variance given by Var[ℓ,i]=[p=1n1Var[ℓ,ip]]-1·\begin{eqnarray} \text{Var}[\hat C_{\ell,i}]=\left[\sum_{p=1}^n {1\over \text{Var}[\hat C_{\ell,ip}]}\right]^{-1}\cdot \end{eqnarray}(13)In order to compute the unbiased power spectrum we devise two methods, as detailed below.

3.1. Difference of power spectra

In the first method, we define the unbiased power spectrum as the unbiased mean power spectrum. We compute the unbiased mean power spectrum by computing the power spectrum Cℓ,ipbias\hbox{$C_{\ell,ip}^{\text{bias}}$} of the mean map δbiasikp=δxijkpbiasj={δ00}kp,shuffle_sn0\hbox{$\delta \bar x_{ikp}^{\text{bias}}=\big<\delta x_{ijkp}^{\text{bias}}\big>_j=\{\delta \bar H_0,\delta \bar q_0, \delta \bar j_0\}_{kp,\text{shuffle\_sn}}$} averaged over the j realizations of the “Shuffle SNe” local estimation and subtracting it from the mean power spectrum Ĉℓ,ipunbiasℓ,ip=ℓ,ipCℓ,ipbias,\begin{eqnarray} \hat C_{\ell,ip}^{\text{unbias}}=\hat C_{\ell,ip}-C_{\ell,ip}^{\text{bias}}, \end{eqnarray}(14)with total variance Var[unbiasℓ,ip]=Var[ℓ,ip]+Var[Cℓ,ipbias]estimator.\begin{eqnarray} \text{Var}[\hat C_{\ell,ip}^{\text{unbias}}] =\text{Var}[\hat C_{\ell,ip}]+\text{Var}[C_{\ell,ip}^{\text{bias}}]^{\text{estimator}}. \end{eqnarray}(15)This was the method suggested in Carvalho & Marques (2016). We then compute the mean unbiased power spectrum unbiasℓ,i\hbox{$\hat C_{\ell,i}^{\text{unbias}}$} by averaging Ĉℓ,ip over the different p pixelations unbiasℓ,i=unbiasℓ,ipp,\begin{eqnarray} \hat C_{\ell,i}^{\text{unbias}}=\big<\hat C_{\ell,ip}^{\text{unbias}}\big>_p, \end{eqnarray}(16)with total variance Var[unbiasℓ,i]=[p=1n1Var[unbiasℓ,ip]]-1·\begin{eqnarray} \text{Var}[\hat C_{\ell,i}^{\text{unbias}}]=\left[\sum_{p=1}^n {1\over \text{Var}[\hat C_{\ell,ip}^{\text{unbias}}]}\right]^{-1}\cdot \end{eqnarray}(17)In the right panels of Fig. 4, we plot the resulting unbiased power spectra, for the different pixelations and for the average over the p pixelations. The variability in Ĉℓ,ip and Cℓ,ipbias\hbox{$C_{\ell,ip}^{\text{bias}}$} from the different pixelations is indicated as colour-shaded regions bordered by dashed and dotted lines respectively. For all parameters, the power spectrum has a maximum at the quadrupole ( = 2). However, given the size of the error, the results are also compatible with a flat spectrum.

For comparison, in the left panels of Fig. 4, using the same line types, we plot the power spectra for the parameters { H0,q0 = ΩM/ 2−ΩΛ,j0 = ΩM + ΩΛ } from the { H0MΛ } estimation in Carvalho & Marques (2016). In this estimation, the unbiased power spectra also follow the same behaviour as the power spectra before the noise bias removal, with the exception of q0 whose subtle maximum at = 2 is erased with the noise bias removal and the power spectrum decreases always with the multipole. Conversely, for H0 the power spectrum increases always with the multipole, and for j0 the power spectrum has a maximum at = 2. Hence, between the { H0MΛ } and the { H0,q0,j0 } estimation, we observe the creation of a peak at = 2 for H0 and q0, and the smoothing of the peak for j0. We also observe that the power spectra in the { H0,q0,j0 } estimation are up to 5 × 102 times the power spectra in the { H0MΛ } estimation, hence yielding an increase of the amplitude.

The power spectrum that accounts for the noise bias Cℓ,ibias\hbox{$C_{\ell,i}^{\text{bias}}$} is up to two orders of magnitude smaller that the mean power spectra Ĉℓ,i, the resulting unbiased mean power spectra unbiasℓ,i\hbox{$\hat C_{\ell,i}^{\text{unbias}}$} following the same behaviour as Ĉℓ,i. Since Cℓ,ibiasℓ,i\hbox{$C_{\ell,i}^{\text{bias}}\ll \hat C_{\ell,i}$} for the { H0,q0,j0 } estimation, this method might not remove entirely the noise bias contribution to the power spectrum. For this reason, we conceived another method.

thumbnail Fig. 4

Power spectrum of the parameters estimated from the JLA type Ia SN sample (method A). The dashed lines are the mean power spectra from the “Cosmic variance” estimation, the dotted lines are the power spectra of the average “Shuffle SNe” estimation and the solid lines are the unbiased power spectra defined as the difference between the former two. Left panels: power spectra from the { H0MΛ } estimation. Right panels: power spectra from the { H0,q0,j0 } estimation.

thumbnail Fig. 5

Power spectrum of the parameters estimated from the JLA type Ia SN sample (method B). The dashed lines are the power spectra of the average “Cosmic variance” estimation, the dotted lines are the power spectra of the average “Shuffle SNe” estimation and the solid lines are the unbiased power spectra defined as the power spectra of the average unbiased maps. Left panels: power spectra from the { H0MΛ } estimation. Right panels: Power spectra from the { H0,q0,j0 } estimation.

3.2. Power spectrum of unbiased map

In the second method, we define the unbiased power spectrum as the power spectrum of the mean unbiased maps. We compute the unbiased power spectrum by computing the power spectrum Cℓ,ipunbias\hbox{$ C_{\ell,ip}^{\text{unbias}}$} of the unbiased mean maps δunbiasikp=(ikpbiasikp)/xfid,\hbox{$\delta \bar x_{ikp}^{\text{unbias}}=(\bar x_{ikp}-\bar x_{ikp}^{\text{bias}})/x^{\text{fid}},$} with total variance Var[Cℓ,ipunbias]=Var[Cℓ,ipunbias]estimator.\begin{eqnarray} \text{Var}[C_{\ell,ip}^{\text{unbias}}]=\text{Var}[C_{\ell,ip}^{\text{unbias}}]^{\text{estimator}} . \end{eqnarray}(18)We then average over the different p pixelations, defining Cℓ,iunbias=Cℓ,ipunbiasp,\begin{eqnarray} C_{\ell,i}^{\text{unbias}}=\big<C_{\ell,ip}^{\text{unbias}}\big>_p, \end{eqnarray}(19)with variance Var[Cℓ,iunbias]=[p=1n1Var[Cℓ,ipunbias]]-1.\begin{eqnarray} \text{Var}[C_{\ell,i}^{\text{unbias}}]=\left[\sum_{p=1}^n {1\over \text{Var}[C_{\ell,ip}^{\text{unbias}}]}\right]^{-1}. \end{eqnarray}(20)In the right panels of Fig. 5, we plot the resulting unbiased power spectra, for the different pixelations and the average over the p pixelations. The variability in Cℓ,ip (the power spectrum of the mean map δikp=δxijkpj={δ00}kp,cosmic_var0\hbox{$\delta \bar x_{ikp}=\big<\delta x_{ijkp}\big>_j=\{\delta \bar H_0,\delta \bar q_0, \delta \bar j_0\}_{kp,\text{cosmic\_var}}$} averaged over the j realizations of the “Cosmic variance” local estimation) and Cℓ,ipbias\hbox{$C_{\ell,ip}^{\text{bias}}$} (the power spectrum of the mean map δbiasikp=δxijkpbiasj={δ00}kp,shuffle_sn0\hbox{$\delta \bar x_{ikp}^{\text{bias}}=\big<\delta x_{ijkp}^{\text{bias}}\big>_j=\{\delta \bar H_0,\delta \bar q_0, \delta \bar j_0\}_{kp,\text{shuffle\_sn}}$} averaged over the j realizations of the “Shuffle SNe” local estimation) from the different pixelations is indicated as colour–shaded regions bordered by dashed and dotted lines respectively and included for reference. For all parameters, the power spectrum has a maximum at = 2. Although this method has smaller errors than the previous method, the results are still compatible with a flat spectrum.

For comparison, in the left panels of Fig. 5, using the same line types, we plot the power spectra for the parameters { H0,q0 = ΩM/ 2−ΩΛ,j0 = ΩM + ΩΛ } from the { H0MΛ } estimation in Carvalho & Marques (2016). In this estimation, both H0 and j0 have a maximum at = 2, whereas q0 decreases always with the multipole. Hence, between the { H0MΛ } and the { H0,q0,j0 } estimation, we observe the creation of a peak at = 2 for q0 and the smoothing of the peak for H0 and j0. We also observe that the power spectra in the { H0,q0,j0 } estimation are up to 1 × 103 times the power spectra in the { H0MΛ } estimation, hence yielding an increase of the amplitude.

The two methods return qualitatively equivalent results; quantitatively, however, the second method might be more efficient at removing the noise bias.

4. Average values

Table 2

Values for the parameters estimated from the JLA type Ia SN sample.

In order to compute the average values of the parameters from the corresponding mean maps \hbox{$\bar x_{ik},$} we average the map of each parameter over the pixel subsamples. When averaging over homogeneous pixels, angular fluctuations in the expansion factor (and consequently in H0) induce a backreaction term in the average deceleration parameter in the form of an extra positive acceleration (Räsänen 2006). By the same reasoning, angular fluctuations in the expansion factor (and consequently in H0 and q0) will also induce a backreaction term in the average jerk parameter. In Carvalho & Marques (2016) we derived the analytical extra positive acceleration for a toy model of an arbitrary number of disjoint, homogeneous regions and computed the overall deceleration parameter assuming a) no backreaction and b) backreaction for the measured angular distribution of H0. Here, for the same toy model, we derive the corresponding extra terms for the time variation in the acceleration and compute the overall jerk parameter assuming a) no backreaction and b) backreaction for the measured angular distribution of H0 and q0.

In the absence of backreaction, the averaging consists in taking the mean weighted by the variance’s inverse \hbox{$w_k=1/\text{Var}[\bar x_{ik}]$} of parameter xi in each pixel k,i=ikk=kNpixelwkikkNpixelwk·\begin{eqnarray} \bar x_i=\big< \bar x_{ik}\big>_{k}={\sum_{k}^{N_{\text{pixel}}}w_k~\bar x_{ik} \over {\sum_{k}^{N_{\text{pixel}}}w_k}}\cdot \end{eqnarray}(21)For the estimated parameters, we find {H0,q0,j0}={70.99,0.497,0.538}±{0.54,0.126,0.214}.\hbox{$\{\overline H_{0},\overline q_0,\overline j_0\} =\{70.99,-0.497,0.538\}\pm \{0.54,0.126,0.214\}.$} After the noise bias removal, we find {H0,q0,j0}unbias={71.32,0.570,0.540}±{0.60,0.162,0.333}.\hbox{$\{\overline H_{0},\overline q_0,\overline j_0\}_{\text{unbias}} =\{71.32,-0.570,0.540\}\pm \{0.60,0.162,0.333\}.$} These values are consistent with the fiducial values. After the noise bias removal, the pixel average values become closer to the corresponding values estimated using the complete sample.

In the presence of backreaction, the averaging consists in taking the mean weighted by the three–volume Vk of each pixel k,i=ikVk=kNpixelVkikkNpixelVk·\begin{eqnarray} \bar x_i=\big< \bar x_{ik}\big>_{V_k}={\sum_{k}^{N_{\text{pixel}}}V_k~\bar x_{ik} \over {\sum_{k}^{N_{\text{pixel}}}V_k}}\cdot \end{eqnarray}(22)Identifying a volume as a pixel and defining vk = Vk/ ∑ kVk, then for Npixel disjoint regions, the average of the Hubble parameter is given by H0Vk=kNpixelvkH0,k.\begin{eqnarray} \left \langle H_0\right\rangle _{V_k}=\sum_{k}^{N_{\text{pixel}}}v_{k} H_{0,k}. \label{eqn:a_dot_Vk} \end{eqnarray}(23)

4.1. Average value of the deceleration parameter

In order to derive the volume average of the acceleration, we take the time derivative of Eq. (23) and find that ¨aaVk=kNpixelvk¨aa+2kNpixell>kNpixelvkvl(H0,kH0,l)2.\begin{eqnarray} \left \langle {\ddot a \over a}\right\rangle _{V_k}= \sum_{k}^{N_{\text{pixel}}}v_{k} {\ddot a \over a} +2\sum_{k}^{N_{\text{pixel}}}\sum_{l>k}^{N_{\text{pixel}}}v_{k}v_{l}\left(H_{0,k}-H_{0,l}\right)^2. \label{eqn:a_dotdot_Vk} \end{eqnarray}(24)Equation (24) decomposes into a linear term in the pixel average of (¨a/a)\hbox{$(\ddot a/a)$} and a quadratic term in differences of H0 between pairs of pixels. The quadratic (backreaction) term generates an acceleration due to the slower regions becoming less represented in the average. In the absence of the quadratic term, the volume average reduces to the pixel average above. Then the volume average of q0 becomes q0Vk=kNpixelvkq0,k\begin{eqnarray} \left \langle q_{0}\right\rangle _{V_k}&=&\sum_{k}^{N_{\text{pixel}}}v_{k}q_{0,k}\nonumber\\ &&-{2\over \left(\sum_{k}^{N_{\text{pixel}}} v_{k}H_{0,k}\right)^2} \sum_{k}^{N_{\text{pixel}}}\sum_{l>k}^{N_{\text{pixel}}}v_{k}v_{l}\left(H_{0,k}-H_{0,l}\right)^2. \label{eqn:q0_Vk} \end{eqnarray}(25)Using the fluctuations in H0 (measured in the “Cosmic variance” local estimation) we find \hbox{$\bar q_{0}=-0.498\pm0.126,$} and after the noise bias removal we find \hbox{$\bar q_{0,\text{unbias}}=-0.570\pm0.162$} (see Table 2).

The quadratic term is of order 7 × 10-4 times the linear term; the corresponding ratio measured in the { H0MΛ } estimation was of order 10-3 (Carvalho & Marques 2016). The difference due to the backreaction is below the standard deviation, hence unobservable. It follows that, for the angular fluctuations in H0 measured with this SN sample, the contribution of the quadratic term in Eq. (25) is insignificant, which renders the volume averaging equivalent to the pixel averaging. These results confirm that, in the context of this toy model of an inhomogeneous space–time and for a kinetic parametrisation, backreaction is not a viable dynamical mechanism to emulate cosmic acceleration.

4.2. Average value of the jerk parameter

In order to derive the volume average of the time variation in the acceleration, we take the time derivative of Eq. (24) and find that ...aaVk=kNpixelvk...akak+2kNpixell>kNpixelvkvl(H0,kH0,l)2mNpixelvmH0,m+2kNpixell>kNpixelvkvl[(q0,kH0,k2)(q0,lH0,l2)]\begin{eqnarray} \left \langle {\dddot a \over a}\right\rangle _{V_k}&=& \sum_{k}^{N_{\text{pixel}}}v_{k} {\dddot a_{k} \over a_{k}}\nonumber\\ &&+2\sum_{k}^{N_{\text{pixel}}}\sum_{l>k}^{N_{\text{pixel}}}v_{k}v_{l}\left(H_{0,k}-H_{0,l}\right)^2 \sum_{m}^{N_{\text{pixel}}}v_{m} H_{0,m}\nonumber\\ &&+2\sum_{k}^{N_{\text{pixel}}}\sum_{l>k}^{N_{\text{pixel}}}v_{k}v_{l} \left[\left(-q_{0,k}H_{0,k}^2\right)-\left(-q_{0,l}H_{0,l}^2\right)\right]\nonumber\\ &&\times\left(H_{0,k}-H_{0,l}\right). \label{eqn:a_dotdotdot_Vk} \end{eqnarray}(26)Equation (26) decomposes into a linear term in the pixel average of (...a/a),\hbox{$(\dddot a/a),$} a quadratic term in differences of H0 between pairs of pixels and a linear term in differences of (q0H02)\hbox{$(-q_0H_{0}^2)$} between pairs of pixels. The quadratic term in differences of H0 generates a contribution that is always positive similar to the quadratic term in Eq. (24). Conversely, the linear term in differences of (q0H02)\hbox{$(-q_0H_{0}^2)$} generates a contribution that can have either sign; in particular, it will be positive in the pair of pixels where (− q0) and H0 vary in the same direction (i.e. both quantities increase or decrease between pixels) and it will be negative in the pair of pixels where (− q0) and H0 vary in opposite directions (i.e. one quantity increases while the other decreases). Since Eq. (26) measures the angular average of the time variation in the acceleration, the difference (backreaction) terms generate an extra jerk that is due to slower regions and/or more slowly varying regions becoming less represented in the average. Similarly to Eq. (24), in the absence of the difference terms, the volume average reduces to the pixel average above. Then the volume average of j0 becomes j0Vk=kNpixelvkj0,k+2(kNpixelvkH0,k)2kNpixell>kNpixelvkvl(H0,kH0,l)22(kNpixelvkH0,k)3\begin{eqnarray} \left \langle j_{0}\right\rangle _{V_k}&=& \sum_{k}^{N_{\text{pixel}}}v_{k} {j_{0,k}}\nonumber\\ &&+{2\over \left(\sum_{k}^{N_{\text{pixel}}} v_{k}H_{0,k}\right)^2} \sum_{k}^{N_{\text{pixel}}}\sum_{l>k}^{N_{\text{pixel}}}v_{k}v_{l}\left(H_{0,k}-H_{0,l}\right)^2\nonumber\\ &&-{2\over \left(\sum_{k}^{N_{\text{pixel}}} v_{k}H_{0,k}\right)^3}\nonumber\\ &&\times\sum_{k}^{N_{\text{pixel}}}\sum_{l>k}^{N_{\text{pixel}}}v_{k}v_{l} \left(q_{0,k}H_{0,k}^2-q_{0,l}H_{0,l}^2\right)\left(H_{0,k}-H_{0,l}\right). \label{eqn:j0_Vk} \end{eqnarray}(27)Using the fluctuations in { H0,q0 } (measured in the “Cosmic variance” local parameter estimation), we find \hbox{$\bar j_{0}=0.546\pm 0.218$} and after the noise bias removal, we find \hbox{$\bar j_{0,\text{unbias}}=0.548\pm 0.345$} (see Table 2).

The quadratic term in differences of H0 is of order 1 × 10-3 times the linear term, whereas the linear term in differences of (q0H02)\hbox{$(-q_0H_{0}^2)$} is of order 1 × 10-2 times the linear term; the corresponding ratios in the { H0MΛ } estimation in Carvalho & Marques (2016) were of order 6 × 10-4 and 2 × 10-3. Since j0 is more poorly constrained than q0 or H0 in the global estimation, the total difference due to the backreaction is still below the standard deviation, hence unobservable. It follows that, for the angular fluctuations in { H0,q0 } measured with this SN sample, the contribution of the backreaction terms in Eq. (27) is insignificant, which renders the pixel averaging equivalent to the pixel averaging. These results imply that an inhomogeneous j0, such that at different pixels the acceleration changed at different times, cannot be distinguished from a globally homogeneous j0, such that the acceleration changed everywhere at the same time.

For comparison, we present the pixel average of { q0 = ΩM/ 2−ΩΛ,j0 = ΩM + ΩΛ } from the { H0MΛ } estimation in Carvalho & Marques (2016), both in the absence and in the presence of backreaction, which we include in Table 2.

5. Conclusions

In this paper we used SN data to fit a kinematic parametrisation of the luminosity distance expressed in terms of time derivatives of the scale factor. This parametrisation records the cosmic expansion without regard to its cause, thus being independent of the cosmological equation of state and consequently adequate to test interpretations of the cosmic acceleration alternative to the cosmological constant. We followed the parameter estimation, first presented in Carvalho & Marques (2016), to fit the parameters { H0,q0,j0 } by performing both a global and a local parameter estimation. From the global parameter estimation, using the complete SN sample, we obtained the fiducial values adopted in this manuscript. From the local parameter estimation, dividing the SNe into subsamples over a pixelated map, we obtained maps of the estimated parameters with the same pixelation as the SN subsamples. We then proceeded to the analysis of this paper’s results as well as to a comparative analysis with the results from Carvalho & Marques (2016) estimated by fitting { H0MΛ } instead.

The measurements of { H0,q0,j0 } are consistent with the measurements of { H0,q0 = ΩM/ 2−ΩΛ,j0 = ΩM + ΩΛ } obtained in Carvalho & Marques (2016). However, whereas the error of q0 is of the same order in both parametrisations (about a 5σ measurement), the error of j0 is significantly larger in the kinematic parametrisation (from a 7σ to a 1σ measurement). This is a consequence of the kinematic parametrisation in part minimising the physical assumptions that enter in the model and in part truncating the Taylor expansion of the luminosity distance.

We measured fluctuations about the average values of order 0.1–5% for H0, 1–150% for q0 and 1–124% for j0. Comparing with the fluctuations measured in Carvalho & Marques (2016), we observe an increase by a factor of two orders of magnitude. We also computed the power spectrum of the corresponding maps of the parameters up to = 3, as determined by the pixel size, finding that all power spectra have a maximum at = 2, regardless of the method used to subtract the noise bias. This observation can be partially ascribed to the absence of objects towards the galactic plane.

Finally, for an analytical toy model of an inhomogeneous ensemble of homogenous pixels, we measured the backreaction term in q0 due to the fluctuations of H0 to be of order 5 × 10-4 the corresponding pixel average in the absence of backreaction, hence of smaller order than that measured in Carvalho & Marques (2016). We also derived the backreaction term in j0 due to the fluctuations of { H0,q0 } and measured it to be of order 1 × 10-2 the corresponding pixel average in the absence of backreaction, hence about 50 times that measured in Carvalho & Marques (2016). However, both backreaction effects are below the corresponding standard deviation and hence rendered unobservable. It follows that backreaction generates insignificant extra acceleration to emulate cosmic acceleration, and insignificant extra jerk to emulate different late–time cosmic histories.

An interesting idea would be to extend this type of analysis to further inhomogeneity studies by combining SN data with other cosmic tracers, such as high–redshift galaxies (Terlevich et al. 2015) and galaxy clusters (Bengaly et al. 2015b; Bolejko et al. 2016).


Acknowledgments

The authors thank A. Afonso and N. Nunes for useful discussions. CSC is funded by Fundação para a Ciência e a Tecnologia (FCT), Grant no. SFRH/BPD/65993/2009.

References

  1. Amendola, L., Appleby, S., Bacon, D., et al. 2013, Liv. Rev. Relat., 16, 6 [Google Scholar]
  2. Aviles, A., Gruber, C., Luongo, O., & Quevedo, H. 2012, Phys. Rev. D, 86, 123516 [NASA ADS] [CrossRef] [Google Scholar]
  3. Bengaly, C. A. P., Bernui, A., & Alcaniz, J. S. 2015a, ApJ, 808, 39 [NASA ADS] [CrossRef] [Google Scholar]
  4. Bengaly, C. A. P., Bernui, A., Alcaniz, J. S., & Ferreira, I. S. 2015b, ArXiv eprints [arXiv:1511.09414] [Google Scholar]
  5. Betoule, M., Kessler, R., Guy, J., et al. 2014, A&A, 568, 22 [Google Scholar]
  6. Blomqvist, M., Enander, J., & Mortsell, E. 2010, Constraining dark energy fluctuations with supernova correlations, JCAP, 10, 018 [NASA ADS] [CrossRef] [Google Scholar]
  7. Bolejko, K., Nazer, M. A., & Wiltshire, D. L. 2016, JCAP, 06, 35 [Google Scholar]
  8. Caldwell, R. R., & Kamionkowski, M. 2004, JCAP, 0409, 009 [NASA ADS] [CrossRef] [Google Scholar]
  9. Caldwell, R. R., & Stebbins, A. 2008, Phys. Rev. Lett., 100, 191302 [NASA ADS] [CrossRef] [MathSciNet] [PubMed] [Google Scholar]
  10. Carvalho, C. S., & Marques, K. 2016, A&A, 592, A102 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Cattoën, C., & Visser, M. 2008, Phys. Rev. D, 78, 063501 [NASA ADS] [CrossRef] [Google Scholar]
  12. Górski, K. M., Hivon, E., Banday, A. J., et al. 2005, ApJ, 622, 759 [NASA ADS] [CrossRef] [Google Scholar]
  13. Heneka, C., Marra, V., & Amendola, L. 2014, MNRAS, 439, 1855 [NASA ADS] [CrossRef] [Google Scholar]
  14. Javanmardi, B., Porciani, C., Kroupa, P., & Pflamm-Altenburg, J. 2015, ApJ, 810, 47 [NASA ADS] [CrossRef] [Google Scholar]
  15. Jiménez, J. B., Salzano, V., & Lazkoz, R. 2015, Phys. Lett. B, 741, 168 [NASA ADS] [CrossRef] [Google Scholar]
  16. Kalus, B., Schwarz, D. J., Seikel, M., & Wiegand, A. 2013, A&A, 553, A56 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Lan, M. X., Li, M., Li, X. D., & Wang, S. 2010, Phys. Rev. D, 82, 023516 [NASA ADS] [CrossRef] [Google Scholar]
  18. Liu, S., & Zhang, T. J. 2014, Phys. Lett. B, 733, 69 [NASA ADS] [CrossRef] [Google Scholar]
  19. Mariano, A., & Perivolaropoulos, L. 2012, Phys. Rev. D, 86, 083517 [NASA ADS] [CrossRef] [Google Scholar]
  20. Marra, V., & Notari, A. 2011, Class. Quant. Grav., 28, 164004 [NASA ADS] [CrossRef] [Google Scholar]
  21. Marra, V., & Pääkkönen, M. 2010, JCAP, 1012, 021 [NASA ADS] [CrossRef] [Google Scholar]
  22. Migkas, K., & Plionis, M. 2016, Rev. Mex. Astron. Astrofis., 52, 133 [Google Scholar]
  23. Moss, A., Zibin, J. P., & Scott, D. 2011, Phys. Rev. D, 83, 103515 [NASA ADS] [CrossRef] [Google Scholar]
  24. Neben, A. R., & Turner, M. S. 2013, ApJ, 769, 133 [NASA ADS] [CrossRef] [Google Scholar]
  25. Planck Collaboration XIII. 2016, A&A, in press, DOI: 10.1051/0004-6361/201525830 [Google Scholar]
  26. Räsänen, S. 2006, JCAP, 11, 003 [Google Scholar]
  27. Riess, A. G., Strolger, L. G., Tonry, J., et al. 2004, ApJ, 607, 665 [NASA ADS] [CrossRef] [Google Scholar]
  28. Riess, A. G., Macri, L., Casterano, S., et al. 2011, ApJ, 730, 119 [NASA ADS] [CrossRef] [Google Scholar]
  29. Schwarz, D. J., & Weinhorst, B. 2007, A&A, 474, 717 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. Terlevich, R., Terlevich, E., Melnick, J., et al. 2015, MNRAS, 451, 3001 [NASA ADS] [CrossRef] [Google Scholar]
  31. Tomita, K., & Inoue, K. T. 2009, Phys. Rev. D, 79, 103505 [NASA ADS] [CrossRef] [Google Scholar]
  32. Turner, M. S., & Riess, A. G. 2002, AJ, 569, 18 [NASA ADS] [CrossRef] [Google Scholar]
  33. Visser, M. 2004, Class. Quant. Grav., 21, 2603 [NASA ADS] [CrossRef] [Google Scholar]
  34. Wiegand, A., & Schwarz, D. J. 2012, A&A, 538, A147 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  35. Zehavi, I., Riess, A. G., Kirshner, R. P., & Dekel, A. 1998, ApJ, 503, 483 [NASA ADS] [CrossRef] [Google Scholar]
  36. Zibin, P., Moss, A., & Scott, D. 2008, Phys. Rev. Lett., 101, 251303 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]

All Tables

Table 1

Values for the parameters estimated from the JLA type Ia SN sample.

Table 2

Values for the parameters estimated from the JLA type Ia SN sample.

All Figures

thumbnail Fig. 1

Relative difference between expressions of the luminosity distance. The solid line shows the relative difference between the luminosity distance dL(z) in Eq. (3) and the theoretical luminosity distance in the ΛCDM model in Eq. (6) for ΩM = 0.308 and H0 = 67.8 km s-1 Mpc-1 (Planck Collaboration XIII 2016). The dashed line shows the relative difference between the luminosity distance dL(y), with yz/ (1 + z), in Cattoën & Visser (2008) and the theoretical luminosity distance in the ΛCDM model. The inset plot shows the corresponding differences in μ, given by Eq. (4), which are indicated by the corresponding line styles.

In the text
thumbnail Fig. 2

Left panel: distance modulus as a function of redshift. The dots are coloured according to the original SN survey, as described in Betoule et al. (2014), and indicate the distance modulus μdata computed from each SN in the sample. Blue: several low-z. Green: SDSS-II. Yellow: SNLS. Red: HST. The black line indicates the distance modulus computed theoretically for the values { H0,q0,j0 } estimated using the complete sample. The black dashed and black dotted lines indicate the distance modulus for two extreme cases. Right panel: Markov chain of one realization. The dots indicate the entries in the Markov chain. The position in the graph indicates the values of (q0,j0) of each entry according to the axes, whereas the colour indicates the value of H0 as detailed in the side colour bar. The dashed line represents j0 = q0(1 + 2q0), which separates between flat and curved space–time.

In the text
thumbnail Fig. 3

Fluctuations of the estimated parameters as a function of the number of SNe per pixel, before (left panel sets) and after (right panel sets) noise bias removal. Each panel set consists of five panels. At each pixel we plot: a) in the first panel, the number of SNe; b) in the second panel, the value of the parameters { q0,j0 } respectively as the solid red and solid blue lines, with the dashed red and dashed blue lines indicating the corresponding fiducial values; c) in the third panel, the standard deviation { σq0,σj0 } respectively as solid red and solid blue lines; d) in the fourth plot, the value of H0 as the solid black line, with the dashed black line marking the fiducial value; e) in the fifth panel, the standard deviation σH0 as the solid black line. Top panel sets: results from the { H0,q0,j0 } estimation for one pixelation. Centre panel sets: results from the { H0,q0,j0 } estimation for the average over the different pixelations. Bottom panel sets: results from the { H0MΛ } estimation for Ωκ = 1−ΩM−ΩΛ.

In the text
thumbnail Fig. 4

Power spectrum of the parameters estimated from the JLA type Ia SN sample (method A). The dashed lines are the mean power spectra from the “Cosmic variance” estimation, the dotted lines are the power spectra of the average “Shuffle SNe” estimation and the solid lines are the unbiased power spectra defined as the difference between the former two. Left panels: power spectra from the { H0MΛ } estimation. Right panels: power spectra from the { H0,q0,j0 } estimation.

In the text
thumbnail Fig. 5

Power spectrum of the parameters estimated from the JLA type Ia SN sample (method B). The dashed lines are the power spectra of the average “Cosmic variance” estimation, the dotted lines are the power spectra of the average “Shuffle SNe” estimation and the solid lines are the unbiased power spectra defined as the power spectra of the average unbiased maps. Left panels: power spectra from the { H0MΛ } estimation. Right panels: Power spectra from the { H0,q0,j0 } estimation.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.