Contents

A&A 407, L1-L4 (2003)
DOI: 10.1051/0004-6361:20030958

A simple empirical redshift indicator for gamma-ray bursts

J.-L. Atteia

Laboratoire d'Astrophysique, Observatoire Midi-Pyrénées, 31400 Toulouse, France

Received 3 April 2003 / Accepted 23 June 203

Abstract
We propose a new empirical redshift indicator for gamma-ray bursts. This indicator is easily computed from the gamma-ray burst spectral parameters and its duration, and it provides "pseudo-redshifts'' accurate to a factor two. Possible applications of this redshift indicator are briefly discussed.

Key words: gamma-rays: bursts

   
1 Introduction

Gamma-ray bursts (GRBs) are huge stellar explosions which have been observed at redshifts ranging from 0.0085 to 4.5. While GRBs are in principle detectable out to very large redshifts ( $z=10\hspace{-0.4mm}-\hspace{-0.7mm}20$, Lamb & Reichart 2000), redshifts measured to date do not exceed 4.5. The method most frequently used to measure GRB redshifts is to find a visible afterglow, and to identify absorption lines in its spectrum, caused by the gas in the GRB host galaxy. The redshift of the host can also be measured at late times from the host emission lines, when the afterglow has faded below detection. Another, less frequent, method uses X-ray lines detected in the X-ray afterglows of some GRBs. The absence of GRB detection beyond z=5 could be explained by the fact that the afterglows of such distant GRBs must be searched for in the infrared, due to the Lyman alpha cutoff. The difficulty to measure spectroscopic redshifts led various authors to propose alternate ways to determine GRB redshifts. Norris et al. (2000) and Reichart et al. (2001) have found empirical luminosity estimators based on GRB light curves. Such luminosity estimators can be used to infer the intrinsic luminosity of individual GRBs, and consequently their redshifts. While these estimators cannot be used to obtain precise redshifts for individual GRBs, they are useful to derive statistical properties of the GRB population.

Redshift estimators based on the gamma-ray data only present two distinctive advantages: they provide redshift estimates for most GRBs detected in gamma-rays, and they do not require extensive follow-up campaigns involving large telescopes on the ground or in space. Important issues can be addressed with moderatly accurate redshifts, like the amount of energy released by GRBs in gamma-rays, the luminosity function of GRBs, or the history of stellar formation at high redshifts.

We propose here a new method to obtain redshift indicators for GRBs from gamma-ray observations. Our method is calibrated with 17 GRBs detected with BeppoSAX (Boella et al. 1997) and HETE (Ricker et al. 2001). In the following we call the redshifts inferred from our redshift indicator "pseudo-redshifts''. Pseudo-redshifts have the advantage of being very easily computed. In addition to the possible applications already mentioned, pseudo-redshifts may become a useful tool to quickly identify high-redshift GRBs.

   
2 An empirical redshift indicator

Finding redshift indicators for GRBs based on the gamma-ray data alone has always faced the problem of the large intrinsic dispersion of GRB properties. This intrinsic dispersion prevents the determination of the redshifts of individual GRBs. With the measure of an increasing number of GRB redshifts it appeared, however, that several properties of GRBs are correlated with the isotropic-equivalent energy radiated in gamma-rays (called  $E_{{\rm rad}}$ in the following). For instance, the correlation of the spectral hardness with  $E_{{\rm rad}}$ has been suspected for a long time (see e.g. Atteia 2000; Lloyd et al. 2000). It has only been demonstrated recently by Amati et al. (2002) for 12 GRBs with known redshifts. The correlation of the duration with  $E_{{\rm rad}}$ is discussed in Lee et al. (2000). These correlations have led some authors to propose using the observed GRB properties to infer  $E_{{\rm rad}}$, and then to deduce the redshift from the comparison of the observed fluence with  $E_{{\rm rad}}$. Norris et al. (2000), for instance, estimate  $E_{{\rm rad}}$ from the magnitude of the time lags between a high energy band and a low energy band. Reichart et al. (2001) estimate  $E_{{\rm rad}}$ from the variability of the light curve. We propose and test here another approach: we search a quantity which depends little on  $E_{{\rm rad}}$, and which has a small intrinsic dispersion which does not blur the dependence on redshift. Starting from empirical considerations, We find such a quantity essentially based on the spectral characteristics of GRBs.

