Long-duration gamma-ray bursts (GRBs) originate from ultra-relativistic jets launched from the collapsing cores of dying massive stars. They are characterised by an initial phase of bright and highly variable radiation in the keV-MeV band that is likely produced within the jet and lasts from milliseconds to minutes, known as the prompt emission. Subsequently, the interaction of the jet with the external medium generates external shock waves, responsible for the afterglow emission, which lasts from days to months, and occurs over a broad energy range, from the radio to the GeV bands. The afterglow emission is generally well explained as synchrotron radiation by electrons accelerated at the external shock. Recently, an intense, long-lasting emission between 0.2 and 1 TeV was observed from the GRB 190114C. Here we present the results of our multi-frequency observational campaign of GRB~190114C, and study the evolution in time of the GRB emission across 17 orders of magnitude in energy, from $5\times10^{-6}$ up to $10^{12}$\,eV. We find that the broadband spectral energy distribution is double-peaked, with the TeV emission constituting a distinct spectral component that has power comparable to the synchrotron component. This component is associated with the afterglow, and is satisfactorily explained by inverse Compton upscattering of synchrotron photons by high-energy electrons. We find that the conditions required to account for the observed TeV component are not atypical, supporting the possibility that inverse Compton emission is commonly produced in GRBs.

The observations of GW170817/GRB170817A have confirmed that the coalescence of a neutron-star binary is the progenitor of a short gamma-ray burst. In the standard picture of a short gamma-ray burst, a collimated highly relativistic outflow is launched after merger and it successfully breaks out from the surrounding ejected matter. Using initial conditions inspired from numerical-relativity binary neutron-star merger simulations, we have performed general-relativistic hydrodynamic (HD) and magnetohydrodynamic (MHD) simulations in which the jet is launched and propagates self-consistently. The complete set of simulations suggests that: (i) MHD jets have an intrinsic energy and velocity polar structure with a ``hollow core'' subtending an angle $\theta_{\rm core}\approx4^{\circ}-5^{\circ}$ and an opening angle of $\theta_{\rm jet}\gtrsim10^{\circ}$; (ii) MHD jets eject significant amounts of matter and two orders of magnitude more than HD jets; (iii) the energy stratification in MHD jets naturally yields the power-law energy scaling $E(>\Gamma\beta)\propto(\Gamma\beta)^{-4.5}$; (iv) MHD jets provide fits to the afterglow data from GRB170817A that are comparatively better than those of the HD jets and without free parameters; (v) finally, both of the best-fit HD/MHD models suggest an observation angle $\theta_{\rm obs} \simeq 21^{\circ}$ for GRB170817A.

GRB 170817A/GW 170817 is the first GRB clearly viewed far from the GRB jet's symmetry axis. Its afterglow was densely monitored over a wide range of frequencies and times. It has been modeled extensively, primarily numerically, and although this endeavour was very fruitful, many of the underlying model parameters remain undetermined. We provide analytic modelling of GRB afterglows observed off-axis, considering jets with a narrow core (of half-opening angle $\theta_{\rm c}$) and power-law wings in energy per unit solid angle ($\epsilon=\epsilon_c\Theta^{-a}$ where $\Theta=[1+(\theta/\theta_{\rm c})^2]^{1/2}$) and initial specific kinetic energy ($\Gamma_0-1=[\Gamma_{\rm c,0}-1]\Theta^{-b}$), as well as briefly discuss Gaussian jets. Our study reveals qualitatively different types of lightcurves that can be viewed in future off-axis GRBs, with either single or double peaks, depending on the jet structure and the viewing angle. Considering the lightcurve shape rather than the absolute normalizations of times and / or fluxes, removes the dependence of the lightcurve on many of the highly degenerate burst parameters. This study can be easily used to determine the underlying jet structure, significantly reduce the effective parameter space for numerical fitting attempts and provide physical insights. As an illustration, we show that for GRB 170817A, there is a strong correlation between the allowed values of $\Gamma_{\rm c,0}$ and $b$, leading to a narrow strip of allowed solutions in the $\Gamma_{\rm c,0}$-$b$ plane above some minimal values $\Gamma_{\rm c,0}\gtrsim 40, b\gtrsim1.2$. Furthermore, the Lorentz factor of the material dominating the early lightcurve can be constrained by three independent techniques to be $\Gamma_{0}(\theta_{\rm min,0})\approx5-7$.

We report on the observations of gamma-ray burst (GRB) 190114C by the Fermi Gamma-ray Space Telescope and the Neil Gehrels Swift Observatory. The early-time observations reveal multiple emission components that evolve independently, with a delayed power-law component that exhibits significant spectral attenuation above 40 MeV in the first few seconds of the burst. This power-law component transitions to a harder spectrum that is consistent with the afterglow emission observed at later times. This afterglow component is clearly identifiable in the GBM and BAT light curves as a slowly fading emission component on which the rest of the prompt emission is superimposed. As a result, we are able to constrain the transition from internal shock to external shock dominated emission. We find that the temporal and spectral evolution of the broadband afterglow emission can be well modeled as synchrotron emission from a forward shock propagating into a wind-like circumstellar environment and find that high-energy photons observed by Fermi LAT are in tension with the theoretical maximum energy that can be achieved through synchrotron emission from a shock. These violations of the maximum synchrotron energy are further compounded by the detection of very high energy (VHE) emission above 300 GeV by MAGIC concurrent with our observations. We conclude that the observations of VHE photons from GRB 190114C necessitates either an additional emission mechanism at very high energies that is hidden in the synchrotron component in the LAT energy range, an acceleration mechanism that imparts energy to the particles at a rate that is faster than the electron synchrotron energy loss rate, or revisions of the fundamental assumptions used in estimating the maximum photon energy attainable through the synchrotron process.

