Next Article in Journal
30 Years of Multifrequency Quasar Variability: A Personal Journey
Next Article in Special Issue
Stable Up-Down Quark Matter Nuggets, Quark Star Crusts, and a New Family of White Dwarfs
Previous Article in Journal
Simulations of Neutrino and Gamma-Ray Production from Relativistic Black-Hole Microquasar Jets
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Low-Energy Spectral Index of Gamma-ray Burst Prompt Emission from Internal Shocks

1
Department of Astronomy, School of Physics, Huazhong University of Science and Technology, Wuhan 430074, China
2
Department of Astronomy, University of Science and Technology of China, Hefei 230026, China
3
School of Astronomy and Space Science, Nanjing University, Nanjing 210023, China
*
Author to whom correspondence should be addressed.
Galaxies 2021, 9(3), 68; https://doi.org/10.3390/galaxies9030068
Submission received: 30 July 2021 / Revised: 9 September 2021 / Accepted: 13 September 2021 / Published: 15 September 2021
(This article belongs to the Special Issue Gamma-Ray Burst Science in 2030)

Abstract

:
The prompt emission of most gamma-ray bursts (GRBs) typically exhibits a non-thermal Band component. The synchrotron radiation in the popular internal shock model is generally put forward to explain such a non-thermal component. However, the low-energy photon index α 1.5 predicted by the synchrotron radiation is inconsistent with the observed value α 1 . Here, we investigate the evolution of a magnetic field during propagation of internal shocks within an ultrarelativistic outflow, and revisit the fast cooling of shock-accelerated electrons via synchrotron radiation for this evolutional magnetic field. We find that the magnetic field is first nearly constant and then decays as B t 1 , which leads to a reasonable range of the low-energy photon index, 3 / 2 < α < 2 / 3 . In addition, if a rising electron injection rate during a GRB is introduced, we find that α reaches 2 / 3 more easily. We thus fit the prompt emission spectra of GRB 080916c and GRB 080825c.

1. Introduction

The prompt radiation mechanism of gamma-ray bursts (GRBs) is still being debated, even though the prompt spectra can usually be fitted well by the Band function [1], which suggests a smoothly jointed broken power law with the low-energy photon index α 1 , the high energy photon index β 2.2 and the peak energy E p 250 keV [2,3]. Currently, neither the possible one-temperature thermal emission from an ultrarelativistic fireball, nor the single synchrotron radiation from shock-accelerated electrons within this fireball, provide an explanation for such a low-energy photon index (see [4,5] for a review).
Generally, there are two mechanisms that explain the low energy photon index α 1 of the GRB prompt emission. The first mechanism is the Comptonized quasi-thermal emission from the photosphere of an ultrarelativistic outflow [6,7,8,9,10,11,12,13,14,15,16,17,18,19,20]. The second mechanism is synchrotron and/or synchrotron self-Compton (SSC) emission in the optically thin region. For fast-cooling synchrotron radiation in the internal shock model, possible solutions include invoking a small-scale rapidly decaying magnetic field [21], a decaying magnetic field with a power-law index in a relativistically-expanding outflow [22,23,24], a decaying magnetic field in a post-shock region [25], Klein–Nishina (KN) cooling [26,27], an adjustable synchrotron self-absorption frequency [28,29], or the acceleration process [30] and other evolutional model parameters [31]. Alternatively, slow cooling was introduced to understand the low-energy photon index [32]. In addition to the internal shock model, the other energy dissipation mechanisms, such as the ICMART model [33], were proposed to solve the low-energy spectral index issue. In some models (e.g., [34]), the observed prompt emission of GRBs is understood to be dominated by the SSC emission, while the synchrotron radiation is in much lower energy bands.
The fast cooling synchrotron radiation in the internal shock model is generally considered to be a straightforward and leading mechanism to explain the GRB prompt emission spectra, and the most important issue in this model is to explain the low-energy spectral index. An underlying assumption in the traditional synchrotron internal shock model is to calculate the electron cooling without considering the evolution of the magnetic field. In other words, the magnetic field is treated as a constant and its effect in the continuity equation of electrons is usually ignored (e.g., [35]). In the fast cooling case, the predicted low-energy photon index α 3 / 2 is much softer than observed. In this paper, we try to alleviate this problem. We calculate the magnetic field in the realistic internal shock model during a collision of two relativistic thick shells and obtain an evolutional form of the magnetic field, B constant before the time δ t that is nearly equal to the ejection time interval of the two shells, and B t 1 after the time δ t . We consider the cooling of electrons accelerated by internal shocks for this evolutional magnetic field, and find the resulting spectral index α 3 / 2 for B constant and α 2 / 3 for B t 1 , by adopting a cooling method similar to that in Ref. [22]. Actually, these two cases may coexist, and the outflow may undergo the first case and then the second case, so theoretically the actual index α will range from 3 / 2 to 2 / 3 . Furthermore, below the peak energy E p there is a gradual process, so that α is only close to 1 . In addition, we consider a rising electron injection rate, leading to a larger α , slightly smaller than 2 / 3 .
This paper is organized as follows. We calculate the dynamics of a collision between two thick shells in Section 2. In Section 3, we investigate the electron cooling and its synchrotron radiation with an evolutional magnetic field and a rising electron injection rate. In the final section, discussions and conclusions are given.

2. Dynamics of Two-Shell Collision