GRB energy spectra are well fit with the so-called GRB model, consisting of two smoothly connected power laws (Band et al. 1993). In the following, $\alpha $ is the index of the low-energy power law, $\beta $ the index of the high-energy power law, and E0 is the break energy. With this parametrization, the peak energy of the $\nu f_{\nu}$ spectrum is $E_{\rm p} = E_0 \times (2+\alpha)$. $E_{\rm p}$ is well defined for $\alpha \ge -2$ and $\beta < -2$.

Our method is based on the recent finding by Amati et al. (2002) of a correlation between the intrinsic (redshift corrected) $E_{\rm p}$ of 12 GRBs with known redshifts, and $E_{{\rm rad}}$, their isotropic-equivalent energy radiated in gamma-rays. According to Amati et al., $E_{\rm p}$ is roughly proportional to the square-root of  $E_{{\rm rad}}$. Since $\alpha $ and $\beta $ do not vary too much from burst to burst, and since the energy radiated in gamma-rays is more or less the product of the number of photons by their typical energy, we make the assumption that the isotropic-equivalent number of photons in a GRB, $N_{\rm\gamma }$, is also roughly proportional to the square-root of $E_{{\rm rad}}$. For this study, we define $N_{\rm\gamma }$ as the number of photons below the break, integrated from  $E_{\rm p}$/100 to  $E_{\rm p}$/2 .


  \begin{figure}
\resizebox{8.8cm}{!}{\includegraphics[clip]{Fd031_f1.eps}}
\end{figure} Figure 1: Dependence of two characteristic GRB quantities on  $E_{{\rm rad}}$, the isotropic-equivalent energy radiated in gamma-rays. The stars show the dependence of the intrinsic peak energy (multiplied here by  $1.4\times 10^{57}$). The circles show the dependence of the isotropic-equivalent number of photons emitted by the source.


  \begin{figure}
\resizebox{8.8cm}{!}{\includegraphics[clip]{Fd031_f2.eps}}
\end{figure} Figure 2: The ratio $N_{\rm\gamma }$/ $E_{\rm p}$ (see text) as a function of $E_{{\rm rad}}$. This figure illustrates the weak dependence of $N_{\rm\gamma }$/ $E_{\rm p}$ with $E_{{\rm rad}}$.


  
Table 1: Observed properties of 17 GRBs with known redshift. The ten columns give the GRB name, the duration T90 in seconds, the three spectral parameters ( $\alpha , \beta $, and E0), the gamma-ray fluence  $S_{\rm\gamma }$ in units of 10-6 erg cm-2, the spectroscopic redshift z, the pseudo-redshift $\hat{z}$, the ratio $\hat{z}$/z, and a reference for the spectral parameters.
\begin{table}\begin{displaymath}
\begin{array}{llllllllll}
\hline\hline
\noali...
....99 & 7 \\
\noalign{\smallskip }
\hline
\end{array}\end{displaymath}\end{table}
$^{{\it a}}$ $\beta $ has been frozen to -2.3 for HETE GRBs 010921 to 030228.
1 Amati et al. (2002). Fluence measured in the range 40-700 keV.
2 Andersen et al. (2000). Fluence measured in the range 28-1800 keV.
3 Barraud et al. (2003). Fluence measured in the range 30-400 keV.
4 Sakamoto et al. (2003). Fluence measured in the range 7-30 keV.
5 Crew et al. (2003). Fluence measured in the range 7-30 keV.
6 Lamb et al. (2003). Fluence measured in the range 30-400 keV.
7 Vanderspek et al. (2003). Fluence measured in the range 7-30 keV.

Figure 1 shows $E_{\rm p}$, the intrinsic peak energy, and $N_{\rm\gamma }$, the isotropic-equivalent number of photons as a function of $E_{{\rm rad}}$ for a sample of 17 GRBs detected by BeppoSAX, BATSE, and HETE. The main characteristics of these GRBs are given in Table 1, along with references for their spectral parameters. The redshifts have been taken from J. Greiner's GRB page at http://www.mpe.mpg.de/~jcg/grbgen.html (except for GRB 020124, which comes from Hjorth et al. 2003). Figure 1 shows that, as we suspected, $E_{\rm p}$ and  $N_{\rm\gamma }$ have roughly the same dependence on  $E_{{\rm rad}}$. We can thus go one step further with our main conjecture: we suppose that the ratio  $N_{\rm\gamma }$/$E_{\rm p}$ is almost independent of  $E_{{\rm rad}}$, and can be used as a redshift indicator. Figure 2 shows that indeed the ratio  $N_{\rm\gamma }$/$E_{\rm p}$ shows very little dependence on  $E_{{\rm rad}}$, confirming our conjecture. This is not sufficient, however, to make it a correct redshift indicator. The critical issue is to find an indicator which has a small dependence on  $E_{{\rm rad}}$, a strong dependence on the redshift, and a not too strong intrinsic dispersion. This issue is discussed in the next section.

   
2.1 Definition of the redshift indicator

