We present the discovery and subarcsecond localization of a new Fast Radio Burst with the Karl G. Jansky Very Large Array and realfast search system. The FRB was discovered on 2019 June 14 with a dispersion measure of 959 pc/cm3. This is the highest DM of any localized FRB and its measured burst fluence of 0.6 Jy ms is less than nearly all other FRBs. The source is not detected to repeat in 15 hours of VLA observing and 153 hours of CHIME/FRB observing. We describe a suite of statistical and data quality tests we used to verify the significance of the event and its localization precision. Follow-up optical/infrared photometry with Keck and Gemini associate the FRB to a pair of galaxies with $\rm{r}\sim23$ mag. The false-alarm rate for radio transients of this significance that are associated with a host galaxy is roughly $3\times10^{-4}\ \rm{hr}^{-1}$. The two putative host galaxies have similar photometric redshifts of $z_{\rm{phot}}\sim0.6$, but different colors and stellar masses. Comparing the host distance to that implied by the dispersion measure suggests a modest (~ 50 pc/cm3) electron column density associated with the FRB environment or host galaxy/galaxies.

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.

We present a high-resolution analysis of the host galaxy of fast radio burst FRB 190608, an SBc galaxy at $z=0.11778$ (hereafter HG 190608), to dissect its local environment and its contributions to the FRB properties. Our Hubble Space Telescope WFC3/UVIS image reveals that the sub-arcsecond localization of FRB 190608 is coincident with a knot of star-formation ($\Sigma_{SFR} = 1.2 \times 10^{-2}~ M_{\odot} \, kpc^{-2}$) in one of the prominent spiral arms of HG 190608. This is confirmed by H$\beta$ emission present in our Keck/KCWI integral field spectrum of the galaxy with a surface brightness of $\mu_{H\beta} = (3.35\pm0.18)\times10^{-17}\;erg\;s^{-1}\;cm^{-2}\;arcsec^{-2}$. We infer an extinction-corrected H$\alpha$ surface brightness and compute a dispersion measure from the interstellar medium of HG 190608 of ${DM}_{Host,ISM} = 82 \pm 35~ pc \, cm^{-3}$. The galaxy rotates with a circular velocity $v_{circ} = 141 \pm 8~ km \, s^{-1}$ at an inclination $i_{gas} = 37 \pm 3^\circ$, giving a dynamical mass $M_{halo}^{dyn} \approx 10^{11.96 \pm 0.08}~ M_{\odot}$. A surface photometric analysis of the galaxy using FORS2 imaging suggests a stellar disk inclination of $i_{stellar} = 26 \pm 3^\circ$. The dynamical mass estimate implies a halo contribution to the dispersion measure of ${DM}_{Host,Halo} = 55 \pm 25\; pc \, cm^{-3}$ subject to assumptions on the density profile and fraction of baryons retained. The relatively high temporal broadening ($\tau = 3.3 \pm 0.2 \; ms$ at 1.28 GHz) and rotation measure ($ RM = 353 \pm 2\; rad \; m^{-2}$) (Day et al. 2020) of FRB 190608 may be attributable to both turbulent gas within the spiral arm and gas local to the FRB progenitor. In contrast to previous high-resolution studies of FRB progenitor environments, we find no evidence for disturbed morphology, emission, nor kinematics for FRB 190608.

We present light curves and spectra of the tidal disruption event (TDE) ASASSN-18pg / AT 2018dyb spanning a period of one year. The event shows a plethora of strong emission lines, including the Balmer series, He II, He I and metal lines of O III $\lambda$3760 and N III $\lambda\lambda$ 4100, 4640 (blended with He II). The latter lines are consistent with originating from the Bowen fluorescence mechanism. By analyzing literature spectra of past events, we conclude that these lines are common in TDEs. The spectral diversity of optical TDEs is thus larger than previously thought and includes N-rich events besides H- and He-rich events. We study how the spectral lines evolve with time, by means of their width, relative strength, and velocity offsets. The velocity width of the lines starts at $\sim$ 13000 km s$^{-1}$ and decreases with time. The ratio of He II to N III increases with time. The same is true for ASASSN-14li, which has a very similar spectrum to AT 2018dyb but its lines are narrower by a factor of $>$2. We estimate a black hole mass of $M_{\rm BH}$ = $3.3^{+5.0}_{-2.0}\times 10^6$ $M_{\odot}$ by using the $M$-$\sigma$ relation. This is consistent with the black hole mass derived using the MOSFiT transient fitting code. The detection of strong Bowen lines in the optical spectrum is an indirect proof for extreme ultraviolet and (reprocessed) X-ray radiation and favors an accretion origin for the TDE optical luminosity. A model where photons escape after multiple scatterings through a super-Eddington thick disk and its optically thick wind, viewed at an angle close to the disk plane, is consistent with the observations.

