We present a self-consistent way of modeling synchrotron self-Compton (SSC) effects in gamma-ray burst afterglows, with and without approximated Klein-Nishina suppressed scattering. We provide an analytic approximation of our results, so that it can be incorporated into the afterglow modeling code \texttt{boxfit}, which is currently based on pure synchrotron emission. We discuss the changes in spectral shape and evolution due to SSC effects, and comment on how these changes affect physical parameters derived from broadband modeling. We show that SSC effects can have a profound impact on the shape of the X-ray light curve using simulations including these effects. This leads to data that cannot be simultaneously fit well in both the X-ray and radio bands when considering synchrotron-only fits, and an inability to recover the correct physical parameters, with some fitted parameters deviating orders of magnitude from the simulated input parameters. This may have a significant impact on the physical parameter distributions based on previous broadband modeling efforts.

We present ZTF20aajnksq (AT2020blt), a fast-fading ($\Delta r=2.4$ mag in $\Delta t=1.3$ days) red ($g-r\approx0.6$ mag) and luminous ($M_{1626}=-25.9$) optical transient at $z=2.9$ discovered by the Zwicky Transient Facility (ZTF). AT2020blt shares several features in common with afterglows to long-duration gamma-ray bursts (GRBs): (1) an optical light curve well-described by a broken power-law with a break at $t_\mathrm{j}=1$ day (observer-frame); (2) a luminous $(L_X = 10^{46}$ $\mathrm{erg}$ $\mathrm{s}^{-1})$ X-ray counterpart; and (3) luminous ($L_\nu = 4 \times 10^{31}$ $\mathrm{erg}$ $\mathrm{sec}^{-1}$ $\mathrm{Hz}^{-1}$ at 10 GHz) radio emission. However, no GRB was detected in the 0.74d between the last ZTF non-detection ($r > 20.64$) and the first ZTF detection ($r = 19.57$), with an upper limit on the isotropic-equivalent gamma-ray energy release of $E_{\gamma,\mathrm{iso}} < 7 \times 10^{52}$ erg. AT2020blt is thus the third afterglow-like transient discovered without a detected GRB counterpart (after PTF11agg and ZTF19abvizsw) and the second (after ZTF19abvizsw) with a redshift measurement. We conclude that the properties of AT2020blt are consistent with a classical (initial Lorentz factor $\Gamma_0 \gtrsim 100$) on-axis GRB that was missed by high-energy satellites. Furthermore, by estimating the rate of transients with light curves similar to that of AT2020blt in ZTF high-cadence data, we agree with previous results that there is no evidence for an afterglow-like phenomenon that is significantly more common than classical GRBs. We conclude by discussing the status and future of fast-transient searches in wide-field high-cadence optical surveys.

The recurrent fast radio burst FRB 180916 was recently shown to exhibit a 16 day period (with possible aliasing) in its bursting activity. Given magnetars as widely considered FRB sources, this period has been attributed to precession of the magnetar spin axis or the orbit of a binary companion. Here, we make the simpler connection to a {\it rotational period}, an idea observationally motivated by the 6.7 hour period of the Galactic magnetar candidate, 1E 161348--5055. We explore three physical mechanisms that could lead to the creation of ultra long period magnetars: (i) enhanced spin-down due to episodic mass-loaded charged particle winds (e.g. as may accompany giant flares), (ii) angular momentum kicks from giant flares and (iii) fallback leading to long lasting accretion disks. We show that particle winds and fallback accretion can potentially lead to a sub-set of the magnetar population with ultra long periods, sufficiently long to accommodate FRB 180916 or 1E 161348--5055. If confirmed, such periods implicate magnetars in relatively mature states (ages $1-10$ kyr) and which possessed large internal magnetic fields at birth $B_{\rm int}\gtrsim 10^{16}$ G. In the low-twist magnetar model for FRBs, such long period magnetars may dominate FRB production for repeaters at lower isotropic-equivalent energies and broaden the energy distribution beyond that expected for a canonical population of magnetars which terminate their magnetic activity at shorter periods $P \lesssim 10$ s.

The Fermi-LAT collaboration presented the second gamma-ray burst (GRB) catalog covering its first 10 years of operations. A significant fraction of afterglow-phase light curves in this catalog cannot be explained by the closure relations of the standard synchrotron forward-shock model, suggesting that there could be an important contribution from another process. In view of the above, we derive the synchrotron self-Compton (SSC) light curves from the reverse shock in the thick- and thin-shell regime for a uniform-density medium. We show that this emission could explain the GeV flares exhibited in some LAT light curves. Additionally, we demonstrate that the passage of the forward shock synchrotron cooling break through the LAT band from jets expanding in a uniform-density environment may be responsible for the late time ($\approx10^2$ s) steepening of LAT GRB afterglow light curves. As a particular case, we model the LAT light curve of GRB 160509A that exhibited a GeV flare together with a break in the long-lasting emission, and also two very high energy photons with energies of 51.9 and 41.5 GeV observed 76.5 and 242 s after the onset of the burst, respectively. Constraining the microphysical parameters and the circumburst density from the afterglow observations, we show that the GeV flare is consistent with a SSC reverse-shock model, the break in the long-lasting emission with the passage of the synchrotron cooling break through the Fermi-LAT band and the very energetic photons with SSC emission from the forward shock when the outflow carries a significant magnetic field ($R_{\rm B} \simeq 30$) and it decelerates in a uniform-density medium with a very low density ($n=4.554^{+1.128}_{-1.121}\times 10^{-4}\,{\rm cm^{-3}}$).

We study turbulent gas diffusion affects on $r$-process abundances in Milky Way stars, by a combination of an analytical approach and a Monte Carlo simulation. Higher $r$-process event rates and faster diffusion, lead to more efficient mixing corresponding to a reduced scatter of $r$-process abundances and causing $r$-process enriched stars to start appearing at lower metallicities. We use three independent observations to constrain the model parameters: (i) the scatter of radioactively stable $r$-process element abundances, (ii) the largest $r$-process enrichment values observed in any solar neighborhood stars and (iii) the isotope abundance ratios of different radioactive $r$-process elements ($^{244}$Pu/$^{238}$U and $^{247}$Cm/$^{238}$U) at the early solar system as compared to their formation. Our results indicate that the Galactic $r$-process rate and the diffusion coefficient are respectively $r<4\times 10^{-5}\mbox{ yr}^{-1}, D>0.1 \mbox{ kpc}^2\mbox{Gyr}^{-1}$ ($r<4\times 10^{-6}\mbox{ yr}^{-1}, D>0.5 \mbox{ kpc}^2\mbox{Gyr}^{-1}$ for collapsars or similarly prolific $r$-process sources) with allowed values satisfying an approximate anti-correlation such that $D\approx r^{-2/3}$, implying that the time between two $r$-process events that enrich the same location in the Galaxy, is $\tau_{\rm mix}\approx 100-200\mbox{ Myr}$. This suggests that a fraction of $\sim 0.8$ ($\sim 0.5$) of the observed $^{247}$Cm ($^{244}$Pu) abundance is dominated by one $r$-process event in the early solar system. Radioactively stable element abundances are dominated by contributions from $\sim 10$ different events in the early solar system. For metal poor stars (with [Fe/H]$\lesssim -2$), their $r$-process abundances are dominated by either a single or several events, depending on the star formation history.