The theoretical considerations in the previous section are based on the study of intrinsic GRB properties. Defining a redshift indicator implies that we do not know the redshifts of the GRBs which are being studied, but only their observed properties. To keep in mind this difference, in all the following we use capital letters for intrinsic quantities, and lower case for observed quantities.

As discussed earlier, the best redshift indicator is not necessarily the one with the smallest intrinsic dispersion, but rather the one which has the best combination of a small intrinsic dispersion and a large dependence on redshift. Relying on the analysis of the previous section, we propose to base our redshift indicator on  $n_{\rm\gamma }$/$e_{\rm p}$  the ratio of the observed number of photons in the GRB on the observed peak energy. We tried various simple combinations of GRB parameters, all involving the ratio  $n_{\rm\gamma }$/$e_{\rm p}$, and found that X $n_{\rm\gamma }$/$e_{\rm p}$/ $\sqrt {{t}_{90}}$ has the right combination of properties for a redshift indicator. In this equation $e_{\rm p}$ is the observed peak energy, $n_{\rm\gamma }$ the observed number of photons between $e_{\rm p}$/100 and $e_{\rm p}$/2, and t90 the observed duration. We do not claim here that X is definitely the best redshift indicator, we nevertheless believe that it is sufficiently good to deserve further discussion.

We derive pseudo-redshifts from the measure of X in the following way: in a first step we compute the theoretical evolution of X with redshift; then we invert this relation to derive a pseudo-redshift from the observed value of X. The evolution of X with redshift can be written as

\begin{displaymath}X = A \times f(z)\end{displaymath}

A is a constant of normalization, and f describes the evolution of X with redshift for a "standard'' GRB ( $\alpha = -1.0$, $\beta = -2.3$, and E0 = 250 keV) in a "standard'' universe (H0 = 65 km s-1 Mpc-1, $\Omega_0 = 0.3$, $\Omega_\Lambda = 0.7$). GRB spectral parameters are not critical here, because we have shown in the previous section that the ratio  $N_{\rm\gamma }$/$E_{\rm p}$ does not vary much from burst to burst. The normalization constant A has been chosen to have about the same number of GRBs below and above the theoretical curve in Fig. 3 (A=60).

Pseudo-redshifts $\hat{z}$ are then defined by

\begin{displaymath}\hat{z} = {1 \over A} \times f^{-1}(X).
\end{displaymath} (1)

Their use as redshift indicators is discussed below.

   
2.2 Evaluation of pseudo-redshifts

Figure 3 shows the values of Xas a function of z, for the 17 GRBs of Table 1. This figure displays a clear anticorrelation between the two quantities. The dotted line indicates the theoretical dependence $X = A \times f(z)$. The coefficient of correlation between z and X is -0.875, corresponding to a correlation significant at the level of 4.9 sigmas using Fisher's Z transformation. We consider that this anticorrelation provides a good support to our intention of using X as a redshift indicator, and we use the Eq. (1) above to compute the pseudo-redshifts of GRBs in Table 1. Table 1 gives the values of z, $\hat{z}$, and their ratio, for the 17 GRBs with known redshift used in our analysis[*]. It shows that $\hat{z}$ is usually within a factor of two of z, except for GRB 000131 (at z=4.5), for which z and $\hat{z}$ differ by a factor 3.3. This discrepancy could be the consequence of the low quality of our redshift indicator for this burst (most probably) or of a problem with the measure of the spectral parameters of this GRB or of its redshift. It might also indicate that the relation between z and $\hat{z}$ is only working (or calibrated) to z=3.5. Because this event is clearly an outlier, we recomputed the coefficient of correlation between z and X without it. We find a coefficient of correlation of -0.927, corresponding to a correlation significant at the level of 6.1 sigmas using Fisher's Z transformation.

We conclude that the intrinsic dispersion of X is not such that it prevents its use as a redshift indicator. We prefer the term redshift indicator than redshift estimator, because the ratio of $\hat{z}$ over z varies too much for a true redshift estimator. In the following we use the name pseudo-redshifts for $\hat{z}$. Because $\hat{z}$ was derived from a purely empirical approach, we expect that an approach based on a physical treatment of GRB emission might provide a better redshift estimator. While pseudo-redshifts are a potentially useful tool, they deserve a word of caution here because we do not know how observational biases affect Fig. 3. We note for instance that GRBs with spectroscopic redshifts certainly represent a biased sample. In addition the relations linking the GRB properties (from which we compute X) to those of their afterglows (from which we measure z) are far from being understood. Anyone using this tool should thus keep in mind that Fig. 3 provides a biased view of the true distribution of GRBs in the (zX) plane.


  \begin{figure}
