Issue 
A&A
Volume 598, February 2017



Article Number  A112  
Number of page(s)  14  
Section  Cosmology (including clusters of galaxies)  
DOI  https://doi.org/10.1051/00046361/201628909  
Published online  09 February 2017 
Cosmology with gammaray bursts
I. The Hubble diagram through the calibrated E_{p,i}–E_{iso} correlation
^{1} Institute for Theoretical Physics, University of Warsaw, Pasteura 5, 02093 Warsaw, Poland
^{2} Department of Astronomy, Williams College, Williamstown, MA 01267, USA
^{3} Dipartimento di Fisica, Università degli Studi di Napoli Federico II, Compl. Univ. Monte S. Angelo, 80126 Napoli, Italy
email: ester@na.infn.it
^{4} I.N.F.N., Sez. di Napoli, Compl. Univ. Monte S. Angelo, Edificio 6, via Cinthia, 80126 Napoli, Italy
^{5} Dipartimento di Fisica e Scienze della Terra, Università degli Studi di Ferrara, 44121 Ferrara, Italy
^{6} Department of Physics, University of Nice Sophia Antipolis, Parc Valrose, 06034 Nice Cedex 2, France
^{7} INAF–IASF, Sezione di Bologna, via Gobetti 101, 40129 Bologna, Italy
Received: 12 May 2016
Accepted: 7 September 2016
Context. Gammaray bursts (GRBs) are the most energetics explosions in the Universe. They are detectable up to very high redshifts. They may therefore be used to study the expansion rate of the Universe and to investigate the observational properties of dark energy, provided that empirical correlations between spectral and intensity properties are appropriately calibrated.
Aims. We used the type Ia supernova (SN) luminosity distances to calibrate the correlation between the peak photon energy, E_{p,i}, and the isotropic equivalent radiated energy, E_{iso} in GRBs. With this correlation, we tested the reliability of applying these phenomena to measure cosmological parameters and to obtain indications on the basic properties and evolution of dark energy.
Methods. Using 162 GRBs with measured redshifts and spectra as of the end of 2013, we applied a local regression technique to calibrate the E_{p,i}–E_{iso} correlation against the type Ia SN data to build a calibrated GRB Hubble diagram. We tested the possible redshift dependence of the correlation and its effect on the Hubble diagram. Finally, we used the GRB Hubble diagram to investigate the dark energy equation of state (EOS). To accomplish this, we focused on the socalled ChevalierPolarskiLinder (CPL) parametrization of the dark energy EOS and implemented the Markov chain Monte Carlo (MCMC) method to efficiently sample the space of cosmological parameters.
Results. Our analysis shows once more that the E_{p,i}–E_{iso} correlation has no significant redshift dependence. Therefore the highredshift GRBs can be used as a cosmological tool to determine the basic cosmological parameters and to test different models of dark energy in the redshift region (z ≥ 3), which is unexplored by the SNIa and baryonic acoustic oscillations data. Our updated calibrated Hubble diagram of GRBs provides some marginal indication (at 1σ level) of an evolving dark energy EOS. A significant enlargement of the GRB sample and improvements in the accuracy of the standardization procedure is needed to confirm or reject, in combination with forthcoming measurements of other cosmological probes, this intriguing and potentially very relevant indication.
Key words: cosmological parameters / dark energy / cosmology: theory
© ESO, 2017
1. Introduction
Starting at the end of the 1990s, observations of highredshift supernovae of type Ia (SNIa) revealed that the Universe is now expanding at an accelerated rate (see e.g. Perlmutter et al. 1998, 1999; Riess et al. 1998, 2007; Schmidt et al. 1998; Astier et al. 2006; Amanullah et al. 2010). This surprising result has been independently confirmed by statistical analyses of observations of smallscale temperature anisotropies of the cosmic microwave background (CMB) radiation Spergel et al. (2007), Planck Collaboration XVI (2014), Planck Collaboration XIII (2016). It is usually assumed that the observed accelerated expansion is caused by the socalled dark energy, a cosmic medium with unusual properties. The pressure of dark energy p_{de} is negative and is related to the positive energy density of dark energy ϵ_{de} by p_{de} = wϵ_{de}, where the proportionality coefficient, that is, the equation of state (EOS), w, is negative (w< −1/3). According to current estimates, about 70% of the matter energy in the Universe is in the form of dark energy, so that today dark energy is the dominant component in the Universe. The nature of dark energy is, however, not known. The models of dark energy proposed so far can be divided into at least three groups: a) a nonzero cosmological constant (see for instance Carroll 2001), in this case w = −1, b) a potential energy of some not yet discovered scalar field (see for instance Sahni et al. 2003; Alam et al. 2003), or c) effects connected with the inhomogeneous distribution of matter and averaging procedures (see for instance Clarkson & Maartens 2010). In the last two cases, in general, w ≠ −1 and is not constant, but depends on the redshift z. To test whether and how w changes with redshift, it is necessary to use more distant objects. It is commonly argued that since the dark energy density term becomes subdominant (with respect to the dark matter) at z ≳ 0.5Riess et al. (2004), it is not important to study its EOS at earlier epochs. However, this argument is unreliable: even within the simplest model, the dark energy contributes nearly ≃10% of the overall cosmic energy density at z ≃ 2 and strongly alters the deceleration parameter with the cosmological constant. Moreover, for several observables the sensitivity to the dark energy equation of state increases at high redshifts. In Fig. 1 we explore this aspect following a simplified approach, considering the modulus of distance μ(z) as observable: we fixed a flat Λ cold dark matter (ΛCDM) fiducial cosmological model, constructing the corresponding μ_{fid}(z,θ), and plot the percentage error on the distance modulus with respect to different corresponding functions evaluated in the framework of a flat CPL quintessence model Chevallier & Polarski (2001), Linder (2003). We selected Ω_{m} = 0.3 and h = 0.7 and varied the dark EOS parameters w_{0} and w_{1}. It is worth noting that a higher sensitivity is reached only for z ≳ 2. Therefore, investigating the cosmic expansion of the Universe also beyond these redshifts remains a fundamental probe of dark energy. In this scenario, new possibilities opened up when gammaray bursts (GRBs) were discovered at higher redshifts. The current record is at z = 9.4 (Tanvir et al. 2009; Salvaterra et al. 2009; Cucchiara 2011). It is worth noting that the photometric redshift on GRB 090429B is quite high, especially on the low side; GRB 090423, for which a spectroscopic redshift is available, is much better determined. GRBs span a redshift range better suited for probing dark energy than the SNIa range, as shown in Fig. 2, where we compare the distribution in redshift of our sample of 162 long GRBs/XRFs with the Union 2.1 SNIa dataset.
Fig. 1
Dependence of the distance modulus on EOS of dark energy. Upper panel: distance modulus μ(z) for different values of the EOS parameters. The black line represents the standard flat ΛCDM model with Ω_{m} = 0.3, h = 0.7, w_{0} = −1, w_{1} = 0. The other curves correspond to different values of the CPL parameters w_{0} and w_{1}. Bottom panel: z dependence of the percentage error in the distance modulus between the fiducial ΛCDM model and different flat CPL models, as in the left panel. 
Fig. 2
Redshift distribution of our sample of 162 long GRBs/XRFs and the Union 2.1 SNIa dataset. 
Gammaray bursts are enigmatic objects, however. First of all, the mechanism that is responsible for releasing the high amounts of energy that a typical GRB emits is not yet completely known, and only some aspects of the progenitor models are well established, in particular that long GRBs are produced by corecollapse events, see for instance Meszaros (2006) and Vedrenne & Atteia (2009). Despite these difficulties, GRBs are promising objects for studying the expansion rate of the Universe at high redshifts. One of the most important aspects of the observational property of long GRBs is that they show several correlations between spectral and intensity properties (luminosity, radiated energy). Here we consider the correlation between the observed photon energy of the peak spectral flux, E_{p,i}, which corresponds to the peak in the νF_{ν} spectra, and the isotropic equivalent radiated energy E_{iso} (e.g., Amati et al. 2002; Amati 2006), (1)where a and b are constants, and E_{p,i} is the spectral peak energy in the GRB cosmological restframe, which can be derived from the observer frame quantity, E_{p}, by E_{p,i} = E_{p}(1 + z). This correlation not only provides constraints for the model of the prompt emission, but also naturally suggests that GRBs can be used as distance indicators. The isotropic equivalent energy E_{iso} can be calculated from the bolometric fluence S_{bolo} as (2)where d_{L} is the luminosity distance and cp denotes the set of parameters that specify the background cosmological model. It is clear that to be able to use GRBs as distance indicators, it is necessary to consistently calibrate this correlation. Unfortunately, owing to the lack of GRBs at very low redshifts, the calibration of GRBs is more difficult than that of SNIa, and several calibration procedures have been proposed so far (see for instance Cardone et al. 2008; Demianski et al. 2011, 2012; Gao et al. 2012; Postnikov et al. 2014; Dai et al. 2004; Liu & Wei 2015; Lin et al. 2015). We here apply a local regression technique to determine the correlation parameters a and b, using the recently updated SNIa sample and to construct a new calibrated GRB Hubble diagram that can be used for cosmological investigations. We then use this calibrated GRB Hubble diagram to investigate the cosmological parameters through the Markov chain Monte Carlo technique (MCMC), which simultaneously computes the full posterior probability density functions of all the parameters. The structure of the paper is as follows. In Sect. 2 we describe the methods used to fit the E_{p,i}–E_{iso} correlation and construct the calibrated GRB Hubble diagram. In Sect. 3 we present our cosmological constraints and investigate the possible redshift dependence and Malmquistlike bias effects. In Sect. 4, as an additional selfconsistency check, we apply the Bayesian method for the noncalibrated E_{p,i}–E_{iso} correlation and simultaneously extract the correlation coefficients and the cosmological parameters of the model. Section 5 is devoted to the discussion of our main results and conclusions.
2. Standardizing GRBs and constructing the Hubble diagram
The GRBs νF_{ν} spectra are well modeled by a phenomenological smoothly broken power law, characterized by a lowenergy index, α, a highenergy index, β, and a break energy E_{0}. They show a peak corresponding to a specific (and observable) value of the photon energy E_{p} = E_{0}(2 + α). For GRBs with measured spectrum and redshift it is possible to evaluate the intrinsic peak energy, E_{p,i} = E_{p}(1 + z) and the isotropic equivalent radiated energy
where N(E) is the Band function: (3)E_{iso} and E_{p,i} span several orders of magnitude (therefore GRBs cannot be considered standard candles), and show distributions approximated by Gaussians plus a tail at low energies. A strong correlation between these two quantities was initially discovered in a small sample ofBeppo SAXGRBs with known redshifts (Amati et al. 2002) and confirmed afterward by HETE2 and Swift observations Lamb et al. (2005), Amati (2006). Several analyses of the E_{p,i}E_{iso} plane of GRBs showed that different classes of GRBs exhibit different behaviors, and while normal long GRBs and Xray flashes (XRF, i.e., particularly soft bursts) follow this correlation, with the exception of the two peculiar subenergetic GRBs 980425 and 031203 see, e.g., Amati et al. (2007), for a discussion on possible explanations, short GRBs do not. These facts may depend on the different emission mechanisms and/or geometry involved in different classes of GRBs and makes this correlation a useful tool to distinguish between them (Amati 2006; Antonelli et al. 2009). Although it was the first strong correlation discovered for the GRB observables, until recent years, the E_{p,i}–E_{iso} correlation was never used for cosmology because it exhibits a significant dispersion around the bestfit power law: the residuals follow a Gaussian with a value of σ_{log Ep,i} ≃ 0.2. This type of additional Poissonian scatter is typically quantified by performing a maximum likelihood analysis that takes the variance and the errors on dependent and independent variables into account. By measuring E_{p,i} in keV and E_{iso} in 10^{52} erg, this method gives σ_{log Ep,i} = 0.19 ± 0.02 (Amati et al. 2008, 2009). However, the recent increase in the efficiency of GRB discoveries combined with the fact that the E_{p,i}–E_{iso} correlation requires only two parameters that are directly inferred from observations (this minimizes the effects of systematics and increases the number of GRBs that can be used) makes this correlation an interesting tool for cosmology. Despite the very large number of bursts consistent with this correlation, its physical origin is still debated. Some authors claimed that it is strongly affected by instrumental selection effects (Nakar & Piran 2005; Band & Preece 2005; Butler et al. 2007; Shahmoradi & Nemiroff 2011). However, many other studies found that such instrumental selection biases, even if they may affect the sample, cannot be responsible for the existence of the spectralenergy correlations (Ghirlanda et al. 2008; Nava et al. 2012). Moreover, a recent timeresolved spectral analysis of GRBs that were detected by the BeppoSAX and Fermi satellite showed that E_{p,i} correlates with the luminosity (e.g., Ghirlanda et al. 2010; Frontera et al. 2012) and radiated energy (e.g., Basak et al. 2013) also during the temporal evolution of the bursts (and the correlation between the spectral peak energy and the evolving flux has been pointed out by, e.g., Golenetskii et al. 1983 based on KonusWIND data). The timeresolved E_{p,i}luminosity and E_{p,i}–E_{iso} correlations within individual GRBs and that their average slope is consistent with that of correlations defined by the timeaveraged spectral properties of different bursts strongly supports the physical origin of the E_{p,i}–E_{iso} correlation. It is very difficult to explain them as a selection or instrumental effect (e.g., Dichiara & Amati 2013; Basak et al. 2013), and the predominant emission mechanism in GRB prompt emission produces a correlation between the spectral peak energy and intensity (which can be characterized either as luminosity, peak luminosity, or radiated energy). In addition to its existence and slope, an important property of the E_{p,i}–E_{iso} correlation is its dispersion. As shown by several works that were based on the socalled jetbreaks, in the optical afterglow light curves of some GRBs (e.g., Ghirlanda et al. 2004, 2006), 50% of the extraPoissonian scatter of the correlation is sometimes probably a result of the distribution of jet opening angles. However, these estimates of jet opening angles are still unconfirmed and modeldependent and can be made only for a small number of GRBs. Other factors that may contribute to the dispersion of the E_{p,i}–E_{iso} and other E_{p,i}intensity correlations include the jet structure, viewing angles, detectors sensitivity, and energy band (but see Amati et al. 2009), the diversity of shock microphysics parameters, and the magnetization within the emitting ejecta. At the current observational and theoretical status, it is very difficult to quantify these single contributions, which seem to act randomly in scattering the data around the best fit powerlaw (e.g., Ghirlanda et al. 2008; Amati et al. 2008). Importantly, it has been shown (e.g., Amati et al. 2009) that a small fraction (5–10%) of the scatter depends on the cosmological model and parameters assumed for computing E_{iso}, which makes this correlation a potential tool for cosmology.
In this section we show how the E_{p,i}–E_{iso} correlation can be calibrated to standardize long GRBs and to build a GRB Hubble diagram, which we use to investigate different cosmological models at very high redshifts (see also Amati & Della Valle 2013; Wang et al. 2015, 2016; Lin et al. 2015, 2016b).
2.1. GRB data
We used a sample of 162 long GRBs/XRFs as of September 2013 taken from the updated compilation of spectral and intensity parameters of GRBs by Sawant & Amati (2016). The redshift distribution of this sample covers a broad range, 0.03 ≤ z ≤ 9.3, which means that it extends far beyond the SNIa range (z ≤ 1.7). These data are of long GRBs/XRFs that are characterized by firm measurements of redshifts and the restframe peak energy E_{p,i}. The main contributions come from the joint detections by Swift/BAT and Fermi/GBM or KonusWIND, except for the small fraction of events for which Swift/BAT can directly provide E_{p,i} when it is in the (15−150) keV interval. For events detected by more than one of these detectors, the uncertainties on the E_{p,i} and E_{iso} take the measurements and uncertainties provided by each individual detector into account. In Table 4 we report for each GRB the redshift, the restframe spectral peak energy, E_{p,i}, and the isotropicequivalent radiated energy, E_{iso}. As detailed in (Sawant & Amati, in prep.), the criterion behind selecting the observations from a particular mission is based on the conditions summarized below.