The central engine in long gamma-ray bursts (GRBs) is thought to be a compact object produced by the core collapse of massive stars, but its exact nature (black hole or millisecond magnetar) is still debatable. Although the central engine of GRB collapsars is hidden to direct observation, its properties may be imprinted on the accompanying electromagnetic signals. We aim to decipher the generic properties of central engines that are consistent with prompt observations of long GRBs detected by the Burst Alert Telescope (BAT) on board the Neil Gehrels Swift Observatory. Adopting a generic model for the central engine, in which the engine power and activity timescale are independent of each other, we perform Monte Carlo simulations of long GRBs produced by jets that successfully breakout from the star. Our simulations consider the dependence of the jet breakout timescale on the engine luminosity and the effects of the detector's flux threshold. The two-dimensional (2D) distribution of simulated detectable bursts in the gamma-ray luminosity versus gamma-ray duration plane is consistent with the observed one for a range of parameter values describing the central engine. The intrinsic 2D distribution of simulated collapsar GRBs peaks at lower gamma-ray luminosities and longer durations than the observed one, a prediction that can be tested in the future with more sensitive detectors. Black-hole accretors, whose power and activity time are set by the large-scale magnetic flux through the progenitor star and stellar structure, respectively, are compatible with the properties of the central engine inferred by our model.

Observational follow-up of well localized short gamma-ray bursts (SGRBs) has left $20-30\%$ of the population without a coincident host galaxy association to deep optical and NIR limits ($\gtrsim 26$ mag). These SGRBs have been classified as observationally hostless due to their lack of strong host associations. It has been argued that these hostless SGRBs could be an indication of the large distances traversed by the binary neutron star system (due to natal kicks) between its formation and its merger (leading to a SGRB). The distances of GRBs from their host galaxies can be indirectly probed by the surrounding circumburst densities. We show that a lower limit on those densities can be obtained from early afterglow lightcurves. We find that $\lesssim16\%$ of short GRBs in our sample took place at densities $\lesssim10^{-4}$ cm$^{-3}$. These densities represent the expected range of values at distances greater than the host galaxy's virial radii. We find that out of the five SGRBs in our sample that have been found to be observationally hostless, none are consistent with having occurred beyond the virial radius of their birth galaxies. This implies one of two scenarios. Either these observationally hostless SGRBs occurred outside of the half-light radius of their host galaxy, but well within the galactic halo, or in host galaxies at moderate to high redshifts ($z\gtrsim 2$) that were missed by follow-up observations.

We present a model for the hydrodynamics of a relativistic jet interacting with the circum-stellar medium (CSM). The shocked CSM and the jet material are assumed to be in an infinitely thin surface, so the original 2D problem is effectively reduced to 1D. From general conservation laws, we derive the equation of motion for each fluid element along this surface, taking into account the deceleration along the surface normal due to newly swept-up mass and lateral expansion due to pressure gradient in the tangential direction. The pressure and energy density of the shocked CSM are given by the jump conditions at the forward shock. The method is implemented with a finite-differencing numerical scheme, along with calculation of synchrotron emission and absorption from shock-accelerated electrons, in a new code $\mathtt{Jedi}$ (for "jet dynamics"). We present a number of test cases, including top-hat jet, power-law structured jet, "boosted fireball" profile, and CSM with density jump at the wind termination shock. Based on the agreement with other analytical and numerical calculations, we conclude that our simplified method provides a good approximation for the hydrodynamics and afterglow emission for a wide variety of jet structures and CSM density profiles. Efficient modeling of the afterglow from e.g., neutron star mergers, will provide important information on the jet energetics, CSM properties, and the viewing angle.

The joint detection of GW~170817/GRB 170817 confirmed the long-standing theory that binary neutron star mergers produce short gamma-ray burst (sGRB) jets that can successfully break out of the surrounding ejecta. At the same time, the association with a kilonova provided unprecedented information regarding the physical properties (such as masses and velocities) of the different ejecta constituents. Combining this knowledge with the observed luminosities and durations of cosmological sGRBs detected by the Burst Alert Telescope (BAT) onboard the Neil Gehrels Swift Observatory, we revisit the breakout conditions of sGRB jets. Assuming self-collimation of sGRB jets does not play a critical role, we find that the time interval between the binary merger and the launching of a typical sGRB jet is $\lesssim0.1$~s. We also show that for a fraction of at least $\sim 30\%$ of sGRBs, the usually adopted assumption of static ejecta is inconsistent with observations, even if the polar ejecta mass is an order of magnitude smaller than the one in GRB 170817. Our results disfavour magnetar central engines for powering cosmological sGRBs, limit the amount of energy deposited in the cocoon prior to breakout, and suggest that the observed delay of $\sim 1.$7~s in GW 170817 /GRB 170817 between the gravitational wave and $\gamma$-ray signals is likely dominated by the propagation time of the jet to the $\gamma$-ray production site.

A luminous radio burst was recently detected in temporal coincidence with a hard X-ray flare from the Galactic magnetar SGR 1935+2154 with a time and frequency structure consistent with cosmological fast radio bursts (FRB) and a fluence within a factor of $\lesssim 10$ of the least energetic extragalactic FRB previously detected. Although active magnetars are commonly invoked FRB sources, several distinct mechanisms have been proposed for generating the radio emission which make different predictions for the accompanying higher frequency radiation. We show that the properties of the coincident radio and X-ray flares from SGR 1935+2154, including their approximate simultaneity and relative fluence $E_{\rm radio}/E_{\rm X} \sim 10^{-5}$, as well as the duration and spectrum of the X-ray emission, are consistent with extant predictions for the synchrotron maser shock model. Rather than arising from the inner magnetosphere, the X-rays are generated by (incoherent) synchrotron radiation from thermal electrons heated at the same shocks which produce the coherent maser emission. Although the rate of SGR 1935+2154-like bursts in the local universe is not sufficient to contribute appreciably to the extragalactic FRB rate, the inclusion of an additional population of more active magnetars with stronger magnetic fields than the Galactic population can explain both the FRB rate as well as the repeating fraction, however only if the population of active magnetars are born at a rate that is at least two-orders of magnitude lower than that of SGR 1935+2154-like magnetars. This may imply that the more active magnetar sources are not younger magnetars formed in a similar way to the Milky Way population (e.g. via ordinary supernovae), but instead through more exotic channels such as superluminous supernovae, accretion-induced collapse or neutron star mergers.

