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 duration of a gamma-ray burst (GRB) is a key indicator of its physics origin, with long bursts perhaps associated with the collapse of massive stars and short bursts with mergers of neutron stars.However, there is substantial overlap in the properties of both short and long GRBs and neither duration nor any other parameter so far considered completely separates the two groups. Here we unambiguously classify every GRB using a machine-learning, dimensionality-reduction algorithm, t-distributed stochastic neighborhood embedding (t-SNE), providing a catalog separating all Swift GRBs into two groups. Although the classification takes place only using prompt emission light curves,every burst with an associated supernova is found in the longer group and bursts with kilonovae in the short, suggesting along with the duration distributions that these two groups are truly long and short GRBs. Two bursts with a clear absence of a supernova belong to the longer class, indicating that these might have been direct-collapse black holes, a proposed phenomenon that may occur in the deaths of more massive stars.

Multi-messenger astronomy received a great boost following the discovery of kilonova AT2017gfo, the optical counterpart of the gravitational wave source GW170817 associated with the short gamma-ray burst GRB 170817A. AT2017gfo was the first kilonova that could be extensively monitored in time both photometrically and spectroscopically. Previously, only few candidates have been observed against the glare of short GRB afterglows. In this work, we aim to search the fingerprints of AT2017gfo-like kilonova emissions in the optical/NIR light curves of 39 short GRBs with known redshift. For the first time, our results allow us to study separately the range of luminosity of the blue and red components of AT2017gfo-like kilonovae in short GRBs. In particular, the red component is similar in luminosity to AT2017gfo, while the blue kilonova can be more than 10 times brighter. Finally, we find further evidence to support all the claimed kilonova detections and we exclude an AT2017gfo-like kilonova in GRBs 050509B and 061201.

We report two secure ($z=3.775, 4.012$) and one tentative ($z\approx3.767$) spectroscopic confirmations of massive and quiescent galaxies through $K$-band observations with Keck/MOSFIRE and VLT/X-Shooter. The stellar continuum emission, the absence of strong nebular emission lines and the lack of significant far-infrared detections confirm the passive nature of these objects, disfavoring the alternative solution of low-redshift dusty star-forming interlopers. We derive stellar masses of $\mathrm{log}(M_{\star}/M_\odot)\sim11$ and ongoing star formation rates placing these galaxies $\gtrsim 1-2$ dex below the main sequence at their redshifts. The adopted parametrization of the star formation history suggests that these sources experienced a strong ($\langle \rm SFR \rangle \sim 1200-3500\,M_\odot\,\mathrm{yr}^{-1}$) and short ($\sim 50$ Myr) burst of star formation, peaking $\sim 150-500$ Myr before the time of observation, all properties reminiscent of the characteristics of sub-millimeter galaxies (SMGs) at $z>4$. We investigate this connection by comparing the comoving number densities and the properties of these two populations. We find a fair agreement only with the deepest sub-mm surveys detecting not only the most extreme starbursts, but also more normal galaxies. We support these findings by further exploring the Illustris-TNG cosmological simulation, retrieving populations of both fully quenched massive galaxies at $z\sim3-4$ and SMGs at $z\sim4-5$, with number densities and properties in agreement with the observations at $z\sim3$, but in increasing tension at higher redshift. Nevertheless, as suggested by the observations, not all the progenitors of quiescent galaxies at these redshifts shine as bright SMGs in their past and, similarly, not all bright SMGs quench by $z\sim3$, both fractions depending on the threshold assumed to define the SMGs themselves.