The amount of cold, molecular gas in high-redshift galaxies is typically inferred from proxies of molecular hydrogen (H$_2$), such as carbon monoxide (CO) or neutral atomic carbon ([CI]) and molecular gas mass conversion factors. The use of these proxies, however, relies on modeling and observations that have not been directly measured outside the local universe. Here, we use recent samples of high-redshift gamma-ray burst (GRB) and quasar molecular gas absorbers to determine this conversion factor $\alpha_{\rm [CI]}$ from the column density of H$_2$, which gives us the mass per unit column, and the [CI]($J=1$) column density, which provides the luminosity per unit column. This technique allows us to make direct measurements of the relative abundances in high-redshift absorption-selected galaxies. Our sample spans redshifts of z=1.9-3.4 and covers two orders of magnitude in gas-phase metallicity. We find that $\alpha_{\rm [CI]}$ scales linearly with the metallicity: $\log \alpha_{\rm [CI]} = -1.13\times \log(Z/Z_{\odot}) + 1.33$, with a scatter of $\sigma_{\alpha_{\rm [CI]}} = 0.2$ dex. Using a sample of emission-selected galaxies at z~0-5, with both [CI] and CO line detections, we apply the $\alpha_{\rm [CI]}$ conversion to derive independent estimates of the molecular gas mass and $\alpha_{\rm CO}$. We find a remarkable agreement between the molecular gas masses inferred from the absorption-derived $\alpha_{\rm [CI]}$ compared to typical $\alpha_{\rm CO}$-based estimates, and an inferred metallicity evolution of $\alpha_{\rm CO}$ that is consistent with $\alpha_{\rm [CI]}$ and previous estimates from the literature. These results thus support the use of the absorption-derived $\alpha_{\rm [CI]}$ conversion factor for emission-selected star-forming galaxies and demonstrate that both methods probe the same universal properties of molecular gas in the local and high-redshift universe.

We report on the discovery of a peculiar Broad Absorption Line (BAL) quasar identified in our Gaia-assisted survey of red quasars. The systemic redshift of this quasar was difficult to establish due to the absence of conspicuous emission lines. Based on deep and broad BAL troughs (at least SiIV, CIV, and AlIII), a redshift of z=2.41 was established under the assumption that the systemic redshift can be inferred from the red edge of the BAL troughs. However, we observe a weak and spatially-extended emission line at 4450 AA most likely due to Lyman-alpha emission, which implies a systemic redshift of z=2.66 if correctly identified. There is also evidence for the onset of Lyman-alpha forest absorption bluewards of 4450 AA and evidence for H-alpha emission in the K-band consistent with a systemic redshift of z=2.66. If this redshift is correct, the quasar is an extreme example of a detached low-ionisation BAL quasar. The BAL lines must originate from material moving with very large velocities ranging from 22000 to 40000 km/s. To our knowledge, this is the first case of a systemic-redshift measurement based on extended Lyman-$\alpha$ emission for a BAL quasar, a method that should also be useful in cases of sufficiently distant BL Lac quasars without systemic-redshift information.

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.

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.

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.

In the early Universe, most of the cold neutral gas that will later form into individual stars and galaxies is practically invisible to us. These neutral gas reservoirs can, however, be illuminated by bright cosmic lightsources such as gamma-ray bursts (GRBs) and quasars. The aim of this thesis is to use these luminous objects as tools to study the environments of intervening or host galaxy absorption systems through cosmic time. Part I is dedicated to examining the gas, dust and metals in the immediate region surrounding GRBs. Part II presents a search for and the study of cold and molecular gas in high-z GRB host galaxy absorption systems. Part III focuses on using quasars to examine gas-rich intervening galaxies in the line of sight, with specific focus on absorption systems rich in dust and metals. This thesis demonstrates the importance of observing large samples of GRB afterglows to 1) allow for statistical studies of the GRB phenomena itself and the associated host galaxy environments and 2) to obtain spectra of peculiar or unusual GRB afterglows, that is only observed rarely. In addition, it highlights that defining a complete and unbiased sample of quasars is vital to fully exploit the potential of quasars as probes of cosmic chemical evolution. In this version of the thesis only the title, author list and abstract for each published paper, presented as individual chapters, are provided. Principal supervisor: Prof. P\'all Jakobsson.

