Search

Gamma-Ray Bursts are among the most powerful events in the Universe. Despite half a century of observations of these transient sources, many open questions remain about their nature. Polarization measurements of the GRB prompt emission have long been theorized to be able to answer most of these questions. With the aim of characterizing the polarization of these prompt emissions, a compact Compton polarimeter, called POLAR, has been launched to space in September 2016. Time integrated polarization analysis of the POLAR GRB catalog have shown that the prompt emission is lowly polarized or fully unpolarized. However, time resolved analysis depicted strong hints of an evolving polarization angle within single pulses, washing out the polarization degree in time integrated analyses. Here we will for the first time present energy resolved polarization measurements with the POLAR data. The novel analysis, performed on several GRBs, will provide new insights and alter our understanding of GRB polarization. The analysis was performed using the 3ML framework to fit polarization parameters versus energy in parallel to the spectral parameters. Although limited by statistics, the results could provide a very relevant input to disentangle between existing theoretical models. In order to gather more statistics per GRB and perform joint time and energy resolved analysis, a successor instrument, called POLAR-2, is under development with a launch window early 2025 to the CSS. After presenting the first energy resolved polarization results of the POLAR mission, we will present the prospects for such measurements with the upcoming POLAR-2 mission., Comment: Proceeding from the 38th International Cosmic Ray Conference (ICRC2023), 9 pages, 6 figures

The POLAR-2 Gamma-Ray Burst (GRB) Polarimetry mission is a follow-up to the successful POLAR mission. POLAR collected six months of data in 2016-2017 on board the Tiangong-2 Chinese Space laboratory. From a polarization study on 14 GRBs, POLAR measured an overall low polarization and a hint for an unexpected complexity in the time evolution of polarization during GRBs. Energy-dependent measurements of the GRB polarization will be presented by N. de Angelis in GA21-09 (August 2nd). These results demonstrate the need for measurements with significantly improved accuracy. Moreover, the recent discovery of gravitational waves and their connection to GRBs justifies a high-precision GRB polarimeter that can provide both high-precision polarimetry and detection of very faint GRBs. The POLAR-2 polarimeter is based on the same Compton scattering measurement principle as POLAR, but with an extended energy range and an order of magnitude increase in total effective area for polarized events. Proposed and developed by a joint effort of Switzerland, China, Poland and Germany, the device was selected for installation on the China Space Station and is scheduled to start operation for at least 2 years in 2025., Comment: Proceeding from the 38th International Cosmic Ray Conference (ICRC2023), 8 pages, 6 figures

`ronswanson` provides a simple-to-use framework for building so-called table or template models for `astromodels`the modeling package for multi-messenger astrophysical data-analysis framework, `3ML`. With `astromodels` and `3ML` one can build the interpolation table of a physical model result of an expensive computer simulation. This then enables efficient reevaluation of the model while, for example, fitting it to a dataset. While `3ML` and `astromodels` provide factories for building table models, the construction of pipelines for models that must be run on high-performance computing (HPC) systems can be cumbersome. `ronswanson` removes this complexity with a simple, reproducible templating system. Users can easily prototype their pipeline on multi-core workstations and then switch to a multi-node HPC system. `ronswanson` automatically generates the required `Python` and `SLURM` scripts to scale the execution of `3ML` with `astromodel`'s table models on an HPC system., Comment: 3 pages. Accepted for publication in JOSS

The last decade has brought about a profound transformation in multimessenger science. Ten years ago, facilities had been built or were under construction that would eventually discover the nature of objects in our universe could be detected through multiple messengers. Nonetheless, multimessenger science was hardly more than a dream. The rewards for our foresight were finally realized through IceCube's discovery of the diffuse astrophysical neutrino flux, the first observation of gravitational waves by LIGO, and the first joint detections in gravitational waves and photons and in neutrinos and photons. Today we live in the dawn of the multimessenger era. The successes of the multimessenger campaigns of the last decade have pushed multimessenger science to the forefront of priority science areas in both the particle physics and the astrophysics communities. Multimessenger science provides new methods of testing fundamental theories about the nature of matter and energy, particularly in conditions that are not reproducible on Earth. This white paper will present the science and facilities that will provide opportunities for the particle physics community renew its commitment and maintain its leadership in multimessenger science., Comment: 174 pages, 12 figures. Contribution to Snowmass 2021. Solicited white paper from CF07. Comments and endorsers welcome. Still accepting contributions (contact editors)

