# Overall Evolution of Realistic Gamma-Ray Burst Remnant and Its Afterglow*footnote **footnote *Supported by the National Natural Science Foundation of China under Grant Nos. 19773007 and 19825109, and the National Climbing Project on Fundamental Researches.

## Abstract

Conventional dynamic model of gamma-ray burst remnants is found to be incorrect for adiabatic blastwaves during the non-relativistic phase. A new model is derived, which is shown to be correct for both radiative and adiabatic blastwaves during both ultra-relativistic and non-relativistic phase. Our model also takes the evolution of the radiative efficiency into account. The importance of the transition from the ultra-relativistic phase to the non-relativistic phase is stressed.

PACS: 95.30.Lz, 98.70.Rz, 97.60.Jd

slugcomment: Published in: Chinese Physics Letters 16 (1999) 775  

The origin of gamma-ray bursts (GRBs) has remained unknown for over 30 years.<sup>1,2</sup> A major breakthrough appeared in early 1997, when the Italian-Dutch BeppoSAX satellite observed X-ray afterglows from GRB 970228 for the first time.<sup>3</sup> Since then, X-ray afterglows have been observed from about 15 GRBs, of which ten events were detected optically and five bursts were also detected in radio wavelengths. The cosmological origin of at least some GRBs is firmly established. The so called fireball model <sup>4,5</sup> is strongly favoured, which is found successful at explaining the major features of the low energy light curves.<sup>6-9</sup>

In the fireball model, low energy afterglows are generated by ultra-relativistic fireballs, which first give birth to GRBs through internal or external shock waves and then decelerate continuously due to collisions with the interstellar medium (ISM). The dynamics of the expansion has been investigated extensively.<sup>6-9</sup> Both analytic solutions and numerical approaches are available. It is a general conception that current models describe the gross features of the process very well and further improvements are possible only by considering some details. However, we find that three serious problems are associated with the popular model.

First, it is usually assumed that the expansion is ultra-relativistic. Then for an adiabatic fireball, the evolution of the bulk Lorentz factor is derived to be:

$$\gamma (200400)E_{51}^{1/8}n_1^{1/8}t^{3/8},$$
(1)   
where $`E_{51}=E_0/(10^{51}\mathrm{erg})`$ with $`E_0`$ the initial fireball energy, $`n_1=n/(1\mathrm{c}\mathrm{m}^3)`$ with $`n`$ the ISM number density, and $`t`$ is observer’s time in unit of s.<sup>6-9</sup> The radius of the blastwave scales as $`Rt^{1/4}`$. Based on Eq.(1), flux density at frequency $`\nu `$ then declines as $`S_\nu t^{3(1p)/4}`$, where $`p`$ is the index characterizing the power-law distribution of the shocked ISM electrons, $`dn_\mathrm{e}^{}/d\gamma _\mathrm{e}\gamma _\mathrm{e}^p`$. These expressions are valid only when $`\gamma 1`$.

However, optical afterglows from GRB 970228 and GRB 970508 were detected for as long as 190 and 260 d respectively, while in Eq.(1), even $`t=30`$ d will lead to $`\gamma 1`$. It is clear that the overall evolution of the postburst fireball can not be regarded as a simple one-phase process, we should pay special attention to the transition from the ultra-relativistic phase to the non-relativistic phase.<sup>10</sup> This is unfortunately ignored in the literature.

Second, the expansion of the fireball might be either adiabatic or highly radiative. Extensive attempts have been made to find a common model applicable for both cases.<sup>11-13</sup> As a result, a differential equation has been proposed by various authors,<sup>12,13</sup>

$$\frac{d\gamma }{dm}=\frac{\gamma ^21}{M},$$
(2)   
where $`m`$ is the rest mass of the swept-up ISM, $`M`$ is the total mass in the co-moving frame, including internal energy $`U`$. Since thermal energy produced during the collisions is $`dE=c^2(\gamma 1)dm`$, usually we assume:<sup>13</sup>

$$dM=\frac{(1ϵ)}{c^2}dE+dm=[(1ϵ)\gamma +ϵ]dm,$$
(3)   
where $`ϵ`$ is defined as the fraction of the shock generated thermal energy (in the co-moving frame) that is radiated. It is putative that Eq.(2) is correct in both ultra-relativistic and non-relativistic phase, for both radiative and adiabatic fireballs.

In the highly radiative case, $`ϵ=1`$, $`dM=dm`$, Eq.(2) reduces to,