\resizebox{8.8cm}{!}{\includegraphics[clip]{Fd031_f3.eps}}
\end{figure} Figure 3: Correlation of X =  $n_{\rm\gamma }$/ $e_{\rm p}$/ $\sqrt {{t}_{90}}$ (see text ) with the measured redshifts of 17 GRBs. The isolated star at z=4.5 is GRB 000131. The dotted line shows the theoretical relation between $\hat{z}$ and z.

   
3 An example of using pseudo-redshifts

In this section we compute the pseudo-redshifts of 18 GRBs detected by HETE/FREGATE, whose spectral parameters are given in Barraud et al. (2003). We compare them with the pseudo-redshifts of 8 GRBs with known redshifts in Table 1 in order to assess the role of the redshift in the non-detection of the afterglows for these GRBs. The pseudo-redshifts of these 18 GRBs are given in Table 2.

The first remark is that short/hard GRBs should probably not be integrated in our framework. GRB 020531 for instance has a low X value, which results in an unrealistically high pseudo-redshift. Having no redshift for short/hard bursts we cannot evaluate, and eventually calibrate, our redshift indicator for these bursts. The two shord/hard GRBs of our sample, GRB 020113 and GRB 020531, are thus excluded from the rest of our analysis.

  
Table 2: Observed properties of 18 GRBs with no measured redshift. The eight columns give the name of the GRB, the time of arrival, the duration T90 in seconds, the three spectral parameters ($\alpha $, $\beta $, and E0), the gamma-ray fluence  $S_{\gamma }$ in units of 10-6 erg cm-2 in the range 30-400 keV, the pseudo-redshift $\hat{z}$, and a comment on the eventual detection of an afterglow (XRR stands for X-Ray Rich GRB, OA, XA, and RA, respectively for Optical Afterglow, X-ray Afterglow, and Radio Afterglow).
\begin{table}\begin{displaymath}
\begin{array}{llllllll}
\hline\hline
\noalign...
...rm XRR} \\
\noalign{\smallskip }
\hline
\end{array}\end{displaymath}\end{table}

The median pseudo-redshift of long GRBs in Table 2 is 1.65, while it is only 0.88 for the 8 FREGATE GRBs with a measured redshift in Table 1. If we believe the correlation between the pseudo-redshifts and the true redshifts, this indicates that the redshift certainly plays a role in the non-detection of the afterglows of FREGATE GRBs, even if this is not the only factor as emphasized by Crew et al. (2003).

While pseudo-redshifts can be useful for statistical analyses, the information they convey is probably not meaningful for individual GRBs. We believe however that pseudo-redshifts could become a useful tool to quickly identify high redshift GRBs from the gamma-ray data alone. A first step in this direction is obviously to prove the validity of pseudo-redshifts for this task. GRB 020127 may appear as a good test case in this context because it has a high pseudo-redshift ( $\hat{z} = 3.9$ in Table 2), a possible X-ray afterglow, a possible radio afterglow, and a candidate host galaxy (Fox et al. 2002a, 2002b, 2002c). If the host candidate is at a redshift of about 3, this would strengthen the validity of pseudo-redshifts as a tool for the quick identification of high redshift GRBs.

   
4 Conclusion

We propose an empirical redshift indicator for GRBs, which is easily computed from the gamma-ray data alone and provides "pseudo-redshifts'' accurate to a factor of two. Despite their moderate accuracy, we believe that their easy computation will make these pseudo-redshifts useful in future GRB studies. Their possible applications include a statistical comparison of the distance distribution of distinct GRB populations, constraints on the star formation rate at high redshifts, and the fast identification of remote GRBs, with redshifts beyond three.

The usefulness of these pseudo-redshifts will ultimately depend on the confirmation of their accuracy, which will be tested as a larger number of GRBs with known redshifts become available.

Acknowledgements

The author is grateful to R. Mochkovitch for insightful discussions on pseudo-redshifts and to C. Barraud for helpful suggestions. JLA acknowledges the use of J. Greiner GRB page at http://www.mpe.mpg.de/~jcg/grbgen.html, and useful comments from the referees.

References



Copyright ESO 2003