Gamma-rays, the most energetic photons, carry information from the far reaches of extragalactic space with minimal interaction or loss of information. They bring messages about particle acceleration in environments so extreme they cannot be reproduced on earth for a closer look. Gamma-ray astrophysics is so complementary with collider work that particle physicists and astroparticle physicists are often one in the same. Gamma-ray instruments, especially the Fermi Gamma-ray Space Telescope, have been pivotal in major multi-messenger discoveries over the past decade. There is presently a great deal of interest and scientific expertise available to push forward new technologies, to plan and build space- and ground-based gamma-ray facilities, and to build multi-messenger networks with gamma rays at their core. It is therefore concerning that before the community comes together for planning exercises again, much of that infrastructure could be lost to a lack of long-term planning for support of gamma-ray astrophysics. Gamma-rays with energies from the MeV to the EeV band are therefore central to multiwavelength and multi-messenger studies to everything from astroparticle physics with compact objects, to dark matter studies with diffuse large scale structure. These goals and new discoveries have generated a wave of new gamma-ray facility proposals and programs. This paper highlights new and proposed gamma-ray technologies and facilities that have each been designed to address specific needs in the measurement of extreme astrophysical sources that probe some of the most pressing questions in fundamental physics for the next decade. The proposed instrumentation would also address the priorities laid out in the recent Astro2020 Decadal Survey, a complementary study by the astrophysics community that provides opportunities also relevant to Snowmass., Comment: Contribution to Snowmass 2021

INTEGRAL/SPI is a coded mask instrument observing since 2002 in the keV to MeV energy range, which covers the peak of the $\nu F\nu$ spectrum of most Gamma-Ray Bursts (GRBs). Since its launch in 2008, Fermi/GBM has been the primary instrument for analyzing GRBs in the energy range between $\approx$ 10 keV to $\approx$ 10 MeV. Herein, we show that SPI, covering a similar energy range, can give equivalently constraining results for some parameters if we use an advanced analysis method. Also, combining the data of both instruments reduces the allowed parameter space in spectral fits. The main advantage of SPI as compared to GBM is the energy resolution of $\approx$ 0.2\% at 1.3 MeV compared to $\approx$ 10\% for GBM. Therefore, SPI is an ideal instrument to precisely measure the curvature of the spectrum. This is important, as it has been shown in recent years that physical models rather than heuristic functions should be fit to GRB data to obtain better insights into their still unknown emission mechanism, and the curvature of the peak is unique to the different physical models. To fit physical models to SPI GRB data and get the maximal amount of information from the data, we developed a new open source analysis software {\tt PySPI}. We apply these new techniques to GRB 120711A in order to validate and showcase {\tt PySPI}'s capabilities. We show that {\tt PySPI} improves the analysis of SPI GRB data compared to the {\tt OSA} analysis. In addition, we demonstrate that the GBM and the SPI data of this GRB can be fitted well with a physical synchrotron model. This evinces that SPI can play an important role in GRB spectral model fitting., Comment: 11 pages, 13 figures, Accepted by A&A

Simulating a survey of fluxes and redshifts (distances) from an astrophysical population is a routine task. \texttt{popsynth} provides a generic, object-oriented framework to produce synthetic surveys from various distributions and luminosity functions, apply selection functions to the observed variables and store them in a portable (HDF5) format. Population synthesis routines can be constructed either using classes or from a serializable YAML format allowing flexibility and portability. Users can not only sample the luminosity and distance of the populations, but they can create auxiliary distributions for parameters which can have arbitrarily complex dependencies on one another. Thus, users can simulate complex astrophysical populations which can be used to calibrate analysis frameworks or quickly test ideas., Comment: Accepted in JOSS