$$\frac{d\gamma }{dm}=\frac{\gamma ^21}{M_{\mathrm{ej}}+m},$$
(4)   
where $`M_{\mathrm{ej}}`$ is the mass ejected from the GRB central engine. Then an analytic solution is available,<sup>11,13</sup> which satisfies $`\gamma R^3`$ when $`\gamma 1`$ and $`vR^3`$ when $`\gamma 1`$, where $`v`$ is the bulk velocity of the material. These scaling laws indicate that Eq.(2) is really correct for highly radiative fireballs. In the adiabatic case, $`ϵ=0`$, $`dM=\gamma dm`$, Eq.(2) also has an analytic solution:<sup>12</sup>

$$M=[M_{\mathrm{ej}}^2+2\gamma _0M_{\mathrm{ej}}m+m^2]^{1/2},$$
(5)   
$$\gamma =\frac{m+\gamma _0M_{\mathrm{ej}}}{M},$$
(6)   
where $`\gamma _0`$ is the initial value of $`\gamma `$. During the ultra-relativistic phase, Eqs.(5) and (6) do produce the familiar power-law $`\gamma R^{3/2}`$, which is often quoted for an adiabatic blastwave decelerating in a uniform medium. In the non-relativistic limit ($`\gamma 1`$, $`m\gamma _0M_{\mathrm{ej}}`$), Chiang and Dermer have derived $`\gamma 1+\gamma _0M_{\mathrm{ej}}/m`$,<sup>12</sup> so that they believe it also agrees with the Sedov solution, i.e., $`vR^{3/2}`$.<sup>14</sup> However we find that their approximation is not accurate,<sup>15</sup> because they have omitted some first-order infinitesimals of $`\gamma _0M_{\mathrm{ej}}/m`$. The correct approximation could be obtained only by retaining all the first and second order infinitesimals, which in fact gives: $`\gamma 1+(\gamma _0M_{\mathrm{ej}}/m)^2/2`$, then we have $`vR^3`$.<sup>15</sup> This is not consistent with the Sedov solution!

The problem is serious: (i) It means that the reliability of Eq.(2) is questionable, although it does correctly reproduce the major features for radiative fireballs and even for adiabatic fireballs in the ultra-relativistic limit. (ii) In the non-relativistic phase of the expansion, the fireball is more likely to be adiabatic rather than highly radiative. However, it is just in this condition that the conventional model fails. So any calculation made according to Eq.(2) will lead to serious deviations in the light curves during the non-relativistic phase.

Third, for simplicity, it is usually assumed that $`ϵ`$ is a constant during the expansion. But in realistic case this is not true. The fireball is expected to be highly radiative ($`ϵ=1`$) at first, due to significant synchrotron radiation. In only one or two days, it will evolve to an adiabatic one ($`ϵ=0`$) gradually. So $`ϵ`$ should evolve with time.<sup>16</sup>

Below, we will construct a new model that is no longer subject to the aforementioned problems.

In the fixed frame, since the total kinetic energy of the fireball is $`E_\mathrm{K}=(\gamma 1)(M_{\mathrm{ej}}+m)c^2+(1ϵ)\gamma U`$,<sup>17</sup> and the radiated thermal energy is $`ϵ\gamma (\gamma 1)dmc^2`$,<sup>11</sup> we have:

$$d[(\gamma 1)(M_{\mathrm{ej}}+m)c^2+(1ϵ)\gamma U]=ϵ\gamma (\gamma 1)c^2dm.$$
(7)   
For the item $`U`$, it is usually assumed: $`dU=c^2(\gamma 1)dm`$.<sup>17</sup> Eq.(2) has been derived just in this way. However, the jump conditions <sup>11</sup> at the forward shock imply that $`U=(\gamma 1)mc^2`$, so we suggest that the correct expression for $`dU`$ should be: $`dU=d[(\gamma 1)mc^2]=(\gamma 1)dmc^2+mc^2d\gamma `$. Here we simply use $`U=(\gamma 1)mc^2`$ and substitute it into Eq.(7), then it is easy to get:<sup>15</sup>

$$\frac{d\gamma }{dm}=\frac{\gamma ^21}{M_{\mathrm{ej}}+ϵm+2(1ϵ)\gamma m}.$$
(8)   
In the highly radiative case ($`ϵ=1`$), Eq.(8) reduces to Eq.(4) exactly. While in the adiabatic case ($`ϵ=0`$), Eq.(8) reduces to:

$$\frac{d\gamma }{dm}=\frac{\gamma ^21}{M_{\mathrm{ej}}+2\gamma m}.$$
(9)   
The analytic solution is:

$$(\gamma 1)M_{\mathrm{ej}}c^2+(\gamma ^21)mc^2E_0.$$
(10)   
Then in the ultra-relativistic limit, we get the familiar relation of $`\gamma R^{3/2}`$, and in the non-relativistic limit, we get $`vR^{3/2}`$ as required by the Sedov solution. From these analysises, we believe that Eq.(8) is really correct for both radiative and adiabatic fireballs, and in both ultra-relativistic and non-relativistic phase.

