Interpreting time-integrated polarization data of gamma-ray burst prompt emission

Aims. With the accumulation of polarization data in the gamma-ray burst (GRB) prompt phase, polarization models can be tested. Methods. We predicted the time-integrated polarizations of 37 GRBs with polarization observation. We used their observed spectral parameters to do this. In the model, the emission mechanism is synchrotron radiation, and the magnetic field configuration in the emission region was assumed to be large-scale ordered. Therefore, the predicted polarization degrees (PDs) are upper limits. Results. For most GRBs detected by the Gamma-ray Burst Polarimeter (GAP), POLAR, and AstroSat, the predicted PD can match the corresponding observed PD. Hence the synchrotron-emission model in a large-scale ordered magnetic field can interpret both the moderately low PDs ($\sim10\%$) detected by POLAR and relatively high PDs ($\sim45\%$) observed by GAP and AstroSat well. Therefore, the magnetic fields in these GRB prompt phases or at least during the peak times are dominated by the ordered component. However, the predicted PDs of GRB 110721A observed by GAP and GRB 180427A observed by AstroSat are both lower than the observed values. Because the synchrotron emission in an ordered magnetic field predicts the upper-limit of the PD for the synchrotron-emission models, PD observations of the two bursts challenge the synchrotron-emission model. Then we predict the PDs of the High-energy Polarimetry Detector (HPD) and Low-energy Polarimetry Detector (LPD) on board the upcoming POLAR-2. In the synchrotron-emission models, the concentrated PD values of the GRBs detected by HPD will be higher than the LPD, which might be different from the predictions of the dissipative photosphere model. Therefore, more accurate multiband polarization observations are highly desired to test models of the GRB prompt phase.


Introduction
Gamma-ray bursts (GRBs) are the most violent high-energy explosions in the Universe.GRBs were divided into two categories, long and short GRBs, based on a rough duration separation of about 2 s.GRB spectra are nonthermal, which can be described by a broken power law with a smooth joint, known as the Band function (Band et al. 1993).The spectrum integrated over the GRB duration can empirically be described by a function with a peak in the ν f ν spectrum, and the peak energy is defined as E p,obs .For the low-energy spectral index α, the typical value for long GRBs is α ∼ −0.92, and short GRBs have a harder low-energy spectral index α ∼ −0.50 (Nava et al. 2011).
Gamma-ray polarization measurements of the prompt emission of GRBs have profound implications for our understanding of the unknown magnetic field configuration and emission mechanism of GRBs prompt emission.With the development of polarimetry, increasingly more GRBs have been measured and can be used for a statistical analysis.Therefore, constraints on the underlying models can be provided (Toma et al. 2009).For the GRB prompt phase, there are two possible emission mechanism, synchrotron radiation and inverse Compton scattering (Chand et al. 2018;Lazzati et al. 2004).Although several thousands of GRBs have been observed to date, few of these have reported polarization detections.The polarization degrees (PDs) of GRB prompt emissions in the measurements so far vary strongly.
The Gamma-ray Burst Polarimeter (GAP) has observed PD values of GRBs 100826A, 110301A, and 110721A, which suggests that GRB prompt emissions are highly polarized (Yonetoku et al. 2011(Yonetoku et al. , 2012)).Then an increasing number of polarimeters became operational.Chattopadhyay et al. (2022) published their updated polarization observational results of 20 GRBs recently, which are the brightest GRBs detected by the Cadmium Zinc Telluride Imager (CZTI) on board AstroSat.The renewed AstroSat data show that most of the bright GRBs are relatively highly polarized (with a typical PD value of ∼45%) in the energy range of 100 keV−600 keV, different from that of their former results with the high polarizations (typical PD was around 60%) in energy range of 100 keV−350 keV (Chattopadhyay et al. 2019;Chand et al. 2019;Gupta et al. 2022).
POLAR is a polarimeter with a comparable energy range as CZTI, which was launched as part of the China Tiangong-2 space laboratory in September 2016.The detection energy range of POLAR is 50 keV−500 keV.During its approximately six months of operation, a total of 55 GRBs were detected (Xiong et al. 2017).Polarization measurements of 5 of these 55 GRBs were reported first, and the results show that they are less polarized than predicted by some popular models (Zhang et al. 2019).Moderate levels of linear polarization were also found in subsequent reports, and the polarization measurements of 9 other GRBs were published next (Kole et al. 2020).Despite the great efforts that have been made in gamma-ray polarimetry, there are still large errors in the