In the popular internal shock model, an ultrarelativistic fireball consisting of a series of shells with different Lorentz factors can produce prompt emission through collisions among these shells. For the dynamics of two-shell collision, we adopt the same approach as the one in [36]. In order to present one GRB prompt emission component (with duration ∼ few seconds), we here consider two thick shell–shell collision to produce a consistent GRB pulse with a duration of few seconds (i.e., the slow pulse) and the fast pulses with a duration of 0.01   s in GRBs may be caused by the density fluctuation of the shell. Under this assumption, a prior slow thick shell A with bulk lorentz factor γ A and kinetic luminosity L k , A , and a posterior fast thick shell B with bulk lorentz factor γ B (where γ B > γ A 1 ) and kinetic luminosity L k , B is adopted. The collision of the two shells begins at radius [36]
R c o l = β B c β A Δ t i n t ( β B β A ) 2 γ A 2 c Δ t i n t 1 ( γ A / γ B ) 2 2 γ A 2 c δ t 5.4 × 10 14 γ A , 2.5 2 δ t , 1 cm ,
where Δ t i n t is the time interval between the two thick shells, and δ t Δ t i n t / [ 1 ( γ A / γ B ) 2 ] is a redefined time interval. For γ A γ B , δ t Δ t i n t . The conventional expression Q , m = Q / 10 m is used. During the collision, there are four regions separated by internal forward-reverse shocks: (1) the unshocked shell A; (2) the shocked shell A; (3) the shocked shell B; and (4) the unshocked shell B, where regions 2 and 3 are separated by a contact discontinuity.
The particle number density of a shell measured in its comoving frame can be calculated as [37]:
n i = L k , i 4 π R 2 γ i 2 m p c 3 ,
where R is the radius of the shell and subscript i can be taken as A or B. As in the literature [36,37,38,39,40], we derive the dynamics of internal forward-reverse shocks. In order to get a high prompt emission luminosity, it is reasonable to assume γ A γ B and L k , A = L k , B L k . Assuming that γ 1 , γ 2 , γ 3 , and γ 4 are Lorentz factors of regions 1, 2, 3 and 4 respectively, we have γ 1 = γ A , γ 4 = γ B , and n 1 n 4 . If a fast shell with low particle number density catches up with a slow shell with high particle number density and then they collide with each other, a Newtonian forward shock (NFS) and a relativistic reverse shock (RRS) may be generated [36,37]. So we can obtain γ 1 γ 2 = γ 3 = γ γ 4 . Then, according to the jump conditions between the two sides of a shock [41], the comoving internal energy densities of the two shocked regions can be calculated following e 2 = ( γ 21 1 ) ( 4 γ 21 + 3 ) n 1 m p c 2 and e 3 = ( γ 34 1 ) ( 4 γ 34 + 3 ) n 4 m p c 2 , where γ 21 = 1 2 ( γ 1 / γ 2 + γ 2 / γ 1 ) and γ 34 = 1 2 ( γ 3 / γ 4 + γ 4 / γ 3 ) are the Lorentz factors of region 2 relative to the unshocked region 1, and region 3 relative to region 4, respectively. It is required that e 2 = e 3 because of the mechanical equilibrium. We have [36,37]
( γ 21 1 ) ( 4 γ 21 + 3 ) ( γ 34 1 ) ( 4 γ 34 + 3 ) = n 4 n 1 = γ 1 γ 4 2 f .
Two relative Lorentz factors can be calculated as γ 21 f γ 4 2 7 γ 1 2 + 1 = 8 7 , and γ 34 = γ 4 2 γ 1 1 . Assuming that t is the observed shell–shell interaction time since the prompt flare onset, the radius of the system during the collision can be written as
R = R c o l + 2 γ 2 c t 2 γ 1 2 c ( t + δ t ) .
During the propagation of the shocks and before the shock crossing time, the instantaneous electron injection numbers (in d t ) in regions 2 and 3 can be calculated as follows [38]:
d N e , 2 = 8 π R 2 n 1 ( γ 21 β 21 / γ β ) γ 2 c d t
and
d N e , 3 = 8 π R 2 n 4 ( γ 34 β 34 / γ β ) γ 2 c d t ,
respectively.

3. Synchrotron Radiation with a Decaying Magnetic Field and a Variable Electron Injection Rate

3.1. Synchrotron Radiation with a Decaying Magnetic Field

