text:
# On the White Dwarf distances to Galactic Globular Clusters

## 1 Introduction

Globular cluster stars are possibly the most suitable objects to be used in order to constrain the age of the universe. Since stars in a given globular cluster are coeval and share the same initial chemical composition, the turn–off brightness derived from the globular cluster colour-magnitude diagram (CMD) provides straightforwardly the cluster age: if the cluster distance is known, one only has to compare the observed cluster turn–off brightness with the correspondent quantity predicted by theoretical isochrones with the appropriate chemical composition in order to derive a firm estimate of its age.

The globular cluster distance scale is however still the subject of an intense debate. Main Sequence (MS) fitting distances using metal poor Hipparcos subdwarfs with well determined parallaxes provide long distances implying globular cluster ages of the order of 12 Gyr (see, e.g., Gratton et al. 1997). The same long distances are derived by Salaris & Weiss (1997, 1998), Mazzitelli, D’Antona & Caloi (1995), Cassisi et al. (1999) when using as standard candles theoretical Zero Age Horizontal Branch (ZAHB) models, while shorter distances are obtained from the ZAHB models of Vandenberg et al. (2000). On the other hand, the calibration of the RR Lyrae stars absolute brightness using the statistical parallax methods (see, for instance, Fernley et al. 1998 and Luri et al. 1998) provides consistently much shorter distances and therefore higher globular cluster ages. Taking these results at face value, the uncertainty of the globular cluster distance scale appears to be still of the order of 0.25–0.30 mag, which translates into an indetermination on the globular cluster age of about 20% — see the discussion in Renzini et al. (1996).

An alternative empirical method to derive globular cluster distances has been recently applied to the galactic globular cluster NGC 6752 (Renzini et al. 1996). The distance indicator is in this case a template local sequence of white dwarf stars with effective temperature ($`T_{\mathrm{eff}}`$) ranging between 10 000 and 20 000 K, and precise parallax measurements, which is fitted to the dereddened cluster white dwarf sequence. The vertical shift applied to the local template sequence in order to fit the sequence of the cluster provides its distance modulus. We shall refer to this technique as the WD–fitting technique, because of its close analogy with the MS–fitting technique.

A key assumption of this method is that the white dwarfs of the local template sequence — which, in the case of the paper by Renzini et al. (1996), have an estimated average mass $`M=0.515M_{}`$ — are totally equivalent to the white dwarfs of the cluster. If this is the case, the main advantage of the WD–fitting technique over the MS–fitting one is that it is in principle independent of the knowledge of the initial chemical composition of the globular cluster, since all white dwarfs have virtually metal free atmospheres, and therefore it avoids the uncertainties introduced in the MS–fitting method by the colour corrections that have to be necessarily applied to the local subdwarfs in order to precisely match the metallicity of the globular cluster under scrutiny.

Another potential advantage of the WD–fitting method with respect to the MS–fitting technique is, according to Renzini et al. (1996), that local white dwarfs are more abundant than metal poor subdwarfs, and therefore it is possible in principle to have a larger sample of local calibrators. However, since white dwarfs are much dimmer than the subdwarfs used in the MS–fitting technique, up to now it is possible to apply the WD–fitting technique only through HST observations of relatively close globular clusters.

Finally it should be mentioned as well that in the $`T_{\mathrm{eff}}`$ range considered by Renzini et al. (1996), local white dwarfs appear in two types, either as the so–called DA spectral type, which are characterized by an envelope made up of pure H (on top of a He layer), or as the non–DA type, which is characterized by an almost pure He envelope and possibly, at least for 20% of them, traces of H with abundances of the order of $`10^4`$ by number, or even less — see the review by Koester & Chanmugam (1990), and references therein. For local white dwarfs the number ratio of DA versus non–DA is of the order of 4:1 in this temperature range. Thus, the heterogeneity of the local sample of white dwarfs could potentially have undesired effects on the determinations of the ages and distances of globular clusters using the WD–fitting method.

Because of the potential advantages offered by the WD–fitting method, we analyze in detail the possible systematic uncertainties involved in this technique due to the poor knowledge of parameters affecting the CMD location of the template and globular cluster white dwarfs, expanding upon the discussion by Renzini et al. (1996). We pay particular attention to the role played by differences in the thickness of the surface H and He layers on the final distance determination — which has not yet been explored —, to the possible differences between white dwarfs in clusters with predominantly blue or red Horizontal Branches (HB), and also to the differences obtained when using DA or non–DA model atmospheres, or evolutionary models of white dwarfs with He cores. All along our analysis we use evolutionary white dwarf models in a purely differential way.

The paper is organized as follows. In Sect. 2 we briefly describe the theoretical models, while in Sect. 3 the influence of the various parameters affecting the CMD of white dwarfs is assessed. An exhaustive estimate of the systematic errors involved in the WD–fitting technique is derived in Sect. 4, which is followed by Sect. 5 were we draw our conclusions.

## 2 The models

The stellar evolution computations have been performed using the same code and the same input physics described in Salaris et al. (2000). We just recall here that in the temperature range we are dealing with ($`T_{\mathrm{eff}}`$ well above 6 000 K) the OPAL radiative opacities (Iglesias & Rogers 1993) with $`Z=0`$ are used, together with the conductive opacities by Itoh et al. (1983) supplemented by the Hubbard & Lampe (1969) ones; the boundary conditions for the integration of the stellar structure have been derived in this case from a grey $`T(\tau )`$ relationship, which is completely adequate at these temperatures — see, for instance, Hansen (1999) and Salaris et al. (2000). The equation of state for the H and He envelopes is that of Saumon, Chabrier & Van Horn (1995), while the equation of state for the carbon-oxygen (CO) core is from Segretain et al. (1994); neutrino energy losses have been taken from Itoh et al. (1996). Bolometric luminosities and effective temperatures have been transformed into $`V`$ magnitudes and colours using the relationships by Bergeron, Wesemael & Beauchamp (1995).

We have computed white dwarf model sequences for several masses, CO stratifications — unless stated otherwise our reference CO stratification is that from Salaris et al. (1997) — and various H and/or He envelope thicknesses, which will be described in the following section. For all the computations an initial model with the selected CO profile and envelope structure was converged at $`\mathrm{log}(L/L_{})2.0`$ and evolved down to sufficiently low temperatures.

## 3 The white dwarf sequence location in the CMD

As already mentioned in §1, the WD–fitting method is based on the fitting of a template white dwarf sequence to the correspondent globular cluster one. The key ingredient to derive reliable distances (and hence ages) is to ensure that the local white dwarfs included in the template sequence are homogeneous with the globular cluster ones.

All stars in a globular cluster are coeval; therefore, the luminosity (and, thus, $`T_{\mathrm{eff}}`$) of a white dwarf of a given mass is constrained by the fact that the sum of its cooling age ($`t_{\mathrm{cool}}`$) plus the evolutionary time of its progenitor ($`t_{\mathrm{prog}}`$) — which, to a good approximation, is equal to its main sequence lifetime, since subsequent evolutionary phases are much shorter — must be equal to the globular cluster age ($`t_{\mathrm{GC}}`$). It is immediately clear that one also has to assume an initial–final mass relationship, that is, a relation between the white dwarf mass and the initial mass of its progenitor. In Fig. 1 we show, as an example, the location of the white dwarfs of an old cluster of initial solar metallicity and ages of 10 and 12 Gyr, respectively, in the $`M_\mathrm{V}/(BV)`$ plane, as predicted by the theoretical models (Salaris et al. 2000). The horizontal lines mark the region where white dwarfs have effective temperatures between 10 000 and 20 000 K, which is the temperature range of the template field white dwarfs used by Renzini et al. (1996) to derive the distance to NGC 6752.

The initial–final mass relationship has been derived considering the CO core masses after the first thermal pulse from the evolutionary models of Salaris et al. (1997). This initial–final mass relationship provides white dwarf masses nearly constant and equal to 0.54 $`M_{}`$ for progenitor masses up to $`\mathrm{\hspace{0.17em}2.5}M_{}`$, and then increasing up to 1.0 $`M_{}`$ when the progenitor mass reaches $`7.0M_{}`$. It can be easily seen in Fig. 1 that the cluster white dwarf sequences show a pronounced turn to the blue at their dimmer end which, as the age of the cluster increases, is located at increasingly larger magnitudes. Until this blue–turn the white dwarf sequence is almost coincident with the cooling track of the $`\mathrm{\hspace{0.17em}0.54}M_{}`$, while the blue–turn is due to the contribution of more massive white dwarfs. It is easy to understand this behaviour by recalling that at each brightness along the cooling sequence the constraint $`t_{\mathrm{GC}}=t_{\mathrm{cool}}+t_{\mathrm{prog}}`$ has to be valid. Since $`t_{\mathrm{cool}}`$ is very short at the bright end of the cooling sequence, and practically negligible with respect to $`t_{\mathrm{GC}}`$, $`t_{\mathrm{prog}}`$ and, hence, the progenitor mass has to be, to a very good approximation, constant for a large magnitude range, and very close to the turn–off mass. On the contrary, towards the dim end of the cooling sequence $`t_{\mathrm{cool}}`$ becomes a sizeable fraction of $`t_{\mathrm{GC}}`$ and, thus, the contribution of the white dwarfs coming from higher mass progenitors (and, consequently, with smaller $`t_{\mathrm{prog}}`$) is quite apparent. A similar feature is expected for the halo white dwarf population (Isern et al. 1998).