1.
We preferred observations with exposure times of at leasttwothirds of the whole event duration. Hence KonusWINDand Fermi/GBM were chosen whenever available. For KonusWIND, the observations were reported in the official literature(see e.g. Ulanov et al. 2005),and also in GCN archives when needed. For Fermi/GBM, theobservations were taken from Gruber et al.(2012), from several other papers (e.g., Ghirlanda etal. 2004, 2005; Friedman& Bloom 2005, etc.), and from GCNs.

2.
The Swift/BAT observations were chosen when no other preferred mission (KonusWIND, Fermi/GBM) was able to provide spectral parameters and the value of E_{p,obs} was within the energy band of this instrument. In particular, the E_{p,i} values derived from BAT spectral analysis were conservatively taken from the results reported by the BAT team (Sakamoto et al. 2008a,b). Other BAT E_{p,i} values reported in the literature were not considered because they were not confirmed by Sakamoto et al. (2008a,b), by a refined analysis (see e.g. Cabrera et al. 2007), or because they were based on speculative methods Butler et al. (2007).
When available, the Band model Band et al. (1993) was considered as the cutoff power law, which sometimes overestimates the value of E_{p,i}. Finally, the error on any value was assumed to be not less than 10% to account for the instrumental capabilities and calibration uncertainties.
2.2. Cosmologically independent calibration: local regression of SNIa
The lack of nearby GRBs creates the socalled circularity problem: GRBs can be used as cosmological tools through the E_{p,i}–E_{iso} correlation, which is based on the cosmological model assumed for the computation of E_{iso}, however. In principle, this problem could be solved in several ways: it is possible, for instance, to simultaneously constrain the calibration parameters (a,b,σ_{int}) ∈ G and the set of cosmological parameters θ ∈ H by considering a likelihood function defined in the space G ⊗ H, which allows simultaneously fitting the parameters (e.g., Diaferio et al. 2011). This procedure is implemented in Sect. 5 and compared with the local regression technique. Alternatively, it has been proposed that the problem might be avoided by considering a sufficiently large number of GRBs within a small redshift bin centered on any zGhirlanda et al. (2006). However, this method, even if quite promising for the future, is currently unrealistic because of the limited number of GRBs with measured redshifts. In this section we implement a procedure for calibrating the E_{p,i}–E_{iso} relation in a way independent of the cosmological model by applying a local regression technique to estimate the distance modulus μ(z) from the recently updated SNIa sample, called Union2.1 (see also Kodama et al. 2008; Liang et al. 2008). Originally implemented by Cleveland & Devlin (1988), this technique combines the simplicity of linear regression with the flexibility of nonlinear regression to localized subsets of the data, and reconstructs a function describing their behavior in the neighborhood of any point z_{0}. A lowdegree polynomial is fitted to a subsample containing a neighborhood of z_{0}, by using weighted leastsquares with a weight function that quickly decreases with the distance from z_{0}. The local regression procedure can be schematically sketched as follows:
 1.