As usual, we assume that fractions ϵ B and ϵ e of the internal energy density in a GRB shock are converted into the energy densities of the magnetic field and electrons, respectively. Thus, using B i = ( 8 π ϵ B e i ) 1 / 2 for i = 2 or 3, we can calculate the strength of the magnetic field before the shock crossing time t c r s by
B 2 = B 3 = ϵ B L k 2 γ 1 6 c 3 ( t + δ t ) 2 1 / 2 ,
and find that the change of the magnetic field before δ t can be ignored (i.e., B i constant ), but after δ t the magnetic field B i decreases linearly with time t (i.e., B i t 1 ). Actually, the evolution of the magnetic field is caused by the expansion of the shocked regions, which is presented in Figure 1. After the shock crossing time t c r s (here, t c r s is comparable with the peak time of the slow pulse in GRBs), the spreading of the hot materials into the vacuum cannot be ignored and the merged shell undergoes an adiabatic cooling. During this phase, the volume of the merged shell is assumed to expand as V i R s , where s is a free parameter and its value is taken to be from 2 to 3. As a result, the particle number density would decrease as n i V i 1 R s , the internal energy density as e i V i 4 / 3 R 4 s / 3 , and the magnetic field strength as B i ( e i ) 1 / 2 R 2 s / 3 t 2 s / 3 . Because no additional shock-accelerated electrons are injected after the shock crossing time t c r s , we only study the prompt emission before t c r s in the remaining part of this paper. What we want to point out is that the redefined time interval δ t is not equal to the shock crossing time ( t c r s ), the latter one is dependent on the thickness of the shells. In this paper, the two shells are must be thick enough so that t c r s δ t .
The electrons accelerated by the shocks are assumed to have a power-law energy distribution, d N e , i / d γ e , i γ e p for γ e , i γ e , m , i , where γ e , m , i is the minimum Lorentz factor of the accelerated electrons. The following electron cooling discussion is not based on the conventional synchrotron and SSC cooling, which always give us the electron distribution, d N e d γ e γ e 2 for γ e < γ e , m and d N e d γ e γ e p 1 for γ e > γ e , m in the fast cooling case, d N e d γ e γ e p for γ e < γ e , m and d N e d γ e γ e p 1 for γ e > γ e , m in the slow cooling case. These electron distributions do not take into account the evolution of the magnetic field. Ref. [22] discussed the electron distribution affected by a decaying magnetic field based on B r b , where r is the fireball radius and b is the magnetic field decaying index. They considered the electron distribution of a group of plasma in a magnetic field with an arbitrary decaying index b, which is called a “toy box model". Here we consider a more physical process, internal shocks, which generate an evolutional magnetic field and a consistent spectrum with the observed Band spectral shape.
In the comoving frame, the evolution of the Lorentz factor of an electron via synchrotron and SSC cooling and adiabatic cooling can be described by [22]
d d t 1 γ e = σ T ( 1 + Y i ) 6 π m e c B i 2 1 3 1 γ e d ln n i d t ,
where Y i [ ( 4 η i ϵ e / ϵ B + 1 ) 1 / 2 1 ] / 2 is the Compton parameter, which is defined by the ratio of the IC to synchrotron luminosity, with η i = min [ 1 , ( γ e , c , i / γ e , m , i ) 2 p ] [42]. γ e , c , i is the cooling Lorentz factor and the comoving time t = 2 γ t .
The minimum Lorentz factor of the accelerated electrons is γ e , m , i = m p m e ( p 2 p 1 ) ϵ e ( γ rel 1 ) (where γ rel = γ 21 or γ 34 in region 2 or 3), so it can be written as:
γ e , m , 3 1.0 × 10 4 g p ϵ e , 1 / 2 γ 4 , 4.5 γ 1 , 2.5 1 ,
γ e , m , 2 30 g p ϵ e , 1 / 2 ,
where g p = 3 ( p 2 ) / ( p 1 ) . Moreover, the cooling Lorentz factor γ e , c , i = 6 π m e c / ( y i σ T B 3 2 γ t ) , can be written as
γ e , c , 3 = γ e , c , 2 3.4 × 10 2 y i , 0 1 ϵ B , 1 / 2 1 L k , 51 1 γ 1 , 2.5 5 ( t + δ t ) , 0 2 t , 0 ,
where y i = 1 + Y i is the ratio of the total luminosity to synchrotron luminosity.
From the electron injection rate based on Equations (5) and (6), one can obtain the injected electrons number between t and t + d t . Assuming the original electron injection distribution d N e , i / d γ e , i γ e p for γ e , i γ e , m , i , the injected electrons number between t and t + d t and between γ e and γ e + d γ e can be derived. So we cut the injected electrons into small pieces in the time space t and the energy space γ e . At the beginning, time t = 0 , a number of electrons d N will be injected into the shocked region in a time interval d t and will be cooled in the initial magnetic field, so one can obtain the change of electron Lorentz factor Δ γ e , 1 based on Equation (8) for the electrons between γ e and γ e + d γ e . In the next time interval d t , these electrons with the Lorentz factor between γ e + Δ γ e , 1 and γ e + d γ e + Δ γ e , 1 will be cooled in the instantaneous magnetic field based on the evolutional magnetic field in Equation (7), and one can obtain another Δ γ e , 2 ( Δ γ e , 2 Δ γ e , 1 ). At the same time, another group electrons are injected and cooled in this instantaneous magnetic field. These processes are continuous before 2 γ t c r s . The shocked electrons are injected as time and all the electrons are cooled in the instantaneous magnetic field. We sum all electrons at time t in the energy space, obtain the electron distributions at time t and present them in Figure 2 ( t o b s = t / 2 γ ). As shown in Figure 2, when t < δ t , the magnetic field does not change significantly (see Figure 1), the electron distribution in the fast cooling case, d N e d γ e γ e 2 for γ e < γ e , m , and d N e d γ e γ e p 1 for γ e > γ e , m , are expected. However, when t > δ t , the electron distribution below γ e , m would be flattened because of the decaying magnetic field. Due to the magnetic field decay, the electrons injected at later times would cool more slowly than the electrons injected at early times (here, all times are before t c r s ). In other words, the cooling efficiency would become smaller due to the decaying magnetic field, which induces more electrons accumulating at γ e , m than in the invariable magnetic field case. When t δ t and t < t c r s , the electron spectral index for γ e < γ e , m is even flattened to zero.
In order to find these new results for t δ t and t < t c r s , we can evaluate the continuity equation of electrons in energy space, t ( d N e , γ e / d γ e ) + γ e [ γ e ˙ ( d N e , γ e / d γ e ) ] = Q ( γ e , t ) , where d N e , γ e / d γ e is the instantaneous electron spectrum at time t , and Q ( γ e , t ) = Q 0 ( t ) ( γ e / γ e , m ) p is the electron injection distribution accelerated by shocks above the minimum injection Lorentz factor γ e , m . By ignoring the inconsequential adiabatic cooling term, we can get d d t ( 1 γ e ) ( 1 + Y i ) B i 2 t 2 , where Y i is assumed to be a constant before the shock crossing time t c r s in the fast cooling case. Then, we can obtain γ e t , and thus γ ˙ e γ e 2 t 2 γ e 0 . For γ e , c < γ e < γ e , m , Q ( γ e , t ) = 0 , to obtain the final and quasi-steady electron spectral shape at the arbitrary time t , by considering a quasi-steady-state system ( / t = 0 ), we can easily find d N e , γ e / d γ e γ e 0 below γ e , m .
Next, the four characteristic frequencies in regions 2 and 3 that can be calculated from ν = q e 2 π m e c B γ e 2 γ are derived as [36]
h ν m , 2 2.1 × 10 4 g p 2 ϵ e , 1 / 2 2 ϵ B , 1 / 2 1 / 2 L k , 51 1 / 2 γ 1 , 2.5 2 ( t + δ t ) , 0 1 keV ,
h ν m , 3 26 g p 2 ϵ e , 1 / 2 2 ϵ B , 1 / 2 1 / 2 L k , 51 1 / 2 γ 4 , 4.5 2 γ 1 , 2.5 4 ( t + δ t ) , 0 1 keV ,
and
h ν c , 2 = h ν c , 3 3.6 × 10 2 y , 0 2 ϵ B , 1 / 2 3 / 2 L k , 51 3 / 2 γ 1 , 2.5 8 ( t + δ t ) , 0 3 t , 0 2 keV .
Here, if γ 1 = 100 and γ 4 = 10,000 , we obtain h ν m , 3 186 keV at time t = 1 s , which is approximatively equal to the typical value of E p of the GRB prompt emission.
We also present the spectrum of region 3 in the top panel of Figure 3 based on the electron distribution shown in Figure 2. However, we do not present the spectrum of region 2 because, from NFS, (1) its photon peak frequency is much smaller than the typical GRB prompt emission E p , (2) the radiation efficiency can not be high enough as a result of slow cooling, and (3) the flux of region 2 is much lower than that of region 3. The last reason can be evaluated from [35,37]
F ν , max , i < N e , i 4 π D L 2 m e c 2 σ T 3 q e B i γ ,
where D L is the luminosity distance of the burst and N e , i is the total number of injected electrons until the time t. Since a portion of the electrons have cooled to a much smaller value than γ b r e a k (where the break Lorentz factor γ b r e a k of an electron distribution, γ b r e a k = γ e , m for fast cooling, and γ b r e a k = γ e , c for slow cooling), the actual number of electrons near γ b r e a k is less than N e , i and thus the actual F ν , max , i is smaller than the right term of inequality Equation (15). So, we can obtain
ν m , 3 F ν , max , 3 < 1.5 × 10 5 g p 2 ϵ e , 1 / 2 2 ϵ B , 1 / 2 L k , 51 2 γ 4 , 4 γ 1 , 2 6 × t , 0 ( t + δ t ) , 0 2 D L , 28 2 keV cm 2 s 1 ,
and
ν c , 2 F ν , max , 2 < 8.7 × 10 1 y , 0 2 ϵ B 1 / 2 1 γ 1 , 2 5 × ( t + δ t ) , 0 2 t , 0 D L , 28 2 keV cm 2 s 1 ,
where γ 1 = 100 and γ 4 = 10,000 are taken.
From Figure 3, we can see that for t < δ t , because of a constant magnetic field, the spectral slope of ν F ν is 1 / 2 as described by [35]. However, when t > δ t , the spectral slope will deviate from 1 / 2 and become a larger value (even 4 / 3 ). If the electron index ( d N e / d γ e γ e u ) is u, the F ν slope of synchrotron radiation ( F ν ν w ) would be w = ( u 1 ) / 2 and the photon spectral index (defined as d N γ / d E γ = E γ α , where E γ is the photon energy, and N γ is the photon number flux) would be α = ( w + 1 ) . Due to the decaying magnetic field, u tends to be zero, and thus w = 1 / 2 and α = 1 / 2 . However, when α > 2 / 3 , because of the overlying effect, the low energy photon index of the electrons with γ e , m is 2 / 3 and will cover the emission of electrons with smaller Lorentz factors. So, due to the effect of the low-energy radiation tail of electrons with Lorentz factor γ e , m , α is at most equal to 2 / 3 , and we can get 3 / 2 < α < 2 / 3 . This is consistent with the observations [2,3], which suggest α 1 .