We present ALMA 97.5 GHz total intensity and linear polarization observations of the mm-band afterglow of GRB 190114C spanning 2.2 to 5.2 hours after the burst. We detect linear polarization at the $\approx 5\,\sigma$ level, decreasing from $\Pi=(0.87\pm0.13)\%$ to $(0.60\pm0.19)\%$, and evolving in polarization position angle from $(10\pm5)^\circ$ to $(-44\pm12)^\circ$ during the course of the observations. This represents the first detection of polarized millimeter emission in a $\gamma$-ray burst. We show that the optical and X-ray observations between $0.03$ days and $\sim0.3$ days are consistent with a fast cooling forward shock expanding into a wind environment. However, the optical observations at $\lesssim0.03$ days, as well as the radio and millimeter observations arise from a separate component, which we interpret as emission from the reverse-shocked ejecta. Using the measured linear polarization, we constrain the coherence scale of tangled magnetic fields in the ejecta to an angular size of $\theta_{\rm B} \approx10^{-3}$ radian, while the rotation of the polarization angle rules out the presence of large scale, ordered axisymmetric magnetic fields, and in particular a large scale toroidal field, in the jet.

Gamma-ray burst afterglows arise from relativistic collisionless shocks in which the postshock tangled magnetic field $\vec{B}$ is produced by the two-stream and/or Weibel instabilities on plasma skin-depth scales $(c/\omega_p)$. The field is expected to be oriented predominantly within the shock plane ($B_{\perp}$; transverse to the shock normal, $\hat{n}_{\rm{sh}}$), and is often approximated to be completely within it ($B_\parallel\equiv\hat{n}_{\rm{sh}}\,\cdot\,\vec{B}=0$). Current 2D/3D particle-in-cell simulations are limited to short timescales and box sizes $\lesssim10^4(c/\omega_p)\ll R/\Gamma_{\rm{sh}}$ much smaller than the shocked region's comoving width, and cannot probe the asymptotic downstream $\vec{B}$ structure. We constrain the latter using the linear polarization upper limit, $|\Pi|<12\%$, on the radio afterglow of GW170817/GRB170817A. Afterglow polarization depends on the jet's angular structure, our viewing angle, and the $\vec{B}$ structure. In GW170817/GRB170817A the latter can be tightly constrained since the former two are constrained from observations. We model $\vec{B}$ as an isotropic field in 3D that is stretched along $\hat{n}_{\rm{sh}}$ by a factor $\xi\equiv{}B_{\parallel}/B_{\perp}$, whose initial value $\xi_f\equiv{}B_{\parallel,f}/B_{\perp,f}$ describes the field that survives downstream on plasma scales $\ll{}R/\Gamma_{\rm{sh}}$. We calculate $\Pi(\xi_f)$ by integrating over the shocked volume for core-dominated structured jets, with a local Blandford-McKee self-similar radial profile. We find that independent of the exact jet structure, $\vec{B}$ has a finite, but initially sub-dominant, parallel component: $0.57\lesssim\xi_f\lesssim0.89$, making it less anisotropic. While this motivates numerical studies of the asymptotic $\vec{B}$ structure in relativistic collisionless shocks, it may be consistent with turbulence amplified magnetic field.

Despite being hard to measure, GRB prompt $\gamma$-ray emission polarization is a valuable probe of the dominant emission mechanism and the outflow's composition and angular structure. During the prompt emission the outflow is ultra-relativistic with Lorentz factors $\Gamma\gg1$. We describe in detail the linear polarization properties of various emission mechanisms: synchrotron radiation from different magnetic field structures (ordered: toroidal $B_{\rm tor}$ or radial $B_\parallel$, and random: normal to the radial direction $B_\perp$), Compton drag, and photospheric emission. We calculate the polarization for different GRB jet angular structures (e.g. top-hat, Gaussian, power-law) and viewing angles $\theta_{\rm obs}$. Synchrotron with $B_\perp$ can produce large polarizations, up to $25\%\lesssim\Pi\lesssim45\%$, for a top-hat jet but only for lines of sight just outside the jet's sharp edge. The same also holds for Compton drag, albeit with a slightly higher overall $\Pi$. Moreover, we demonstrate how $\Gamma$-variations during the GRB or smoother jet edges would significantly reduce $\Pi$. We construct a semi-analytic model for non-dissipative photospheric emission from structured jets. Such emission can produce up to $\Pi\lesssim15\%$ with reasonably high fluences, but this requires steep gradients in $\Gamma(\theta)$. A polarization of $50\%\lesssim\Pi\lesssim65\%$ can robustly be produced only by synchrotron emission from a transverse magnetic field ordered on angles $\gtrsim\!1/\Gamma$ around our line of sight (like a global toroidal field). Therefore, such a model would be strongly favored even by a single secure measurement within this range. We find that such a model would also be favored if $\Pi\gtrsim20\%$ is measured in most GRBs within a large enough sample, by deriving the polarization distribution for our different emission and jet models.

