Gamma-ray bursts are the strongest explosions in the Universe since the Big Bang, believed to be produced either in forming black holes at the end of massive star evolution [1, 2, 3] or merging of compact objects [4]. Spectral and timing properties of gamma-ray bursts suggest that the observed bright gamma-rays are produced in the most relativistic jets in the Universe [4]; however, the physical properties, especially the structure and magnetic topologies in the jets are still not well known, despite several decades of studies. It is widely believed that precise measurements of the polarization properties of gamma-ray bursts should provide crucial information on the highly relativistic jets [5]. As a result there have been many reports of gamma-ray burst polarization measurements with diverse results, see [1], however many such measurements suffered from substantial uncertainties, mostly systematic [7, and references therein]. After the first successful measurements by the GAP and COSI instruments [2, 3, 4], here we report a statistically meaningful sample of precise polarization measurements, obtained with the dedicated gamma-ray burst polarimeter, POLAR onboard China's Tiangong-2 spacelab. Our results suggest that the gamma-ray emission is at most polarized at a level lower than some popular models have predicted; although our results also show intrapulse evolution of the polarization angle. This indicates that the low polarization degrees could be due to an evolving polarization angle during a gamma-ray burst.POLAR is a dedicated Gamma-ray Burst (GRB) polarization detection experiment onboard China's Tiangong-2 spacelab [11], launched on Sept. 15th, 2016 and stopped operation on March 31, 2017. POLAR detected 55 GRBs with high significance. In order to make statistically significant GRB polarization measurements and yet with negligible systematic errors, we select a subsample of five GRBs for detailed analysis of their polarization properties; please refer to the supplementary information (SI) for the sample selection criteria and the properties of the five selected GRBs. We employ a straight forward χ 2 based analysis, similar to that successfully employed in [4], to study the polarization properties of the five GRBs, while a Bayesian method is employed to accurately determine the credible regions of the measurements. The studies rely on extensive ground and in-orbit calibration data and Monte-Carlo simulations matching the calibration data [12, 13]. Please refer to the methods section for details of the methodology and analysis.In Figure 1, we show the measured modulation curves of the five GRBs integrated over the whole GRB duration, together with the best fitting simulated modulation curves from linear polarization and fitting residuals. All fittings are statistically acceptable with no significant systematic deviations. In Figure 2, we show the 2-D posterior distributions of the five GRBs, i.e., the posterior probability as functions of both polarization angle (PA) and degree (PD). Clearly the measured P...

Context. Despite over 50 years of research, many open questions remain about the origin and nature of gamma-ray bursts (GRBs). Linear polarization measurements of the prompt emission of these extreme phenomena have long been thought to be key to answering a range of these questions. The POLAR detector was designed to produce the first set of detailed and reliable linear polarization measurements in the 50 − 500 keV energy range. During late 2016 and early 2017, POLAR detected a total of 55 GRBs. The analysis results of 5 of these GRBs have been reported, and were found to be consistent with a low or unpolarized flux. However, previous reports by other collaborations found high levels of linear polarization, including some as high as 90%. Aims. We study the linear polarization for the 14 GRBs observed by POLAR for which statistically robust inferences are possible. Additionally, time-resolved polarization studies are performed on GRBs with sufficient apparent flux. Methods. A publicly available polarization analysis tool, developed within the Multi-Mission Maximum Likelihood framework (3ML), was used to produce statistically robust results. The method allows spectral and polarimetric data from POLAR to be combined with spectral data from the Fermi Gamma-ray Burst Monitor (Fermi-GBM) and the Neil Gehrels Swift Observatory (hereafter Swift) to jointly model the spectral and polarimetric parameters. Results. The time-integrated analysis finds all results to be compatible with low or zero polarization with the caveat that, when timeresolved analysis is possible within individual pulses, we observe moderate linear polarization with a rapidly changing polarization angle. Therefore, time-integrated polarization results, while pointing to lower polarization, are potentially an artifact of summing over the changing polarization signal and thus washing out the true moderate polarization. We therefore caution against overinterpretation of any time-integrated results inferred herein and encourage the community to wait for more detailed polarization measurements from forthcoming missions such as POLAR-2 and LEAP.

Context. Simultaneous γ-ray measurements of γ-ray burst spectra and polarization offer a unique way to determine the underlying emission mechanism(s) in these objects, as well as probing the particle acceleration mechanism(s) that lead to the observed γ-ray emission. Aims. We examine the jointly observed data from POLAR and Fermi-GBM of GRB 170114A to determine its spectral and polarization properties, and seek to understand the emission processes that generate these observations. We aim to develop an extensible and statistically sound framework for these types of measurements applicable to other instruments. Methods. We leveraged the existing 3ML analysis framework to develop a new analysis pipeline for simultaneously modeling the spectral and polarization data. We derived the proper Poisson likelihood for γ-ray polarization measurements in the presence of background. The developed framework is publicly available for similar measurements with other γ-ray polarimeters. The data are analyzed within a Bayesian probabilistic context and the spectral data from both instruments are simultaneously modeled with a physical, numerical synchrotron code. Results. The spectral modeling of the data is consistent with a synchrotron photon model as has been found in a majority of similarly analyzed single-pulse gamma-ray bursts. The polarization results reveal a slight trend of growing polarization in time reaching values of ∼30% at the temporal peak of the emission. We also observed that the polarization angle evolves with time throughout the emission. These results suggest a synchrotron origin of the emission but further observations of many GRBs are required to verify these evolutionary trends. Furthermore, we encourage the development of time-resolved polarization models for the prompt emission of gamma-ray bursts as the current models are not predictive enough to enable a full modeling of our current data.