We suggest that neutron star mergers eject an ultra-relativistic envelope of mass $m\sim 10^{-7}M_\odot$, which helps explain the gamma-ray burst from GW170817. One ejection mechanism is the ablation of the neutron star surface by the burst of neutrinos in the first $30\,\mu$s of the merger. Another, more efficient, mechanism for inflating the ultra-relativistic envelope is an internal shock in the massive ejecta from the merger. A strong shock is expected if the merger product is a magnetar, which emits a centrifugally accelerated wind. The shock propagates outward through the ejecta and accelerates in its outer layers at radii $r\sim 10^9-10^{10}$cm, launching an ultra-relativistic opaque envelope filled with $\sim 10^4$ photons per nucleon. The Lorentz factor profile of the envelope rises outward and determines its homologous expansion, which adiabatically cools the trapped photons. Once the magnetar loses its differential rotation and collapses into a black hole, a powerful jet forms. It drives a blast wave into the envelope, chasing its outer layers and eventually catching up with the envelope photosphere at $r\sim 10^{12}$cm. The ultra-relativistic photospheric breakout of the delayed blast wave emits a gamma-ray burst in a broad solid angle around the merger axis. This model explains the gamma-ray pulse from merger GW170817 with luminosity $L_\gamma\sim 10^{47}$erg/s, duration $\Delta t_{\rm obs}\sim 0.5$s, and characteristic photon energy $\sim 100$keV. The blast wave Lorentz factor at the envelope photosphere is consistent with $\Gamma\geq 5$ that we derive from the observed light curve of the burst. We suggest future tests of the model.

It has been suggested that the prompt emission in gamma-ray bursts (GRBs) could be described by radiation from the photosphere in a hot fireball. Such models must be tested by directly fitting them to data. In this work we use data from the Fermi Gamma-ray Space Telescope and consider a specific photospheric model, in which kinetic energy of a low-magnetisation outflow is dissipated locally by internal shocks below the photosphere. We construct a table model with a physically motivated parameter space and fit it to time-resolved spectra of the 36 brightest Fermi GRBs with known redshift. We find that about two thirds of the examined spectra cannot be described by the model, as it typically under-predicts the observed flux. However, since the sample is strongly biased towards bright GRBs, we argue that this fraction will be significantly lowered when considering the full population. From the successful fits we find that the model can reproduce the full range of spectral slopes present in the sample. For these cases we also find that the dissipation consistently occurs at a radius of $\sim 10^{12}$ cm and that only a few percent efficiency is required. Furthermore, we find a positive correlation between the fireball luminosity and the Lorentz factor. Such a correlation has been previously reported by independent methods. We conclude that if GRB spectra are due to photospheric emission, the dissipation cannot only be the specific scenario we consider here.

The emission mechanism during the prompt phase in gamma-ray bursts (GRBs) can be investigated through correlations between spectral properties. Here, we revisit the correlation relating the instantaneous flux, $F$, and the photon index below the spectral break, $\alpha$, in individual emission pulses, by studying the 38 most prominent pulses in the Fermi/GBM GRB catalogue. First, we search for signatures of the bias in the determination of $\alpha$ due to the limited spectral coverage (window effect) expected in the synchrotron case. The absence of such a characteristic signature argues against the simplest synchrotron models. We instead find that the observed correlation between $F$ and $\alpha$ can, in general, be described by the relation $F(t) \propto {\rm e}^{k\,\alpha(t)}$, for which the median $k = 3$. We suggest that this correlation is a manifestation of subphotospheric heating in a flow with a varying entropy. Around the peak of the light curve, a large entropy causes the photosphere to approach the saturation radius, leading to an intense emission with a narrow spectrum. As the entropy decreases the photosphere secedes from the saturation radius, and weaker emission with a broader spectrum is expected. This simple scenario naturally leads to a correlated variation of the intensity and spectral shape, covering the observed range.