The afterglow of GRB$\,$170817A/GW$\,$170817 was very unusual, slowly rising as $F_\nu\propto{}t_{\rm{}obs}^{0.8}\nu^{-0.6}$, peaking at $t_{\rm{obs,pk}}\sim\,150\;$days, and sharply decaying as $\sim{}t_{\rm{}obs}^{-2.2}$. VLBI observations revealed an unresolved radio afterglow image whose flux centroid moved superluminally with $v_{\rm{app}}\approx4c$, clearly indicating that the afterglow was dominated by a relativistic jet's compact core. Different jet angular structures explained the afterglow lightcurves: Gaussian and steep power-law profiles with narrow core angles $\theta_c\lesssim5^\circ$ and larger viewing angles $\theta_{\rm{}obs}/\theta_c\sim3-5$. However, a top-hat jet (with sharp edges at $\theta=\theta_0$) was ruled out since it appeared to produce an early flux rise much steeper than observed. Using 2D relativistic hydrodynamic simulations we show that the initial steep flux rise is an artifact caused by the simulation's finite start time, $t_0$, missing its flux contributions from $t<t_0$ and sometimes "compensated" using an analytic top-hat jet. While an initially top-hat jet is not very physical, such simulations are particularly useful at $t_{\rm{}obs}\gtrsim{}t_{\rm{obs,pk}}$ when the afterglow emission is dominated by the jet's core and becomes insensitive to its exact initial angular profile if it drops off sharply outside of the core. We demonstrate that an initially top-hat jet fits GW$\,$170817/GRB$\,$170817A's afterglow lightcurves and flux centroid motion at $t_{\rm{}obs}\gtrsim{}t_{\rm{obs,pk}}$, for $\theta_{\rm{}obs}/\theta_0\approx3$ and may also fit the earlier lightcurves for $\Gamma_0=\Gamma(t_0)\gtrsim10^{2.5}$. We analytically express the degeneracies between the model parameters, and find a minimal jet energy of $E_{\rm{}min}\approx5.3\times10^{48}\;$erg and circum-burst medium density of $n_{\min}\approx5.3\times10^{-6}~{\rm cm}^{-3}$.

The main hard pulse of prompt gamma-ray emission in GRB$\,$170817A had a duration of $\sim0.5\,{\rm s}$ and its onset was delayed with respect to the gravitational-wave chirp signal by $t_{\rm del} \approx 1.74\,{\rm s}$. Detailed follow-up of the subsequent broadband kilonova emission revealed a two-component ejecta -- a lanthanide-poor ejecta with mass $M_{\rm ej,blue}\approx0.025\,M_\odot$ that powered the early but rapidly fading blue emission and a lanthanide-rich ejecta with mass $M_{\rm ej,red}\approx 0.04\,M_\odot$ that powered the longer lasting redder emission. Both the prompt gamma-ray onset delay and the existence of the blue ejecta with modest electron fraction, $0.2\lesssim Y_e\lesssim0.3$, can be explained if the collapse to a black hole was delayed by the formation of a hypermassive neutron star (HMNS). Here, we determine the survival time of the merger remnant by combining two different constraints, namely, the time needed to produce the requisite blue-ejecta mass and that necessary for the relativistic jet to bore its way out of the expanding ejecta. In this way, we determine that the remnant of GW170817 must have collapsed to a black hole after $t_{\rm coll}=0.98_{-0.26}^{+0.31}\,{\rm s}$. We also discuss how future detections and the delays between the gravitational and electromagnetic emissions can be used to constrain the properties of the merged object.

In this White Paper we present the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to Observatory Science targets. These include flaring stars, supernova remnants, accreting white dwarfs, low and high mass X-ray binaries, radio quiet and radio loud active galactic nuclei, tidal disruption events, and gamma-ray bursts. eXTP will be excellently suited to study one common aspect of these objects: their often transient nature. Developed by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Science, the eXTP mission is expected to be launched in the mid 2020s.

We analyze the prompt emission of GRB 100724B and GRB 160509A, two of the brightest Gamma-Ray Bursts (GRBs) observed by Fermi at $\lesssim{\rm MeV}$ energies but surprisingly faint at $\gtrsim100\;$MeV energies. Time-resolved spectroscopy reveals a sharp high-energy cutoff at energies $E_c\sim20-60\;$MeV for GRB~100724B and $E_c\sim80-150\;$MeV for GRB~160509A. We first characterize phenomenologically the cutoff and its time evolution. We then fit the data to two models where the high-energy cutoff arises from intrinsic opacity to pair production within the source ($\tau_{\gamma\gamma}$): (i) a Band spectrum with $\tau_{\gamma\gamma}$ from the internal-shocks motivated model of Granot et al. (2008), and (ii) the photospheric model of Gill & Thompson (2014). Alternative explanations for the cutoff, such as an intrinsic cutoff in the emitting electron energy distribution, appear to be less natural. Both models provide a good fit to the data with very reasonable physical parameters, providing a direct estimate of bulk Lorentz factors in the range $\Gamma\sim 100-400$, on the lower end of what is generally observed in Fermi GRBs. Surprisingly, their lower cutoff energies $E_c$ compared to other Fermi-LAT GRBs arise not predominantly from the lower Lorentz factors, but also at a comparable level from differences in variability time, luminosity, and high-energy photon index. Finally, particularly low $E_c$ values may prevent detection by Fermi-LAT, thus introducing a bias in the Fermi-LAT GRB sample against GRBs with low Lorentz factors or variability times.

