There is no complete description of the emission physics during the prompt phase in gamma-ray bursts. Spectral analyses, however, indicate that many spectra are narrower than what is expected for non-thermal emission models. Here, we reanalyse the sample of 37 bursts in \citet{Yu2019}, by fitting the narrowest time-resolved spectrum in each burst. We perform model comparison between a photospheric and a synchrotron emission model based on Bayesian evidence. We choose to compare the shape of the narrowest expected spectra: emission from the photosphere in a non-dissipative flow and slow-cooled synchrotron emission from a narrow electron distribution. We find that the photospheric spectral shape is preferred by $54 \pm 8 \%$ of the spectra (20/37), while $38 \pm 8 \%$ of the spectra (14/37) prefer the synchrotron spectral shape; three spectra are inconclusive. We hence conclude that GRB spectra are indeed very narrow and that more than half of the bursts have a photospheric emission episode. We also find that a third of all analysed spectra, not only prefer, but are also compatible with a non-dissipative photosphere, confirming previous similar findings. Furthermore, we notice that the spectra, that prefer the photospheric model, all have a low-energy power-law indices $\alpha > -0.5$. This means that $\alpha$ is a good estimator of which model is preferred by the data. Finally, we argue that the spectra which statistically prefer the synchrotron model, could equally well be caused by subphotospheric dissipation. If that is the case, photospheric emission during the early, prompt phase would be even more dominant.

A large fraction of gamma-ray burst (GRB) spectra are very hard below the peak. Indeed, the observed distribution of sub-peak power-law indices, $\alpha$, has been used as an argument for a photospheric origin of GRB spectra. Here, we investigate what fraction of GRBs have spectra that are consistent with emission from a photopshere in a non-dissipative outflow. This is the simplest possible photospheric emission scenario. We create synthetic spectra, with a range of peak energies, by folding the theoretical predictions through the detector response of the FERMI/GBM detector. These simulated spectral data are fit with typically employed empirical models. We find that the low-energy photon indices obtain values ranging $-0.4 < \alpha < 0.0$, peaking at around $-0.1$, thus covering a non-negligible fraction of observed values. These values are significantly softer than the asymptotic value of the theoretical spectrum of $\alpha \sim 0.4$. The reason for the $\alpha$-values to be much softer than expected, is the limitation of the empirical functions to capture the true curvature of the theoretical spectrum. We conclude that more than a $1/4$ of the bursts in the GBM catalogue have at least one time-resolved spectrum, whose $\alpha$-values are consistent with a non-dissipative outflow, releasing its thermal energy at the photosphere. The fraction of spectra consistent with emission from the photosphere will increase even more if dissipation of kinetic energy in the flow occurs below the photosphere.

Many different physical processes have been suggested to explain the prompt gamma-ray emission in gamma-ray bursts (GRBs). Although there are examples of both bursts with photospheric and synchrotron emission origins, these distinct spectral appearances have not been generalized to large samples of GRBs. Here, we search for signatures of the different emission mechanisms in the full Fermi Gamma-ray Space Telescope GBM catalogue. We use Gaussian Mixture Models to cluster bursts according to their parameters from the Band function ($\alpha$, $\beta$, and $E_{pk}$) as well as their fluence and $T_{90}$. We find five distinct clusters. We further argue that these clusters can be divided into bursts of photospheric origin (2/3 of all bursts, divided into 3 clusters) and bursts of synchrotron origin (1/3 of all bursts, divided into 2 clusters). For instance, the cluster that contains predominantly short bursts is consistent of photospheric emission origin. We discuss several reasons that can determine which cluster a burst belongs to: jet dissipation pattern and/or the jet content, or viewing angle.

We study the narrowest spectra expected from gamma-ray bursts. We present an analytical function for the spectrum that is emitted from the photosphere of a radiation-dominated flow that is under acceleration. This is the narrowest possible spectrum and it differs from a Planck function. We also present numerical spectra from photospheres occurring during the transition into the coasting phase of the flow. Using these spectral models, we reanalyse Fermi observations of GRB100507 and GRB101219, which both have been reported to have very narrow spectra. The bursts can be fitted by the spectral models: For GRB101219 the spectrum is consistent with the photosphere occurring below or close to the saturation radius, while for GRB100507 the photosphere position relative to the saturation radius can be determined as a function of time. In the latter case, we find that the photosphere initially occurs in the acceleration phase and thereafter transitions into the coasting phase. We also find that this transition occurs at the same time as the change in observed cooling behaviour: the temperature is close to constant before the break and decays after. We argue that such a transition can be explained by an increasing mass outflow rate. Both analysed bursts thus give strong evidence that the jets are (initially) radiation dominated.

We present timing analysis of {\emph{RXTE}}-PCA and {\emph{INTEGRAL}}-ISGRI observations of X Per between 1998 and 2010. All pulse arrival times obtained from the {\emph{RXTE}}-PCA observations are phase connected and a timing solution is obtained using these arrival times. We update the long-term pulse frequency history of the source by measuring its pulse frequencies using {\emph{RXTE}}-PCA and {\emph{INTEGRAL}}-ISGRI data. From the {\emph{RXTE}}-PCA data, the relation between frequency derivative and X-ray flux suggests accretion via the companion's stellar wind. On the other hand, detection of the transient QPO feature peaking at $\sim 0.2$ Hz suggests the existence of an accretion disc. We find that double break models fit the average power spectra well, which suggests that the source has at least two different accretion flow components dominating the overall flow. From the power spectrum of frequency derivatives, we measure a power law index of $\sim -1$ which implies that on short time scales disc accretion dominates over noise, while on time scales longer than the viscous time scales the noise dominates. From pulse profiles, we find a correlation between pulse fraction and count rate of the source.