This feature can also be seen in Fig. 2 (upper panel) where we show the run of the initial masses as a function of the $`M_\mathrm{V}`$ magnitude along the cluster white dwarf sequences shown in Fig. 1. In the $`T_{\mathrm{eff}}`$ range we are dealing with (marked, as in Fig. 1, by vertical dashed lines), the initial mass is almost constant — as well as, due to our selected initial–final mass relationship, the actual white dwarf mass shown in the lower panel. It is remarkable that even two magnitudes below the lower limit of this effective temperature range the initial mass of white dwarf progenitors is still constant. This result (the constancy of the progenitor mass of white dwarfs contributing to the effective temperature range, not its actual value) is completely general, that is, independent of the adopted initial–final mass relationship, since for all possible white dwarf masses the evolutionary times are very fast at these high temperatures; moreover, this result is also independent of the initial metallicity of the white dwarf progenitors and therefore valid also in the case of white dwarf progenitors with lower metallicities, typical of galactic globular clusters.

Since the progenitor mass is constant, one can then reasonably expect that, even in the case of a completely different initial–final mass relationship, the mean white dwarf mass and the spread around this value is constant, and equal to the values attained at the beginning of the cooling sequence. This behaviour is very different from what happens in field disk white dwarfs. In this case, because of the ongoing star formation processes, very different white dwarf masses can populate the same $`T_{\mathrm{eff}}`$ range, depending on the age of their progenitors.

There are in principle four quantities which can introduce a systematic error in the white dwarf distance determination to globular clusters, provided that there are significant differences between the properties of the local template white dwarf sequence and the globular cluster one, namely:

1. The mass of the white dwarf cooling sequence.
2. The chemical composition of the envelope (DA, non–DA).
3. The thickness (in mass) of the envelope.
4. The chemical stratification of the core (different CO profiles or, even, He–core white dwarfs).

In the following we are going to discuss the influence of these four parameters on the CMD location of white dwarf cooling sequences in the $`M_\mathrm{V}/(BV)`$ (hereinafter $`BV`$) and $`M_\mathrm{V}/(VI)`$ (hereinafter $`VI`$) planes, which are actually the most popular CMDs employed for studying cluster white dwarfs — see, e.g., Renzini et al. (1996), Richer et al. (1997), and Von Hippel et al. (2000). We consider as a reference the $`T_{\mathrm{eff}}`$ range between 10 000 and 20 000 K, already selected by Renzini et al. (1996) for their distance determination. In this temperature range the white dwarf sequence is sufficiently populated and bright enough to be detected with a reasonable photometric error in close globular clusters, using the HST or giant ground based telescopes.

### 3.1 The mass of the white dwarf sequence

Stars in a globular cluster lose their mass along the Red Giant Branch (RGB) due to stellar winds. The amount of mass lost during this phase is not the same for each red giant star, and this, in turn, is the origin of the extended horizontal branches (HB) observed in Galactic globular clusters; this also means that the mass evolving along the RGB (which at the beginning of the red giant phase is almost equal to the mass of the star at the cluster turn–off) can give birth to HB stars with different values of the total mass, but with the same initial He core mass. These different HB stars will all end up as white dwarfs with a mass range which is constrained mainly by the size of the helium core at the beginning of the He-burning phase and by the mass loss processes along the subsequent Asymptotic Giant Branch (AGB) phase.

Renzini et al. (1996) summarize the results about the few semiempirical determinations of white dwarf masses in galactic globular clusters, providing a very narrow range of values, namely $`0.53\pm \mathrm{\hspace{0.17em}0.02}M_{}`$. This range of values nicely overlaps with current determinations of the CO core mass at the first thermal pulse derived from stellar evolutionary models — see, e.g., Wagenhuber (1996) and Salaris et al. (1997) — for initial masses of about $`0.81.0M_{}`$. Moreover, it is also in agreement with determinations of the mass distribution of nearby field white dwarfs, which peaks at about $`0.55M_{}`$ (Bragaglia et al. 1995; Reid 1996). However, regarding this last result, the reader should keep also in mind the recent results of Bergeron, Leggett & Ruiz (2000), who find a mean mass of field white dwarfs about $`0.10M_{}`$ higher at lower $`T_{\mathrm{eff}}`$ values.

Very recently however, Alves, Bond & Livio (2000) derived semiempirically the mass of the central star in a planetary nebula of the globular cluster M 15, obtaining a value $`0.60\pm \mathrm{\hspace{0.17em}0.02}M_{}`$. It is not clear yet if this result is an indication of a globular cluster initial–final mass relationship with a larger dispersion than that derived by Renzini et al. (1996), or it is due (as discussed by the authors) to an increment of the mass through an interaction in a close binary system. We just recall here that there is also another possible indication of a large dispersion for the initial–final mass relationship in the results obtained by Reid (1996) when considering the white dwarf population of the Hyades open cluster.

Another important piece of evidence must also be taken into account. The determination of the mass range $`0.53\pm \mathrm{\hspace{0.17em}0.02}M_{}`$ — and the results of Alves et al. (2000), as well — is mainly based on the maximum brightness of AGB stars and on the luminosity of post-AGB stars, which can be directly related to the mass of the degenerate CO core prior to the beginning of the cooling phase. However, there are several globular clusters with an horizontal branch morphology showing the presence of an extended blue tail, as in the case of NGC 6752. Regardless of the physical mechanism(s) producing the stars which populate the hot side of the HB in globular clusters, it is well known that they are stars which have lost a large amount of their envelope during the RGB phase. At the end of the He–burning phase these stars, depending on the mass of their residual H–rich envelope, can behave as post–Early AGB structures or as AGB-manqué ones — see, for instance, Greggio & Renzini (1990) — and do not experience the thermally pulsing phase on the AGB as most massive HB stars do. On theoretical grounds, the minimum (initial) HB mass which does not experience the AGB thermal pulses is a quite robust prediction (Dorman et al. 1993, Bono et al. 1997), being equal to $`\mathrm{\hspace{0.17em}0.52}M_{}`$. Since the mass of the CO core at the end of the HB phase is of the order of $`0.45M_{}`$, globular clusters with extended blue tails can produce also white dwarfs with masses in the range $`0.450.52M_{}`$, on average $`\mathrm{\hspace{0.17em}0.05}M_{}`$ less massive than the white dwarf progeny of stars climbing up the AGB. Thus, it is worth considering the possibility of different masses for the globular cluster cooling sequences.

Turning now our attention to the CMD location of white dwarfs of different masses, it is well known from relatively simple physical considerations that the larger the white dwarf mass is, the smaller its radius is. In Fig. 3 we display in the $`BV`$ and $`VI`$ planes, for effective temperatures between 10 000 and 20 000 K, two representative white dwarf cooling tracks with masses equal to 0.54 and $`0.61M_{}`$. The higher mass is shifted to higher $`M_\mathrm{V}`$ values. The derivative $`\mathrm{\Delta }M_\mathrm{V}/\mathrm{\Delta }(M/M_{})`$ equals to 2.3 for masses between 0.45 and 0.60 $`M_{}`$.

### 3.2 The chemical composition of the envelope

As already mentioned in the introduction, field white dwarfs with $`T_{\mathrm{eff}}`$ between $``$ 20 000 and 10 000 K can be present either in the DA or non–DA types, with a number ratio of DA to non–DA equal to 4. Our selected standard envelope thicknesses (see next subsection) are $`\mathrm{log}q(\mathrm{H})=4.0`$ and $`\mathrm{log}q(\mathrm{He})=2.0`$ for DA white dwarfs, and $`\mathrm{log}q(\mathrm{He})=3.5`$ for the non–DA spectral type, where $`q`$(H) and $`q`$(He) indicate the ratio of the mass contained in the H and He envelope layers to the total white dwarf mass, respectively.