We examine the possible role of turbulence in feeding the emission of gamma-ray bursts (GRBs). Turbulence may develop in a GRB jet as the result of hydrodynamic or current-driven instabilities. The jet carries dense radiation and the turbulence cascade can be damped by Compton drag, passing kinetic fluid energy to photons through scattering. We identify two regimes of turbulence dissipation: (1) "Viscous" - the turbulence cascade is Compton damped on a scale $\ell_{\rm damp}$ greater than the photon mean free path $\ell_\star$. Then turbulence energy is passed to photons via bulk Comptonization by smooth shear flows on scale $\ell_\star<\ell_{\rm damp}$. (2) "Collisionless" - the cascade avoids Compton damping and extends to microscopic plasma scales much smaller than $\ell_\star$. The collisionless dissipation energizes plasma particles, which radiate the received energy; how the dissipated power is partitioned between particles needs further investigation with kinetic simulations. We show that the dissipation regime switches from viscous to collisionless during the jet expansion, at a critical value of the jet optical depth which depends on the amplitude of turbulence. Turbulent GRB jets are expected to emit nonthermal photospheric radiation. Our analysis also suggests revisions of turbulent Comptonization in black hole accretion disks discussed in previous works.

Gamma-ray bursts (GRBs) are exceptionally bright electromagnetic events occurring daily on the sky. The prompt emission is dominated by X-/$\gamma$-rays. Since their discovery over 50 years ago, GRBs are primarily studied through spectral and temporal measurements. The properties of the emission jets and underlying processes are not well understood. A promising way forward is the development of missions capable of characterising the linear polarisation of the high-energy emission. For this reason, the SPHiNX mission has been developed for a small-satellite platform. The polarisation properties of incident high-energy radiation (50-600 keV) are determined by reconstructing Compton scattering interactions in a segmented array of plastic and Gd$_3$Al$_2$Ga$_3$O$_{12}$(Ce) (GAGG(Ce)) scintillators. During a two-year mission, $\sim$200 GRBs will be observed, with $\sim$50 yielding measurements where the polarisation fraction is determined with a relative error $\leq$10%. This is a significant improvement compared to contemporary missions. This performance, combined with the ability to reconstruct GRB localisation and spectral properties, will allow discrimination between leading classes of emission models.

Internal shocks provide a plausible heating mechanism in the jets of gamma-ray bursts (GRBs). Shocks occurring below the jet photosphere are mediated by radiation. It was previously found that radiation mediated shocks (RMSs) inside GRB jets are inefficient photon producers, and the photons that mediate the RMS must originate from an earlier stage of the explosion. We show that this conclusion is valid only for non-magnetized jets. RMSs that propagate in moderately magnetized plasma develop a collisionless subshock which locally heats the plasma to a relativistic temperature, and the hot electrons emit copious synchrotron photons inside the RMS. We find that this mechanism is generally effective for mildly relativistic shocks and may be the main source of photons observed in GRBs. We derive a simple analytical formula for the generated photon number per proton, $Z$, which gives $Z=10^5$-$10^6$, consistent with observations. The photons are initially injected with low energies, well below the observed GRB peak. Their number is controlled by two main factors: (1) the abundance of electron-positron pairs created in the shock, which is self-consistently calculated, and (2) the upper limit on the brightness temperature of soft radiation set by induced downscattering. The injected soft photons that survive induced downscattering gain energy in the RMS via bulk Comptonization and shape its nonthermal spectrum.