The X-ray to radio afterglow emission of GRB 170817A/GW 170817 so far scales as $F_\nu\propto\nu^{-0.6}t^{0.8}$ with observed frequency and time, consistent with a single power-law segment of the synchrotron spectrum from the external shock going into the ambient medium. This requires the effective isotropic equivalent afterglow shock energy in the visible region to increase as $\sim t^{1.7}$. The two main channels for such an energy increase are (i) \emph{radial}: more energy carried by slower material (in the visible region) gradually catches up with the afterglow shock and energizes it, and (ii) \emph{angular}: more energy in relativistic outflow moving at different angles to our line of sight, whose radiation is initially beamed away from us but its beaming cone gradually reaches our line of sight as it decelerates. One cannot distinguish between these explanations (or combinations of them) using only the X-ray to radio $F_\nu(t)$. Here we demonstrate that the most promising way to break this degeneracy is through afterglow imaging and polarization, by calculating the predicted evolution of the afterglow image (its size, shape and flux centroid) and linear polarization $\Pi(t)$ for different angular and/or radial outflow structures that fit $F_\nu(t)$. We consider two angular profiles -- a Gaussian and a narrow core with power-law wings in energy per solid angle, as well as a (cocoon motivated) (quasi-) spherical flow with radial velocity profile. For a jet viewed off-axis (and a magnetic field produced in the afterglow shock) $\Pi(t)$ peaks when the jet's core becomes visible, at $\approx2t_p$ where the lightcurve peaks at $t_p$, and the image can be elongated with aspect ratios$\;\gtrsim2$. A quasi-spherical flow has an almost circular image and a much lower $\Pi(t)$ (peaking at $\approx t_p$) and flux centroid displacement.

The highly luminous and variable prompt emission in Gamma-Ray Bursts (GRBs) arises in an ultra-relativistic outflow. The exact underlying radiative mechanism shaping its non-thermal spectrum is still uncertain, making it hard to determine the outflow's bulk Lorentz factor $\Gamma$. GRBs with spectral cutoff due to pair production ($\gamma\gamma\to e^+e^-$) at energies $E_c\gtrsim10\;$MeV are extremely useful for inferring $\Gamma$. We find that when the emission region has a high enough compactness, then as it becomes optically thick to scattering, Compton downscattering by non-relativistic $e^\pm$-pairs can shift the spectral cutoff energy well below the self-annihilation threshold, $E_{\rm sa}=\Gamma m_ec^2/(1+z)$. We treat this effect numerically and show that $\Gamma$ obtained assuming $E_c=E_{\rm sa}$ can under-predict its true value by as much as an order of magnitude.

The short-duration ($\lesssim2\;$s) GRB 170817A in the nearby ($D=40\;$Mpc) elliptical galaxy NGC 4993 is the first electromagnetic counterpart of the first gravitational wave (GW) detection of a binary neutron-star (NS-NS) merger. It was followed by optical, IR, and UV emission from half a day up to weeks after the event, as well as late time X-ray and radio emission. The early UV, optical, and IR emission showed a quasi-thermal spectrum suggestive of radioactive-decay powered kilonova-like emission. Comparison to kilonova models favors the formation of a short-lived ($\sim1\;$s) hypermassive NS, which is also supported by the $\Delta t\approx1.74\;$s delay between the GW chirp signal and the prompt GRB onset. However, the late onset of the X-ray (8.9$\;$days) and radio (16.4$\;$days) emission, together with the low isotropic equivalent $\gamma$-ray energy output ($E_{\rm\gamma,iso}\approx5\times10^{46}\;$erg), strongly suggest emission from a narrow relativistic jet viewed off-axis. Here we set up a general framework for off-axis GRB jet afterglow emission, comparing analytic and numerical approaches, and showing their general predictions for short-hard GRBs that accompany binary NS mergers. The prompt GRB emission suggests a viewing angle well outside the jet's core, and we compare the afterglow lightcurves expected in such a case to the X-ray to radio emission from GRB 170817A. We fit an afterglow off-axis jet model to the X-ray and radio data and find that the observations are explained by a viewing angle $\theta_{\rm obs}\approx16^\circ-26^\circ$, GRB jet energy $E\sim10^{48.5}-10^{49.5}~{\rm erg}$, and external density $n\sim10^{-5}-10^{-1}~{\rm cm}^{-3}$ for a $\xi_e\sim 0.1$ non-thermal electron acceleration efficiency.