3.2. The Effect of the Variable Electron Injection Rate

Although theoretically the low-energy photon spectral index α can reach 2 / 3 caused by a decaying magnetic field, since this is a gradual process, α would be softer than 2 / 3 for E E p , which can be seen from Figure 3. In fact, α 1 can be fitted easily, but fitting α slightly smaller than 2 / 3 is difficult. We here consider a variable electron injection rate, which could induce α 2 / 3 . The variable electron injection rate may be suggested by that the actual GRB shell is not homogeneous and presents a density profile, for example, a Gaussian density profile, inducing a rising electron injection rate. Nonetheless, we do not know its growing method clearly. Ref. [22] discussed this effect in their “toy box model”, and suggested, because of a rising electron injection rate, α goes from 0.82 to 1.03 , which is dependent on the growing power-law index q (where the injection rate t q with q = 1 , 2 or 3). Here we adopt similar expressions of the rising electron injection rate,
d N e , 2 = 8 π R 2 n 1 ( γ 21 β 21 / γ β ) γ 2 c d t × t t 0 q
and
d N e , 3 = 8 π R 2 n 4 ( γ 34 β 34 / γ β ) γ 2 c d t × t t 0 q ,
where the factor ( t t 0 ) is to maintain the same electron injection number in the interval t 0 as that for the constant injection rate q = 0 .
We show the electron distribution for a rising electron injection rate in the top panel of Figure 4. The rising electron injection would increase the electrons injected later, which would cool in a weaker magnetic field and pile up at γ e , m . This can result in a harder electron distribution and a relative spectrum. The slopes of the spectra are presented in the bottom panel of Figure 4. We can see that the slopes of the spectra tend to reach 4 / 3 more easily than in the constant electron injection case. In addition, a larger q would generate a harder low-energy spectral index.

4. Application to the Actual GRB Spectra

In order to compare with the actual GRB spectrum, we select the broad band spectrum of GRB 080916c in the interval “b” detected by Gamma-ray Burst Monitor (GBM) and the Large Area Telescope (LAT) aboard the Fermi satellite (see Ref. [43]), from 3.58 s to 7.68 s since the lightcurve during this period is presented as a single and pure pulse. Moreover, its low energy photon index is close to the typical value of the low energy photon index of the GRB, that is, α 1 , harder than the expectation of synchrotron fast cooling ( 1.5 ). The “b” spectrum of GRB 080916c can be well fitted by the Band function with the low energy photon index of α = 1.02 ± 0.02 , the high energy photon spectral slope β = 2.21 ± 0.03 and the peak energy E p = 1170 ± 140 keV [43]. Since the observational data can be well fitted by this Band function with a very small error range, the Band function is precise enough to represent the actual GRB emission. We select some representative points (black points in Figure 5) in this Band function to present the tendency of the actual GRB emission. In addition, more black points around the peak energy in the figure are taken to present the gradual change in behavior there. In Figure 5, by using a time-averaged energy spectrum from t = 0 s to t = 3 s , the emission of GRB 080916c can be fitted well in our model with the proper parameters, which have been listed in Table 1.
We select the single pulse spectrum of GRB 080825c in the interval “a” detected by Fermi GBM and LAT (see Ref. [44]), from 0.0 s to 2.7 s , which has a harder photon index, α 0.76 . The “a” spectrum of GRB 080825c can be well fitted by the Band function with the low energy photon index α = 0.76 ± 0.05 , the high energy photon index β = 2 . 54 0.17 + 0.11 and the peak energy E p = 291 22 + 25 keV [44]. Such a hard photon index could not be approached easily for a constant electron injection rate, that is, q = 0 , so we consider a rising electron injection rate as suggested in Section 3.2. Some representative points (black points in Figure 6) in this Band function are selected to present the tendency of the actual GRB emission as the same as the treatment for the GRB 080916c. The observational spectrum can be reproduced well in Figure 6 phenomenally by using a time-averaged energy spectrum from t = 0 s to t = 3 s with an index of rising electron injection rate q = 2 and other reasonable parameters (all parameters are listed in Table 1).
The main parameters to effect the final spectrum are listed in Table 1. The dependence of the break energy of the spectrum on the listed parameters could be found in Equation (13) and for the magnitude of peak flux the dependence could be derived roughly in Equation (16). During the model fitting, for simplification, the energy equipartition factors for electrons and the magnetic field, that is, ϵ e and ϵ B , and γ 4 are fixed, and then the Lorentz factor γ 1 and kinetic luminosity L k are adjusted to match the observational peak energy and the peak flux. The electron injection index p is determined by the observational high-energy photon index since the relation between them is β ( p 2 ) / 2 , 1 suggested by the synchrotron radiation. The shock cross time is comparably adopted with the typical duration of the slow pulse of the GRB, namely, 3 s . Different values of δ t could affect the evolutional form of the magnetic field (as shown in Figure 1) and adjust the weight of the cooling in a constant magnetic field and the cooling in a decaying magnetic field. In other words, a smaller δ t could make it so that the electron synchroton cooling mainly takes place in a decaying magnetic field and leads the photon index to be harder, while for a larger δ t , the electrons are mainly cooling in a constant magnetic field and generating a photon index close to 1.5 . As a result, for GRB 080916c with a photon index 1 , a relatively small δ t = 0.1 s is adopted. A harder photon index 0.76 for GRB 080825c and a rising electron injection rate with an index q = 2 are taken into account as suggested in Section 3.2. Therefore, a certain range of a low-energy photon index from 3 / 2 to 2 / 3 could be approached through the adjustment of δ t and the index of the electron injection rate q. However, for a low-energy photon index harder than 2 / 3 , this model would become invalid.

5. Discussions and Conclusions

