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.

We present the results of our year-long afterglow monitoring of GW170817, the first binary neutron star (NS) merger detected by advanced LIGO and advanced Virgo. New observations with the Australian Telescope Compact Array (ATCA) and the Chandra X-ray Telescope were used to constrain its late-time behavior. The broadband emission, from radio to X-rays, is well-described by a simple power-law spectrum with index ~0.585 at all epochs. After an initial shallow rise ~t^0.9, the afterglow displayed a smooth turn-over, reaching a peak X-ray luminosity of ~5e39 erg/s at 160 d, and has now entered a phase of rapid decline ~t^(-2). The latest temporal trend challenges most models of choked jet/cocoon systems, and is instead consistent with the emergence of a relativistic structured jet seen at an angle of ~22 deg from its axis. Within such model, the properties of the explosion (such as its blastwave energy E_K~2E50 erg, jet width theta_c~4 deg, and ambient density n~3E-3 cm^(-3)) fit well within the range of properties of cosmological short GRBs.

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.

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.

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.

We report a framework for spectroscopic follow-up design for optimizing supernova photometric classification. The strategy accounts for the unavoidable mismatch between spectroscopic and photometric samples, and can be used even in the beginning of a new survey -- without any initial training set. The framework falls under the umbrella of active learning (AL), a class of algorithms that aims to minimize labelling costs by identifying a few, carefully chosen, objects which have high potential in improving the classifier predictions. As a proof of concept, we use the simulated data released after the Supernova Photometric Classification Challenge (SNPCC) and a random forest classifier. Our results show that, using only 12\% the number of training objects in the SNPCC spectroscopic sample, this approach is able to double purity results. Moreover, in order to take into account multiple spectroscopic observations in the same night, we propose a semi-supervised batch-mode AL algorithm which selects a set of $N=5$ most informative objects at each night. In comparison with the initial state using the traditional approach, our method achieves 2.3 times higher purity and comparable figure of merit results after only 180 days of observation, or 800 queries (73% of the SNPCC spectroscopic sample size). Such results were obtained using the same amount of spectroscopic time necessary to observe the original SNPCC spectroscopic sample, showing that this type of strategy is feasible with current available spectroscopic resources. The code used in this work is available in the COINtoolbox: https://github.com/COINtoolbox/ActSNClass .

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.

Here we built up a sample of 22 GRBs at redshifts $z > 2$ observed with X-shooter to determine the abundances of hydrogen, metals, dust, and molecular species. This allows us to study the metallicity and dust depletion effects in the neutral ISM at high redshift and to answer the question whether (and why) there might be a lack of H$_2$ in GRB-DLAs. We fit absorption lines and measure the column densities of different metal species as well as atomic and molecular hydrogen. The derived relative abundances are used to fit dust depletion sequences and determine the dust-to-metals ratio and the host-galaxy intrinsic visual extinction. There is no lack of H$_2$-bearing GRB-DLAs. We detect absorption lines from H$_2$ in 6 out of 22 GRB afterglow spectra, with molecular fractions ranging between $f\simeq 5\cdot10^{-5}$ and $f\simeq 0.04$, and claim tentative detections in three other cases. The GRB-DLAs in the present sample have on average low metallicities ($\mathrm{[X/H]}\approx -1.3$), comparable to the rare population of QSO-ESDLAs (log N(HI) $> 21.5$). H$_2$-bearing GRB-DLAs are found to be associated with significant dust extinction, $A_V > 0.1$ mag, and have dust-to-metals ratios DTM$ > 0.4$. All of these systems exhibit column densities of log N(HI) $> 21.7$. The overall fraction of H$_2$ detections is $\ge 27$% (41% including tentative detections), which is three times larger than in the general population of QSO-DLAs. For $2<z<4$, and for log N(HI) $> 21.7$, the H$_2$ detection fraction is 60-80% in GRB-DLAs as well as in extremely strong QSO-DLAs. This is likely a consequence of the fact that both GRB- and QSO-DLAs with high N(HI) probe sight-lines with small impact parameters that indicate that the absorbing gas is associated with the inner regions of the absorbing galaxy, where the gas pressure is higher and the conversion of HI to H$_2$ takes place.

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.

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.