We present a detailed analysis of a large sample of spectroscopically confirmed ultra-massive quiescent galaxies (${\rm{log}}(M_{\ast}/M_{\odot})\sim11.5$) at $z\gtrsim2$. This sample comprises 15 galaxies selected in the COSMOS and UDS fields by their bright K-band magnitudes and followed up with VLT/X-shooter spectroscopy and HST/WFC3 $H_{F160W}$ imaging. These observations allow us to unambiguously confirm their redshifts ascertain their quiescent nature and stellar ages, and to reliably assess their internal kinematics and effective radii. We find that these galaxies are compact, consistent with the high mass end of the mass-size relation for quiescent galaxies at $z=2$. Moreover, the distribution of the measured stellar velocity dispersions of the sample is consistent with the most massive local early-type galaxies from the MASSIVE Survey showing that evolution in these galaxies, is dominated by changes in size. The HST images reveal, as surprisingly high, that $40\ \%$ of the sample have tidal features suggestive of mergers and companions in close proximity, including three galaxies experiencing ongoing major mergers. The absence of velocity dispersion evolution from $z=2$ to $0$, coupled with a doubling of the stellar mass, with a factor of four size increase and the observed disturbed stellar morphologies support dry minor mergers as the primary drivers of the evolution of the massive quiescent galaxies over the last 10 billion years.

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.

Half of all the elements in the universe heavier than iron were created by rapid neutron capture. The theory for this astrophysical `$r$-process' was worked out six decades ago and requires an enormous neutron flux to make the bulk of these elements. Where this happens is still debated. A key piece of missing evidence is the identification of freshly-synthesised $r$-process elements in an astrophysical site. Current models and circumstantial evidence point to neutron star mergers as a probable $r$-process site, with the optical/infrared `kilonova' emerging in the days after the merger a likely place to detect the spectral signatures of newly-created neutron-capture elements. The kilonova, AT2017gfo, emerging from the gravitational-wave--discovered neutron star merger, GW170817, was the first kilonova where detailed spectra were recorded. When these spectra were first reported it was argued that they were broadly consonant with an outflow of radioactive heavy elements, however, there was no robust identification of any element. Here we report the identification of the neutron-capture element strontium in a re-analysis of these spectra. The detection of a neutron-capture element associated with the collision of two extreme-density stars establishes the origin of $r$-process elements in neutron star mergers, and demonstrates that neutron stars contain neutron-rich matter.

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].

Pulsars are well studied all over the electromagnetic spectrum, and the Crab pulsar may be the most studied object in the sky. Nevertheless, a high-quality optical to near-infrared spectrum of the Crab or any other pulsar has not been published to date. Obtaining a properly flux-calibrated spectrum enables us to measure the spectral index of the pulsar emission, without many of the caveats from previous studies. This was the main aim of this project, but we could also detect absorption and emission features from the pulsar and nebula over an unprecedentedly wide wavelength range. A spectrum was obtained with the X-shooter spectrograph on the Very Large Telescope. Particular care was given to the flux-calibration of these data. A high signal-to-noise spectrum of the Crab pulsar was obtained from 300 to 2400nm. The spectral index fitted to this spectrum is flat with alpha_nu=0.16 +- 0.07. For the emission lines we measure a maximum velocity of 1600 km/s, whereas the absorption lines from the material between us and the pulsar is unresolved at the 50 km/s resolution. A number of Diffuse Interstellar Bands and a few near-IR emission lines that have previously not been reported from the Crab are highlighted.

We report observations and analysis of the nearby gamma-ray burst GRB\,161219B (redshift $z=0.1475$) and the associated Type Ic supernova (SN) 2016jca. GRB\,161219B had an isotropic gamma-ray energy of $\sim 1.6 \times 10^{50}$\,erg. Its afterglow is likely refreshed at an epoch preceding the first photometric points (0.6\,d), which slows down the decay rates. Combined analysis of the SN light curve and multiwavelength observations of the afterglow suggest that the GRB jet was broad during the afterglow phase (full opening angle $\sim 42^\circ \pm 3^\circ$). Our spectral series shows broad absorption lines typical of GRB supernovae (SNe), which testify to the presence of material with velocities up to $\sim 0.25$c. The spectrum at 3.73\,d allows for the very early identification of a SN associated with a GRB. Reproducing it requires a large photospheric velocity ($35,000 \pm 7000$\,\kms). The kinetic energy of the SN is estimated through models to be \KE $\approx 4 \times 10^{52}$\,erg in spherical symmetry. The ejected mass in the explosion was \Mej $\approx 6.5 \pm 1.5$\,\Msun, much less than that of other GRB-SNe, demonstrating diversity among these events. The total amount of \Nifs\ in the explosion was $0.27 \pm 0.05$\,\Msun. The observed spectra require the presence of freshly synthesised \Nifs\ at the highest velocities, at least 3 times more than a standard GRB-SN. We also find evidence for a decreasing \Nifs\ abundance as a function of decreasing velocity. This suggests that SN\,2016jca was a highly aspherical explosion viewed close to on-axis, powered by a compact remnant. Applying a typical correction for asymmetry, the energy of SN\,2016jca was $\sim$ (1--3) $\times 10^{52}$\,erg, confirming that most of the energy produced by GRB-SNe goes into the kinetic energy of the SN ejecta.