Currently alleviating the tension between the expectation of synchrotron and observations in the GRB prompt regime is a more and more important issue. Two classes of model have been proposed to explain the low-energy photon index of GRB prompt emission, Comptonized quasi-thermal emission from the photosphere within a relativistic outflow and synchrotron and/or SSC emission in the optically thin region. These models can experience difficulties. For Comptonized quasi-thermal emission, the most significant effect to obtain α 1 is the equal arrival time effect in this model, which is relevant to the end time of central engine activity, but may not be applicable during the prompt emission phase when a continuous wind is ejected from the central engine [45]. Synchrotron slow cooling in internal shocks may not provide a high radiative efficiency, and a dominant SSC component usually predicts an even more dominant 2nd-order SSC component, which significantly exceeds the total energy budget of GRBs [46,47]. Thus, some evolutional parameters, such as the magnetic field, the fraction of the accelerated electrons, and the energy equipartition factors, were suggested to explain the low-energy index.
In this paper, we have considered a straightforward model, that is, the fast cooling synchrotron radiation in internal shocks. We obtain the magnetic field evolutional form in a practical shell–shell collision, B constant before δ t and B t 1 after this time, and recalculate the electron distribution for this evolutional magnetic field. When t < δ t , the magnetic field is nearly constant, and the fraction of cooling electrons in the invariable magnetic field is high enough so that d N e d γ e γ e 2 for γ e < γ e , m is expected. However, when t δ t but t < t c r s , the fraction of cooling electrons in the evolutional magnetic field is higher than in the invariable magnetic field, so that d N e d γ e will be gradually proportional to γ e 0 . 2 γ 2 c δ t and 2 γ 2 c t indicate roughly the collision radius and the propagation distance of a relativistic outflow after the collision takes place but before the shock crossing time t c r s , respectively. In other words, if the propagation distance of the outflow is smaller than the collision radius before the shock crossing time, the magnetic field can be treated as a constant and it is not necessary to consider the evolution of the magnetic field when calculating the electron cooling. However, if the propagating distance of the outflow is larger than the collision radius before the shock crossing time, we have to consider the evolution of the magnetic field and can obtain a different electron distribution. Actually, the outflow may undergo the first case and then the second case, so we can obtain a reasonable range of the low-energy photon index α , from 3 / 2 to 2 / 3 theoretically. Since d N e d γ e proportional to γ e 0 is a gradual process below E p , it is usually difficult to get α to be exactly equal to 2 / 3 , but this index is only slightly smaller than 2 / 3 . Moreover, we also consider a rising electron injection rate, which may exacerbate this situation, inducing α to be closer to 2 / 3 .
Ref. [22] considered a decaying magnetic field varying with the distance from the central engine to explore the range of a low energy photon index in the GRB prompt regime. They discussed the radiation spectra of a cloud of plasma in a decaying magnetic field with an arbitrary decaying index b for a simplified model, which is called a “toy box model". Different from their work, we adopt a more physical case for the internal shock by considering the collision of two shells and inducing the decaying form of the magnetic field. As a result, a time-dependent magnetic field is derived (as shown in Figure 1). In fact, a time-dependent magnetic field could be translated to a distance-dependent form due to the propagation of relativistic outflow. For the evolutional magnetic field form obtained from the practical internal shock, we study the influence on the spectral index. In addition to the detailed treatment of shell–shell collision, the kinetic luminosity and the energy equipartition parameters, ϵ B and ϵ e are taken into account to obtain the radiation spectra, comparing them with the actual GRB spectra for GRB 080916c and GRB 080825c.
In our model, in order to obtain the high prompt emission luminosity, we assume that γ 4 γ 1 . This assumption is reasonable. This is because estimates based on four methods by Ref. [48] show that the mean observed value of the bulk Lorentz factors of GRB outflows is a few hundred, corresponding to γ 2 = γ 3 γ 1 100 in our model. Furthermore, within the framework of the collapsar model, a prior relativistic jet-like shell (e.g., shell A) first has to propagate through the envelope of a massive star and clean up almost all of the baryons along the propagation direction of this shell, leaving behind a clean passage for a posterior jet-like shell (e.g., shell B). This, therefore, leads to a reasonable possibility that the Lorentz factor of shell B is much greater than that of shell A.
Usually, we have γ 1 100 , so γ 4 10 4 γ 1 2 is a universal relationship to obtain the high prompt emission luminosity. Ref. [37] also mentioned that, when γ 4 γ 1 2 , the highest luminosity from internal shocks is expected. In fact, this assumption is not a special case. When collisions among a series of shells with different Lorentz factors occur, the highest luminosity from one collision will cover the others. In other words, we always see the brightest. According to Equation (13), if deeming that γ 4 γ 1 2 does not vary significantly among bursts, we can easily obtain the so-called “Yonetoku Relation”, E p L i s o 1 / 2 [49], and the “Amati Relation”, E p E i s o 1 / 2 [50]. However, this model is also confronted with some issues, for example, the spectrum is somewhat broad near E p in contrast to the observed data or the Band function [51], which can be seen in Figure 6.
It is important that we should beware of the empirical Band function. The thermal components and more spectral structures are found in the prompt regimes of some GRBs, which deviate from the so-called Band function [52,53]. The thermal emission generated by the photosphere is a natural prediction of the generic fireball scenario. The relative strength of thermal emission and non-thermal emission may depend on the various environments [54,55]. Ref. [53] also claimed that the GRB spectra below the peak energy may present an extra break energy around a few keV, inducing a consistent spectral shape with expectation from the classical synchrotron radiation. More spectral structures of GRBs may make the simple Band function become invalid, and result in an incorrect low-energy spectral index if one forcibly fits them using a Band function. Although our model can present a consistent low-energy spectral index with observations in a certain range, due to the complexities of GRB prompt spectra, more detailed studies are needed.

Author Contributions

All authors have contributed equally to the preparation of this manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Key Research and Development Program of China (grant No. 2017YFA0402600), the National SKA Program of China (grant No. 2020SKA0120300), the National Natural Science Foundation of China under grants No. 11833003, 12003007 and the Fundamental Research Funds for the Central Universities (No. 2020kfyXJJS039).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We thank the anonymous referees for the helpful suggestions and comments.

Conflicts of Interest

The author declares no conflict of interest.

Note

1
The high-energy electron distribution above the break electron energy is d N e / d γ e γ e ( p 1 ) .