set the redshift z where μ(z) has to be reconstructed;
 2.
sort the SNIa Union2.1 sample by increasing value of  z−z_{i}  and select the first , where α is a userselected value and the total number of SNIa;
 3.
apply the weight function (4)where u =  z−z_{i}  /Δ and Δ the highest value of the  z−z_{i}  over the previously selected subset;
 4.
fit a firstorder polynomial to the data previously selected and weighted, and use the zerothorder term as the bestfit value of the modulus of distance μ(z);
 5.
evaluate the error σ_{μ} as the root mean square of the weighted residuals with respect to the bestfit value.
It is worth stressing that both the choice of the weight function and the order of the fitting polynomial are somewhat arbitrary. Similarly, the value of n to be used must not be too low to make up a statistically valuable sample, but also not too high to prevent the use of a loworder polynomial. To check our local regression routine, we simulated a large catalog with the same redshift and error distribution as the Union2.1 survey. We adopted a quintessence cosmological model with a constant EOS, w, and fixed values of the cosmological parameters (Ω_{M},w,h). For each redshift value in the Union2.1 sample, we extracted the corresponding modulus of distance from a Gaussian distribution centered on the theoretically predicted value and with the standard deviation σ = 0.15. To this value, we added the value of the error, corresponding to the same relative uncertainty as the data in the Union sample. This simulated catalog was used as input to the local regression routine, and the reconstructed μ(z) values were compared to the input ones.
Fig. 3
Bestfit curves for the E_{p,i}–E_{iso} correlation relation superimposed on the data. The solid and dashed lines refer to the results obtained with the maximum likelihood (Reichart likelihood) and weighted χ^{2} estimator, respectively. 
Defining the percentage deviation with μ_{rec.} and μ_{w} the local regression estimate and the input values, respectively, and averaging over 500 simulations, we found that the choice α = 0.02 gives (δμ/μ)_{rms} ≃ 0.3% with  ϵ  ≤ 1% independently on the redshift z. This result implies that the local regression method allows correctly reconstructing the underlying distance modulus regardless of redshift z< 1.4 from the Union SNIa sample. We also tested this results in different cosmological backgrounds by adopting an evolving EOS and averaging over five realizations of the mock catalog. With this efficient way of estimating μ(z) at redshift z in a modelindependent way, we can now fit the E_{p,i}–E_{iso} correlation relation, using the reconstructed μ(z). We considered only GRBs with z ≤ 1.414 to cover the same redshift range as is spanned by the SNIa data. To standardize the E_{p,i}–E_{iso} relation as expressed by Eq. (1), we need to fit a data array { x_{i},y_{i} } with uncertainties { σ_{x,i},σ_{y,i} }, to a straight line (5)and determine the parameters (a,b). We expect a certain amount of intrinsic additional Poissonian scatter, σ_{int}, around the bestfit line that has to be taken into account and determined together with (a,b) by the fitting procedure. We used a likelihood, implemented by Reichart et al. (2001), that solved this problem of fitting data that are affected by extrinsic scatter in addition to the intrinsic uncertainties along both axes: (6)where the sum is over the objects in the sample. We note that this maximization was performed in the twoparameter space (a,σ_{int}) since b may be calculated analytically by solving the equation , we obtained (7)To quantitatively estimate the goodness of this fit, we used the median and root mean square of the bestfit residuals, defined as δ = y_{obs}−y_{fit}. To quantify the uncertainties of some fit parameter p_{i}, we evaluated the marginalized likelihood ℒ_{i}(p_{i}) by integrating over the other parameters. The median value for the parameter p_{i} was then found by solving (8)The 68% (95%) confidence range (p_{i,l},p_{i,h}) was then found by solving e.g., D’Agostini (2005)(9)with ε = 0.68 and ε = 0.95 for the 68% and 95% confidence level. The maximum likelihood values of a and σ_{int} are and . In Fig. 3 we show the correlation log E_{p,i}−log E_{iso}. The solid line is the best fit obtained using the Reichart likelihood, and the dashed line is the best fit obtained by the weighted χ^{2} method. The marginalized likelihood functions are shown in Fig. 4.
Fig. 4
Marginalized likelihood functions: the likelihood function ℒ_{a} is obtained by marginalizing over σ_{int}; and the likelihood function, ℒ_{σint}, is obtained by marginalizing over a. 
As noted above, b can be evaluated analytically. We obtained b = 52.53 ± 0.02. It is worth noting that in the literature results for the inverse correlation are commonly reported, that is, the correlation E_{iso}–E_{p,i}: using the local regression technique and the Reichart likelihood, we also obtained this inverse calibration. We obtained , . In Fig. 5 we plot the bestfit curves for the E_{iso}–E_{p,i} correlation relation superimposed on the data.
Fig. 5
Bestfit curves for the E_{iso}–E_{p,i} correlation relation superimposed on the data. The solid and dashed lines refer to the results obtained with the maximum likelihood (Reichart likelihood) and weighted χ^{2} estimator, respectively. 
Fig. 6
Calibrated GRB Hubble diagram (black filled diamond) up to very high values of redshift, as constructed by applying the local regression technique with the Reichart likelihood: we show overplotted the behavior of the theoretical distance modulus μ(z) = 25 + 5log d_{L}(z) corresponding to the favored fitted CPL model (red line), with w_{0} = −0.29, w_{1} = −0.12, h = 0.74, Ω_{m} = 0.24, and the standard ΛCDM model (gray line) with Ω_{m} = 0.33, and h = 0.74. They are defined by Eqs. (12) and (21). The blue filled points at lower redshift are the SNIa data. 
2.2.1. Constructing the Hubble diagram
After fitting the correlation and estimating its parameters, we used them to construct the GRB Hubble diagram. We recall that the luminosity distance of a GRB with the redshift z may be computed as (10)The uncertainty of d_{L}(z) was then estimated through the propagation of the measurement errors on the involved quantities. In particular, recalling that our correlation relation can be written as a linear relation, as in Eq. (5), where y = E_{iso} is the distance dependent quantity, while x is not, the error on the distance dependent quantity y was estimated as (11)where σ_{b} is properly evaluated through the Eq. (27), which implicitly defines b as a function of a and σ_{int}, and is then added in quadrature to the uncertainties of the other terms entering Eq. (10) to obtain the total uncertainty. The distance modulus μ(z) is easily obtained from its standard definition (12)with its uncertainty obtained again by error propagation: (13)We finally estimated the distance modulus for each ith GRB in the sample at redshift z_{i} to build the Hubble diagram plotted in Fig. 6. We refer to this data set as the calibrated GRB Hubble diagram below since to compute the distances Hubble diagram we relied on the calibration based on the SNIa Hubble diagram. The derived distance moduli are divided into two subsets, listed in Tables 4 and 5. Table 4 lists GRBs with z ≤ 1.46, the same redshift range as for known SNIa, and Table 5 lists GRBs with z ≥ 1.47. In Fig. 7 we finally compare the GRB Hubble diagram (black points) with the SNIa Hubble diagram (blue points) and with BAO data (red points).
Fig. 7
GRB Hubble diagram (black points) is compared with the SNIa Hubble diagram (blue points) and with BAO data (red points). 
3. Cosmological constraints derived from the calibrated GRB Hubble diagram
Here we illustrate the possibilities of using the GRB Hubble diagram to constrain the cosmological models and investigate the dark energy EOS. Within the standard homogeneous and isotropic cosmology, the dark energy appears in the Friedman equations: where a is the scale factor, H = ȧ/a the Hubble parameter, ρ_{M} the density of matter, ρ_{de} the density of dark energy, p_{de} its pressure, and the dot denotes the time derivative. The continuity equation for any component of the cosmological fluid is (16)where the energy density is ρ_{i}, the pressure p_{i}, and the EOS of the ith component is defined by . The standard nonrelativistic matter has w = 0, and the cosmological constant has w = −1. The dark energy EOS and other constituents of the Universe determine the Hubble function H(z) and any derivations of it that are needed to obtain the observable quantities. When only matter and dark energy are present, the Hubble function is given by (17)where , and w(z,θ) describes the dark energy EOS, characterized by n parameters θ = (θ_{1},θ_{2},...,θ_{n}). We limit our analysis below to the CPL parametrization, (18)where w_{0} and w_{1} are constant parameters and they represent the w(z) present value and its overall time evolution, respectively Chevallier & Polarski (2001), Linder (2003). If we introduce the dimensionless Hubble parameter E(z,θ): (19)we can define the luminosity distance and the modulus of distance. Actually the luminosity distance id given by: (20)where , d_{M}(z,θ) is the transverse comoving distance and it is defined as (21)being d_{C}(z,θ) the comoving distance: (22)therefore we can define the modulus of distance μ_{th}(z,θ): (23)The standard ΛCDM model corresponds to w_{0} = −1, w_{1} = 0.
To constrain the parameters specifying different cosmological models, we maximized the likelihood function ℒ(θ) ∝ exp [−χ^{2}(θ)/2], where θ indicates the set of cosmological parameters and the χ^{2}(θ) was defined as usual by (24)Here, μ_{obs} and μ_{th} are the observed and theoretically predicted values of the distance modulus. The parameter space is efficiently sampled by using the MCMC method, thus running five parallel chains and using the GelmanRubin test to check the convergence Gelman & Rubin (1992). The confidence levels are estimated from the likelihood quantiles. We recall that we performed the analysis assuming a nonzero spatial curvature (not flat ΛCDM), and only in this case did we take w = −1. To alleviate the strong degeneracy of the curvature parameters and any EOS parameters, we added a Gaussian prior on Ω_{k}, centered on the value provided by the Planck Collaboration XVI (2014), , and with a dispersion ten times of , where . We investigated the CPL parameters assuming a flat universe. In a forthcoming paper we will present a full cosmological analysis using the highredshift GRB Hubble diagram to test different cosmological models, where several parametrizations of the dark energy EOS will be used and also different dark energy scenarios, for instance the scalar field quintessence. In Table 1 we summarize the results of our analysis. There are indications that the dark energy EOS is evolving.
Constraints on the cosmological parameters.
Fig. 8
Regions of confidence for the marginalized likelihood function ℒ(w_{0},w_{1}). On the axes we also plot the boxandwhisker diagrams for the w_{0} and w_{1} parameters: the bottom and top of the diagrams are, as commonly done in the boxplots, the 25th and 75th percentile (the lower and upper quartiles, respectively), and the band near the middle of the box is the 50th percentile (the median). It is evident that the ΛCDM model, corresponding to w_{0} = −1 and w_{1} = 0 is not favored. 
Constraints on the cosmological parameters.
The joint probability for different sets of parameters that characterize the CPL EOS is shown in Fig. 8.
4. E_{iso} and E_{p,i} correlation
Since we here used the calibrated GRB Hubble diagram to perform the cosmological investigation described above, we discuss some arguments about the reliability of using the E_{iso}−E_{p,i} relation for cosmological tasks. For instance, the calibration method we implemented so far relies on the underlying assumption that the calibration parameters do not evolve with redshift. It is worth noting that the problem of the redshift dependence of the GRBs correlations is widely debated in the literature, with different conclusions. This problem is intimately connected to the problem of the influence of possible selection effects or biases on the observed correlations, see, for instance, Li et al. (2008), Dainotti et al. (2013a,b), Liu & Wei (2015), Wang et al. (2016), Lin et al. (2016a,b). Any answer to these fundamental questions is far from being settled until more data with known redshifts are available. In this section, however, we revisit this question from an observational point of view: we test the validity of this commonly adopted working hypothesis and search for any evidence of such a redshift dependence. We also investigate possible effects on the GRB Hubble diagram. As a first simplified approach we considered two subsamples with a comparable number of bursts, divided according to redshift: a lower redshift sample of 97 bursts with z ≤ 2, and a higher redshift sample of 67 burst with z > 2. We estimated the cosmological parameters for flat ΛCDM and CPL dark energy EOS cosmological models, considering the two subsamples separately, and compared the results. Even when the bursts belonging to these two samples experience different environment conditions, we did not find significant indication that a spurious zevolution of the slope affects the cosmological fit: the values of the cosmological parameters derived from these two samples are statistically consistent within 1σ, as shown in Table 2. This result indirectly shows that redshift evolution dependence, even if it exists, does not undermine the reliability of the GRBs as probes of the cosmological expansion. Moreover, to investigate this redshift dependence in more detail, that is, how it affects E_{iso} and/or E_{p,i} (which we generically denote by y), we used two different approaches. First, we evaluated the Spearman rank correlation coefficient, C(z,y), taking into account that since our sample is not too large, a few points could dominate the final value of the rank correlation, which would introduce a bias. We applied a jacknife resampling method by evaluating C(z,y) for N−1 samples obtained by excluding one GRB at a time, and we adopted the mean value and the standard deviation to estimate C(z,y). We found C(z,E_{p,i}) = 0.299 ± 0.004, and C(z,E_{iso}) = 0.278 ± 0.004, which indicates a moderate evolution of the correlation. The Spearman rank correlation coefficient, however, does not include the errors of E_{iso} and E_{p,i}, so that we also implemented an alternative, and completely different, approach to determine whether and how strongly these variables are correlated with redshift. We assumed that the evolutionary functions can be parametrized by a simple powerlaw functions: g_{iso}(z) = ^{(}1 + z^{)}^{kiso} and g_{p}(z) = ^{(}1 + z^{)}^{kp} (see for instanceDainotti et al. 2013a), so that /g_{iso}(z) and /g_{p}(z) are the deevolved quantities. In this case, the effective E_{iso}−E_{p,i} correlation can be written as a 3D correlation:(25)Calibrating this 3D relation means determining the coefficients (a, b, k_{iso}, and k_{p}) plus the intrinsic scatter σ_{int}. Low values for k_{iso} and k_{p} would indicate a lack of evolution, or at least negligible evolutionary effects. To fit these coefficients, we constructed a 3D Reichart likelihood, but we consider no error on the redshift of each GRB: (26)where α = ak_{p}. As in the 2D case, we maximized our likelihood in the space (a, k_{iso}, and α) since b was evaluated analytically by solving the equation , we obtain(27)We also used the MCMC method and ran five parallel chains and the GelmanRubin convergence test, as previously explained. We finally studied the median and 68% confidence range of k_{iso} and α to test whether the correlation evolves and noted that a null value for these parameters is strong evidence for a lack of any evolution. We found that , k_{iso} = −0.04 ± 0.1; α = −0.02 ± 0.2; σ_{int} = 0.35 ± 0.03, so that . We can safely conclude that the E_{iso} and E_{p,i} correlation shows, at this stage, weak indications of evolution. In Fig. 9 we plot both the deevolved and evolved/original correlation, and the deevolved and evolved/original Hubble diagram: we do not see any signs of evolution.
5. Fully Bayesian analysis
Fig. 9
Evolutionary effects on the E_{iso} – E_{p,i} correlation and the Hubble flow. Upper panel: deevolved (red points) and evolved/original correlation (green poins); there is no discernable evolution. Bottom panel: deevolved (red points) and evolved/original (green points) Hubble diagram; evolutionary effects of the correlation do not affect the GRB Hubble diagram. 
In this section we simultaneously constrain the calibration parameters (a,b,σ_{int}) and the set of cosmological parameters by considering the same likelihood function as in Eq. (7). Our task is to determine the multidimensional probability distribution function (PDF) of the parameters { a,b,σ_{int},p }, where p is the dimensional vector of the cosmological parameters. The Amati correlation can be written in the form (28)We note that the bestfit zero point b can be analytically expressed as a function of (a,σ_{int}) when the cosmological parameter p_{c} is specified. A more computationally efficient strategy is to let b free and add it to the list of quantities to be determined, which thus sums up to . To efficiently sample the dimensional parameter space, we used the MCMC method and ran five parallel chains and used the GelmanRubin convergence test, as described in the previous section. Since we are mainly interested in the calibration problem and not in constraining the cosmological parameters (Ω_{M},Ω_{Λ},h), we considered here only the particular case of a flat ΛCDM model, as already done in the previous analysis. The analysis determines, at the same time, the cosmological parameters and the correlation coefficients, which are listed in Table 3.
Constraints on the calibration parameters.
Although the calibration procedure is different (since we now fit for the cosmological parameters as well), it is nevertheless worth comparing this determination of (a,b,σ_{int}) with the one obtained in the previous analysis based on the SNIa sample. It is evident that although the median values change, the 95% confidence levels are in full agreement so that we cannot find any statistically significant difference. Figure 10 shows the marginalized probability distribution function of the correlation coefficients of the relation E_{iso}−E_{p,i}.
As far as the cosmological parameters are concerned, it turns out that , the range of confidence at 1σ is (0.13,0.54), and , the range of confidence at 1σ is (0.50,0.87). This result implies that , and that the range of confidence is (− 0.00730,0.00605), in a good agreement with results derived by using the calibrated GRBs Hubble diagram.
Also in this case we finally estimate the distance modulus for each ith GRB in our sample at redshift z_{i}, to build the fiducial Hubble diagram, by using Eqs. (1) and (10). It turns out that the fiducial and calibrated Hubble diagrams are fully statistically consistent, as shown in Fig. 11.
Fig. 10
Regions of confidence for the marginalized likelihood function ℒ(a,σ), obtained by marginalizing over b and the cosmological parameters. The gray regions indicate the 1σ (full zone) and 2σ (empty zone) regions of confidence. On the axes we also plot the boxandwhisker diagrams for the a and σ_{int} parameters: the bottom and top of the diagrams are the 25th and 75th percentile (the lower and upper quartiles, respectively), and the band near the middle of the box is the 50th percentile (the median). 
6. Discussion and conclusions
The E_{p,i}–E_{iso} correlation is one of the most intriguing properties of the long GRBs, with significant implications for the use of GRBs as cosmological probes. Here we explored the Amati relation in a way independent of the cosmological model. Using the recently updated data set of 162 highredshift GRBs, we applied a local regression technique to estimate the distance modulus using the recent Union SNIa sample (Union2.1). The derived calibration parameters are statistically fully consistent with the results of our previous work Demianski et al. (2011). Moreover, we tested the validity of the commonly adopted working hypothesis that the GRB Hubble diagram is slightly affected by redshift dependence of the E_{p,i}–E_{iso} correlation. As a first qualitative and simplified approach we considered a lower redshift sample of 97 bursts with z ≤ 2, and a higher redshift sample of 67 burst with z > 2. We estimated the cosmological parameters for flat ΛCDM and CPL dark energy EOS cosmological models, considering the two subsamples separately, and compared the results. Even when the bursts belonging to these two samples had different environment conditions, we found no significant indications that a spurious zevolution of the slope affected the cosmological fit. Moreover, to quantify this redshift dependence, we used two different approaches. First, we evaluated the Spearman rank correlation coefficient, C(z,y) by applying a jacknife resampling method by evaluating C(z,y) for N−1 samples obtained by excluding one GRB at a time, and we adopted the mean value and the standard deviation to estimate C(z,y). C(z,E_{p,i}) = 0.299 ± 0.004, and C(z,E_{iso}) = 0.278 ± 0.004, which indicates a negligible evolution of the correlation. Moreover, we also implemented an alternative method, assuming that the redshift evolution can be parametrized by simple powerlaw functions: g_{iso}(z) = ^{(}1 + z^{)}^{kiso} and g_{p}(z) = ^{(}1 + z^{)}^{kp} and that the correlation holds for the deevolved quantities and . In this case, we rewrote an effective 3D E_{iso}−E_{p,i} correlation, with included the evolutionary terms. Since we were interested in the implications of possible evolutionary effects of the GRB Hubble diagram and to simplify the fit, we introduced an auxiliary variable α = ak_{p}. A null value for k_{iso} and α is strong evidence for a lack of any evolution.
Fig. 11
Fiducial GRB Hubble diagram, superimposed on the calibrated Hubble diagram. 
To fit the coefficient of our 3D correlation we constructed a 3D Reichart likelihood and used a MCMC method running five parallel chains and using the GelmanRubin convergence test. The E_{iso} and E_{p,i} correlation shows at this stage only weak indications of evolution. The derived calibration parameters were used to construct an updated calibrated GRB Hubble diagram, which we adopted as a tool to constrain the cosmological models and to investigate the dark energy EOS. In particular, we searched for any indications that w(z) ≠ −1, which reflects the possibility of a deviation from the ΛCDM cosmological model. To accomplish this task, we focused on the CPL parametrization as an example. To efficiently sample the cosmological parameter space, we again used a MCMC method. At 1σ level we found indications for a time evolution of the EOS in the considered parametrization; we conclude that the ΛCDM model is not favored even though it is not ruled out by these observations so far. Moreover, for w = −1 we also performed our analysis assuming a nonzero spatial curvature, adding a Gaussian prior on Ω_{k} centered on the value provided by the Planck data Planck Collaboration XVI (2014). The GRB Hubble diagram alone provides for the range of confidence at 1σ(− 0.007,0.0064). Finally, to investigate the reliability of the E_{p,i}–E_{iso} correlation in greater detail, we also used a different method to simultaneously extract the correlation coefficients and the cosmological parameters of the model from the observed quantities. To illustrate this method we assumed here, by way of an example, only the particular case of a nonflat ΛCDM model, as already done in the previous analysis. This analysis simultaneously determines the cosmological parameters and the correlation coefficients. Although the calibration procedure is different (since we now fit for the cosmological parameters as well), it is nevertheless worth comparing this determination of (a,b,σ_{int}) with the one obtained in the previous analysis based on the SNIa sample. It is evident that although the median values change, the 95% confidence levels are in full agreement so that we cannot find any statistically significant difference. This means that the highredshift GRBs can be used as cosmological probes, mainly in a redshift region, , which is unexplored by SNIa and BAO samples, and that both the calibration technique based on a local regression with SNIa and a fully Bayesian approach are reliable. As a final remark we note that our results for the cosmological parameters are statistically consistent with those previously obtained by other teams, as shown in Table 6, where we display some of the most recent measurements of cosmological parameters obtained with GRBs, even if following a slightly different approach. The main peculiarity in our analysis consists of the procedure used to build up the dataset, consisting of 162 objects, as discussed in Sect. 2.1 (see also Sawant & Amati, in prep.), even more than the specific statistical analysis, which follows a Bayesian approach through the implementation of the MCMC. Because the reliability of GRBs as distance indicators has not yet been clearly proved and because it constitutes an independent topic in the field of GRB research, it is not meaningless to investigate the reliability of the E_{p,i}–E_{iso} correlation whenever new and improved datasets are available. In the near future we intend to enhance the cosmological analysis by using the highredshift GRB Hubble diagram to test different cosmological models, where several parametrizations of the dark energy EOS are used, but a cosmographic approach will also be implemented, which will update the analysis performed in Demianski et al. (2012) to check whether with this new updated dataset the estimates of the jerk and the dark energy parameters will confirm deviations from ΛCDM cosmological model, as has been indicated in our previous analysis.
Derived distance moduli for the low redshift subsample, consisting of the GRBs whose redshifts cross the range spanned by the SNIa.
Derived distance moduli for the high redshift subsample, consisting of the GRBs whose redshifts pass through the range spanned by the SNIa.
Some of the most recent constraints of cosmological parameters by GRBs.
Acknowledgments
M.D. is grateful to the INFN for financial support through the Fondi FAI GrIV. E.P. acknowledges the support of INFN Sez. di Napoli (Iniziative Specifica QGSKY and TEONGRAV). L.A. acknowledges support by the Italian Ministry for Education, University and Research through PRIN MIUR 2009 project on “Gamma ray bursts: from progenitors to the physics of the prompt emission process”. (Prot. 2009 ERC3HT.)
References
 Alam, U., Sahni, V., Saini, T. D., & Starobinsky, A. A. 2003, MNRAS, 344, 1057 [NASA ADS] [CrossRef] [Google Scholar]
 Amati, L. 2006, MNRAS, 372, 233 [NASA ADS] [CrossRef] [Google Scholar]
 Amati, L., Frontera, F., Tavani, M., et al. 2002, A&A, 390, 81 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Amati, L., Della Valle, D., Frontera, F., et al. 2007, A&A, 463, 913 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Amati, L., Guidorzi, C., Frontera, F., et al. 2008, MNRAS, 391, 577 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Amati, L., Frontera, F., & Guidorzi, C. 2009, A&A, 508, 173 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Amati, L., & Della Valle, M. 2013, Int. J. Mod. Phys. D, 22, 30028 [Google Scholar]
 Amanullah, R., Lidman, C., Rubin, D., et al. 2010, ApJ, 71, 712 [NASA ADS] [CrossRef] [Google Scholar]
 Antonelli, A. L., D’Avanzo, P., Perna, R., et al. 2009, A&A, 507, L45 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Astier, P., Guy, J., Regnault, N., et al. 2006, A&A, 447, 31 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Band, D. L., & Preece, R. D. 2005, ApJ, 627, 319 [NASA ADS] [CrossRef] [Google Scholar]
 Band, D., Matteson, J., Ford, L., et al. 1993, ApJ, 413, 281 [NASA ADS] [CrossRef] [Google Scholar]
 Basak, R., & Rao, A. R. 2013, MNRAS, 436, 3082 [NASA ADS] [CrossRef] [Google Scholar]
 Butler, N. R., Kocevski, D., Bloom, J. S., et al. 2007, ApJ, 671, 656 [NASA ADS] [CrossRef] [Google Scholar]
 Cabrera, J. I., Firmani, C., AvilaReese, V., et al. 2007, MNRAS, 382, 342 [NASA ADS] [CrossRef] [Google Scholar]
 Cardone, V. F., Capozziello, S., & Dainotti, M. G. 2008, MNRAS, 391, L79 [NASA ADS] [CrossRef] [Google Scholar]
 Carroll, S. M. 2001, Liv. Rev. Relativ., 4, 1 [Google Scholar]
 Chevallier, M., & Polarski, D. 2001, Int. J. Mod. Phys. D, 10, 213 [NASA ADS] [CrossRef] [Google Scholar]
 Clarkson, C., & Maartens, R. 2010, Class. Quant. Gravity, 27, 124008 [Google Scholar]
 Cucchiara, A., Levan, A. J., Fox, D. B., et al. 2011, ApJ, 736, 7 [NASA ADS] [CrossRef] [Google Scholar]
 D’Agostini, G. 2005, [arXiv:physics/0511182] [Google Scholar]
 Dai, Z. G., Liang, E. W., & Xu, D. 2004, ApJ, 612, L101 [NASA ADS] [CrossRef] [Google Scholar]
 Dainotti, M. G., Petrosian, V., Singal, J., & Ostrowski, M. 2013a, ApJ, 774, 157 [NASA ADS] [CrossRef] [Google Scholar]
 Dainotti, M. G., Cardone, V. F., Piedipalumbo, E., & Capozziello, S. 2013b, MNRAS, 436, 82 [NASA ADS] [CrossRef] [Google Scholar]
 Diaferio, A., Ostorero, L., & Cardone, V. F. 2011, JCAP, 10, 008 [NASA ADS] [CrossRef] [Google Scholar]
 Dichiara, S., & Amati, L. 2013, Acta Politechnica, 53, 687 [Google Scholar]
 Demianski, M., Piedipalumbo, E., & Rubano, C. 2011, MNRAS, 411, 1213 [NASA ADS] [CrossRef] [Google Scholar]
 Demianski, M., Piedipalumbo, E., Rubano, C., & Scudellaro, P. 2012, MNRAS, 426, 13961415 [NASA ADS] [CrossRef] [Google Scholar]
 Friedman, A. S., & Bloom, J. S. 2005, ApJ, 627, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Frontera, F., Amati, L., Guidorzi, C., Landi, R., & in’t Zand, J. 2012, ApJ, 754, 138 [NASA ADS] [CrossRef] [Google Scholar]
 Gao, H., Liang, N., & Zhu, Z.H. 2012, Int. J. Mod. Phys. D, 21, 12500161 [NASA ADS] [CrossRef] [Google Scholar]
 Gelman, A., & Rubin, D. B. 1992, Stat. Sci., 7, 457 [Google Scholar]
 Ghirlanda, G., Ghisellini, G., & Lazzati, D. 2004, ApJ, 616, 331 [NASA ADS] [CrossRef] [Google Scholar]
 Ghirlanda, G., Ghisellini, G., Firmani, C., Celotti, A., & Bosnjak, Z. 2005, MNRAS, 360, L45 [NASA ADS] [CrossRef] [Google Scholar]
 Ghirlanda, G., Ghisellini, G., & Firmani, C. 2006, New J. Phys., 8, 123 [NASA ADS] [CrossRef] [Google Scholar]
 Ghirlanda, G., Nava, L., Ghisellini, G., Firmani, C., & Cabrera, J. I. 2008, MNRAS, 387, 319 [NASA ADS] [CrossRef] [Google Scholar]
 Ghirlanda, G., Nava, L., & Ghisellini, G. 2010, A&A, 511, A43 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Golenetskii, S. V., Mazets, E. P., Aptekar, R. L., & Ilinskii, V. N. 1983, Nature, 306, 451 [NASA ADS] [CrossRef] [Google Scholar]
 Kodama, Y., Yonetoku, D., Murakami, T., et al. 2008, MNRAS, 391, L1 [NASA ADS] [CrossRef] [Google Scholar]
 Lamb, D. Q., Donaghy, T. Q., & Graziani, C. 2005, ApJ, 620, 355 [NASA ADS] [CrossRef] [Google Scholar]
 Izzo, L., Muccino, M., Zaninoni, E., Amati, L., & Della Valle, M. 2015, A&A, 582, A115 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Liang, N., Xiao, W. K., Liu, Y., & Zhang, S. N. 2008, ApJ, 685, 354 [NASA ADS] [CrossRef] [Google Scholar]
 Lin, H.N., Li, X., Liu, Y., Wang, S., & Chang, Z. 2015, MNRAS, 453, 128 [NASA ADS] [CrossRef] [Google Scholar]
 Lin, H.N., Li, X., & Chang, Z. 2016a, MNRAS, 455, 2131 [NASA ADS] [CrossRef] [Google Scholar]
 Lin, H.N., Li, X., & Chang, Z. 2016b, MNRAS, 459, 2501 [NASA ADS] [CrossRef] [Google Scholar]
 Liu, J., & Wei, H. 2015, General Relativity and Gravitation, 47, 141 [NASA ADS] [CrossRef] [Google Scholar]
 Linder, E. V. 2003, Phys. Rev. Lett., 90, 091301 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Meszaros, P. 2006, Rep. Prog. Phys., 69, 2259 [NASA ADS] [CrossRef] [Google Scholar]
 Nakar, E., & Piran, T. 2005, MNRAS, 360, L36 [Google Scholar]
 Nava, L., Salvaterra, R., Ghirlanda, G., et al. 2012, MNRAS, 421, 1256 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Perlmutter, S., Aldering, G., della Valle, M., et al. 1998, Nature, 391, 51 [Google Scholar]
 Perlmutter, S., Aldering, G., Goldhaber, G., et al. 1999, ApJ, 517, 565 [NASA ADS] [CrossRef] [Google Scholar]
 Planck Collaboration XVI. 2014, A&A, 571, A16 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XIII. 2016, A&A, 594, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Postnikov, V. S., Dainotti, M. G., Hernandez, X., & Capozziello, S. 2014, ApJ, 783, 126 [NASA ADS] [CrossRef] [Google Scholar]
 Reichart, D. E., Lamb, D. Q., Fenimore, E. E., et al. 2001, ApJ, 552, 57 [NASA ADS] [CrossRef] [Google Scholar]
 Riess, A. G., Filippenko, A. V., Challis, P., et al. 1998, AJ, 116, 1009 [NASA ADS] [CrossRef] [Google Scholar]
 Riess, A. G., Strolger, L.G., Tonry, J., et al. 2004, ApJ, 607, 665 [NASA ADS] [CrossRef] [Google Scholar]
 Riess, A. G., Strolger, L. G., Casertano, S., et al. 2007, ApJ, 659, 98 [NASA ADS] [CrossRef] [Google Scholar]
 Sahni, V., Saini, T. D., Starobinsky, A. A., & Alam, U. 2003, JETP Lett., 77, 201 [Google Scholar]
 Sakamoto, T., Hullinger, D., Sato, G., Yamazaki, R., & Barbier, L. 2008a, ApJ, 679, 570 [NASA ADS] [CrossRef] [Google Scholar]
 Sakamoto, T., Barthelmy, S. D., Barbier, L., et al. 2008b, ApJS, 175, 179 [NASA ADS] [CrossRef] [Google Scholar]
 Salvaterra, R., Della Valle, M., Campana, S., et al. 2009, Nature, 461, 1258 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Schmidt, B. P., Suntzeff, N. B., Phillips, M. M., et al. 1998, ApJ, 507, 46 [Google Scholar]
 Shahmoradi, A., & Nemiroff, R. J. 2011, MNRAS, 411, 1843 [NASA ADS] [CrossRef] [Google Scholar]
 Spergel, D. N., Bean, R., Doré, O., et al. 2007, ApJS, 170, 377 [NASA ADS] [CrossRef] [Google Scholar]
 Tanvir, N. R., Fox, D. B., Levan, A. J., et al. 2009, Nature, 461, 1254 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Ulanov, M. V., Golenetskii, S. V., Frederiks, D. D., et al. 2005, Nuovo Cimento C, 28, 351 [Google Scholar]
 Wang, F. Y., Dai, Z. G., & Liang, E. W. 2015, New Astron. Rev., 67, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Wang, J. S., Wang, F. Y., Cheng, K. S., & Dai, Z. G. 2016, A&A, 585, A68 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