Sub-photospheric shock dissipation is one of the main proposed mechanisms for producing the prompt gamma-ray burst (GRB) emission. Such shocks are mediated by scattering of radiation. We introduce a time dependent, special relativistic code which dynamically couples Monte Carlo radiative transfer to the flow hydrodynamics. The code also self-consistently implements electron-positron pair production and annihilation. We simulate shocks with properties relevant for GRBs and study the steady-state solutions, which are accurate deep below the jet photosphere. The shock generates a power-law photon spectrum through the first-order Fermi mechanism, extending upwards from the typical upstream photon energy. Strong shocks (for which the downstream pressure is much larger than the upstream pressure) have rising $\nu F_\nu$ shock spectra. The spectrum extends up to $\epsilon_{max} \equiv E_{max}/m_e c^2 \sim v^2$ for non-relativistic shocks, where $m_e$ is the electron rest mass and $v$ is the relative speed between the upstream and downstream in units of the speed of light $c$. For mildly relativistic shocks the power law softens at $\epsilon \gtrsim 10^{-1}$ due to Klein-Nishina effects, and shocks with $v\gamma \gtrsim 1$, where $\gamma \equiv (1-v^2)^{-1/2}$, produce electron-positron pairs. As an example, a strong shock with $v\gamma = 3$ and a photon-to-proton ratio of $n_\gamma/n_p = 2 \times 10^5$ has a peak pair-to-proton ratio of $Z_\pm \approx 225$. The main effect of pairs in a steady-state shock is to decrease its spatial width by a factor of $\sim Z_\pm$. The post-shock spectrum thermalizes in the downstream. In absence of emission and absorption processes, kinetic equilibrium at temperature $\theta_d \equiv kT_d/m_e c^2 \approx \epsilon_d/3$ is reached at an optical depth of $\tau \gg \theta_d^{-1}$ behind the shock, where $\epsilon_d$ is the average downstream photon energy.

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.

The MeV spectral peak of gamma-ray bursts (GRBs) is best explained as photospheric emission from a dissipative relativistic jet. The observed non-blackbody spectrum shows that sub-photospheric dissipation involves both thermal plasma heating and injection of nonthermal particles, which quickly cool through inverse Compton scattering and emission of synchrotron radiation. Synchrotron photons emitted around and above the photosphere are predicted to dominate the low-energy part of the GRB spectrum, starting from roughly a decade in energy below the MeV peak. We show that this leads to a unique polarization signature: a rise in GRB polarization toward lower energies. We compute the polarization degree of GRB radiation as a function of photon energy for a generic jet model, and show the predictions for GRBs 990123, 090902B and 110721A. The expected polarization is significant in the X-ray band, in particular for bursts similar to GRB 090902B. Radiation in the MeV peak (and at higher energies) is unpolarized as long as the jet is approximately uniform on angular scales $\delta\theta \gtrsim \Gamma^{-1}$ where $\Gamma$ is the bulk Lorentz factor of the jet.

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 consider the polarization properties of photospheric emission originating in jets consisting of a highly relativistic core of opening angle theta_j and Lorentz factor Gamma_0, and a surrounding shear layer where the Lorentz factor is decreasing as a power law of index p with angle from the jet axis. We find significant degrees of linear polarization for observers located at viewing angles theta_v \gtrsim theta_j. In particular, the polarization degree of emission from narrow jets (theta_j \approx 1/Gamma_0) with steep Lorentz factor gradients (p \gtrsim 4) reaches \sim 40%. The angle of polarization may shift by pi/2 for time-variable jets. The spectrum below the thermal peak of the polarized emission appears non-thermal due to aberration of light, without the need for additional radiative processes or energy dissipation. Furthermore, above the thermal peak a power law of photons forms due to Comptonization of photons that repeatedly scatter between regions of different Lorentz factor before escaping. We show that polarization degrees of a few tens of percent and broken power law spectra are natural in the context of photospheric emission from structured jets. Applying the model to gamma-ray bursts, we discuss expected correlations between the spectral shape and the polarization degree of the prompt emission.