We study the linear and non-linear development of the Kruskal-Schwarzchild Instability in a relativisitically expanding striped wind. This instability is the generalization of Rayleigh-Taylor instability in the presence of a magnetic field. It has been suggested to produce a self-sustained acceleration mechanism in strongly magnetized outflows found in active galactic nuclei, gamma-ray bursts, and micro-quasars. The instability leads to magnetic reconnection, but in contrast with steady-state Sweet-Parker reconnection, the dissipation rate is not limited by the current layer's small aspect ratio. We performed two-dimensional (2D) relativistic magneto-hydrodynamic (RMHD) simulations featuring two cold and highly magnetized ($1\leq\sigma\leq10^{3}$) plasma layers with an anti-parallel magnetic field separated by a thin layer of relativistically hot plasma with a local effective gravity induced by the outflow's acceleration. Our simulations show how the heavier relativistically hot plasma in the reconnecting layer drips out and allows oppositely oriented magnetic field lines to reconnect. The instability's growth rate in the linear regime matches the predictions of linear stability analysis. We find turbulence rather than an ordered bulk flow near the reconnection region, with turbulent velocities up to $\sim0.1$c, largely independent of model parameters. However, the magnetic energy dissipation rate is found to be much slower, corresponding to an effective ordered bulk velocity inflow into the reconnection region $v_{\rm in}=\beta_{\rm in}c$, of $10^{-3}\lesssim\beta_{\rm in}\lesssim 5\times10^{-3}$. This occurs due to the slow evacuation of hot plasma from the current layer, largely because of the Kelvin-Helmholtz instability experienced by the dripping plasma. 3D RMHD simulations are needed to further investigate the non-linear regime.

The first, long awaited, detection of a gravitational wave (GW) signal from the merger of a binary neutron-star (NS-NS) system was finally achieved (GW$\,$170817), and was also accompanied by an electromagnetic counterpart -- the short-duration GRB 170817A. It occurred in the nearby ($D\approx40\;$Mpc) elliptical galaxy NGC$\,$4993, and showed optical, IR and UV emission from half a day up to weeks after the event, as well as late time X-ray (at $\geq 8.9\;$days) and radio (at $\geq 16.4\;$days) emission. There was a delay of $\Delta t \approx 1.74\;$s between the GW merger chirp signal and the prompt-GRB emission onset, and an upper limit of $\theta_{\rm obs}<28^\circ$ was set on the viewing angle w.r.t the jet's symmetry axis from the GW signal. In this letter we examine some of the implications of these groundbreaking observations. The delay $\Delta t$ sets an upper limit on the prompt-GRB emission radius, $R_\gamma\lesssim 2c\Delta t/(\theta_{\rm obs}-\theta_0)^2$, for a jet with sharp edges at an angle $\theta_0<\theta_{\rm obs}$. GRB 170817A's relatively low isotropic equivalent $\gamma$-ray energy-output may suggest a viewing angle slightly outside the jet's sharp edge, $\theta_{\rm{}obs}-\theta_0\sim(0.05-0.1)(\Gamma/100)^{-1}$, but its peak $\nu{}F_\nu$ photon energy and afterglow emission suggest instead that the jet does not have sharp edges and the prompt emission was dominated by less energetic material along our line of sight, at $\theta_{\rm{}obs}\gtrsim 2\theta_0$. Finally, we consider the type of remnant that is produced by the NS-NS merger and find that a relatively long-lived ($>2\;$s) massive NS is strongly disfavored, while a hyper-massive NS of lifetime $\sim1\;$s appears to be somewhat favored over the direct formation of a black hole.

We analyzed broad-band X-ray and radio data of the magnetar SGR J1935+2154 taken in the aftermath of its 2014, 2015, and 2016 outbursts. The source soft X-ray spectrum <10 keV is well described with a BB+PL or 2BB model during all three outbursts. NuSTAR observations revealed a hard X-ray tail, $\Gamma=0.9$, extending up to 79 keV, with flux larger than the one detected <10 keV. Imaging analysis of Chandra data did not reveal small-scale extended emission around the source. Following the outbursts, the total 0.5-10 keV flux from SGR J1935+2154 increased in concordance to its bursting activity, with the flux at activation onset increasing by a factor of $\sim7$ following its strongest June 2016 outburst. A Swift/XRT observation taken 1.5 days prior to the onset of this outburst showed a flux level consistent with quiescence. We show that the flux increase is due to the PL or hot BB component, which increased by a factor of $25$ compared to quiescence, while the cold BB component $kT=0.47$ keV remained more or less constant. The 2014 and 2015 outbursts decayed quasi-exponentially with time-scales of $\sim40$ days, while the stronger May and June 2016 outbursts showed a quick short-term decay with time-scales of $\sim4$ days. Our Arecibo radio observations set the deepest limits on the radio emission from a magnetar, with a maximum flux density limit of 14 $\mu$Jy for the 4.6 GHz observations and 7 $\mu$Jy for the 1.4 GHz observations. We discuss these results in the framework of the current magnetar theoretical models.

The origin of prompt emission from gamma ray bursts remains to be an open question. Correlated prompt optical and gamma-ray emission observed in a handful of GRBs strongly suggests a common emission region, but failure to adequately fit the broadband GRB spectrum prompted the hypothesis of different emission mechanisms for the low- and high-energy radiations. We demonstrate that our multi-component model for GRB gamma-ray prompt emission provides an excellent fit to GRB 110205A from optical to gamma-ray energies. Our results show that the optical and highest gamma-ray emissions have the same spatial and spectral origin, which is different from the bulk of the X- and softest gamma-ray radiation. Finally, our accurate redshift estimate for GRB 110205A demonstrates promise for using GRBs as cosmological standard candles.