References

  1. Band, D.; Matteson, J.; Ford, L.; Schaefer, B.; Palmer, D.; Teegarden, B.; Cline, T.; Briggs, M.; Paciesas, W.; Pendleton, G.; et al. BATSE Observations of Gamma-Ray Burst Spectra. I. Spectral Diversity. Astrophys. J. 1993, 413, 281. [Google Scholar] [CrossRef]
  2. Gruber, D.; Goldstein, A.; von Ahlefeld, V.W.; Bhat, P.N.; Bissaldi, E.; Briggs, M.S.; Byrne, D.; Clevel, W.H.; Connaughton, V.; Diehl, R.; et al. The Fermi GBM Gamma-Ray Burst Spectral Catalog: Four Years of Data. Astrophys. J. Suppl. 2014, 211, 12. [Google Scholar] [CrossRef] [Green Version]
  3. Preece, R.D.; Briggs, M.S.; Mallozzi, R.S.; Pendleton, G.N.; Paciesas, W.S.; Band, D.L. The BATSE Gamma-Ray Burst Spectral Catalog. I. High Time Resolution Spectroscopy of Bright Bursts Using High Energy Resolution Data. Astrophys. J. 2000, 126, 19. [Google Scholar] [CrossRef] [Green Version]
  4. Zhang, B. Gamma-ray Burst Prompt Emission. IJMPD 2014, 23, 1430002. [Google Scholar] [CrossRef] [Green Version]
  5. Kumar, P.; Zhang, B. The Physics of Gamma-Ray Bursts and Relativistic Jets. Phys. Rep. 2014, 561, 1. [Google Scholar] [CrossRef] [Green Version]
  6. Beloborodov, A.M. Collisional Mechanism for Gamma-Ray Burst Emission. Mon. Not. R. Astron. Soc. 2010, 407, 1033. [Google Scholar] [CrossRef] [Green Version]
  7. Deng, W.; Zhang, B. Cosmological Implications of Fast Radio Burst/Gamma-Ray Burst Associations. Astrophys. J. 2014, 783, 35. [Google Scholar] [CrossRef] [Green Version]
  8. Giannios, D. Prompt GRB emission from gradual energy dissipation. Astron. Astrophys. 2008, 480, 305. [Google Scholar] [CrossRef] [Green Version]
  9. Lazzati, D.; Begelman, M.C. Non-Thermal Emission from the Photospheres of Gamma-Ray Burst Outflows. I. High-Frequency Tails. Astrophys. J. 2010, 725, 1137. [Google Scholar] [CrossRef] [Green Version]
  10. Lazzati, D.; Morsony, B.J.; Margutti, R.; Begelman, M.C. Photospheric Emission As the Dominant Radiation Mechanism in Long-Duration Gamma-Ray Bursts. Astrophys. J. 2013, 765, 103. [Google Scholar] [CrossRef] [Green Version]
  11. Lundman, C.; Pe’er, A.; Ryde, F. A Theory of Photospheric Emission from Relativistic, Collimated Outflows. Mon. Not. R. Astron. Soc. 2013, 428, 2430. [Google Scholar] [CrossRef] [Green Version]
  12. Meszaros, P.; Ramirez-ruiz, E.; Rees, M.J.; Zhang, B. X-Ray-rich Gamma-Ray Bursts, Photospheres, and Variability. Astrophys. J. 2002, 578, 812. [Google Scholar] [CrossRef]
  13. Meszaros, P.; Rees, M.J. Steep Slopes and Preferred Breaks in Gamma-Ray Burst Spectra: The Role of Photospheres and Comptonization. Astrophys. J. 2000, 530, 292. [Google Scholar] [CrossRef] [Green Version]
  14. Pe’er, A. Temporal Evolution of Thermal Emission from Relativistically Expanding Plasma. Astrophys. J. 2008, 682, 463. [Google Scholar] [CrossRef]
  15. Pe’er, A.; Meszaros, P.; Rees, M.J. The Observable Effects of a Photospheric Component on GRB and XRF Prompt Emission Spectrum. Astrophys. J. 2006, 642, 995. [Google Scholar] [CrossRef]
  16. Pe’er, A.; Ryde, F. A Theory of Multicolor Blackbody Emission from Relativistically Expanding Plasmas. Astrophys. J. 2011, 732, 49. [Google Scholar] [CrossRef]
  17. Rees, M.J.; Meszaros, P. Dissipative Photosphere Models of Gamma-Ray Bursts and X-Ray Flashes. Astrophys. J. 2005, 628, 847. [Google Scholar] [CrossRef]
  18. Thompson, C. A Model of Gamma-Ray Bursts. Mon. Not. R. Astron. Soc. 1994, 270, 480. [Google Scholar] [CrossRef]
  19. Thompson, C.; Meszaros, P.; Rees, M.J. Thermalization in Relativistic Outflows and the Correlation between Spectral Hardness and Apparent Luminosity in Gamma-Ray Bursts. Astrophys. J. 2007, 666, 1012. [Google Scholar] [CrossRef] [Green Version]
  20. Deng, W.; Zhang, B. Low Energy Spectral Index and Ep Evolution of Quasi-thermal Photosphere Emission of Gamma-Ray Bursts. Astrophys. J. 2014, 785, 112. [Google Scholar] [CrossRef] [Green Version]
  21. Pe’er, A.; Zhang, B. Synchrotron Emission in Small-Scale Magnetic Fields as a Possible Explanation for Prompt Emission Spectra of Gamma-Ray Bursts. Astrophys. J. 2006, 653, 454. [Google Scholar] [CrossRef] [Green Version]
  22. Uhm, Z.L.; Zhang, B. Fast-Cooling Synchrotron Radiation in a Decaying Magnetic Field and Γ-Ray Burst Emission Mechanism. Nat. Phys. 2014, 10, 351. [Google Scholar] [CrossRef]
  23. Zhang, B.B.; Uhm, Z.L.; Connaughton, V.; Briggs, M.S.; Zhang, B. Synchrotron Origin of the Typical GRB Band Function—A Case Study of GRB 130606B. Astrophys. J. 2016, 816, 72. [Google Scholar] [CrossRef]
  24. Geng, J.J.; Huang, Y.F.; Wu, X.F.; Zhang, B.; Zong, H.S. Low-energy Spectra of Gamma-Ray Bursts from Cooling Electrons. Astrophys. J. 2018, 234, 3. [Google Scholar] [CrossRef]
  25. Zhao, X.; Li, Z.; Liu, X.; Zhang, B.B.; Bai, J.; Meszaros, P. Gamma-Ray Burst Spectrum with Decaying Magnetic Field. Astrophys. J. 2014, 780, 12. [Google Scholar] [CrossRef] [Green Version]
  26. Daigne, F.; Bosnjak, Z.; Dubus, G. Reconciling observed gamma-ray burst prompt spectra with synchrotron radiation? Astron. Astrophys. 2011, 526, 110. [Google Scholar] [CrossRef] [Green Version]
  27. Wang, X.Y.; Li, Z.; Dai, Z.G.; Meszaros, P. GRB 080916C: On the Radiation Origin of the Prompt Emission from keV/MeV TO GeV. Astrophys. J. 2009, 698, 98. [Google Scholar] [CrossRef] [Green Version]
  28. Lloyd, N.M.; Petrosian, V. Synchrotron Radiation as the Source of Gamma-Ray Burst Spectra. Astrophys. J. 2000, 543, 722. [Google Scholar] [CrossRef]
  29. Shen, R.F.; Zhang, B. Prompt Optical Emission and Synchrotron Self-Absorption Constraints on Emission Site of GRBs. Mon. Not. R. Astron. Soc. 2009, 398, 1936. [Google Scholar] [CrossRef] [Green Version]
  30. Xu, S.Y.; Yang, Y.P.; Zhang, B. On the Synchrotron Spectrum of GRB Prompt Emission. Astrophys. J. 2018, 853, 43. [Google Scholar] [CrossRef] [Green Version]
  31. Bosnjak, Z.; Daigne, F. Spectral Evolution in Gamma-Ray Bursts: Predictions of the Internal Shock Model and Comparison to Observations. Astron. Astrophys. 2014, 568, A45. [Google Scholar] [CrossRef] [Green Version]
  32. Burgess, J.M. Time-resolved Analysis of Fermi Gamma-Ray Bursts with Fast- and Slow-cooled Synchrotron Photon Models. Astrophys. J. 2014, 784, 17. [Google Scholar] [CrossRef] [Green Version]
  33. Zhang, B.; Yan, H. The Internal-collision-induced Magnetic Reconnection and Turbulence (ICMART) Model of Gamma-ray Bursts. Astrophys. J. 2011, 726, 90. [Google Scholar] [CrossRef] [Green Version]
  34. Panaitescu, A.; Meszaros, P. Gamma-Ray Bursts from Upscattered Self-absorbed Synchrotron Emission. Astrophys. J. 2000, 544, 17. [Google Scholar] [CrossRef] [Green Version]
  35. Sari, R.; Piran, T.; Narayan, R. Spectra and Light Curves of Gamma-Ray Burst Afterglows. Astrophys. J. 1998, 497, L17. [Google Scholar] [CrossRef]
  36. Wang, K.; Dai, Z.G. GeV Emission during X-Ray Flares from Late Internal Shocks: Application to GRB 100728A. Astrophys. J. 2013, 772, 152. [Google Scholar] [CrossRef] [Green Version]
  37. Yu, Y.W.; Dai, Z.G. X-Ray and High Energy Flares from Late Internal Shocks of Gamma-Ray Bursts. Astrophys. J. 2009, 692, 133. [Google Scholar] [CrossRef]
  38. Dai, Z.G.; Lu, T. Hydrodynamics of Relativistic Blast Waves in a Density-Jump Medium and Their Emission Signature. Astrophys. J. 2002, 565, L87. [Google Scholar] [CrossRef] [Green Version]
  39. Kumar, P.; Piran, T. Some Observational Consequences of Gamma-Ray Burst Shock Models. Astrophys. J. 2000, 532, 286. [Google Scholar] [CrossRef] [Green Version]
  40. Sari, R.; Piran, T. Hydrodynamic Timescales and Temporal Structure of Gamma-Ray Bursts. Astrophys. J. 1995, 455, L143. [Google Scholar] [CrossRef]
  41. Blandford, R.; McKee, C. Fluid Dynamics of Relativistic Blast Waves. Phys. Fluids 1976, 19, 1130. [Google Scholar] [CrossRef]
  42. Sari, R.; Esin, A.A. On The Synchrotron Self-Compton Emission from Relativistic Shocks and Its Implications for Gamma-Ray Burst Afterglows. Astrophys. J. 2001, 548, 787. [Google Scholar] [CrossRef] [Green Version]
  43. Abdo, A.A.; Ackermann, M.; Arimoto, M.; Asano, K.; Atwood, W.B.; Axelsson, M.; Baldini, L.; Ballet, J.; Band, D.L.; Barbiellini, G.; et al. Fermi Observations of High-Energy Gamma-Ray Emission from GRB 080916C. Science 2009, 323, 1688. [Google Scholar] [PubMed]
  44. Abdo, A.A.; Ackermann, M.; Asano, K.; Atwood, W.B.; Axelsson, M.; Baldini, L.; Ballet, J.; Band, D.L.; Barbiellini, G.; Bastieri, D.; et al. Fermi observations of high-energy gamma-ray emission from GRB 080825C. Astrophys. J. 2009, 707, 580. [Google Scholar] [CrossRef] [Green Version]
  45. Zhang, B. Open Questions in GRB Physics. Comptes Rendus Phys. 2011, 12, 206. [Google Scholar] [CrossRef] [Green Version]
  46. Derishev, E.V.; Kocharovsky, V.V.; Kocharovsky, V.V. Physical Parameters and Emission Mechanism in Gamma-Ray Bursts. Astron. Astrophys. 2011, 372, 1071. [Google Scholar] [CrossRef] [Green Version]
  47. Piran, T.; Sari, R.; Zou, Y.C. Observational Limits on Inverse Compton Processes in Gamma-Ray Bursts. Mon. Not. R. Astron. Soc. 2009, 393, 1107. [Google Scholar] [CrossRef] [Green Version]
  48. Racusin, J.L.; Oates, S.R.; Schady, P.; Burrows, D.N.; De Pasquale, M.; Donato, D.; Gehrels, N.; Koch, S.; McEnery, J.; Piran, T.; et al. Fermi and Swift Gamma-ray Burst Afterglow Population Studies. Astrophys. J. 2011, 738, 138. [Google Scholar] [CrossRef] [Green Version]
  49. Yonetoku, D.; Murakami, T.; Nakamura, T.; Yamazaki, R.; Inoue, A.K.; Ioka, K. Gamma-Ray Burst Formation Rate Inferred from the Spectral Peak Energy-Peak Luminosity Relation. Astrophys. J. 2004, 609, 935. [Google Scholar] [CrossRef] [Green Version]
  50. Amati, L.; Frontera, F.; Tavani, M.; Antonelli, A.; Costa, E.; Feroci, M.; Guidorzi, C.; Heise, J.; Masetti, N.; Montanari, E.; et al. Intrinsic Spectra and Energetics of BeppoSAX Gamma–Ray Bursts with Known Redshifts. Astron. Astrophys. 2002, 390, 81. [Google Scholar] [CrossRef]
  51. Yu, H.F.; Van Eerten, H.J.; Greiner, J.; Sari, R.E.; Bhat, P.N.; Von Kienlin, A.; Paciesas, W.S.; Preece, R.D. The Sharpness of Gamma-Ray Burst Prompt Emission Spectra. Astron. Astrophys. 2015, 583, 129. [Google Scholar] [CrossRef] [Green Version]
  52. Guiriec, S.; Connaughton, V.; Briggs, M.S.; Burgess, M.; Ryde, F.; Daigne, F.; Meszaros, P.; Goldstein, A.; McEnery, J.; Omodei, N.; et al. Detection of a Thermal Spectral Component in the Prompt Emission of GRB 100724B. Astrophys. J. Lett. 2011, 727, L33. [Google Scholar] [CrossRef]
  53. Oganesyan, G.; Nava, L.; Ghirlanda, G.; Celotti, A. Detection of Low-energy Breaks in Gamma-Ray Burst Prompt Emission Spectra. Astrophys. J. 2017, 846, 137. [Google Scholar] [CrossRef] [Green Version]
  54. Daigne, F. Mochkovitch, R.The Expected Thermal Precursors of Gamma-Ray Bursts in the Internal Shock Model. Mon. Not. R. Astron. Soc. 2002, 336, 1271. [Google Scholar] [CrossRef] [Green Version]
  55. Ryde, F. Is Thermal Emission in Gamma-Ray Bursts Ubiquitous? Astrophys. J. Lett. 2005, 625, L95. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The magnetic field as a function of time. The two blue vertical dotted line represents the redefined interval δ t = 0.1 s and the shock crossing time t c r s = 3 s , respectively. After the shock crossing time, the merged shell expands adiabatically and s = 3 is assumed. The dynamics parameters L k = 10 51 erg s 1 , γ 1 = 300 , γ 4 = 30,000 , p = 2.5 , ϵ e = 0.3 , ϵ B = 0.3 , and z = 1 are taken from numerical calculations.