Context. Gamma-ray bursts can be located via arrival time signal triangulation using gamma-ray detectors in orbit throughout the solar system. The classical approach based on cross-correlations of binned light curves ignores the Poisson nature of the time-series data, and is unable to model the full complexity of the problem. Aims. To present a statistically proper and robust GRB timing/triangulation algorithm as a modern update to the original procedures used for the Interplanetary Network (IPN). Methods. A hierarchical Bayesian forward model for the unknown temporal signal evolution is learned via random Fourier features (RFF) and fitted to each detector's time-series data with time-differences that correspond to GRB's position on the sky via the appropriate Poisson likelihood. Results. Our novel method can robustly estimate the position of a GRB as verified via simulations. The uncertainties generated by the method are robust and in many cases more precise compared to the classical method. Thus, we have a method that can become a valuable tool for gravitational wave follow-up. All software and analysis scripts are made publicly available here (https://github.com/grburgess/nazgul) for the purpose of replication., Comment: 18 pages, 30 figures. Submitted to A&A. Criticism gladly welcome

Despite over 50 years of research, many open questions remain about the origin and nature of GRBs. Polarization measurements of the prompt emission of these extreme phenomena have long been thought to be the key to answering a range of these questions. The POLAR detector was designed to produce the first set of detailed and reliable polarization measurements in an energy range of approximately 50-500 keV. During late 2016 and early 2017, POLAR detected a total of 55 GRBs. Analysis results of 5 of these GRBs have been reported in the past. The results were found to be consistent with a low or unpolarized flux. However, previous reports by other collaborations found high levels of polarization. We study the polarization for all 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. A publicly available polarization analysis tool, developed within the 3ML framework, was used to produce statistically robust results. The method allows to combine spectral and polarimetric data from POLAR with spectral data from the Fermi GBM and Swift-BAT to jointly model the spectral and polarimetric parameters. The time integrated analysis finds all results to be compatible with a low or zero polarization with the caveat that, when time-resolved analysis is possible within individual pulses, we observe moderate polarization with a rapidly changing polarization angle. Thus, 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. Therefore, we caution against over interpretation of any time-integrated results and encourage one to wait for more detailed polarization measurements from forthcoming missions such as POLAR-2 and LEAP., Comment: 43 pages, 78 figures, Accepted for publication in A&A

We present the first physically motivated background model for the Gamma-Ray Burst Monitor (GBM) onboard the Fermi satellite. Such a physically motivated background model has the potential to significantly improve the scientific output of Fermi/GBM, as it can be used to improve the background estimate for spectral analysis and localization of Gamma-Ray Bursts (GRBs) and other sources. Additionally, it can also lead to detections of new transient events, since long/weak or slowly rising ones do not activate one of the existing trigger algorithms. In this paper we show the derivation of such a physically motivated background model, which includes the modeling of the different background sources and the correct handling of the response of GBM. While the goal of the paper is to introduce the model rather than developing a transient search algorithm, we demonstrate the ability of the model to fit the background seen by GBM by showing four applications, namely (1) for a canonical GRB, (2) for the ultra-long GRB 091024, (3) for the V404 Cygni outburst in June 2015, and (4) the ultra-long GRB 130925A., Comment: Accepted for publication in A&A, 15 pages, 22 figures

The localizations of gamma-ray bursts (GRBs) detected with the Gamma-ray Burst Monitor (GBM) onboard the Fermi satellite are known to be affected by significant systematic errors of 3-15 degrees. This is primarily due to mismatch of the employed Band function templates and the actual GRB spectrum. This problem can be avoided by simultaneously fitting for the location and the spectrum of a GRB, as demonstrated with an advanced localization code, BALROG (arXiv:1610.07385). Here, we analyze in a systematic way a sample of 105 bright GBM-detected GRBs for which accurate reference localizations are available from the Swift observatory. We show that the remaining systematic error can be reduced to $\sim$1-2 degrees.

Simultaneous $\gamma$-ray measurements of gamma-ray burst (GRB) 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 $\gamma$-ray emission. Herein we examine the jointly-observed data from POLAR and 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. We leverage the existing 3ML analysis framework to develop a new analysis pipeline for simultaneously modeling the spectral and polarization data. We derive the proper Poisson likelihood for $\gamma$-ray polarization measurements in the presence of background. The developed framework is publicly available for similar measurements with other $\gamma$-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. 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 GRBs. The polarization results reveal a slight trend of growing polarization in time reaching values of ~30% at the temporal peak of the emission. Additionally, it is observed that the polarization angle evolves with time throughout the emission. These results hint at 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 GRBs as the current models are not predictive enough to enable a full modeling of our current data., Comment: 22 pages, 23 figures. Published to A&A

Gamma-ray bursts are the most energetic electromagnetic sources in the Universe. Their prompt gamma-ray radiation corresponds to an energy release of 1E42-1E47J. Fifty years after their discovery and several dedicated space-based instruments, the physical origin of this emission is still unknown. Synchrotron emission has been one of the early contenders but was criticized because spectral fits of empirical models (such as a smoothly-connected broken power law or a cut-off power law) suggest too hard a slope of the low-energy power law, violating the so-called synchrotron line-of-death. We perform time-resolved gamma-ray spectroscopy of single-peaked GRBs as measured with Fermi/GBM. We demonstrate that idealized synchrotron emission, when properly incorporating time-dependent cooling of the electrons, is capable of fitting ~95% of all these GBM spectra. The comparison with spectral fit results based on previous empirical models demonstrates that the past exclusion of synchrotron radiation as an emission mechanism derived via the line-of-death was misleading. Our analysis probes the physics of these ultra-relativistic outflows and the related microphysical processes, and for the first time provides estimates of magnetic field strength and Lorentz factors of the emitting region directly from spectral fits. Our modeling of the Fermi/GBM observations provides evidence that GRBs are produced by moderately magnetized jets in which relativistic mini-jets emit optically-thin synchrotron radiation at large emission radii., Comment: submitted

With the confirmed detection of short gamma-ray burst (GRB) in association with a gravitational wave signal, we present the first fully Bayesian {\it Fermi}-GBM short GRB spectral catalog. Both peak flux and time-resolved spectral results are presented. Additionally, we release the full posterior distributions and reduced data from our sample. Following our previous study, we introduce three variability classes based of the observed light curve structure., Comment: Submitted to ApJL

The unexpected nearby gamma-ray burst GRB 170817A associated with the LIGO binary neutron star merger event GW170817 presents a challenge to the current understanding of the emission physics of short gamma-ray bursts (GRBs). The event's low luminosity but similar peak energy compared to standard short GRBs are difficult to explain with current models, challenging our understanding of the GRB emission process. Emission models invoking synchrotron radiation from electrons accelerated in shocks and photospheric emission are particularly challenging explanations for this burst., Comment: Submitted to ApJL on the 22nd of October

The recent coincident detection of gravitational waves (GW) from a binary neutron star merger with aLIGO/Virgo and short-lived gamma-ray emission with Fermi/GBM (called GW 170817) is a milestone for the establishment of multi-messenger astronomy. Merging neutron stars (NS) represent the standard scenario for short-duration (< 2 sec) gamma-ray bursts (GRBs) which are produced in a collimated, relativistically expanding jet with an opening angle of a few degrees and a bulk Lorentz factor of 300-1000. While the present aLIGO detection is consistent with predictions, the measured faint gamma-ray emission from GW 170817A, if associated to the merger event at a distance of 40 Mpc, is about 1000x less luminous than known short-duration GRBs (sGRBs). Hence, the presence of this sGRB in the local Universe is either a very rare event, or points to a dramatic ignorance of the emission properties of sGRBs outside their narrow jets. Here we show that the majority of previously detected faint sGRBs are local, at redshift smaller than 0.1, seen off-axis. In contrast, the brighter sGRBs are seen on-axis, and therefore out to larger distances, consistent with the measured redshift distribution. Examining the observer-frame parameter space of all Fermi/GBM sGRBs shows that the sGRB associated with GW 170817A is extreme in its combination of flux, spectral softness and temporal structure. We identify a group of similar GRBs, one of which has been associated to a bright galaxy at 75 Mpc. We incorporate off-axis emission in the estimate of the rates of sGRBs, and predict that the majority of future GW-detections of NS-NS mergers will be accompanied by faint gamma-ray emission, contrary to previous thinking. The much more frequent off-axis emission of sGRBs also implies a much higher deadly rate of gamma-rays for extraterrestrial life in the Universe., Comment: originally submitted September 13, 2017 and forced to withdraw, now re-submitted in its original form

The spectral width and sharpness of unfolded, observed GRB spectra have been presented as a new tool to infer physical properties about GRB emission via spectral fitting of empirical models. Following the tradition of the 'line-of-death', the spectral width has been used to rule out synchrotron emission in a majority of GRBs. This claim is investigated via reexamination of previously reported width measures. Then, a sample of peak-flux GRB spectra are fit with an idealized, physical synchrotron model. It is found that many spectra can be adequately fit by this model even when the width measures would reject it. Thus, the results advocate for fitting a physical model to be the sole tool for testing that model. Finally, a smoothly-broken power law is fit to these spectra allowing for the spectral curvature to vary during the fitting process in order to understand why the previous width measures poorly predict the spectra. It is found that the failing of previous width measures is due to a combination of inferring physical parameters from unfolded spectra as well as the presence of multiple widths in the data beyond what the Band function can model., Comment: Accepted in A&A

The accurate spatial location of gamma-ray bursts (GRBs) is crucial for both producing a detector response matrix (DRM) and follow-up observations by other instruments. The Fermi Gamma-ray Burst Monitor (GBM) has the largest field of view (FOV) for detecting GRBs as it views the entire unocculted sky, but as a non-imaging instrument it relies on the relative count rates observed in each of its 14 detectors to localize transients. Improving its ability to accurately locate GRBs and other transients is vital to the paradigm of multi-messenger astronomy, including the electromagnetic follow-up of gravitational wave signals. Here we present the BAyesian Location Reconstruction Of GRBs ({\tt BALROG}) method for localizing and characterising GBM transients. Our approach eliminates the systematics of previous approaches by simultaneously fitting for the location and spectrum of a source. It also correctly incorporates the uncertainties in the location of a transient into the spectral parameters and produces reliable positional uncertainties for both well-localized sources and those for which the GBM data cannot effectively constrain the position. While computationally expensive, {\tt BALROG} can be implemented to enable quick follow-up of all GBM transient signals. Also, we identify possible response problems that require attention as well as caution when using standard, public GBM DRMs. Finally, we examine the effects of including the variance in location on the spectral parameters of GRB 080916C. We find that spectral parameters change and no extra components are required when these effects are included in contrast to when we use a fixed location. This finding has the potential to alter both the GRB spectral catalogs as well as the reported spectral composition of some well-known GRBs., Comment: 22 pages, 23 figures, submitted to MNRAS

Since its launch in 2008, the Fermi Gamma-ray Burst Monitor (GBM) has triggered and located on average approximately two gamma-ray bursts (GRB) every three days. Here we present the third of a series of catalogs of GRBs detected by GBM, extending the second catalog by two more years, through the middle of July 2014. The resulting list includes 1405 triggers identified as GRBs. The intention of the GBM GRB catalog is to provide information to the community on the most important observables of the GBM detected GRBs. For each GRB the location and main characteristics of the prompt emission, the duration, peak flux and fluence are derived. The latter two quantities are calculated for the 50-300~keV energy band, where the maximum energy release of GRBs in the instrument reference system is observed, and also for a broader energy band from 10-1000 keV, exploiting the full energy range of GBM's low-energy NaI(Tl) detectors. Using statistical methods to assess clustering, we find that the hardness and duration of GRBs are better fitted by a two-component model with short-hard and long-soft bursts, than by a model with three components. Furthermore, information is provided on the settings and modifications of the triggering criteria and exceptional operational conditions during years five and six in the mission. This third catalog is an official product of the Fermi GBM science team, and the data files containing the complete results are available from the High-Energy Astrophysics Science Archive Research Center (HEASARC)., Comment: 225 pages, 13 figures and 8 tables. Accepted for publication in Astrophysical Journal Supplement 2016

Ultra-relativistic motion and efficient conversion of kinetic energy to radiation are required by gamma-ray burst (GRB) observations, yet they are difficult to simultaneously achieve. Three leading mechanisms have been proposed to explain the observed emission emanating from GRB outflows: radiation from either relativistic internal or external shocks, or thermal emission from a photosphere. Previous works were dedicated to independently treating these three mechanisms and arguing for a sole, unique origin of the prompt emission of gamma-ray bursts. In contrast, herein, we first explain why all three models are valid mechanisms and that a contribution from each of them is expected in the prompt phase. Additionally, we show that a single parameter, the dimensionless entropy of the GRB outflow, determines which mechanism contributes the most to the emission. More specifically, internal shocks dominate for low values of the dimensionless entropy, external shocks for intermediate values and finally, photospheric emission for large values. We present a unified framework for the emission mechanisms of GRBs with easily testable predictions for each process., Comment: 6 pages, 1 fig, accepted by ApJ

Gamma-ray bursts (GRBs) are characterised by a strong correlation between the instantaneous luminosity and the spectral peak energy within a burst. This correlation, which is known as the hardness-intensity correlation or the Golenetskii correlation, not only holds important clues to the physics of GRBs but is thought to have the potential to determine redshifts of bursts. In this paper, I use a hierarchical Bayesian model to study the universality of the rest-frame Golenetskii correlation and in particular I assess its use as a redshift estimator for GRBs. I find that, using a power-law prescription of the correlation, the power-law indices cluster near a common value, but have a broader variance than previously reported ($\sim 1-2$). Furthermore, I find evidence that there is spread in intrinsic rest-frame correlation normalizations for the GRBs in our sample ($\sim 10^{51}-10^{53}$ erg s$^{-1}$). This points towards variable physical settings of the emission (magnetic field strength, number of emitting electrons, photospheric radius, viewing angle, etc.). Subsequently, these results eliminate the Golenetskii correlation as a useful tool for redshift determination and hence a cosmological probe. Nevertheless, the Bayesian method introduced in this paper allows for a better determination of the rest frame properties of the correlation, which in turn allows for more stringent limitations for physical models of the emission to be set., Comment: submitted to MNRAS

The prompt emission of the long, smooth, and single-pulsed gamma-ray burst, GRB $\textit{141028A}$, is analyzed under the guise of an external shock model. First, we fit the $\gamma$-ray spectrum with a two-component photon model, namely synchrotron+blackbody, and then fit the recovered evolution of the synchrotron $\nu F_{\nu}$ peak to an analytic model derived considering the emission of a relativistic blast-wave expanding into an external medium. The prediction of the model for the $\nu F_{\nu}$ peak evolution matches well with the observations. We observe the blast-wave transitioning into the deceleration phase. Further we assume the expansion of the blast-wave to be nearly adiabatic, motivated by the low magnetic field deduced from the observations. This allows us to recover within an order of magnitude the flux density at the $\nu F_{\nu}$ peak, which is remarkable considering the simplicity of the analytic model. Across all wavelengths, synchrotron emission from a single forward shock provides a sufficient solution for the observations. Under this scenario we argue that the distinction between $\textit{prompt}$ and $\textit{ afterglow}$ emission is superfluous as both early and late time emission emanate from the same source. While the external shock model is clearly not a universal solution, this analysis opens the possibility that at least some fraction of GRBs can be explained with an external shock origin of their prompt phase.

The analysis of gamma-ray burst (GRB) spectra with multi-component emission models has become an important part of the field. In particular, multi-component analysis where one component is a blackbody representing emission from a photosphere has enabled both a more detailed understanding of the energy content of the jet as well as the ability to examine the dynamic structure of the outflow. While the existence of a blackbody-like component has been shown to be significant and not a byproduct of background fluctuations, it is very possible that it can be an artifact of spectral evolution of a single component that is being poorly resolved in time. Herein, this possibility is tested by simulating a single component evolving in time and then folding the spectra through the $Fermi$ detector response to generate time-tagged event Gamma-ray Burst Monitor (GBM) data. We then fit both the time integrated and resolved generated spectral data with a multi-component model using standard tools. It is found that in {\it time-integrated} spectra, a blackbody can be falsely identified due to the spectral curvature introduced by the spectral evolution. However, in time-resolved analysis defined by time bins that can resolve the evolution of the spectra, the significance of the falsely identified blackbody is very low. Additionally, the evolution of the artificial blackbody parameters does not match the recurring behavior that has been identified in the actual observations. These results reinforce the existence of the blackbody found in {\it time-resolved} analysis of GRBs and stress the point that caution should be taken when using time-integrated spectral analysis for identifying physical properties of GRBs.

The long standing problem of identifying the emission mechanism operating in gamma-ray bursts (GRBs) has produced a myriad of possible models that have the potential of explaining the observations. Generally, the empirical Band function is fit to the observed gamma-ray data and the fit parameters are used to infer which radiative mechanisms are at work in GRB outflows. In particular, the distribution of the Band function's low-energy power law index, $\alpha$, has led to the so-called synchrotron "line-of-death" (LOD) which is a statement that the distribution cannot be explained by the simplest of synchrotron models alone. As an alternatively fitting model, a combination of a blackbody in addition to the Band function is used, which in many cases provide a better or equally good fit. It has been suggested that such fits would be able to alleviate the LOD problem for synchrotron emission in GRBs. However, these conclusions rely on the Band function's ability to fit a synchrotron spectrum within the observed energy band. In order to investigate if this is the case, we simulate synchrotron and synchrotron+blackbody spectra and fold them through the instrumental response of the $Fermi$ Gamma-ray Burst Monitor (GBM). We then perform a standard data analysis by fitting the simulated data with both Band and Band+blackbody models. We find two important results: the synchrotron LOD is actually more severe than the original predictions: $\alpha_{\bf LOD} \sim -0.8$. Moreover, we find that intrinsic synchrotron+blackbody emission is insufficient to account for the entire observed $\alpha$ distribution. This implies that some other emission mechanism(s) are required to explain a large fraction of observed GRBs.

We study the time-resolved spectra of eight GRBs observed by Fermi GBM in its first five years of mission, with 1 keV - 1 MeV fluence $f>1.0\times10^{-4}$ erg cm$^{-2}$ and signal-to-noise level $\text{S/N}\geq10.0$ above 900 keV. We aim to constrain in detail the spectral properties of GRB prompt emission on a time-resolved basis and to discuss the theoretical implications of the fitting results in the context of various prompt emission models. We perform time-resolved spectral analysis using a variable temporal binning technique according to optimal S/N criteria, resulting in a total of 299 time-resolved spectra. We fit the Band function to all spectra and obtain the distributions for the low-energy power-law index $\alpha$, the high-energy power-law index $\beta$, the peak energy in the observed $\nu F_\nu$ spectrum $E_\text{p}$, and the difference between the low- and high-energy power-law indices $\Delta s=\alpha-\beta$. Using the distributions of $\Delta s$ and $\beta$, the electron population index $p$ is found to be consistent with the "moderately fast" scenario which fast- and slow-cooling scenarios cannot be distinguished. We also apply a physically motivated synchrotron model, which is a triple power-law with constrained power-law indices and a blackbody component, to test for consistency with a synchrotron origin for the prompt emission and obtain the distributions for the two break energies $E_\text{b,1}$ and $E_\text{b,2}$, the middle segment power-law index $\beta$, and the Planck function temperature $kT$. A synchrotron model is found consistent with the majority of time-resolved spectra for these eight energetic Fermi GBM bursts with good high-energy photon statistics, as long as both the cooling and injection break are included and the leftmost spectral slope is lifted either by inclusion of a thermal component or when an evolving magnetic field is accounted for., Comment: 20 pages, 7 figures, 8 tables, accepted for publication in A&A

The understanding of the prompt $\gamma$-ray spectra of gamma-ray bursts (GRBs) is of great importance to correctly interpreting the physical mechanisms that produce the underlying event as well as the structure of the relativistic jet from which the emission emanates. Time-resolved analysis of these spectra is the main method of extracting information from the data. In this work, several techniques for temporal binning of GRB spectra are examined to understand the systematics associated with each with the goal of finding the best method(s) to bin lightcurves for analysis. The following binning methods are examined: constant cadence (CC), Bayesian blocks (BB), signal-to-noise (S/N) and Knuth bins (KB). I find that both the KB and BB methods reconstruct the intrinsic spectral evolution accurately while the S/N method fails in most cases. The CC method is accurate when the cadence is not too coarse but does not necessarily bin the data based on the true source variability. Additionally, the integrated pulse properties are investigated and compared to the time-resolved properties. If intrinsic spectral evolution is present then the integrated properties are not useful in identifying physical and cosmological properties of GRBs without knowing the physical emission mechanism and its evolution.

Recent observations by the $Fermi$ Gamma-ray Space Telescope have confirmed the existence of thermal and non-thermal components in the prompt photon spectra of some Gamma-ray bursts (GRBs). Through an analysis of six bright Fermi GRBs, we have discovered a correlation between the observed photospheric and non-thermal $\gamma$-ray emission components of several GRBs using a physical model that has previously been shown to be a good fit to the Fermi data. From the spectral parameters of these fits we find that the characteristic energies, $E_{\rm p}$ and $kT$, of these two components are correlated via the relation $E_{\rm p} \propto T^{\alpha}$ which varies from GRB to GRB. We present an interpretation in which the value of index $\alpha$ indicates whether the jet is dominated by kinetic or magnetic energy. To date, this jet composition parameter has been assumed in the modeling of GRB outflows rather than derived from the data.

This is the second of a series of catalogs of gamma-ray bursts (GRBs) observed with the Fermi Gamma-ray Burst Monitor (GBM). It extends the first two-year catalog by two more years, resulting in an overall list of 953 GBM triggered GRBs. The intention of the GBM GRB catalog is to provide information to the community on the most important observables of the GBM detected GRBs. For each GRB the location and main characteristics of the prompt emission, the duration, peak flux and fluence are derived. The latter two quantities are calculated for the 50 - 300 keV energy band, where the maximum energy release of GRBs in the instrument reference system is observed and also for a broader energy band from 10 - 1000 keV, exploiting the full energy range of GBMs low-energy detectors. Furthermore, information is given on the settings and modifications of the triggering criteria and exceptional operational conditions during years three and four in the mission. This second catalog is an official product of the Fermi GBM science team, and the data files containing the complete results are available from the High-Energy Astrophysics Science Archive Research Center (HEASARC)., Comment: 151 pages, 10 figures, 9 tables; to be published in ApJS

Gamma-ray burst (GRB) 130427A is one of the most energetic GRBs ever observed. The initial pulse up to 2.5 s is possibly the brightest well-isolated pulse observed to date. A fine time resolution spectral analysis shows power-law decays of the peak energy from the onset of the pulse, consistent with models of internal synchrotron shock pulses. However, a strongly correlated power-law behavior is observed between the luminosity and the spectral peak energy that is inconsistent with curvature effects arising in the relativistic outflow. It is difficult for any of the existing models to account for all of the observed spectral and temporal behaviors simultaneously., Comment: This is the author's version of the work. It is posted here by permission of the AAAS for personal use, not for redistribution. The definitive version was published online in Science Express on 21 November 2013 [DOI:10.1126/science.1242302] 15 pages, 4 figures and 1 table; includes supplementary online materials

A Band function has become the standard spectral function used to describe the prompt emission spectra of gamma-ray bursts (GRBs). However, deviations from this function have previously been observed in GRBs detected by BATSE and in individual GRBs from the \textit{Fermi} era. We present a systematic and rigorous search for spectral deviations from a Band function at low energies in a sample of the first two years of high fluence, long bursts detected by the \textit{Fermi} Gamma-Ray Burst Monitor (GBM). The sample contains 45 bursts with a fluence greater than 2$\times10^{-5}$ erg / cm$^{2}$ (10 - 1000 keV). An extrapolated fit method is used to search for low-energy spectral anomalies, whereby a Band function is fit above a variable low-energy threshold and then the best fit function is extrapolated to lower energy data. Deviations are quantified by examining residuals derived from the extrapolated function and the data and their significance is determined via comprehensive simulations which account for the instrument response. This method was employed for both time-integrated burst spectra and time-resolved bins defined by a signal to noise ratio of 25 $\sigma$ and 50 $\sigma$. Significant deviations are evident in 3 bursts (GRB\,081215A, GRB\,090424 and GRB\,090902B) in the time-integrated sample ($\sim$ 7%) and 5 bursts (GRB\,090323, GRB\,090424, GRB\,090820, GRB\,090902B and GRB\,090926A) in the time-resolved sample ($\sim$ 11%).} The advantage of the systematic, blind search analysis is that it can demonstrate the requirement for an additional spectral component without any prior knowledge of the nature of that extra component. Deviations are found in a large fraction of high fluence GRBs; fainter GRBs may not have sufficient statistics for deviations to be found using this method.