All Tables
Derived distance moduli for the low redshift subsample, consisting of the GRBs whose redshifts cross the range spanned by the SNIa.
Derived distance moduli for the high redshift subsample, consisting of the GRBs whose redshifts pass through the range spanned by the SNIa.
All Figures
Fig. 1
Dependence of the distance modulus on EOS of dark energy. Upper panel: distance modulus μ(z) for different values of the EOS parameters. The black line represents the standard flat ΛCDM model with Ω_{m} = 0.3, h = 0.7, w_{0} = −1, w_{1} = 0. The other curves correspond to different values of the CPL parameters w_{0} and w_{1}. Bottom panel: z dependence of the percentage error in the distance modulus between the fiducial ΛCDM model and different flat CPL models, as in the left panel. 

In the text 
Fig. 2
Redshift distribution of our sample of 162 long GRBs/XRFs and the Union 2.1 SNIa dataset. 

In the text 
Fig. 3
Bestfit curves for the E_{p,i}–E_{iso} correlation relation superimposed on the data. The solid and dashed lines refer to the results obtained with the maximum likelihood (Reichart likelihood) and weighted χ^{2} estimator, respectively. 

In the text 
Fig. 4
Marginalized likelihood functions: the likelihood function ℒ_{a} is obtained by marginalizing over σ_{int}; and the likelihood function, ℒ_{σint}, is obtained by marginalizing over a. 