Binary quasars are extremely rare objects, used to investigate clustering on very small scales at different redshifts. The cases where the two quasar components are gravitationally bound, known as physical binary quasars, can also exhibit enhanced astrophysical activity and therefore are of particular scientific interest. Here we present the serendipitous discovery of a physical pair of quasars with an angular separation of $\Delta\theta = (8.76 \pm 0.11)$ arcsec. The redshifts of the two quasars are consistent within the errors and measured as $z = (1.76\pm 0.01)$. Under the motivated assumption that the pair does not arise from a single gravitationally lensed quasar, the resulting projected physical separation was estimated as $(76 \pm 1)$ kpc. For both targets we detected Si VI, C VI, C III], and Mg II emission lines. However, the two quasars show significantly different optical colours, one being among the most reddened quasars at $z > 1.5$ and the other with colours consistent with typical quasar colours at the same redshift. Therefore it is ruled out that the sources are a lensed system. This is our second serendipitous discovery of a pair of two quasars with different colours, having a separation $\lesssim 10$ arcsec, which extends the very limited catalogue of known quasar pairs. We ultimately argue that the number of binary quasars may have been significantly underestimated in previous photometric surveys, due to the bias arising from paired quasars with very different colours.

Depletion studies provide a way to understand the chemical composition of interstellar dust grains. We here examine 23 gamma-ray bursts (GRB) optical afterglow spectra (spanning 0.6<z<5.0) and compare their silicon and iron dust-phase column densities with different extinction curve parameters to study the composition of the interstellar dust grains in these high-redshift GRB host galaxies. The majority of our sample (87%) show featureless extinction curves and only vary in shape. We observe strong correlations (with >96% significance) between the total-to-selective extinction, R_V, and the dust-phase column densities of Si and Fe. Since a large fraction of interstellar iron is locked in silicate grains, this indicates that high Si and Fe depletion leads to an increase in the fraction of large silicate grains and vice versa. This suggests that silicates play a vital role to induce the entire extinction at any wavelength. On the other hand, the far-UV extinction is usually attributed to the presence of small silicates. However, we find no trend between the far-UV parameter of the extinction curve, c_4, and the abundance of Si and Fe in the dust phase. We, therefore, propose that the far-UV extinction could be a combined effect of small (probably nanoparticles) dust grains from various species.

Quasar lines of sight intersect intervening galaxy discs or circum-galactic environments at random impact parameters and potential well depths. Absorption line velocity widths ($\Delta v_{90}$) are known to scale with host galaxy stellar masses, and inversely with the projected separation from the quasar line of sight. Its dependence on stellar mass can be eliminated by normalising with the emission-line widths of the host galaxies, $\sigma_{em}$, so that absorbers with a range of $\Delta v_{90}$ values can be compared directly. Using a sample of DLA systems at 0.2 < z < 3.2 with spectroscopically confirmed host galaxies, we find that the velocity ratio $\Delta v_{90}/\sigma_{em}$ decreases with projected distances from the hosts. We compare the data with expectations of line-of-sight velocity dispersions derived for different dark matter halo mass distributions, and find that models with steeper radial dark matter profiles provide a better fit to the observations, although the scatter remains large. Gas outflows from the galaxies may cause an increased scatter, or scale radii of dark matter halo models may not be representative for the galaxies. We demonstrate by computing virial velocities, that metal-rich DLAs that belong to massive galaxy halos (M$_{halo} \approx 10^{12}$ M$_{\odot}$) mostly remain gravitationally bound to the halos.

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 report our identification of the optical afterglow and host galaxy of the short-duration gamma-ray burst GRB 160821B. The spectroscopic redshift of the host is $z=0.162$, making it one of the lowest redshift sGRBs identified by Swift. Our intensive follow-up campaign using a range of ground-based facilities as well as HST, XMM and Swift, shows evidence for a late-time excess of optical and near-infrared emission in addition to a complex afterglow. The afterglow light-curve at X-ray frequencies reveals a narrow jet, $\theta_j\sim1.9^{+0.10}_{-0.03}$ deg, that is refreshed at $>1$ day post-burst by a slower outflow with significantly more energy than the initial outflow that produced the main GRB. Observations of the 5 GHz radio afterglow shows a reverse shock into a mildly magnetised shell. The optical and near-infrared excess is fainter than AT2017gfo associated with GW170817, and is well explained by a kilonova with dynamic ejecta mass $M_{\rm dyn}=(1.0\pm0.6)\times10^{-3}$ M$_{\odot}$ and a secular (postmerger) ejecta mass with $M_{\rm pm}=(1.0\pm0.6)\times10^{-2}$ M$_\odot$, consistent with a binary neutron star merger resulting in a short-lived massive neutron star. This optical and near-infrared dataset provides the best-sampled kilonova light-curve without a gravitational wave trigger to date.

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.