Large surveys of galaxy clusters with the Hubble and Spitzer Space Telescopes, including CLASH and the Frontier Fields, have demonstrated the power of strong gravitational lensing to efficiently deliver large samples of high-redshift galaxies. We extend this strategy through a wider, shallower survey named RELICS, the Reionization Lensing Cluster Survey. This survey, described here, was designed primarily to deliver the best and brightest high-redshift candidates from the first billion years after the Big Bang. RELICS observed 41 massive galaxy clusters with Hubble and Spitzer at 0.4-1.7um and 3.0-5.0um, respectively. We selected 21 clusters based on Planck PSZ2 mass estimates and the other 20 based on observed or inferred lensing strength. Our 188-orbit Hubble Treasury Program obtained the first high-resolution near-infrared images of these clusters to efficiently search for lensed high-redshift galaxies. We observed 46 WFC3/IR pointings (~200 arcmin^2) with two orbits divided among four filters (F105W, F125W, F140W, and F160W) and ACS imaging as needed to achieve single-orbit depth in each of three filters (F435W, F606W, and F814W). As previously reported by Salmon et al., we discovered 322 z ~ 6 - 10 candidates, including the brightest known at z ~ 6, and the most distant spatially-resolved lensed arc known at z ~ 10. Spitzer IRAC imaging (945 hours awarded, plus 100 archival) has crucially enabled us to distinguish z ~ 10 candidates from z ~ 2 interlopers. For each cluster, two HST observing epochs were staggered by about a month, enabling us to discover 11 supernovae, including 3 lensed supernovae, which we followed up with 20 orbits from our program. We delivered reduced HST images and catalogs of all clusters to the public via MAST and reduced Spitzer images via IRSA. We have also begun delivering lens models of all clusters, to be completed before the JWST GO call for proposals.

Individual highly magnified stars have been recently discovered at lookback times of more than half the age of the Universe, in lensed galaxies that straddle the critical curves of massive galaxy clusters. Having confirmed their detectability, it is now important to carry out systematic searches for them in order to establish their frequency, and in turn learn about the statistical properties of high-redshift stars and of the granularity of matter in the foreground deflector. Here we report the discovery of a highly magnified star at redshift $z = 0.94$ in a strongly lensed arc behind a Hubble Frontier Field galaxy cluster, MACS J0416.1-2403, discovered as part of a systematic archival search. The bright transient (dubbed "Warhol") was discovered in Hubble Space Telescope data taken on 2014 September 15 and 16. This single image faded over a period of two weeks, and observations taken on 2014 September 1 show that the duration of the microlensing event was at most four weeks in total. The light curve may also exhibit slow changes over a period of years consistent with the level of microlensing expected from stars responsible for the intracluster light (ICL) of the cluster. Optical and infrared observations taken near peak brightness can be fit by a stellar spectrum with moderate host-galaxy extinction. A blue supergiant matches the measured spectral energy distribution near peak, implying a temporary magnification of at least several thousand. While the spectrum of an O-type star would also fit the transient's spectral energy distribution, extremely luminous O-type stars are much less common than blue supergiants. The short timescale of the event and the estimated effective temperature indicate that the lensed source is an extremely magnified star.

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.

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.