−
Notes. (* ) The spectral parameters are obtained from Kole et al. (2020), who performed a joint fit using an external spectrum with POLAR data based on the Multi-Mission Maximum Likelihood (3ML) framework (Vianello et al. 2015).
current data, which allow us to present only preliminary constraints on the various models of the GRB prompt phase.It is encouraging that more detailed polarization measurements will become available from forthcoming missions such as POLAR-2 (de Angelis & Polar-2 Collaboration 2022), which will help us to understand the magnetic field configuration and emission mechanism of GRBs.
In this paper, we have numerically calculated (Toma et al. 2009) the ranges of theoretical PDs of 37 GRBs detected by GAP, AstroSat, and POLAR based on the values of the observed spectral parameters.The paper is arranged as follows.In Sect. 2 we present our data.The model and numerical results are described in Sect.3. Finally, we give our conclusions and discussion in Sect. 4. Yonetoku et al. (2011Yonetoku et al. ( , 2012) ) reported polarization observations of the prompt emission of GRB 100826A, GRB 110301A, and GRB 110721A with GAP.For GRB 100826A, an averaged polarization of 27 ± 11% with a confidence level of 99.4% (2.9σ) was reported, with systematic errors being considered for the first time in their analysis.For GRB 110301A and GRB 110721A, the observed linear polarizations are 70 ± 22% and 84 +16 −28 % with confidence levels of 3.7σ and 3.3σ, respectively.Berger (2011) reported a redshift value of 0.382 for GRB 110721A.The spectral parameters used in our calculations for all three GRBs are from the Fermi-GBM catalog in the energy range of 50 keV−300 keV 1 (von Kienlin et al. 2014,  2020; Gruber et al. 2014; Narayana Bhat et al. 2016), and are presented in Table 1.

