On 2019 August 14, the LIGO and Virgo interferometers detected a high-significance event labelled S190814bv. Preliminary analysis of the GW data suggests that the event was likely due to the merger of a compact binary system formed by a BH and a NS. ElectromagNetic counterparts of GRAvitational wave sources at the VEry Large Telescope (ENGRAVE) collaboration members carried out an intensive multi-epoch, multi-instrument observational campaign to identify the possible optical/near infrared counterpart of the event. In addition, the ATLAS, GOTO, GRAWITA-VST, Pan-STARRS and VINROUGE projects also carried out a search on this event. Our observations allow us to place limits on the presence of any counterpart and discuss the implications for the kilonova (KN) possibly generated by this NS-BH merger, and for the strategy of future searches. Altogether, our observations allow us to exclude a KN with large ejecta mass $M\gtrsim 0.1\,\mathrm{M_\odot}$ to a high ($>90\%$) confidence, and we can exclude much smaller masses in a subsample of our observations. This disfavours the tidal disruption of the neutron star during the merger. Despite the sensitive instruments involved in the campaign, given the distance of S190814bv we could not reach sufficiently deep limits to constrain a KN comparable in luminosity to AT 2017gfo on a large fraction of the localisation probability. This suggests that future (likely common) events at a few hundreds Mpc will be detected only by large facilities with both high sensitivity and large field of view. Galaxy-targeted observations can reach the needed depth over a relevant portion of the localisation probability with a smaller investment of resources, but the number of galaxies to be targeted in order to get a fairly complete coverage is large, even in the case of a localisation as good as that of this event.

We present the observations of Lyman continuum (LyC) emission in the afterglow spectra of GRB 191004B at $z=3.5055$, together with those of the other two previously known LyC-emitting long gamma-ray bursts (LGRB) (GRB 050908 at $z=3.3467$, and GRB 060607A at $z=3.0749$), to determine their LyC escape fraction and compare their properties. From the afterglow spectrum of GRB 191004B we determine a neutral hydrogen column density at the LGRB redshift of $\log(N_{\rm HI}/cm^{-2})= 17.2 \pm 0.15$, and negligible extinction ($A_{\rm V}=0.03 \pm 0.02$ mag). The only metal absorption lines detected are CIV and SiIV. In contrast to GRB 050908 and GRB 060607A, the host galaxy of GRB 191004B displays significant Ly$\alpha$ emission. From its Ly$\alpha$ emission and the non-detection of Balmer emission lines we constrain its star-formation rate (SFR) to $1 \leq$ SFR $\leq 4.7$ M$_{\odot}\ yr^{-1}$. We fit the Ly$\alpha$ emission with a shell model and find parameters values consistent with the observed ones. The absolute LyC escape fractions we find for GRB 191004B, GRB 050908 and GRB 060607A are of $0.35^{+0.10}_{-0.11}$, $0.08^{+0.05}_{-0.04}$ and $0.20^{+0.05}_{-0.05}$, respectively. We compare the LyC escape fraction of LGRBs to the values of other LyC emitters found from the literature, showing that LGRB afterglows can be powerful tools to study LyC escape for faint high-redshift star-forming galaxies. Indeed we could push LyC leakage studies to much higher absolute magnitudes. The host galaxies of the three LGRB presented here have all $M_{\rm 1600} > -19.5$ mag, with the GRB 060607A host at $M_{\rm 1600} > -16$ mag. LGRB hosts may therefore be particularly suitable for exploring the ionizing escape fraction in galaxies that are too faint or distant for conventional techniques. Furthermore the time investment is very small compared to galaxy studies. [Abridged]

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 process of massive star ($M\geq8~M_\odot$) formation is still poorly understood. Observations of massive young stellar objects (MYSOs) are challenging due to their rarity, short formation timescale, large distances, and high circumstellar extinction. Here, we present the results of a spectroscopic analysis of a population of MYSOs in the Large Magellanic Cloud (LMC). We took advantage of the spectral resolution and wavelength coverage of X-shooter (300-2500 nm), mounted on the European Southern Observatory Very Large Telescope, to detect characteristic spectral features in a dozen MYSO candidates near 30 Doradus, the largest starburst region in the Local Group hosting the most massive stars known. The X-shooter spectra are strongly contaminated by nebular emission. We used a scaling method to subtract the nebular contamination from our objects. We detect H$\alpha,\beta$, [O I] 630.0 nm, Ca II infrared triplet, [Fe II] 1643.5 nm, fluorescent Fe II 1687.8 nm, H$_2$ 2121.8 nm, Br$\gamma$, and CO bandhead emission in the spectra of multiple candidates. This leads to the spectroscopic confirmation of 10 candidates as bona fide MYSOs. We compare our observations with photometric observations from the literature and find all MYSOs to have a strong near-infrared excess. We compute lower limits to the brightness and luminosity of the MYSO candidates, confirming the near-infrared excess and the massive nature of the objects. No clear correlation is seen between the Br$\gamma$ luminosity and metallicity. Combining our sample with other LMC samples results in a combined detection rate of disk features such as fluorescent Fe II and CO bandheads which is consistent with the Galactic rate (40\%). Most of our MYSOs show outflow features.

GRB 190114C is the first GRB for which the detection of very-high energy emission up to the TeV range has been reported. It is still unclear whether environmental properties might have contributed to the production of these very high-energy photons, or if it is solely related to the released GRB emission. The relatively low redshift of the GRB (z=0.425) allows us to study the host galaxy of this event in detail, and to potentially identify idiosyncrasies that could point to progenitor characteristics or environmental properties responsible for such a unique event. We use ultraviolet, optical, infrared and submillimetre imaging and spectroscopy obtained with HST, VLT and ALMA to obtain an extensive dataset on which the analysis of the host galaxy is based. The host system is composed of a close pair of interacting galaxies (Delta v = 50 km s^-1), both of which are well-detected by ALMA in CO(3-2). The GRB occurred within the nuclear region (~170 pc from the centre) of the less massive but more star-forming galaxy of the pair. The host is more massive (log(M/M_odot)=9.3) than average GRB hosts at that redshift and the location of the GRB is rather unique. The enhanced star-formation rate was probably triggered by tidal interactions between the two galaxies. Our ALMA observations indicate that both host galaxy and companion have a high molecular gas fraction, as has been observed before in interacting galaxy pairs. The location of the GRB within the core of an interacting galaxy with an extinguished line-of-sight is indicative of a denser environment than typically observed for GRBs and could have been crucial for the generation of the very-high-energy photons that were observed.

Mapping of the large-scale structure through cosmic time has numerous applications in the studies of cosmology and galaxy evolution. At $z > 2$, the structure can be traced by the neutral intergalactic medium (IGM) by way of observing the Ly$\alpha$, forest towards densely-sampled lines-of-sight of bright background sources, such as quasars and star forming galaxies. We investigate the scientific potential of MOSAIC, a planned multi-object spectrograph on the European Extremely Large Telescope (ELT), for the 3D mapping of the IGM at $z \gtrsim 3$. We simulate a survey of $3 \lesssim z \lesssim 4$ galaxies down to a limiting magnitude of $m_{r}\sim 25.5$ mag in an area of 1 degree$^2$ in the sky. Galaxies and their spectra (including the line-of-sight Ly$\alpha$ absorption) are taken from the lightcone extracted from the Horizon-AGN cosmological hydrodynamical simulation. The quality of the reconstruction of the original density field is studied for different spectral resolutions and signal-to-noise ratios of the spectra. We demonstrate that the minimum $S/N$ (per resolution element) of the faintest galaxies that such survey has to reach is $S/N = 4$. We show that a survey with such sensitivity enables a robust extraction of cosmic filaments and the detection of the theoretically-predicted galaxy stellar mass and star-formation rate gradients towards filaments. By simulating the realistic performance of MOSAIC we obtain $S/N(T_{\rm obs}, R, m_{r})$ scaling relations. We estimate that $\lesssim 35~(65)$ nights of observation time are required to carry out the survey with the instrument's high multiplex mode and with the spectral resolution of $R=1000~(2000)$. A survey with a MOSAIC-concept instrument on the ELT is found to enable the mapping of the IGM at $z > 3$ on Mpc scales, and as such will be complementary to and competitive with other planned IGM tomography surveys. [abridged]