Very-high-energy (VHE) emission is usually interpreted in the synchrotron-self Compton (SSC) scenario, and expected from the low-redshift and high-luminosity gamma-ray bursts (GRBs), as GRB 180720B and GRB 190114C. Recently, VHE emission was detected by the H.E.S.S. telescopes from one of the closest burst GRB 190829A which was associated with the supernova (SN) 2019oyw. In this paper, we present a temporal and spectral analysis from optical bands to Fermi-LAT energy range over multiple observational periods beginning just after the BAT trigger time and extending for almost three months. We show that the X-ray and optical observations are consistent with synchrotron forward-shock emission evolving between the characteristic and cooling spectral breaks during the early and late afterglow in a uniform-density medium. Modeling the light curves together with its spectral energy distribution, it is shown that the outflow expands with an initial bulk Lorentz factor of $\Gamma\sim 30$, which is high for a low-luminosity GRBs and low for a high-luminosity GRBs. The values of the initial bulk Lorentz factor and the isotropic equivalent energy suggest that GRB 190829A is classified as an intermediate-luminosity burst and consequently, it becomes the first burst of this class in being detected in the VHE gamma-ray band by an imaging atmospheric Cherenkov telescope, and, in turn, the first event without being simultaneously observed by the Fermi-LAT instrument. Analyzing the intermediate-luminosity bursts with $z\lesssim 0.2$ such as GRB 130702A, we show that bursts with intermediate luminosity are potential candidates to be detected in very-high energies.

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 present a systematic search for short-duration gamma-ray bursts (GRBs) in the local Universe based on 14 years of observations with the Neil Gehrels Swift Observatory. We cross-correlate the GRB positions with the GLADE catalogue of nearby galaxies, and find no event at a distance $\lesssim$100 Mpc and four plausible candidates in the range 100 Mpc$\lesssim$$D$$\lesssim$200 Mpc. Although affected by low statistics, this number is higher than the one expected for chance alignments to random galaxies, and possibly suggests a physical association between these bursts and nearby galaxies. By assuming a local origin, we use these events to constrain the range of properties for X-ray counterparts of neutron star mergers. Optical upper limits place tight constraints on the onset of a blue kilonova, and imply either low masses ($\lesssim10^{-3}\,M_{\odot}$) of lanthanide-poor ejecta or unfavorable orientations ($\theta_{obs}\gtrsim$30 deg). Finally, we derive that the all-sky rate of detectable short GRBs within 200 Mpc is $1.3^{+1.7}_{-0.8}$ yr$^{-1}$ (68\% confidence interval), and discuss the implications for the GRB outflow structure. If these candidates are instead of cosmological origin, we set a upper limit of $\lesssim$2.0 yr$^{-1}$ (90\% confidence interval) to the rate of nearby events detectable with operating gamma-ray observatories, such as Swift and Fermi.

Using multiple observational arguments, recent work has shown that cosmological GRBs are typically viewed at angles within, or close to the cores of their relativistic jets. One of those arguments relied on the lack of tens-of-days-long periods of very shallow evolution that would be seen in the afterglow light-curves of GRBs viewed at large angles. Motivated by these results, we consider that GRBs efficiently produce $\gamma$-rays only within a narrow region around the core. We show that, on these near-core lines-of-sight, structured jets naturally produce shallow phases in the X-ray afterglow of GRBs. These plateaus would be seen by a large fraction of observers and would last between $10^2-10^5$ s. They naturally reproduce the observed distributions of time-scales and luminosities as well as the inter-correlations between plateau duration, plateau luminosity and prompt $\gamma$-ray energy. An advantage of this interpretation is that it involves no late time energy injection which would be both challenging from the point of view of the central engine and, as we show here, less natural given the observed correlations between plateau and prompt properties.

We explore the burst energy distribution of fast radio bursts (FRBs) in the low-twist magnetar model of Wadiasingh and Timokhin (2019). Motivated by the power-law fluence distributions of FRB 121102, we propose an elementary model for the FRB luminosity function of individual repeaters with an inversion protocol which directly relates the power-law distribution index of magnetar short burst fluences to that for FRBs. The protocol indicates the FRB energy scales virtually linearly with crust/field dislocation amplitude, if magnetar short bursts prevail in the magnetoelastic regime. Charge starvation in the magnetosphere during bursts (required in WT19) for individual repeaters implies the predicted burst fluence distribution is narrow, $\lesssim 3$ decades for yielding strains and oscillation frequencies feasible in magnetar crusts. Requiring magnetic confinement and charge starvation, we obtain a death line for FRBs which segregates magnetars from the normal pulsar population, suggesting only the former will host recurrent FRBs. We convolve the burst energy distribution for individual magnetars to define the distribution of luminosities in evolved magnetar populations. The broken power-law luminosity function's low energy character depends on the population model, while the high energy index traces that of individual repeaters. Independent of the evolved population, the broken power-law isotropic-equivalent energy/luminosity function peaks at $\sim10^{37}-10^{40}$ erg with a low-energy cutoff at $\sim 10^{37}$ erg. Lastly, we consider the local fluence distribution of FRBs, and find that it can constrain the subset of FRB-producing magnetar progenitors. Our model suggests that improvements in sensitivity may reveal flattening of the global FRB fluence distribution and saturation in FRB rates.

The up-to-now only binary neutron star merger gravitational wave event with detected electromagnetic counterparts, GRB170817A, occurred in a rarefied medium with a density smaller than $10^{-3}-10^{-2}~{\rm cm}^{-3}$. As neutron star binaries are imparted kicks upon formation, and due to their long delay-times before merger, such low-density circum-merger media are generally expected. However, there is some indirect evidence for a class of faster-merging binaries, which would coalesce in denser environments. Nonetheless, astronomical data is largely inconclusive on the possibility of these high-density mergers. We describe a method to directly probe this hypothetical population of high-density mergers through multi-messenger observations of binary neutron star merger afterglows, exploiting the sharp sensitivity of these to the circum-merger medium density. This method is based on a yet-to-be-collected sample of merger afterglows. Its constraining power is large even with a small sample of events. We discuss the method's limitations and applicability. In the upcoming era of 3$^{\rm rd}$-generation gravitational wave detectors, this method's potential will be fully realized as it will allow to probe mergers having occurred soon after the peak of cosmic star formation.

GRB 190114C, a long and luminous burst, was detected by several satellites and ground-based telescopes from radio wavelengths to GeV gamma-rays. In the GeV gamma-rays, the Fermi LAT detected 48 photons above 1 GeV during the first hundred seconds after the trigger time, and the MAGIC telescopes observed for more than one thousand seconds very-high-energy (VHE) emission above 300 GeV. Previous analysis of the multi-wavelength observations showed that although these are consistent with the synchrotron forward-shock model that evolves from a stratified stellar-wind to homogeneous ISM-like medium, photons above few GeVs can hardly be interpreted in the synchrotron framework. In the context of the synchrotron forward-shock model, we derive the light curves and spectra of the synchrotron self-Compton (SSC) model in the stratified and homogeneous medium. In particular, we study the evolution of these light curves during the stratified-to-homogeneous afterglow transition. Using the best-fit parameters reported for GRB 190114C we interpret the photons beyond the synchrotron limit in the SSC framework and model its spectral energy distribution. We conclude that low-redshift GRBs described under a favourable set of parameters as found in the early afterglow of GRB 190114C could be detected at hundreds of GeVs, and also afterglow transitions would allow that VHE emission could be observed for longer periods.

MAXI J1621-501 is the first Swift/XRT Deep Galactic Plane Survey transient that was followed up with a multitude of space missions (NuSTAR, Swift, Chandra, NICER, INTEGRAL, and MAXI) and ground-based observatories (Gemini, IRSF, and ATCA). The source was discovered with MAXI on 2017 October 19 as a new, unidentified transient. Further observations with NuSTAR revealed 2 Type I X-ray bursts, identifying MAXI J1621-501 as a Low Mass X-ray Binary (LMXB) with a neutron star primary. Overall, 24 Type I bursts were detected from the source during a 15 month period. At energies below 10 keV, the source spectrum was best fit with three components: an absorbed blackbody with kT = 2.3 keV, a cutoff power law with index $\Gamma{}$ = 0.7, and an emission line centered on 6.3 keV. Timing analysis of the X-ray persistent emission and burst data has not revealed coherent pulsations from the source or an orbital period. We identified, however, a super-orbital period $\sim{}$78 days in the source X-ray light curve. This period agrees very well with the theoretically predicted radiative precession period of $\sim{}$82 days. Thus, MAXI J1621-501 joins a small group of sources characterized with super-orbital periods.

