The fate and observable properties of gamma-ray burst jets depend crucially on their interaction with the progenitor material that surrounds the central engine. We present a semi-analytical model of such interaction, which builds upon several previous analytical and numerical works, aimed at predicting the angular distribution of jet and cocoon energy and Lorentz factor after breakout, given the properties of the ambient material and of the jet at launch. Using this model, we construct synthetic populations of structured jets, assuming either a collapsar (for long gamma-ray bursts -- LGRBs) or a binary neutron star merger (for short gamma-ray bursts -- SGRBs) as progenitor. We assume all progenitors to be identical, and we allow little variability in the jet properties at launch: our populations therefore feature a quasi-universal structure. These populations are able to reproduce the main features of the observed LGRB and SGRB luminosity functions, although several uncertainties and caveats remain to be addressed.

The Gamma Ray Burst (GRB) 180720B is one of the brightest events detected by the Fermi satellite and the first GRB detected by the H.E.S.S. telescope above 100 GeV. We analyse the Fermi (GBM and LAT) and Swift (XRT and BAT) data and describe the evolution of the burst spectral energy distribution in the 0.5 keV - 10 GeV energy range over the first 500 seconds of emission. We reveal a smooth transition from the prompt phase, dominated by synchrotron emission in a moderately fast cooling regime, to the afterglow phase whose emission has been observed from the radio to the GeV energy range. The LAT (0.1 - 100 GeV) light curve initially rises ($F_{\rm LAT}\propto t^{2.4}$), peaks at $\sim$78 s, and falls steeply ($F_{\rm LAT}\propto t^{-2.2}$) afterwards. The peak, which we interpret as the onset of the fireball deceleration, allows us to estimate the bulk Lorentz factor $\Gamma_{0}\sim 150 \ (300)$ under the assumption of a wind-like (homogeneous) circum-burst medium density. We derive a flux upper limit in the LAT energy range at the time of H.E.S.S. detection, but this does not allow us to unveil the nature of the high energy component observed by H.E.S.S. We fit the prompt spectrum with a physical model of synchrotron emission from a non-thermal population of electrons. The 0 - 35 s spectrum after its $E F(E)$ peak (at 1 - 2 MeV) is a steep power law extending to hundreds of MeV. We derive a steep slope of the injected electron energy distribution $N(\gamma)\propto \gamma^{-5}$. Our fit parameters point towards a very low magnetic field ($B'\sim 1 $ G) in the emission region.

We explore the host galaxies of compact-object binaries (black hole--black hole binaries, BHBs; neutron star--black hole binaries, NSBHs; double-neutron stars; DNSs) across cosmic time, by means of population-synthesis simulations combined with the Illustris cosmological simulation. At high redshift ($z\gtrsim{}4$) the host galaxies of BHBs, NSBHs and DNSs are very similar and are predominantly low-mass galaxies (stellar mass $M<10^{11}$ M$_\odot$). If $z\gtrsim{}4$ most compact objects form and merge in the same galaxy, with a short delay time. At low redshift ($z\leq{}2$), the host galaxy populations of DNSs differ significantly from the host galaxies of both BHBs and NSBHs. DNSs merging at low redshift tend to form and merge in the same galaxy, with relatively short delay time. The stellar mass of DNS hosts peaks around $\sim{}10^{10}-10^{11}$ M$_\odot$. In contrast, BHBs and NSBHs merging at low redshift tend to form in rather small galaxies at high redshift and then to merge in larger galaxies with long delay times. This difference between DNSs and black hole binaries is a consequence of their profoundly different metallicity dependence.

We investigate the host galaxies of compact objects merging in the local Universe, by combining the results of binary population-synthesis simulations with the Illustris cosmological box. Double neutron stars (DNSs) merging in the local Universe tend to form in massive galaxies (with stellar mass $>10^{9}$ M$_\odot$) and to merge in the same galaxy where they formed, with a short delay time between the formation of the progenitor stars and the DNS merger. In contrast, double black holes (DBHs) and black hole $-$ neutron star binaries (BHNSs) form preferentially in small galaxies (with stellar mass $<10^{10}$ M$_\odot$) and merge either in small or in larger galaxies, with a long delay time. This result is an effect of metallicity: merging DBHs and BHNSs form preferentially from metal-poor progenitors ($Z\leq{}0.1$ Z$_\odot$), which are more common in high-redshift galaxies and in local dwarf galaxies, whereas merging DNSs are only mildly sensitive to progenitor's metallicity and thus are more abundant in massive galaxies nowadays. The mass range of DNS hosts we predict in this work is consistent with the mass range of short gamma-ray burst hosts.