The number of supernovae known to be connected with long-duration gamma-ray bursts is increasing and the link between these events is no longer exclusively found at low redshift ($z \lesssim 0.3$) but is well established also at larger distances. We present a new case of such a liaison at $z = 0.33$ between GRB\,171010A and SN\,2017htp. It is the second closest GRB with an associated supernova of only three events detected by Fermi-LAT. The supernova is one of the few higher redshift cases where spectroscopic observations were possible and shows spectral similarities with the well-studied SN\,1998bw, having produced a similar Ni mass ($M_{\rm Ni}=0.33\pm0.02 ~\rm{M_{\odot}}$) with slightly lower ejected mass ($M_{\rm ej}=4.1\pm0.7~\rm{M_{\odot}}$) and kinetic energy ($E_{\rm K} = 8.1\pm2.5 \times 10^{51} ~\rm{erg}$). The host-galaxy is bigger in size than typical GRB host galaxies, but the analysis of the region hosting the GRB revealed spectral properties typically observed in GRB hosts and showed that the progenitor of this event was located in a very bright HII region of its face-on host galaxy, at a projected distance of $\sim$ 10 kpc from its galactic centre. The star-formation rate (SFR$_{GRB} \sim$ 0.2 M$_{\odot}$~yr$^{-1}$) and metallicity (12 + log(O/H) $\sim 8.15 \pm 0.10$) of the GRB star-forming region are consistent with those of the host galaxies of previously studied GRB-SN systems.

We report the detections of molecular hydrogen (H$_2$), vibrationally-excited H$_2$ (H$^*_2$), and neutral atomic carbon (CI), in two new afterglow spectra of GRBs\,181020A ($z=2.938$) and 190114A ($z=3.376$), observed with X-shooter at the Very Large Telescope (VLT). Both host-galaxy absorption systems are characterized by strong damped Lyman-$\alpha$ absorbers (DLAs) and substantial amounts of molecular hydrogen with $\log N$(HI, H$_2$) = $22.20\pm 0.05,~20.40\pm 0.04$ (GRB\,181020A) and $\log N$(HI, H$_2$) = $22.15\pm 0.05,~19.44\pm 0.04$ (GRB\,190114A). The DLA metallicites, depletion levels and dust extinctions are [Zn/H] = $-1.57\pm 0.06$, [Zn/Fe] = $0.67\pm 0.03$, and $A_V = 0.27\pm 0.02$\,mag (GRB\,181020A) and [Zn/H] = $-1.23\pm 0.07$, [Zn/Fe] = $1.06\pm 0.08$, and $A_V = 0.36\pm 0.02$\,mag (GRB\,190114A). We then examine the molecular gas content of all known H$_2$-bearing GRB-DLAs and explore the physical conditions and characteristics of these systems. We confirm that H$_2$ is detected in all CI- and H$^*_2$-bearing GRB absorption systems, but that these rarer features are not necessarily detected in all GRB H$_2$ absorbers. We find that a large molecular fraction of $f_{\rm H_2} \gtrsim 10^{-3}$ is required for CI to be detected. The defining characteristic for H$^*_2$ to be present is less clear, though a large H$_2$ column density is an essential factor. We then derive the H$_2$ excitation temperatures of the molecular gas and find that they are relatively low with $T_{\rm ex} \approx 100 - 300$\,K, however, there could be evidence of warmer components populating the high-$J$ H$_2$ levels in GRBs\,181020A and 190114A. Finally, we demonstrate that the otherwise successful X-shooter GRB afterglow campaign is hampered by a significant dust bias excluding the most dust-obscured H$_2$ absorbers from identification [Abridged].

We present a study of the extinction and depletion-derived dust properties of gamma-ray burst (GRB) absorbers at $1<z<3$ showing the presence of neutral carbon (\ion{C}{I}). By modelling their parametric extinction laws, we discover a broad range of dust models characterizing the GRB \ion{C}{I} absorption systems. In addition to the already well-established correlation between the amount of \ion{C}{I} and visual extinction, $A_V$, we also observe a correlation with the total-to-selective reddening, $R_V$. All three quantities are also found to be connected to the presence and strength of the 2175\,{\AA} dust extinction feature. While the amount of \ion{C}{I} is found to be correlated with the SED-derived dust properties, we do not find any evidence for a connection with the depletion-derived dust content as measured from [Zn/Fe] and $N$(Fe)$_{\rm dust}$. To reconcile this, we discuss a scenario where the observed extinction is dominated by the composition of dust particles confined in the molecular gas-phase of the ISM. We argue that since the depletion level trace non-carbonaceous dust in the ISM, the observed extinction in GRB \ion{C}{I} absorbers is primarily produced by carbon-rich dust in the molecular cloud and is therefore only observable in the extinction curves and not in the depletion patterns. This also indicates that the 2175\,{\AA} dust extinction feature is caused by dust and molecules in the cold and molecular gas-phase. This scenario provides a possible resolution to the discrepancy between the depletion- and SED-derived amounts of dust in high-$z$ absorbers.

(Abridged) Long gamma-ray bursts (LGRB) have been suggested as promising tracers of star formation owing to their association with the core-collapse of massive stars. The goal of this work is to characterise the population of host galaxies of LGRBs at 1 < z < 2, investigate the conditions in which LGRBs form at these redshifts and assess their use as tracers of star formation. We perform a spectro-photometric analysis to determine the stellar mass, star formation rate, specific star formation rate and metallicity of the complete, unbiased host galaxy sample of the Swift/BAT6 LGRB sample at 1 < z < 2. We compare the distribution of these properties to the ones of typical star-forming galaxies from the MOSDEF and COSMOS2015 Ultra Deep surveys, within the same redshift range. We find that, similarly to z < 1, LGRBs do not directly trace star formation at 1 < z < 2, and they tend to avoid high-mass, high-metallicity host galaxies. We also find evidence for an enhanced fraction of starbursts among the LGRB host sample with respect to the star-forming population of galaxies. Nonetheless we demonstrate that the driving factor ruling the LGRB efficiency is metallicity. The LGRB host distributions can be reconciled with the ones expected from galaxy surveys by imposing a metallicity upper limit of 12+logOH ~ 8.55. Metallicity rules the LGRB production efficiency, which is stifled at Z > 0.7 Zsun. Under this hypothesis we can expect LGRBs to trace star formation at z > 3, once the bulk of the star forming galaxy population are characterised by metallicities below this limit. The moderately high metallicity threshold found is in agreement with the conditions necessary to rapidly produce a fast-rotating Wolf-Rayet star a in close binary system, and could be accommodated by single star models under chemically homogeneous mixing with very rapid rotation and weak magnetic coupling.

Long gamma-ray bursts mark the death of massive stars, as revealed by their association with energetic broad-lined stripped-envelope supernovae. The scarcity of nearby events and the brightness of the GRB afterglow, dominating the first days of emission, have so far prevented the study of the very early stages of the GRB-SN evolution. Here we present detailed, multi-epoch spectroscopic observations of SN 2017iuk, associated with GRB 171205A which display features at extremely high expansion velocities of $\sim$ 100,000 km s$^{-1}$ within the first day after the burst. These high-velocity components are characterized by chemical abundances different from those observed in the ejecta of SN 2017iuk at later times. Using spectral synthesis models developed for the SN 2017iuk, we explain these early features as originating not from the supernova ejecta, but from a hot cocoon generated by the energy injection of a mildly-relativistic GRB jet expanding into the medium surrounding the progenitor star. This cocoon becomes rapidly transparent and is outshone by the supernova emission which starts dominating three days after the burst. These results proves that the jet plays an important role not only in powering the GRB event but also its associated supernova.

We present a survey for neutral atomic-carbon (CI) along gamma-ray burst (GRB) sightlines, which probes the shielded neutral gas-phase in the interstellar medium (ISM) of GRB host galaxies at high redshift. We compile a sample of 29 medium- to high-resolution GRB optical afterglow spectra spanning a redshift range through most of cosmic time from $1 < z < 6$. We find that seven ($\approx 25\%$) of the GRBs entering our statistical sample have CI detected in absorption. It is evident that there is a strong excess of cold gas in GRB hosts compared to absorbers in quasar sightlines. We investigate the dust properties of the GRB CI absorbers and find that the amount of neutral carbon is positively correlated with the visual extinction, $A_V$, and the strength of the 2175 \AA\ dust extinction feature, $A_{\mathrm{bump}}$. GRBs with CI detected in absorption are all observed above a certain threshold of $\log N$(HI)$/\mathrm{cm}^{-2}$ + [X/H] > 20.7 and a dust-phase iron column density of $\log N$(Fe)$_{\mathrm{dust}}/\mathrm{cm}^{-2}$ > 16.2. In contrast to the SED-derived dust properties, the strength of the CI absorption does not correlate with the depletion-derived dust properties. This indicates that the GRB CI absorbers trace dusty systems where the dust composition is dominated by carbon-rich dust grains. The observed higher metal and dust column densities of the GRB CI absorbers compared to H$_2$- and CI-bearing quasar absorbers is mainly a consequence of how the two absorber populations are selected, but is also required in the presence of intense UV radiation fields in actively star-forming galaxies.