We present a systematic study of the impact of a dust bias on samples of damped Lyman-$\alpha$ absorbers (DLAs). This bias arises as an effect of the magnitude and colour criteria utilized in the Sloan Digital Sky Survey (SDSS) quasar target selection up until data release 7 (DR7). The bias has previously been quantified assuming only a contribution from the dust obscuration. In this work, we apply the full set of magnitude and colour criteria used up until SDSS DR7 in order to quantify the full impact of dust biasing against dusty and metal-rich DLAs. We apply the quasar target selection algorithm on a modelled population of intrinsic colours, and by exploring the parameter space consisting of redshift, ($z_{\rm QSO}$ and $z_{\rm abs}$), optical extinction, and HI column density, we demonstrate how the selection probability depends on these variables. We quantify the dust bias on the following properties derived for DLAs at z=3: the incidence rate, the mass density of neutral hydrogen and metals, the average metallicity. We find that all quantities are significantly affected. When considering all uncertainties, the mass density of neutral hydrogen is underestimated by 10 to 50%, and the mass density in metals is underestimated by 30 to 200%. Lastly, we find that the bias depends on redshift. At redshift z=2.2, the mass density of neutral hydrogen and metals might be underestimated by up to a factor of 2 and 5, respectively. Characterizing such a bias is crucial in order to accurately interpret and model the properties and metallicity evolution of absorption-selected galaxies.

Damped Lyman-{\alpha} Absorbers (DLAs) as a class of QSO absorption-line systems are currently our most important source of detailed information on the cosmic chemical evolution of galaxies. However, the degree to which this information is biased by dust remains to be understood. One strategy is to specifically search for QSOs reddened by metal-rich and dusty foreground absorbers. In this Letter we present the discovery of a z=2.60 QSO strongly reddened by dust in an intervening extremely-strong DLA at z=2.226. This QSO was identified through a novel selection combining the astrometric measurements from ESA's Gaia satellite with extent optical and near/mid-infrared photometry. We infer a total neutral atomic-hydrogen column density of log N(HI)=21.95{\pm}0.15 and a lower limit on the gas-phase metallicity of [Zn/H]>-0.96. This DLA is also remarkable in exhibiting shielded neutral gas witnessed in CI and tentative detections of CO molecular bands. The Spectral Energy Distribution (SED) of the QSO is well-accounted for by a normal QSO-SED reddened by dust from a DLA with a 10%-of-Solar metallicity, dust extinction of A_V=0.82{\pm}0.02mag, and LMC-like extinction curve including the characteristic 2175{\AA} extinction feature. Such QSO absorption-line systems have shown to be very rare in previous surveys, which have mostly revealed sight-lines with low extinction. The present case therefore suggests that previous samples have under-represented the fraction of dusty absorbers. Building a complete sample of such systems is needed to assess the significance of this effect.

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.