Binary neutron stars (BNS) mergers are prime sites for $r$-process nucleosynthesis. Their rate determines the chemical evolution of heavy elements in the Milky Way. The merger rate of BNS is a convolution of their birth rate and the gravitational radiation spiral-in delay time. Using the observed population of Galactic BNS we show here that the lifetimes of pulsars in observed BNSs are sufficiently short that the ages of BNSs have little to no effect on the observed merger time distribution. We find that at late times ($t\gtrsim 1$ Gyr) the gravitational wave delay time distribution (DTD) follows the expected $ t^{-1}$. However, a significant excess of rapidly merging systems (between $40-60\%$ of the entire population) is apparent at shorter times. Although the exact shape of the DTD cannot be determined with the existing data, in all models that adequately describe the data we find at least $40\%$ of BNSs with merger times less than 1Gyr. This population of rapid mergers implies a declining deposition rate of $r$-process materials that is consistent with several independent observations of heavy element abundances in the Milky Way. At the same time this population that requires initial binary separations of roughly one solar radius clearly indicates that these binaries had common envelope progenitors. Our results suggest that a significant fraction of future LIGO/Virgo BNS mergers would reside in star forming galaxies.

We have performed a series of Chandra X-ray Observatory observations of the gamma-ray binary LS 2883, which is comprised of a young pulsar (PSR B1259-63) orbiting a massive Be star with a period of 1236.7 days. The system was observed in 5 epochs, spanning a range from 352 to 1175 days after the periastron passage on 2014 May 4. The observations confirmed the recurrent nature of the high-speed ejecta that appear as an extended X-ray structure (clump) moving away from the binary. Compared to the results of the previous monitoring campaign (between the 2010 and 2014 periastron passages), this time we find evidence suggesting that the clump is accelerated to a projected velocity $ v_{\perp}\approx0.15c$ with an acceleration $a_{\perp}=47\pm2$ cm s$^{-2}$ (for uniformly accelerated motion), assuming that it was launched near periastron passage. The observed X-ray properties of the clump are consistent with synchrotron emission from pulsar wind particles accelerated at the interface between the pulsar wind and the clump. We have also performed contemporaneous observations with the {\sl Hubble Space Telescope}, which are used to set an upper limit on the optical flux of the extended emission.

We constrain the formation rate of Galactic magnetars directly from observations. Combining spin-down rates, magnetic activity, and association with supernova remnants, we put a 2$\sigma$ limit on their Galactic formation rate at $2.3-20\mbox{kyr}^{-1}$. This leads to a fraction $0.4_{-0.28}^{+0.6}$ of neutron stars being born as magnetars. We study evolutionary channels that can account for this rate as well as for the periods, period derivatives and luminosities of the observed population. We find that their typical magnetic fields at birth are $3\times 10^{14}-10^{15}$G, and that those decay on a time-scale of $\sim 10^4$years, implying a maximal magnetar period of $P_{\rm max}\approx 13$s. A sizable fraction of the magnetars' energy is released in outbursts. Giant Flares with $E\geq 10^{46}$ erg are expected to occur in the Galaxy at a rate of $\sim 5\mbox{kyr}^{-1}$. Outside our Galaxy, such flares remain observable by {\it Swift} up to a distance of $\sim 100$~Mpc, implying a detection rate of $\sim 5\mbox{ yr}^{-1}$. The specific form of magnetic energy decay is shown to be strongly tied to the total number of observable magnetars in the Galaxy. A systematic survey searching for magnetars could determine the former and inform physical models of magnetic field decay.

The field of time-domain astrophysics has entered the era of Multi-messenger Astronomy (MMA). One key science goal for the next decade (and beyond) will be to characterize gravitational wave (GW) and neutrino sources using the next generation of Extremely Large Telescopes (ELTs). These studies will have a broad impact across astrophysics, informing our knowledge of the production and enrichment history of the heaviest chemical elements, constrain the dense matter equation of state, provide independent constraints on cosmology, increase our understanding of particle acceleration in shocks and jets, and study the lives of black holes in the universe. Future GW detectors will greatly improve their sensitivity during the coming decade, as will near-infrared telescopes capable of independently finding kilonovae from neutron star mergers. However, the electromagnetic counterparts to high-frequency (LIGO/Virgo band) GW sources will be distant and faint and thus demand ELT capabilities for characterization. ELTs will be important and necessary contributors to an advanced and complete multi-messenger network.

The neutron star merger, GW170817, was followed by an optical-infrared transient (a kilonova) which indicated that a substantial ejection of mass at trans-relativistic velocities occurred during the merger. Modeling of the kilonova is able to constrain the kinetic energy of the ejecta and its characteristic velocity, but not the high-velocity distribution of the ejecta. Yet, this distribution contains crucial information on the merger dynamics. In this work, we assume a power-law distribution of the form $E(>\beta\Gamma)\propto(\beta\Gamma)^{-\alpha}$ for the energy of the kilonova ejecta and calculate the non-thermal signatures produced by the interaction of the ejecta with the ambient gas. We find that ejecta with minimum velocity $\beta_0\simeq 0.3$ and energy $E\sim 10^{51}$ erg, as inferred from kilonova modeling, has a detectable radio, and possibly X-ray, afterglow for a broad range of parameter space. This afterglow component is expected to dominate the observed emission on a timescale of a few years post merger and peak around a decade later. Its light curve can be used to determine properties of the kilonova ejecta and in particular the ejecta velocity distribution $\alpha$, the minimum velocity $\beta_0$ and its total kinetic energy $E$. We also predict that an afterglow rebrightening, that is associated with the kilonova component, will be accompanied by a shift of the centroid of the radio source towards the initial position of the explosion.

The joint detection of gravitational waves (GWs) and $\gamma$-rays from a binary neutron-star (NS) merger provided a unique view of off-axis GRBs and an independent measurement of the NS merger rate. Comparing the observations of GRB170817 with those of the regular population of short GRBs (sGRBs), we show that an order unity fraction of NS mergers result in sGRB jets that breakout of the surrounding ejecta. We argue that the luminosity function of sGRBs, peaking at $\approx 2\times 10^{52}\, \mbox{erg s}^{-1}$, is likely an intrinsic property of the sGRB central engine and that sGRB jets are typically narrow with opening angles $\theta_0 \approx 0.1$. We perform Monte Carlo simulations to examine models for the structure and efficiency of the prompt emission in off axis sGRBs. We find that only a small fraction ($\sim 0.01-0.1$) of NS mergers detectable by LIGO/VIRGO in GWs is expected to be also detected in prompt $\gamma$-rays and that GW170817-like events are very rare. For a NS merger rate of $\sim 1500$ Gpc$^{-3}$ yr$^{-1}$, as inferred from GW170817, we expect within the next decade up to $\sim 12$ joint detections with off-axis GRBs for structured jet models and just $\sim 1$ for quasi-spherical cocoon models where $\gamma$-rays are the result of shock breakout. Given several joint detections and the rates of their discoveries, the different structure models can be distinguished. In addition the existence of a cocoon with a reservoir of thermal energy may be observed directly in the UV, given a sufficiently rapid localisation of the GW source.