Here we built up a sample of 22 GRBs at redshifts $z > 2$ observed with X-shooter to determine the abundances of hydrogen, metals, dust, and molecular species. This allows us to study the metallicity and dust depletion effects in the neutral ISM at high redshift and to answer the question whether (and why) there might be a lack of H$_2$ in GRB-DLAs. We fit absorption lines and measure the column densities of different metal species as well as atomic and molecular hydrogen. The derived relative abundances are used to fit dust depletion sequences and determine the dust-to-metals ratio and the host-galaxy intrinsic visual extinction. There is no lack of H$_2$-bearing GRB-DLAs. We detect absorption lines from H$_2$ in 6 out of 22 GRB afterglow spectra, with molecular fractions ranging between $f\simeq 5\cdot10^{-5}$ and $f\simeq 0.04$, and claim tentative detections in three other cases. The GRB-DLAs in the present sample have on average low metallicities ($\mathrm{[X/H]}\approx -1.3$), comparable to the rare population of QSO-ESDLAs (log N(HI) $> 21.5$). H$_2$-bearing GRB-DLAs are found to be associated with significant dust extinction, $A_V > 0.1$ mag, and have dust-to-metals ratios DTM$ > 0.4$. All of these systems exhibit column densities of log N(HI) $> 21.7$. The overall fraction of H$_2$ detections is $\ge 27$% (41% including tentative detections), which is three times larger than in the general population of QSO-DLAs. For $2<z<4$, and for log N(HI) $> 21.7$, the H$_2$ detection fraction is 60-80% in GRB-DLAs as well as in extremely strong QSO-DLAs. This is likely a consequence of the fact that both GRB- and QSO-DLAs with high N(HI) probe sight-lines with small impact parameters that indicate that the absorbing gas is associated with the inner regions of the absorbing galaxy, where the gas pressure is higher and the conversion of HI to H$_2$ takes place.

Evolved asymptotic giant branch (AGB) stars and Type Ia supernovae (SNe) are important contributors to the elements that form dust in the interstellar medium of galaxies, in particular, carbon and iron. However, they require at least a Gyr to start producing these elements, therefore, a change in dust quantity or properties may appear at high redshifts. In this work, we use extinction of gamma-ray burst (GRB) afterglows as a tool to look for variations in dust properties at z>3. We use a spectroscopically selected sample of GRB afterglows observed with the VLT/X-shooter instrument to determine extinction curves out to high redshifts. We present ten new z>3 X-shooter GRBs of which six are dusty. Combining these with individual extinction curves of three previously known z>3 GRBs, we find an average extinction curve consistent with the SMC-Bar. A comparison with spectroscopically selected GRBs at all redshifts indicates a drop in visual extinction (A_V) at z>3.5 with no moderate or high extinction bursts. We check for observational bias using template spectra and find that GRBs up to z~8 are detectable with X-shooter up to A_V~0.3 mag. Although other biases are noted, a uniformly low dust content above z>3.5 indicates a real drop, suggesting a transition in dust properties and/or available dust building blocks. The remarkable increase in dust content at z<3.5 could arise due to carbon and possibly iron production by the first carbon-rich AGB and Type Ia SNe, respectively. Alternatively, z>3.5 dust drop could be the result of low stellar masses of GRB host galaxies.

The gamma-ray bursts hosts (GRBHs) are excellent targets to study the extinction properties of dust and its effects on the global emission of distant galaxies. The dust extinction curve is measured along the GRB afterglow line of sight and the analysis of the spectral energy distribution (SED) of the host galaxy gives access to the global dust attenuation of the stellar light. We selected a sample of 30 GRBs for which the extinction curve along the GRB afterglow line-of-sight (l.o.s) is measured in the rest-frame ultraviolet (UV) up to optical and we analysed the properties of the extinction curve as a function of the host galaxy properties. From these 30 GRBs, we selected seven GRBHs with a good rest-frame UV to near infra-red spectral coverage for the host. The attenuation curve was derived by fitting the SEDs of the GRBH sample with the CIGALE SED fitting code. Different star formation histories were studied to recover the star formation rates (SFR) derived using H$_{\alpha}$ luminosities. The most extinguished GRBs are preferentially found in the more massive hosts and the UV bump is preferentially found in the most extinguished GRB l.o.s. Five out of seven hosts are best fitted with a recent burst of star formation, leading to lower stellar mass estimates than previously found. The average attenuation in the host galaxies is about 70% of the amount of extinction along the GRB l.o.s. We find a great variety in the derived attenuation curves of GRBHs, the UV slope can be similar, flatter or even steeper than the extinction curve slope. We find that the flatter (steeper) attenuation curves are found in galaxies with the highest (lowest) SFR and stellar masses. The comparison of our results with radiative transfer simulations leads to a uniform distribution of dust and stars in a very clumpy ISM for half the GRBHs and various dust-stars geometries for the second half of the sample.

In the era of Extremely Large Telescopes, the current generation of 8-10m facilities are likely to remain competitive at far-blue visible wavelengths for the foreseeable future. High-efficiency (>20%) observations of the ground UV (300-400 nm) at medium resolving power (R~20,000) are required to address a number of exciting topics in stellar astrophysics, while also providing new insights in extragalactic science. Anticipating strong demand to better exploit this diagnostic-rich wavelength region, we revisit the science case and instrument requirements previously assembled for the CUBES concept for the Very Large Telescope.

Broad-line Ic supernovae (SNe Ic-BL) are a very rare class of core-collapse supernovae exhibiting high ejecta velocities and high kinetic energies. They are the only type of SNe that accompany long gamma-ray burst (GRB) explosions. Systematic differences found in the spectra of SNe Ic-BL with and without GRBs (GRB-SNe and SNe Ic-BL, respectively) suggest that either the progenitor or/and explosion mechanism of SNe Ic-BL with and without a GRB differ, or the difference could be only due to the viewing angle of the observer with respect to the orientation of the collimated explosion. We present the systematic comparison of the host galaxies of broad-lined SNe Ic with and without a detected GRB, the latter being detected in untargeted surveys, with the aim to find out whether there are any systematic differences between the environments in which these two classes of SNe preferentially explode. We study photometric properties of the host galaxies of a sample of 8 GRB-SNe and a sample of 28 SNe Ic-BL at z < 0.2. The two galaxy samples have indistinguishable luminosity and proper size distribution. We find indications that GRB-SNe on average occur closer to the centres of their host galaxies, i.e. the samples have a different distribution of projected offsets, normalized by the galaxy sizes. In addition we compare gas-phase metallicities of the GRB-SNe and SNe Ic-BL host samples and find that a larger fraction of super-solar metallicity hosts are found among the SNe Ic-BL without a GRB. Our results are indicative of a genuine difference between the two types of explosions and suggest that the viewing angle is not the main source of difference in the spectra of the two classes. We discuss the implications our results have on our understanding of progenitors of SNe Ic-BL with and without a GRB.

The UV extinction feature at 2175 \AA\ is ubiquitously observed in the Galaxy but is rarely detected at high redshifts. Here we report the spectroscopic detection of the 2175 \AA\ bump on the sightline to the \gamma-ray burst (GRB) afterglow GRB 180325A at z=2.2486, the only unambiguous detection over the past ten years of GRB follow-up, at four different epochs with the Nordic Optical Telescope (NOT) and the Very Large Telescope (VLT)/X-shooter. Additional photometric observations of the afterglow are obtained with the Gamma-Ray burst Optical and Near-Infrared Detector (GROND). We construct the near-infrared to X-ray spectral energy distributions (SEDs) at four spectroscopic epochs. The SEDs are well-described by a single power-law and an extinction law with R_V~4.4, A_V~1.5, and the 2175 \AA\ extinction feature. The bump strength and extinction curve are shallower than the average Galactic extinction curve. We determine a metallicity of [Zn/H]>-0.98 from the VLT/X-shooter spectrum. We detect strong neutral carbon associated with the GRB with an equivalent width of Wr(\lambda 1656) = 0.85+/-0.05. We also detect optical emission lines from the host galaxy. Based on the H\alpha emission line flux, the derived dust-corrected star-formation rate is ~46+/-4 M_sun/yr and the predicted stellar mass is log M*/M_sun~9.3+/-0.4, suggesting the host galaxy is amongst the main-sequence star-forming galaxies.