In the text 
Fig. 5
Bestfit curves for the E_{iso}–E_{p,i} correlation relation superimposed on the data. The solid and dashed lines refer to the results obtained with the maximum likelihood (Reichart likelihood) and weighted χ^{2} estimator, respectively. 

In the text 
Fig. 6
Calibrated GRB Hubble diagram (black filled diamond) up to very high values of redshift, as constructed by applying the local regression technique with the Reichart likelihood: we show overplotted the behavior of the theoretical distance modulus μ(z) = 25 + 5log d_{L}(z) corresponding to the favored fitted CPL model (red line), with w_{0} = −0.29, w_{1} = −0.12, h = 0.74, Ω_{m} = 0.24, and the standard ΛCDM model (gray line) with Ω_{m} = 0.33, and h = 0.74. They are defined by Eqs. (12) and (21). The blue filled points at lower redshift are the SNIa data. 

In the text 
Fig. 7
GRB Hubble diagram (black points) is compared with the SNIa Hubble diagram (blue points) and with BAO data (red points). 

In the text 
Fig. 8
Regions of confidence for the marginalized likelihood function ℒ(w_{0},w_{1}). On the axes we also plot the boxandwhisker diagrams for the w_{0} and w_{1} parameters: the bottom and top of the diagrams are, as commonly done in the boxplots, the 25th and 75th percentile (the lower and upper quartiles, respectively), and the band near the middle of the box is the 50th percentile (the median). It is evident that the ΛCDM model, corresponding to w_{0} = −1 and w_{1} = 0 is not favored. 

In the text 
Fig. 9
Evolutionary effects on the E_{iso} – E_{p,i} correlation and the Hubble flow. Upper panel: deevolved (red points) and evolved/original correlation (green poins); there is no discernable evolution. Bottom panel: deevolved (red points) and evolved/original (green points) Hubble diagram; evolutionary effects of the correlation do not affect the GRB Hubble diagram. 

In the text 
Fig. 10
Regions of confidence for the marginalized likelihood function ℒ(a,σ), obtained by marginalizing over b and the cosmological parameters. The gray regions indicate the 1σ (full zone) and 2σ (empty zone) regions of confidence. On the axes we also plot the boxandwhisker diagrams for the a and σ_{int} parameters: the bottom and top of the diagrams are the 25th and 75th percentile (the lower and upper quartiles, respectively), and the band near the middle of the box is the 50th percentile (the median). 

In the text 
Fig. 11
Fiducial GRB Hubble diagram, superimposed on the calibrated Hubble diagram. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext 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 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.