The Fermi Gamma-ray Burst Monitor (GBM) is designed to enhance the scientific return from Fermi in studying gamma-ray bursts (GRBs). In its first two years of operation GBM triggered on 491 GRBs. We summarize the criteria used for triggering and quantify the general characteristics of the triggered GRBs, including their locations, durations, peak flux, and fluence. This catalog is an official product of the Fermi GBM science team, and the data files containing the complete results are available from the High-Energy Astrophysics Science Archive Research Center (HEASARC)., Comment: 92 pages, 11 figures, 5 tables; to be published in ApJS

We present systematic spectral analyses of GRBs detected by the Fermi Gamma-Ray Burst Monitor (GBM) during its first two years of operation. This catalog contains two types of spectra extracted from 487 GRBs, and by fitting four different spectral models, this results in a compendium of over 3800 spectra. The models were selected based on their empirical importance to the spectral shape of many GRBs, and the analysis performed was devised to be as thorough and objective as possible. We describe in detail our procedure and criteria for the analyses, and present the bulk results in the form of parameter distributions. This catalog should be considered an official product from the Fermi GBM Science Team, and the data files containing the complete results are available from the High-Energy Astrophysics Science Archive Research Center (HEASARC)., Comment: 43 pages, 24 Figures, accepted to ApJ Supplements; http://heasarc.gsfc.nasa.gov/W3Browse/fermi/fermigbrst.html