The envelope chemical composition strongly affects not only the white dwarf cooling times, but also their location in the CMD, and the precise shape of the cooling track, as clearly shown in Fig. 4 where we show the cooling tracks for a $`0.54M_{}`$ DA white dwarf (solid line) and a non–DA white dwarf of the same mass. It is evident that in the $`BV`$ plane DA white dwarfs are brighter at a fixed colour. This is clearly confirmed by observations of local DA white dwarfs, see e.g. panel 2 of Fig. 1 in Renzini et al. (1996), while the reverse happens in the $`VI`$ plane. Moreover, the separation between DA and non–DA white dwarfs is larger in the $`BV`$ plane, making possibly more evident their identification in a given globular cluster. $`M_\mathrm{V}`$ differences of up to $``$ 1 mag at a fixed $`BV`$ value are possible in the $`BV`$ plane, while they are a factor $``$ 5 smaller in the $`VI`$ plane. Also the shape of the DA and non–DA cooling sequences is different, especially in the $`BV`$ plane.

So far, we have considered DA and non–DA white dwarfs having their very outer layers made of either pure hydrogen or pure helium. However, observations suggest that in this temperature range DA white dwarfs can have a small amount of helium in their spectra, of the order of $`10^5`$ by number, and that 20% of non–DA white dwarfs show some hydrogen in their envelopes, with abundances of the order of $`10^4`$ by number at most. We have tested the effect of these small admixtures of helium in the outer envelope of DA white dwarfs, and of hydrogen in the envelope of non–DA white dwarfs, by computing appropriate models for a $`0.61M_{}`$ white dwarf, and using the results of Bergeron, Saumon & Wesemael (1995) about the influence of H–He mixtures on the derived white dwarf colours. We have found negligible variations of the CMD location with respect to the case of our standard He–free DA white dwarf envelopes and H–free non–DA atmospheres.

### 3.3 The thickness of the envelope

The thicknesses of the hydrogen and helium envelopes of white dwarfs has been the subject of many investigations during the last decade. The thickness of these layers is a key ingredient to determine the evolution of the white dwarf since, due to their opacity, they basically regulate the energy loss rate of the isothermal, highly conductive electron degenerate core. Moreover, the envelope thickness also affects the radius of the white dwarf at a given effective temperature. As discussed by D’Antona & Mazzitelli (1990), due to the unknown details of the mass loss process during the AGB phase and the planetary nebula ejection, theoretical evolutionary models cannot yet provide too stringent predictions about the thickness of the hydrogen and helium layers surrounding the degenerate CO core.

Observational constraints based on spectroscopic (Barstow et al. 1993) as well as on asteroseismological (Clemens 1995) analyses suggest that local field DA white dwarfs have typically hydrogen layers of about $`\mathrm{log}q(\mathrm{H})=4.0`$. On the other hand, the study of the mass–radius relationship for a sample of field white dwarfs with known Hipparcos parallaxes provides indications that the thickness of the hydrogen layers spans a range of values within $`\mathrm{log}q(\mathrm{H})4.0`$ and $`\mathrm{log}q(\mathrm{H})7.0`$ (Provencal et al. 1998). As for the thickness of the helium layer below the hydrogen envelope the assumed reference value comes basically from stellar evolution constraints, and it is $`\mathrm{log}q(\mathrm{He})2.0`$ (Hansen 1999 and references therein). We have verified, by computing white dwarf cooling tracks with $`\mathrm{log}q(\mathrm{He})`$ decreased by 1 dex (keeping the total white dwarf mass constant) that the location in the CMD of DA white dwarfs is basically unchanged. As for non–DA white dwarfs, estimates of the envelope thickness range between $`10^4M_{}`$ and $`10^2M_{}`$ (Pelletier et al. 1986; MacDonald, Hernanz & José 1998).

All these results are however for field white dwarfs; no indications exist yet about $`\mathrm{log}q(\mathrm{H})`$ and $`\mathrm{log}q(\mathrm{He})`$ in globular clusters, apart from the fact that theoretical evolutionary models of white dwarfs coming from blue HB progenitors — in the hypothesis that mass loss during the He-burning phases is negligible — predict an upper limit to the thickness of the very outer H layers ranging from $`\mathrm{log}q(\mathrm{H})3.5`$ to $`4.0`$ (Castellani et al. 1994a).

In Fig. 5 we show the effect of varying the thickness of the external H and He layers for, respectively, a DA and and a non–DA white dwarf of $`0.54M_{}`$. In the case of non–DA white dwarfs our reference value is $`\mathrm{log}q(\mathrm{He})=3.5`$, which we have changed by $`\pm `$1.0 dex, obtaining virtually no variation in the location of the cooling track (the three tracks perfectly overlap in Fig. 5). In the case of H envelopes the situation is different, since the thickness of the H layers affects appreciably the location of the track. In Fig. 5 representative models with $`\mathrm{log}q(\mathrm{H})=4.0`$ and $`\mathrm{log}q(\mathrm{H})=6.0`$ are displayed; a reduction of the thickness of the H envelope shifts the track — at a constant value of the colour — towards higher $`M_\mathrm{V}`$ values, with a derivative $`\mathrm{\Delta }M_\mathrm{V}/\mathrm{\Delta }\mathrm{log}q(\mathrm{H})0.035`$ for $`\mathrm{log}q(\mathrm{H})`$ ranging between $`4.0`$ and $`7.0`$.

### 3.4 The chemical stratification of the CO core

The chemical stratification of the CO core may potentially affect as well the CMD location of the white dwarf track since the the mass-radius relationship also depends on the electron mean molecular weight of the electron-degenerate core. To this regard, it is important to realize that the value of the CO ratio along the white dwarf core is subject to some uncertainties due to our poor knowledge of the value of the $`{}_{}{}^{12}\mathrm{C}(\alpha ,\gamma )^{16}\mathrm{O}`$ reaction rate — see the detailed discussions in Salaris et al. (1997) or Cassisi et al. (2001), and references therein. Moreover, for a given value of the $`{}_{}{}^{12}\mathrm{C}(\alpha ,\gamma )^{16}\mathrm{O}`$ reaction rate, the CO profile is possibly a function of the initial metallicity of the white dwarf progenitor (Umeda et al. 1999).

As a numerical test, and in order to mimic these effects, we have computed two cooling tracks for a $`0.54M_{}`$ white dwarf, where, for the first case, we have substituted our reference CO stratification shown in Fig. 6 with a flat profile with 50% of C and 50% of O by mass, whereas in the second case we have inverted the role of C and O, preserving the shape of the initial chemical profile. The resulting cooling tracks are however negligibly affected by these changes in the effective temperature range we are dealing with. This means that, regardless of the still large uncertainty affecting the $`{}_{}{}^{12}\mathrm{C}(\alpha ,\gamma )^{16}\mathrm{O}`$ reaction rate, and regardless of the possible effects on the CO profile of the initial progenitor metallicity, the CMD location of the bright portion of the white dwarf cooling sequence is not appreciably affected.

### 3.5 He–core white dwarfs

He–core white dwarfs are the byproduct of strong mass loss along the red giant phase (due for example to the interaction with a binary companion), which strips out the H envelope on top of the H burning shell before the degenerate He core reaches the critical mass to experience the He flash and the subsequent quiescent He burning phase (Kippenhahn, Kohl & Weigert 1967; Castellani, Luridiana & Romaniello 1994b; Hansen & Phinney 1998). Their mass can range from a value close to the degenerate He core–mass at the red giant tip (about $`0.50M_{}`$ for the most metal poor globular clusters as M 92, decreasing down to about $`0.48M_{}`$ for the metal rich ones as 47 Tuc) down to masses of the order of $`0.20M_{}`$. There is possibly at least one direct observational indication that He–core white dwarfs contribute to the cooling sequence of globular clusters (Moehler et al. 2000).

In Fig. 7 we show the CMD location of a $`0.45M_{}`$ He–core white dwarf with a H surface layer of $`\mathrm{log}q(\mathrm{H})=4.0`$. In the same plot our reference DA and non–DA model sequences of $`0.54M_{}`$ are also shown. The He–core white dwarf mass is close to the upper possible value, lower masses being shifted to the red side of the CMD. As expected purely on the base of the mass difference with respect to the DA white dwarf with a CO core plotted in the figure, the He–core white dwarf is shifted to higher brightnesses at a fixed colour. This feature is reminiscent of the two bright white dwarfs discarded by Renzini et al. (1996) in their fitting procedure — compare the $`BV`$ panel in Fig. 7 with Fig. 1 of Renzini et al. (1996) — because they were clearly located to the right side of the main DA cluster sequence. In the $`VI`$ plane the He–core white dwarf sequence overlaps with the non–DA one at the bright end of the $`T_{\mathrm{eff}}`$ interval. In the $`BV`$ plane, due to the steeper slope of the H–envelope cooling sequences, the He–core white dwarf is on average closer in colour to the CO one than in the $`VI`$ plane.

## 4 Systematic errors on the WD fitting distances