We have compiled the most comprehensive burst sample from magnetar 4U 0142+61, comprising 27 bursts from its three burst-active episodes in 2011, 2012 and the latest one in 2015 observed with Swift/BAT and Fermi/GBM. Bursts from 4U 0142+61 morphologically resemble typical short bursts from other magnetars. However, 4U 0142+61 bursts are less energetic compared to the bulk of magnetar bursts. We uncovered an extended tail emission following a burst on 2015 February 28, with a thermal nature, cooling over a time-scale of several minutes. During this tail emission, we also uncovered pulse peak phase aligned X-ray bursts , which could originate from the same underlying mechanism as that of the extended burst tail, or an associated and spatially coincident but different mechanism.

The first wind nebula around a magnetar was recently discovered in X-rays around Swift~J1834.9$-$0846. We study this magnetar's global energetics and the properties of its particle wind or outflows. At a distance of $\sim4\;$kpc, Swift~J1834.9$-$0846 is located at the center of the supernova remnant (SNR) W41 whose radius is $\sim 19\;$pc, an order of magnitude larger than that of the X-ray nebula ($\sim2\;$pc). The association with SNR W41 suggests a common age of $\sim5-100\;$kyr, while its spin-down age is $4.9$~kyr. A small natal kick velocity may partly explain why a wind nebula was detected around this magnetar but not around other magnetars, most of which appear to have larger kick velocities and may have exited their birth SNR. We find that the GeV and TeV source detected by Fermi/LAT and H.E.S.S., respectively, of radius $\sim11\;$pc is most likely of hadronic origin. The dynamics and internal structure of the nebula are examined analytically to explain the nebula's current properties. Its size may naturally correspond to the diffusion-dominated cooling length of the X-ray emitting $e^+e^-$ pairs. This may also account for the spectral softening of the X-ray emission from the nebula's inner to outer parts. Analysis of the X-ray synchrotron nebula implies that (i) the nebular magnetic field is $\gtrsim 11\;\mu$G (and likely $\lesssim30\;\mu$G), and (ii) the nebula is not powered predominantly by the magnetar's quiescent spin-down-powered MHD wind, but by other outflows that contribute most of its energy. The latter are most likely associated with the magnetar's bursting activity, and possibly dominated by outflows associated with its past giant flares. The energy source for the required outflows cannot be the decay of the magnetar's dipole field alone, and is most likely the decay of its much stronger internal magnetic field.

We use measurements from the Rosetta plasma consortium (RPC) Langmuir probe (LAP) and mutual impedance probe (MIP) to study the spatial distribution of low-energy plasma in the near-nucleus coma of comet 67P/Churyumov-Gerasimenko. The spatial distribution is highly structured with the highest density in the summer hemisphere and above the region connecting the two main lobes of the comet, i.e. the neck region. There is a clear correlation with the neutral density and the plasma to neutral density ratio is found to be about 1-2x10^-6, at a cometocentric distance of 10 km and at 3.1 AU from the sun. A clear 6.2 h modulation of the plasma is seen as the neck is exposed twice per rotation. The electron density of the collisonless plasma within 260 km from the nucleus falls of with radial distance as about 1/r. The spatial structure indicates that local ionization of neutral gas is the dominant source of low-energy plasma around the comet.

We present the Fermi Gamma-ray Burst Monitor (GBM) and Large Area Telescope (LAT) observations of the LIGO binary black hole merger event GW151226 and candi- date LVT151012. No candidate electromagnetic counterparts were detected by either the GBM or LAT. We present a detailed analysis of the GBM and LAT data over a range of timescales from seconds to years, using automated pipelines and new techniques for char- acterizing the upper limits across a large area of the sky. Due to the partial GBM and LAT coverage of the large LIGO localization regions at the trigger times for both events, dif- ferences in source distances and masses, as well as the uncertain degree to which emission from these sources could be beamed, these non-detections cannot be used to constrain the variety of theoretical models recently applied to explain the candidate GBM counterpart to GW150914.

We report on the analysis of two deep XMM-Newton observations of the magnetar Swift J1834.9-0846 and its surrounding extended emission taken in March 2014 and October 2014, 2.5 and 3.1 years after the source went into outburst. The magnetar is only weakly detected in the first observation with an absorption corrected flux $F_{\rm 0.5-10 keV}\approx4\times10^{-14}$ erg s$^{-1}$ cm$^{-2}$, and a $3\sigma$ upper limit during the second observation of about $3\times10^{-14}$ erg s$^{-1}$ cm$^{-2}$. This flux level is more than 3 orders of magnitude lower than the flux measured at the outburst onset on September 2011. The extended emission, centered at the magnetar position and elongated towards the south-west, is clearly seen in both observations; it is best fit by a highly absorbed power-law (PL), with a hydrogen column density of $N_{\rm H}=8.0\times10^{22}$ cm$^{-2}$ and PL photon index $\Gamma=2.2\pm0.2$. Its flux is constant between the two observations at $F_{\rm 0.5-10 keV}=1.3\times10^{-12}$ erg s$^{-1}$ cm$^{-2}$. We find no statistically significant changes in the spectral shape or the flux of this extended emission over a period of 9 years from 2005 to 2014. These new results strongly support the extended emission nature as a wind nebula and firmly establish Swift J1834.9-0846 as the first magnetar to show a surrounding wind nebula. Further, our results imply that such nebulae are no longer exclusive to rotation-powered pulsars and narrow the gap between these two sub-populations of isolated neutron stars. The size and spectrum of the nebula are compatible with those of pulsar-wind nebulae but its radiative efficiency $\eta_{\rm X}=L_{\rm X}/\dot{E}\approx0.1$ is markedly high, possibly pointing to an additional wind component in Swift J1834.9-0846.