A long-standing paradigm in astrophysics is that collisions- or mergers- of two neutron stars (NSs) form highly relativistic and collimated outflows (jets) powering gamma-ray bursts (GRBs) of short (< 2 s) duration. However, the observational support for this model is only indirect. A hitherto outstanding prediction is that gravitational wave (GW) events from such mergers should be associated with GRBs, and that a majority of these GRBs should be off-axis, that is, they should point away from the Earth. Here we report the discovery of the X-ray counterpart associated with the GW event GW170817. While the electromagnetic counterpart at optical and infrared frequencies is dominated by the radioactive glow from freshly synthesized r-process material in the merger ejecta, known as kilonova, observations at X-ray and, later, radio frequencies exhibit the behavior of a short GRB viewed off-axis. Our detection of X-ray emission at a location coincident with the kilonova transient provides the missing observational link between short GRBs and GWs from NS mergers, and gives independent confirmation of the collimated nature of the GRB emission.

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.

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.

In view of the recent report by Connaughton we analyse continuous TTE data of Fermi-GBM around the time of the gravitational wave event GW 150914. We find that after proper accounting for low count statistics, the GBM transient event at 0.4 s after GW 150914 is likely not due to an astrophysical source, but consistent with a background fluctuation, removing the tension between the INTEGRAL/ACS non-detection and GBM. Additionally, reanalysis of other short GRBs shows that without proper statistical modeling the fluence of faint events is over-predicted, as verified for some joint GBM-ACS detections of short GRBs. We detail the statistical procedure to correct these biases. As a result, faint short GRBs, verified by ACS detections, with significances in the broad-band light curve even smaller than that of the GBM-GW150914 event are recovered as proper non-zero source, while the GBM-GW150914 event is consistent with zero fluence.

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).

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.

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.

It has been suggested that the prompt emission in gamma-ray bursts consists of several components giving rise to the observed spectral shape. Here we examine a sample of the 8 brightest, single pulsed {\it Fermi} bursts whose spectra are modelled by using synchrotron emission as one of the components. Five of these bursts require an additional photospheric component (blackbody). In particular, we investigate the inferred properties of the jet and the physical requirements set by the observed components for these five bursts, in the context of a baryonic dominated outflow, motivated by the strong photospheric component. We find similar jet properties for all five bursts: the bulk Lorentz factor decreases monotonously over the pulses and lies between 1000 and 100. This evolution is robust and can neither be explained by a varying radiative efficiency nor a varying magnetisation of the jet assuming the photosphere radius is above the coasting radius. Such a behaviour challenges several dissipation mechanisms, e.g., the internal shocks. Furthermore, in all 8 cases the data clearly reject a fast-cooled synchrotron spectrum (in which a significant fraction of the emitting electrons have cooled to energies below the minimum injection energy), inferring a typical electron Lorentz factor of $10^4 - 10^7$. Such values are much higher than what is typically expected in internal shocks. Therefore, while the synchrotron scenario is not rejected by the data, the interpretation does present several limitations that need to be addressed. Finally, we point out and discuss alternative interpretations.

Astrophysical sources are now observed by many different instruments at different wavelengths, from radio to high-energy gamma-rays, with an unprecedented quality. Putting all these data together to form a coherent view, however, is a very difficult task. Each instrument has its own data format, software and analysis procedure, which are difficult to combine. It is for example very challenging to perform a broadband fit of the energy spectrum of the source. The Multi-Mission Maximum Likelihood framework (3ML) aims to solve this issue, providing a common framework which allows for a coherent modeling of sources using all the available data, independent of their origin. At the same time, thanks to its architecture based on plug-ins, 3ML uses the existing official software of each instrument for the corresponding data in a way which is transparent to the user. 3ML is based on the likelihood formalism, in which a model summarizing our knowledge about a particular region of the sky is convolved with the instrument response and compared to the corresponding data. The user can choose between a frequentist analysis, and a Bayesian analysis. In the former, parameters of the model are optimized in order to obtain the best match to the data (i.e., the maximum of the likelihood). In the latter, the priors specified by the user are used to build the posterior distribution, which is then sampled with Markov Chain Monte Carlo or Multinest. Our implementation of this idea is very flexible, allowing the study of point sources as well as extended sources with arbitrary spectra. We will review the problem we aim to solve, the 3ML concepts and its innovative potential.

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.