Motivated by GW170817 we examine constraints that observations put on the angular structure of long gamma-ray burst (GRB) jets. First, the relatively narrow observed distribution of $E_{X}/E_{\gamma}$ (the isotropic equivalent early X-ray afterglow to prompt $\gamma$-ray energy ratio) implies that at any angle that $\gamma$-rays are emitted the Lorentz factor must be high. Specifically, the Lorentz factor of $\gamma$-ray emitting material cannot drop rapidly with angle, and must be $\Gamma(\theta)\gtrsim 50$ even if there are angles for which the gamma-ray received energy is lower by three orders of magnitude compared to the jet core. Second, jets with an angular structure of the $\gamma$-ray emission that over-produce events with a $\gamma$-ray luminosity below the peak of the observed luminosity function are ruled-out. This eliminates models in which the $\gamma$-ray energy angular distribution isn't sufficiently steep and the Lorentz factor distribution isn't sufficiently shallow. Finally, models with a steep structure (e.g. Gaussian) which are detected away from the jet core generate afterglow light-curves that were never observed. We conclude that even if the jet kinetic energy distribution drops continuously with latitude, efficient $\gamma$-ray emission seems to be restricted to material with $\Gamma\gtrsim 50$ and is most likely confined to a narrow region around the core. While our study is confined to long GRBs, where the observed sample is larger and more complete, there are indications that similar conclusions may be applicable also to short GRBs. We discuss the possible implications to the $\gamma$-rays observed in GRB 170817A.

We develop stellar population models to predict the number of binaries with a single luminous member is Gaia and Hipparcos. Our models yield dozens of detections of black hole - luminous companion binaries (BHLC) and hundreds to thousands of neutron star - luminous companion binaries (NSLC) with Gaia. Interestingly, our models also yield a single detection of BHLC binary with Hipparcos, and a few NSLC binaries. We also show how the statistical distribution of detected binaries with a single luminous companion can be used to constrain the formation process of neutron stars and black holes.

The synthesis of $r$-process elements is known to involve extremely energetic explosions. At the same time, recent observations find significant $r$-process enrichment even in extremely small ultra-faint dwarf (UFD) galaxies. This raises the question of retainment of those elements within their hosts. We estimate the retainment fraction and find that it is large $\sim 0.9$, unless the $r$-process event is very energetic ($\gtrsim 10^{52}$erg) and / or the host has lost a large fraction of its gas prior to the event. We estimate the $r$-process mass per event and rate as implied by abundances in UFDs, taking into account imperfect retainment and different models of UFD evolution. The results are consistent with previous estimates (Beniamini et al. 2016) and with the constraints from the recently detected macronova accompanying a neutron star merger (GW170817). We also estimate the distribution of abundances predicted by these models. We find that $\sim 0.07$ of UFDs should have $r$-process enrichment. The results are consistent with both the mean values and the fluctuations of [Eu/Fe] in galactic metal poor stars, supporting the possibility that UFDs are the main 'building blocks' of the galactic halo population.

The young energetic pulsar J1747-2958 ($\tau=26$ kyr, $\dot{E}=2.5\times 10^{36}$ erg s$^{-1}$) powers the Mouse pulsar wind nebula (PWN), famous for its spectacular tail spanning 45$"$ in X-rays and 12$'$ in radio ($d\sim5$ kpc). We present the results of Chandra observations of the PWN and the analysis of archival lower-frequency data. The Chandra HRC image reveals a point-like source at the pulsar position, $\approx$1$'$ behind the bow shock apex of the PWN. The flattened appearance of the compact nebula is consistent with an equatorial outflow deformed by the ram pressure, implying that the angle between the pulsar's spin axis and line of sight is $\sim70^\circ$ (in agreement with the radio and $\gamma$-ray pulse profiles). The spatially-resolved spectroscopy with Chandra ACIS shows that the power-law (PL) spectrum steepens from $\Gamma=1.65\pm0.02$ to $3.0\pm0.1$ over the 45$'$ extent of the X-ray tail. However, the tail's X-ray spectrum integrated over its 45$"$ length fits a single absorbed PL with $\Gamma=2.09\pm0.03$. We also found the Mouse PWN in 150 MHz GMRT data, and a possible counterpart in 24 ${\rm \mu}$m Spitzer data. The multiwavelength data suggest that, at low frequencies, the spectrum of the X-ray-emitting tail region can be described by a broken PL with at least one turnover between radio and X-rays. This is consistent with synchrotron cooling of electrons injected at the termination shock (with an SED slope of 2.2) in an equipartition magnetic field $B\sim200$ ${\rm \mu G}$ and a bulk flow speed $v\sim4000$ km s$^{-1}$.

The recent discovery of gravitational waves from GW170817, associated with a short Gamma-Ray Burst (GRB) at a distance of $40$Mpc, has demonstrated that short GRBs can occur locally and at a reasonable rate. Furthermore, gravitational waves enable us to detect close by GRBs, even when we are observing at latitudes far from the jet's axis. We consider here Compton echoes, the scattered light from the prompt and afterglow emission. Compton echoes, an as yet undetected counterpart of GRBs, peak in X-rays and maintain a roughly constant flux for hundreds to thousands of years after the burst. Though too faint to be detected in typical cosmological GRBs, a fraction of close by bursts with a sufficiently large energy output in X-rays, and for which the surrounding medium is sufficiently dense, may indeed be observed in this way. The detection of a Compton echo could provide unique insight into the burst properties and the environment's density structure. In particular, it could potentially determine whether or not there was a successful jet that broke through the compact binary merger ejecta. We discuss here the properties and expectations from Compton echoes and suggest methods for detectability.

Rapidly spinning, strongly magnetized proto-neutron stars ("millisecond proto-magnetars") are candidate central engines of long-duration gamma-ray bursts (GRB), superluminous supernovae (SLSNe), and binary neutron star mergers. Magnetar birth may be accompanied by the fall-back of stellar debris, lasting for seconds or longer following the explosion. Accretion alters the magnetar evolution by (1) providing an additional source of rotational energy (or a potential sink, if the propeller mechanism operates); (2) enhancing the spin-down luminosity above the dipole rate by compressing the magnetosphere and expanding the polar cap region of open magnetic field lines; (3) supplying an additional accretion-powered neutrino luminosity that sustains the wind baryon-loading, even after the magnetar's internal neutrino luminosity has subsided. The more complex evolution of the jet power and magnetization of an accreting magnetar more readily accounts for the high 56Ni yields GRB SNe and irregular time evolution of some GRB light curves (e.g.~bursts with precursors followed by a long quiescent interval before the main emission episode). Additional baryon-loading from accretion-powered neutrino irradiation of the magnetar polar cap lengthens the timeframe over which the jet magnetization is in the requisite range sigma ~< 1e3 for efficient gamma-ray emission, thereby accommodating GRBs with ultra-long durations. Though accretion does not significantly raise the maximum energy budget from the limit of <~ few 1e52 ergs for an isolated magnetar, it greatly expands the range of magnetic field strengths and birth spin periods capable of powering GRB jets, reducing the differences between the magnetar properties normally invoked to explain GRBs versus SLSNe.