We present the consolidated scientific case for multi-object spectroscopy with the MOSAIC concept on the European ELT. The cases span the full range of ELT science and require either 'high multiplex' or 'high definition' observations to best exploit the excellent sensitivity and wide field-of-view of the telescope. Following scientific prioritisation by the Science Team during the recent Phase A study of the MOSAIC concept, we highlight four key surveys designed for the instrument using detailed simulations of its scientific performance. We discuss future ways to optimise the conceptual design of MOSAIC in Phase B, and illustrate its competitiveness and unique capabilities by comparison with other facilities that will be available in the 2020s.

We present here a survey of high-ionization absorption lines in the afterglow spectra of long-duration gamma-ray bursts (GRBs) obtained with the VLT/X-shooter spectrograph. Our main goal is to investigate the circumburst medium in the natal regions of GRBs. Our primary focus is on the NV 1238,1242 line transitions, but we also discuss other high-ionization lines such as OVI, CIV and SiIV. We find no correlation between the column density of NV and the neutral gas properties such as metallicity, HI column density and dust depletion, however the relative velocity of NV, typically a blueshift with respect to the neutral gas, is found to be correlated with the column density of HI. This may be explained if the NV gas is part of an HII region hosting the GRB, where the region's expansion is confined by dense, neutral gas in the GRB's host galaxy. We find tentative evidence (at 2-sigma significance) that the X-ray derived column density, N_H,X, may be correlated with the column density of NV, which would indicate that both measurements are sensitive to the column density of the gas located in the vicinity of the GRB. We investigate the scenario where NV (and also OVI) is produced by recombination after the corresponding atoms have been stripped entirely of their electrons by the initial prompt emission, in contrast to previous models where highly-ionized gas is produced by photoionization from the GRB afterglow.

The part played by stars in the ionization of the intergalactic medium remains an open question. A key issue is the proportion of the stellar ionizing radiation that escapes the galaxies in which it is produced. Spectroscopy of gamma-ray burst afterglows can be used to determine the neutral hydrogen column-density in their host galaxies and hence the opacity to extreme ultra-violet radiation along the lines-of-sight to the bursts. Thus, making the reasonable assumption that long-duration GRB locations are representative of the sites of massive stars that dominate EUV production, one can calculate an average escape fraction of ionizing radiation in a way that is independent of galaxy size, luminosity or underlying spectrum. Here we present a sample of NH measures for 138 GRBs in the range 1.6<z<6.7 and use it to establish an average escape fraction at the Lyman limit of <fesc>~0.005, with a 98% confidence upper limit of ~0.015. This analysis suggests that stars provide a small contribution to the ionizing radiation budget of the IGM at z<5, where the bulk of the bursts lie. At higher redshifts, z>5, firm conclusions are limited by the small size of the GRB sample, but any decline in average HI column-density seems to be modest. We also find no indication of a significant correlation of NH with galaxy UV luminosity or host stellar mass, for the subset of events for which these are available. We discuss in some detail a number of selection effects and potential biases. Drawing on a range of evidence we argue that such effects, while not negligible, are unlikely to produce systematic errors of more than a factor ~2, and so would not affect the primary conclusions. Given that many GRB hosts are low metallicity, high specific star-formation rate, dwarf galaxies, these results present a particular problem for the hypothesis that such galaxies dominated the reionization of the universe.

The extinction profiles in Gamma-Ray Burst (GRB) afterglow spectral energy distributions (SEDs) are usually described by the Small Magellanic Cloud (SMC)-type extinction curve. In different empirical extinction laws, the total-to-selective extinction, RV, is an important quantity because of its relation to dust grain sizes and compositions. We here analyse a sample of 17 GRBs (0.34<z<7.84) where the ultraviolet to near-infrared spectroscopic observations are available through the VLT/X-shooter instrument, giving us an opportunity to fit individual extinction curves of GRBs for the first time. Our sample is compiled on the basis that multi-band photometry is available around the X-shooter observations. The X-shooter data are combined with the Swift X-ray data and a single or broken power-law together with a parametric extinction law is used to model the individual SEDs. We find 10 cases with significant dust, where the derived extinction, AV, ranges from 0.1-1.0mag. In four of those, the inferred extinction curves are consistent with the SMC curve. The GRB individual extinction curves have a flat RV distribution with an optimal weighted combined value of RV = 2.61+/-0.08 (for seven broad coverage cases). The 'average GRB extinction curve' is similar to, but slightly steeper than the typical SMC, and consistent with the SMC Bar extinction curve at ~95% confidence level. The resultant steeper extinction curves imply populations of small grains, where large dust grains may be destroyed due to GRB activity. Another possibility could be that young age and/or lower metallicities of GRBs environments are responsible for the steeper curves.

In this work we present spectra of all $\gamma$-ray burst (GRB) afterglows that have been promptly observed with the X-shooter spectrograph until 31-03-2017. In total, we obtained spectroscopic observations of 103 individual GRBs observed within 48 hours of the GRB trigger. Redshifts have been measured for 97 per cent of these, covering a redshift range from 0.059 to 7.84. Based on a set of observational selection criteria that minimize biases with regards to intrinsic properties of the GRBs, the follow-up effort has been focused on producing a homogeneous sample of 93 afterglow spectra for GRBs discovered by the Swift satellite. We here provide a public release of all the reduced spectra, including continuum estimates and telluric absorption corrections. For completeness, we also provide reductions for the 18 late-time observations of the underlying host galaxies. We provide an assessment of the degree of completeness with respect to the parent GRB population, in terms of the X-ray properties of the bursts in the sample and find that the sample presented here is representative of the full Swift sample. We constrain the fraction of dark bursts to be < 28 per cent and we confirm previous results that higher optical darkness is correlated with increased X-ray absorption. For the 42 bursts for which it is possible, we provide a measurement of the neutral hydrogen column density, increasing the total number of published HI column density measurements by $\sim$ 33 per cent. This dataset provides a unique resource to study the ISM across cosmic time, from the local progenitor surroundings to the intervening universe.

We report the detection of the CO(3-2) emission line from the host galaxy of Gamma Ray Burst (GRB) 080207 at $z$ = 2.086. This is the first detection of molecular gas in emission from a GRB host galaxy beyond redshift 1. We find this galaxy to be rich in molecular gas with a mass of $1.1 \times 10^{11}\,\rm M_{\odot}$ assuming $\alpha_{\rm CO}=$ 4.36 $\rm M_{\odot}(\rm K\,km\,s^{-1}\,pc^2)^{-1}$. The molecular gas mass fraction of the galaxy is $\sim$ 0.5, typical of star forming galaxies (SFGs) with similar stellar masses and redshifts. With a $\rm SFR_{FIR}$ of 260 $\rm M_{\odot}\,yr^{-1}$, we measure a molecular-gas-depletion timescale of 0.43 Gyr, near the peak of the depletion timescale distribution of SFGs at similar redshifts. Our findings are therefore in contradiction with the proposed molecular gas deficiency in GRB host galaxies. We argue that the reported molecular gas deficiency for GRB hosts could be the artifact of improper comparisons or neglecting the effect of the typical low metallicities of GRB hosts on the CO-to-molecular-gas conversion factor. We also compare the kinematics of the CO(3-2) emission line to that of the H$\alpha$ emission line from the host galaxy. We find the H$\alpha$ emission to have contributions from two separate components, a narrow and a broad one. The narrow component matches the CO emission well in velocity space. The broad component, with a FWHM of $\sim$ 1100 $\rm km\,s^{-1}$, is separated by $+390$ $\rm km\,s^{-1}$ in velocity space from the narrow component. We speculate this broad component to be associated with a powerful outflow in the host galaxy or in an interacting system.

We present a comprehensive temporal and spectral analysis of the long Swift GRB 120327A afterglow data to investigate the possible causes of the observed early time colour variations. We collected data from various instruments/telescopes in different bands (X-rays, ultra- violet, optical and near-infrared) and determined the shapes of the afterglow early-time light curves. We studied the overall temporal behaviour and the spectral energy distributions from early to late times. The ultra-violet, optical, and near-infrared light curves can be modelled with a single power-law component between 200 and 2e4 s after the burst event. The X-ray light curve shows a canonical steep-shallow-steep behaviour, typical of long gamma-ray bursts. At early times a colour variation is observed in the ultra-violet/optical bands, while at very late times a hint of a re-brightening is visible. The observed early time colour change can be explained as a variation in the intrinsic optical spectral index, rather than an evolution of the optical extinction.