In this section we quantify the possible systematic errors on the distances derived from the WD–fitting technique taking into account the results previously discussed. We study the cases in which, due to our lack of precise empirical or theoretical determinations, the parameters which determine the position of the white dwarf sequence on the CMD are possibly different in the local template and in the globular cluster sequence, but within the range allowed by independent observations. Because the ratio between the number of DA and non–DA white dwarfs in the field is 4:1 (and presumably this ratio is the same in globular clusters) one is forced to use the cooling sequences of DA white dwarfs in order to build the template sequence. This, in turn, implies that an important role is played by the thickness of the H surface layers. To this date, we do not have any empirical indication about the value of $`q`$(H) for white dwarfs in globular clusters, and about its dependence on the HB morphology. In principle, there could be sizeable systematic differences due, for example, to the fact that in globular clusters with blue HB (like M 92 and NGC 6752) the main component of the white dwarf population could be produced by stars not experiencing the AGB phase, in contrast with white dwarfs populating red HB globular clusters (like 47 Tuc). We are forced, therefore, to treat the thickness of the envelope layers as a free parameter, in the hypothesis that the possible range of values is within the range spanned by field white dwarfs.

### 4.1 Uncertainties on the template sequence

Since the mass range for white dwarfs in globular clusters appears to be mainly in the range $`0.53\pm \mathrm{\hspace{0.17em}0.02}M_{}`$, it is clear that the template white dwarf sequence should be made by stars within this mass range. However, there are almost no field white dwarfs suitable for the WD–fitting method, for which the value of the mass is derived empirically (like, for instance, the case of white dwarfs in close binary systems, as Sirius B). In order to build the template DA white dwarf sequence, one needs therefore to employ mass determinations obtained from semiempirical methods, as discussed in Bragaglia et al. (1995). It is possible to determine spectroscopical values for the gravity ($`g`$) and $`T_{\mathrm{eff}}`$ of a sample of local white dwarfs, and then derive the masses employing theoretical $`\mathrm{log}gT_{\mathrm{eff}}`$ relationships. However, the assumed thicknesses of the envelope layers of DA white dwarfs affect the $`\mathrm{log}gT_{\mathrm{eff}}`$ relationship. In fact, it can be shown that a uniform distribution of $`\mathrm{log}q(\mathrm{H})`$ between $`4.0`$ and $`7.0`$ produces a mass range dispersion of about $`0.04M_{}`$ for a given couple of $`\mathrm{log}g`$ and $`T_{\mathrm{eff}}`$ values.

One could speculate if, after deriving a certain mass range for the template white dwarfs assuming a given spread in the $`\mathrm{log}q`$(H) values, it is possible to judge if the results and the assumptions about the envelope thickness are consistent with the distribution of the stars in the CMD. As we are going to show now, this is not an easy task when considering realistic errors on the parallax of the objects. We determined, using our cooling tracks, three template sequences of about 100 DA white dwarfs by means of a simplified version of our Monte–Carlo simulator (García–Berro et al. 1999), assuming a negligible error on their colours, and allowing for an error of $`\pm \mathrm{\hspace{0.17em}0.10}`$ mag in $`M_\mathrm{V}`$, which is approximately the average error of the brightness of the template white dwarfs of Renzini et al. (1996) due to the error on their parallaxes (which is of the order of 5%). The first sequence (sequence A) comprises masses in the range between 0.51 and $`0.55M_{}`$, and $`\mathrm{log}q`$(H) uniformly distributed between $`4`$ and $`7`$; sequence B has the same mass range but the value of the thickness of the hydrogen outer layer was kept constant at $`\mathrm{log}q(\mathrm{H})=4`$; finally, sequence C is characterized by $`M=0.55M_{}`$ and $`\mathrm{log}q(\mathrm{H})=4`$. In Fig. 8 we show the three sequences in the $`BV`$ CMD (for $`T_{\mathrm{eff}}`$ between 10 000 and 20 000 K), after shifting vertically sequence B by +0.04 mag, and sequence C by +0.015 mag, in order to reproduce the average brightness of sequence A. The average brightness of the three sequences is different, but the dispersion around the average looks very much the same, and dominated by the parallax error, in spite of the fact that the mass and/or envelope thickness ranges are different.

The main result of this exercise is that uncertainties in the value of $`\mathrm{log}q`$(H) produce unavoidable uncertainties on the precise value of the masses of the template white dwarfs and, thus, on the precise location of the template sequence, which may contribute to the error budget of the WD–fitting technique by amounts of the order of less than 0.05 mag. For the sake of conciseness, in the following we will consider as our reference template a sequence (determined from the Monte–Carlo simulations) made of about 100 DA white dwarfs with masses in the range between 0.51 and $`0.55M_{}`$, $`\mathrm{log}q`$(H) uniformly distributed between $``$4 and $``$7 and effective temperatures in the range between 10 000 and 20 000 K.

### 4.2 Uncertainties on the white dwarf mass and envelope thickness. Clusters with a red HB

We define here as red HB globular clusters those clusters in which the mass of the HB stars is always larger than $`\mathrm{\hspace{0.17em}0.52}M_{}`$, corresponding to a colour of HB stars larger than $`(BV)`$0.21. This means that all the white dwarfs populating red HB clusters are a product of AGB evolution (see the discussion in our previous section). In this case, empirical estimates of the white dwarf masses show a range between 0.51 and $`0.55M_{}`$. Thus, if the unknown thickness of the H layers of the DA cluster white dwarfs is similar to that of the local white dwarfs, in principle there should be no systematic error in the distances derived using our reference template sequence.

To obtain a reasonable estimate of the maximum possible systematic error, our reference sequence has been fitted to a cluster DA sequence with, respectively, $`M=0.51M_{}`$ and $`\mathrm{log}q(\mathrm{H})=4`$, and $`M=0.55M_{}`$ and $`\mathrm{log}q(\mathrm{H})=7`$. The systematic error in the distance modulus is of $`0.10`$ mag in the first case and of $`+0.10`$ mag in the second case, with a $`\mathrm{\hspace{0.17em}60}`$% contribution due to the effect of the envelope thickness. It is obvious that in the case that the cluster white dwarfs have an intrinsic spread in mass and/or $`\mathrm{log}q(\mathrm{H})`$ within the mentioned range, the systematic error on their distance modulus is smaller than $`\pm \mathrm{\hspace{0.17em}0.10}`$ mag.

If one accepts the results of Alves et al. (2000) about M 15 as due to a genuine spread in the initial–final mass relationship, the mass range of globular cluster white dwarfs could possibly extend up to $`0.60M_{}`$. By repeating the previous exercise considering as an extreme case $`M=0.60M_{}`$ and $`\mathrm{log}q(\mathrm{H})=7`$ for the cluster white dwarfs, one gets systematic errors of +0.20 mag in the cluster distance modulus. Of course this is a very extreme case. Probably, it si more realistic to consider a mass distribution between 0.51 and $`0.60M_{}`$ and $`\mathrm{log}q(\mathrm{H})`$ distributed between $`4`$ and $`7`$. With this sample of 100 cluster white dwarfs and using the same reference template white dwarf sequence, we determine a systematic error on the distance modulus of +0.05 mag. Should all the globular cluster white dwarfs have the same H layers thickness, the systematic errors would be then $`+0.01`$ mag if $`\mathrm{log}q(\mathrm{H})=4`$ and $`+0.11`$ mag if $`\mathrm{log}q(\mathrm{H})=7`$.

### 4.3 Uncertainties on the white dwarf mass and envelope thickness. Clusters with a blue HB

We define as blue HB globular clusters those clusters in which the mass of the HB stars extends below $`\mathrm{\hspace{0.17em}0.52}M_{}`$. In general, even globular clusters with long blue tails always have a fraction of their HB population located at colours larger than $`(BV)0.21`$, and, therefore, white dwarfs populating clusters with a blue HB are the product of both AGB evolution and AGB-manqué or post-Early AGB stars. As discussed before, one expects that this kind of globular clusters is populated by white dwarfs with masses down to $`\mathrm{\hspace{0.17em}0.45}M_{}`$. It appears that blue HB clusters should therefore show a larger mass spread for the white dwarf population, with white dwarf masses ranging between approximately 0.45 and $`0.55M_{}`$.

By fitting our template sequence to a cluster sequence with a mass range between 0.45 and $`0.55M_{}`$ and a range of $`\mathrm{log}q`$(H) as in the template one, we obtained a systematic error on the distance modulus of $``$0.06 mag. If we accept a larger mass range of up to $`0.60M_{}`$, the systematic error goes down to $`0.04`$ mag, since the higher mass objects tend to compensate for the presence of masses lower than the template ones. This is of course just an estimate of the average systematic error in case of globular clusters with a blue HB. To obtain more precise evaluations one should compute synthetic CMDs for each given cluster in order to determine, from the observed HB morphology, the mass distribution along the HB and therefore, on the base of the discussion in Sect. 3.1, infer the possible distribution of white dwarf masses.