Gamma-ray bursts (GRBs) show a dramatic pulse structure that requires bulk relativistic motion, but whose physical origin has remained murky. We focus on a hot, magnetized jet that is emitted by a black hole and interacts with a confining medium. Strongly relativistic expansion of the magnetic field, as limited by a corrugation instability, may commence only after it forms a thin shell. Then the observed $T_{90}$ burst duration is dominated by the curvature delay, and null periods arise from angular inhomogeneities, not the duty cycle of the engine. We associate the $O(1)$ s timescale observed in the pulse width distribution of long GRBs with the collapse of the central 2.5-3$M_\odot$ of a massive stellar core. A fraction of the baryons are shown to be embedded in the magnetized outflow by the hyper-Eddington radiation flux; they strongly disturb the magnetic field after the compactness drops below $\sim 4\times 10^3(Y_e/0.5)^{-1}$. The high-energy photons so created have a compressed pulse structure. Delayed breakout of magnetic field from heavier baryon shells is also a promising approach to X-ray flares. In the second part of the paper, we calculate the imprint of an expanding, scattering photosphere on pulse evolution. Two models for generating the high-energy spectral tail are contrasted: i) pair breakdown due to reheating of an optically thin pair plasma embedded in a thermal radiation field; and ii) continuous heating extending from large to small scattering depth. The second model is strongly inconsistent with the observed hard-to-soft evolution in GRB pulses. The first shows some quantitative differences if the emission is purely spherical, but we show that finite shell width, mild departures from spherical curvature, and latitudinal Lorentz factor gradients have interesting effects.

We consider delayed, volumetric heating in a magnetized outflow that has broken out of a confining medium and expanded to a high Lorentz factor ($\Gamma \sim 10^2-10^3$) and low optical depth to scattering ($\tau_{\rm T} \sim 10^{-3}-10^{-2}$). The energy flux at breakout is dominated by the magnetic field, with a modest contribution from quasi-thermal gamma rays whose spectrum was calculated in Paper I. We focus on the case of extreme baryon depletion in the magnetized material, but allow for a separate baryonic component that is entrained from a confining medium. Dissipation is driven by relativistic motion between these two components, which develops once the photon compactness drops below $ 4\times 10^3(Y_e/0.5)^{-1}$. We first calculate the acceleration of the magnetized component following breakout, showing that embedded MHD turbulence provides significant inertia, the neglect of which leads to unrealistically high estimates of flow Lorentz factor. After re-heating begins, the pair and photon distributions are evolved self-consistently using a one-zone kinetic code that incorporates an exact treatment of Compton scattering, pair production and annihilation, and Coulomb scattering. Heating leads to a surge in pair creation, and the scattering depth saturates at $\tau_{\rm T} \sim$ 1-4. The plasma maintains a very low ratio of particle to magnetic pressure, and can support strong anisotropy in the charged particle distribution, with cooling dominated by Compton scattering. High-energy power-law spectra with photon indices in the range observed in GRBs ($-3 < \beta < -3/2$) are obtained by varying the ratio of heat input to the seed energy in quasi-thermal photons. We contrast our results with those for continuous heating across an expanding photosphere, and show that the latter model produces soft-hard evolution that is inconsistent with observations of GRBs.

The evolution of a dilute electron-positron fireball is calculated in the regime of strong magnetization and very high compactness (l ~10^3-10^8). Heating is applied at a low effective temperature (< 25 keV), and the fireball is allowed to expand, so that the formation of a black-body spectral distribution is inhibited by pair annihilation. The diffusion equation for Compton scattering is coupled to a single-temperature pair gas and an exact (trans-relativistic) cyclo-synchrotron photon source. We find that the photon spectrum develops a quasi-thermal peak, with a power-law slope below it that is characteristic of gamma-ray bursts. The formation of a thermal high-frequency spectrum is checked using the full kinetic equations at l ~ 10^3. These results have several implications for the central engine of GRBs, and the mechanism of energy transport. 1. Baryon rest mass carries less than ~ 10^{-5} of the energy flux at jet breakout inside ~ 10^{12} cm from the engine, with most carried by the magnetic field. 2. This degree of baryon purity points to the presence of an event horizon in the engine, and neutrons play a negligible role in the prompt emission mechanism. 3. X-ray flashes are emitted by outflows carrying enough baryons that the photosphere is pair-depleted, which we show results in faster thermalization. 4. The relation between observed peak frequency and burst luminosity is bounded below by the observed Amati et al. relation if jet Lorentz factor ~ 1/(opening angle) at breakout. 5. Stellar models are used to demonstrate an inconsistency between the highest observed GRB energies, and a hydrodynamic nozzle: magnetic collimation is required. 6. The magnetized pair gas is dilute enough that high-frequency Alfven waves may become charge starved. Finally, we suggest that limitations on magnetic reconnection from plasma collisionality have been overestimated.