Super-luminous supernovae (SLSNe) are rare events defined as being significantly more luminous than normal terminal stellar explosions. The source of the extra powering needed to achieve such luminosities is still unclear. Discoveries in the local Universe (i.e. $z<0.1$) are scarce, but afford dense multi-wavelength observations. Additional low-redshift objects are therefore extremely valuable. We present early-time observations of the type I SLSN ASASSN-18km/SN~2018bsz. These data are used to characterise the event and compare to literature SLSNe and spectral models. Host galaxy properties are also analysed. Optical and near-IR photometry and spectroscopy were analysed. Early-time ATLAS photometry was used to constrain the rising light curve. We identified a number of spectral features in optical-wavelength spectra and tracked their time evolution. Finally, we used archival host galaxy photometry together with HII region spectra to constrain the host environment. ASASSN-18km/SN~2018bsz is found to be a type I SLSN in a galaxy at a redshift of 0.0267 (111 Mpc), making it the lowest-redshift event discovered to date. Strong CII lines are identified in the spectra. Spectral models produced by exploding a Wolf-Rayet progenitor and injecting a magnetar power source are shown to be qualitatively similar to ASASSN-18km/SN~2018bsz, contrary to most SLSNe-I that display weak/non-existent CII lines. ASASSN-18km/SN~2018bsz displays a long, slowly rising, red 'plateau' of $>$26 days, before a steeper, faster rise to maximum. The host has an absolute magnitude of --19.8 mag ($r$), a mass of M$_{*}$ = 1.5$^{+0.08}_{-0.33}$ $\times$10$^{9}$ M$_{\odot}$ , and a star formation rate of = 0.50$^{+2.22}_{-0.19}$ M$_{\odot}$ yr$^{-1}$. A nearby HII region has an oxygen abundance (O3N2) of 8.31$\pm$0.01 dex.

Long gamma-ray bursts are produced during the dramatic deaths of massive stars with very short lifetimes, meaning that they explode close to the birth place of their progenitors. During a short period they become the most luminous objects observable in the Universe, being perfect beacons to study high-redshift star-forming regions. To use the afterglow of GRB 161023A at a redshift $z=2.710$ as a background source to study the environment of the explosion and the intervening systems along its line-of-sight. r the first time, we complement UV/Optical/NIR spectroscopy with millimetre spectroscopy using ALMA, which allows us to probe the molecular content of the host galaxy. The X-shooter spectrum shows a plethora of absorption features including fine-structure and metastable transitions of Fe, Ni, Si, C and O. We present photometry ranging from 43 s to over 500 days after the burst. We infer a host-galaxy metallicity of [Zn/H] $=-1.11\pm0.07$, which corrected for dust depletion results in [X/H] $=-0.94\pm0.08$. We do not detect molecular features in the ALMA data, but we derive limits on the molecular content of $log(N_{CO}/cm^{-2})<15.7$ and $log(N_{HCO+}/cm^{-2})<13.2$, which are consistent with those that we obtain from the optical spectra, $log(N_{H_2}/cm^{-2})<15.2$ and $log(N_{CO}/cm^{-2})<14.5$. Within the host galaxy we detect three velocity systems through UV/Optical/NIR absorption spectroscopy, all with levels that were excited by the GRB afterglow. We determine the distance from these systems to the GRB to be in the range between 0.7 and 1.0 kpc. The sight-line to GRB 161023A shows 9 independent intervening systems, most of them with multiple components. (Abridged)

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

Galaxy-cluster gravitational lenses can magnify background galaxies by a total factor of up to ~50. Here we report an image of an individual star at redshift z=1.49 (dubbed "MACS J1149 Lensed Star 1 (LS1)") magnified by >2000. A separate image, detected briefly 0.26 arcseconds from LS1, is likely a counterimage of the first star demagnified for multiple years by a >~3 solar-mass object in the cluster. For reasonable assumptions about the lensing system, microlensing fluctuations in the stars' light curves can yield evidence about the mass function of intracluster stars and compact objects, including binary fractions and specific stellar evolution and supernova models. Dark-matter subhalos or massive compact objects may help to account for the two images' long-term brightness ratio.

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 here present the spectroscopic follow-up observations with VLT/X-shooter of the Swift long-duration gamma-ray burst GRB 160804A at z = 0.737. Typically, GRBs are found in low-mass, metal-poor galaxies which constitute the sub-luminous population of star-forming galaxies. For the host galaxy of the GRB presented here we derive a stellar mass of $\log(M_*/M_{\odot}) = 9.80\pm 0.07$, a roughly solar metallicity (12+log(O/H) = $8.74\pm 0.12$) based on emission line diagnostics, and an infrared luminosity of $M_{3.6/(1+z)} = -21.94$ mag, but find it to be dust-poor ($E(B-V) < 0.05$ mag). This establishes the galaxy hosting GRB 160804A as one of the most luminous, massive and metal-rich GRB hosts at z < 1.5. Furthermore, the gas-phase metallicity is found to be representative of the physical conditions of the gas close to the explosion site of the burst. The high metallicity of the host galaxy is also observed in absorption, where we detect several strong FeII transitions as well as MgII and MgI. While host galaxy absorption features are common in GRB afterglow spectra, we detect absorption from strong metal lines directly in the host continuum (at a time when the afterglow was contributing to < 15%). Finally, we discuss the possibility that the geometry and state of the absorbing and emitting gas is indicative of a galactic scale outflow expelled at the final stage of two merging galaxies.