Previous studies have shown that the incidence rate of intervening strong MgII absorbers towards GRBs were a factor of 2 - 4 higher than towards quasars. Exploring the similar sized and uniformly selected legacy data sets XQ-100 and XSGRB, each consisting of 100 quasar and 81 GRB afterglow spectra obtained with a single instrument (VLT/X-shooter), we demonstrate that there is no disagreement in the number density of strong MgII absorbers with rest-frame equivalent widths $W_r^{2796} >$ 1 {\AA} towards GRBs and quasars in the redshift range 0.1 < z < 5. With large and similar sample sizes, and path length coverages of $\Delta$z = 57.8 and 254.4 for GRBs and quasars, respectively, the incidences of intervening absorbers are consistent within 1 sigma uncertainty levels at all redshifts. For absorbers at z < 2.3 the incidence towards GRBs is a factor of 1.5$\pm$0.4 higher than the expected number of strong MgII absorbers in SDSS quasar spectra, while for quasar absorbers observed with X-shooter we find an excess factor of 1.4$\pm$0.2 relative to SDSS quasars. Conversely, the incidence rates agree at all redshifts with reported high spectral resolution quasar data, and no excess is found. The only remaining discrepancy in incidences is between SDSS MgII catalogues and high spectral resolution studies. The rest-frame equivalent width distribution also agrees to within 1 sigma uncertainty levels between the GRB and quasar samples. Intervening strong MgII absorbers towards GRBs are therefore neither unusually frequent, nor unusually strong.

It is notoriously difficult to localize short $\gamma$-ray bursts (sGRBs) and their hosts to measure their redshifts. These measurements, however, are critical to constrain the nature of sGRB progenitors, their redshift distribution and the $r$-process element enrichment history of the universe. Here, we present spectroscopy of the host galaxy of GRB111117A and measure its redshift to be $z = 2.211$. This makes GRB111117A the most distant high-confidence short duration GRB detected to date. Our spectroscopic redshift supersedes a lower, previously estimated photometric redshift value for this burst. We use the spectroscopic redshift, as well as new imaging data to constrain the nature of the host galaxy and the physical parameters of the GRB. The rest-frame X-ray derived hydrogen column density, for example, is the highest compared to a complete sample of sGRBs and seems to follow the evolution with redshift as traced by the hosts of long GRBs (lGRBs). The host lies in the brighter end of the expected sGRB host brightness distribution at $z = 2.211$, and is actively forming stars. Using the host as a benchmark for redshift determination, we find that between 43 and 71 per cent of all sGRB redshifts should be missed due to host faintness for hosts at $z\sim2$. The high redshift of GRB111117A is evidence against a lognormal delay-time model for sGRBs through the predicted redshift distribution of sGRBs, which is very sensitive to high-$z$ sGRBs. From the age of the universe at the time of GRB explosion, an initial neutron star (NS) separation of $a_0 < 3.2~R_\odot$ is required in the case where the progenitor system is a circular pair of inspiralling NSs. This constraint excludes some of the longest sGRB formation channels for this burst.

We present a catalog of early-time (~10^2-10^4s) photometry and polarimetry of all Gamma-Ray Burst (GRB) optical afterglows observed with RINGO2 imaging polarimeter on the Liverpool Telescope. For the 19 optical afterglows observed, the following 9 were bright enough to perform photometry and attempt polarimetry: GRB 100805A, GRB 101112A, GRB 110205A, GRB 110726A, GRB 120119A, GRB 120308A, GRB 120311A, GRB 120326A and GRB 120327A. We present multi-wavelength light curves for these 9 GRBs, together with estimates of their optical polarization degrees and/or limits. We carry out a thorough investigation of detection probabilities, instrumental properties and systematics. Using two independent methods, we confirm previous reports of significant polarization in GRB 110205A and 120308A, and report new detection of $P=6^{+3}_{-2}\%$ in GRB101112A. We discuss the results for the sample in the context of the reverse and forward shock afterglow scenario, and show that GRBs with detectable optical polarization at early time have clearly identifiable signatures of reverse-shock emission in their optical light curves. This supports the idea that GRB ejecta contain large-scale magnetic fields and highlights the importance of rapid-response polarimetry.

Star forming galaxies have long been considered the dominant sources of the cosmic ultraviolet background radiation at early epochs. However, observing and characterizing the galaxy population with significant ionizing emission has proven to be challenging. In particular, the fraction of ionizing radiation that escapes the local environment to the intergalactic medium is poorly known. We investigate the relation between the escape fraction and galaxy luminosity. We combine the deep ultraviolet observations of Hubble Ultra Deep Field (UVUDF) with the deep Multi Unit Spectroscopic Explorer (MUSE) observations of the same field, collecting a sample of 165 faint star forming galaxies in the $3 < z < 4$ redshift range with deep rest-frame observations of the Lyman continuum. In our sample, we do not find any galaxy with significant emission of LyC radiation. We bin the galaxies in various redshift and brightness intervals and stack their images. From stacked images we estimate the relative escape fraction upper limits as a function of the luminosity. Thanks to the depth of the sample we measure meaningful 1$\sigma$ upper limits of $f_{esc,rel} < 0.07, 0.2$ and 0.6 at $L \sim L_{\rm z=3}^{*}, 0.5L_{\rm z=3}^{*}$ and $0.1L_{\rm z=3}^{*}$, respectively. We use our estimates and theoretical predictions from the literature to study a possible dependence of the escape fraction on galaxy luminosity by modelling the ionizing background with different prescriptions of $f_{\rm esc} (M_{\rm UV})$. We show that the understanding of the luminosity dependence hinges on the ability to constrain the escape fraction down to $M_{\rm UV} \sim -18$ mag in the future.

We investigate the existence of a metallicity threshold for the production of long gamma-ray bursts (LGRBs). We used the host galaxies of the Swift/BAT6 sample of LGRBs. We considered the stellar mass, star formation rate (SFR), and metallicity determined from the host galaxy photometry and spectroscopy up to z = 2 and used them to compare the distribution of host galaxies to that of field galaxies in the mass-metallicity and fundamental metallicity relation plane. We find that although LGRBs also form in galaxies with relatively large stellar masses, the large majority of host galaxies have metallicities below log(O=H)~8.6. The extension to z = 2 results in a good sampling of stellar masses also above Log(Mstar/Msun)~9.5 and provides evidence that LGRB host galaxies do not follow the fundamental metallicity relation. As shown by the comparison with dedicated numerical simulations of LGRB host galaxy population, these results are naturally explained by the existence of a mild (~0.7 Zsun) threshold for the LGRB formation. The present statistics does not allow us to discriminate between different shapes of the metallicity cutoff, but the relatively high metallicity threshold found in this work is somewhat in disagreement to most of the standard single-star models for LGRB progenitors.

We present the extensive follow-up campaign on the afterglow of GRB 110715A at 17 different wavelengths, from X-ray to radio bands, starting 81 seconds after the burst and extending up to 74 days later. We performed for the first time a GRB afterglow observation with the ALMA observatory. We find that the afterglow of GRB 110715A is very bright at optical and radio wavelengths. We use optical and near infrared spectroscopy to provide further information about the progenitor's environment and its host galaxy. The spectrum shows weak absorption features at a redshift z = 0.8225, which reveal a host galaxy environment with low ionization, column density and dynamical activity. Late deep imaging shows a very faint galaxy, consistent with the spectroscopic results. The broadband afterglow emission is modelled with synchrotron radiation using a numerical algorithm and we determine the best fit parameters using Bayesian inference in order to constrain the physical parameters of the jet and the medium in which the relativistic shock propagates. We fitted our data with a variety of models, including different density profiles and energy injections. Although the general behaviour can be roughly described by these models, none of them are able to fully explain all data points simultaneously. GRB 110715A shows the complexity of reproducing extensive multi-wavelength broadband afterglow observations, and the need of good sampling in wavelength and time and more complex models to accurately constrain the physics of GRB afterglows.

Long gamma-ray bursts (LGRBs) and superluminous supernovae (SLSNe) are both explosive transients with very massive progenitor stars. Clues about the nature of the progenitors can be found by investigating environments in which such transients occur. While studies of LGRB host galaxies have a long history, dedicated observational campaigns have only recently resulted in a high enough number of photometrically and spectroscopically observed SLSN hosts to allow statistically significant analysis of their properties. In this paper we make a comparison of the host galaxies of hydrogen-poor (H-poor) SLSNe and the Swift/BAT6 sample of LGRBs. In contrast to previous studies we use a complete sample of LGRBs and we address a special attention to the comparison methodology and the selection of SLSN sample whose data have been compiled from the available literature. At intermediate redshifts (0.3 < z < 0.7) the two classes of transients select galaxies whose properties (stellar mass, luminosity, star-formation rate, specific star-formation rate and metallicity) do not differ on average significantly. Moreover, the host galaxies of both classes of objects follow the fundamental metallicity relation and the fundamental plane of metallicity. In contrast to previous studies we show that at intermediate redshifts the emission line equivalent widths of the two populations are essentially the same and that the previous claims regarding the higher fraction of SLSN hosts among the extreme emission line galaxies with respect to LGRBs are mostly due to a larger fraction of strong-line emitters among SLSN hosts at z < 0.3, where samples of LGRB hosts are small and poorly defined.