Considering now the possibility of different envelope thicknesses between template and cluster white dwarfs, if the white dwarf mass range is between 0.45 and $`0.55M_{}`$ and $`\mathrm{log}q(\mathrm{H})=4`$ in the cluster, the systematic error on the derived distance modulus is of about $`0.12`$ mag. If we adopt the same mass distribution but we change $`\mathrm{log}q`$(H) to $`7`$, the error is negligible, since the different envelope thickness compensates for the different mass range between the template and cluster sequence.

### 4.4 Contamination of the DA sample by non–DA white dwarfs

Up to know, we have made the assumption that cluster DA white dwarfs can be distinguished from non–DA ones. In absence of spectroscopical identification, the only way to discriminate between DA and non–DA white dwarfs is by checking the relative CMD location of the white dwarf sample (see Fig. 4). Here we estimate the maximum photometric error allowing for a clear distinction between DA and non–DA objects in the $`BV`$ and $`VI`$ CMDs. Moreover, we estimate the systematic errors introduced in the globular cluster distances when using a DA template sequence fitted to a cluster white dwarf sample made of both DA and non–DA stars. We assume that the cluster DA/non–DA number ratio in the $`T_{\mathrm{eff}}`$ interval we are dealing with is 4:1 as in the field. Notice, however, that below $`T_{\mathrm{eff}}=\mathrm{10\hspace{0.17em}000}`$ K the number of non–DAs is similar or even larger than the number of DAs. The most reasonable explanation for this is that the outer convective region mixes the H layer into the He envelope (Bergeron et al. 2000).

Figs. 9 and 10 show a sequence of about 80 DA and 20 non–DA white dwarfs ($`M=0.54M_{}`$, $`\mathrm{log}q(\mathrm{H})=4`$ for the DAs and $`\mathrm{log}q(\mathrm{He})=3.5`$ for the non–DAs) computed including 1$`\sigma `$ photometric errors in $`B`$, $`V`$, and $`I`$ of, respectively, 0.02 and 0.05 mag. Even in the case of the smallest errors, the non–DA and DA sequences are not well separated in the $`VI`$ plane, while the DA sequence is clearly distinguishable from the non–DA one in the $`BV`$ plane. With an error bar of the order of 0.05 mag, non–DA white dwarfs start to be mixed up with DA ones also in the $`BV`$ plane. In order to be more precise, when fitting a template DA sequence to a cluster white dwarf sequence including both DA and non–DA objects, one underestimates the globular cluster distance by only 0.03 mag if the fitting is performed in the $`VI`$ plane, whereas in the case in which the fitting is performed in the $`BV`$ plane the distance is overestimated by $`\mathrm{\hspace{0.17em}0.20}`$ mag. These figures are of course reduced if the ratio DA/non–DA is smaller than the value we have adopted here.

### 4.5 He-core white dwarfs

Another potential source of systematic errors in the WD–fitting distances is the presence of He–core white dwarfs (see Fig. 7), and although it is hard to give a precise quantitative assessment of this error source, we can however provide an argument to give some hints about the magnitude of this effect. We will show that the effect would be possibly small or even negligible. White dwarfs with cores made of He are the product of binary evolution. Actual estimates of the binary frequency in globular clusters provide a value of the order of 10% or less (Hut et al. 1992). If one accepts this estimate, we have verified that, by computing a cluster cooling sequence containing about 130 CO–core white dwarfs with DA/non-DA ratio 4:1 and constant envelope thickness ($`\mathrm{log}q(\mathrm{H})=4`$ for the DAs and $`\mathrm{log}q(\mathrm{He})=3.5`$ for the non–DAs), plus about 14 (about 10% of the total sample) He–core DA objects, even if all the He–core white dwarfs have masses of the order of $`0.45M_{}`$, their presence causes a systematic error of the order of only $``$0.01 mag. In the case that there is a spectrum of masses for He–core white dwarfs, ranging between $`\mathrm{\hspace{0.17em}0.20}`$ and $`0.45M_{}`$, then even less objects will be close enough to the CO cooling sequences for the relevant effective temperature range to be confused with them, while the others will lay at larger colours, clearly separated from the more populated CO sequences (see Fig. 11).

### 4.6 Cluster reddening

The knowledge of the reddening of the globular cluster is fundamental in order to derive the cluster distance using the WD–fitting technique. We tested the sensitivity of the derived distances to the uncertainty on the cluster reddening, by considering two DA cooling tracks with $`M=0.54M_{}`$ and $`\mathrm{log}q(\mathrm{H})=4`$, shifted in colour one with respect to the other; we then fitted the two sequences one of top of the other for different amounts of the relative shift, obtaining in the $`BV`$ and $`VI`$ planes a value for the derivative of the apparent distance modulus with respect to the $`(BV)`$ colour excess, $`\mathrm{\Delta }(mM)_\mathrm{V}/\mathrm{\Delta }E(BV)5.5`$.

## 5 Summary and conclusions

In this paper we have thoroughly discussed many possible sources of systematic errors on the globular clusters distances obtained using the WD–fitting technique. In order to do so we have first investigated in detail the changes produced in the $`BV`$ and $`VI`$ CMDs when varying the adopted mass of the template white dwarf cooling sequence, the influence of the assumed thicknesses of their envelopes on the location of the cooling tracks, as well as the consequences of changing the chemical composition of their envelopes and of their cores. We have then applied these results in order to derive, by means of Monte–Carlo simulations, a realistic estimate of the systematic errors involved in the application of this technique to galactic globular clusters. Our main results can be summarized as follows:

1. The unknown thickness of the H layers in cluster DA white dwarfs plays a non–negligible role, comparable to the role played by uncertainties on the white dwarf masses. For reasonable assumptions (derived from observations of field white dwarfs and constraints from globular cluster CMDs) about the unknown mass and $`\mathrm{log}q(\mathrm{H})`$ values in cluster DA white dwarfs, a realistic estimate of the maximum systematic error on the derived distance moduli is within $`\pm \mathrm{\hspace{0.17em}0.10}`$ mag. However, one should be aware that particular combinations of white dwarf masses and envelope thicknesses — still allowed by current weak or non-existent observational constraints — could produce larger errors.
2. A photometric precision better than $`\mathrm{\hspace{0.17em}0.05}`$ mag is needed in order to distinguish DA from non–DA white dwarfs in the $`BV`$ plane. An even better precision is needed when using the $`VI`$ plane. In the case of larger observational errors and no clear (spectroscopical) distinction between DA and non–DA cluster white dwarfs, fitting a template DA sequence to a cluster sequence made of a mixture of DA and non–DA stars introduces a very small systematic error $`0.03`$ mag in the $`VI`$ plane, but this error amounts to $`+0.20`$ mag, in the $`BV`$ plane. The $`VI`$ plane looks therefore better suited for the application of the WD–fitting method, since it permits to get rid of the uncertainty due to a possible contamination of the DA sample by non–DA white dwarfs.
3. Contamination by He–core white dwarfs should not influence appreciably the WD–fitting distances.
4. Due to the steep slopes of the white dwarf cooling curves in the CMD, the distance derived from the WD–fitting technique has a non negligible dependence on the adopted cluster reddening. We obtained a derivative $`\mathrm{\Delta }(mM)_\mathrm{V}/\mathrm{\Delta }E(BV)5.5`$, a dependence which is similar to the one of the MS–fitting technique in the $`BV`$ CMD.

###### Acknowledgements.

We thank M. Zoccali for useful comments to a preliminary version of the manuscript. Part of this work was supported by the Spanish DGES project numbers PB98–1183–C03–02, AYA2000-1785 and ESP98-1348, by the CIRIT and by Sun MicroSystems under the Academic Equipment Grant AEG–7824–990325–SP. One of us (S.C.) has been supported by MURST-Cofin2000- under the scientific project ”Stellar Observables of Cosmological Relevance”.


format: text

text:
# Broadband Observations of the Afterglow of GRB 000926: Observing the Effect of Inverse Compton Scattering

## 1 Introduction

Broadband observations of gamma-ray burst afterglows can in principle be used to constrain fundamental physical parameters of the explosion. In the fireball model, a relativistic blast wave expands into the surrounding medium, its hydrodymanical evolution being strongly influenced by the density structure of the medium as well as by the energy content and geometry (in particular collimation) of the ejecta. The temporal behavior of the afterglow emission which arises from the shocked gas depends on the shock evolution, and the partition of energy between the magnetic field and relativistic electrons, and can therefore probe these physical parameters given data of sufficient quality.