Cosmological simulations predict that an intergalactic magnetic field (IGMF) pervades the large scale structure (LSS) of the Universe. Measuring the IGMF is important to determine its origin (i.e. primordial or otherwise). Using data from the LOFAR Two Metre Sky Survey (LoTSS), we present the Faraday rotation measure (RM) and depolarisation properties of the giant radio galaxy J1235+5317, at a redshift of $z = 0.34$ and 3.38 Mpc in size. We find a mean RM difference between the lobes of $2.5\pm0.1$ rad/m$^2$ , in addition to small scale RM variations of ~0.1 rad/m$^2$ . From a catalogue of LSS filaments based on optical spectroscopic observations in the local universe, we find an excess of filaments intersecting the line of sight to only one of the lobes. Associating the entire RM difference to these LSS filaments leads to a gas density-weighted IGMF strength of ~0.3 {\mu}G. However, direct comparison with cosmological simulations of the RM contribution from LSS filaments gives a low probability (~5%) for an RM contribution as large as 2.5 rad/m$^2$ , for the case of IGMF strengths of 10 to 50 nG. It is likely that variations in the RM from the Milky Way (on 11' scales) contribute significantly to the mean RM difference, and a denser RM grid is required to better constrain this contribution. In general, this work demonstrates the potential of the LOFAR telescope to probe the weak signature of the IGMF. Future studies, with thousands of sources with high accuracy RMs from LoTSS, will enable more stringent constraints on the nature of the IGMF.

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.

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 report on the results of a search for CO(3-2) emission from the galaxy counterpart of a high-metallicity Damped Ly-alpha Absorber (DLA) at z=2.5832 towards the quasar Q0918+1636. We do not detect CO emission from the previously identified DLA galaxy counterpart. The limit we infer on M_gas / M_star is in the low end of the range found for DLA galaxies, but is still consistent with what is found for other star-forming galaxies at similar redshifts. Instead we detect CO(3-2) emission from another intensely star-forming galaxy at an impact parameter of 117 kpc from the line-of-sight to the quasar and 131 km s^-1 redshifted relative to the velocity centroid of the DLA in the quasar spectrum. In the velocity profile of the low- and high-ionisation absorption lines of the DLA there is an absorption component consistent with the redshift of this CO-emitting galaxy. It is plausible that this component is physically associated with a strong outflow in the plane of the sky from the CO-emitting galaxy. If true, this would be further evidence, in addition to what is already known from studies of Lyman-break galaxies, that galactic outflows can be traced beyond 100 kpc from star-forming galaxies. The case of this z=2.583 structure is an illustration of this in a group environment.

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.

We here apply a novel technique selecting quasar candidates purely as sources with zero proper motions in the Gaia data release 2 (DR2). We demonstrate that this approach is highly efficient toward high Galactic latitudes with < 25% contamination from stellar sources. Such a selection technique offers a very pure sample completeness, since all cosmological point sources are selected regardless of their intrinsic spectral properties within the limiting magnitude of Gaia. We carry out a pilot-study by defining a sample compiled by including all Gaia-DR2 sources within one degree of the North Galactic Pole (NGP) selected to have proper motions consistent with zero within 2-sigma uncertainty. By cross-matching the sample to the optical Sloan Digital Sky Survey (SDSS) and the mid-infrared AllWISE photometric catalogues we investigate the colours of each of our sources. Together with already spectroscopically confirmed quasars we are therefore able to determine the efficiency of our selection. The majority of the zero proper motion sources have optical to mid-infrared colours consistent with known quasars. The remaining population may be contaminating stellar sources, but some may also be quasars with colours similar to stars. Spectroscopic follow-up of the zero proper motion sources is needed to unveil such a hitherto hidden quasar population. This approach has the potential to allow substantial progress on many important questions concerning quasars such as determining the fraction of dust-obscured quasars, the fraction of broad absorption line (BAL) quasars, and the metallicity distribution of damped Lyman-$\alpha$ absorbers. The technique could also potentially reveal new types of quasars or even new classes of cosmological point sources.

The cosmic chemical enrichment as measured from damped Ly$\alpha$ absorbers (DLAs) will be underestimated if dusty and metal-rich absorbers have evaded identification. Here we report the discovery and present the spectroscopic observations of a quasar, KV-RQ\,1500-0031, at $z=2.520$ reddened by a likely dusty DLA at $z=2.428$ and a strong MgII absorber at $z=1.603$. This quasar was identified as part of the KiDS-VIKING Red Quasar (KV-RQ) survey, specifically aimed at targeting dusty absorbers which may cause the background quasars to escape the optical selection of e.g. the SDSS quasar survey. For the DLA we find an HI column density of $\log N$(HI) = $21.2\pm 0.1$ and a metallicity of [X/H] = $-0.90\pm 0.20$ derived from an empirical relation based on the equivalent width of SiII$\lambda$1526. We observe a total visual extinction of $A_V=0.16$ mag induced by both absorbers. We compile a sample of 17 additional dusty ($A_V > 0.1$ mag) DLAs toward quasars (QSO-DLAs) from the literature for which we characterize the properties of HI column density, metallicity and dust. From this sample we also estimate a correction factor to the overall DLA metallicity budget. We demonstrate that the dusty QSO-DLAs have high metal column densities ($\log N$(HI) + [X/H]) and are more similar to gamma-ray burst (GRB)-selected DLAs (GRB-DLAs) than regular QSO-DLAs. We evaluate the effect of dust reddening in DLAs as well as illustrate how the induced color excess of the underlying quasars can be significant (up to $\sim 1$ mag in various optical bands), even for low to moderate extinction values ($A_V \lesssim 0.6$ mag). Finally we discuss the direct and indirect implications of a significant dust bias in both QSO- and GRB-DLA samples. [Abridged]

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.

Gravitational waves were discovered with the detection of binary black hole mergers and they should also be detectable from lower mass neutron star mergers. These are predicted to eject material rich in heavy radioactive isotopes that can power an electromagnetic signal called a kilonova. The gravitational wave source GW170817 arose from a binary neutron star merger in the nearby Universe with a relatively well confined sky position and distance estimate. Here we report observations and physical modelling of a rapidly fading electromagnetic transient in the galaxy NGC4993, which is spatially coincident with GW170817 and a weak short gamma-ray burst. The transient has physical parameters broadly matching the theoretical predictions of blue kilonovae from neutron star mergers. The emitted electromagnetic radiation can be explained with an ejected mass of 0.04 +/- 0.01 Msol, with an opacity of kappa <= 0.5 cm2/gm at a velocity of 0.2 +/- 0.1c. The power source is constrained to have a power law slope of beta = -1.2 +/- 0.3, consistent with radioactive powering from r-process nuclides. We identify line features in the spectra that are consistent with light r-process elements (90 < A < 140). As it fades, the transient rapidly becomes red, and emission may have contribution by a higher opacity, lanthanide-rich ejecta component. This indicates that neutron star mergers produce gravitational waves, radioactively powered kilonovae, and are a nucleosynthetic source of the r-process elements.

We present a comprehensive study of the relations between gas kinematics, metallicity, and stellar mass in a sample of 82 GRB-selected galaxies using absorption and emission methods. We find the velocity widths of both emission and absorption profiles to be a proxy of stellar mass. We also investigate the velocity-metallicity correlation and its evolution with redshift and find the correlation derived from emission lines to have a significantly smaller scatter compared to that found using absorption lines. Using 33 GRB hosts with measured stellar mass and metallicitiy, we study the mass-metallicity relation for GRB host galaxies in a stellar mass range of $10^{8.2} M_{\odot}$ to $10^{11.1} M_{\odot}$ and a redshift range of $ z\sim 0.3-3.4$. The GRB-selected galaxies appear to track the mass-metallicity relation of star forming galaxies but with an offset of 0.15 towards lower metallicities. This offset is comparable with the average error-bar on the metallicity measurements of the GRB sample and also the scatter on the MZ relation of the general population. It is hard to decide whether this relatively small offset is due to systematic effects or the intrinsic nature of GRB hosts. We also investigate the possibility of using absorption-line metallicity measurements of GRB hosts to study the mass-metallicity relation at high redshifts. Our analysis shows that the metallicity measurements from absorption methods can significantly differ from emission metallicities and assuming identical measurements from the two methods may result in erroneous conclusions.

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.

To fully exploit the potential of quasars as probes of cosmic chemical evolution and the internal gas dynamics of galaxies it is important to understand the selection effects behind the quasar samples and in particular if the selection criteria exclude foreground galaxies with certain properties (most importantly a high dust content). Here we present spectroscopic follow-up from the 10.4-m GTC telescope of a dust-reddened quasar, eHAQ0111+0641, from the extended High A_V Quasar (HAQ) survey. We find that the z=3.21 quasar has a foreground Damped Lyman-alpha Absorber (DLA) at z=2.027 along the line of sight. The DLA has very strong metal lines due to a moderately high metallicity (with an inferred lower limit of 25% of the solar metallicity), but a very large gas column density along the line-of-sight in its host galaxy. This discovery is further evidence that there is a dust bias affecting the census of metals, caused by the combined effect of dust obscuration and reddening, in existing samples of z>2 DLAs. The case of eHAQ0111+0641 illustrates that dust bias is not only caused by dust obscuration, but also dust reddening.

Since the first discovery of a broad-lined type Ic supernova (SN) with a long-duration gamma-ray burst (GRB) in 1998, fewer than fifty gamma-ray burst supernovae (GRB-SNe) have been discovered. The intermediate-luminosity Swift GRB 161219B and its associated supernova SN 2016jca, which occurred at a redshift of z=0.1475, represents only the seventh GRB-SN to have been discovered within 1 Gpc, and hence provides an excellent opportunity to investigate the observational and physical properties of these very elusive and rare type of SN. As such, we present optical to near-infrared photometry and optical spectroscopy of GRB 161219B and SN 2016jca, spanning the first three months since its discovery. GRB 161219B exploded in the disk of an edge-on spiral galaxy at a projected distance of 3.4 kpc from the galactic centre. GRB 161219B itself is an outlier of the Amati relation, while SN 2016jca had a rest-frame, peak absolute V-band magnitude of M_V = -19.0, which it reached after 12.5 rest-frame days. We find that the bolometric properties of SN 2016jca are inconsistent with being powered solely by a magnetar central engine, as proposed by other authors, and demonstrate that it was likely powered exclusively by energy deposited by the radioactive decay of nickel and cobalt into their daughter products, which were nucleosynthesized when its progenitor underwent core collapse. We find that 0.22 solar masses of nickel is required to reproduce the peak luminosity of SN 2016jca, and we constrain an ejecta mass of 5.8 solar masses and a kinetic energy of ~5 x 10^52 erg. Finally, we report on a chromatic, pre-maximum bump in the g-band light curve, and discuss its possible origin. [Abridged]

Gamma-ray bursts (GRBs) are powerful probes of early stars and galaxies, during and potentially even before the era of reionization. Although the number of GRBs identified at z>6 remains small, they provide a unique window on typical star-forming galaxies at that time, and thus are complementary to deep field observations. We report the identification of the optical drop-out afterglow of Swift GRB 120923A in near-infrared Gemini-North imaging, and derive a redshift of z=7.84_{-0.12}^{+0.06} from VLT/X-shooter spectroscopy. At this redshift the peak 15-150 keV luminosity of the burst was 3.2x10^52 erg/s, and in fact the burst was close to the Swift/BAT detection threshold. The X-ray and near-infrared afterglow were also faint, and in this sense it was a rather typical long-duration GRB in terms of rest-frame luminosity. We present ground- and space-based follow-up observations spanning from X-ray to radio, and find that a standard external shock model with a constant-density circumburst environment with density, n~4x10^-2 cm^-3 gives a good fit to the data. The near-infrared light curve exhibits a sharp break at t~3.4 days in the observer frame, which if interpreted as being due to a jet corresponds to an opening angle of ~5 degrees. The beaming corrected gamma-ray energy is then E_gamma~2x10^50 erg, while the beaming-corrected kinetic energy is lower, E_K~10^49 erg, suggesting that GRB 120923A was a comparatively low kinetic energy event. We discuss the implications of this event for our understanding of the high-redshift population of GRBs and their identification.

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.

We present the results of a new spectroscopic survey for dusty intervening absorption systems, particularly damped Ly$\alpha$ absorbers (DLAs), towards reddened quasars. The candidate quasars are selected from mid-infrared photometry from the Wide-field Infrared Survey Explorer combined with optical and near-infrared photometry. Out of 1073 candidates, we secure low-resolution spectra for 108 using the Nordic Optical Telescope on La Palma, Spain. Based on the spectra, we are able to classify 100 of the 108 targets as quasars. A large fraction (50 %) is observed to have broad absorption lines (BALs). Moreover, we find 6 quasars with strange breaks in their spectra, which are not consistent with regular dust reddening. Using template fitting we infer the amount of reddening along each line of sight ranging from A(V)$\approx$0.1 mag to 1.2 mag (assuming an SMC extinction curve). In four cases, the reddening is consistent with dust exhibiting the 2175{\AA} feature caused by an intervening absorber, and for two of these, a MgII absorption system is observed at the best-fit absorption redshift. In the rest of the cases, the reddening is most likely intrinsic to the quasar. We observe no evidence for dusty DLAs in this survey. However, the large fraction of BAL quasars hampers the detection of absorption systems. Out of the 50 non-BAL quasars only 28 have sufficiently high redshift to detect Ly$\alpha$ in absorption.

The sub-population of quasars reddened by intrinsic or intervening clouds of dust are known to be underrepresented in optical quasar surveys. By defining a complete parent sample of the brightest and spatially unresolved quasars in the COSMOS field, we quantify to which extent this sub-population is fundamental to our understanding of the true population of quasars. By using the available multiwavelength data of various surveys in the COSMOS field, we built a parent sample of 33 quasars brighter than $J=20$ mag, identified by reliable X-ray to radio wavelength selection techniques. Spectroscopic follow-up with the NOT/ALFOSC was carried out for four candidate quasars that had not been targeted previously to obtain a 100\% redshift completeness of the sample. The population of high $A_V$ quasars (HAQs), a specific sub-population of quasars selected from optical/near-infrared photometry, is found to contribute $21\%^{+9}_{-5}$ of the parent sample. The full population of bright spatially unresolved quasars represented by our parent sample consists of $39\%^{+9}_{-8}$ reddened quasars defined by having $A_V>0.1$, and $21\%^{+9}_{-5}$ of the sample having $E(B-V)>0.1$ assuming the extinction curve of the Small Magellanic Cloud. We show that the HAQ selection works well for selecting reddened quasars, but some are missed because their optical spectra are too blue to pass the $g-r$ color cut in the HAQ selection. This is either due to a low degree of dust reddening or anomalous spectra. We find that the fraction of quasars with contributing light from the host galaxy is most dominant at $z \lesssim 1$. At higher redshifts the population of spatially unresolved quasars selected by our parent sample is found to be representative of the full population at $J<20$ mag. This work quantifies the bias against reddened quasars in studies that are based solely on optical surveys.

We present the serendipitous discovery of a projected pair of quasi-stellar objects (QSOs) with an angular separation of $\Delta\theta =4.50$ arcsec. The redshifts of the two QSOs are widely different: one, our programme target, is a QSO with a spectrum consistent with being a narrow line Seyfert 1 AGN at $z=2.05$. For this target we detect Lyman-$\alpha$, \ion{C}{4}, and \ion{C}{3]}. The other QSO, which by chance was included on the spectroscopic slit, is a Type 1 QSO at a redshift of $z=1.68$, for which we detect \ion{C}{4}, \ion{C}{3]} and \ion{Mg}{2}. We compare this system to previously detected projected QSO pairs and find that only about a dozen previously known pairs have smaller angular separation.

