A classical model based on a power law assumption for the radius-time relationship in the expansion of a supernova (SN) allows to derive an analytical expression for the flow of mechanical kinetic energy and the time duration of gamma-ray burst (GRB). A random process based on the ratio of two truncated lognormal distributions for luminosity and luminosity distance allows to derive the statistical distribution for time duration of GRBs. The high velocities involved in the first phase of expansion of a SN require a relativistic treatment. The circumstellar medium is assumed to follow a density profile of Plummer type with eta = 6. A series solution for the relativistic flow of kinetic energy allows to derive in a numerical way the duration time for GRBs. Here we analyze two cosmologies: the standard cosmology and the plasma cosmology.
The theoretical efforts for gamma-ray bursts (GRBs) analyze: 1) the different predictions between cannonball and fireball model, see [
This paper reviews in Section 2 the current status of the observations for the time duration of GRBs. Section 3 develops a simple classic model for theoretical time duration for GRBs. Section 4 contains a simulation for time duration of GRBs based on the random generation of luminosity and distance luminosity. Section 5 derives a relativistic equation of motion for SN, deduces the relativistic flow of energy and finally evaluates the time duration for GRBS in a relativistic environment.
This section reviews the observations of the Fermi Gamma-ray Burst Monitor (GBM) for the time duration and fits the data with the truncated lognormal PDF.
Most of the typical GRBs have two stages: the initial prompt emission phase observed in the gamma-ray band and the following afterglow phase with multibands’ emission. The observational time duration is measured during the prompt emission phase with the typical energy band between tens of keV to hundreds of keV. Observationally, the light curve of the afterglow phase can be fitted by a broken power-law, but the light curve during the prompt emission phase is highly variable and completely different from the afterglow phase. The time duration of GRBs has been monitored by GBM in the 50 - 300 keV energy range, see the third catalog [
A simple test on the isotropy of the arrival direction of GRBs can be performed dividing the the surface area of a sphere of unit radius in N × N boxes of equal area. We then choose N in order to have » 7 theoretical events, n t h , assuming a uniform distribution on the total surface area, which means N = 14 . We now evaluate the observed averaged number of GRBs in each box, n, which turns out to be n = 7.16 ± 2.7 . The goodness of the isotropy is evaluated through the percentage error, δ , which is
δ = | n t h − n | n t h × 100 = 2.4 % . (1)
The low value of the percentage error for the isotropy allows to state that:
• the GRBs have an extra-galactic origin;
• the spatial distribution of GRBs is isotropic.
The histogram of time distribution for GRBs for the case of T 50 in
Scheme 1 and Scheme 2 report the main statistical parameters of T 50 and T 90 respectively.
The high number of decades, » 4, covered by the distribution of T 50 and T 90 suggests to fit the data with the truncated lognormal (TL) PDF
T L ( x ; m , σ , x l , x u ) = 2 e − 1 2 1 σ 2 ( ln ( x m ) ) 2 π σ ( − erf ( 1 2 2 σ ln ( x l m ) ) + erf ( 1 2 2 σ ln ( x u m ) ) ) x , (2)
where x is the random variable, x l is the lower bound, x u is the upper bound, m is the scale parameter and σ is the shape parameter, see [
The DF is
D F ( x ; m , σ , x l , x u ) = − erf ( 1 2 2 σ ln ( x m ) ) + erf ( 1 2 2 σ ln ( x l m ) ) erf ( 1 2 2 σ ln ( x l m ) ) − erf ( 1 2 2 σ ln ( x u m ) ) . (3)
Scheme 1. The sample parameters of T 50 for GRBs in seconds.
Scheme 2. The sample parameters of T 90 for GRBs in seconds.
Scheme 3. The TL parameters of T 50 for GRBs.
Scheme 4. The TL PDF and lognormal PDF parameters of T 90 for GRB.
A random generation of the variate X can be found by solving the following non linear equation
R = D F ( X ; m , σ , x l , x u ) , (4)
where R is the unit rectangular variate.
We assume that the observed radius-time relationship, R ( t ) , for SN has a power law dependence of the type
R ( t ) = C × t α , (5)
where C is a constant and α a parameter which can be fixed by the observation of the temporal evolution of the radius of a SN, in the case of SN 1993J α = 0.82 , see [
R 0 ( t ) = C × t 0 α . (6)
The velocity, V ( t ) , is
V ( t ) = C × α t ( α − 1 ) , (7)
and the velocity at time t 0 is
V 0 ( t ) = C × α t 0 ( α − 1 ) . (8)
In classical physics the density of kinetic energy, K, is
K = 1 2 ρ V 2 , (9)
where ρ is the density and V the velocity. In presence of an area A and when the velocity is perpendicular to that area, the mechanical flux of kinetic energy L m is
L m = 1 2 ρ A V 3 , (10)
which in SI is measured in W and in CGS in erg s−1 see Formula (A28) in [
L m = 1 2 ρ 4 π R 2 V 3 . (11)
The density in the advancing layer as function of the radius, ρ ( R ) , is assumed to scale as
ρ ( R ) = ρ 0 ( R R 0 ) − d , (12)
where ρ 0 is the density at radius R 0 . According to previous formula, the scaling law for the mechanical luminosity as function of the time is
L m = L m ,0 ( t t 0 ) 5 α − d α − 3 , (13)
where L m ,0 is the mechanical luminosity at t = t 0 . The observed luminosity at a given frequency ν , L , is assumed to be proportional to the mechanical luminosity
L = ϵ L m , (14)
where ϵ is a constant of conversion. The luminosity at a given range of energy is expressed in erg/s in CGS and W in SI. As a useful example the astrophysical version of the luminosity at t = t 0 is
L = 2.68521 n 0 ( R 0, p c ) 2 ( V 0 c ) 3 10 32 W , (15)
where R 0, p c is the radius at t = t 0 in pc, V 0 is the velocity at t = t 0 in km/s, c is the light velocity in km/s and n 0 is the number density expressed in particles cm−3 (density ρ 0 = n 0 m , where m = m H ).
The observed flux, F , as a function of the luminosity distance, D L , is
F = L 4 π D 2 , (16)
or
F ( L , D L ) = L 0 ( t t 0 ) 5 α − d α − 3 4 π D L 2 , (17)
where L 0 = ϵ L m , 0 . The observed flux at a given range of energy is expressed in erg/scm2 in CGS and W/m2 in SI. An on line collection of light curves for GRBs is made by the Swift GRB Mission and can be found at http://swift.gsfc.nasa.gov/index.html, see also [
observations is (0.3 - 10) kev. As a practical example we processed GRB161214B, see [
Scheme 5. Some parameters of GRB161214B.
A first numerical analysis of the observed flux, F 0 , versus time relationship can be done by assuming a power law dependence
F 0 ( t ) = C k t k erg ⋅ s − 1 , (18)
where the two parameters C and k are found from a numerical analysis of the data, see Scheme 5. We now match the observed value of k with the theoretical exponent given by Equation (17)
k = 5 α − d α − 3. (19)
The above equation allows to deduce d once α and k are given: as an example, when α = 0.82 and k = − 2.54 , d = 4.44 . The decreasing flux will be visible until a minimum value of threshold in the flux F T is reached. As an example the instrument Windowed Timing (WT) of X-ray telescope (XRT) on the Swift satellite has a threshold value of F T ≈ 3.1 fW / m 2 . Therefore the GRB will be visible up to a maximum value in time of t max
t max = t 0 ( F T 4 π D L 2 L 0 ) 1 5 α − d α − 3 , (20)
or
t max = t 0 ( F T F 0 ) 1 5 α − d α − 3 , (21)
where F 0 is the observed flux at t = t 0 .
A simulation of the duration time for GRBs according to Equation (20) requires the ratio of two two lognormal PDFs (one for L 0 and the other one for D L ) and a cosmological environment.
We now evaluate the mean and the variance of the ratio of two lognormal
distributions, X Y . From the fact that ln ( X Y ) = ln ( X ) − ln ( Y ) and X and Y
are lognormally distributed, it turns out that ln ( X ) and ln ( Y ) are distributed as a normal PDF. We assume that ln ( X ) and ln ( Y ) have means μ X and μ Y , variances σ X 2 and σ Y 2 , and covariance σ X Y 2 (equal to zero if X and Y are independent) and are jointly normally distributed. The difference Z is then distributed as a normal PDF with mean μ Z = μ X − μ Y and variance σ Z 2 = σ X 2 + σ Y 2 − 2 σ X Y 2 .
Note that X Y = exp ( Z ) , which means that X Y is distributed as a lognormal
PDF with parameters μ Z and σ Z 2 .
We deal with two cosmologies: the LCDM cosmology and the plasma cosmology. The the LCDM cosmology is characterized by three parameters which are the Hubble constant, H 0 , expressed in kms−1Mpc−1 and the two numbers Ω M and Ω Λ ,, see Scheme 6.
A simple expression for the luminosity distance, D l , in LCDM is obtained with the minimax approximation ( p = 3 , q = 2 )
D L , 3 , 2 = 0.3597252600 + 5.612031882 z + 5.627811123 z 2 + 0.05479466285 z 3 0.0105878216 + 0.1375418627 z + 0.1159043801 z 2 . (22)
More details on the analytical derivation of the luminosity distance in terms of Padé approximant can be found in [
The plasma cosmology is characterized by one parameter which is H 0 and by a simple formula for the luminosity distance which is the same of the distance, d p ,
D L ( z ; H 0 ) ≡ d p ( z ; H 0 ) = 0.359725 + 5.61203 z + 5.62781 z 2 + 0.0547946 z 3 0.010587 + 0.137541 z + 0.115904 z 2 , (23)
see [
A careful examination of Equation (20) which gives the time duration for a GRB isolates five fixed basic parameters which are α , d , R 0 , V 0 , t 0 and two random parameters are , L , D L . In order to allow the simulation the five fixed parameters are reported in Scheme 7.
The variable L is generated in a random way which follows a TL PDF with the following meanings: L l lower luminosity, L u upper luminosity and L * the scale, see Scheme 8.
Scheme 6. Numerical values for parameters of the two cosmologies.
Scheme 7. The fixed parameters of the simple model.
Scheme 8. The parameters of the simulation for L in LCDM and plasma cosmology.
The theoretical luminosity at t = t 0 , L 0 , is obtained inserting in Equation (15) the number density n 0 for which L = L 0 .
The variable D L , in the case of LCDM or plasma cosmology, is generated in a random way which follows a TL PDF with the following meanings: D L , l lower luminosity distance, L L , u upper luminosity distance and m the scale , see Scheme 9.
Once a sequence of theoretical times of duration is obtained according to Equation (20), see results in Scheme 10, we compare the observed and simulated data, see
Scheme 9. The parameters of the simulation for D L in LCDM and plasma cosmology.
Scheme 10. Theoretical time duration for GRBs in seconds.
The theoretical reason which allows to fit the ratio of two lognormal PDFs with another lognormal is outlined in Section 4.1.
The density, ρ , of the ISM at a distance r from the SN is here modeled by a Plummer-like profile, see [
ρ ( r ; R f l a t ) = ρ c ( R f l a t ( R f l a t 2 + r 2 ) 1 / 2 ) η . (24)
where r is the distance from the center, ρ is the density, ρ c is the density at the center, R f l a t is the distance before which the density is nearly constant, and η is the power law exponent at large values of r. The following transformation, R f l a t = 3 b , gives
ρ ( r ; b ) = ρ c ( 1 1 + 1 3 r 2 b 2 ) η / 2 . (25)
The total mass M ( r ; b ) comprised between 0 and r, when η = 6 , is
M ( r ; b ) = ∫ 0 r 4 π r 2 ρ ( r ; b ) d r = 3 2 ρ _ c , π b 3 ( 3 b 2 + r 2 ) 2 ( 9 arctan ( 1 3 r 3 b ) 3 b 4 + 6 a r c t a n ( 1 3 r 3 b ) 3 b 2 r 2 + a r c t a n ( 1 3 r 3 b ) 3 r 4 − 9 b 3 r + 3 b r 3 ) . (26)
The relativistic conservation of momentum, see [
M ( R 0 ; b ) γ 0 β 0 = M ( r ; b ) γ β , (27)
with
γ 0 = 1 1 − β 0 2 ; γ = 1 1 − β 2 , (28)
and
β 0 = V 0 c ; β = v c , (29)
where R 0 is the initial radius of the advancing sphere, V 0 is the initial velocity at R 0 and c is the light velocity.
The relativistic conservation of momentum for a Plummer profile with η = 6 is
A N A D = B N B D , (30)
A N = 3 π b 3 ( 9 arctan ( 1 3 r ( t ) 3 b ) 3 b 4 + 6 arctan ( 1 3 r ( t ) 3 b ) 3 b 2 ( r ( t ) ) 2 + arctan ( 1 3 r ( t ) 3 b ) 3 ( r ( t ) ) 4 − 9 b 3 r ( t ) + 3 b ( r ( t ) ) 3 ) d d t r ( t ) (31)
A D = 2 ( 3 b 2 + ( r ( t ) ) 2 ) 2 c 1 − ( d d t r ( t ) ) 2 c 2 , (32)
B N = 3 π b 3 ( 9 a r c t a n ( 1 3 R 0 3 b ) 3 b 4 + 6 a r c t a n ( 1 3 R 0 3 b ) 3 b 2 R 0 2 + a r c t a n ( 1 3 R 0 3 b ) 3 R 0 4 − 9 b 3 R 0 + 3 b R 0 3 ) β 0 , (33)
B D = 2 ( 3 b 2 + R 0 2 ) 2 1 − β 0 2 . (34)
The relativistic transfer of energy through a surface, A , is
L m , r = A γ 2 ( ρ c 2 + p ) v , (35)
where p is the pressure here assumed to be p = 0 ; see Equation (A31) in [
The astrophysical version of the relativistic transfer of energy (the luminosity) at t = t 0 is
L = 5.37 n 0 R 0, p c 2 β 0 1 − β 0 2 10 32 W , (36)
where R 0, p c is the radius at t = t 0 in pc.
In the case of a spherical cold expansion
L m , r = 4 π r ( t ) 2 1 1 − β ( t ) 2 ρ ( t ) c 3 β ( t ) . (37)
We now assume the following power law behavior for the density in the advancing thin layer
ρ ( t ) = ρ 0 ( t 0 t ) d , (38)
and we obtain
L m , r = 4 π r ( t ) 2 1 1 − β ( t ) 2 ρ 0 ( t 0 t ) d c 3 β ( t ) . (39)
We can now derive L m , r in two ways: 1) from a numerical evaluation of r(t) and v(t); 2) from a Taylor series of L m , r ( t ) of the type
L m , r ( t ) = ∑ n = 0 3 a n ( t − t 0 ) n . (40)
The first two coefficients are
a 0 = − ρ 0 4 π R 0 2 ( t 0 t ) d c 3 β 0
a 1 = ρ 0 A 1 N A 1 D ,
where
A 1 N = − 8 π R 0 ( t 0 t ) d c 4 β 0 2 ( 27 3 b 6 arctan ( 1 3 r 0 3 b ) + 27 3 b 4 R 0 2 arctan ( 1 3 R 0 3 b ) + 9 3 b 2 R 0 4 arctan ( 1 3 R 0 3 b ) + 3 R 0 6 arctan ( 1 3 R 0 3 b ) − 36 R 0 3 β 0 2 b 3 − 27 b 5 R 0 − 36 R 0 3 b 3 + 3 b R 0 5 ) , (41)
and
A 1 D = ( 27 3 b 6 arctan ( 1 3 R 0 3 b ) + 27 3 b 4 R 0 2 arctan ( 1 3 R 0 3 b ) + 9 3 b 2 R 0 4 arctan ( 1 3 R 0 3 b ) + 3 R 0 6 arctan ( 1 3 R 0 3 b ) − 27 b 5 R 0 + 3 b R 0 5 ) × ( 1 − β 0 2 ) . (42)
The relativistic theory is now applied to GRB050814 at 0.3 - 10 kev in the time interval 10−5 - 3 days, see [
The theoretical flux without absorption is given by Equation (16) and
The presence of the absorption can be modeled as
F = L 4 π D 2 ( 1 − e − τ ν ( t ) ) , (43)
where τ n u ( t ) is the optical thickness here assumed to be dependent from the time. As a model for τ ν as function of time we select a logarithmic polynomial approximation, of degree 5 , see [
In the relativistic case the theoretical luminosity is provided by a series solution of order 7, see Equation (36), The fixed parameters adopted for the simulation of duration time are reported in Scheme 12.
The random luminosity L is generated according to a TL PDF, see Scheme 13.
The value of n 0 which allows to generate L 0 in a random way is deduced by L = L 0 with L as evaluated in Equation (32). The distance luminosity, D L , in the case of LCDM or plasma cosmology, i randomly generated according to a TL PDF, see Scheme 14.
The sequence of theoretical times of duration is obtained in a numerical way, see results in Scheme 15. The observed and simulated data are displayed in
Scheme 11. Numerical values of the parameters used in the Plummer relativistic solution.
Scheme 12. Numerical values of the parameters used in the Plummer relativistic simulation for duration time of GRBs.
Scheme 13. The parameters of the relativistic simulation for L in LCDM and plasma cosmology.
Scheme 14. The parameters of the relativistic simulation for DL in LCDM and plasma cosmology.
Scheme 15. Theoretical time duration for GRBs in seconds in a relativistic framework.
Extra-galactic versus Galactic origin
The direction of arrival of GRBs shows an isotropic universe with a percentage error of 2.4%. This means that the GRBs are originated in external galaxies in a random way as suggested by [
Truncated lognormal distribution
The lognormal PDF is usually adopted to model the duration time of GRBs, see [
A simple model
The theoretical time of duration can be derived from the light curve for GRBs. A power law approximation for the time of expansion of a shell, see Equation (5), coupled with a power law behavior for the density as function of the radius, see Equation (12), allows to derive a formula for the theoretical luminosity, see Equation (14). A theoretical time of duration is obtained, see Equation (20) which contains an evaluation of the luminosity and the luminosity distance. The cosmological evaluation of the luminosity distance is given both in LCDM cosmology and plasma cosmology. The results are given in Scheme 10,
Relativistic model
The temporal evolution of a SN in a medium of the Plummer type, η = 6 , can be found by applying the conservation of relativistic momentum in the thin layer approximation. This relativistic invariant is evaluated as a differential
equation of the first order, see Equation (26). Two different relativistic solutions for the theoretical luminosity as a function of time are derived: 1) a numerical solution, see Equation (35); 2) a series solution, see Equation (36), which has a limited range of validity, 10 − 1 s < t < 10 6 s .
The coupling of the previous series with a logarithmic polynomial approximation allows to model fine details such as the oscillation in LC visible at »1000 s in GRB161214B , see
The relativistic time of duration is reported in in Scheme 15,
An enlightening example of such relativistic simulation is absence of bursts at t ≈ 15.8 s visible in
This work made use of data supplied by the UK Swift Science Data Centre at the University of Leicester.
Zaninetti, L. (2017) Classical and Relativistic Models for Time Duration of Gamma-Ray Bursts. Journal of High Energy Physics, Gravitation and Cosmology, 3, 808-827. https://doi.org/10.4236/jhepgc.2017.34059