In this paper, we report the synthesized results from our multi-frequency followup campaign on the relatively bright GRB 000926. This campaign was aimed at studying the evolution of the afterglow to constrain the model parameters described above. Price et al. (2001) have reported the results from our multi-band ($`BVRI`$) optical monitoring. We combine these data with 4 epochs taken with the Hubble Space Telescope WFPC2, with Chandra X-ray Observatory (CXO) target of opportunity (TOO) observations, and with multi-frequency radio monitoring from the Very Large Array (VLA)111The NRAO is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. NRAO operated the VLA, the Ryle Telescope, and the Owens Valley Radio Observatory. We interpret the resulting broadband lightcurve in the context of a theoretical afterglow model.

## 2 Observations and Data Reduction

The Interplanetary Network discovered GRB 000926 on 2000 Sep 26.993 UT (Hurley et al., 2000). The afterglow of this 25 s long event was identified less than a day later (Gorosabel et al., 2000; Dall et al., 2000). The redshift, measured from optical absorption features, is $`2.0369\pm 0.0007`$ (Fynbo et al., 2000; Castro et al., 2000). The afterglow was well-monitored in the optical (Price et al., 2001; Fynbo et al., 2001), and was detected in the IR (DiPaola et al., 2000; Fynbo et al., 2001). Here we describe Hubble Space Telescope (HST), Chandra X-ray Observatory (CXO) and radio observations.

### 2.1 HST Observations

As part of an HST cycle 9 program we observed GRB 000926 at four epochs with the Wide Field Planetary Camera 2 (WFPC2), with the OT placed on WFPC CCD#3. In the first three epochs we observed at three passbands corresponding to the F450W, F606W, and F814W filters, and in the final epoch we used only the F606W and F814W filters. These observations took place between Oct 7.25 (10.26 days after the GRB) and Dec 16.9 (81.9 days after).

Table 1 shows a log of the HST observations, along with the magnitude derived for each filter for the 2-pixel radius region surrounding the OT. We determined the aperture correction using a 2 – 5 pixel radius, and we quote the corresponding 5 pixel radius magnitude. We calibrated the zeropoints, and converted the WFPC2 filters to Johnson Cousins magnitudes using the color transformations from Holtzman et al. (1995). We estimate the associated calibration uncertainty to be about 0.10 magnitude in $`B,V,`$ and $`R`$ bands, and 0.20 magnitudes in the $`I`$ band. To perform the transformations, we have interpolated the 3-filter WFPC2 data to obtain photometric points in $`B,V,R`$ and $`I`$. Measurements in only three of the bands are, therefore, truly independent. We tied the WFPC2 calibration to the ground-based data of Price et al. (2001) using their tertiary star at position RA $`12^h04^m11^s.75`$ DEC $`+51^{}46^{}56^{\prime \prime }.1`$ (J2000). The magnitudes we derive for this star agree with those obtained by Price et al. (2001) to within 0.10 mag in $`B,V`$, and $`R`$, and to within $`0.20`$ mag in $`I`$. Galama et al. (2001) provide further details on the image reduction and magnitude calibration.

From flattening in the late-time ($`t10`$ days) optical lightcurve we infer that the HST data are contaminated by a contribution from the host galaxy. We do not actually resolve the host in our HST images, and we observe flattening even with our small (0.14′′) extraction radius. This implies that the host itself is compact, or that the constant emission is from a bright knot or filament associated with the host. Future HST observations scheduled for late May 2001 may resolve this feature, and we will present these results elsewhere. To determine the flux from the OT we subtract this constant contribution.

To derive the magnitude for the constant term, we fit the ground-based optical lightcurve from Price et al. (2001) combined with our HST data to a decaying function with a constant added:

$$F(t,\nu )=F_0\nu ^\beta \left((t/t_{})^{\alpha _1s}+(t/t_{})^{\alpha _2s}\right)^{1/s}+F_{\mathrm{host}}(\nu ).$$
(1)
The decaying function, taken from Beuermann et al. (1999), has no physical significance, but provides a simple and general parametric description of the data. Here $`\alpha _1`$ and $`\alpha _2`$ are the early and late time asymptotic temporal slopes respectively, $`t_{}`$ is the time of the temporal slope break, $`s`$ is a parameter that determines the sharpness of the transition, and $`\beta `$ is the spectral slope. From the combined ground-based and HST data we derive $`\alpha _1=1.46\pm 0.11`$, $`\alpha _2=2.38\pm .07`$, $`t_{}=1.8\pm 0.1`$ days, $`\beta =1.49\pm 0.07`$. All errors in this paper are 1-$`\sigma `$. The best-fit value for the constant is $`F_{host}(B)=25.77\pm 0.14`$, $`F_{host}(V)=25.60\pm 0.1`$, $`F_{host}\left(R\right)=25.38\pm 0.1`$, $`F_{host}\left(I\right)=25.12\pm 0.10`$. We note that our fit values are similar to those derived by Price et al. (2001), but the late-time, high-resolution HST data allow us to determine the flux from the compact region near the OT.

We derive the OT magnitudes by subtracting the constant term, $`F_{host}`$, given above from each band, correcting for the foreground Galactic extinction of $`E_{BV}=0.0235`$ mag (Schlegel, Finkbeiner & Davis, 1998). Table 2 contains the resulting dereddened, host-subtracted OT flux values converted to $`\mu `$Jy. The measurement errors include a contribution (1–$`\sigma `$) of 0.026(B), 0.016 (V), 0.018 (R), 0.020 (I) $`\mu `$Jy added in quadrature to the statistical and systematic calibration errors to reflect the uncertainty in the constant host term derived from the fit.

### 2.2 CXO Observations

The BeppoSAX Medium Energy Concentrating Spectrometer (MECS) instrument (sensitive from 1.6 – 10 keV) discovered the X-ray afterglow of GRB 000926 in an observation made on Sep 29.03 – 29.52 (Piro et al., 2000). Garmire et al. (2000) observed the source for 10 ksec as part of a Cycle 1 CXO program using the ACIS S3 backside-illuminated chip on Sep 29.674 – 29.851 (referred to here as CXO-P1). Supporting our HST cycle 9 program, CXO again observed GRB 000926 in a 33-ksec long TOO taken Oct 10.176–10.76 UT (CXO-P2), also with the transient placed on the ACIS S3 backside-illuminated chip. The afterglow was clearly detected in each of these observations. The transient was not detected in the BeppoSAX Low Energy Concentrating Spectrometer (co-aligned with the MECS, and sensitive in the 0.1 – 10 keV band), however the effective exposure time was only 5 ksec. Piro et al. (2001) present a detailed spectral analysis of the X-ray data, including combined fits to all pointings.

We analyzed the data from CXO-P1 and CXO-P2 using software provided by the Chandra X-ray Observatory Center (CXC) to filter the events, extract source counts and subtract background. For the latter we used an annular region surrounding the source. We fit the two Chandra pointings separately using a powerlaw model plus absorption. The best fit photon spectral index for the CXO-P1 data is $`\alpha =1.9\pm 0.21`$ and for CXO-P2 is $`\alpha =2.23\pm 0.34`$, with best-fit hydrogen column depth $`N_H=(4.8\pm 3)\times 10^{20}`$ cm-2 (CXO-P1) and $`N_H=(3.0\pm 2.5)\times 10^{20}`$ cm-2 (CXO-P2). We find no evidence for additional absorption in the X-ray spectrum at the redshift of the host ($`z=2.04`$). This is consistent with the host $`A(V)`$ of 0.1 mag derived from model fits to the optical data (see below), if we adopt typical dust to gas ratios. For an LMC-like extinction curve, for example, this $`A(V)`$ corresponds to column depth of $`N_H=2\times 10^{20}`$ cm-2 for a host at $`z=2.04`$. Given the host redshift, such a low column is not detectable in the X-ray spectrum, and so no correction for absorption in the GRB host is necessary.

The detection statistics in the X-ray are limited, so for the purposes of our modeling we divided the energy range into soft (0.2 – 1.5 keV) and hard (1.5 – 8 keV) bands. We converted counts to flux using exposure maps weighted in energy using a photon powerlaw spectral index of $`\alpha =2`$ and $`N_H=2.7\times 10^{20}`$ cm-2, corresponding to the value for our own Galaxy, as determined by W3nH222W3nH is available at http://heasarc.gsfc.nasa.gov.. We adopt this value since it is consistent with that derived from our spectral modeling, and we find no significant evidence for additional absoprtion in the host. For the SAX observation, we used the best fit spectral model to determine the hard band flux (the MECS response does not allow a soft-band flux to be determined).

