An analytic derivation of the empirical correlations of gamma-ray bursts

Empirical correlations between various key parameters have been extensively explored ever since the discovery of gamma-ray bursts (GRBs) and have been widely used as standard candles to probe the Universe. The Amati relation and the Yonetoku relation are two good examples that enjoyed special attention. The former reflects the connection between the peak photon energy (Ep) and the isotropic $\gamma$-ray energy release (Eiso), while the latter links Ep with the isotropic peak luminosity (Lp), both in the form of a power-law function. Most GRBs are found to follow these correlations well, but a theoretical interpretation is still lacking. Some obvious outliers may be off-axis GRBs and may follow correlations that are different from those at the on-axis. Here we present a simple analytical derivation for the Amati relation and the Yonetoku relation in the framework of the standard fireball model, the correctness of which is then confirmed by numerical simulations. The off-axis Amati relation and Yonetoku relation are also derived. They differ markedly from the corresponding on-axis relation. Our results reveal the intrinsic physics behind the radiation processes of GRBs, and they highlight the importance of the viewing angle in the empirical correlations of GRBs.


Introduction
After about five decades of research, some insights into gammaray bursts (GRBs) have been obtained.Generally speaking, there are two different phases in GRBs: the prompt emission, and the afterglow.The afterglow can be relatively well interpreted by the external shock model (Mészáros & Rees 1997;Sari et al. 1998;Huang et al. 2006;Geng et al. 2016;Lazzati et al. 2018;Xu et al. 2022).However, the radiation process of the prompt emission is still under debate.Different models have been proposed to explain the complicated prompt emission, such as the internal shock model (Rees & Meszaros 1994;Kobayashi et al. 1997;Bošnjak et al. 2009), the dissipative photosphere model (Rees & Mészáros 2005;Pe'er & Ryde 2011), and the internalcollision-induced magnetic reconnection and turbulence (IC-MART) model (Zhang & Yan 2011;Zhang & Zhang 2014).The most commonly discussed model is the internal shock model, which is naturally expected for a highly variable central engine.In this model, the collision and merger of shells create relativistic shocks to accelerate particles.Then the accelerated particles will cause the observed GRB prompt emission.
The spectrum of GRB prompt emission has traditionally been described by the empirical Band function (Band et al. 1993).It has been suggested that synchrotron radiation may be responsible for the non-thermal component (Rees & Meszaros 1994;Sari et al. 1998).However, for a standard fast-cooling spectrum, the theoretical low-energy spectral index is too soft (Sari et al. 1998).Several studies have shown that this problem can be mitigated when a decaying magnetic field and the e-mail: hyf@nju.edu.cndetailed cooling process are considered (Uhm & Zhang 2014;Zhang et al. 2016;Geng et al. 2018).Some authors argued that the low-energy index of the empirical Band function may be misleading (Burgess et al. 2020).They proposed that the synchrotron emission mechanism can interpret most spectra and can satisfactorily fit the observational data (Zhang et al. 2020).
In addition to the spectrum, the diversity of GRB energy is another mystery.For a typical long GRB (duration longer than 2 s), the isotropic energy is around 10 52 − 10 54 erg.However, there also exist some low-luminosity GRBs (LLGRBs) with energies down to 10 48 erg (GRB 980425) for long GRBs and 10 46 erg (GRB 170817A) for short GRBs (duration shorter than 2 s).Some authors claimed that LLGRBs may form a distinct population of GRBs (Liang et al. 2007).However, a later study reported no clear separation between LLGRBs and standard highluminosity GRBs in a larger GRB sample (Sun et al. 2015).On the other hand, LLGRBs can naturally be explained by off-axis jets.Due to the relativistic beaming effect, a GRB will become dimmer when the viewing angle (θ obs ) is larger than the jet opening angle (θ jet ) (Granot et al. 2002;Huang et al. 2002;Yamazaki et al. 2003a).The interesting short GRB 170817A clearly proves that at least some LLGRBs are viewed off-axis (Granot et al. 2017).
The GRB empirical relations can help us understand the physical nature of GRBs.The most famous relation is the socalled Amati relation (Amati et al. 2002).This relation connects the isotropic energy (E iso ) and the peak photon energy (E p ).The index of the Amati relation (E p − E iso ) is about 0.5 (Amati 2006;Nava et al. 2012;Demianski et al. 2017;Minaev & Pozanenko 2020).However, it was found that this relation is not followed by some LLGRBs.These outliers of the Amati relation appear to follow a flatter track on the E p − E iso plane (Farinelli et al. 2021).Previous studies showed that the off-axis effect may play a role in this phenomenon (Ramirez-Ruiz et al. 2005;Dado & Dar 2012).
Analytical derivation of the on-axis and off-axis Amati relation indices has been attempted in several researches (Granot et al. 2002;Eichler & Levinson 2004;Ramirez-Ruiz et al. 2005;Dado & Dar 2012).Eichler & Levinson (2004) derived the index of the Amati relation as 0.5 by considering a uniform, axisymmetric jet with a hole cut out of it (i.e., a ring-shaped fireball).Dado & Dar (2012) argued that the Amati relation should have an index of 1/2 ± 1/6 in the framework of the so-called cannonball model.Here, we show that the standard fireball model will also lead to an index of 0.5 for the on-axis Amati relation.For the off-axis Amati relation, previous studies showed that the index should be 1/3 based on the effect of the viewing angle alone (Granot et al. 2002;Ramirez-Ruiz et al. 2005).However, the effect of the Lorentz factor was usually ignored or not fully considered for the off-axis cases.Here we argue that the Lorentz factor is equally responsible for determining the index of the offaxis Amati relation.We suggest that this index should be between 1/4 and 4/13 after fully considering both the effect of the viewing angle and the Lorentz factor.
Researchers have also tried to reproduce the Amati relation by means of numerical calculations (Yamazaki et al. 2004;Kocevski 2012;Mochkovitch & Nava 2015).Recently, Farinelli et al. (2021) used an empirical comoving-frame spectrum (a spectrum described by a smoothly broken power-law function) to simulate the prompt emission.Their radiation flux was calculated by averaging over the pulse duration.They found that the Amati relation should be E p ∝ E 0.5 iso for on-axis and E p ∝ E 0.25 iso for off-axis cases.Here, we consider the detailed synchrotron spectrum instead of an empirical one in our study, so that we can gain more insights into the detailed physics.Furthermore, the isotropic peak luminosity L p can be precisely calculated in our model.As a result, in addition to deriving the Amati relation, we can also examine the Yonetoku relation between E p and L p (Yonetoku et al. 2004), which was previously discussed only for the on-axis cases (Zhang et al. 2009;Ito et al. 2019).
Our paper is organized as follows.In Section 2 we briefly describe our model.Then, in Section 3, we present an analytic derivation for the quantities of E p , E iso , and L p for both on-axis and off-axis cases.The relations between these parameters are also derived.Numerical results are presented in Section 4. Next, in Section 5, Monte Carlo simulations are performed to confirm the theoretically derived correlations.The theoretical results are compared with the observational data in Section 6. Conclusions and discussion are presented in Section 7.

Jet model
We mainly focus on long GRBs.We used a simple top-hat jet model similar to that of Farinelli et al. (2021).For simplicity, the radiation was assumed to only last for a very short time interval (δt) in the local burst frame.In other words, photons are emitted from the shell at almost the same time.Hence the photon arrival time is largely decided by the curvature effect.A photon emitted at a radius r with a polar angle θ will reach the observer at where c is the speed of light, and z is the redshift of the source.θ 1 refers to the angle from which the photons first arrive.We have θ 1 = 0 for an on-axis jet (θ jet ≥ θ obs ), and θ 1 = θ obs − θ jet for an off-axis jet (θ jet < θ obs ).The emission is produced through internal shocks, which can naturally dissipate the kinetic energy of a baryonic fireball (Zhang 2018).As a result, we have r ∼ γ 2 d, where γ is the bulk Lorentz factor of the shell, and d is the initial separation between the clumps ejected by the central engine (Kobayashi et al. 1997).A schematic illustration of our model is presented in Figure 1.Usually, the light curve of a single GRB pulse is composed of a fast rising phase and an exponential decay phase (Li & Zhang 2021).Complicated processes may be involved in the prompt GRB phase, such as the hydrodynamics of the outflow and the cooling of electrons (Zhang et al. 2007).Analytical solutions will be too difficult to derive if these ingredients are included.Especially in the context of the internal shock model, our assumption of a small δt implies that the pulse is simply produced by the collision of two very thin shells.In realistic cases, a reverse shock may form and propagate backward during the collision.The pulse profile could then be affected by the hydrodynamics of the flow, not only by the curvature of the emitting shell.Depending on the initial distribution of the Lorentz factor, the effect can remain moderate, however, and the two-shell model we considered is approximately valid.
Electrons accelerated by internal shocks should be in the fast-cooling regime.They will follow a distribution of (Geng et al. 2018) with where γ e is the Lorentz factor of electrons, and N e is their total number.The spectral index p usually ranges from 2 -3 (Huang et al. 2000).γ m and γ max are the minimum and maximum electron Lorentz factor, respectively.γ m is less than the initial Lorentz factor of the injected electrons (γ inj ) in the fast-cooling regime (Burgess et al. 2020), while γ max can be calculated approximately as γ max ∼ 10 8 ( B 1G ) −0.5 (Dai & Lu 1999;Huang et al. 2000), where B is the magnetic field strength in the comoving frame.
In the comoving frame, the synchrotron radiation power of these electrons at frequency ν is (Rybicki & Lightman 1979) where e is the electron charge, and m e is the electron mass.ν c = 3γ 2 e eB /4πm e c and F(x) = x +∞ x K 5/3 (k)dk, where x = ν ν c and K 5/3 (k) is the Bessel function.
In the local burst frame, the angular distribution of the radiation power is (Rybicki & Lightman 1979) where times.This is called the curvature effect.Summing up the contribution from the whole jet, we can obtain the observed γ-ray light curve.We first consider a time interval of ∆t in the local burst frame.It corresponds to a thin ring of [θ, θ + ∆θ].When the electrons are uniformly distributed in the shell, the number of electrons in the ring is where f b = 1 − cos θ jet is the beaming factor, N tot is the total number of electrons in the shell, and φ(θ) is the annular angle of the ring.N tot can be derived from the isotropic-equivalent mass of the shell m sh as N tot = 2π f b m sh /m p , where m p is the mass of the proton.The annular angle φ(θ) takes the form of Since the radiation process only lasts for a short interval of δt, the energy emitted into one unit solid angle is This energy corresponds to a duration of ∆t = r sin θ∆θ/c in the local burster frame.The luminosity is then L ν = ∆E ν /∆t.For an observer at distance D L , the observed flux at t obs (θ) should be The total isotropic energy can be calculated as where is the end time of the pulse in the observer frame.The isotropic peak luminosity is where θ p is the angle corresponding to the peak time of the light curve.The peak photon energy E p can be derived from the timeintegrated spectrum.In the next section, we present an analytic derivation for the prompt γ-ray emission, paying special attention to the three parameters E iso , L p , and E p .

Analytic derivation of E p , E iso , and L p
We first consider the on-axis cases.From Equations 1 and 10, we can derive E iso as where θ 2 = θ jet + θ obs .Further combining Equation 4, we have We note that C 1 in Equation 3 is largely dependent on γ m and γ inj .In this study, a simple case of γ m = γ inj /10 is assumed, which naturally leads to C 1 ∝ γ inj .Then with Equations 2, 3, and 6, we can further obtain For an on-axis observer, most of the observed photons are emitted by electrons within a small angle around the line of sight.According to Equation 7, it is safe for us to simply take φ(θ) as 2π.Changing the integral order in Equation 14, we obtain e dγ e dθ. (15) In the last step above, the integral of F ν ν c has been simplified as a constant (k 0 ).Equation 15 can be further reduced as where k 1 and k 2 are integration constants.When θ << 1 and β ∼ 1, we have 1/(1 − β cos θ) ∼ 2γ 2 /(θ 2 γ 2 + 1).The integral in Equation 16 can therefore be approximated as Finally, combining Equations 16 and 17, it is easy to obtain E iso ∝ γ 2 inj m sh B 2 γ 2 .Now we derive the peak luminosity of L p .Since the main difference between Equations 10 and 11 is the integral of time, we can first consider the luminosity corresponding to a particular angle θ as The flux peaks at θ p ∼ 0, thus the peak luminosity is The peak photon energy in the observer frame can be derived by considering the standard synchrotron emission mechanism, Therefore, the peak photon energy in the comoving frame is Next, we consider the off-axis cases.Again, we first focus on L p .For an off-axis jet, the main difference is that φ(θ) can no longer be taken as 2π .Instead, it should be calculated as φ(θ) ≡ 2 arccos( cos θ jet −cos θ obs cos θ sin θ obs sin θ ). Equation 18 now becomes Here θ ∈ (θ obs −θ jet , θ obs +θ jet ).We define R(θ) = cos θ jet −cos θ obs cos θ sin θ obs sin θ .Then it is obvious that R(θ) ∈ (0, 1), which leads to the approximation arccos We note that R(θ) can be simplified as When θ > 1 γ , we have Combining Equations 20 and 21, we obtain In this case, L p will peak at θ p = 1 15 (16θ obs − 255θ 2 jet + θ 2 obs ) ∼ 16 15 (θ obs − θ jet ).We note that the condition of θ p > 1 γ leads to 16 15 (θ obs −θ jet ) > 1 γ .Under this condition, L p can be further written as The last factor is . Its effect is not important compared with (θ obs − θ jet ) −8.5 , and therefore, we ignore this factor for simplicity.Finally, we obtain L p ∝ γ 2 inj m sh d −1 B 2 γ −6 (θ obs − θ jet ) −8.5 .Now we continue to calculate E iso .Since most of the γ-ray energy is released in the decay stage of the light curve, we only need to integrate the energy over a θ ranging from θ p to θ obs +θ jet , that is, Here the value of I(θ obs , θ jet , θ) will not change significantly in the range from θ p to θ obs +θ jet .Hence it can be taken as a constant for simplicity.Since θ p ∼ 16 15 (θ obs − θ jet ), we obtain the result as E iso ∝ γ 2 inj m sh B 2 γ −4 (θ obs − θ jet ) −6.5 .Finally, we calculate E p .Similar to Equation 19, the peak energy in the rest frame is We see that E p ∝ γ 2 inj B γ −1 (θ obs − θ jet ) −2 .To summarize, in the on-axis cases, we have while in the off-axis cases (θ obs − θ jet > 1/γ), we obtain The equations derived above are very simple but intriguing.There are eight input parameters in our model: γ, B , γ inj , θ jet , θ obs , m sh c 2 , d, and the electron spectral index p.Seven of them are closely related to the prompt emission parameters of E p , E iso , and L p , while p is largely irrelevant.Of all the input parameters, γ plays the most important role in affecting the observed spectral peak energies and fluxes for on-axis cases (Ghirlanda et al. 2018).Usually, γ may vary in a wide range for different GRBs, from lower than 100 to more than 1000.In contrast, γ inj and m sh only vary in much narrower ranges.Taking γ inj and m sh approximately as constants, from Equation 26, we can easily derive the on-axis Amati relation as E p ∝ E 0.5 iso and the on-axis Yonetoku relation as E p ∝ L 0.5 p .In both relations, the power-law indices are 0.5, which is largely consistent with observations (Nava et al. 2012;Demianski et al. 2017).
In the off-axis cases, previous studies mainly focused on the effect of the variation of the viewing angle among different GRBs.Consequently, the power-law index of the off-axis Amati relation is derived as 1/3 based on the sharp-edge homogeneous jet geometry (Granot et al. 2002;Ramirez-Ruiz et al. 2005).Here, from our Equation 27, we see that E iso and L p sensitively depend on both the viewing angle and the Lorentz factor.This indicates that both γ and θ obs are important parameters that will affect the slope of the Amati relation and the Yonetoku relation.In Equation 27, if only the variation of the viewing angle is considered (i.e., the indices of (θ obs − θ jet ) in Equation 27 are taken into account, but the item of γ is omitted), an index of 4/13 will be derived for the off-axis Amati relation, which is very close to the previous result of 1/3.When the combined effect of varying γ and θ obs is included, the off-axis Amati relation is derived as E p ∝ E 1/4∼4/13 iso , and the corresponding off-axis Yonetoku relation is E p ∝ L 1/6∼4/17 p .The power-law indices become smaller than in the on-axis cases.

Numerical results
Some approximations have been made in deriving Equations 26 and 27.In this section, we carry out numerical simulations to confirm whether the above analytical derivations are correct.For convenience, a set of standard values were taken for the eight input parameters in our simulations, as shown in Table 1.When studying the effect of one particular parameter, we only changed this parameter, but fixed all other parameters at the standard values.Notes. (a) We set the standard value as θ obs = 0.0 for an on-axis jet and θ obs = 0.2 for an off-axis jet. (b) The value of m sh here is the isotropic equivalent mass.
To examine the correctness of our analytical derivations, we chose one input parameter as a variable, and standard values were taken for all other input parameters so that the dependence of E p , E iso , and L p on that particular parameter can be illustrated.
We first consider the effect of the bulk Lorentz factor γ for both on-axis and off-axis cases.Here we let γ vary between 50 -1000.Figure 2 shows the time-integrated spectrum obtained for different γ.The peak photon energy is correspondingly marked with a black star on the curve.
For on-axis cases (θ jet ≥ θ obs ), the numerical results are plotted in the left panel of Figure 2. When the value of γ increases, the peak photon energy and the flux increase simultaneously.However, a different trend is found for the off-axis cases (θ jet < θ obs ), as shown in the right panel of Figure 2.Both the flux and E p decrease as γ increases.The flux is significantly lower than that of the on-axis cases.The results are consistent with those of Farinelli et al. (2021).
Figure 2 clearly shows that both E p and the flux are positively correlated with γ in the on-axis cases, but they are negatively correlated with γ in the off-axis cases.Next, we wish to determine the precise relation between the prompt emission pa-rameters and the input parameters.Hence, we plot E p , E iso , and L p as functions of γ, B , γ inj , |θ obs − θ jet |, and p in Figure 3.
Generally, the simulation results are well consistent with our analytical results of Equations 26 and 27.Especially E p , E iso , and L p have a positive dependence on the Lorentz factor γ in the on-axis cases, but they have a negative dependence on γ in the off-axis cases.It is also interesting to note that the prompt emission parameters are nearly independent of θ obs − θ jet in the on-axis cases, since the line of sight is always within the jet cone.However, in the off-axis cases, the dependence of the prompt emission parameters on θ obs − θ jet can be described as a twophase behavior.When θ obs − θ jet 1/γ, the prompt emission parameters are independent of θ obs − θ jet .When θ obs − θ jet > 1/γ, however, they decrease sharply with the increase of θ obs − θ jet , as indicated in Equation 27.Finally, the last column of Figure 3 shows that the parameter p has little impact on the prompt emission parameters.

Monte Carlo simulations
In this section, we perform Monte Carlo simulations to generate a large number of mock GRBs to further test the existence of the Amati relation and the Yonetoku relation.For this purpose, we performed Monte Carlo simulations as follows.First, a group of eight input parameters is generated randomly, assuming that each parameter follows a particular distribution.This group of parameters defines a mock GRB.Second, our model is applied to numerically calculate the corresponding values of E p , E iso , and L p for this mock GRB.The above two steps are repeated until we obtain a large sample of mock GRBs.Finally, the mock bursts are plotted on the E p − E iso plane and the E p − L p plane, and a best-fit correlation is obtained for them.
For convenience, we designate the distribution function of a particular parameter a as D(a).Here a refers to the input parameters (i.e., γ, B , γ inj , θ jet , θ obs , p, m sh c 2 , and d).Three kinds of distributions are adopted in this study: (i) A normal Gaussian distribution, which is noted as N(µ, σ), where µ and σ refer to the mean value and standard deviation, respectively.(ii) A uniform distribution, which is noted as U(min, max), where min and max define the range of the parameter.Finally, (iii) a log-normal distribution, noted as D(log(a)), which means that log(a) follows a normal Gaussian distribution of N(µ, σ).The detailed distribution for each of the input parameters adopted in our Monte Carlo simulations is described below.
The ranges of some input parameters, such as γ, θ jet , and θ obs , can be inferred from observations.For the bulk Lorentz factor, a distribution of D(log(γ)) = N(2.2,0.8) was assumed, which is consistent with observational constraints (Liang et al. 2010;Ghirlanda et al. 2018).As for the angle parameters, we set θ jet as a uniform distribution of D(θ jet ) = U(0.04,0.2) for the onand off-axis cases (Frail et al. 2001;Wang et al. 2018).The distribution of θ obs was taken as D(θ obs ) = U(0.0,0.2) in the onaxis cases, together with the restriction of θ jet ≥ θ obs .On the other hand, in the off-axis cases, θ obs was taken as D(θ obs ) = U(0.1,0.5), together with θ jet < θ obs .
Other parameters cannot be directly measured from observations.Their distributions are then taken based on some theoretical assumptions.The isotropic ejecta mass is usually thought to range from 10 −3 M to 10 −1 M for a typical long GRB.Hence we assumed that m sh c 2 has a log-normal distribution.We set its mean value as 10 −2 M c 2 and its 2σ range as 10 −3 M c 2 to 10 −1 M c 2 , that is, D(log(m sh c 2 /erg)) = N(52.2,0.5).The distribution of d was taken as D(log(d/cm)) = N(10, 0.2), which, b a Some parameters (B , γ inj , p) concern the microphysics of relativistic shocks.The comoving magnetic field strength (B ) is hard to estimate from observations (Zhang & Yan 2011;Burgess et al. 2020).Here, a log-normal distribution with a mean value of 10 G and a standard deviation of 0.4 G was assumed for it, that is, D(log(B /G)) = N(1, 0.4).The mean value of γ inj was taken as 10 5 , with a distribution of D(log(γ inj )) = N(5, 0.2) (Burgess et al. 2020).For the parameter p, we simply assumed a normal distribution of D(p) = N(2.5, 0.4).
With the above distribution functions of the input parameters, two samples of mock GRBs were generated through Monte Carlo simulations.One sample included 200 on-axis events, and the other sample included 200 off-axis cases.E p , E iso , and L p were calculated for each mock event.
The distribution of the simulated bursts on the E p − E iso and E p − L p planes is shown in Figure 4.In the on-axis cases, the best-fit slope is 0.495 ± 0.015 for the E p − E iso correlation, and it is 0.511±0.016for the E p − L p correlation.These two indices are very close to our theoretical results of 0.5, proving credibility for our analytical derivations.The intrinsic scatters of both relations are about 0.21, indicating that the correlations are rather tight.
In the off-axis cases, the best-fit slope is 0.282 ± 0.006 for the E p − E iso correlation, which is consistent with our expected range of 1/4 ∼ 4/13 for the off-axis Amati relation.Similarly, on the E p − L p plane, the best-fit slope is 0.206 ± 0.006, which also agrees well with our theoretical range of 1/6 ∼ 4/17 for the off-axis Yonetoku relation.Another interesting point in Figure 4c and 4d is that most of the simulated off-axis GRBs have an E iso lower than 10 48 erg and an L p lower than 10 46 erg s −1 .Only about 5% of our simulated off-axis GRBs have greater E iso and L p .We thus argue that off-axis GRBs are usually LLGRBs.It also indicates that the majority of currently observed normal GRBs should be on-axis GRBs.In Figure 5, all the simulated on-axis and off-axis bursts are plotted in the same panel.Only a very small number of off-axis events are strong enough to be detected.Their distribution obviously deviates from that of the on-axis events on both the E p − E iso and E p − L p plane.
Using the mock bursts, we can also test the Ghirlanda relation (Ghirlanda et al. 2004), which links E p with the collimationcorrected energy of E γ = (1 − cos θ jet )E iso .The distribution of the simulated bursts on the E p − E γ plane is shown in Figure 6.The data points are best fit by a power-law function of E p ∝ E 0.459±0.018γ for the on-axis case (Figure 6a), with an intrinsic scatter of 0.263 ± 0.031.This agrees well with the recently updated Ghirlanda relation of E p ∝ E 0.44±0.07γ for a sample of 55 observed GRBs (Wang et al. 2018).For the off-axis case (Figure 6b), we find E p ∝ E 0.269±0.008γ , with a larger scatter of 0.272 ± 0.029.This indicates that the Ghirlanda relation may be largely connected with the Amati relation, but compared with the latter, the slope of the Ghirlanda relation is slightly flatter, and the data points are also obviously more scattered.It is difficult task to compare the off-axis Ghirlanda relation with observations because it needs the precise determination of both the beaming angle and the viewing angle for a number of off-axis GRBs, which themselves are rather dim.

Comparison with observations
In Section 5 we illustrated the simulated E p − E iso and E p − L p relations for both on-axis and off-axis cases.We now proceed to compare our results with observational data.
A sample containing 172 long GRBs was used for this purpose, 162 of which were collected from the previous study by Demianski et al. (2017) and Nava et al. (2012).E p and E iso parameters are available for all these 162 events.The L p parameter is available for only 45 GRBs (Nava et al. 2012).Additionally, we collected another 10 bursts that were argued to be possible outliers of the normal Amati relation in previous studies.It has been suggested that they might be off-axis GRBs (Yamazaki et al. 2003b;Ramirez-Ruiz et al. 2005).The observational data of these 10 events are listed in Table 2.
All the observed GRB samples are plotted on the E p − E iso and E p − L p planes in Figure 7.Most of the GRBs are well consistent with our on-axis results, thus they both follow the Amati relation and the Yonetoku relation.
On the E p − E iso plane (Figure 7a), four of the ten possible outliers are located between our on-and off-axis lines.It is thus difficult to judge whether these events are on-axis or off-axis GRBs simply from the Amati relation.However, the remaining six bursts can only be matched by our off-axis line.They should indeed be off-axis events.Figure 4c and 4d clearly show that most off-axis GRBs should be LLGRBs with E iso < 10 48 erg, with only a small portion falling in the range of 10 48 − 10 52 erg (see also Figure 5a).This is easy to understand.The input parameters should take some extreme values for an off-axis GRB to be strong enough.Especially these off-axis events should generally have a low θ obs − θ jet value (see Equation 27), which means that they are still slightly off-axis.In Figure 7a, only a small number of GRBs are outliers of the on-axis Amati relation.This low event rate is also consistent with our simulation results.
For the E p − L p correlation (Figure 7b), we note that almost all the ten possible outliers clearly deviate from the on-axis Yonetoku relation.From combining the Amati relation and the Yonetoku relation, we therefore argue that the ambiguous sample should also be off-axis GRBs.Furthermore, Figure 7 seems to indicate that the Yonetoku relation may be a better tool for distinguishing between on-and off-axis GRBs.

Conclusions and discussion
The prompt emission of GRBs was studied for both on-and off-axis cases.Especially the three prompt emission parameters E p , E iso , and L p were considered.Their dependence on the input model parameters was obtained via both analytical derivations and numerical simulations, which are consistent with each other.We confirmed that E p , E iso , and L p are independent of θ obs and θ jet as long as our line of sight is within the homogeneous jet cone (i.e., the on-axis cases).However, they strongly depend on the value of θ obs − θ jet when the jet is viewed offaxis.Additionally, their dependence on γ is very different for onand off-axis cases, as shown in Equations 26 and 27.Through Monte Carlo simulations, we found that the Amati relation is E p ∝ E 0.495±0.015iso in on-axis cases.Correspondingly, the Yonetoku relation is E p ∝ L 0.511±0.016p .On the other hand, in off-axis cases, the Amati relation is E p ∝ E 0.282±0.006iso and the Yonetoku relation is E p ∝ L 0.206±0.006p .The simulated samples were directly compared with the observational samples.We found that they are well consistent with each other in the slopes and intrinsic scatters.The E p − E γ relation was also tested and was found to be E p ∝ E 0.459±0.018γ for on-axis GRBs and E p ∝ E 0.269±0.008γ for off-axis GRBs.
We have focused on the empirical correlations of long GRBs in the above analysis.Here we present some discussion of short GRBs.Generally speaking, long GRBs usually contain many pulses (tens or even up to hundreds) in their light curves, while short GRBs contain much fewer pulses.For short GRBs, we can therefore perform a similar analysis, but only when the number of pulses considered is reduced significantly.In our model, the total energy release (E iso ) is proportional to the number of pulses, while the peak luminosity (L p ) and the spectral parameter E p are nearly irrelevant.In other words, a short GRB will have a much smaller E iso , but its L p and E p are largely unchanged.As a result, short GRBs will follow a parallel track on the E p − E iso plane as compared with long GRBs, the only difference is that their E iso are significantly smaller.Alternatively speaking, the power-law index of the Amati relation should be the same for short GRBs and long GRBs, consistent with currently available observations (Zhang et al. 2009(Zhang et al. , 2018)).As for the Yonetoku relation on the E p − L p plane, we argue that it should be almost identical for both long and short GRBs.However, the number of currently available short GRBs that is suitable for exploring the Yonetoku A simple top-hat jet geometry was adopted in this study.It should be noted that the GRB outflow might be a structured jet (Ioka & Nakamura 2018).A recent study suggested that E p and the flux will not change much, regardless of the choice of a tophat jet or a not-too-complicated structured jet (Farinelli et al. 2021).However, in some more complicated scenarios, structured jets may further have an angle-dependent energy density (Lamb et al. 2021).They may be choked jets or even in jet-cocoon systems (Ioka & Nakamura 2018;Mooley et al. 2018;Troja et al. 2019).The empirical correlations of GRBs in these situations are beyond the scope of this study and could be further considered in the future.
In our calculations, the prompt emission was assumed to be due to nonthermal radiation from internal shocks.Other mechanisms such as photosphere radiation and magnetic reconnection above the photosphere may also contribute to the observed flux of GRBs.Generally, there are two types of photosphere models: nondissipative photosphere models, and dissipative photosphere models.The former can produce a narrow quasi-Planckian component and could account for the thermal component observed in some GRBs (Pe'er 2008;Beloborodov 2011), although this component does not play a dominant role in most GRBs (Guiriec et al. 2011;Axelsson et al. 2012).In the dissipative photosphere models, various subphotosphere dissipation processes (Rees & Mészáros 2005), especially Comptonization, may affect the spectrum (Thompson 1994;Pe'er et al. 2006;Pe'er & Ryde 2011;Veres & Mészáros 2012), which may lead to a nonthermal spectrum.Some authors have even used this model to explain the Band spectrum of GRBs (Thompson 1994;Beloborodov 2013;Lundman et al. 2013).Reconnection above the photosphere has also been considered, especially in the view of the ICMART model, focusing on Poynting-flux-dominated outflows (Zhang & Yan 2011;Zhang & Zhang 2014).It is still unclear how commonly photosphere radiation and magnetic reconnection are involved in the prompt phase of GRBs.If these radiation mechanisms are included, the derivation will become much more complicated, but it deserves a trial in the future.
The Amati and Yonetoku relations are essential for us to understand the nature of GRBs.They can also help us probe the high-redshift universe (Xu et al. 2021;Hu et al. 2021;Zhao & Xia 2022;Jia et al. 2022;Deng et al. 2023).Our study shows that they are due to the ultra-relativistic effect of highly collimated jets.To be more specific, they mainly result from the variation in b a b a Fig. 5. Distribution of simulated off-axis bursts as compared with that of on-axis ones.a, All mock bursts on the E p − E iso plane.b, All mock bursts on the E p − L p plane.The dash-dotted line and the dotted line correspond to the best-fit result for on-and off-axis samples, respectively.The dashed lines represent the corresponding 3σ confidence level.Most of the off-axis bursts are too weak to appear in the two panels here.the Lorentz factor in different events.It may provide useful insights for a better understanding of these empirical correlations.

Fig. 1 .
Fig. 1.Geometry of our jet model.The upper panel shows the case for an on-axis jet, where θ jet ≥ θ obs .The lower panel shows an off-axis jet with θ jet < θ obs .

Fig. 2 .Fig. 3 .
Fig. 2. Time-integrated spectrum for different γ values (marked in the figure).Panel a Spectrum of on-axis cases.Panel b Spectrum of off-axis cases.The black star on each curve marks the position of the peak energy (E p ).

Fig. 4 .
Fig. 4. Monte Carlo simulation results.(a,b), Distribution of simulated bursts on the E p − E iso plane and the E p − L p plane for the on-axis cases.(c,d),Distribution of simulated bursts on the E p − E iso plane and the E p − L p plane for off-axis cases.The solid lines show the best-fitting results for the simulated data points, and the dashed lines represent the 3σ confidence level.To conclude, for on-axis bursts, the best-fit Amati relation is E p ∝ E 0.495±0.015iso

Fig. 6 .Fig. 7 .
Fig. 6.Distribution of the mock GRBs on the E p − E γ plane for the on-axis case.Panel a On-axis case.Panel b Off-axis case.The solid line shows the best-fit result for the simulated data points, and the dashed lines represent the 3σ confidence level.The on-axis Ghirlanda relation is E p ∝ E 0.459±0.018γ

Table 1 .
Standard values assumed for the eight input parameters