The merger of two neutron stars is predicted to give rise to three major detectable phenomena: a short burst of gamma-rays, a gravitational wave signal, and a transient optical/near-infrared source powered by the synthesis of large amounts of very heavy elements via rapid neutron capture (the r-process). Such transients, named "macronovae" or "kilonovae", are believed to be centres of production of rare elements such as gold and platinum. The most compelling evidence so far for a kilonova was a very faint near-infrared rebrightening in the afterglow of a short gamma-ray burst at z = 0.356, although findings indicating bluer events have been reported. Here we report the spectral identification and describe the physical properties of a bright kilonova associated with the gravitational wave source GW 170817 and gamma-ray burst GRB 170817A associated with a galaxy at a distance of 40 Mpc from Earth. Using a series of spectra from ground-based observatories covering the wavelength range from the ultraviolet to the near-infrared, we find that the kilonova is characterized by rapidly expanding ejecta with spectral features similar to those predicted by current models. The ejecta is optically thick early on, with a velocity of about 0.2 times light speed, and reaches a radius of about 50 astronomical units in only 1.5 days. As the ejecta expands, broad absorption-like lines appear on the spectral continuum indicating atomic species produced by nucleosynthesis that occurs in the post-merger fast-moving dynamical ejecta and in two slower (0.05 times light speed) wind regions. Comparison with spectral models suggests that the merger ejected 0.03-0.05 solar masses of material, including high-opacity lanthanides.

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.

The spectacular success of type Ia supernovae (SNe Ia) in SN-cosmology is based on the assumption that their photometric and spectroscopic properties are invariant with redshift. However, this fundamental assumption needs to be tested with observations of high-z SNe Ia. To date, the majority of SNe Ia observed at moderate to large redshifts (0.4 < z < 1.0) are faint, and the resultant analyses are based on observations with modest signal-to-noise ratios that impart a degree of ambiguity in their determined properties. In rare cases however, the Universe offers a helping hand: to date a few SNe Ia have been observed that have had their luminosities magnified by intervening galaxies and galaxy clusters acting as gravitational lenses. In this paper we present long-slit spectroscopy of the lensed SNe Ia 2016geu, which occurred at a redshift of z=0.409, and was magnified by a factor of ~55 by a galaxy located at z=0.216. We compared our spectra, which were obtained a couple weeks to a couple months past peak light, with the spectroscopic properties of well-observed, nearby SNe Ia, finding that SN 2016geu's properties are commensurate with those of SNe Ia in the local universe. Based primarily on the velocity and strength of the Si II 6355 absorption feature, we find that SN 2016geu can be classified as a high-velocity, high-velocity gradient and "core-normal" SN Ia. The strength of various features (measured though their pseudo-equivalent widths) argue against SN 2016geu being a faint, broad-lined, cool or shallow-silicon SN Ia. We conclude that the spectroscopic properties of SN 2016geu imply that it is a normal SN Ia, and when taking previous results by other authors into consideration, there is very little, if any, evolution in the observational properties of SNe Ia up to z~0.4. [Abridged]