Data list
Recently, Kole et al. (2020) published a detailed polarization catalog reporting the polarization properties of 14 GRBs observed with POLAR.We searched the parameters of the Band function for all of these GRBs and list them in Table 2 along with the instruments that provide them.Among these spectral parameters, all those from Konus were measured in the energy range of 20 keV−15 MeV (Frederiks et al. 2016;Tsvetkova et al. 2017;Svinkin et al. 2017), while those from Fermi-GBM were obtained in the energy range of 50 keV−300 keV (von Kienlin & Roberts 2017; Roberts & Mailyan 2017;Stanbro et al. 2017).
The updated polarization measurements and the corresponding spectral parameters of 20 GRBs observed by the CZTI on board AstroSat were also reported by Chattopadhyay et al. (2022) recently.In Table 3 we list the detailed polarization information and spectral properties for them (Chattopadhyay et al. 2022).In addition, the redshift values have been reported for 4 of 20 GRBs (Malesani et al. 2016;de Ugarte Postigo et al. 2017;Gupta et al. 2022;D'Avanzo et al. 2018), which provides more precise parameters for our calculations.Notes. (a) Fitting this spectrum with the Band function only presents a lower limit on β s of 1.48.

Model and numerical results
An ultrarelativistic jet is assumed to be an optically thin shell to gamma-rays with an emitting region of radius r, located at redshift z, and a source with a luminosity distance d L .Its fluence can be expressed as follows (Toma et al. 2009;Granot et al. 1999;Woods & Loeb 1999;Ioka & Nakamura 2001): In the above equation, θ V is the viewing angle of the observer, θ j is the half-opening angle of the jet, and θ is the angle between the line of sight and the local direction of the fluid velocity.The physical quantities that are primed and unprimed are in the comoving and observer frame, respectively.For example, ν = ν obs (1 + z)γ(1 − β 0 cos θ) is the observational frequency in comoving frame, with the bulk Lorentz factor γ and the velocity of jet β 0 in units of the speed of light.ν obs is the observational frequency in the observer frame.ν 1 and ν 2 are the energy ranges of the corresponding detectors (e.g., ν 1 = 50 keV and ν 2 = 500 keV for POLAR).φ is the angle between the projection of the jet axis and the projection of the local fluid velocity direction on the sky plane.More information about ∆φ can be obtained from Toma et al. (2009).E p,obs can be converted into the comoving frame by E p = E p,obs (1 + z)/2/γ.We adopted the following form for the spectrum of GRB prompt emission described by the Band function (Band et al. 1993): (2) α s and β s are the low-energy and high-energy spectral indices, respectively.ν 0 = E p /h is the comoving break energy of the Band spectrum.h is the Planck constant.In this paper, α s and β s are the spectral indices of the flux density F ν .In our calculation, the source was assumed to be at a redshift of 1 unless its redshift value has been reported.We assumed an aligned largescale ordered magnetic field in the emission region with an orientation of δ = π/6 (Lan et al. 2016).Other fixed parameters are θ j = 0.1 rad, θ V = 0 rad, and γ = 100.We then calculated the PDs of the GRBs with the polarization observation, using the observed spectral parameters (including α s , β s , and E p,obs ) and redshift values as well as the energy range of the polarimeters.In general, the calculated PD (PD cal ) of a GRB consists of a typical value and its upper and lower limits.In our calculations, we used the redshift value, the detector energy range, and the typical values of E p,obs , α s , and β s to calculate a typical value of PD cal .For the same GRB (i.e., the redshift value and the upper and lower limits of the detector energy range are fixed), the upper limit of its PD cal was taken when its α s and β s took the maximum values and E p,obs took the minimum value; conversely, the minimum values of α s and β s and the maximum value of E p,obs determine the lower limit of the PD cal .We compare the calculated PDs (PD cal ) and the observed PDs (PD obs ) in Figs.1-3.
Figure 1 shows a comparison of PD cal and PD obs observed with GAP.The energy ranges of GAP (for polarization observations) and of the Fermi-GBM (for spectra observations) overlap exactly (50 keV−300 keV).For GRB 100826A detected with GAP, the polarization evolutions of this burst were simulated with the collision-induced magnetic reconnection model (Deng et al. 2016), and the results can reproduce the A160, page 3 of 6 time-resolved polarizations, especially the 90 • polarization angle (PA) change between the two pulses.The observed PD of GRB 110721A is larger than the predicted one.We also calculated the PD cal ranges and compared them with the PD obs observed by POLAR, as shown in Fig. 2. PD cal for 10 of the 13 GRBs in the light blue region overlaps with their corresponding PD obs .The PD cal of the remaining 4 GRBs is significantly higher than the PD obs ranges.All these PD cal show higher polarization degrees.
In Fig. 3 we numerically calculated the PD cal ranges of the GRBs observed by the AstroSat and found that the results match most of the observations, with a distribution around ∼40%.The only burst with observed PD larger than the predicted value is GRB 180427A.Our integrated energy range of Stokes parameters for AstroSat is 100 keV−600 keV.The range of Π 0 is [0, 1].This requires that the spectral index (α s or β s ) should be higher than −1 according to the local polarization as shown below (Toma et al. 2009).
The PDs of the POLAR bursts are concentrated around 10%, while they are around 40%−50% for AstroSat bursts.To interpret this discrepancy, we plot the spectral indices with peak energy in Fig. 4. The typical values of high-energy spectral indices are similar for POLAR bursts and AstroSat bursts.However, the typical value of low-energy spectral index is higher for AstroSat bursts (typically α s ∼ 0.0) than the POLAR bursts (typically α s ∼ −0.5), resulting in a higher PD cal for AstroSat bursts.And the integrated energy range (100 keV−600 keV) of AstroSat bursts, compared with that (50 keV−500 keV) of POLAR bursts, is shifted to the higher energy range.For the bursts with similar spectral parameters the contributions from the high-energy photons (with larger local PD) will be larger for AstroSat burst, which will lead to a higher energy-integrated PD cal .These might be the main reasons for this discrepancy.
Because the energy range of the polarimeter also affects the observed polarization properties, we numerically predict the PDs of the long and short GRBs measured by two detectors, Low-energy Polarimetry Detector (LPD) and Highenergy Polarimetry Detector (HPD) (whose energy ranges are 2 keV−30 keV and 30 keV−800 keV, respectively) on board POLAR-2 (de Angelis & Polar-2 Collaboration 2022), based on the typical values and distribution of their spectral parameters (Nava et al. 2011;Preece et al. 2000).We present the results in Fig. 5, where the gray area and the light blue area denote the energy ranges of LPD and HPD, respectively.The typical PD values of GRBs detected by LPD and HPD are shown as black diamonds and red points.Because the typical PD values for different detectors are calculated under the typical spectral parameters and the number of GRBs with typical spectral parameters is maximum, the observed PD value with the maximum number of GRBs for each detector is concentrated around the typical PD values predicted by the model.

Discussion and conclusion
Polarization properties of GRBs are essential for diagnosing the magnetic field configuration and the geometry of the emission A160, page 4 of 6 Long GRBs Short GRBs PD (%) region and the observation.We used the observed energy spectrum to calculate the corresponding GRB polarization properties within the synchrotron-emission model and compared them with the observed time-integrated PDs.In our model, we used the large-scale aligned magnetic field in the emission region.Therefore, the predicted PDs give upper limits for the synchrotronemission models.
For GAP, POLAR, and AstroSat the predicted PDs of the model can match most of the corresponding observed PDs, indicating that in the GRB prompt phase or at least during the peak time of the burst, the magnetic field configuration is approximately large-scale ordered with ξ B > 1 (Lan et al. 2019(Lan et al. , 2021)).The large-scale ordered magnetic field in the GRB emission region may originate from its central engine.In the scenario of the internal shock in a fireball (Paczynski & Xu 1994;Rees & Meszaros 1994), the magnetic field may be mixed with a low ξ B value (ξ B < 1), so that this model is not favored by the current PD observations.For the internal shock with an ordered magnetic field (Fan et al. 2004), the magnetization parameter σ is required to be smaller than 1.The observed PD values require that it cannot be too small, however, otherwise, turbulence will develop and destroy the ordered magnetic field (Deng et al. 2017).For the ICMART model (Zhang & Yan 2011), the magnetic field becomes less ordered with the magnetic reconnection during the burst (i.e., a decrease in ξ B ).The observed data indicate that at the peak time of these bursts, the ξ B values of the magnetic fields are still higher than 1 (i.e., the magnetic fields are dominated by the ordered component at the peak times of the bursts).
For POLAR, 10% of the observed PDs can also be interpreted as synchrotron emission in an ordered magnetic field with a small low-energy spectral index.However, PD cal of four GRBs is still higher than PD obs .For these four GRBs, the magnetic field configurations in the emission regions may be mixed (Lan et al. 2019), or the PAs show rotation or an abrupt change during the bursts.Future time-resolved polarization observations will enable us to distinguish the two scenarios.For AstroSat, the predicted PDs are concentrated around 40% and can interpret the measurements of all GRBs except one (GRB 180427A).There is a discrepancy between the moderately low PDs (∼10%) detected with POLAR and the relatively high PDs (about 40%−50%) observed with AstroSat.The main reasons for this difference may originate from both the higher low-energy spectral indices and higher integrated energy range for AstroSat bursts.
The PD data of GRB 180427A detected by AstroSat and of GRB 110721A detected by GAP are both higher than the predicted values.Therefore the two PD observations challenge the models invoking synchrotron radiation in an ordered magnetic field.Because the synchrotron radiation in an ordered magnetic field gives the upper limit of the PD of the synchrotron-emission models with a mixed magnetic field for on-axis observations (Lan & Dai 2020), and the GRBs selected for the polarization analysis are usually bright (indicating on-axis observations), the PD data of GRB 180427A and GRB 110721A finally challenge the synchrotron-emission models for GRB prompt phase.
With co-observations of the HPD and LPD on board POLAR-2 (de Angelis & Polar-2 Collaboration 2022), the polarization spectrum will be obtained in the near future.We predict that the concentrated PD values of the GRBs detected by the HPD will be higher than the LPD in the synchrotronemission model.A reversed polarization spectrum was predicted by the dissipative photosphere model, however, that is, the concentrated PD values detected by the HPD will be lower than the LPD.The two models can be tested in this way with polarization observations of POLAR-2.The emission mechanism at the highenergy gamma-ray band is multiple inverse-Compton scattering for the dissipative photosphere model (Lundman et al. 2018), which is different from the synchrotron-emission model.With the observations of the polarization spectrum of the POLAR-2, the emission mechanism in the high-energy gamma-ray band can therefore be determined.A160, page 5 of 6

Fig. 1 .
Fig. 1.Ranges of PD cal and PD obs observed by GAP.The black diamond and red points correspond to the observed and calculated PD values, respectively.The area in light blue in the figure indicates the energy range of GAP, which is 50 keV−300 keV.

Fig. 4 .
Fig. 4. Comparison of the spectral parameters of the bursts observed with POLAR and AstroSat.The black squares and circles represent lowand high-energy spectral indices for bursts detected by POLAR, the red diamonds and stars correspond to low-and high-energy spectral indices for bursts observed by AstroSat.

Fig. 5 .
Fig. 5. Comparison of polarization measurements of the long and short GRBs observed with HPD and LPD.The concentrated PD values of LPD and the HPD are shown with black diamonds and red points.

Table 1 .
Spectral parameters and polarization properties of the three GRBs observed with GAP.

Table 2 .
Spectral parameters and polarization properties of the 14 GRBs observed with POLAR.

Table 3 .
Spectral parameters and polarization properties of the 20 GRBs observed with AstroSat.