Selection of extragalactic point sources, e.g. QSOs, is often hampered by significant selection effects causing existing samples to have rather complex selection functions. We explore whether a purely astrometric selection of extragalactic point sources, e.g. QSOs, is feasible with the ongoing Gaia mission. Such a selection would be interesting as it would be unbiased in terms of colours of the targets and hence would allow selection also with colours in the stellar sequence. We have analyzed a total of 18 representative regions of the sky by using GUMS, the simulator prepared for ESAs Gaia mission, both in the range of $12\le G \le 20$ mag and $12\le G \le 18$ mag. For each region we determine the density of apparently stationary stellar sources, i.e. sources for which Gaia cannot measure a significant proper motion. The density is contrasted with the density of extragalactic point sources, e.g. QSOs, in order to establish in which celestial directions a pure astrometric selection is feasible. When targeting regions at galactic latitude $|b| \ge 30^\mathrm{o}$ the ratio of QSOs to apparently stationary stars is above 50\% and when observing towards the poles the fraction of QSOs goes up to about $\sim80$\%. We show that the proper motions from the proposed Gaia successor mission in about 20 years would dramatically improve these results at all latitudes. Detection of QSOs solely from zero proper motion, unbiased by any assumptions on spectra, might lead to the discovery of new types of QSOs or new classes of extragalactic point sources.