Much evidence points towards that the photosphere in the relativistic outflow in GRBs plays an important role in shaping the observed MeV spectrum. However, it is unclear whether the spectrum is fully produced by the photosphere or whether a substantial part of the spectrum is added by processes far above the photosphere. Here we make a detailed study of the $\gamma-$ray emission from single pulse GRB110920A which has a spectrum that becomes extremely narrow towards the end of the burst. We show that the emission can be interpreted as Comptonisation of thermal photons by cold electrons in an unmagnetised outflow at an optical depth of $\tau \sim 20$. The electrons receive their energy by a local dissipation occurring close to the saturation radius. The main spectral component of GRB110920A and its evolution is thus, in this interpretation, fully explained by the emission from the photosphere including localised dissipation at high optical depths.

We present observations by the Fermi Gamma-Ray Space Telescope Gamma-Ray Burst Monitor (GBM) of the nearby (z=0.55) GRB 101219B. This burst is a long GRB, with an associated supernova and with a blackbody component detected in the early afterglow observed by the Swift X-ray Telescope (XRT). Here we show that the prompt gamma-ray emission has a blackbody spectrum, making this the second such burst observed by Fermi GBM. The properties of the blackbody, together with the redshift and our estimate of the radiative efficiency, makes it possible to calculate the absolute values of the properties of the outflow. We obtain an initial Lorentz factor Gamma=138\pm 8, a photospheric radius r_phot=4.4\pm 1.9 \times 10^{11} cm and a launch radius r_0=2.7\pm 1.6 \times 10^{7} cm. The latter value is close to the black hole and suggests that the jet has a relatively unobstructed path through the star. There is no smooth connection between the blackbody components seen by GBM and XRT, ruling out the scenario that the late emission is due to high-latitude effects. In the interpretation that the XRT blackbody is prompt emission due to late central engine activity, the jet either has to be very wide or have a clumpy structure where the emission originates from a small patch. Other explanations for this component, such as emission from a cocoon surrounding the jet, are also possible.

The Fermi Gamma-ray Burst Monitor (GBM) has detected over 1400 Gamma-Ray Bursts (GRBs) since it began science operations in July, 2008. We use a subset of over 300 GRBs localized by instruments such as Swift, the Fermi Large Area Telescope, INTEGRAL, and MAXI, or through triangulations from the InterPlanetary Network (IPN), to analyze the accuracy of GBM GRB localizations. We find that the reported statistical uncertainties on GBM localizations, which can be as small as 1 degree, underestimate the distance of the GBM positions to the true GRB locations and we attribute this to systematic uncertainties. The distribution of systematic uncertainties is well represented (68% confidence level) by a 3.7 degree Gaussian with a non-Gaussian tail that contains about 10% of GBM-detected GRBs and extends to approximately 14 degrees. A more complex model suggests that there is a dependence of the systematic uncertainty on the position of the GRB in spacecraft coordinates, with GRBs in the quadrants on the Y-axis better localized than those on the X-axis.