Figure 1. The magnetic field as a function of time. The two blue vertical dotted line represents the redefined interval δ t = 0.1 s and the shock crossing time t c r s = 3 s , respectively. After the shock crossing time, the merged shell expands adiabatically and s = 3 is assumed. The dynamics parameters L k = 10 51 erg s 1 , γ 1 = 300 , γ 4 = 30,000 , p = 2.5 , ϵ e = 0.3 , ϵ B = 0.3 , and z = 1 are taken from numerical calculations.
Galaxies 09 00068 g001
Figure 2. The electron distribution in energy space after cooling time t in the evolutional magnetic field in Figure 1. The same δ t , t c r s , and dynamics parameters as in Figure 1 are taken in numerical calculations.
Figure 2. The electron distribution in energy space after cooling time t in the evolutional magnetic field in Figure 1. The same δ t , t c r s , and dynamics parameters as in Figure 1 are taken in numerical calculations.
Galaxies 09 00068 g002
Figure 3. The top panel corresponds to time-resolved spectra in four different t as in Figure 2 and the bottom panel shows corresponding synchrotron spectral slopes. The same δ t , t c r s , and dynamics parameters as in Figure 1 are taken in numerical calculations.
Figure 3. The top panel corresponds to time-resolved spectra in four different t as in Figure 2 and the bottom panel shows corresponding synchrotron spectral slopes. The same δ t , t c r s , and dynamics parameters as in Figure 1 are taken in numerical calculations.
Galaxies 09 00068 g003
Figure 4. The top panel shows the electron distributions in evolutional magnetic fields and different electron injection rising indices. We adopt the electron injection rising index, q = 0 (solid line), q = 1 (dotted line) and q = 2 (dashed line). The bottom panel shows the corresponding synchrotron spectral slopes for these electron distributions. The same parameters as in Figure 1 are taken for numerical calculations.
Figure 4. The top panel shows the electron distributions in evolutional magnetic fields and different electron injection rising indices. We adopt the electron injection rising index, q = 0 (solid line), q = 1 (dotted line) and q = 2 (dashed line). The bottom panel shows the corresponding synchrotron spectral slopes for these electron distributions. The same parameters as in Figure 1 are taken for numerical calculations.
Galaxies 09 00068 g004
Figure 5. The time-averaged spectrum to fit the interval “b” of GRB 080916c. The black points are selected from the Band function with the low-energy photon index α = 1.02 ± 0.02 , the high-energy photon index β = 2.21 ± 0.03 and the peak energy E p = 1170 ± 140 keV provided in Ref. [43], which are precise enough to present the tendency of the actual GRB emission. This spectral duration is from 3.58 s to 7.68 s and we fit it by adopting the time-averaged spectrum from t = 0 s to t = 3 s ( t m a x t c r s ). The fitting parameters are listed in Table 1.
Figure 5. The time-averaged spectrum to fit the interval “b” of GRB 080916c. The black points are selected from the Band function with the low-energy photon index α = 1.02 ± 0.02 , the high-energy photon index β = 2.21 ± 0.03 and the peak energy E p = 1170 ± 140 keV provided in Ref. [43], which are precise enough to present the tendency of the actual GRB emission. This spectral duration is from 3.58 s to 7.68 s and we fit it by adopting the time-averaged spectrum from t = 0 s to t = 3 s ( t m a x t c r s ). The fitting parameters are listed in Table 1.
Galaxies 09 00068 g005
Figure 6. The time-averaged spectrum to fit the interval “a” of GRB 080825c. The black points are selected from the Band function with the low-energy photon index α = 0.76 ± 0.05 , the high-energy photon index β = 2.54 0.17 + 0.11 and the peak energy E p = 291 22 + 25 keV provided in Ref. [44], which are precise enough to present the tendency of the actual GRB emission. This spectral duration is from 0.0 s to 2.7 s and we fit it by adopting the time-averaged spectrum from t = 0 s to t = 3 s ( t m a x t c r s ). The fitting parameters are listed in Table 1.
Figure 6. The time-averaged spectrum to fit the interval “a” of GRB 080825c. The black points are selected from the Band function with the low-energy photon index α = 0.76 ± 0.05 , the high-energy photon index β = 2.54 0.17 + 0.11 and the peak energy E p = 291 22 + 25 keV provided in Ref. [44], which are precise enough to present the tendency of the actual GRB emission. This spectral duration is from 0.0 s to 2.7 s and we fit it by adopting the time-averaged spectrum from t = 0 s to t = 3 s ( t m a x t c r s ). The fitting parameters are listed in Table 1.
Galaxies 09 00068 g006
Table 1. The parameters adopted to fit the spectra of GRB 080916c and GRB 080825c.
Table 1. The parameters adopted to fit the spectra of GRB 080916c and GRB 080825c.
ParametersSymbolGRB 080916cGRB 080825c
Redshiftz4.351
Index of electron injection rateq02
redefined time interval (s) δ t 0.10.1
Shock cross time (s) t c r s 33
Kenetic luminosity (erg/s) L k 3.3 × 10 53 1.2 × 10 51
Bulk Lorentz factor of region 1 γ 1 146255
Bulk Lorentz factor of region 2 γ 4 3 × 10 4 3 × 10 4
Electron injection indexp2.53.2
Electron equipartition factor ϵ e 0.30.3
Magnetic equipartition factor ϵ B 0.30.3
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, K.; Dai, Z.-G. The Low-Energy Spectral Index of Gamma-ray Burst Prompt Emission from Internal Shocks. Galaxies 2021, 9, 68. https://doi.org/10.3390/galaxies9030068

AMA Style

Wang K, Dai Z-G. The Low-Energy Spectral Index of Gamma-ray Burst Prompt Emission from Internal Shocks. Galaxies. 2021; 9(3):68. https://doi.org/10.3390/galaxies9030068

Chicago/Turabian Style

Wang, Kai, and Zi-Gao Dai. 2021. "The Low-Energy Spectral Index of Gamma-ray Burst Prompt Emission from Internal Shocks" Galaxies 9, no. 3: 68. https://doi.org/10.3390/galaxies9030068

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop