# The Afterglow of GRB 990123 and a Dense Medium

## 1. Introduction

The gamma-ray burst (GRB) 990123 was an extraordinary event. It was the brightest burst yet detected with the Wide Field Camera on the BeppoSAX satellite (Feroci et al. 1999), and had a total gamma-ray fluence of $`5\times 10^4\mathrm{erg}\mathrm{cm}^2`$, which is in the top 0.3% of all bursts. It was the first burst to be simultaneous detected in the optical band. Optical emission with peak magnitude of $`V9`$ was discovered by the Robotic Optical Transient Search Experiment (ROTSE) during the burst and was found to have rapidly faded down immediately after the gamma-ray emission (Akerlof et al. 1999). The detection of the redshift showed that the burst appears at $`z1.6`$ (Andersen et al. 1999; Kulkarni et al. 1999a). This implies that if the GRB emission was directed isotropically, the inferred energy release is $`1.6\times 10^{54}\mathrm{ergs}`$ (Kulkarni et al. 1999a; Briggs et al. 1999).

The burst’s afterglow was detected and monitored at X-ray, optical and radio bands. It was the brightest of all GRB X-ray afterglows observed until now. The BeppoSAX detected the flux of the afterglow at 2-10 keV six hours after the gamma-ray trigger to be $`1.1\times 10^{11}\mathrm{erg}\mathrm{cm}^2\mathrm{s}^1`$ and the subsequent temporal decay index to be $`\alpha _X=1.44\pm 0.07`$ (Heise et al. 1999a, b). The R-band optical afterglow about 3.5 hours after the burst showed a power-law decay with index $`\alpha _{1R}=1.1\pm 0.03`$ (Kulkarni et al. 1999a; Castro-Tirado et al. 1999; Fruchter et al. 1999). This law continued until about $`2.04\pm 0.46`$ days after the burst. Then the optical emission began to decline based on another power law with index $`\alpha _{2R}=1.65\pm 0.06`$ (Kulkarni et al. 1999a) or $`1.75\pm 0.11`$ (Castro-Tirado et al. 1999) or $`1.8`$ (Fruchter et al. 1999). In addition, a radio flare was also detected about 1 day after the burst (Kulkarni et al. 1999b; Galama et al. 1999).

A scenario has been proposed to explain these observations. If the burst is assumed to be produced from a jet, the steepening of the late optical afterglow decay is due to the possibility that this jet has undergone the transition from a spherical-like phase to a sideways-expansion phase (Rhoads 1997, 1999; Kulkarni et al. 1999a; Fruchter et al. 1999; Sari, Piran & Halpern 1999) or that we have observed the edge of the jet (Panaitescu & Mészáros 1998; Mészáros & Rees 1999).

In this Letter we propose another possible scenario, in which the steepening of the late optical afterglow decay is due to the shock which has evolved from a relativistic phase to a nonrelativistic phase in a dense medium. According to the standard afterglow shock model (for a review see Piran 1998), the afterglow is produced by synchrotron radiation or inverse Compton scattering in the external forward wave (blast wave) of the GRB fireball expanding in a homogeneous medium. The external reverse shock of the fireball may lead to a prompt optical flash (Sari & Piran 1999). As more and more ambient matter is swept up, the forward shock gradually decelerates and eventually enters a nonrelativistic phase. In the meantime, the emission from such a shock fades down, dominating at the beginning in X-rays and progressively at optical to radio energy band. There are two limiting cases (adiabatic and highly radiative) for the hydrodynamical evolution of the shock. These cases have been well studied both analytically (e.g., Mészáros & Rees 1997; Wijers, Rees & Mészáros 1997; Waxman 1997a, b; Reichart 1997; Sari 1997; Vietri 1997; Katz & Piran 1997; Mészáros, Rees & Wijers 1998; Dai & Lu 1998a; Sari, Piran & Narayan 1998; etc) and numerically (e.g., Panaitescu, Mészáros & Rees 1998; Huang et al. 1998; Huang, Dai & Lu 1998). A partially radiative (intermediate) case has been investigated (Chiang & Dermer 1998; Cohen, Piran & Sari 1998; Dai, Huang & Lu 1999). Here we only consider the limiting cases. In the highly radiative model, since all shock-heated electrons cool faster than the age of the shock, the optical afterglow should have the same temporal decay index as the X-ray afterglow (Sari et al. 1998), incompatible with the observations (Kulkarni et al. 1999a). In the adiabatic model, however, the difference in the decay index between optical and X-ray afterglows is found to be likely $`1/4`$, which is consistent with the observational result $`\mathrm{\Delta }\alpha =\alpha _{1R}\alpha _X0.3`$. This implies that the shock producing the afterglow of GRB 990123 has evolved adiabatically. This is the starting point of our analysis. For an adiabatic shock, the time at which it enters a nonrelativistic phase $`n^{1/3}`$, where $`n`$ is the baryon number density of the medium. Therefore, this time for a shock expanding in a dense medium with density of $`n10^6\mathrm{cm}^3`$ is two orders of magnitude smaller than that for a shock with the same energy in a thin medium with density of $`n1\mathrm{cm}^3`$. Furthermore, as given in Section 2, the afterglow at the nonrelativistic phase decays faster than at the relativistic phase. It is natural to expect that this effect can provide an explanation for the steepening feature of the afterglow from GRB 990123.

Dense media have been discussed in the context of GRBs. First, Katz (1994) suggested collisions of relativistic nucleons with a dense cloud as an explanation of the delayed hard photons from GRB 940217. Second, to explain the radio flare of GRB 990123, Shi & Gyuk (1999) speculated that a relativistic shock may have ploughed into a dense medium off the line of sight. Third, Piro et al. (1999) and Yoshida et al. (1999) have reported an iron emission line in the X-ray afterglow spectrum of GRB 970508 and GRB 970828 respectively. The observed line intensity requires a dense medium with a large iron mass concentrated in the vicinity of the burst (Lazzati, Campana & Ghisellini 1999). Finally, dense media (e.g., clouds or ejecta) may appear in the context of some energy source models, e.g., failed supernovae (Woosley 1993), hypernovae (Paczyński 1998), supranovae (Vietri & Stella 1998), phase transition of neutron stars to strange stars (Dai & Lu 1998b), baryon decay of neutron stars (Pen & Loeb 1998), etc.

## 2. The Evolution of a Shock in a Dense Medium

### 2.1. Relativistic Phase

Now we consider an adiabatic relativistic shock expanding in a dense medium. The Blandford-McKee (1976) self similar solution gives the Lorentz factor of the shock,

$`\gamma `$ $`=`$ $`{\displaystyle \frac{1}{4}}\left[{\displaystyle \frac{17E(1+z)^3}{\pi nm_pc^5t_{}^3}}\right]^{1/8}`$ (1)   
$`=`$ $`2E_{54}^{1/8}n_5^{1/8}t_{\mathrm{day}}^{3/8}[(1+z)/2.6]^{3/8},`$   
where $`E=E_{54}\times 10^{54}\mathrm{ergs}`$ is the total isotropic energy, $`n_5=n/10^5\mathrm{cm}^3`$, $`t_{}=t_{\mathrm{day}}\times 1\mathrm{day}`$ is the observer’s time since the gamma-ray trigger, $`z`$ is the the redshift of the source generating this shock, and $`m_p`$ is the proton mass.

In analyzing the spectrum and light curve of synchrotron radiation from the shock, one needs to know two crucial frequencies: the synchrotron radiation peak frequency ($`\nu _m`$) and the cooling frequency ($`\nu _c`$). In the standard afterglow shock picture, the electrons heated by the shock are assumed to have a power-law distribution: $`dN_e/d\gamma _e\gamma _e^p`$ for $`\gamma _e\gamma _{em}`$, where $`\gamma _e`$ is the electron Lorentz factor and the minimum Lorentz factor $`\gamma _{em}=610ϵ_e\gamma `$. The power-law index $`p2.56`$ by fitting the spectrum and light curve of the observed afterglow of GRB 990123 (see below). We further assume that $`ϵ_e`$ and $`ϵ_B`$ are ratios of the electron and magnetic energy densities to the thermal energy density of the shocked medium respectively. Based on these assumptions, the synchrotron radiation peak frequency in the observer’s frame can be written as

$`\nu _m`$ $`=`$ $`{\displaystyle \frac{\gamma \gamma _{em}^2}{1+z}}{\displaystyle \frac{eB^{}}{2\pi m_ec}}`$ (2)   
$`=`$ $`8.0\times 10^{11}ϵ_e^2ϵ_{B,6}^{1/2}E_{54}^{1/2}t_{\mathrm{day}}^{3/2}`$   
$`\times [(1+z)/2.6]^{1/2}\mathrm{Hz},`$   
where $`ϵ_{B,6}=ϵ_B/10^6`$ and $`B^{}=(32\pi ϵ_B\gamma ^2nm_pc^2)^{1/2}`$ is the internal magnetic field strength of the shocked medium. According to Sari et al. (1998), the cooling frequency, the frequency of electrons with Lorentz factor of $`\gamma _c`$ that cool on the dynamical time of the shock, is given by

$`\nu _c`$ $`=`$ $`{\displaystyle \frac{\gamma \gamma _c^2}{1+z}}{\displaystyle \frac{eB^{}}{2\pi m_ec}}={\displaystyle \frac{18\pi em_ec(1+z)}{\sigma _T^2B^3\gamma t_{}^2}}`$ (3)   
$`=`$ $`1.9\times 10^{16}ϵ_{B,6}^{3/2}E_{54}^{1/2}n_5^1t_{\mathrm{day}}^{1/2}`$   
$`\times [(1+z)/2.6]^{1/2}\mathrm{Hz},`$   
where $`\sigma _T`$ is the Thompson scattering cross section. From equations (2) and (3), Sari et al. (1998) have further defined two critical times, when the breaking frequencies $`\nu _m`$ and $`\nu _c`$ cross the observed frequency $`\nu =\nu _{15}\times 10^{15}\mathrm{Hz}`$: $`t_m=8.6\times 10^3ϵ_e^{4/3}ϵ_{B,6}^{1/3}E_{54}^{1/3}[(1+z)/2.6]^{1/3}\nu _{15}^{2/3}\mathrm{days}`$, and $`t_c=380ϵ_{B,6}^3E_{54}^1n_5^2[(1+z)/2.6]^1\nu _{15}^2\mathrm{days}`$. Therefore we see that for $`E_{54}1.6`$, $`ϵ_e0.1`$, $`ϵ_{B,6}0.02`$, and $`n_530`$ inferred in the next section, the optical afterglow in several days after the burst should result from those slowly-cooling electrons and the X-ray afterglow from those fastly-cooling electrons.

The observed synchrotron radiation peak flux can be obtained by

$`F_{\nu _m}`$ $`=`$ $`{\displaystyle \frac{N_e\gamma P_{\nu _m}^{}(1+z)}{4\pi D_L^2}}`$ (4)   
$`=`$ $`4.2ϵ_{B,6}^{1/2}E_{54}n_5^{1/2}[(1+z)/2.6]D_{L,28}^2\mathrm{Jy},`$   
where $`N_e`$ is the total number of swept-up electrons, $`P_{\nu _m}^{}=m_ec^2\sigma _TB^{}/(3e)`$ is the radiated power per electron per unit frequency in the frame comoving with the shocked medium, and $`D_L=D_{L,28}\times 10^{28}\mathrm{cm}`$ is the distance to the source. In the light of equations (2)-(4), one can easily find the spectrum and light curve of the afterglow,