Table 3 shows the flux values, not corrected for Galactic absorption, for the three observations. To determine the center energy for the band, we took the mean, weighted using a photon index of $`\alpha =2`$. We determined flux errors by adding in quadrature the statistical error and the error due to the uncertainty in spectral slope. Since the flux is decaying during each observation interval, we weight the time of the observation by $`t^2`$ and average to determine the mean epoch. Table 4 lists the center frequencies and flux values for the X-ray transient (converted to $`\mu `$Jy), where we have corrected for Galactic absorption.

### 2.3 Radio Observations

We obtained radio observations at frequencies from 1.43 to 100 GHz at a number of facilities: the Owens Valley Radio Observatory Interferometer (OVRO), the Ryle Telescope, and the VLA. Table 5 summarizes these observations, organized by frequency. The 98.48 GHz data were taken at OVRO, the 15 GHz data at Ryle, and all other frequencies were observed using the VLA. In performing the observations, reducing the data and deriving flux errors we adopted the methodology described in detail in Kulkarni et al. (1999), Frail et al. (2000), and Berger et al. (2000). For the OVRO observations, we used the quasar 1726+455 for phase calibration and Uranus for absolute flux calibration.

## 3 Afterglow Model

To interpret the broadband lightcurve we employ a model for the emission from the relativistic shock (resulting from the GRB explosion) as it expands into the surrounding medium. Our modeling is quite comprehensive. It allows for either a constant density medium (which we refer to as ISM), or for a density depending on radius as $`\rho =Ar^2`$ referred to as the wind model, since this profile would be produced by a stellar progenitor’s wind (Chevalier & Li, 1999). We also allow for either isotropic or collimated ejecta, and we fit for extinction in the host galaxy. We include the effects of inverse Compton scattering on the evolution of the synchrotron spectrum, as well as the contribution of IC emission to the observed spectrum. IC scattering is not included in most models used to fit GRB afterglows, but if the ambient medium is of sufficiently high density, it becomes an important effect. Compton emission can dominate the total afterglow cooling rate for months after the event, and can even be directly observed in the X-ray (Panaitescu & Kumar, 2000; Sari & Esin, 2001).

We calculate the shock emission at a given time and frequency from a number of fundamental parameters: the initial isotropic-equivalent energy, $`E`$ in the shock, the electron powerlaw slope, $`p`$, the electron and magnetic field energy fractions, $`ϵ_e`$ and $`ϵ_B`$, and the density of the circumburst medium. We employ the equations given by Sari, Piran & Narayan (1998) and Granot, Piran & Sari (1999a, b) to determine the evolution of the synchrotron spectrum for the constant density medium, and the equations given by Chevalier & Li (1999, 2000) for the wind model. We include the effect of inverse Compton scattering using the treatment given by Sari & Esin (2001). We calculate the time at which collimation becomes evident ($`\theta =1/\mathrm{\Gamma }`$, where $`\mathrm{\Gamma }`$ is the Lorentz factor of the ejecta) and the subsequent spectral evolution using Sari, Piran & Halpern (1999). For the host extinction we use the parameterization of A(V) of Cardelli et. al. (1989) and Fitzpatrick & Massa (1998), joined using the formula from Reichart (1999) in the rest-frame optical/UV, and absorption cross-sections given by Morrison & McCammon (1983) in the rest-frame X-ray. We consider both LMC or SMC-like extinction laws. To fit the broadband lightcurve we vary the fundamental parameters, iterating to minimize $`\chi ^2`$.

At any given time the synchrotron emission, which dominates at most energies, has a spectrum characterized by a number of breaks: a low-frequency roll over due to self-absorption at frequency $`\nu _a`$; a peak at frequency $`\nu _m`$ due to the peak energy in the electron distribution; and a cooling break at $`\nu _c`$, where radiative energy losses become a significant fraction of the electron energy. These breaks evolve in time as the shock evolves, producing frequency-dependent breaks in the lightcurve. For the case of collimated ejecta, the lightcurve will begin to steepen when $`\theta `$ is smaller than the angle into which the emission is beamed due to relativistic effects ($`1/\mathrm{\Gamma }`$, where $`\mathrm{\Gamma }`$ is the Lorentz factor). We use a sharp transition for this “jet break”, and employ spherical evolution laws up until the time of the break. After the jet transition we assume the ejecta expand laterally, and we employ the appropriate asymptotic formulae until the time the ejecta become non-relativistic. The transition to the non-relativistic evolution is again treated as sharp.

The inverse Compton component is also characterized by a series of spectral breaks. With the synchrotron spectrum as the photon source, Sari & Esin (2001) find that the upscattered component is characterized by three breaks at $`\nu _a^{IC}2\mathrm{m}\mathrm{i}\mathrm{n}(\gamma _c,\gamma _m)^2\nu _a`$, $`\nu _m^{IC}2\gamma _m^2\nu _m`$, $`\nu _c^{IC}2\gamma _c^2\nu _c`$ where $`\gamma _m`$ is the characteristic Lorentz factor of the electrons emitting with peak frequency $`\nu _m`$. We adopt a spectral shape for the IC component similar to the synchrotron spectrum, but with the breaks given by the values in Sari & Esin (2001). This provides a good approximation to the IC component spectral shape except for at $`\nu <\nu _a^{IC}`$, but for our data this is not a concern (see below). Between $`\nu _m^{IC}`$ and $`\nu _c^{IC}`$ there is only a logarithmic correction, which we ignore.

## 4 Fits to the Broadband GRB 000926 Lightcurve

We fit the broadband GRB 000926 lightcurve employing the afterglow model described above. In addition to the data in Tables 1 – 4, we included the ground-based optical points given in Price et.al. (2001), as well as $`J`$ band data from DiPaola et al. (2000) and $`J,H,`$ $`K^{}`$ and $`K`$-band points from Fynbo et al. (2001). We converted the $`K^{}`$ points to $`K`$ using the prescription in Wainscoat & Cowie (1992). All the data are corrected for absorption in our Galaxy, and in the case of the optical data we subtracted the contribution from the host. The uncertainty in the radio flux is dominated by interstellar scintillation over the timescales of interest here. In calculating the $`\chi ^2`$ for the fit, we add an uncertainty to the radio fluxes based on the scintillation envelopes calculated using the fireball size derived from the model. We have excluded the first 8.46 data point from the fits, as it is likely associated with a separate component from the reverse shock (see below).

We considered several possible cases in fitting the data; constant density as well as a density gradient for the surrounding medium, and two cases for the IC emission. The need to consider two cases arises from the fact that for the same sychrotron emission spectrum, Sari & Esin (2001) find two possible IC solutions. These correspond to the two limits $`\eta ϵ_e/ϵ_Bf<<1`$, and $`f>>1`$, where $`ϵ_e`$ and $`ϵ_B`$ are the fractions of the total explosion energy that go into accelerating shocked electrons and amplifying the post-shock magnetic field respectively. Here $`\eta `$ is the fraction of electron energy radiated away, so describes whether cooling is dominant. If $`f<<1`$ (low IC), the IC cooling rate is unimportant compared to that of synchrotron, whereas if $`f>>1`$ (high IC) inverse Compton cooling dominates the total emission.

Table 6 summarizes the physical parameters corresponding to the best fit for each model. All models we considered required the ejecta to be collimated in order to explain the break in the optical lightcurve at 1–2 days. In the optical, this break is roughly independent of frequency. Extinction of $`A(V)0.10.2`$ mag, presumably due to the host galaxy, is required to fit the optical spectrum in three out of four cases. We employed an LMC-like curve for the models, however the magnitude of the extinction is relatively insensitive to the extinction law. Price et al. (2001) derived similar values for collimation and extinction based on fits to the optical data alone.

The model which best represents the broadband data is a constant density ISM with the high-IC solution. Figures 1a and 1b show the broadband spectrum around day 2 and day 10, with the best-fit model overplotted. The dot-dashed lines show the separate contributions from the synchrotron and IC components. Synchrotron emission dominates at optical frequencies and below, however on both day 2 and day 10 the model predicts an equal or larger contribution to the X-ray data from IC scattering. From Figure 1 it is evident that the approximation to the IC contribution made by assuming the same shape as the synchrotron is not an issue in fitting the data. As mentioned previously, the error in this approximation is only significant below $`\nu _a^{IC}`$, and for the parameters we derive the synchrotron dominates the total flux in this region by several orders of magnitude. Figures 2, 3, and 4 show the X-ray, optical and radio lightcurves. The $`\chi ^2`$ for the fit is 124 for 114 degrees of freedom.

The first column of data in Table 6 lists the physical parameters for the best fit high-IC solution. The derived parameters are all reasonable, although we note that $`ϵ_e`$ of $`0.30`$ is only marginally consistent with the assumption of adiabatic evolution inherent to our model. While this will not alter our basic conclusions, in future work we will expand the model to include partially radiative blast waves (see Sari (1997) and Cohen, Piran & Sari (1998) for treatments of the afterglow evolution in this case).