In realistic fireballs, $`ϵ`$ is a variable dependent on the ratio of synchrotron-radiation-induced to expansion-induced loss rate of energy.<sup>16</sup> As usual, we assume that in the co-moving frame the magnetic field energy density is a fraction $`\xi _\mathrm{B}^2`$ of the energy density $`e^{}`$, $`B^2/(8\pi )=\xi _\mathrm{B}^2e^{}`$, and that the electron carries a fraction $`\xi _\mathrm{e}`$ of the energy, $`\gamma _{\mathrm{e},\mathrm{min}}=\xi _\mathrm{e}(\gamma 1)m_\mathrm{p}/m_\mathrm{e}+1`$, where $`m_\mathrm{p}`$ and $`m_\mathrm{e}`$ are proton and electron masses, respectively. The co-moving frame expansion time is $`t_{\mathrm{ex}}^{}=R/(\gamma c)`$, and the synchrotron cooling time is $`t_{\mathrm{syn}}^{}=6\pi m_\mathrm{e}c/(\sigma _\mathrm{T}B^2\gamma _{\mathrm{min},\mathrm{e}})`$, where $`\sigma _\mathrm{T}`$ is the Thompson cross section. Then we have:<sup>16</sup>

$$ϵ=\xi _\mathrm{e}\frac{t_{\mathrm{syn}}^1}{t_{\mathrm{syn}}^1+t_{\mathrm{ex}}^1}.$$
(11)   
We have evaluated Eqs.(8) and (11) numerically, bearing in mind that:<sup>18</sup>

$$dm=4\pi R^2nm_\mathrm{p}dR,$$
(12)   
$$dR=v\gamma (\gamma +\sqrt{\gamma ^21})dt.$$
(13)   
We take $`E_0=10^{52}`$ erg, $`n=1`$ cm<sup>-3</sup>, $`M_{\mathrm{ej}}=2\times 10^5`$ M. Figs.(1) and (2) illustrate the evolution of $`\gamma `$ and $`R`$, where full, dotted, and dashed lines correspond to constant $`ϵ`$ values of 0, 0.5, and 1 respectively. Dash-dotted lines are plot by allowing $`ϵ`$ to vary according to Eq.(11). It is clearly shown that our new model overcomes the shortcomings of Eq.(2). For example, for highly radiative expansion, the dashed lines in these figures approximately satisfy $`\gamma t^{3/7}`$, $`Rt^{1/7}`$, $`\gamma R^3`$ when $`\gamma 1`$, and $`vt^{3/4}`$, $`Rt^{1/4}`$, $`vR^3`$ when $`\gamma 1`$. While for adiabatic expansion, the full lines satisfy $`\gamma t^{3/8}`$, $`Rt^{1/4}`$, $`\gamma R^{3/2}`$ when $`\gamma 1`$, and satisfy $`vt^{3/5}`$, $`Rt^{2/5}`$, $`vR^{3/2}`$ when $`\gamma 1`$.

In order to compare with observations, we have also calculated the synchrotron radiation from the shocked ISM. Fig.(3) illustrates R band afterglows from the realistic fireball. We see that after entering the non-relativistic phase, the light curve becomes steeper only slightly, consistent with the prediction made by Wijers et al..<sup>7</sup> In contrast, Eq.(2) generally leads to a much sharper decline.<sup>18</sup> In our new model optical afterglows from GRB 970228 are generally well fitted.

To conclude, current researches are mainly concentrated on the ultra-relativistic phase of the expansion of GRB remnants. The popular dynamic model is in fact incorrect for adiabatic fireballs during the non-relativistic phase. This is completely unnoticed in the literature. We have revised the model. Our new model has been shown to be correct in both ultra-relativistic and non-relativistic phase. The revision is of great importance, taking account of the following facts: (i) Optical afterglows lasting for more than 100 – 200 d have been observed from some GRBs, the advent of the non-relativistic phase seems inevitable. (ii) Beaming effects also lead to a steepening in the optical light curve, non-relativistic effects should be considered carefully to tell whether GRB ejecta are beamed or not, which is crucial in understanding the GRB origin. (iii) HI supershells might be highly evolved GRB remnants,<sup>19,20</sup> to address this question in detail, we should deal with non-relativistic blastwaves. Additionally we suggest that at very late stages, GRB remnants might become highly radiative again, just in the same way that supernova remnants do.<sup>14</sup> This might occur when the bulk velocity is just several tens kilometers per second.

Figure Caption