A massive galaxy cluster can serve as a magnifying glass for distant stellar populations, with strong gravitational lensing exposing details in the lensed background galaxies that would otherwise be undetectable. The MACS J0416.1-2403 cluster (hereafter MACS0416) is one of the most efficient lenses in the sky, and in 2014 it was observed with high-cadence imaging from the Hubble Space Telescope (HST). Here we describe two unusual transient events that appeared behind MACS0416 in a strongly lensed galaxy at redshift $z = 1.0054 \pm 0.0002$. These transients---designated HFF14Spo-NW and HFF14Spo-SE and collectively nicknamed "Spock"---were faster and fainter than any supernova (SN), but significantly more luminous than a classical nova. They reached peak luminosities of $\sim10^{41}$ erg s$^{-1}$ ($M_{\rm AB} < -14$ mag) in 5 rest-frame days, then faded below detectability in roughly the same time span. Models of the cluster lens suggest that these events may be spatially coincident at the source plane, but are most likely not temporally coincident. We find that HFF14Spo can be explained as a luminous blue variable (LBV), a recurrent nova (RN), or a pair of stellar microlensing events. To distinguish between these hypotheses will require a clarification of the positions of nearby critical curves, along with high-cadence monitoring of the field that could detect new transient episodes in the host galaxy.

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 the spectroscopic and photometric late-time follow-up of the host galaxy of the long-duration Swift gamma-ray burst GRB 140506A at z = 0.889. The optical and near-infrared afterglow of this GRB had a peculiar spectral energy distribution (SED) with a strong flux-drop at 8000 {\AA} (4000 {\AA} rest-frame) suggesting an unusually steep extinction curve. By analyzing the contribution and physical properties of the host galaxy, we here aim at providing additional information on the properties and origin of this steep, non-standard extinction. We find that the strong flux-drop in the GRB afterglow spectrum at < 8000 {\AA} and rise at < 4000 {\AA} is well explained by the combination of a steep extinction curve along the GRB line of sight and contamination by the host galaxy light so that the scenario with an extreme 2175 {\AA} extinction bump can be excluded. We localise the GRB to be at a projected distance of approximately 4 kpc from the centre of the host galaxy. Based on emission-line diagnostics of the four detected nebular lines, Halpha, Hbeta, [O II] and [O III], we find the host to be a modestly star forming (SFR = 1.34 +/- 0.04 Msun yr^-1) and relatively metal poor (Z = 0.35^{+0.15}_{-0.11} Zsun) galaxy with a large dust content, characterized by a measured visual attenuation of A_V = 1.74 +/- 0.41 mag, thus unexceptional in all its physical properties. We model the extinction curve of the host-corrected afterglow and show that the standard dust properties causing the reddening seen in the Local Group are inadequate in describing the steep drop. We conclude that the steep extinction curve seen in the afterglow towards the GRB is of exotic origin, is sightline-dependent only and thus solely a consequence of the circumburst environment.

In Hubble Space Telescope (HST) imaging taken on 10 November 2014, four images of supernova (SN) "Refsdal" (redshift z=1.49) appeared in an Einstein-cross--like configuration (images S1-S4) around an early-type galaxy in the cluster MACS J1149.5+2223 (z=0.54). Almost all lens models of the cluster have predicted that the SN should reappear within a year in a second host-galaxy image created by the cluster's potential. In HST observations taken on 11 December 2015, we find a new source at the predicted position of the new image of SN Refsdal approximately 8" from the previous images S1-S4. This marks the first time the appearance of a SN at a particular time and location in the sky was successfully predicted in advance! We use these data and the light curve from the first four observed images of SN Refsdal to place constraints on the relative time delay and magnification of the new image (SX), compared to images S1-S4. This enables us, for the first time, to test "blind" lens model predictions of both magnifications and time delays for a lensed SN. We find that the timing and brightness of the new image are consistent with the blind predictions of a fraction of the models. The reappearance illustrates the discriminatory power of this blind test and its utility to uncover sources of systematic uncertainty. From planned HST photometry, we expect to reach a precision of 1-2% on the time delay between S1-S4 and SX.