The cooling theory of neutron stars is corroborated by its comparison with observations of thermally emitting isolated neutron stars and accreting neutron stars in binary systems. An important ingredient for such an analysis is the age of the object, which, typically, is obtained from the spin-down history. This age is highly uncertain if the object's magnetic field varies appreciably over time. Other age estimators, such as supernova remnant ages and kinematic ages, only apply to few handful of neutron stars. We conduct a population synthesis study of the nearby isolated thermal emitters and obtain their ages statistically from the observed luminosity function of these objects. We argue that a more sensitive blind scan of the galactic disk with the upcoming space telescopes can help to constrain the ages to higher accuracy.

We examine an external trigger mechanism that gives rise to the intense soft gamma-ray repeater (SGR) giant flares. Out of the three giant flares, two showcased the existence of a precursor, which we show to have had initiated the main flare. We develop a reconnection model based on the hypothesis that shearing motion of the footpoints causes the materialization of a Sweet-Parker current layer in the magnetosphere. The thinning of this macroscopic layer due to the development of an embedded super-hot turbulent current layer switches on the impulsive Hall reconnection, which powers the giant flare. We show that the thinning time is on the order of the pre-flare quiescent time.

The light pseudoscalar particle, dubbed the axion, borne out of the Peccei-Quinn solution to the strong CP problem in QCD remains elusive. One avenue of inferring its existence is through its coupling to electromagnetic radiation. So far, laboratory experiments have dedicated all efforts to detect the axion in the mass range $10^{-6} < m_a < 10^{-3}$ eV with a photon-axion coupling strength $g_{a\gamma\gamma} < 10^{-10}GeV^{-1}$, where the limits are derived from astrophysical considerations. In this study, we present a novel way of constraining $g_{a\gamma\gamma}$ by looking at the level of linear polarization in the radiation emerging from magnetic white dwarfs (mWDs). We find that photon-axion oscillations in WD magnetospheres can enhance the degree of linear polarization. Observing that most mWDs show only 5% linear polarization, we derive upper limits on $g_{a\gamma\gamma}$ for different axion masses.

We examine two trigger mechanisms, one internal and the other external to the neutron star, that give rise to the intense soft gamma-ray repeater (SGR) giant flares. So far, three giant flares have been observed from the three out of the seven confirmed SGRs on March 5, 1979, August 27, 1998, and December 27, 2004. The last two events were found to be much more powerful than the first, and both showcased the existence of a precursor, that we show to have had initiated the main flare. In the internal mechanism, we propose that the strongly wound up poloidal magnetic field develops tangential discontinuities and dissipates its torsional energy in heating the crust. The timescale for the instability to develop coincides with the duration of the quiescent state that followed the precursor. Alternatively, we develop a reconnection model based on the hypothesis that shearing motion of the footpoints causes the materialization of a Sweet-Parker current layer in the magnetosphere. The thinning of this macroscopic layer due to the development of an embedded super-hot turbulent current layer switches on the impulsive Hall reconnection, which powers the giant flare. Again, we show that the thinning time is on the order of the preflare quiescent time. This model naturally explains the origin of the observed nonthermal radiation during the flares, as well as the post flare radio afterglows.

We compare the expected abundance of cosmic-ray electrons and positrons from pulsars and magnetars. We assume that the distribution of infant pulsars and magnetars follows that of high-mass stars in the Milky Way and that the production rate of cosmic rays is proportional to the spin-down and magnetic-decay power of pulsars and magnetars, respectively. In combination with primary and secondary cosmic-ray leptons from other sources (especially supernova remnants), we find that both magnetars and pulsars can easily account for the observed cosmic-ray spectrum, in particular the dip seen by HESS at several TeV and the increase in positron fraction found by PAMELA.

A fully relativistic treatment of Bernstein waves in an electron-positron pair plasma has remained too formidable a task owing to the very complex nature of the problem. In this article, we perform contour integration of the dielectric response function and numerically compute the dispersion curves for a uniform, magnetized, relativistic electron-positron pair plasma. The behavior of the dispersion solution for several cases with different plasma temperatures is highlighted. In particular, we find two wave modes that exist only for large wavelengths and frequencies similar to the cyclotron frequency in a moderately relativistic pair plasma. The results presented here have important implications for the study of those objects where a hot magnetized electron-positron plasma plays a fundamental role in generating the observed radiation.

Magnetars, neutron stars with ultra strong magnetic fields ($B\sim 10^{14} - 10^{15}$G), manifest their exotic nature in the form of soft gamma-ray repeaters and anomalous X-ray pulsars. This study estimates the birthrate of magnetars to be $\sim$ 0.22 per century with a galactic population comprising of $\sim$17 objects. A population synthesis was carried out based on the five anomalous X-ray pulsars detected in the ROSAT All-Sky Survey by comparing their number to that of massive OB stars in a well defined volume. Additionally, the group of seven X-ray dim isolated neutron stars detected in the same survey were found to have a birthrate of $\sim$ 2 per century with a galactic population of $\sim$ 22,000 objects.