Previous studies have considered synchrotron as the emission mechanism for prompt Gamma-Ray Bursts (GRBs). These works have shown that the electrons must cool on a timescale comparable to the dynamic time at the source in order to satisfy spectral constraints while maintaining high radiative efficiency. We focus on conditions where synchrotron cooling is balanced by a continuous source of heating, and in which these constraints are naturally satisfied. Assuming that a majority of the electrons in the emitting region are contributing to the observed peak, we find that the energy per electron has to be $E\gtrsim 20$ GeV and that the Lorentz factor of the emitting material has to be very large $10^3\lesssim \Gamma_{\rm em} \lesssim 10^4$, well in excess of the bulk Lorentz factor of the jet inferred from GRB afterglows. A number of independent constraints then indicate that the emitters must be moving relativistically, with $\Gamma'\approx 10$, relative to the bulk frame of the jet and that the jet must be highly magnetized upstream of the emission region, $\sigma_{\rm up}\gtrsim 30$. The emission radius is also strongly constrained in this model to $R\gtrsim 10^{16}$cm. These values are consistent with magnetic jet models where the dissipation is driven by magnetic reconnection that takes place far away from the base of the jet.

Neutron star mergers have been long considered as promising sites of heavy $r$-process nucleosynthesis. We overview observational evidence supporting this scenario including: the total amount of $r$-process elements in the Galaxy, extreme metal poor stars, geological radioactive elemental abundances, dwarf galaxies, and short gamma-ray bursts (sGRBs). Recently, the advanced LIGO and Virgo observatories discovered a gravitational-wave signal of a neutron star merger, GW170817, as well as accompanying multi-wavelength electromagnetic (EM) counterparts. The ultra-violet, optical, and near infrared observations point to $r$-process elements that have been synthesized in the merger ejecta. The rate and ejected mass inferred from GW170817 and the EM counterparts are consistent with other observations. We find however that, within simple one zone chemical evolution models (based on merger rates with reasonable delay time distributions as expected from evolutionary models, or from observations of sGRBs), it is difficult to reconcile the current observations of the europium abundance history of Galactic stars for [Fe/H] $\gtrsim -1$. This implies that to account for the role of mergers in the Galactic chemical evolution, we need a Galactic model with multiple populations that have different spatial distributions and/or varying formation rates.

Gamma-ray burst (GRB) afterglows have been observed across the electromagnetic spectrum, and physical parameters of GRB jets and their surroundings have been derived using broadband modeling. While well-sampled lightcurves across the broadband spectrum are necessary to constrain all the physical parameters, some can be strongly constrained by the right combination of just a few observables, almost independently of the other unknowns. We present a method involving the peaks of radio lightcurves to constrain the fraction of shock energy that resides in electrons, $\epsilon_e$. This parameter is an important ingredient for understanding the microphysics of relativistic shocks; Based on a sample of 36 radio afterglows, we find $\epsilon_e$ has a narrow distribution centered around $0.13-0.15$. Our method is suggested as a diagnostic tool for determining $\epsilon_e$, and to help constrain the broadband modeling of GRB afterglows. Some earlier measurements of the spreads in parameter values for $\epsilon_e$, the kinetic energy of the shock, and the density of the circumburst medium, based on broadband modeling across the entire spectrum, are at odds with our analysis of radio peaks. This could be due to different modeling methods and assumptions, and possibly missing ingredients in past and current modeling efforts. Furthermore, we show that observations at $\gtrsim10$~GHz performed $0.3-30$ days after the GRB trigger, are best suited for pinpointing the synchrotron peak frequency, and consequently $\epsilon_e$. At the same time, observations at lower radio frequencies can pin down the synchrotron self-absorption frequency and help constrain the other physical parameters of GRB afterglows.

We examine millisecond magnetars as central engines of Gamma Ray Bursts' (GRB) prompt emission. Using the proto-magnetar wind model of Metzger et al. 2011, we estimate the temporal evolution of the magnetization and power injection at the base of the GRB jet and apply these to different prompt emission models to make predictions for the GRB energetics, spectra and lightcurves. We investigate both shock and magnetic reconnection models for the particle acceleration, as well as the effects of energy dissipation across optically thick and thin regions of the jet. The magnetization at the base of the jet, $\sigma_0$, is the main parameter driving the GRB evolution in the magnetar model and the emission is typically released for $100\lesssim \sigma_0 \lesssim 3000$. Given the rapid increase in $\sigma_0$ as the proto-magnetar cools and its neutrino-driven mass loss subsides, the GRB duration is typically limited to $\lesssim 100$ s. This low baryon loading at late times challenges magnetar models for ultra-long GRBs, though black hole models likely run into similar difficulties without substantial entrainment from the jet walls. The maximum radiated gamma-ray energy is $\lesssim 5 \times 10^{51}$erg, significantly less than the magnetar's total initial rotational energy and in strong tension with the high end of the observed GRB energy distribution. However, the gradual magnetic dissipation model (Beniamini & Giannios 2017) applied to a magnetar central engine, naturally explains several key observables of typical GRBs, including energetics, durations, stable peak energies, spectral slopes and a hard to soft evolution during the burst.

It has been proposed that an intermediate mass black hole can produce compact molecular clouds with large velocity dispersions. This model is called into question due to the discovery of several such short lived clouds, which suggests they are common, and hence it is unlikely that their formation involves an exotic object such as an intermediate mass black hole. In this paper we propose an alternative model, where supernovae produce such clouds. We apply our model to a such a cloud called CO-0.40-0.22. A previous work pointed out the fact that the kinetic energy of this cloud is comparable to that from a supernova, but disqualified the supernova scenario due to lack of thermal emission. We explain the lack of thermal emission as a consequence of the large density of the medium. In such a medium, the supernova shock is radiation dominated, and hence the initial temperatures are lower by many orders of magnitude compared with a typical supernova. Finally, we discuss some of the consequences of having a significant population of intermediate mass black holes in the galactic centre.

Context: The origin of GRBs' prompt emission is highly debated. Proposed scenarios involve dissipation processes above or below the photosphere of an ultra-relativistic outflow. Aims: We search for observational features that would favour one scenario over the others by constraining the dissipation radius, the outflow magnetization or by indicating the presence of shocks. Bursts showing peculiarities can emphasize the role of a specific physical ingredient, which becomes more apparent under certain circumstances. Methods: We study GRB 070110, which exhibited several remarkable features during its early afterglow: a very flat plateau terminated by an extremely steep drop and immediately followed by a bump. We model the plateau as photospheric emission from a long lasting outflow of moderate Lorentz factor ($\Gamma\sim 20$) which lags behind an ultra-relativistic ($\Gamma> 100$) ejecta responsible for the prompt emission. We compute the dissipation of energy in the forward and reverse shocks resulting from this ejecta's deceleration by the external medium. Results: Photospheric emission from the long-lasting outflow can account for the plateau properties (luminosity and spectrum) assuming some dissipation takes place in the flow. The geometrical timescale at the photospheric radius is so short that the observed decline at the end of the plateau likely corresponds to the shut-down of the central engine. The following bump results from dissipated power in the reverse shock, which develops when the slower material catches up with the initially fast component, after it had been decelerated. Conclusions: Our interpretation suggests that the prompt phase resulted from dissipation above the photosphere while the plateau had a photospheric origin. If the bump is produced by the reverse shock, it implies an upper limit ($\sigma \lesssim 0.1$) on the magnetization of the slower material.