Swift/BAT detected the first burst from 1E 1841-045 in May 2010 with intermittent burst activity recorded through at least July 2011. Here we present Swift and Fermi/GBM observations of this burst activity and search for correlated changes to the persistent X-ray emission of the source. The T90 durations of the bursts range between 18-140 ms, comparable to other magnetar burst durations, while the energy released in each burst ranges between (0.8 - 25)E38 erg, which is in the low side of SGR bursts. We find that the bursting activity did not have a significant effect on the persistent flux level of the source. We argue that the mechanism leading to this sporadic burst activity in 1E 1841-045 might not involve large scale restructuring (either crustal or magnetospheric) as seen in other magnetar sources., Comment: accepted for publication in The Astrophysical Journal Letters

Discerning the radiative dissipation mechanism for prompt emission in Gamma-Ray Bursts (GRBs) requires detailed spectroscopic modeling that straddles the $\nu F_{\nu}$ peak in the 100 keV - 1 MeV range. Historically, empirical fits such as the popular Band function have been employed with considerable success in interpreting the observations. While extrapolations of the Band parameters can provide some physical insight into the emission mechanisms responsible for GRBs, these inferences do not provide a unique way of discerning between models. By fitting physical models directly this degeneracy can be broken, eliminating the need for empirical functions; our analysis here offers a first step in this direction. One of the oldest, and leading, theoretical ideas for the production of the prompt signal is the synchrotron shock model (SSM). Here we explore the applicability of this model to a bright {\it Fermi} GBM burst with a simple temporal structure, GRB {\it 090820}A. Our investigation implements, for the first time, thermal and non-thermal synchrotron emissivities in the RMFIT forward-folding spectral analysis software often used in GBM burst studies. We find that these synchrotron emissivities, together with a blackbody shape, provide at least as good a match with the data as the Band GRB spectral fitting function. This success is achieved in both time-integrated and time-resolved spectral fits.