$$F_\nu =\{\begin{array}{cccc}(\nu /\nu _m)^{(p1)/2}F_{\nu _m}\hfill & & & \\ \nu ^{(p1)/2}t_{}^{3(1p)/4}\mathrm{if}\nu _m<\nu <\nu _c;\hfill & & & \\ (\nu _c/\nu _m)^{(p1)/2}(\nu /\nu _c)^{p/2}F_{\nu _m}\hfill & & & \\ \nu ^{p/2}t_{}^{(23p)/4}\mathrm{if}\nu >\nu _c,\hfill & & & \end{array}$$
(5)   
where the low-frequency radiation component has not been considered (Sari et al. 1998). In the GRB 990123 case, we require $`\nu _m<\nu <\nu _c`$ for the optical afterglow and $`\nu >\nu _c`$ for the X-ray afterglow. Thus, the R-band afterglow decay index $`\alpha _R=3(1p)/4`$ and the X-ray decay index $`\alpha _X=(23p)/4`$, which are well consistent with the observational results $`\alpha _{1R}=1.1\pm 0.03`$ and $`\alpha _X=1.44\pm 0.07`$ if $`p2.56`$.

### 2.2. Nonrelativistic Phase

As it sweeps up sufficient ambient matter, the shock will eventually go into a nonrelativistic phase. During such a phase, the shock’s velocity $`vt_{}^{3/5}`$, its radius $`rt_{}^{2/5}`$, the internal field strength $`B^{}t_{}^{3/5}`$ and the typical electron Lorentz factor $`\gamma _{em}t_{}^{6/5}`$. Thus, we obtain the synchrotron peak frequency $`\nu _m\gamma _{em}^2B^{}t_{}^3`$, the cooling frequency $`\nu _cB^3t_{}^2t_{}^{1/5}`$ and the peak flux $`F_{\nu _m}N_eP_{\nu _m}^{}r^3B^{}t_{}^{3/5}`$. According to these scaling laws, we further derive the spectrum and light curve at the nonrelativistic stage:

$$F_\nu =\{\begin{array}{cccc}(\nu /\nu _m)^{(p1)/2}F_{\nu _m}\hfill & & & \\ \nu ^{(p1)/2}t_{}^{(2115p)/10}\mathrm{if}\nu _m<\nu <\nu _c;\hfill & & & \\ (\nu _c/\nu _m)^{(p1)/2}(\nu /\nu _c)^{p/2}F_{\nu _m}\hfill & & & \\ \nu ^{p/2}t_{}^{(43p)/2}\mathrm{if}\nu >\nu _c.\hfill & & & \end{array}$$
(6)   
From equation (6), we can see the R-band decay index $`\alpha _R=(2115p)/10`$ for radiation from slowly-cooling electrons or $`\alpha _R=(43p)/2`$ for radiation from rapidly-cooling electrons. If $`p2.56`$, then $`\alpha _R1.74`$ or $`1.84`$, in excellent agreement with the observations in the time interval of 2.5 days to 20 days after the burst (Kulkarni et al. 1999a; Fruchter et al. 1999; Castro-Tirado et al. 1999).

## 3. Constraints on Parameters

In the above section, we show that as an adiabatic shock expands in a dense medium from an ultrarelativistic phase to a nonrelativistic phase, the decay of the radiation from such a shock will steepen. This effect may fit the observed steepening better than the alternative interpretation — jet sideways expansion. In the latter interpretation, the temporal decay of a late afterglow is very likely to be $`t_{}^p`$ (Rhoads 1997, 1999; Sari et al. 1999). We further analyze our effect and infer some parameters of the model.

According to the analysis on the R-band light curve of the GRB 990123 afterglow (Kulkarni et al. 1999a; Fruchter et al. 1999; Castro-Tirado et al. 1999), the observed break occurred at $`t_{}=2.04\pm 0.46`$days. This implies $`\gamma 1`$ at $`t_{\mathrm{day}}2.5`$. From equation (1), therefore, we find $`n_516E_{54}`$, where the redshift $`z=1.6`$ has been used. We now continue to consider two observational results. First, on January 23.577 UT, the Palomar 60-inch telescope detected the R-band magnitude $`R=18.65\pm 0.04`$, corresponding to the flux $`F_R100\mu \mathrm{Jy}`$ at $`t_{\mathrm{day}}0.17`$ (Kulkarni et al. 1999a). Considering this result in equation (5) together with equations (2) and (4), we can derive

$$ϵ_e^{p1}ϵ_{B,6}^{(p+1)/4}E_{54}^{(p+3)/4}n_5^{1/2}0.01,$$
(7)   
where the right number has been obtained by taking $`p2.56`$ and $`D_{L,28}3.7`$. Second, on January 24.65 UT, the BeppoSAX observed the X-ray (2-10 keV) flux $`F_X5\times 10^2\mu `$Jy (Heise et al. 1999a, b). Combining this result with equations (2)-(5), we can also derive

$$ϵ_e^{p1}ϵ_{B,6}^{(p2)/4}E_{54}^{(p+2)/4}0.03.$$
(8)   
Since $`E_{54}1.6`$ (Briggs et al. 1999; Kulkarni et al. 1999a), the medium density $`n_530`$ and the solution of equations (7) and (8) is $`ϵ_e0.1`$ and $`ϵ_{B,6}0.02`$. Our inferred value of $`ϵ_e`$ is near the equipartition value, in agreement with the result of Wijers & Galama (1998) and Granot, Piran & Sari (1998), while our $`ϵ_B`$ is about six orders of magnitude smaller than the value inferred from the afterglow of GRB 970508. Of course, the field density for GRB 971214 has been estimated to be less than $`10^5`$ times the equipartition value (Wijers & Galama 1998). As suggested by Galama et al. (1999), such differences in field strength may reflect differences in energy flow from the central engine.

## 4. Discussion and Conclusion

In the above section, we find the medium density $`n3\times 10^6\mathrm{cm}^3`$ for our model to fit the observed optical and X-ray afterglow of GRB 990123. Now we show that even in the the presence of such a dense medium, the optical and X-ray radiations from the forward shock were neither self absorbed in the shocked medium nor scattered in the unshocked medium. First, the self-absorption frequency of the shocked medium is (Wijers & Galama 1998; Granot et al. 1998) $`\nu _a10^3\mathrm{GHz}(ϵ_e/0.1)^1(ϵ_{B,6}/0.01)^{1/5}E_{54}^{1/5}(n_5/10)^{3/5}`$. This estimate should be the upper limit because of the presence of a possible low-energy electron population (Waxman 1997b). Clearly, $`\nu _a`$ is much less than the optical frequency, implying that the self absorption in the shocked medium didn’t affect the optical and X-ray afterglow. In fact, this estimate is valid only for $`\nu _a<\nu _m`$. When $`\nu _a>\nu _m`$, $`\nu _a`$ must have decayed. As a result, the flux at 8.46 GHz first increased as $`t_{}^{1.25}`$ and then declined as $`t_{}^{1.74}`$ for $`\nu _a<8.46`$ GHz during the nonrelativistic phase. This might provide an explanation for the observed radio flare. Second, a photon emitted from the shock may be scattered by the electrons in the unshocked medium. The scattering optical depth $`\tau \sigma _TnR`$ (where $`R`$ is the typical radius of the medium). If the medium was distributed isotropically and homogeneously and its mass $`M10M_{}`$ (the typical mass of a supernova ejecta), then $`\tau 0.05(M/10M_{})^{1/3}(n_5/10)^{2/3}1`$. This implies that the afterglow from the shock was hardly affected by the medium.

For other well-studied afterglows, e.g., GRB 970228 and GRB 970508, their ambient densities must be very low for three reasons: (i) In these bursts there was no observed break in the optical light curve as long as the afterglow could be observed (Fruchter et al. 1998; Zharikov et al. 1998). (ii) The fluctuation appearing in the radio afterglow light curve of GRB 970508 requires the shock had been relativistic for several weeks (Waxman, Kulkarni & Frail 1998). (iii) The analysis of the afterglow spectrum of GRB 970508 leads to a low ambient density $`n<10\mathrm{cm}^3`$ (Wijers & Galama 1998; Granot et al. 1998). However, the observed iron emission line in the X-ray afterglow spectrum of GRB 970508 indeed requires a dense medium with density $`10^9\mathrm{cm}^3`$ (Lazzati et al. 1999). The only way to reconcile a monthly lasting power-law afterglow with iron line emission is through a particular geometry, in which the line of sight is devoid of the dense medium. In contrast to this idea, we suggest that for GRB 990123 a dense medium of $`n3\times 10^6\mathrm{cm}^3`$ appears at least at the line of sight or perhaps isotropically.

How was the dense medium produced? One possibility was a cloud and another possibility was an ejecta from the GRB site. There have been several source models (mentioned in Introduction) in the literature which may lead to massive ejecta. Here we want to discuss one of them in detail. Timmes, Woosley & Weaver (1996) showed that Type II supernovae may produce a kind of neutron star with $`1.73M_{}`$. If these massive neutron stars have very short periods at birth, they may subsequently convert into strange stars due to rapid loss of angular momenta (Cheng & Dai 1998), and perhaps the strange stars are differentially rotating (Dai & Lu 1998b). Even though this model is somewhat similar to the supranova model of Vietri & Stella (1998), resultant compact objects are strange stars in our model and black holes in the supranova model. We further discuss implications of our model. First, the model leads to low-mass loading matter because of thin baryonic crusts of the strange stars. Second, such stars result in GRBs with spiky light curves, being consistent with the analytical result from the observed data of GRB 990123 (Fenimore, Ramirez-Ruiz & Wu 1999). The third advantage of this model is to be able to explain well the property of the early afterglow of GRB 970508 by considering energy injection from the central pulsar (Dai & Lu 1998b, c). Finally, a dense medium, the supernova ejecta, appears naturally.

Our scenario proposed in this Letter requires a dense medium with density $`3\times 10^6\mathrm{cm}^3`$ to explain the steepening in the temporal decay of the R-band afterglow about 2.5 days after GRB 990123. We also suggest that this medium could be a supernova/supranova/hypernova ejecta. Thus, if the mass of the medium is assumed to be $`M10M_{}`$, its radius can be estimated to be $`R3\times 10^{17}\mathrm{cm}(M/10M_{})^{1/3}(n_5/10)^{1/3}`$. According to equation (1), we can integrate $`dr=2\gamma ^2cdt_{}`$ and thus find that the postburst 2.5-day time in the observer’s frame corresponds to about 20 days in the unshocked medium’s frame. This implies that the radius at which the shock entered a nonrelativistic phase is about $`5\times 10^{16}`$ cm. This radius is much less than that of the medium. Therefore, the medium discussed here was so wide and dense that the ultrarelativistic shock must have become nonrelativistic about 2.5 days after the burst.

In summary, a simple explanation for the “steepening” observed in the temporal decay of the late R-band afterglow of GRB 990123 is that a shock expanding in a dense medium with density of $`3\times 10^6\mathrm{cm}^3`$ has evolved from a relativistic phase to a nonrelativistic phase. We find that this scenario not only explains well the optical afterglow but also accounts for the observed X-ray afterglow quantatitively.

We would like to thank J. I. Katz, S. R. Kulkarni, A. Mitra and the anonymous referee for invaluable suggestions, and Y. F. Huang and D. M. Wei for helpful discussions. This work was supported by the National Natural Science Foundation of China (grants 19825109 and 19773007).