Using VLT/X-shooter we searched for emission line galaxies associated to four damped Lyman-$\alpha$ systems (DLAs) and one sub-DLA at 2.73<=z<=3.25 towards QSO J2358+0149. We detect [O III] emission from a "low-cool" DLA at z_abs = 2.9791 (having log N(HI)=21.69+\-0.10, [Zn/H] = -1.83+\-0.18) at an impact parameter of, $\rho$ ~12 kpc. The associated galaxy is compact with a dynamical mass of (1-6)x10^9 M_solar, very high excitation ([O III]/[O II] and [O III]/[H$\beta$] both greater than 10), 12+[O/H]<=8.5 and moderate star formation rate (SFR <=2 M_solar yr^{-1}). Such properties are typically seen in the low-z extreme blue compact dwarf galaxies. The kinematics of the gas is inconsistent with that of an extended disk and the gas is part of either a large scale wind or cold accretion. We detect Ly$\alpha$ emission from the z_abs = 3.2477 DLA (having log N(HI)=21.12+\-0.10 and [Zn/H]=-0.97+\-0.13).The Ly$\alpha$ emission is redshifted with respect to the metal absorption lines by 320 km s^{-1}, consistent with the location of the red hump expected in radiative transport models. We derive SFR ~0.2-1.7 M_solar yr^{-1} and Ly$\alpha$ escape fraction of >=10 per cent. No other emission line is detected from this system. Because the DLA has a small velocity separation from the quasar (~500 km s^{-1}) and the DLA emission is located within a small projected distance ($\rho<5$ kpc), we also explore the possibility that the Ly$\alpha$ emission is being induced by the QSO itself. QSO induced Ly$\alpha$ fluorescence is possible if the DLA is within a physical separation of 340 kpc to the QSO. Detection of stellar continuum light and/or the oxygen emission lines would disfavor this possibility. We do not detect any emission line from the remaining three systems.

Long gamma-ray bursts (LGRBs) are associated with the deaths of massive stars and could thus be a potentially powerful tool to trace cosmic star formation. However, especially at low redshifts (z < 1.5) LGRBs seem to prefer particular types of environment. Our aim is to study the host galaxies of a complete sample of bright LGRBs to investigate the impact of the environment on GRB formation. We study host galaxy spectra of the Swift/BAT6 complete sample of 14 z < 1 bright LGRBs. We use the detected nebular emission lines to measure the dust extinction, star formation rate (SFR) and nebular metallicity (Z) of the hosts and supplement the data set with previously measured stellar masses M$_{\star}$. The distributions of the obtained properties and their interrelations (e.g. mass-metallicity and SFR-M$_{\star}$ relations) are compared to samples of field star-forming galaxies.We find that LGRB hosts at z < 1 have on average lower SFRs than if they were direct star-formation tracers. By directly comparing metallicity distributions of LGRB hosts and star-forming galaxies, we find a good match between the two populations up to Z $\sim 8.4-8.5$, after which the paucity of metal-rich LGRB hosts becomes apparent. The LGRB host galaxies of our complete sample are not inconsistent with the mass-metallicity relation at similar mean redshift and stellar masses. The cutoff against high metallicities (and high masses) can explain the low SFR values of LGRB hosts. We find a hint of increased incidence of starburst galaxies in the Swift/BAT6 z < 1 sample with respect to that of a field star-forming population. Given that the SFRs are low on average, the latter is ascribed to low stellar masses. Nevertheless the limits on the completeness and metallicity availability of current surveys, coupled with the limited number of LGRB host galaxies, prevent us from reaching more quantitative conclusions.

[Abridged] We present data and initial results from VLT/X-Shooter emission-line spectroscopy of 96 GRB-selected galaxies at 0.1<z<3.6, the largest sample of GRB host spectroscopy available to date. Most of our GRBs were detected by Swift and 76% are at 0.5<z<2.5 with a median z~1.6. Based on Balmer and/or forbidden lines of oxygen, nitrogen, and neon, we measure systemic redshifts, star formation rates (SFRs), visual attenuations, oxygen abundances (12+log(O/H)), and emission-line widths. We find a strong change of the typical physical properties of GRB hosts with redshift. The median SFR, for example, increases from ~0.6 M_sun/yr at z~0.6 up to ~15 M_sun/yr at z~2. A higher ratio of [OIII]/[OII] at higher redshifts leads to an increasing distance of GRB-selected galaxies to the locus of local galaxies in the BPT diagram. Oxygen abundances of the galaxies are distributed between 12+log(O/H)=7.9 and 12+log(O/H)=9.0 with a median of 12+log(O/H)~8.5. The fraction of GRB-selected galaxies with super-solar metallicities is around 20% at z<1 in the adopted metallicity scale. This is significantly less than the fraction of star formation in similar galaxies, illustrating that GRBs are scarce in high-metallicity environments. At z~3, sensitivity limits us to probing only the most luminous GRB hosts for which we derive metallicities of Z ~< 0.5 Z_sun. Together with a high incidence of galaxies with similar metallicity in our sample at z~1.5, this indicates that the metallicity dependence observed at low redshift will not be dominant at z~3.

Gamma-ray burst GRB 140430A was detected by the Swift satellite and observed promptly with the imaging polarimeter RINGO3 mounted on the Liverpool Telescope, with observations beginning while the prompt $\gamma$-ray emission was still ongoing. In this paper, we present densely sampled (10-second temporal resolution) early optical light curves in 3 optical bands and limits to the degree of optical polarization. We compare optical, X-ray and gamma-ray properties and present an analysis of the optical emission during a period of high-energy flaring. The complex optical light curve cannot be explained merely with a combination of forward and reverse shock emission from a standard external shock, implying additional contribution of emission from internal shock dissipation. We estimate an upper limit for time averaged optical polarization during the prompt phase to be as low as P < 12% (1$\sigma$). This suggests that the optical flares and early afterglow emission in this GRB are not highly polarized. Alternatively, time averaging could mask the presence of otherwise polarized components of distinct origin at different polarization position angles.

LGRBs are associated with massive stars and are therefore linked to star formation. The conditions necessary to produce LGRBs can affect the relation between the LGRB rate and star formation. By using the power of a complete LGRB sample, our aim is to understand whether such a bias exists and, if it does, what is its origin. In this first paper, we build the SED of the z<1 host galaxies of the BAT6 LGRB sample, and determine their stellar masses from SED fitting. We compare the resulting stellar mass distribution (i) with star-forming galaxies observed in deep surveys (UltraVISTA); (ii) with semi-analitical models of the z<1 star forming galaxy population and (iii) with numerical simulations of LGRB hosts having different metallicity thresholds for the progenitor star environment. We find that at z<1 LGRBs tend to avoid massive galaxies and are powerful in selecting faint low-mass star-forming galaxies. The stellar mass distribution of the hosts is not consistent with that of the UltraVISTA star-forming galaxies weighted for their SFR. This implies that, at least at z<1, LGRBs are not unbiased tracers of star formation. To make the two distributions consistent, a much steeper faint-end of the mass function would be required, or a very shallow SFR-Mass relation for the low mass galaxy population. GRB host galaxy simulations indicates that, to reproduce the stellar mass distribution, a metallicity threshold of the order of Z_th=0.3-0.5Z_sun is necessary. Models without a metallicity threshold or with an extreme threshold of Z_th = 0.1Z_sun are excluded at z<1. The use of the BAT6 complete sample makes this result not affected by possible biases which could have influenced past results based on incomplete samples. The preference for low metallicities (Z<~0.5Z_sun) can be a consequence of the particular conditions needed for the progenitor star to produce a GRB. (Abridged)

Long gamma-ray bursts (LGRBs) are associated with the death of massive stars. Their host galaxies therefore represent a unique class of objects tracing star formation across the observable Universe. Indeed, recently accumulated evidence shows that GRB hosts do not differ substantially from general population of galaxies at high (z > 2) redshifts. However, it has been long recognised that the properties of z < 1.5 hosts, compared to general star-forming population, are unusual. To better understand the reasons for the supposed difference in LGRB hosts properties at z < 1.5, we obtained VLT/X- Shooter spectra of six hosts lying in the redshift range of 0.8 < z < 1.3. Some of these hosts have been observed before, yet we still lack well constrained information on their characteristics such as metallicity, dust extinction and star formation rate. We search for emission lines in the VLT/X-Shooter spectra of the hosts and measure their fluxes. We perform a detailed analysis, estimating host average extinction, star-formation rates, metallicities and electron densities where possible. Measured quantities of our hosts are compared to a larger sample of previously observed GRB hosts at z < 2. Star-formation rates and metallicities are measured for all the hosts analyzed in this paper and metallicities are well determined for 4 hosts. The mass-metallicity relation, the fundamental metallicity relation and SFRs derived from our hosts occupy similar parameter space as other host galaxies investigated so-far at the same redshift. We therefore conclude that GRB hosts in our sample support the found discrepancy between the properties of low-redshift GRB hosts and the general population of star- forming galaxies.