We present here our results of the temporal and spectral analysis of a sample of 52 bright and hard gamma-ray bursts (GRBs) observed with the Fermi Gamma-ray Burst Monitor (GBM) during its first year of operation (July 2008-July 2009). Our sample was selected from a total of 253 GBM GRBs based on each event peak count rate measured between 0.2 and 40MeV. The final sample comprised 34 long and 18 short GRBs. These numbers show that the GBM sample contains a much larger fraction of short GRBs, than the CGRO/BATSE data set, which we explain as the result of our (different) selection criteria and the improved GBM trigger algorithms, which favor collection of short, bright GRBs over BATSE. A first by-product of our selection methodology is the determination of a detection threshold from the GBM data alone, above which GRBs most likely will be detected in the MeV/GeV range with the Large Area Telescope (LAT) onboard Fermi. This predictor will be very useful for future multiwavelength GRB follow ups with ground and space based observatories. Further we have estimated the burst durations up to 10MeV and for the first time expanded the duration-energy relationship in the GRB light curves to high energies. We confirm that GRB durations decline with energy as a power law with index approximately -0.4, as was found earlier with the BATSE data and we also notice evidence of a possible cutoff or break at higher energies. Finally, we performed time-integrated spectral analysis of all 52 bursts and compared their spectral parameters with those obtained with the larger data sample of the BATSE data. We find that the two parameter data sets are similar and confirm that short GRBs are in general harder than longer ones., Comment: 40 pages, 11 figures, 3 tables, Submitted to ApJ