{\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.

Relativistic outflows in the form of jets are common in many astrophysical objects. By their very nature, jets have angle dependent velocity profiles, Gamma = Gamma(r, theta, phi), where Gamma is the outflow Lorentz factor. In this work we consider photospheric emission from non-dissipative jets with various Lorentz factor profiles, of the approximate form Gamma \approx Gamma_0/[(theta/theta_j)^p + 1], were theta_j is the characteristic jet opening angle. In collimated jets, the observed spectrum depends on the viewing angle, theta_v. We show that for narrow jets (theta_j Gamma_0 \lesssim few), the obtained low energy photon index is alpha \approx -1 (dN/dE \propto E^alpha), independent of viewing angle, and weakly dependent on the Lorentz factor gradient (p). A similar result is obtained for wider jets observed at theta_v \approx theta_j. This result is surprisingly similar to the average low energy photon index seen in gamma-ray bursts. For wide jets (theta_j Gamma_0 \gtrsim few) observed at theta_v \ll theta_j, a multicolor blackbody spectrum is obtained. We discuss the consequences of this theory on our understanding of the prompt emission in gamma-ray bursts.

We analyze the spectral evolution of GRB 090902B and show that subphotospheric dissipation can explain both the spectra and the spectral evolution. The emission from a GRB photosphere can give rise to a variety of spectral shapes. The spectrum can have a shape close to that of a Planck function (as is observed during the first half of GRB090902B) or be broadened, resembling a typical Band function (as is observed during the second half of GRB090902B). The shape mainly depends on the strength and location of the dissipation in the jet, the ratio of the energy densities of thermal photons and of the electrons at the dissipation site, as well as on the strength of the magnetic field. We further discuss numerical models of the dissipation and relate these to the observed spectra.

The emission from a gamma-ray burst (GRB) photosphere can give rise to a variety of spectral shapes. The spectrum can retain the shape of a Planck function or it can be broadened and have the shape of a Band function. This fact is best illustrated by studying GRB090902B: The main gamma-ray spectral component is initially close to a Planck function, which can only be explained by emission from the jet photosphere. Later, the same component evolves into a broader Band function. This burst thus provides observational evidence that the photosphere can give rise to a non-thermal spectrum. We show that such a broadening is most naturally explained by subphotospheric dissipation in the jet. The broadening mainly depends on the strength and location of the dissipation, on the magnetic field strength, and on the relation between the energy densities of thermal photons and of the electrons. We suggest that the evolution in spectral shape observed in GRB090902B is due to a decrease of the bulk Lorentz factor of the flow, leading to the main dissipation becoming subphotospheric. Such a change in the flow parameters can also explain the correlation observed between the peak energy of the spectrum and low-energy power law slope, alpha, a correlation commonly observed in GRBs. We conclude that photospheric emission could indeed be a ubiquitous feature during the prompt phase in GRBs and play a decisive role in creating the diverse spectral shapes and spectral evolutions that are observed.

We present a time-resolved spectral analysis of the bright, long GRB 061007 (z=1.261) using Swift BAT and Suzaku WAM data. We find that the prompt emission of GRB 061007 can be equally well explained by a photospheric component together with a power law as by a Band function, and we explore the implications of the former model. The photospheric component, which we model with a multicolour blackbody, dominates the emission and has a very stable shape throughout the burst. This component provides a natural explanation for the hardness-intensity correlation seen within the burst and also allows us to estimate the bulk Lorentz factor and the radius of the photosphere. The power-law component dominates the fit at high energies and has a nearly constant slope of -1.5. We discuss the possibility that this component is of the same origin as the high-energy power laws recently observed in some Fermi LAT bursts.

The Fermi Gamma-ray Space Telescope observed the bright and long GRB090902B, lying at a redshift of z = 1.822. Together the Large Area Telescope (LAT) and the Gamma-ray Burst Monitor (GBM) cover the spectral range from 8 keV to >300 GeV is covered. Here we show that the prompt burst spectrum is consistent with emission from the jet photosphere combined with non-thermal emission described by a single power-law with photon index -1.9. The photosphere gives rise to a strong quasi-blackbody spectrum which is somewhat broader than a single Planck function and has a characteristic temperature of ~290keV. We derive the photospheric radius Rph = (1.1 \pm 0.3) x 10^12 Y^{1/4} cm and the bulk Lorentz factor of the flow, which is found to vary by a factor of two and has a maximal value of Gamma = 750 Y^{1/4}. Here Y is the ratio between the total fireball energy and the energy emitted in the gamma-rays. We find that during the first quarter of the prompt phase the photospheric emission dominates, which explains the delayed onset of the observed flux in the LAT compared to the GBM. We model the photospheric emission with a multi-color blackbody and its shape indicates that the photospheric radius increases at higher latitudes. We interpret the broad band emission as synchrotron emission at R ~ 4x10^15 cm. Our analysis emphasize the importance of having high temporal resolution when performing spectral analysis on GRBs, since there is strong spectral evolution.