High-redshift gamma-ray bursts have several advantages for the study of the distant universe, providing unique information about the structure and properties of the galaxies in which they exploded. Spectroscopic identification with large ground-based telescopes has improved our knowledge of the class of such distant events. We present the multi-wavelength analysis of the high-$z$ Swift gamma-ray burst GRB140515A ($z = 6.327$). The best estimate of the neutral hydrogen fraction of the intergalactic medium (IGM) towards the burst is $x_{HI} \leq 0.002$. The spectral absorption lines detected for this event are the weakest lines ever observed in gamma-ray burst afterglows, suggesting that GRB140515A exploded in a very low density environment. Its circum-burst medium is characterised by an average extinction (A$_{\rm V} \sim 0.1$) that seems to be typical of $z \ge 6$ events. The observed multi-band light curves are explained either with a very flat injected spectrum ($p = 1.7$) or with a multi-component emission ($p = 2.1$). In the second case a long-lasting central engine activity is needed in order to explain the late time X-ray emission. The possible origin of GRB140515A from a Pop III (or from a Pop II stars with local environment enriched by Pop III) massive star is unlikely.

Gamma Ray Bursts are detectable in the gamma-ray band if their jets are oriented towards the observer. However, for each GRB with a typical theta_jet, there should be ~2/theta_jet^2 bursts whose emission cone is oriented elsewhere in space. These off-axis bursts can be eventually detected when, due to the deceleration of their relativistic jets, the beaming angle becomes comparable to the viewing angle. Orphan Afterglows (OA) should outnumber the current population of bursts detected in the gamma-ray band even if they have not been conclusively observed so far at any frequency. We compute the expected flux of the population of orphan afterglows in the mm, optical and X-ray bands through a population synthesis code of GRBs and the standard afterglow emission model. We estimate the detection rate of OA by on-going and forthcoming surveys. The average duration of OA as transients above a given limiting flux is derived and described with analytical expressions: in general OA should appear as daily transients in optical surveys and as monthly/yearly transients in the mm/radio band. We find that ~ 2 OA yr^-1 could already be detected by Gaia and up to 20 OA yr^-1 could be observed by the ZTF survey. A larger number of 50 OA yr^-1 should be detected by LSST in the optical band. For the X-ray band, ~ 26 OA yr^-1 could be detected by the eROSITA. For the large population of OA detectable by LSST, the X-ray and optical follow up of the light curve (for the brightest cases) and/or the extensive follow up of their emission in the mm and radio band could be the key to disentangle their GRB nature from other extragalactic transients of comparable flux density.

We present predictions of centimeter and millimeter radio emission from reverse shocks in the early afterglows of gamma-ray bursts with the goal of determining their detectability with current and future radio facilities. Using a range of GRB properties, such as peak optical brightness and time, isotropic equivalent gamma-ray energy and redshift, we simulate radio light curves in a framework generalized for any circumburst medium structure and including a parametrization of the shell thickness regime that is more realistic than the simple assumption of thick- or thin-shell approximations. Building on earlier work by Mundell et al. (2007) and Melandri et al. (2010) in which the typical frequency of the reverse shock was suggested to lie at radio, rather than optical wavelengths at early times, we show that the brightest and most distinct reverse-shock radio signatures are detectable up to 0.1 -- 1 day after the burst, emphasizing the need for rapid radio follow-up. Detection is easier for bursts with later optical peaks, high isotropic energies, lower circumburst medium densities, and at observing frequencies that are less prone to synchrotron self-absorption effects - typically above a few GHz. Given recent detections of polarized prompt gamma-ray and optical reverse-shock emission, we suggest that detection of polarized radio/mm emission will unambiguously confirm the presence of low-frequency reverse shocks at early time.

In this work we use gamma-ray burst (GRB) afterglow spectra observed with the VLT/X-shooter spectrograph to measure rest-frame extinction in GRB lines-of-sight by modeling the broadband near-infrared (NIR) to X-ray afterglow spectral energy distributions (SEDs). Our sample consists of nine Swift GRBs, eight of them belonging to the long-duration and one to the short-duration class. Dust is modeled using the average extinction curves of the Milky Way and the two Magellanic Clouds. We derive the rest-frame extinction of the entire sample, which fall in the range $0 \lesssim {\it A}_{\rm V} \lesssim 1.2$. Moreover, the SMC extinction curve is the preferred extinction curve template for the majority of our sample, a result which is in agreement with those commonly observed in GRB lines-of-sights. In one analysed case (GRB 120119A), the common extinction curve templates fail to reproduce the observed extinction. To illustrate the advantage of using the high-quality X-shooter afterglow SEDs over the photometric SEDs, we repeat the modeling using the broadband SEDs with the NIR-to-UV photometric measurements instead of the spectra. The main result is that the spectroscopic data, thanks to a combination of excellent resolution and coverage of the blue part of the SED, are more successful in constraining the extinction curves and therefore the dust properties in GRB hosts with respect to photometric measurements. In all cases but one the extinction curve of one template is preferred over the others. We show that the modeled values of the extinction and the spectral slope, obtained through spectroscopic and photometric SED analysis, can differ significantly for individual events. Finally we stress that, regardless of the resolution of the optical-to-NIR data, the SED modeling gives reliable results only when the fit is performed on a SED covering a broader spectral region.

The long ${\it Swift}$ gamma-ray burst GRB 120326A at redshift $z=1.798$ exhibited a multi-band light curve with a striking feature: a late-time, long-lasting achromatic rebrightening, rarely seen in such events. Peaking in optical and X-ray bands $\sim 35$ ks ($\sim 12.5$ ks in the GRB rest frame) after the 70-s GRB prompt burst, the feature brightens nearly two orders of magnitude above the underlying optical power-law decay. Modelling the multiwavelength light curves, we investigate possible causes of the rebrightening in the context of the standard fireball model. We exclude a range of scenarios for the origin of this feature: reverse-shock flash, late-time forward shock peak due to the passage of the maximal synchrotron frequency through the optical band, late central engine optical/X-ray flares, interaction between the expanding blast wave and a density enhancement in the circumburst medium and gravitational microlensing. Instead we conclude that the achromatic rebrightening may be caused by a refreshed forward shock or a geometrical effect. In addition, we identify an additional component after the end of the prompt emission, that shapes the observed X-ray and optical light curves differently, ruling out a single overall emission component to explain the observed early time emission.

Context. Gamma-ray burst (GRBs) afterglows probe sightlines to star-forming regions in distant star-forming galaxies. Here we present a study of the peculiar afterglow spectrum of the z = 0.889 Swift GRB 140506A. Aims. Our aim is to understand the origin of the very unusual properties of the absorption along the line-of-sight. Methods. We analyse spectroscopic observations obtained with the X-shooter spectrograph mounted on the ESO/VLT at two epochs 8.8 h and 33 h after the burst as well as imaging from the GROND instrument. We also present imaging and spectroscopy of the host galaxy obtained with the Magellan telescope. Results. The underlying afterglow appears to be a typical afterglow of a long-duration GRB. However, the material along the line-of- sight has imprinted very unusual features on the spectrum. Firstly, there is a very broad and strong flux drop below 8000 AA (4000 AA in the rest frame), which seems to be variable between the two spectroscopic epochs. We can reproduce the flux-drops both as a giant 2175 AA extinction bump and as an effect of multiple scattering on dust grains in a dense environment. Secondly, we detect absorption lines from excited H i and He i. We also detect molecular absorption from CH+ . Conclusions. We interpret the unusual properties of these spectra as reflecting the presence of three distinct regions along the line-of-sight: the excited He i absorption originates from an H ii-region, whereas the Balmer absorption must originate from an associated photodissociation region. The strong metal line and molecular absorption and the dust extinction must originate from a third, cooler region along the line-of-sight. The presence of (at least) three separate regions is reflected in the fact that the different absorption components have different velocities relative to the systemic redshift of the host galaxy.