We considered a model for the prompt phase of Gamma-Ray Burst (GRB) emission arising from a magnetized jet undergoing gradual energy dissipation due to magnetic reconnection. The dissipated magnetic energy is translated to bulk kinetic energy and to acceleration of particles. The energy in these particles is released via synchrotron radiation as they gyrate around the strong magnetic fields in the jet. At small radii, the optical depth is large, and the radiation is reprocessed through Comptonization into a narrow, strongly peaked, component. At larger distances the optical depth becomes small and radiation escapes the jet with a non-thermal distribution. The obtained spectra typically peak around $\approx 300$keV (as observed) and with spectral indices below and above the peak that are, for a broad range of the model parameters, close to the observed values. The small radius of dissipation causes the emission to become self absorbed at a few keV and can sufficiently suppress the optical and X-ray fluxes within the limits required by observations (Beniamini & Piran 2014).

Recent observations of $r$-process material in ultra-faint dwarf galaxies (UFDs) shed light on the sources of these elements. Strong upper limits on the Eu mass in some UFDs combined with detections of much larger masses in a UFD, Reticulum II, and other dwarf galaxies imply that Eu production is dominated by rare events, and that the minimal Eu mass observed in any UFD is approximately the amount of Eu mass produced per event. This is consistent with other independent observations in the Galaxy. We estimate, using a model independent likelihood analysis, the rate and Eu (Fe) mass produced per $r$-process (Fe production) event in dwarf galaxies including classical dwarfs and UFDs. The mass and rate of the Fe production events are consistent with the normal core-collapse supernova~(ccSN) scenario. The Eu mass per event is $3\times 10^{-5}M_{\odot}<\tilde{m}_{\rm Eu}<2\times 10^{-4}M_{\odot}$, corresponding to a total $r$-process mass per event of $6\times 10^{-3}M_{\odot}<\tilde{m}_{r-process}<4\times 10^{-2}M_{\odot}$. The rate of $r$-process events is $2.5\times 10^{-4}<R_{rp/SN}<1.4\times 10^{-3}$ as compared with the ccSNe rate. These values are consistent with the total Eu mass observed in our own Galaxy, suggesting that the same mechanism is behind the production of $r$-process events in both dwarf galaxies and the Milky Way, and that it may be the dominant mechanism for production of $r$-process elements in the Universe. The results are consistent with neutron star mergers estimates, but cannot rule out other rare core collapse scenarios, provided that they produce a significant amounts of $r$-process material per event.

Merging neutron star binaries are prime candidate sources for heavy r-process nucleosynthesis. The amount of heavy r-process material is consistent with the mass ejection and rates of mergers, and abundances of relic radioactive materials suggest that heavy r-process material is produced in rare events. Observations of possible macronovae provide further support for this model. Still, some concerns remain. One is the observation of heavy r-process elements in Ultra Faint Dwarf (UFD) galaxies. The escape velocities from UFDs are so small that the natal kicks, taking place at neutron stars birth, might eject such binaries from UFDs. Furthermore the old stellar populations of UFDs requires that r-process nucleosynthesis must have taken place very early on, while it may take several Gyr for compact binaries to merge. This last problem arises also within the Milky Way where heavy r-process materials has been observed in some low metallicity stars. We show here that since a significant fraction of neutron star binaries form with a very small proper motion those won't be ejected even from a UFD. Furthermore, a significant fraction of these binaries merge within a few hundred Myr and in some cases the merger time is even shorter than a hundred Myr. This population of "rapid mergers" explains the appearance of heavy r-process material in both UFDs and in the Milky Way.

The prompt Gamma-Ray Bursts' (GRBs) efficiency is an important clue on the emission mechanism producing the $\gamma$-rays. Previous estimates of the kinetic energy of the blast waves, based on the X-ray afterglow luminosity $L_X$, suggested that this efficiency is large, with values above 90\% in some cases. This poses a problem to emission mechanisms and in particular to the internal shocks model. These estimates are based, however, on the assumption that the X-ray emitting electrons are fast cooling and that their Inverse Compton (IC) losses are negligible. The observed correlations between $L_X$ (and hence the blast wave energy) and $E_{\gamma\rm ,iso}$, the isotropic equivalent energy in the prompt emission, has been considered as observational evidence supporting this analysis. It is reasonable that the prompt gamma-ray energy and the blast wave kinetic energy are correlated and the observed correlation corroborates, therefore, the notion $L_X$ is indeed a valid proxy for the latter. Recent findings suggest that the magnetic field in the afterglow shocks is significantly weaker than was earlier thought and its equipartition fraction, $\epsilon_B$, could be as low as $10^{-4}$ or even lower. Motivated by these findings we reconsider the problem, taking now IC cooling into account. We find that the observed $L_X-E_{\gamma\rm ,iso}$ correlation is recovered also when IC losses are significant. For small $\epsilon_B$ values the blast wave must be more energetic and we find that the corresponding prompt efficiency is significantly smaller than previously thought. For example, for $\epsilon_B\sim10^{-4}$ we infer a typical prompt efficiency of $\sim15\%$.

The energy dissipation mechanism within Gamma-Ray Burst (GRB) outflows, driving their extremely luminous prompt $\gamma$-ray emission is still uncertain. The leading candidates are internal shocks and magnetic reconnection. While the emission from internal shocks has been extensively studied, that from reconnection still has few quantitative predictions. We study the expected prompt-GRB emission from magnetic reconnection and compare its temporal and spectral properties to observations. The main difference from internal shocks is that for reconnection one expects relativistic bulk motions with Lorentz factors $\Gamma'\gtrsim\;$a few in the jet's bulk frame. We consider such motions of the emitting material in two anti-parallel directions (e.g. of the reconnecting magnetic-field lines) within an ultra-relativistic (with $\Gamma\gg1$) thin spherical reconnection layer. The emission's relativistic beaming in the jet's frame greatly affects the light-curves. For emission at radii $R_0<R<R_0+\Delta{}R$ (with $\Gamma=\rm{const}$) the observed pulse width is $\Delta{}T\sim(R_0/2c\Gamma^2)\max(1/\Gamma',\,\Delta{}R/R_0)$, i.e. up to $\sim\Gamma'$ times shorter than for isotropic emission in the jet's frame. We consider two possible magnetic reconnection modes: a quasi steady-state with continuous plasma flow into and out of the reconnection layer, and sporadic reconnection in relativistic turbulence that produces relativistic plasmoids. Both of these modes can account for many observed prompt-GRB properties: variability, pulse asymmetry, the very rapid declines at their end, and pulse evolutions that are either hard to soft (for $\Gamma'\lesssim2$) or intensity tracking (for $\Gamma'>2$). However, the relativistic turbulence mode is more likely to be relevant for the prompt sub-MeV emission and can naturally account also for the peak luminosity$\,$--$\,$peak frequency correlation.

Double Neutron Stars (DNS) have to survive two supernovae and still remain bound. This sets strong limits on the nature of the second collapse in these systems. We consider the masses and orbital parameters of the DNS population and constrain the two distributions of mass ejection and kick velocities directly from observations with no a-priori assumptions regarding evolutionary models and/or the types of the supernovae involved. We show that there is strong evidence for two distinct types of supernovae in these systems, where the second collapse in the majority of the observed systems involved small mass ejection ($\Delta M\lesssim 0.5M_{\odot}$) and a corresponding low-kick velocity ($v_{k}\lesssim 30 $km\,s$^{-1}$). This formation scenario is compatible, for example, with an electron capture supernova. Only a minority of the systems have formed via the standard SN scenario involving larger mass ejection of $\sim 2.2 M_{\odot}$ and kick velocities of up to $400$km\,s$^{-1}$. Due to the typically small kicks in most DNS (which are reflected by rather low proper motion), we predict that most of these systems reside close to the galactic disc. In particular, this implies that more NS-NS mergers occur close to the galactic plane. This may have non-trivial implications to the estimated merger rates of DNS and to the rate of LIGO / VIRGO detections.

Observations of X-ray flares from Gamma Ray Bursts (GRBs) imply strong constraints on possible physical models. We provide a general discussion of these. In particular, we show that in order to account for the relatively flat and weak optical flux during the X-ray flares, the size of the emitting region should be $\lesssim 3\times 10^{14}$cm. The bolometric luminosity of flares also strongly constrain the energy budget, and are inconsistent with late time activity of a central engine powered by the spin-down of a magnetar. We provide a simple toy model according to which flares are produced by an outflow of modest Lorentz factor (a few tens instead of hundreds) that is launched more or less simultaneously with the highly relativistic jet which produced the prompt gamma-ray emission. The "slower" moving outflow produces the flare as it reaches its photosphere. The existence of such a component would naturally resolve the observational challenges imposed by flares, outlined in this work.

We consider a sample of ten GRBs with long lasting ($\gtrsim10^2\rm\,sec$) emission detected by Fermi/LAT and for which X-ray data around $1\,$day are also available. We assume that both the X-rays and the GeV emission are produced by electrons accelerated at the external forward shock, and show that the X-ray and the GeV fluxes lead to very different estimates of the initial kinetic energy of the blast wave. The energy estimated from GeV is on average $\sim50$ times larger than the one estimated from X-rays. We model the data (accounting also for optical detections around $1\,$day, if available) to unveil the reason for this discrepancy and find that good modelling within the forward shock model is always possible and leads to two possibilities: either the X-ray emitting electrons (unlike the GeV emitting electrons) are in the slow cooling regime or ii) the X-ray synchrotron flux is strongly suppressed by Compton cooling, whereas, due to the Klein-Nishina suppression, this effect is much smaller at GeV energies. In both cases the X-ray flux is no longer a robust proxy for the blast wave kinetic energy. On average, both cases require weak magnetic fields ($10^{-6}\lesssim \epsilon_B \lesssim 10^{-3}$) and relatively large isotropic kinetic blast wave energies $10^{53}\rm\,erg<E_{0,kin}<10^{55}\rm\,erg$ corresponding to large lower limits on the collimated energies, in the range $10^{52}\rm\,erg<E_{\theta,kin}<5\times10^{52}\rm\,erg$ for an ISM environment with $n\sim 1\mbox{cm}^{-3}$ and $10^{52}\rm\,erg<E_{\theta,kin}<10^{53}\rm\,erg$ for a wind environment with $A_* \sim 1$. These energies are larger than those estimated from the X-ray flux alone, and imply smaller inferred values of the prompt efficiency mechanism, reducing the efficiency requirements on the still uncertain mechanism responsible for prompt emission.