We present the first year of Hubble Space Telescope imaging of the unique supernova (SN) 'Refsdal', a gravitationally lensed SN at z=1.488$\pm$0.001 with multiple images behind the galaxy cluster MACS J1149.6+2223. The first four observed images of SN Refsdal (images S1-S4) exhibited a slow rise (over ~150 days) to reach a broad peak brightness around 20 April, 2015. Using a set of light curve templates constructed from SN 1987A-like peculiar Type II SNe, we measure time delays for the four images relative to S1 of 4$\pm$4 (for S2), 2$\pm$5 (S3), and 24$\pm$7 days (S4). The measured magnification ratios relative to S1 are 1.15$\pm$0.05 (S2), 1.01$\pm$0.04 (S3), and 0.34$\pm$0.02 (S4). None of the template light curves fully captures the photometric behavior of SN Refsdal, so we also derive complementary measurements for these parameters using polynomials to represent the intrinsic light curve shape. These more flexible fits deliver fully consistent time delays of 7$\pm$2 (S2), 0.6$\pm$3 (S3), and 27$\pm$8 days (S4). The lensing magnification ratios are similarly consistent, measured as 1.17$\pm$0.02 (S2), 1.00$\pm$0.01 (S3), and 0.38$\pm$0.02 (S4). We compare these measurements against published predictions from lens models, and find that the majority of model predictions are in very good agreement with our measurements. Finally, we discuss avenues for future improvement of time delay measurements -- both for SN Refsdal and for other strongly lensed SNe yet to come.

We have acquired Hubble Space Telescope (HST) and Very Large Telescope near-infrared spectra and images of supernova (SN) Refsdal after its discovery as an Einstein cross in Fall 2014. The HST light curve of SN Refsdal matches the distinctive, slowly rising light curves of SN 1987A-like supernovae (SNe), and we find strong evidence for a broad H-alpha P-Cygni profile in the HST grism spectrum at the redshift (z = 1.49) of the spiral host galaxy. SNe IIn, powered by circumstellar interaction, could provide a good match to the light curve of SN Refsdal, but the spectrum of a SN IIn would not show broad and strong H-alpha absorption. From the grism spectrum, we measure an H-alpha expansion velocity consistent with those of SN 1987A-like SNe at a similar phase. The luminosity, evolution, and Gaussian profile of the H-alpha emission of the WFC3 and X-shooter spectra, separated by ~2.5 months in the rest frame, provide additional evidence that supports the SN 1987A-like classification. In comparison with other examples of SN 1987A-like SNe, SN Refsdal has a blue B-V color and a high luminosity for the assumed range of potential magnifications. If SN Refsdal can be modeled as a scaled version of SN 1987A, we estimate it would have an ejecta mass of 20+-5 solar masses. The evolution of the light curve at late times will provide additional evidence about the potential existence of any substantial circumstellar material (CSM). Using MOSFIRE and X-shooter spectra, we estimate a subsolar host-galaxy metallicity (8.3+-0.1 dex and <8.4 dex, respectively) near the explosion site.

Quasi-stellar object (QSO) spectral templates are important both to QSO physics and for investigations that use QSOs as probes of intervening gas and dust. However, combinations of various QSO samples obtained at different times and with different instruments so as to expand a composite and to cover a wider rest frame wavelength region may create systematic effects, and the contribution from QSO hosts may contaminate the composite. We have constructed a composite spectrum from luminous blue QSOs at 1 < z < 2.1 selected from the Sloan Digital Sky Survey (SDSS). The observations with X-shooter simultaneously cover ultraviolet (UV) to near- infrared (NIR) light, which ensures that the composite spectrum covers the full rest-frame range from Ly$\beta$ to 11350 $\AA$ without any significant host contamination. Assuming a power-law continuum for the composite we find a spectral slope of $\alpha_\lambda$ = 1.70+/-0.01, which is steeper than previously found in the literature. We attribute the differences to our broader spectral wavelength coverage, which allows us to effectively avoid fitting any regions that are affected either by strong QSO emissions lines (e.g., Balmer lines and complex [Fe II] blends) or by intrinsic host galaxy emission. Finally, we demonstrate the application of the QSO composite spectrum for evaluating the reddening in other QSOs.

[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.

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.