Long-duration gamma-ray bursts (LGRBs) are produced by the collapse of very massive stars. Due to the short lifetime of their progenitors, LGRBs pinpoint star-forming galaxies. We present here a multi-band search for the host galaxy of the long dark GRB 050219A within the enhanced Swift/XRT error circle. We used spectroscopic observations acquired with VLT/X-shooter to determine the redshift and star-formation rate of the putative host galaxy. We compared the results with the optical/IR spectral energy distribution obtained with different facilities. Surprisingly, the host galaxy is a old and quiescent early-type galaxy at z = 0.211 characterised by an unprecedentedly low specific star-formation rate. It is the first LGRB host to be also an early-type post-starburst galaxy. This is further evidence that GRBs can explode in all kind of galaxies, with the only requirement being an episode of star-formation.

We use a parent sample of 118 gamma-ray burst (GRB) afterglows, with known redshift and host galaxy extinction, to separate afterglows with and without signatures of dominant reverse-shock emission and to determine which physical conditions lead to a prominent reverse-shock emission. We identify 10 GRBs with reverse shock signatures - GRBs 990123, 021004, 021211, 060908, 061126, 080319B, 081007, 090102, 090424 and 130427A. By modeling their optical afterglows with reverse and forward shock analytic light curves and using Monte Carlo simulations, we estimate the parameter space of the physical quantities describing the ejecta and circumburst medium. We find that physical properties cover a wide parameter space and do not seem to cluster around any preferential values. Comparing the rest-frame optical, X-ray and high-energy properties of the larger sample of non-RS-dominated GRBs, we show that the early-time ($<$ 1ks) optical spectral luminosity, X-ray afterglow luminosity and $\gamma$-ray energy output of our reverse-shock dominated sample do not differ significantly from the general population at early times. However, the GRBs with dominant reverse shock emission have fainter than average optical forward-shock emission at late time ($>$ 10 ks). We find that GRBs with an identifiable reverse shock component show high magnetization parameter $R_{\mathrm{B}} = \varepsilon_{\rm B,r}/\varepsilon_{\rm B,f} \sim 2 - 10^4$. Our results are in agreement with the mildly magnetized baryonic jet model of GRBs.

Long-duration Gamma-Ray Bursts (GRBs) are an extremely rare outcome of the collapse of massive stars, and are typically found in the distant Universe. Because of its intrinsic luminosity ($L\sim 3 \times 10^{53}$ erg s$^{-1}$) and its relative proximity ($z=0.34$), GRB 130427A was a unique event that reached the highest fluence observed in the gamma-ray band. Here we present a comprehensive multiwavelength view of GRB 130427A with Swift, the 2-m Liverpool and Faulkes telescopes and by other ground-based facilities, highlighting the evolution of the burst emission from the prompt to the afterglow phase. The properties of GRB 130427A are similar to those of the most luminous, high-redshift GRBs, suggesting that a common central engine is responsible for producing GRBs in both the contemporary and the early Universe and over the full range of GRB isotropic energies.

We use high--quality, multi-band observations of Swift GRB120404A, from gamma-ray to radio frequencies, together with the new hydrodynamics code of van Eerten et al. (2012) to test the standard synchrotron shock model. The evolution of the radio and optical afterglow, with its prominent optical rebrightening at t_rest 260-2600 s, is remarkably well modelled by a decelerating jet viewed close to the jet edge, combined with some early re-energization of the shock. We thus constrain the geometry of the jet with half-opening and viewing angles of 23 and 21 deg respectively and suggest that wide jets viewed off-axis are more common in GRBs than previously thought. We also derive the fireball microphysics parameters epsilon_B=2.4e-4 and epsilon_e=9.3e-2 and a circumburst density of n=240 cm^-3. The ability to self-consistently model the microphysics parameters and jet geometry in this way offers an alternative to trying to identify elusive canonical jet breaks at late times. The mismatch between the observed and model-predicted X-ray fluxes is explained by the local rather than the global cooling approximation in the synchrotron radiation model, constraining the microphysics of particle acceleration taking place in a relativistic shock and, in turn, emphasising the need for a more realistic treatment of cooling in future developments of theoretical models. Finally, our interpretation of the optical peak as due to the passage of the forward shock synchrotron frequency highlights the importance of high quality multi-band data to prevent some optical peaks from being erroneously attributed to the onset of fireball deceleration.

We present a broadband study of gamma-ray burst (GRB) 091024A within the context of other ultra-long-duration GRBs. An unusually long burst detected by Konus-Wind, Swift, and Fermi, GRB 091024A has prompt emission episodes covering ~1300 s, accompanied by bright and highly structured optical emission captured by various rapid-response facilities, including the 2-m autonomous robotic Faulkes North and Liverpool Telescopes, KAIT, S-LOTIS, and SRO. We also observed the burst with 8- and 10-m class telescopes and determine the redshift to be z = 1.0924 \pm 0.0004. We find no correlation between the optical and gamma-ray peaks and interpret the optical light curve as being of external origin, caused by the reverse and forward shock of a highly magnetized jet (R_B ~ 100-200). Low-level emission is detected throughout the near-background quiescent period between the first two emission episodes of the Konus-Wind data, suggesting continued central-engine activity; we discuss the implications of this ongoing emission and its impact on the afterglow evolution and predictions. We summarize the varied sample of historical GRBs with exceptionally long durations in gamma-rays (>~ 1000 s) and discuss the likelihood of these events being from a separate population; we suggest ultra-long GRBs represent the tail of the duration distribution of the long GRB population.

We present a multi-wavelength analysis of gamma-ray burst GRB 090727, for which optical emission was detected during the prompt gamma-ray emission by the 2-m autonomous robotic Liverpool Telescope and subsequently monitored for a further two days with the Liverpool and Faulkes telescopes. Within the context of the standard fireball model, we rule out a reverse shock origin for the early time optical emission in GRB 090727 and instead conclude that the early time optical flash likely corresponds to emission from an internal dissipation processes. Putting GRB 090727 into a broader observational and theoretical context, we build a sample of 36 gamma-ray bursts (GRBs) with contemporaneous early time optical and gamma-ray detections. From these GRBs, we extract a sub-sample of 18 GRBs, which show optical peaks during prompt gamma-ray emission, and perform detailed temporal and spectral analysis in gamma-ray, X-ray, and optical bands. We find that in most cases early time optical emission shows sharp and steep behaviour, and notice a rich diversity of spectral properties. Using a simple internal shock dissipation model, we show that the emission during prompt GRB phase can occur at very different frequencies via synchrotron radiation. Based on the results obtained from observations and simulation, we conclude that the standard external shock interpretation for early time optical emission is disfavoured in most cases due to sharp peaks ($\Delta t/t < 1$) and steep rise/decay indices, and that internal dissipation can explain the properties of GRBs with optical peaks during gamma-ray emission.

Because of their cosmological origin, gamma-ray burst (GRB) optical afterglows are attenuated when they pass intergalactic absorbers in the GRB line-of-sight. Without the knowledge of the number of absorbers and their physical properties, the effect of absorption on the observed magnitudes can not be determined precisely. Different methods have been applied in order to correct for this effect statistically, either using semi-analytical calculations or numerical simulations. We follow these works and present the expected magnitude corrections as a function of redshift for a set of filters most commonly used in the scientific community. The results are publically available on the web (http://igmac.fmf.uni-lj.si).

In this paper we present a sample of 10 short gamma-ray bursts (GRBs) with a robust redshift determination, discovered by the Swift satellite up to January 2011. We measure their X-ray absorbing column densities and collect data on the host galaxy offsets. We find evidence for intrinsic absorption and no correlation between the intrinsic absorbing column density and the projected offset of the GRB from its host galaxy center. We find that the properties in the gamma regime (T90, fluence and 1-s peak photon flux) of short GRBs with "bright" and "faint" X-ray afterglow likely disfavour different prompt emission mechanisms. Host galaxy offset and GRB duration (T90) do not correlate. Instead, there is a hint of anti-correlation between the effective radius normalised host galaxy offset and T90. Finally, we examine the properties of short GRBs with short-lived and long-lived X-ray afterglows, finding that some short GRBs with short-lived X-ray afterglows have their optical afterglow detected. In light of this, the X-ray afterglow duration does not seem to be an unique indicator of a specific progenitor and/or environment for short GRBs.

With the launch of Gaia satellite, detection of many different types of transient sources will be possible, one of them being optical afterglows of gamma-ray bursts (GRBs). Using the knowledge of the satellites dynamics and properties of GRB optical afterglows we performed a simulation in order to estimate an average GRB detection rate with Gaia. Here we present the simulation results for two types of GRB optical afterglows, differing in the observer's line-of-sight compared to a GRB jet axis: regular (on-axis) and orphan afterglows. Results show that for on-axis GRBs, less than 10 detections in five years of foreseen Gaia operational time are expected. The orphan afterglows simulation results are more promising, giving a more optimistic number of several tens of detections in five years.