Simultaneous Swift and Fermi observations of gamma-ray bursts (GRBs) offer a unique broadband view of their afterglow emission, spanning more than ten decades in energy. We present the sample of X-ray flares observed by both Swift and Fermi during the first three years of Fermi operations. While bright in the X-ray band, X-ray flares are often undetected at lower (optical), and higher (MeV to GeV) energies. We show that this disfavors synchrotron self-Compton processes as origin of the observed X-ray emission. We compare the broadband properties of X-ray flares with the standard late internal shock model, and find that, in this scenario, X-ray flares can be produced by a late-time relativistic (Gamma>50) outflow at radii R~10^13-10^14 cm. This conclusion holds only if the variability timescale is significantly shorter than the observed flare duration, and implies that X-ray flares can directly probe the activity of the GRB central engine.

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.

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.

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 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).

Time-resolved spectroscopy is performed on eight bright, long gamma-ray bursts (GRBs) dominated by single emission pulses that were observed with the {\it Fermi Gamma-ray Space Telescope}. Fitting the prompt radiation of GRBs by empirical spectral forms such as the Band function leads to ambiguous conclusions about the physical model for the prompt radiation. Moreover, the Band function is often inadequate to fit the data. The GRB spectrum is therefore modeled with two emission components consisting of optically thin nonthermal synchrotron radiation from relativistic electrons and, when significant, thermal emission from a jet photosphere, which is represented by a blackbody spectrum. To produce an acceptable fit, the addition of a blackbody component is required in 5 out of the 8 cases. We also find that the low-energy spectral index \alpha is consistent with a synchrotron component with \alpha = -0.81\pm 0.1. This value lies between the limiting values of \alpha = -2/3 and \alpha = -3/2 for electrons in the slow and fast-cooling regimes, respectively, suggesting ongoing acceleration at the emission site. The blackbody component can be more significant when using a physical synchrotron model instead of the Band function, illustrating that the Band function does not serve as a good proxy for a nonthermal synchrotron emission component. The temperature and characteristic emission-region size of the blackbody component are found to, respectively, decrease and increase as power laws with time during the prompt phase. In addition, we find that the blackbody and nonthermal components have separate temporal behaviors.

We present systematic spectral analyses of GRBs detected with the Burst and Transient Source Experiment (BATSE) onboard the Compton Gamma-Ray Observatory (CGRO) during its entire nine years of operation. This catalog contains two types of spectra extracted from 2145 GRBs and fitted with five different spectral models resulting in a compendium of over 19000 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 procedures 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 BATSE Science Team, and the data files containing the complete results are available from the High-Energy Astrophysics Science Archive Research Center (HEASARC).

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.

{\it Fermi Gamma-ray Space Telescope} observations of GRB110721A have revealed two emission components from the relativistic jet: emission from the photosphere, peaking at $\sim 100$ keV and a non-thermal component, which peaks at $\sim 1000$ keV. We use the photospheric component to calculate the properties of the relativistic outflow. We find a strong evolution in the flow properties: the Lorentz factor decreases with time during the bursts from $\Gamma \sim 1000$ to $\sim 150$ (assuming a redshift $z=2$; the values are only weakly dependent on unknown efficiency parameters). Such a decrease is contrary to the expectations from the internal shocks and the isolated magnetar birth models. Moreover, the position of the flow nozzle measured from the central engine, $r_0$, increases by more than two orders of magnitude. Assuming a moderately magnetised outflow we estimate that $r_0$ varies from $10^6$ cm to $\sim 10^9$ cm during the burst. We suggest that the maximal value reflects the size of the progenitor core. Finally, we show that these jet properties naturally explain the observed broken power-law decay of the temperature which has been reported as a characteristic for GRB pulses.