The jet opening angle of gamma-ray bursts (GRBs) is an important parameter for determining the characteristics of the progenitor, and the information contained in the opening angle gives insight into the relativistic outflow and the total energy that is contained in the burst. Unfortunately, a confident inference of the jet opening angle usually requires broadband measurement of the afterglow of the GRB, from the X-ray down to the radio and from minutes to days after the prompt gamma-ray emission, which may be difficult to obtain. For this reason, very few of all detected GRBs have constrained jet angles. We present an alternative approach to derive jet opening angles from the prompt emission of the GRB, given that the GRB has a measurable Epeak and fluence, and which does not require any afterglow measurements. We present the distribution of derived jet opening angles for the first two years of the Fermi Gamma-ray Burst Monitor (GBM) operation, and we compare a number of our derived opening angles to the reported opening angles using the traditional afterglow method. We derive the collimation-corrected gamma-ray energy, E_{\gamma}, for GRBs with redshift and find that some of the GRBs in our sample are inconsistent with a proto-magnetar progenitor. Finally, we show that the use of the derived jet opening angles results in a tighter correlation between the rest-frame Epeak and E_{\gamma} than has previously been presented, which places long GRBs and short GRBs onto one empirical power law., Comment: 23 pages, 6 figures, 1 table; submitted to ApJ

In three years of observations since the beginning of nominal science operations in 2008 August, the Large Area Telescope (LAT) on board the Fermi Gamma-Ray Space Telescope has observed high-energy great than (20 MeV) gamma-ray emission from 35 gamma-ray bursts (GRBs). Among these, 28 GRBs have been detected above 100 MeV and 7 GRBs above approximately 20 MeV. The first Fermi-LAT catalog of GRBs is a compilation of these detections and provides a systematic study of high-energy emission from GRBs for the first time. To generate the catalog, we examined 733 GRBs detected by the Gamma-Ray Burst Monitor (GBM) on Fermi and processed each of them using the same analysis sequence. Details of the methodology followed by the LAT collaboration for the GRB analysis are provided. We summarize the temporal and spectral properties of the LAT-detected GRBs. We also discuss characteristics of LAT-detected emission such as its delayed onset and longer duration compared with emission detected by the GBM, its power-law temporal decay at late times, and the fact that it is dominated by a power-law spectral component that appears in addition to the usual Band model.