The isotropic-equivalent energy in the blast wave ($`E_{iso}`$ in Table 6) of $`1.8\times 10^{53}`$ erg is valid for $`t\mathrm{}>1`$ day. We use an adiabatic model to fit the data, which assumes the total energy is constant in time. While this is approximately valid at $`t\mathrm{}>12`$ days (i.e. for the majority of our data), the energy we derive will not represent the total initial blastwave energy, since significant energy is lost to radiation at early times. In particular, for the high-IC case, the shock wave radiates a significant amount of energy early on, emitted as MeV gamma-rays. Using the equations in Cohen, Piran & Sari (1998) we estimate the energy in the initial fireball to be about five times higher than that associated with the adiabatically expanding blast wave at a few days, or $`E_{\mathrm{initial}}1\times 10^{54}`$ erg. The observed isotropic energy release in the gamma-ray burst itself was $`E_\gamma =3\times 10^{53}`$ erg (Hurley et al., 2000). Recent theoretical work suggests the internal shock process producing the GRB should radiate less energy than the blast wave, except under extreme physical conditions (Beloborodov, 2000; Kobayahi & Sari, 2001). The value we derive for the fireball energy is, therefore, reasonable.

To fit the low-IC ISM model, we initially allowed all of the physical parameters to vary without limitation. For the low-IC case, the observed emission at all frequencies is from synchrotron radiation alone. The fitting program adjusts the physical parameters in order to reproduce all data with this component. Reproducing the broadband observations requires a large fraction of energy in the post-shock magnetic field. The best-fit value of $`ϵ_B=24`$ is unphysical, and so we have fixed $`ϵ_B`$ at unity – the maximum value physically possible. For either value of $`ϵ_B`$ the low-IC model provides a poor fit to the data. The $`\chi ^2`$ is 187 for 114 datapoints for $`ϵ_B=1`$ and 178 for 114 datapoints for $`ϵ_B=24`$. The largest contribution to the $`\chi ^2`$ comes from the optical data, where the syncrotron component does not fit the data as well as that derived for the high-IC case.

Table 6 provides the physical parameters associated with the best-fit low-IC ISM model for $`ϵ_B`$ fixed at unity. In addition to the large $`ϵ_B`$, the isotropic-equivalent energy in the blast wave of $`10^{53}`$ erg for the low-IC case is smaller than the observed gamma-ray energy release. Radiative corrections for this model are much smaller than for the high-IC case, so that the adiabatic value more closely represents the energy in the initial fireball. As pointed out above, $`E_\gamma \mathrm{}>E_{\mathrm{fireball}}`$ is hard to accomodate except under extreme conditions.

Figure 2 shows the X-ray data with both high and low-IC ISM models overplotted. The high-IC case provides a somewhat better fit to the data, as the low-IC case systematically underpredicts the measurements. The X-ray observations alone, however, are not sufficient in this case to indicate the presence of an inverse Compton component. This requires the results from the broadband fit. Our calculations show that the IC scattering hardens the X-ray spectrum over some time intervals compared to the prediction low-IC model. There is only marginal evidence for this in these data; the observed X-ray spectral photon index is $`1.9\pm 0.21`$ compared with the $`p/2+1=2.2`$ associated with the synchrotron component. This emphasizes the role of future, higher-significance X-ray observations in confirming the role if IC scattering in GRB afterglows.

We also investigated a wind density profile in both high- and low-IC regimes. Neither of these cases reproduce the data as well as the high-IC ISM model. The low-IC wind solution has a poor $`\chi ^2`$ of 244 for 114 degrees of freedom. Again we have fixed $`ϵ_B=1.0`$, since the best-fit value of 174 is unphysical. The high-IC wind case also has a large $`\chi ^2`$ of 167 for 114 d.o.f, and under-predicts the 5 GHz radio data by a factor 2 – 3. It does, however, reproduce the X-ray data reasonably well. We note that early-time optical or high-frequency radio data is particularly important in ruling out the wind model. This is because the lightcurve evolution after the jet break is similar for ISM and wind. Directly observing the cooling break evolution provides the greatest leverage, since this evolution is significantly different for wind and ISM models. Generally, this requires optical observations or high-frequency radio observations at early epochs (minutes to hours after the event). Unfortunately optical coverage of this event began rather late ($`t0.8`$ days), but future early-time observations of other events will provide a more sensitive probe for density structure in the surrounding medium.

Finally, we point out that even for our best fit (high-IC ISM) solution, the first 8 GHz data point exceeds the model’s expectation by a factor of 3.3 (7.9 $`\sigma `$). This discrepancy is worse for the alternative models. Scintillation gains this large are extremely improbable, and it is therefore likely that we are seeing evidence for an additional component due to the reverse shock emission, as was observed in early radio observations of GRB 990123 (Kulkarni et al., 1999).

## 5 Testing the High Density Model for GRB000926

Piro et al. (2001) analyzed the X-ray and optical lightcurves of GRB 000926, and noted that the spectral and temporal properties cannot be explained using a standard jet model with moderate ISM density. By invoking inverse Compton cooling, we were able to fit the data with a moderate density medium and relatively standard fireball parameters. In order to describe the X-ray and optical lightcurves, Piro et al. invoke moderate collimation of the fireball (opening angle $`\theta =25^{}`$), and a dense ($`n=4\times 10^4`$ cm-3) medium. This high density results in rapid deceleration of the fireball, and an early transition to the non-relativistic regime. These features seemingly reproduce the general characteristics of the high-frequency observations. The high density found by Piro et al. would suggest the event took place in a molecular cloud, and further the low extinction measured in the optical would indicate that the GRB destroyed much of the dust in its vicinity.

The radio data provide an important means of distinguishing these two interpretations. If the fireball expands in a very dense medium, the ejecta decelerate much more rapidly, and at a given time the fireball radius is significantly smaller than for the moderate densities derived in this paper. The self-absorption frequency, easily observed as a cutoff in the radio spectrum, is correspondingly much higher at a given time for the high-density case.

We have fit a model with the ISM density fixed at $`n=4\times 10^4`$ cm-3, and the jet opening angle of $`\theta =25^{}`$ to the broadband data set. We investigated both high-IC and low-IC constant-density solutions, and find we obtain the best fit using the high-IC model. Figure 5 shows the broadband spectrum on day 2, with a comparison of our best fit with the data. The radio data quite clearly exclude the high-density possibility, since the self-absorption frequency is more than a decade higher than we measure. On day two, the high density model predicts an 8 GHz radio flux of $`0.2\mu `$Jy, while we measure a flux of 200$`\mu `$Jy. This points out the importance of broadband coverage, including radio, in constraining models for the GRB environment.

## 6 Conclusions

GRB 000926 has one of the highest-quality broadband afterglow lightcurves studied to-date. We can explain all the general features using a model of a relativistic shock produced by ejecta collimated to a $`8^{}`$ opening angle. If we correct the isotropic-equivalent gamma-ray energy release for this jet angle we get $`E_\gamma =E_{\gamma ,iso}\times \theta ^2/2=3\times 10^{51}`$ erg. This is reasonably accounted for in most currently-popular progenitor models. Our observations require extinction in the host galaxy, with a corresponding hydrogen column of $`2\times 10^{20}`$ cm-2 if the dust to gas ratio is typical of the LMC or SMC. This is consistent with the event taking place in a galactic disk.

We find that an explosion occurring in a spatially homogeneous medium best describes the broadband lightcurves, provided that the ratio $`\eta ϵ_e/ϵ_B>>1`$. This condition implies that Compton emission dominates the total cooling rate, and that IC scattering is potentially directly observable in the afterglow spectrum. We do, in fact, find that for our best-fit model IC emission dominates over the synchrotron contribution at 1 keV on the 2 – 10 day timescales spanned by our X-ray data. Detecting the IC emission in the X-ray band directly implies a lower-limit on the density (Panaitescu & Kumar, 2000; Sari & Esin, 2001), in this case $`n\mathrm{}>10`$ cm-3, and our best fit yields $`n=30`$ cm-3. This value for $`n`$ is higher than the average ISM density in a typical galaxy, and is consistent with a diffuse interstellar cloud, such as those commonly found in star forming regions. Combining the density with the host $`N_H`$ derived from our optical observations (again assuming typical dust to gas ratios) implies a scale size of 2 pc for the cloud.

We wish to thank Harvey Tananbaum and the CXO operations staff for facilitating the Chandra TOO observations. We thank Steve Beckwith and the HST operations staff for facilitating the WFPC2 observations. FAH acknowledges support from a Presidential Early Career award. SRK and SGD thank NSF for support of their ground-based GRB programs. RS is grateful for support from a NASA ATP grant. RS & TJG acknowledge support from the Sherman Fairchild Foundation.