The short GRB 120323A had the highest flux ever detected with the Fermi/GBM. Here we study its remarkable spectral properties and their evolution using two spectral models: (i) a single emission component scenario, where the spectrum is modeled by the empirical Band function, and (ii) a two component scenario, where thermal (Planck-like) emission is observed simultaneously with a non-thermal component (a Band function). We find that the latter model fits the integrated burst spectrum significantly better than the former, and that their respective spectral parameters are dramatically different: when fit with a Band function only, the Epeak of the event is unusually soft for a short GRB, while adding a thermal component leads to more typical short GRB values. Our time-resolved spectral analysis produces similar results. We argue here that the two-component model is the preferred interpretation for GRB 120323A, based on: (i) the values and evolution of the Band function parameters of the two component scenario, which are more typical for a short GRB, and (ii) the appearance in the data of a significant hardness-intensity correlation, commonly found in GRBs, when we employee two-component model fits; the correlation is non-existent in the Band-only fits. GRB 110721A, a long burst with an intense photospheric emission, exhibits the exact same behavior. We conclude that GRB 120323A has a strong photospheric emission contribution, first time observed in a short GRB. Magnetic dissipation models are difficult to reconcile with these results, which instead favor photospheric thermal emission and fast cooling synchrotron radiation from internal shocks. Finally, we derive a possibly universal hardness-luminosity relation in the source frame using a larger set of GRBs L,i=(1.59+/-0.84).10^50 (Epeak,i)^(1.33+/-0.07) erg/s), which could be used as a possible redshift estimator for cosmology.

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).

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).

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.

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.

Aims: In this paper we study the main spectral and temporal properties of gamma-ray bursts (GRBs) observed by Fermi/GBM. We investigate these key properties of GRBs in the rest-frame of the progenitor and test for possible intra-parameter correlations to better understand the intrinsic nature of these events. Methods: Our sample comprises 32 GRBs with measured redshift that were observed by GBM until August 2010. 28 of them belong to the long-duration population and 4 events were classified as short/hard bursts. For all of these events we derive, where possible, the intrinsic peak energy in the $\nu F_{\nu}$ spectrum (\eprest), the duration in the rest-frame, defined as the time in which 90% of the burst fluence was observed (\tninetyrest) and the isotropic equivalent bolometric energy (\eiso). Results: The distribution of \eprest has mean and median values of 1.1 MeV and 750 keV, respectively. A log-normal fit to the sample of long bursts peaks at ~800 keV. No high-\ep population is found but the distribution is biased against low \ep values. We find the lowest possible \ep that GBM can recover to be ~ 15 keV. The \tninetyrest distribution of long GRBs peaks at ~10 s. The distribution of \eiso has mean and median values of $8.9\times 10^{52}$ erg and $8.2 \times 10^{52}$ erg, respectively. We confirm the tight correlation between \eprest and \eiso (Amati relation) and the one between \eprest and the 1-s peak luminosity ($L_p$) (Yonetoku relation). Additionally, we observe a parameter reconstruction effect, i.e. the low-energy power law index $\alpha$ gets softer when \ep is located at the lower end of the detector energy range. Moreover, we do not find any significant cosmic evolution of neither \eprest nor \tninetyrest.

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.

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.

In this paper we examine gamma-ray and optical data of GRB 091024, a gamma-ray burst (GRB) with an extremely long duration of T90~1020 s, as observed with the Fermi Gamma-Ray Burst Monitor (GBM). We present spectral analysis of all three distinct emission episodes using data from Fermi/GBM. Because of the long nature of this event, many ground-based optical telescopes slewed to its location within a few minutes and thus were able to observe the GRB during its active period. We compare the optical and gamma-ray light curves. Furthermore, we estimate a lower limit on the bulk Lorentz factor from the variability and spectrum of the GBM light curve and compare it with that obtained from the peak time of the forward shock of the optical afterglow. From the spectral analysis we note that, despite its unusually long duration, this burst is similar to other long GRBs, i.e. there is spectral evolution (both the peak energy and the spectral index vary with time) and spectral lags are measured. We find that the optical light curve is highly anti-correlated to the prompt gamma-ray emission, with the optical emission reaching the maximum during an epoch of quiescence in the prompt emission. We interpret this behavior as the reverse shock (optical flash), expected in the internal-external shock model of GRB emission but observed only in a handful of GRBs so far. The lower limit on the initial Lorentz factor deduced from the variability time scale ($\Gamma_{min}=195_{-110}^+{90}$)is consistent within the error to the one obtained using the peak time of the forward shock ($\Gamma_0=120$) and is also consistent with Lorentz factors of other long GRBs.