We examine GRBs with both Fermi-LAT and X-ray afterglow data. Assuming that the 100MeV (LAT) emission is radiation from cooled electrons accelerated by external shocks, we show that the kinetic energy of the blast wave estimated from the 100MeV flux is 50 times larger than the one estimated from the X-ray flux. This can be explained if either: i) electrons radiating at X-rays are significantly cooled by SSC (suppressing the synchrotron flux above the cooling frequency) or ii) if the X-ray emitting electrons, unlike those emitting at 100MeV energies, are in the slow cooling regime. In both cases the X-ray flux is no longer an immediate proxy of the blast wave kinetic energy. We model the LAT, X-ray and optical data and show that in general these possibilities are consistent with the data, and explain the apparent disagreement between X-ray and LAT observations. All possible solutions require weak magnetic fields: $10^{-6}< \epsilon_B < 10^{-3}$ (where $\epsilon_B$ is the fraction of shocked plasma energy in magnetic fields). Using the LAT emission as a proxy for the blast wave kinetic energy we find that the derived prompt efficiencies are of order 15%. This is considerably lower compared with previous estimates (87% and higher for the same bursts). This provides at least a partial solution to the "prompt high efficiency paradox".

We consider the conditions within a Gamma-Ray Burst (GRB) emission region that is Poynting flux dominated. Due to the enormous magnetic energy density, relativistic electrons will cool in such a region extremely rapidly via synchrotron. As there is no known mechanism that can compete with synchrotron it must be the source of the prompt sub-MeV emission. This sets strong limits on the size and Lorentz factor of the outflow. Furthermore, synchrotron cooling is too efficient. It overproduces optical and X-ray as compared with the observations. This overproduction of low energy emission can be avoided if the electrons are re-accelerated many times ($\gtrsim 5\times 10^4$)during each pulse (or are continuously heated) or if they escape the emitting region before cooling down. We explore the limitations of both models practically ruling out the later and demonstrating that the former requires two different acceleration mechanisms as well as an extremely large magnetic energy to Baryonic energy ratio. To be viable, any GRB model based on an emission region that is Poynting flux dominated must demonstrate how these conditions are met. We conclude that if GRB jets are launched magnetically dominated they must dissipate somehow most of their magnetic energy before they reach the emission region.

We reexamine the general synchrotron model for GRBs' prompt emission and determine the regime in the parameter phase space in which it is viable. We characterize a typical GRB pulse in terms of its peak energy, peak flux and duration and use the latest Fermi observations to constrain the high energy part of the spectrum. We solve for the intrinsic parameters at the emission region and find the possible parameter phase space for synchrotron emission. Our approach is general and it does not depend on a specific energy dissipation mechanism. Reasonable synchrotron solutions are found with energy ratios of $10^{-4}<\epsilon_B/\epsilon_e<10$, bulk Lorentz factor values of $300<\Gamma<3000$, typical electrons' Lorentz factor values of $3\times 10^3<\gamma_e<10^5$ and emission radii of the order $10^{15}cm<R<10^{17}$cm. Most remarkable among those are the rather large values of the emission radius and the electron's Lorentz factor. We find that soft (with peak energy less than 100KeV) but luminous (isotropic luminosity of $~1.5 \times 10^{53}$) pulses are inefficient. This may explain the lack of strong soft bursts. In cases when most of the energy is carried out by the kinetic energy of the flow, such as in the internal shocks, the synchrotron solution requires that only a small fraction of the electrons are accelerated to relativistic velocities by the shocks. We show that future observations of very high energy photons from GRBs by CTA, could possibly determine all parameters of the synchrotron model or rule it out altogether.

The Large Area Telescope (LAT) on board of the Fermi satellite detected emission above 30 MeV only in a small fraction of the long gamma-ray bursts (GRBs) detected by the Fermi Gamma-ray Burst Monitor (GBM) at 8 keV - 10 MeV. Those bursts that were detected by the LAT were among the brightest GBM bursts. We examine a sample of the most luminous GBM bursts with no LAT detection and obtain upper limits on their high energy fluence. We find an average upper limit of LAT/GBM fluence ratio of 0.13 for GeV fluence during $T_{90}$ and an average upper limit ratio of 0.45 for GeV fluence during the first 600 seconds after the trigger. These ratios strongly constrain various emission models and in particular rule out SSC models for the prompt emission. In about a third of both LAT detected and LAT non-detected bursts, we find that the extrapolation of the MeV range Band spectrum to the GeV range is larger than the observed GeV fluence (or its upper limit). While this excess is not highly significant for any specific burst, the overall excess in a large fraction of the bursts suggests a decline in the high energy spectral slope in at least some of these bursts. Possibly an evidence for the long sought after pair creation limit.