Arxiv:1802.07727V1 [Astro-Ph.HE] 21 Feb 2018 Tion Systems to Standard Candles in Cosmology (E.G., Wijers Et Al

Home , Galactic plane, Galactic quadrant, GRB 100621A, GRB 160625B, Interacting galaxy

Astronomy & Astrophysics manuscript no. XSGRB_sample_arxiv c ESO 2018 2018-02-23

The X-shooter GRB afterglow legacy sample (XS-GRB)?

J. Selsing1,??, D. Malesani1, 2, 3,†, P. Goldoni4,†, J. P. U. Fynbo1, 2,†, T. Krühler5,†, L. A. Antonelli6,†, M. Arabsalmani7, 8, J. Bolmer5, 9,†, Z. Cano10,†, L. Christensen1, S. Covino11,†, P. D’Avanzo11,†, V. D’Elia12,†, A. De Cia13, A. de Ugarte Postigo1, 10,†, H. Flores14,†, M. Friis15, 16, A. Gomboc17, J. Greiner5, P. Groot18, F. Hammer14, O.E. Hartoog19,†, K. E. Heintz1, 2, 20,†, J. Hjorth1,†, P. Jakobsson20,†, J. Japelj19,†, D. A. Kann10,†, L. Kaper19, C. Ledoux9, G. Leloudas1, A.J. Levan21,†, E. Maiorano22, A. Melandri11,†, B. Milvang-Jensen1, 2, E. Palazzi22, J. T. Palmerio23,†, D. A. Perley24,†, E. Pian22, S. Piranomonte6,†, G. Pugliese19,†, R. Sánchez-Ramírez25,†, S. Savaglio26, P. Schady5, S. Schulze27,†, J. Sollerman28, M. Sparre29,†, G. Tagliaferri11, N. R. Tanvir30,†, C. C. Thöne10, S.D. Vergani14,†, P. Vreeswijk18, 26,†, D. Watson1, 2,†, K. Wiersema21, 30,†, R. Wijers19, D. Xu31,†, and T. Zafar32

(Aﬃliations can be found after the references)

Received/ accepted

ABSTRACT

In this work we present spectra of all γ-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 ∼ 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. Key words. Gamma-ray burst: general — Galaxies: high-redshift — ISM: general — Techniques: spectroscopic — Galaxies: star formation — catalogs

1. Introduction ati et al. 2008; Vergani et al. 2009; Prochaska et al. 2009; Hjorth & Bloom 2011; Rowlinson et al. 2017; Christensen et al. 2017). Gamma-ray bursts (GRBs) are bright glimpses of electromagnetic radiation that pierce through the universe, all the way from The Neil Gehrels Swift Observatory (Swift) satellite (Gehrels the edges of the observable universe. They provide constraints et al. 2004, 2009), which was launched in 2004, has made it pos- on a very wide range of topics in astrophysics. Examples range sible to harvest much of the great potential in using GRBs as from small-scale phenomena relating to magnetars, properties probes of the intergalactic medium, which was already hinted of highly relativistic jets, hyper/supernova explosions, the inter- at by results from earlier missions (e.g., van Paradijs et al. stellar medium, dust extinction curves, starburst galaxies, chem- 2000; Ricker & Team 2004). With the three on-board instru- ical and molecular abundances, escape of ionizing radiation, the ments, the Burst Alert Telescope (BAT; Barthelmy et al. 2005), ionization state of the intergalactic medium, intervening absorp- the X-Ray Telescope (XRT; Burrows et al. 2005), and the Ul- arXiv:1802.07727v1 [astro-ph.HE] 21 Feb 2018 tion systems to standard candles in cosmology (e.g., Wijers et al. traViolet and Optical Telescope (UVOT; Roming et al. 2005), 1998; Savaglio 2006; Ghirlanda 2007; Molinari et al. 2007; Am- Swift is an ideal observatory for GRB hunting. A crucial as- pect of the success of the Swift mission has been the extensive

? ground-based follow-up observations of the afterglows and of Based on observations collected at the European Southern Obser- the host galaxies of the GRBs, involving a large community vatory, Paranal, Chile, Program IDs: 084.A-0260, 084.D-0265, 085.A- of researchers. This fruitful collaboration has been facilitated 0009, 086.A-0073, 087.A-0055, 088.A-0051, 089.A-0067, 090.A- 0088, 091.A-0877, 0091.C-0934, 092.A-0124, 092.D-0056, 092.D- by the open data access policy of the Swift mission. The close 0633, 093.A-0069, 094.A-0134, 095.A-0045, 095.B-0811, 096.A-0079, collaboration between detection facilities and electromagnetic 097.A-0036, 098.A-0136, and 098.A-0055. follow-up campaigns continue to be immensely rewarding, as ?? e-mail: [email protected] recently highlighted by the simultaneous detection of gravita- † On-call observer tional waves and light from the neutron star merger in the shape

Article number, page 1 of 43 A&A proofs: manuscript no. XSGRB_sample_arxiv of GW170817/GRB170817A/AT2017gfo (Abbott et al. 2017; tribution. When investigating a specific distribution function, a LIGO Scientific Collaboration et al. 2017). high degree of completeness is desired (e.g., Perley et al. 2016b). In the beginning of the Swift era most of the follow-up af- In defining the selection criteria, we simultaneously aim to terglow spectroscopy was secured using low-resolution spectro- minimize any biases against intrinsic astrophysical conditions graphs (typically R = λ/∆λ<1000, e.g., Fynbo et al. 2009). while at the same time maximizing the likelihood of success- Spectroscopy is powerful as it allows us to secure information ful observations, hence allowing us to obtain a higher degree of even for very faint targets (Krühler et al. 2012). This allows completeness. By restricting the selection criteria to conditions the measurement of a number of important parameters such as local to the Milky Way and therefore independent of intrinsic redshifts, spectral slopes, and extinction. For a handful of very GRB properties, the aim is that the collected sample represents bright afterglows high-resolution (typically R>20000) spectra the underlying distribution of GRBs in a fair way. The selection have been secured, and for these events much more informa- criteria used here are based on previous, similar studies (Jakobs- tion about conditions inside the host galaxies were extracted son et al. 2006; Fynbo et al. 2009; Hjorth et al. 2012). We char- (e.g., Fiore et al. 2005; Thöne et al. 2007; Prochaska et al. 2007; acterize the sample completeness in Sect. 4.2. Vreeswijk et al. 2007; D’Elia et al. 2009; Castro-Tirado et al. The selection criteria that define a GRB as part of our initial 2010). statistical sample are: The X-shooter spectrograph (Vernet et al. 2011) was the first of the second generation instruments at the ESO Very Large 1. GRB trigger by BAT onboard the Swift satellite Telescope (VLT). It was designed very much with transient 2. XRT started observing within 10 minutes after the GRB; an follow-up in mind as the fading luminosities of such sources XRT position must be distributed within 12 hr. makes it urgent to secure as extensive wavelength coverage as 3. The target must be visible from Cerro Paranal for at least 60 minutes, 30 degrees above the horizon, with the Sun below possible in the shortest possible time. At the same time, the res- 1 olution was designed to be in the range 4000–9000 in order to −12 degrees . be able to get a large useful spectral range between the many 4. Galactic AV ≤ 0.5 mag according to the maps of Schlegel sky-background emission lines in the red and near-IR spectral et al.(1998). ranges. The near-IR spectral coverage allows for spectroscopic 5. No bright, nearby stars within 1.8 + 0.4 × exp[(20 − R)/2.05] observations of the highest redshift GRBs. arcsec, where R is the USNO magnitude of the star). In this paper, we present the results of a dedicated effort over Our ability to observe GRB afterglows is strongly depen- the years 2009 – 2017 to use the X-shooter spectrograph to se- dent on the timing and the precision of the target positions de- cure spectroscopic observations of afterglows and host galaxies livered by the triggering facilities. By selecting only bursts that of GRBs detected by Swift. We here make all the data result- have been triggered on board the Swift space telescope (Gehrels ing from the survey publicly available in reduced form (see Sect et al. 2004), based on the BAT, we start out with a sample where 3.7). burst characteristics are delivered immediately, allowing for an The paper is organized in the following way: In Sect.2 we informed follow-up strategy. However, the BAT sensitivity varies describe the sample including the sample selection criteria and across its field of view, so selection is not entirely homogeneous. the observational strategy. In Sect. 2.4, we describe the obser- Despite the complexity of the triggering mechanism on board vations and the instrumental setups, and in Sect.3 we detail the Swift (Band 2006; Coward et al. 2013), attempts at inferring methodological strategies adopted in the data reduction process properties of the underlying GRB population based on the detec- and auxiliary material. In Sect.4 we describe the results of the tion thresholds and triggering algorithms have been made (Lien survey, i.e. the efficiency of the follow-up effort and the char- et al. 2014; Graff et al. 2016). Restricting the follow-up effort to acteristics of the observed bursts. We also assess the complete- bursts detected by Swift, we therefore ensure that the limitations ness of the realized sample. Finally, we offer our conclusions in of the parent sample are well studied. Sect.5. We use the ΛCDM cosmology parameters provided by Because the localization accuracy of BAT is 1 – 4 arcmin- the Planck Collaboration et al.(2016) in which the universe is −1 −1 utes (Barthelmy et al. 2005), an afterglow identification based flat with H0 = 67.7 km s Mpc and Ωm = 0.307. on BAT alone would be harder, and host association impossible. We therefore additionally require an X-ray position from XRT to be distributed to the GCN network (Barthelmy 2000) within 2. Sample selection criteria and observations 12 hours and to account for observing constraints on Swift that XRT began observations within 10 minutes. The additional tim- 2.1. Sample selection criteria ing requirement of the XRT follow-up means that all bursts in Being of transient nature, it is difficult to impose strong sam- our sample have detected X-ray afterglows. Because the XRT ple selection criteria on GRBs without hampering the follow-up completeness is very high for promptly-repointed GRBs (Bur- effort. Many natural follow-up restrictions exist already, being rows et al. 2007), this cut should not alter the parent sample sig- it weather conditions, pointing restrictions of the telescope or nificantly. unconstrained burst localizations as reported by the alerting fa- To ensure a minimum of observability, we require that the cility. To maximize the return of the follow-up campaign, we GRB is visible from the telescope site at Cerro Paranal, Chile, − have chosen a few selection criteria that facilitate an unbiased for a least 1 hour after the trigger with the sun below 12 de- selection of bursts, while at the same time allowing for a high grees. This secures time for the spectroscopic observations to be follow-up success rate. The importance of defining unbiased se- completed. Since the GRB population is isotropically distributed lection criteria has been highlighted previously (Jakobsson et al. on the sky as seen from Earth, and because the GRB properties 2006; Salvaterra et al. 2012; Hjorth et al. 2012; Vergani et al. 1 Note that in the P84 proposal the criteria have been stated differ- 2015; Perley et al. 2016a), when trying to address the intrinsic ently, the visibility constraint being replaced by a declination + Sun underlying distribution functions such as the redshift distribu- angle constraint. The criteria listed above are, however, those defining tion, host metallicity distribution, or afterglow brightness dis- the sample.

Article number, page 2 of 43 Selsing et al.: The X-shooter GRB afterglow legacy sample do not depend on position on the sky (Meegan et al. 1992; Briggs commence within minutes of the burst. One limitation is that no et al. 1996; Ukwatta & Wózniak 2016), this cut does not influ- instrument change is allowed in RRM, which lowers the number ence our ability to fairly sample the underlying GRB population. of successful triggers (X-shooter shares the telescope with two The same arguments apply to the requirement that there are no other instruments). In case of a promptly visible GRB, a robotic nearby foreground bright stars. We additionally require that the trigger is sent to the telescope if at the time of the GCN notice Galactic extinction is below AV . 0.5 mag, based on the extinc- the GRB fulfills the following criteria: tion maps by Schlegel et al.(1998) 2. Choosing low-extinction sightlines also reduces the problem of field crowding and the 1. The GRB triggered Swift onboard. contamination from Galactic outbursts posing as GRB impos- 2. The X-ray position must be available less than 1 hr after the GRB. tors. These additional cuts should not influence the optical prop- ◦ erties of the bursts themselves, only our ability to successfully 3. The elevation of the source in the sky is > 22 (both at trigger time and 15 min after). secure the observations that allow us to investigate the spectro- ◦ scopic properties of GRBs. 4. The Sun elevation from Paranl is < −12 (both at trigger time Bursts that fulfill these selection criteria are what we de- and 15 min after). fine as our initial statistical sample. To reach a higher degree 5. Galactic AV ≤ 1.0 mag according to the maps of Schlegel of sample completeness, we impose additional cuts that increase et al.(1998). our ability to infer population properties. These are presented in 6. The Swift-circulated tags: KNOWN_SOURCE, COSMIC_RAY, DEF_NOT_GRB are set to false Sect. 4.2.1. Using this sample we will be able to address pop- 3 ulation properties of Swift-detected bursts. We further discuss (PROB_NOT_GRB can be true) . the effect of these selection criteria and their implication for the The criteria to trigger RRM are looser than those that define completeness of the sample in Sect. 4.2. the statistical sample. This is both because of the expected larger brightness of GRB counterparts soon after the explosion, and 2.2. Follow-up procedure because of the rarity of RRM triggers. The use of RRM is unique as it allows to sample a long loga- Our collaboration has set up a procedure to promptly react to rithmic time span in the GRB lifetime, and it exploits the extreme alerts from Swift and other satellites (such as Fermi). We pro- brightness of early afterglows. A pivotal exampl is the study of cess automatically GCN notices in order to flag those events be- temporal variability of GRB afterglow absorption systems due longing to our statistical sample. Two people are permanently on to effect of the GRB itself on the sorrounding medium (e.g., see alert to manually supervise each event, for example to recognize Dessauges-Zavadsky et al. 2006; Vreeswijk et al. 2007; D’Elia events of special interest beyond their inclusion in the sample. et al. 2009; Vreeswijk et al. 2013). Members of our collaboration have access to numerous facilities There are nine GRBs that have been observed with X-shooter spread throughout the world, among which the William Herschel in RRM mode. One of the RRM triggers is outside the statistical Telescope (WHT), Nordic Optical Telescope (NOT), the Gran sample, and two of the RRM triggers are on short GRBs. Our Telescopio Canarias (GTC), the Telescopio Nazionale Galileo fastest response (between the Swift GRB trigger time and the (TNG), the Gamma-Ray Burst Optical Near-infrared Detector beginning of spectroscopic integration) was for GRB 160410A, (GROND Greiner et al. 2008), and the Xinlong observatory. for which the delay was only 8.4 minutes. Thanks to this network, we can often conduct searches for op- In many cases, unfortunately, we could not use RRM even tical/NIR afterglows ahead of the time of the X-shooter spec- for promptly visible events, e.g. because of the unavailability troscopy, thus helping to plan the spectroscopic observations. If of X-shooter. Potentially, the proximity of the telescope (Cerro this is not possible, we use the built-in acquisition camera of X- Paranal, Chile) to the South Atlantic Anomaly (SAA) affects the shooter to search for counterparts, and adapt the observing strat- rate of RRM triggers (Greiner et al. 2011). BAT switches off egy on the fly. In many cases we could interact directly with the when the satellite goes through the SSA and that means we have staff at the telescope to aid the observation. Raw data are usually fewer events immediately observable, compared to the rest of the available within minutes in the ESO archive and are promptly Swift orbit. reduced by members of our collaboration, often using archival calibration data which are readily available. This allows us in most cases to report the preliminary results (redshift, identifi- 2.4. Observations cation of the most prominent emission and absorption feature, The observations obtained for this sample have been secured etc.) within a few hours after the beginning of the observations with the cross-dispersed echelle spectrograph, X-shooter (Ver- and plan additional X-shooter observations if deemed necessary. net et al. 2011), mounted on one of two Unit Telescopes at ESO/VLT, UT2 (Kueyen) and UT3 (Melipal) during the 2.3. Rapid response mode (RRM) observations duration of this follow-up campaign. The observations have been taken during a period of eight years corresponding to Under rare circumstances, the use of the ESO rapid response the ESO observing periods P84 through P98 under the fol- mode (RRM; Vreeswijk et al. 2010) has been possible. The lowing programme IDs: 084.A-0260, 085.A-0009, 086.A-0073, RRM is a system to automatically override ongoing observations 087.A-0055, 088.A-0051, 089.A-0067, 090.A-0088, 092.A- at the telescope. This allows the shortest possible delay between 0124, 093.A-0069, 094.A-0134, 095.A-0045, 096.A-0079, the GRB trigger and the initiation of observations, where spec- 097.A-0036, and 098.A-0055 (PI: Fynbo) and 0091.C-0934 (PI: troscopic integration of the rapidly fading, optical transient can Kaper). These proposals were initiated on Guaranteed Time. We have included a few additional bursts, from the programmes 2 We use the updated values of Schlafly & Finkbeiner(2011) to correct 084.D-0265 (PI: Benetti), 091.A-0877 (PI: Schady), 092.D-0056 for foreground extinction, but the sample criterion is based on the old Schlegel et al.(1998) values for consistency with the first semesters of 3 These tags are distributed as part of the GCN notices. Explanations our program. for all tags are available online.

Article number, page 3 of 43 A&A proofs: manuscript no. XSGRB_sample_arxiv

0.0 0.1 0.2 0.3 0.4 0.5 AV

Fig. 1. Mollweide projection in Galactic coordinates of the full sky showing the positions on the sky of the bursts presented in this work. The equator is the Galactic plane. The empty stars/circles are the positions of all the 1266 Swift bursts detected until 31/03/2017. Stars indicate bursts with measured redshifts and circles indicate those without. Blue stars indicate the position of the 104 bursts fulfilling the sample criteria specified in Sect. 2.1 that have a measured redshift. Red outlines are added to the 84 GRBs that enter our statistical sample with both X-shooter spectroscopy and a measured redshift. Red outlines of empty stars represent bursts which have been observed with X-shooter and has a measured redshift, but is not a part of our statistical sample. The blue dots show the positions of the 61 GRBs of our statistical sample that lack redshift measurements. A red outline is added around the six bursts in our statistical sample which were observed with X-shooter, but did not yield a redshift measurement. Two empty circles with red outlines indicate bursts outside the statistical sample that were followed up with X-shooter, but without a redshift measurement. The different samples are compared in Sect.4. The background shows the dust maps presented in Schlegel et al.(1998). Please note that we removed the background where the sample criterion is violated (AV > 0.5 mag) and replaced it with a white background. The grey scale ◦ bar below indicates the value of AV on the plot. The dotted lines indicate intervals of 30 in longitude and latitude.

(PI: Rau), 092.D-0633, 098.A-0136 (PI: Greiner), and 095.B- X-shooter covers the spectral wavelength region from 300 0811 (PI: Levan). The total collection of spectra represents all nm to 2480 nm in a single exposure, by splitting the light GRB afterglows that have been followed up by X-shooter up into three separate spectroscopic arms through the use of two to 31/03/2017, which marks the end of the XS-GRB legacy dichroics. The ultraviolet blue (UVB) arm covers 300 - 550 nm, follow-up program. the visual (VIS) arm covers 550 - 1020 nm, and the near-infrared (NIR) arm covers 1020 - 2480 nm. For some of the observations, The first GRB followed up was GRB 090313 (de Ugarte we have applied a K-band blocking filter, cutting the coverage Postigo et al. 2010), observed on the 15th of March, 2009, dur- of the NIR arm at 2100 nm. The K-band blocking filter is only ing the commissioning of X-shooter on the second Unit Tele- used after 2012, where it was installed. This is done to reduce scope (UT2) of the VLT . The bursts observed during the com- the amount of scattered background light from the thermal in- missioning or the science verification process (GRB 090313, frared. For the majority of observations, a nodding observing GRB 090530, GRB 090809, GRB 090926A) are not a part of the scheme has been employed, with a nodding throw of 500. Each sample we use to investigate the statistical properties of GRB af- nodding observation has typically been carried out in a standard terglows, due to the different criteria for their selection. The first ABBA pattern. For some cases, conditions during the follow-up burst observed after science verification and the mounting of X- (either technical or weather), have necessitated alterations to this shooter on UT2, was GRB 091018, which is the first burst en- scheme as described in App.B. For RRM triggers, a slightly dif- tering our statistically homogeneous sample. For all bursts that ferent observing strategy was employed. Starting as rapidly as fulfill our sample selection criteria, described in Sect. 2.1, spec- possible, a simple stare mode sequence was started, with 5 spectroscopic follow-up have been attempted with X-shooter. Vari- troscopic integrations with increasing exposure times. ous conditions can affect our ability to follow up a given burst, and a discussion of these conditions and their consequences for For the majority of the bursts, we have observed with a slit the sample is included in Sect. 4.2.1. width of 100. 0, 000. 9, and 000. 9 for the UVB, VIS, and NIR-arm

Article number, page 4 of 43 Selsing et al.: The X-shooter GRB afterglow legacy sample respectively. This sets a lower limit on the delivered resolving the slit(see Sect. 3.4). This means that for these bursts, the initial power of the spectra based on the tabulated values of the deliv- sky-estimate was made from a limited number of pixels in the ered resolutions, which is 4350, 7450, and 5300 for the UVB, spatial direction. The subtraction of the sky background on the VIS and NIR-arm respectively4. For accurate measurements in- un-rectified image ensures that the bulk of the sky background is volving line profiles, knowledge of the precise instrumental res- not redistributed by the rectification process. olution is required. The spectral resolution becomes better than X-shooter is an echelle spectrograph, and therefore the indi- the nominal one, when the delivered seeing is smaller than the vidual echelle orders are curved across image space. The indi- projected width of the slit on the sky. We discuss how we deter- vidual orders therefore need to be rectified. In order to transform mine the effective instrumental resolution in Sect. 3.3. the image space (pixels) into a physical (wavelength-slit) space, Due to a mechanical failure, the atmospheric dispersion cor- the image pixels are resampled onto a physical grid, while prop- rector (ADC) was disabled from 1st of August 2012 until the end agating the pixel uncertainties derived by the pipeline. This rec- of this program. Only GRB 100728B was affected by the fail- tification process correlates neighboring pixels and in order to ing ADC prior to disablement, resulting in a lower-than-nominal minimize the degree of correlation, we need to choose a physical throughput. To avoid chromatic slit losses due to atmospheric sampling that matches the pixel sampling. We rectified the im- dispersion, nearly all subsequent observations have been carried age onto an equidistant grid with a dispersion sampling of 0.02 out at parallactic angle. A consequence of this is that for all ob- nm/pixel and a 000.16 per pixel spatial sampling for the UVB and servations following 1st of August 2012, the centroid of the trace VIS arm and 0.06 nm/pixel with a 000.21 per pixel in the NIR arm. of the source changes position across the spatial direction of the Because the tabulated resolution is a lower limit to the delivered slit as a function of wavelength. This effect has been modeled in resolution, we choose a sampling of 0.02 nm/pixel to ensure that the extraction procedure, as described in Sect. 3.4. the lowest wavelength part of neither of the arms have a sampling We provide an overview of all the observations in Table1 lower than the Nyquist sampling rate of 2 pixels per resolution and plot the positions of all the bursts on the celestial sphere FWHM. in Galactic coordinates in Fig.1. Away from the central zone of avoidance, due to high Galactic extinction cutoff (marked in white), the GRB positions have an isotropic distribution, except 3.1. Post-processing in the upper left quadrant which cannot be observed from Paranal due to the declination constraints of the telescope. For a typical observation, each of the exposures in the nodding sequence have been reduced as a single observation and then Thirty percent of the spectra presented here, primarily host subsequently combined to form a single image. We employ this observations, have already been published in Krühler et al. strategy so that we can reject outliers in the stack and weight by (2015). Single bursts have additionally been published, based an averaged measure of the inverse variance of the background. on unusual properties in their afterglows. We present individ- When weighting images, where the noise in each pixel is dom- ual notes on all the X-shooter spectra and their previous use inated by Poisson noise, it is important to estimate the back- in App.B. We include independent reductions of them here for ground variance in a large enough region, so that any correlation completeness. between the signal and the weights are removed. To this end, the weight map is generated by a running median window over the 3. Data processing variance map produced by the pipeline, where the trace has been masked and the width of the window is chosen to be wide enough In this section we describe how the final data products for the median variance to be calculated on the basis of several are produced and subsequently post-processed. All post- hundred pixels. This weighting scheme automatically also op- processing scripts developed for this dataset are made pub- timally combines images of different exposure times or images licly available at https://github.com/jselsing/XSGRB_ where the background is varying, which is often the case when a reduction_scripts. burst has been observed close to twilight. Before any reductions are initiated, the raw object images An additional sky-subtraction procedure is run on all recti- are run through the cosmic-ray removal algorithm (van Dokkum fied 2D spectra. This is done to remove residual sky, still present. 5 2001) implementation, Astro-SCRAPPY , where a wide clipping At each pixel in the dispersion direction, the spatial dimension is radius was used around detected cosmic ray hits to ensure that fit with a low-order polynomial, after all sources are masked out. edge residuals are robustly rejected. This low order polynomial is then convolved with a few pixel The basis for the reductions is the VLT/X-shooter pipeline, wide Gaussian filter in the dispersion direction and subtracted version 2.7.1 or newer (Goldoni et al. 2006; Modigliani from the entire 2D image. et al. 2010). The pipeline is managed with the Reflex interface In the NIR arm, where the background is very bright and (Freudling et al. 2013) and is used for subtraction of bias level, there are a high number of bright sky-lines, an alternative ap- flat-fielding, tracing of the echelle orders, wavelength calibra- proach to sky subtraction has been employed. When there are no tions with the use of arc-line lamps, flux calibration using spec- contaminating sources in the slit, the sky has been put back on trophotometric standards (Vernet et al. 2009; Hamuy et al. 1994), the images and the images are combined in pairs of two before mirror flexure compensation(see Sect. 3.2), sky-subtraction and extraction, subtracting the two from each other while keeping lastly the rectification and merging of the orders. Errors and bad the WCS static. Due to the nodding offsets used between obser- pixel maps are propagated throughout the extraction. For the ini- vations, this conserves the source flux while removing the sky at tial sky-subtraction, the background has been estimated by a run- √ the expense of a decrease in signal-to-noise by a factor of 2. ning median in regions adjacent to the object trace clear of con- This amounts to the regular nodding reduction, only we can re- taminating sources. Due to the broken ADC, for some objects ject outliers and weight by the averaged inverse variance map. there is curvature in the object trace along the dispersion axis of Reducing the images as single observations for all exposures, 4 See this URL for the nominal spectral resolutions. we additionally get a spectrum of the sky which we can use to 5 https://github.com/astropy/astroscrappy recalibrate the wavelength solution in the post-processing steps.

Article number, page 5 of 43 A&A proofs: manuscript no. XSGRB_sample_arxiv

3.2. Correction for offsets in the wavelength calibration NIR arm Since X-shooter is installed at the VLT Cassegrain focus, it is VIS arm prone to ﬂexures. The ﬂexures modify the projection of the slit 60 on the detector with respect to the one obtained in daytime cal-

) 50 1

ibration. This requires a modiﬁcation of the wavelength solu- − tion in order to correctly process the night-time data. Part of this km s correction is performed by the pipeline using the frames taken 40 during the X-shooter Active Flexure Compensation procedure6. 30

The remaining offset is corrected by cross-correlating the ob- Resolution ( served sky spectrum with a synthetic sky spectrum (Noll et al. 2012; Jones et al. 2013) after the continuum, estimated as the 20 mode of all flux values, has been subtracted. To get the correct seeing point-spread-function (PSF) with which to convolve the 10 synthetic sky, an initial refinement of the wavelength solution 0.0 0.5 1.0 1.5 2.0 2.5 has been obtained by cross-correlating the observed sky with an Spatial FWHM (”) unconvolved synthetic sky. This preliminary wavelength calibration is applied to the observed sky. The synthetic spectrum is Fig. 2. Blue data points show the FWHM (km/s) of Gaussian fits to un- then convolved with an increasing seeing PSF and the width that resolved telluric absorption lines in the VIS spectra, as a function of the minimizes χ2 with the updated observed sky is chosen to be the airmass corrected DIMM seeing. The red data points show a subsample effective sky-PSF. Using the synthetic sky with the matched res- of measurements obtained for NIR spectra. The colored lines show the corresponding linear fits to the data points. The colored bounds contain olution, a final wavelength calibration can then be calculated by 68 per cent of the best-fit probability mass. The red and blue dashed cross-correlating the observed sky with the correctly broadened lines indicate the nominal resolution for the slit-widths used for the tel- sky spectrum as a function of a velocity offset. Both a multiplica- luric standard observations. As can be seen, the effective resolution is tive and an additive offset to the wavelength calibration has been in many cases superior to the nominal one. tested, but in terms of χ2, the model with only a multiplicative offset is preferred. The resulting offsets, which were smaller than 0.01 nm in the UVB and VIS data and smaller than 0.05 nm in spectra, usually just after the science data, as part of the ESO the NIR spectra, but changing over short periods of time were X-shooter calibration plan. These spectra have been reduced to- applied to the corresponding spectra7. Using the convolved syn- gether with the afterglow spectra, using identical pipeline set- thetic sky, the pixels containing the brightest sky lines have been tings with the same version of the pipeline. To get the resolution flagged as such in the bad pixel map. for each observation, we select a series of atmospheric transitions that are resolved multiples which should be intrinsically unresolved, and are in areas with well defined continuum flux. 3.3. Spectral resolution We then fit the lines with Voigt-profiles and calculate a corre- The afterglow spectra described in this paper are obtained in sponding FWHM, which we can then convert to a delivered res- Target-of-Opportunity mode. In most cases, there is little possi- olution using the wavelength of the chosen transitions. To get the bility to tweak slit widths to the seeing at the time of observations resolution as a function of ambient conditions, for each obser- (i.e. to optimize spectral resolution and signal-to-noise), and al- vations we also calculate the airmass-corrected DIMM seeing, most all our data are therefore taken with a fixed set of slit widths which is measured at 500 nm. and binning, described above. In a fair number of cases, the see- The resulting distribution of spectral FWHM (km/s) as a ing full width at half maximum (FWHM) is considerably smaller function of spatial FWHM at 500 nm is fairly well described than the slit width, and the delivered spectral resolution will then in the VIS arm by a linear relation a + b ∗ x, with x the spatial be determined by the seeing rather than by the slit width, as af- FWHM in arcseconds, a = 25.2±0.8 km/s, b = 3.8±0.7 (see Fig. terglows are point sources (this is evidently not the case for ex- 2). For the NIR arm the corresponding relation is a = 36.0 ± 1.8 tended sources, e.g for host galaxies). The delivered resolution km/s, b = 8.4 ± 1.8. We use these linear relations as a way to es- for slit-width dominated spectra post-reduction and extraction timate the spectral resolution for medium to low signal-to-noise can easily be determined from the bright sky emission lines. For afterglow spectra. To extend this to the UVB arm, we calculate afterglow spectra with very high signal-to-noise, the delivered the ratio between the VIS and the NIR arm resolutions and find spectral resolution can at times be determined from the science that the resulting distribution is consistent with a simple scaling data themselves. However, in the presence of multiple velocity of the VIS arm relation by the ratio of resolutions of the NIR and components in absorption, other forms of line broadening, and VIS arm for unresolved, slit filling, sources as given on the ESO a lack of lines at some redshifts, this is difficult to do at lower instrument website. The UVB arm contains no suitable absorp- signal-to-noise ratios (the majority of spectra in our sample). A tion lines to use, and we therefore use a scaled value. This simple broad starting value for the expected resolution will help fitting analysis gives a sufficiently accurate estimate for the analysis of of these spectra, and can be important in upper limit determina- the low signal-to-noise science spectra. In all cases the deter- tion, and for this reason we construct a crude relation between mined resolution is written to the header with the "RESOLU- the seeing and the delivered resolution at our slit width, binning, TION" keyword. and reduction pipeline settings. To this end we use observations of telluric standard stars that are taken with identical instrument settings as our afterglow 3.4. Spectral extraction 6 X-shooter User Manual available at ESO website. To extract the afterglow spectrum from the rectified 2D-image, 7 The wavelength shifts have been extensively studied by ESO staff in several techniques have been employed based on the brightness this document. of the afterglow and the complexity of the objects entering the

Article number, page 6 of 43 Selsing et al.: The X-shooter GRB afterglow legacy sample slit. Due to the malfunctioning ADC (see Sect. 2.4), the spectral the update in Schlafly & Finkbeiner(2011) 8, and the extinction trace changes position across the slit in the spatial direction as curve by Cardelli et al.(1989) with a total to selective extinc- a function of wavelength. For a large fraction of the observed tion RV = 3.1. The wavelengths of the extracted 1D-spectra are bursts, using a single aperture for the spectral extraction is in- wavelength recalibrated (described in Sect. 3.2), moved to vac- adequate due to the large amount of background that would then uum, and corrected for barycentric motion. Pixels with pixel-to- enter the slit. To optimally select the extraction regions we there- pixel variation larger than 50σ are additionally added to the bad fore need to model the trace position. pixel map. To get the shape and the position of the spectral PSF as a function of location on the image, we need to chose a model 3.5. Telluric correction which can represent how the light falls on the slit. We know from Trujillo et al.(2001) that a Mo ffat function (Moffat 1969) with For all Earth-based telescopes, the light first has to pass through an index of β = 4.765 adequately describes an imaging PSF Earth’s atmosphere, where the atmospheric content and condi- due to atmospheric turbulence, but because of aberrations in the tions make an imprint on the received spectrum. These telluric optical dispersion elements and the rectification process, the PSF features can be corrected for in a multitude of ways. We employ we are trying to model is different from this profile. To allow for a prioritized list of methods here, depending on the availability flexibility in the model, we have chosen the Voigt function as a of the chosen methods. Since the observations are often taken at model for the spectral PSF and we describe how this is evaluated odd times under varying conditions, this prioritized list ensures in App.A. Since the host galaxy could also give a contribution that we are always doing the best possible correction. to the image profile, this choice allows for the required freedom The highest priority method is using the GRB afterglow con- if additional flux is in the wings of the profile. tinuum itself, where the atmospheric conditions have directly been imprinted on the spectrum. The telluric features can di- To guide the estimated position of the trace on the slit as a rectly be fit with an atmospheric model (Molecfit; Smette et al. function of wavelength, we have used the analytic prescription 2015; Kausch et al. 2015)9, which can then be used to correct for for the trace position described in Filippenko(1982), where the the absorption. The accuracy of the correction depends on the header keywords of the observations have been queried for the signal-to-noise per pixel of the target spectrum, where we have ambient conditions which controls the degree to which the trace chosen the requirement that the afterglow continuum spectrum changes position in the spatial direction. This analytic approach has a median signal-to-noise higher than a value of 10. is only valid for a plane-parallel atmosphere, but because the If the afterglow is not sufficiently bright, telluric standard final position is refined in the fit, it is adequate for our purposes. stars observed close in time to the GRB can be used as a proxy Based on the signal-to-noise of the afterglow continuum, the for the atmospheric condition during the GRB observation. Here 2D-image has been binned in the spectral direction to a num- we employ the telluric correction method that has been devel- ber of elements that allows for an accurate tracing of the PSF, oped in Selsing et al.(2015), where a library of synthetic tem- typically 200 bins for moderate signal-to-noise. For each of the plates is fit to the observed telluric standard. bins, using the analytically guided guess position, the spectral In the last case, where the object is neither bright enough, or PSF has been fit using the unweighted chi-squared minimiza- there for some reason a telluric standard have not been observed, tion algorithm implemented in scipy.optimize.curve_fit we rely on a synthetic sky model (Noll et al. 2012; Jones et al. (Oliphant 2007). Since we know that the trace varies slowly as 2013) for which we generate a synthetic transmission spectrum, a function of wavelength, we have then fitted a low-order poly- where the ambient parameters for the observations have been nomial to the fit parameters as a function of wavelength, which used. allows us to evaluate the spectral PSF at all wavelengths and in this way accurately model the entire spectral PSF. 3.6. Continuum estimate Equipped with a model for how the light is distributed across We additionally provide an estimate of the continuum for all the entire dispersion direction, we can employ the optimal ex- the spectra presented here. For this, we have developed an al- traction algorithm (Horne 1986), which weights the extraction gorithmic approach that attempts to automatically estimate the aperture by the spectral profile, or alternatively sum all pixels continuum placement along with the error on the continuum es- within 1 FWHM of the modeled profile. Where possible, we timate through an iterative procedure. The method is entirely have used the optimal extraction. In cases where the trace is very data-driven and does not rely on any physical assumptions. The weak, even in the binned images, an aperture has been selected method is applied on each arm separately for each spectrum, to manually which covers all emission lines, if present, and when allow the widest possible wavelength range of the spectral shape nothing is immediately visible, the entire nodding window. The to guide the normalization. error- and bad pixel maps are in all cases propagated throughout To estimate the continuum, a number of points (typically the extraction. on the order of 100) are inserted at random positions along the In cases where multiple traces are visible in the slit, addi- wavelength direction, and the flux-density of each point is deter- tional components for the profile are used in the optimal extrac- mined by the median value of the spectrum in a small region (∼ tion. The additional components do not share the PSF parame- 1 Å) surrounding each point. The points are fitted with a low or- ters and in cases where the additional component is an extended der polynomial (we use numpy.polynomial.chebyshev) and object, the fits have been inspected to ensure that the additional iteratively, the point furthest away is removed until the polyno- component does not skew the fit towards a different PSF. The mial fit differed from the points by less than ∼ 5 per cent. This fil- additional components are not used for the weights in the ex- tractions. 8 Queried from http://irsa.ipac.caltech.edu/applications/DUST/index.html using astroquery (Ginsburg et al. 2016). The spectra are corrected for Galactic extinction using the 9 http://www.eso.org/sci/software/pipelines/skytools/ E(B − V) value from the dust maps of Schlegel et al.(1998) with molecfit

Article number, page 7 of 43 A&A proofs: manuscript no. XSGRB_sample_arxiv

2 ∗ ∗ β α II II IV IV VI VI II iII iII S C S Ly Ly C O O Si Si S 1

0 3200 3400 3600 3800 4000 4200 4400 4600 4800 2 ∗ ∗∗ II IV IV iII iII S Fe Al C C S 1

0 5000 5200 5400 5600 5800 6000 6200 6400

2 ∗ ∗ ∗ II II II II II II II II Fe Fe Cr Zn Zn Fe Fe Fe 1

0 6600 6800 7000 7200 7400 7600 7800 8000 8200 2 I II II II II II II Mg Fe Fe Mg Mg Mn Mn

λ 1 F

0 8400 8600 8800 9000 9200 9400 9600 9800 ]

2 ] ] II II γ II II III H Ca Ca [O [O [O 1

0 10000 11000 12000 13000 14000 ] ] ]

2 ] ] ] II II I II β I II III III III H [S [O He Na Na [N MgH [O [O [O 1

0 15000 16000 17000 18000 19000

2 ] ] ] ] ] α II II II II II H [S [S [N [N [N 1

0 20000 21000 22000 23000 24000 Observed Wavelength (A)˚

Fig. 3. Telluric corrected, normalized spectrum of GRB 121024A at z = 2.300 that illustrates the typical data quality. The continuum estimate is shown in dashed red and the error spectrum in solid blue. The acquisition magnitude is R = 20, meaning it is in the brighter end of the sample presented here, but not the brightest. The spectrum is rich in absorption and emission lines, including absorption from molecular H2. The absorption trough visible at ∼ 4000 Å is due to Lyα in the host. We have marked the most prominent lines seen in GRB afterglows from Christensen et al. (2011). The regions of most severe telluric absorption are highlighted by grey-shading the background. Additionally, three intervening systems are seen in this sightline. This spectrum has been analyzed in detail in Friis et al.(2015). Article number, page 8 of 43 Selsing et al.: The X-shooter GRB afterglow legacy sample tering is used because the intrinsic afterglow continuum of GRBs – WAVE - Observed wavelength in vacuum, corrected for can be modeled by power laws (Piran 2005), and removing barycentric motion and drifts in the wavelength solution (Å). points that differ significantly from a smooth continuum shape – FLUX - Observed flux density (ergs−1cm−2Å−1). will guide the continuum estimate to a shape more reasonable for – ERR - Associated flux density error (ergs−1cm−2Å−1). GRB afterglows. The reason for the non-physical model for the – QUAL - Bad pixel map, where a value different from zero normalization is that it has the flexibility to capture instrumental indicates a bad pixel. variations of the continuum level that are not easily contained in – CONTINUUM - Continuum estimate based on Sect. 3.6 a more physically motivated model. Additionally, points spaced (ergs−1cm−2Å−1). closer than 1 per cent of the total spectral coverage are pruned. – CONTINUUM_ERR - Relative error on the continuum esti- The remaining filtered and pruned points are then spline in- mate. terpolated using scipy.interpolate.splrep.splev, which – TELL_CORR - Inverse transmission spectrum. Multiply serves as a first estimate of the continuum placement. To prevent FLUX and ERR column with this column to correct for tel- the spline from diverging at the edges, the spline-based contin- luric absorption. uum is tapered with a low order polynomial. An attempt to identify absorption and emission regions is then carried out, where they are marked as such if the difference between the estimated 4. Results continuum and the observed spectrum is larger than 3 - 5 times the associated error spectrum. All regions marked as affected are In this section, we describe the efficiency of the follow-up effort then masked. and the characteristics of the observed bursts. We also assess the degree to which the obtained sample is representative for the full Using the masked spectrum, the entire process is then re- Swift sample. An important note is that here we provide a release peated 500 times where the final continuum estimate is the mean for all GRBs after 14/03/2009, that have been observed with X- of the continuum realizations and the associated error estimate shooter, while only a subset of these constitutes our statistical is the standard deviation. This error reflects the stability of the sample. The statistical sample is based on the selection crite- algorithm across the spectrum. An example of the performance ria described in Sect. 2.1. Some bursts not fulfilling the sample α is shown in Fig.3. In the Ly forest, where there is very lit- criteria have been followed up due to interesting characteristics, tle flux at the continuum level, the performance of the nor- e.g., curios properties of their light curves, their brightness, etc. malization algorithm depends on the continuum coverage red- These bursts are not discussed as part of the investigation of the α wards of the Ly line. In some cases, very little continuum is statistical properties of the GRB population. A prime example contained in a single arm and a manual continuum estimate of a spectacular burst outside the statistical sample is the bright is provided, similar to what is done in López et al.(2016). In INTEGRAL burst GRB 161023A (de Ugarte Postigo et al., in these cases the continuum error is set to 10 per cent. The code prep), that contains at least 15 intervening absorption systems for the continuum estimate is released along with the paper at (See B.117). https://github.com/jselsing/XSGRB-sample-paper.

4.1. Follow-up timing and afterglow brightness 3.7. Science data products Redshift determination of bursts for which the host is too faint All the spectra are made available as a single ZIP file, through for a spectroscopic redshift measurement relies on the detection of absorption lines imprinted on the GRB afterglow continuum. http://grbspec.iaa.es (de Ugarte Postigo et al. 2014a), and −1 additionally through the ESO archive in the form of phase 3 Because the optical afterglow rapidly fades (typically as ∼ t ) material. This release includes both prompt afterglow observa- a rapid follow-up is essential. In Fig.4 we plot the delay from tions as well as late time observations of the associated hosts, the BAT trigger to the start of the spectroscopic observation. The and represents all afterglow spectra of GRBs carried out by the shortest delays are observed in RRM-mode. The fastest follow- X-shooter spectrograph since the commissioning of the instru- up between BAT trigger and start of spectroscopic observations ment, 14/03/2009, and until the end of the last period of the pro- for any observation is for the short, z = 1.717, GRB 160410A gram 098.A-0055, 31/03/2017 and thus constitutes eight years for which spectroscopic integration was initiated only 8.4 min- of GRB afterglow observations with X-shooter. An overview of utes after the BAT trigger. To illustrate the importance of the all the spectra and their observational setups is given in Table follow-up delay for the redshift completeness, we plot the red- 1. For each burst, each individual observation is provided in a shift completeness as a function of delay time in Fig.4 for all separate reduction, and in cases where observations have been the bursts we have followed up, including the ones outside the repeated for an increased signal-to-noise or to follow the tempo- statistical sample. As can be seen from the figure, the fraction of ral evolution, a combined spectrum is also provided. No attempt GRBs with a redshift determination decreases with follow-up de- has been made to join the spectroscopic arms, so for each obser- lay. The redshift completeness for bursts that we have followed vation, three spectra are provided in separate files. up is 94 per cent. This degree of completeness in the followed bursts, illustrates the efficiency of VLT/X-shooter in redshift de- All spectra are released in the ESO Science Data Prod- termination. Not shown in the figure are an additional 12 bursts uct (SDP) format (Micol et al. 2016), and formatted as bi- that have redshift determinations based on late-time host obser- nary FITS files. The naming convention is based on the GRB vations with delay times longer than ∼10 days. name and the observation number, and follow the scheme GRBxxxxxxx_OBxarm.fits. For example, the visual arm of the third observation of GRB 151021A, observed in RRM mode (see 4.2. Sample completeness GRB151021A_OB3VIS.fits Sect. 2.3), is named . Of all the BAT-triggered bursts, a total of 165 bursts fulfill the Each file contains 7 columns with the following contents and sample criteria specified in Sect. 2.1, since the commissioning descriptions: of VLT/X-shooter. This sample constitutes the statistical sample

Article number, page 9 of 43 A&A proofs: manuscript no. XSGRB_sample_arxiv

vided in Table2. It can be seen from the values that the three Redshift 1 2 3 4 5 samples have very similar distributions in terms of the point estimates chosen. This suggests that our selection criteria are unbiased compared to the Swift-sample and that additionally, the 16 follow-up eﬀort conserves the distributions of the intrinsic GRB 100 properties (except perhaps for NHI,X, see Sect. 4.4). 17 Additionally, using a 2-sided Kolmogorov-Smirnov test (KS- -band) 95 test), we can assess the degree to which the null hypothesis, that R, i 18 the two distributions are drawn from the same parent distribu- 19 tion, can be rejected. We show a graphic representation of the 90 test statistics in Fig.6. A high p-value indicates little evidence 20 against the null hypothesis . The distribution of NHI,X exhibits 85 the highest degree of dissimilarity, but the two distribution are 21 still consistent with being drawn from the same underlying distribution. 80 22 Redshift completeness (%) We therefore conclude that the statistics of the sample presented here, conserves the intrinsic properties of the GRBs in 23 75 the full Swift sample - at least in terms of BAT ﬂuence and X-ray

Acquisition camera24 magnitude ( ﬂux at 11 hours.

0.1 1.0 10.0 100.0 4.2.1. Properties of rejected triggers Follow-up delay (hours) Out of the 165 bursts meeting our initial selection criteria for the statistical sample, 36 (22 per cent) were not observed due to rea- Fig. 4. Afterglow magnitude at the start of observation and redshift completeness as a function of follow-up delay for all the afterglows that have sons unrelated to the GRB or afterglow properties. The reasons been followed up. The points have been colored based on the redshift include unavailability of the telescope due to technical mainte- of the corresponding burst. Red symbols indicate GRBs without a mea- nance (e.g., mirror re-coating), a visiting observer rejecting the sured redshift and arrows indicate bursts for which the afterglow was ToO trigger, or bad weather. Because this cut is unrelated to the not detected in the acquisition image. In red is shown the redshift com- GRB properties, it will not change the statistical properties of pleteness as a function of follow-up delay. the full sample. Removing these bursts from the statistical sample, dramatically improves the redshift completeness from 61 per cent to 88 per cent. The remaining burst not followed up from which we will derive statistical properties of the GRB pop- already had a redshift from other instruments or were very faint ulation. The redshift completeness of the full statistical sample and without a host association, thus observations were unlikely is 61 per cent. We return to the question of redshift completeness to yield a redshift measurement. In the remainder of the text, we in Sect. 4.3. From this sample, 93 GRBs have been spectroscop- consider the 129 bursts our statistical sample. ically followed up with X-shooter. In order to assess whether the subset of bursts followed up are representative of the underlying GRB parent population, we compare intrinsic properties of 4.3. On the redshift distribution of GRBs GRBs in our sample to GRBs in the full sample followed up by One of the objectives of our follow-up campaign is to measure Swift. We show the comparison between the BAT (15 – 150 keV) the redshift distribution for a well-defined, observationally un- fluence, the XRT flux (0.3 – 10 keV) at 11 hours, and the in- biased and statistically useful sample of GRBs. The imposed trinsic X-ray derived equivalent hydrogen column density at the selection criteria (see Sect. 2.1) ensure that the GRBs entering redshift of the GRB, in excess of the Galactic X-ray absorption our homogeneous sample, fairly represent the underlying popu- column, NHI,X, in Fig.5. For the latter, we can only use values lation. The redshift distribution of such a sample holds valuable of NHI,X derived for bursts with a measured redshift, excluding information about the occurrence of GRBs through cosmic time ∼ 75 per cent of the full Swift sample. We return to the last point (Jakobsson et al. 2012; Perley et al. 2016a). in Sect. 4.4. Using the observational characteristics of the 1266 bursts observed until 31/03/2017 by Swift, and the derived NHI,X (Evans Sample Nbursts zcompleteness zmean zmed +1.83 et al. 2009), we can quantify the degree to which our sample is XS-GRB 129 88 % 1.89 1.52−0.91 biased relative to the overall Swift sample. The values are queried SHOALS 119 99 % 2.18 2.06+1.27 from the online Swift database10,11. Three samples are of inter- −1.20 BAT6 58 97 % 1.90 1.71+1.31 est in order to assess the completeness of the follow-up cam- −1.04 . +1.42 paign (Fig.5); the full Swift sample consisting of all the bursts TOUGH 69 87 % 2.20 2 11−1.36 +1.28 observed by Swift (blue), all the bursts that fulfill the selection Fynbo09 146 49 % 2.2 2.1−1.23 criteria imposed in Sect. 2.1 (green), and the bursts actually fol- Table 3. Comparison between the redshift distributions of previous lowed up with X-shooter (red). complete samples. The SHOALS redshift characteristics are taken from For each of the samples, we calculate the median, 16th, and Perley et al.(2016a), the BAT6 redshifts are from Salvaterra et al.(2012) 84th percentiles of each of the distributions, which can be used with the update from Pescalli et al.(2015), TOUGH is from Hjorth et al. as point estimates for the population distribution. These are pro- (2012) with the update from Schulze et al.(2015), and Fynbo09 are from Fynbo et al.(2009). The errors shown on the median redshift con- 10 http://swift.gsfc.nasa.gov/archive/grb_table/ tain 68 per cent of the probability mass. 11 http://www.swift.ac.uk/xrt_live_cat/

Article number, page 10 of 43 Selsing et al.: The X-shooter GRB afterglow legacy sample

Swift 1.0 1.0 Statistical 1.0 Observed

0.8 0.8 0.8

0.6 0.6 0.6

0.4 0.4 0.4

0.2 0.2 0.2

0.0 0.0 0.0

-8.0 -6.0 -4.0 17.5 15.0 12.5 10 15 20 2 − − − 2 1 2 log(15 150 keV Fluence/erg cm− ) log(0.3 10 keV Flux/erg cm− s− ) log(NH,X/cm− ) − −

Fig. 5. Comparison between the burst properties of all bursts observed with Swift (blue), the statistical sample that fullfill the criteria specified in Sect. 2.1 (green), and the subset that has been observed as part of the statistical sample (red). The left plot shows the fluence in the 15-150 keV band as observed by BAT. The middle panel shows the 0.3 - 10 keV flux, 11 hours after the bursts as measured by XRT. In the right most panel, we show NHI,X based on the XRT spectrum (Evans et al. 2009).

Full Swift sample Statistical sample Followed up bursts

NBAT 981 163 92 − − +0.7 − +0.7 − +0.7 log(15 150 keV fluence) 5.9−0.6 5.9−0.6 5.9−0.7 NXRT 902 160 90 +0.7 +0.7 +0.9 log(0.3 − 10 keV flux) −12.3−0.8 −12.4−0.8 −12.4−0.7 NHIx 248 99 79 +0.6 +0.7 +0.7 log(NHIx ) 21.7−0.9 21.5−3.4 21.6−4.5 Table 2. Population properties (median and 16th and 84th percentiles as the error intervals) for the Swift sample and the subset of bursts fulfilling the sample criteria. The population characteristics of the three samples are very similar, which shows that our selection criteria effectively conserve the statistical properties of the underlying population, as least for these parameters. Notice that not all bursts have measurements of the quantities we compare.

In Fig.7 we show the redshift distribution of all the observed anism of Swift, this sensitivity limit should be interpreted with GRBs. In the top panel we show a histogram for the full sample some caution as the effective limit depends on the light curve of and the statistical sample and in the main panel, we show the the prompt emission signal. Due to the dilution of light with dis- redshifts of the individual bursts as a function of GRB energy in tance, the Swift GRB luminosity detection limit is almost an or- the observed 15 – 150 keV band. To calculate the energy, EBAT, der of magnitude brighter at z = 2 than at z = 1. At z ≥ 4 we are we follow a similar procedure as Lien et al.(2016) and define only able to observe GRBs that are ∼ hundred times brighter than 2 −1 EBAT = Fγ 4π dL (1 + z) , where Fγ is the observed BAT flu- the faintest bursts at z = 1 and below. The effect of GRB redshift ence in the 15-150 keV band and dL is the luminosity distance on the Swift triggering criteria have previously been studied in to the burst at the given redshift. Note that this measure of lumi- detail (Littlejohns et al. 2013). nosity does not include any k-correction. As an indication of the effect of the Swift sensitivity limit on the redshift distribution, We compare the point estimates for the redshift distributions also in the figure, we have shown the so called ∼ 1 s flux BAT of previous complete samples of GRBs in Tab.3. We see that sensitivity limit (∼ 3 × 10−8 erg s−1 cm−2; Baumgartner et al. when we compare to other complete samples, the XS-GRB pre- 2013; Lien et al. 2016). Due to the complex triggering mech- sented here has the lowest average redshift. However, the other samples also exhibit a large spread in the redshift distributions.

Article number, page 11 of 43 A&A proofs: manuscript no. XSGRB_sample_arxiv

that the extinction of the optical afterglows is primarily driven by the presence of dust in the host galaxies and not solely by unfor- tunate placement of the synchrotron spectral break frequencies. Hjorth et al.(2012) find that systems with no optical afterglow have higher NH,X, irrespective of the nature of the host – which however also turn out to be redder. Additionally, the ISM absorption lines in dark sight-lines are found to be stronger compared to optically brighter bursts (Christensen et al. 2011), which is consistent with the dark bursts being found in more metal-rich and dustier galaxies. For all bursts with follow-up within 100 hours we calculate the "darkness"-parameter, βOX (Jakobsson et al. 2004). This requires the simultaneous measurement of the X-ray flux density and the optical flux density which is in practice possible, but in Fig. 6. Relational graph showing the respective p-values. They all rep- reality extremely rarely available. As a proxy, we use the mea- resent the degree that the different samples are drawn from the same sured acquisition camera magnitude reported in Table1 to get the underlying distribution. The arrows represent the comparison direction with each of the p-values they are listed next to. Only in the HI column optical flux density at the beginning of the spectroscopic integra- density distribution is there mild evidence against the null hypothesis, tion. Because we know the delay between the follow-up and the but the discrepancy is mainly driven by a relatively larger fractional Swift trigger, we can use the measured XRT lightcurve (Evans contribution from low-column hosts in the statistical sample. et al. 2007, 2009)12 to infer the corresponding X-ray flux density at the time of the optical observation. This is done by either linearly interpolating between temporally neighboring XRT A 2-sided KS test reveals that the XS-GRB sample is consis- measurements or by extrapolating the last few X-ray data points tent with being drawn from the same parent sample with the fol- to the time of the spectroscopic observation. When the afterglow lowing p-values: SHOALS (p-value = 0.13), BAT6 (p-value = is not detected in the acquisition camera, an upper limit of > 24 0.95), Fynbo09 (p-value = 0.08) and TOUGH (p-value = 0.09). mags is used, which propagates into an upper limit on βOX. As a small note, the redshift distribution of BAT6 are not ex- In Fig.8 we compare the βOX - NH,X distribution with the pected to be identical to the other complete samples due to the one presented in Fynbo et al.(2009). We take the NHI,X val- additional cut on the GRB peak flux in the BAT band. ues from the XRT spectral fits(Evans et al. 2009). The values Because the redshift completeness of our statistical sample is from Fynbo et al.(2009) have been treated as detections, mean- 88 per cent, making an inference of the true redshift distribution ing that we artificially bias the distribution towards higher βOX- of GRBs based on this sample is impossible. For instance, only values. The two distributions exhibit a large degree of overlap. the brightest GRBs are seen above redshift z & 1 as shown in We confirm the result by Fynbo et al.(2009), that dark bursts, −2 Fig.7. As described in detail in Hjorth et al.(2012) and Perley βOX < 0.5, have higher log(NH,X/cm ). Specifically, for bursts et al.(2016b), bursts for which the redshift is measured from the with measured redshift either from the afterglow or the host afterglow are systematically found in host galaxies with a lower galaxy, we find the following: For bursts with βOX ≥ 0.5 we luminosity than bursts for which the redshift is measured from −2 +0.7 find log(NH,X/cm ) = 21.4−1.0, whereas for βOX < 0.5 we find the host galaxy. Only a few GRBs hosted in galaxies, with stellar −2 +0.5 10 log(NH,X/cm ) 21.8−0.9 where 68 per cent of the probability masses more than 10 M have the redshift measured based on mass is contained within the error intervals. A 2-sided KS test the afterglow continuum. This is likely related to the presence fails to reject the null hypothesis that they are drawn from the of higher contents of dust in more massive galaxies, leading to a same distribution with p = 0.11, meaning that there is no strong larger fraction of extincted afterglows. evidence for a discrepancy. A Kendall’s τ test, however, suggest a statistically significant, low degree of negative correlation 4.4. Sample darkness (Γ = −0.21 at a p-value = 0.01). Using the table maintained by J. Greiner13, we can see how A fraction of all GRBs exhibit no detectable or very faint opti- the presence of an optical afterglow affects the follow-up statis- cal afterglows (Groot et al. 1998; Djorgovski et al. 2001; Fynbo tics. 50.5 per cent of all Swift-triggered bursts in this list do et al. 2001). The degree of optical extinction relative to the X- not have a detected optical afterglow. This number also includes ray brightness has been parametrized in terms of their optical bursts where no optical observations were available, so the real darkness, using the measurement of, or limit on, the optical to number is likely to be lower. For the bursts that enter our statis- X-ray spectral index βOX (Jakobsson et al. 2004; Rol et al. 2005; tical sample, the percentage of bursts without a detected optical van der Horst et al. 2009). The X-ray properties of such bursts afterglow is 28 per cent, close to the upper limit on the fraction of have previously been investigated (De Pasquale et al. 2003; dark bursts found in a complete sample with a very high degree Fynbo et al. 2009; Melandri et al. 2012) and there are some indi- of redshift completeness (Melandri et al. 2012). Of the bursts cations that dark bursts have somewhat higher X-ray luminosity for which follow-up has actually been attempted, this number and NH,X compared to the optically bright bursts (Campana et al. is 23 per cent, suggesting a slight bias against bursts without a 2012; Watson & Jakobsson 2012). The X-ray column density has detected optical afterglow in the spectroscopic sample. been shown to be roughly correlated with the gas column den- However, the fraction of dark bursts for which we have mea- sity, which, for a given metallicity, also correlates with the dust sured redshifts is lower than the ones with a detected optical extinction (Watson et al. 2013; Covino et al. 2013a), though the afterglow. For the afterglows we have observed as part of the range in metallicity introduces a large additional scatter in the statistical sample that do not have a detected optical afterglow, correlation between the extinction and the NH,X. This indicates along with investigations of host galaxy properties (Greiner et al. 12 http://www.swift.ac.uk/xrt_curves/trigger_numer/flux.qdp 2011; Krühler et al. 2011; Hjorth et al. 2012; Perley et al. 2016b), 13 http://www.mpe.mpg.de/ jcg/grbgen.html

Article number, page 12 of 43 Selsing et al.: The X-shooter GRB afterglow legacy sample

Lookback time [Gyr]

54 0 10.5 12.3 12.9 13.2 13.3

GRB 090423 53 z = 8.2

GN−z11 z = 11.1 EGSY8p7 52 ULAS J1342+0928 z = 8.7 z = 7.54 QSO / erg) BAT

−2 51 −8 erg cm )

log (E Swift/BAT detection limit (3 x 10

Sample XS−GRBs

Epoch of reionization All XS−GRBs 49 0 2 4 6 8 10 Redshift (z)

Fig. 7. Redshift distribution as a function of intrinsic BAT γ-ray energy, EBAT. Bursts that are a part of the statistical sample are marked by blue stars whereas black dots show all GRBs observed with X-shooter. All Swift GRBs with measured redshifts are shown in gray. For comparison, we overplot with red stars the GRB (Tanvir et al. 2009; Salvaterra et al. 2009), quasar (Bañados et al. 2017), and galaxies (Zitrin et al. 2015; Oesch et al. 2016) with the highest spectroscopically conﬁrmed redshifts, the latter three shown at arbitrary EBAT. The blue solid line represents the so-called 1 s BAT sensitivity limit described in the text. The estimated epoch of reionization is shown by the black dot-dashed line, with the uncertainty shown as the black dashed lines, from the most recent measurement by Planck Collaboration et al.(2016). On the top of the plot is shown the marginalization of the redshift. Again the blue histogram represents the bursts than enter our sample and the white histogram the full GRB population.

the redshift completeness is 53 per cent. For comparison, for the 2003; Nysewander et al. 2009), if the measured NH,X is primarily afterglows in the statistical sample we have observed with an driven by the gas column in the neighborhood of the burst. Be- optical afterglow detected, the redshift completeness is 92 per cause we only use values for NH,X in the comparison for which cent. This also shows that the lack of redshift completeness in the GRB has a measured redshift, this measure is likely biased the sample presented here is in part due to the increased diffi- toward optically brighter bursts (e.g. Watson & Jakobsson 2012). culty of measuring a redshift for bursts without an afterglow. To measure a redshift, we either need a detected afterglow to obtain a spectrum or to locate the host galaxy and determine the 4.5. Hydrogen column densities redshift from there. It is more difficult to correctly associate a galaxy with a burst when there is no detected optical afterglow The locations of long GRBs are associated with intensely star- and hence a correct redshift measurement is more difficult to forming regions (Hogg & Fruchter 1999; Bloom et al. 2002; make, see Jakobsson et al.(2005); Levesque et al.(2010) and Fruchter et al. 2006; Lyman et al. 2017). Because a significant Perley et al.(2017). fraction of the hydrogen along the line of sight in these regions has not been ionized, the optical depth at the wavelength of Regardless the fraction of dark bursts being lower in the Lyα is very high, saturating not only the line center, but also the observed sample, compared to the statistical sample, the X-ray damping wings. This causes a strong absorption system from the properties do not differ significantly, as shown in Sect. 4.2. This Lyα-transition to appear in the afterglow continuum. For bursts is despite spectroscopic follow-up only being carried out in cases with z & 1.7, the position of Lyα moves into the spectroscopic where either a detectable optical afterglow or a clear counterpart coverage of X-shooter, meaning that we are able to detect this are seen, which naively should be biased against dark bursts oc- absorption trough due to Lyα. The exceptionally good UV re- curring in more obscured galaxies, which is shown to exhibit dif- sponse of X-shooter, allows us to robustly measure Hi-columns ferent galactic properties (Perley et al. 2009; Krühler et al. 2011; - also at these relatively low redshifts. Rossi et al. 2012; Perley et al. 2013, 2015). That the decreased Due to the stochastic nature of the Lyα-forest, the blue fraction of dark bursts in observed sample does not alter the ob- wing of the Lyman-α absorption line is randomly superposed served prompt X-ray brightness distribution, potentially reflects with Lyman-α forest systems, along with strong absorption from the independence of the X-ray brightness on the density of the Mn ii and Si iii, making it notoriously difficult to model. Addi- circumburst medium (Freedman & Waxman 2001; Berger et al. tionally, the red wing has ISM signatures imprinted on it, es-

Article number, page 13 of 43 A&A proofs: manuscript no. XSGRB_sample_arxiv

24.0 1.2 Fynbo et al. 2009 Tanvir et al. 2017 This work 23.5 Noterdaeme et al. 2012b This work 1.0 23.0 ) 0.8 2 22.5 − cm / OX 22.0 β 0.6 HI N

log( 21.5 βOX = 0.5 0.4 21.0 0.2 20.5

0.0 20.0 20 21 22 23 1 2 3 4 5 6 7 2 log(NH,X/cm− ) Redshift

Fig. 8. βOX against X-ray derived hydrogen column density. In red is shown the sample presented in Fynbo et al.(2009) where the lines indi- Fig. 9. Distributions of hydrogen column densities for absorbers found cate the kernel density estimate of the distribution. For the kernel den- in quasar absorption lines, from Noterdaeme et al.(2012) in blue. Over- sity estimate, the limits have been replaced with values. We note that plotted in green is the kernel density estimate of absorbers in GRB sight- contrary to what is presented in Fynbo et al.(2009), we exclude all lines. Values are taken from the compilation in Tanvir et al. (submitted), bursts for which we do not have a redshift. Darker colors represent a along with the new values presented in this sample. We also show in the higher density of points. In blue are the points for the bursts presented red the values derived in this work. The marginal distributions for the here along with the marginal histograms. Limits on βOX are shown by three samples are also shown along the left side and on the top, where downwards facing arrows. The corresponding marginal distributions are the different environments probed are clearly visible in the hydrogen shown along the edges of the plot. A Kendall’s τ test, indicates a sta- column densities, as previously also noted in Fynbo et al.(2009). tistically significant, low degree of negative correlation (Γ = −0.21 at a p-value = 0.01) dian, the 16th, and 84th percentiles of the two distributions. The −2 +0.3 sample presented has a log(NHI/cm ) = 21.8−0.8 and the rest −2 +0.4 pecially strong absorption due to Si ii,S ii and N v, which can of the literature values has log(NHI/cm ) = 21.5−1.5. We see exhibit significant velocity structure. Along with instrumental that the two distributions have a large degree of overlap due to effects, the generative model for the data that we would use in the large width of the distributions, but we find a slightly higher a likelihood-based analysis would be very complicated. We have median value for the new sample presented here. A 2-sided KS therefore decided not to make formal χ2 fitting of the hydrogen test gives a p-value of p = 0.006, meaning relatively strong evi- column densities, but instead use a more subjective visual mea- dence against the null hypothesis that the two samples are drawn surement to the absorption profile. Using an analytic approxi- from the same underlying distribution. Because the bursts that mation to the absorption profile from Tepper García(2006), we have measurements of the hydrogen column density are selected overplot a synthetic absorption line with a specified column den- solely based on our ability to infer a column, it is difficult to sity on our observed, normalised spectrum. By tuning the value make any strong conclusions about the population statistics in of the hydrogen column density until the synthetic absorption terms of gas content. line matches the spectrum, we can thereby infer the actual col- In Fig.9, we also show the column density distribution for −2 umn density of the GRB sight line in a manual way. Similarly, the 12081 quasar absorbers with log(NHI/cm ) > 20 from No- the uncertainty on the hydrogen column can be estimated by ad- terdaeme et al.(2012). The fact that GRBs are systematically lo- −2 justing the error, until the confidence bounds contain the con- cated behind the highest log(NHI/cm ), previously noted (e.g., tinuum variation. We show the results of this procedure for all Prochaska et al. 2007; Fynbo et al. 2009), is very clear in this fig- bursts where possible in Fig. 10 and the inferred hydrogen col- ure. The reason for this is that quasar sample sight-lines through umn densities in Table4. galaxies that are cross-section selected, whereas GRB sight-lines probe the dense, star-forming regions in their hosts. 12 of the NHI measurements for these spectra have previously been presented in Cucchiara et al.(2015) (See Table4). We provide new measurements here for completeness. In the compila- 5. Discussion and conclusions tion of NHI measurements towards GRBs in Tanvir et al. (submitted) there are 93 published NHI values, excluding the measure- In this paper we have presented the results of a dedicated effort ments provided here. We here provide 30 new neutral hydrogen over the years 2009 – 2017 to use the X-shooter spectrograph on column density measurements – an increase of the number of the ESO-VLT to secure spectroscopic observations of afterglows optically derived hydrogen column densities of ∼ 33 per cent. and host galaxies of GRBs detected by Swift. This work was We show the two distributions in Fig.9. We compare the me- initiated by the consortium that built X-shooter and included this

Article number, page 14 of 43 Selsing et al.: The X-shooter GRB afterglow legacy sample project as a key part of the GTO program, but over the years the Table 4. Hydrogen column densities for the 41 bursts exhibiting Lyα ab- project continued in open time. sorption in the spectral coverage of X-shooter. The corresponding measurements are shown in Fig. 10. The sample presented here includes spectroscopic observations of 93 systems fulfilling our sample criteria, including 18 spectra that are late-time observations of the underlying host GRB Hydrogen Column −2 galaxies. All spectra have been made publicly available in the log(NHI/cm ) reduced form used in this paper. GRB 090809Ac 21.7 ± 0.2 Our sample serves the purpose to further characterize the en- GRB 090926A 21.55 ± 0.10 vironments of GRBs that was also much advanced by the previ- GRB 100219Ac 21.2 ± 0.2 ous surveys based primarily on lower-resolution spectroscopy. GRB 100425Ac 21.0 ± 0.2 GRB afterglow sight-lines are unique in the sense that only GRB 100728B 21.2 ± 0.5 after observing more than 12000 damped Lyman-α absorbers GRB 110128A 21.90 ± 0.15 (DLAs) towards about 105 quasars, a handful systems with −2 GRB 110818A 21.9 ± 0.4 log(NHI/cm ) > 22 have been identified (e.g., five in Noter- GRB 111008Ac 22.40 ± 0.10 daeme et al. 2012). Long GRB afterglow spectra, by contrast, GRB 111107Ac 21.0 ± 0.2 reveal such systems routinely (Jakobsson et al. 2006; Fynbo et al. GRB 120119A 22.6 ± 0.2 2009; Cucchiara et al. 2015, and this work). With afterglow spec- GRB 120327Ac 22.00 ± 0.05 troscopy (throughout the electromagnetic spectrum from X-rays GRB 120404A 20.7 ± 0.3 to the sub-mm) we are able to characterize the properties of GRB 120712A 19.95 ± 0.15 star-forming galaxies over cosmic history in terms of redshifts, GRB 120716Ac 21.80 ± 0.15 metallicities, molecular contents, ISM temperatures, UV-flux GRB 120815Ac 22.10 ± 0.10 densities, extinction curves, etc. A number of independent pa- GRB 120909Ac 21.75 ± 0.10 pers have been published or submitted for publication focusing GRB 121024Ac 21.85 ± 0.10 on many of these specific issues of our sample such as extinction GRB 121027A 22.8 ± 0.3 curves (Japelj et al. 2015, Zafar et al. submitted, see also Fynbo GRB 121201Aa c 22.0 ± 0.3 et al. 2014; Heintz et al. 2017), emission lines from the under- GRB 121229A 21.8 ± 0.2 lying host galaxies (Krühler et al. 2015), the frequency of inter- GRB 130408Ac 21.80 ± 0.10 ii vening Mg absorbers (Christensen et al. 2017), Arabsalmani GRB 130427B 21.9 ± 0.3 et al.(2018) on the metallicity-scaling relations, and escape of GRB 130606A 19.91 ± 0.05 ionizing radiation (Tanvir et al., submitted). A number of addi- GRB 130612A 22.1 ± 0.2 tional companion papers are also planned, investigating the de- GRB 131011A 22.0 ± 0.3 tailed properties of the sample presented here, including equiv- GRB 131117A 20.0 ± 0.3 alent width distributions (de Ugarte Postigo et al., in prepara- GRB 140311A 22.40 ± 0.15 tion), metallicities and kinematics (Thöne et al., in preparation), GRB 140430A 21.8 ± 0.3 high ionization lines (Heintz et al., in submitted), molecular lines GRB 140515A 19.0 ± 0.5 (Bolmer et al., in preparation), fine-structure lines (Vreeswijk GRB 140614A 21.6 ± 0.3 et al., in preparation), and composite GRB afterglow spectrum GRB 141028A 20.60 ± 0.15 (Selsing et al., in preparation). GRB 141109A 22.10 ± 0.10 The potential of using GRB sightlines as probes is far from GRB 150206A 21.7 ± 0.4 fully harvested. The sample of sightlines probed by our spectra GRB 150403A 21.8 ± 0.2 are not representative for all GRB sightlines as we have shown GRB 150915Aa 21.2 ± 0.3 and consistent with earlier findings from samples based on low- GRB 151021Aa 22.3 ± 0.2 resolution spectroscopy (e.g., Fynbo et al. 2009) and from stud- GRB 151027B 20.5 ± 0.2 ies of complete samples of GRB host galaxies (Hjorth et al. GRB 160203A 21.75 ± 0.10 2012; Covino et al. 2013b; Perley et al. 2016a). Krühler et al. GRB 160410Ab 21.2 ± 0.2 (2013) argue, that very rich sightlines like that probed by the GRB 161014A 21.4 ± 0.3 remarkable GRB 080607 (Prochaska et al. 2009; Sheffer et al. GRB 161023A 20.96 ± 0.05 2009; Perley et al. 2011) are probably significantly more fre- GRB 170202Aa 21.55 ± 0.10 quent than in the sightlines sampled by our spectra. However, with current instrumentation, these sightlines are out of reach Notes. (a) α (b) (c) except under very fortunate circumstances as in the case of Has Ly emission in the trough. Short burst. Previously published in Cucchiara et al.(2015). GRB 080607 when the afterglow could be observed only a few minutes after the burst with a 10-m class telescope. Observations of such sightlines with X-shooter-like spectrographs on the next port from NOVA and NWO-FAPESP grant for advanced instrumentation in as- generation of 20–40-m telescopes is likely to be very rewarding, tronomy. KEH and PJ acknowledge support by a Project Grant (162948–051) given that a suitable GRB detector will be available. from The Icelandic Research Fund. AG acknowledges the financial support from the Slovenian Research Agency (research core funding No. P1-0031 and project Acknowledgements. Much of the work done here would not have been done grant No. J1-8136). CT acknowledges support from a Spanish National Research without Neil Gehrels, who has now passed away. We owe him a large debt for Grant of Excellence under project AYA 2014-58381-P and funding associated his invaluable work with the Swift satellite. We also want to remember Javier to a Ramón y Cajál fellowship under grant number RyC-2012-09984. AdUP Gorosabel, who was taken from us too early. JPUF, BMJ and DX acknowledge acknowledges support from a Ramón y Cajal fellowship, a BBVA Foundation support from the ERC-StG grant EGGS-278202. The Dark Cosmology Centre Grant for Researchers and Cultural Creators, and the Spanish Ministry of Econ- was funded by the Danish National Research Foundation. This work was sup- omy and Competitiveness through project AYA2014-58381-P. ZC acknowledges ported by a VILLUM FONDEN Investigator grant to JH (project number 16599). support from the Spanish research project AYA 2014-58381-P and support from TK acknowledges support by the European Commission under the Marie Curie Juan de la Cierva Incorporación fellowships IJCI-2014-21669. DAK acknowl- Intra-European Fellowship Programme in FP7. LK and JJ acknowledges sup- edges support from the Spanish research project AYA 2014-58381-P and sup-

Article number, page 15 of 43 A&A proofs: manuscript no. XSGRB_sample_arxiv port from Juan de la Cierva Incorporación fellowships IJCI-2015- 26153. RSR Christensen, L., Vergani, S. D., Schulze, S., et al. 2017, A&A, 608, A84, acknowledges AdUP’s BBVA Foundation Grant for Researchers and Cultural arXiv:1709.01084 Creators and support from ASI (Italian Space Agency) through the Contract n. Cobb, B. E., Bloom, J. S., Perley, D. A., et al. 2010, ApJ, 718, L150 2015-046-R.0 and from European Union Horizon 2020 Programme under the Covino, S., Melandri, A., Salvaterra, R., et al. 2013a, MNRAS, 432, 1231, AHEAD project (grant agreement n. 654215). GL is supported by a research arXiv:1303.4743 grant (19054) from VILLUM FONDEN. SDV acknowledges the support of the Covino, S., Melandri, A., Salvaterra, R., et al. 2013b, MNRAS, 432, 1231, French National Research Agency (ANR) under contract ANR-16-CE31-0003 arXiv:1303.4743 BEaPro DM acknowledges support from the Instrument Center for Danish As- Coward, D. M., Howell, E. J., Branchesi, M., et al. 2013, MNRAS, 432, 2141 trophysics (IDA). This work made use of data supplied by the UK Swift Science Cucchiara, A., Fumagalli, M., Rafelski, M., et al. 2015, ApJ, 804, 51, Data Centre at the University of Leicester. Finally, it is with pleasure that we ac- arXiv:1408.3578 knowledge expert support from the ESO staﬀ at the Paranal and La Silla obser- D’Avanzo, P., Sparre, M., Watson, D., et al. 2011a, GRB Coordinates Network, vatories in obtaining these target of opportunity data. This research made use of Circular Service, No. 12284, #1 (2011), 12284 Astropy, a community-developed core Python package for Astronomy (The As- D’Avanzo, P., Vergani, S. D., Flores, H., et al. 2011b, GRB Coordinates Network, tropy Collaboration et al. 2013). The analysis and plotting was achieved using the Circular Service, No. 12542, #1 (2011), 12542 Python-based packages Matplotlib (Hunter 2007), Numpy, and Scipy (Oliphant De Pasquale, M., Piro, L., Perna, R., et al. 2003, ApJ, 592, 1018 2007; Van Der Walt et al. 2011), along with other community-developed pack- de Ugarte Postigo, A., Blazek, M., Janout, P., et al. 2014a, in Proc. SPIE, Vol. ages. 9152, id. 91520B 9 pp. (2014)., ed. G. Chiozzi & N. M. Radziwill, Vol. 9152, 91520B de Ugarte Postigo, A., Cano, Z., Izzo, L., et al. 2016a, GRB Coordinates Net- work, Circular Service, No. 20342, #1 (2016), 20342 References de Ugarte Postigo, A., D’Odorico, S., Vernet, J., et al. 2009a, GRB Coordinates Network, 9015 Abbott, B. P., Abbott, R., Abbott, T. D., et al. 2017, Phys. Rev. Lett., 119, 161101, de Ugarte Postigo, A., Fynbo, J. P. U., Jakobsson, P., et al. 2011a, GRB Coordi- arXiv:1710.05836 nates Network, Circular Service, No. 12258, #1 (2011), 12258 Abrarov, S. M. & Quine, B. M. 2015a, J. Math. Res., 7, 163 de Ugarte Postigo, A., Goldoni, P., Milvang-Jensen, B., et al. 2011b, GRB Coor- Abrarov, S. M. & Quine, B. M. 2015b, J. Math. Res., 7, 44, arXiv:1411.1024v2 dinates Network, Circular Service, No. 11579, #1 (2011), 11579 Ackermann, M., Ajello, M., Asano, K., et al. 2014, Sci, 343, 42, arXiv:1311.5623 de Ugarte Postigo, A., Goldoni, P., Thöne, C. C., et al. 2009b, GRB Coordinates Amati, L., Guidorzi, C., Frontera, F., et al. 2008, MNRAS, 391, 577, Network, Circular Service, No. 1042, #1 (2009), 10042 arXiv:0805.0377 de Ugarte Postigo, A., Goldoni, P., Thöne, C. C., et al. 2010, A&A, 513, A42, Antonelli, L. A., Maund, J. R., Palazzi, E., et al. 2010, GRB Coordinates Net- arXiv:0911.3901 work, Circular Service, No. 10620, #1 (2010), 10620 de Ugarte Postigo, A., Kann, D. A., Izzo, L., et al. 2016b, GRB Coordinates Arabsalmani, M., Møller, P., Perley, D. A., et al. 2018, MNRAS, 473, 3312 Network, Circular Service, No. 20014, #1 (2016), 20014 Ashall, C., Pian, E., Mazzali, P. A., et al. 2017, eprint arXiv:1702.04339, 1, de Ugarte Postigo, A., Krühler, T., Flores, H., & Fynbo, J. P. U. 2015a, GRB arXiv:1702.04339 Coordinates Network, Circular Service, No. 17523, #1 (2015), 17523 Axelsson, M., Baldini, L., Barbiellini, G., et al. 2012, Astrophys. J. Lett., 757, de Ugarte Postigo, A., Malesani, D., Krühler, T., et al. 2016c, GRB Coordinates L31, arXiv:1207.6109v1 Network, Circular Service, No. 19154, #1 (2016), 19154 Bañados, E., Venemans, B. P., Mazzucchelli, C., et al. 2017, Nature, 553, 473, de Ugarte Postigo, A., Pugliese, G., Xu, D., & Malesani, D. 2016d, GRB Coor- arXiv:1712.01860 dinates Network, Circular Service, No. 18886, #1 (2016), 18886 Band, D. L. 2006, ApJ, 644, 378, arXiv:0602267 de Ugarte Postigo, A., Thöne, C., Lombardi, G., & Perez, A. 2015b, GRB Coor- Barthelmy, S. D. 2000, in AIP Conf. Proc., Vol. 526 (AIP), 731–735 dinates Network, Circular Service, No. 18274, #1 (2015), 18274 Barthelmy, S. D., Barbier, L. M., Cummings, J. R., et al. 2005, SSRv, 120, 143, De Ugarte Postigo, A., Thöne, C. C., Goldoni, P., & Fynbo, J. P. U. 2011, AN, arXiv:0507410 332, 297 Baumgartner, W. H., Tueller, J., Markwardt, C. B., et al. 2013, Astrophys. Jour- de Ugarte Postigo, A., Thöne, C. C., Rowlinson, A., et al. 2014b, A&A, 563, nal, Suppl. Ser., 207, 19, arXiv:1212.3336 A62, arXiv:1308.2984 Benner, D. C., Rinsland, C. P., Devi, V. M., Smith, M. A. H., & Atkins, D. 1995, de Ugarte Postigo, A., Xu, D., Malesani, D., & Tanvir, N. R. 2015c, GRB Coor- JQSRT, 53, 705 dinates Network, Circular Service, No. 17822, #1 (2015), 17822 Berger, E. 2011, GRB Coordinates Network, Circular Service, No. 12193, #1 D’Elia, V., Fiore, F., Perna, R., et al. 2009, ApJ, 694, 332, arXiv:0804.2141 (2011), 12193 D’Elia, V., Fynbo, J. P. U., Covino, S., et al. 2010, A&A, 523, A36 Berger, E., Chornock, R., Holmes, T. R., et al. 2011, ApJ, 743, arXiv:1106.3073 D’Elia, V., Fynbo, J. P. U., Goldoni, P., et al. 2014, A&A, 564, A38 Berger, E., Fong, W., & Chornock, R. 2013, Astrophys. J. Lett., 774, L23, D’Elia, V., Goldoni, P., Malesani, D., et al. 2012a, GRB Coordinates Network, arXiv:1306.3960 Circular Service, No. 13227, #1 (2012), 13227 Berger, E., Kulkarni, S., & Frail, D. 2003, ApJ, 590, 379 D’Elia, V., Goldoni, P., Xu, D., et al. 2012b, GRB Coordinates Network, Circular Bloom, J. S., Kulkarni, S. R., & Djorgovski, S. G. 2002, AJ, 123, 1111 Service, No. 13494, #1 (2012), 13494 Bouwens, R. J., Illingworth, G. D., Oesch, P. A., et al. 2015, ApJ, 803, 1, D’Elia, V., Krühler, T., Wiersema, K., et al. 2015a, GRB Coordinates Network, arXiv:1403.4295 Circular Service, No. 18187, #1 (2015), 18187 Briggs, M. S., Paciesas, W. S., Pendleton, G. N., et al. 1996, ApJ, 459, 40, D’Elia, V., Krühler, T., Xu, D., et al. 2015b, GRB Coordinates Network, Circular arXiv:9509078 Service, No. 18318, #1 (2015), 18318 Bufano, F., Bufano, F., Benetti, S., et al. 2010, GRB Coordinates Network, Cir- D’Elia, V. & Stratta, G. 2011, A&A, 532, A48, arXiv:1108.1077 cular Service, No. 10543, #1 (2010), 10543 D’Elia, V., Xu, D., de Ugarte Postigo, A., et al. 2012c, GRB Coordinates Net- Bufano, F., Pian, E., Sollerman, J., et al. 2012, ApJ, 753, 67, arXiv:1111.4527 work, Circular Service, No. 13507, #1 (2012), 13507 Burrows, D. N., Hill, J. E., Nousek, J. A., et al. 2005, SSRv, 120, 165, D’Elia, V.,Xu, D., Malesani, D., et al. 2016, GRB Coordinates Network, Circular arXiv:0508071 Service, No. 19192, #1 (2016), 19192 Burrows, D. N., Kennea, J. A., Abbey, A. F., et al. 2007in (International Society Dessauges-Zavadsky, M., Chen, H.-W., Prochaska, J. X., Bloom, J. S., & Barth, for Optics and Photonics), 668607 A. J. 2006, ApJ, 648, L89 Campana, S., Salvaterra, R., Melandri, A., et al. 2012, MNRAS, 421, 1697, Djorgovski, S. G., Frail, D. A., Kulkarni, S. R., et al. 2001, ApJ, 562, 654, arXiv:1112.5111 arXiv:0107539 Cano, Z., Bersier, D., Guidorzi, C., et al. 2011, ApJ, 740, arXiv:1104.5141 Elliott, J., Schmidl, S., Greiner, J., et al. 2013, GRB Coordinates Network, Cir- Cano, Z., Izzo, L., de Ugarte Postigo, A., et al. 2017, A&A, 605, A107, cular Service, No. 14898, #1 (2013), 14898 arXiv:1704.05401 Evans, P. A., Beardmore, A. P., Page, K. L., et al. 2009, MNRAS, 397, 1177, Cardelli, J. a., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 345, 245, arXiv:0812.3662 arXiv:0706.1275 Evans, P. A., Beardmore, A. P., Page, K. L., et al. 2007, A&A, 469, 379, Castro-Tirado, A. J., Møller, P., García-Segura, G., et al. 2010, A&A, 517, A61, arXiv:0704.0128 arXiv:1004.1131v1 Filippenko, A. V. 1982, PASP, 94, 715 Chen, T.-W., Klose, S., Guelbenzu, A. N., et al. 2016, GRB Coordinates Net- Fiore, F., D’Elia, V., Lazzati, D., et al. 2005, ApJ, 624, 853, arXiv:0409717 work, Circular Service, No. 19975, #1 (2016), 19975 Flores, H., Covino, S., de Ugarte Postigo, A., et al. 2013a, GRB Coordinates Chornock, R. & Berger, E. 2011, GRB Coordinates Network, Circular Service, Network, Circular Service, No. 14493, #1 (2013), 14493 No. 11518, #1 (2011), 11518 Flores, H., Covino, S., Xu, D., et al. 2013b, GRB Coordinates Network, Circular Christensen, L., Fynbo, J. P. U., Prochaska, J. X., et al. 2011, ApJ, 727, 73 Service, No. 14491, #1 (2013), 14491

Article number, page 16 of 43 Selsing et al.: The X-shooter GRB afterglow legacy sample

Flores, H., Fynbo, J. P. U., de Ugarte Postigo, A., et al. 2010, GRB Coordinates Jakobsson, P., Hjorth, J., Fynbo, J. P. U., et al. 2004, ApJ, 617, L21, Network, Circular Service, No. 11317, #1 (2010), 11317 arXiv:0411036 Fong, W., Berger, E., Chornock, R., et al. 2013, ApJ, 769, 56, arXiv:1302.3221 Jakobsson, P., Hjorth, J., Malesani, D., et al. 2012, ApJ, 752, 62 Freedman, D. L. & Waxman, E. 2001, ApJ, 547, 922, arXiv:9912214v1 Jakobsson, P., Levan, A., Fynbo, J. P. U., et al. 2006, A&A, 447, 897, Freudling, W., Romaniello, M., Bramich, D. M., et al. 2013, A&A, 559, A96, arXiv:0509888 arXiv:1311.5411v1 Japelj, J., Covino, S., Gomboc, A., et al. 2015, A&A, 579, A74 Friis, M., De Cia, A., Kruhler, T., et al. 2015, MNRAS, 451, 167, Jones, A., Noll, S., Kausch, W., Szyszka, C., & Kimeswenger, S. 2013, A&A, arXiv:1409.6315v2 560, A91, arXiv:1310.7030 Fruchter, A. S., Levan, A. J., Strolger, L., et al. 2006, Nature, 441, 463, Kann, D. A., Schady, P., E., F. O., et al. 2017, eprint arXiv:1706.00601, arXiv:0603537 arXiv:1706.00601 Fynbo, J. P. U., de Ugarte Postigo, A., D’Elia, V., Xu, D., & Malesani, D. 2012a, Kausch, W., Noll, S., Smette, A., et al. 2015, A&A, 576, A78, GRB Coordinates Network, Circular Service, No. 13477, #1 (2012), 13477 arXiv:1501.07239v1 Fynbo, J. P. U., Jakobsson, P., Prochaska, J. X., et al. 2009, ApJS, 185, 175, Knust, F., Kann, D. A., Krühler, T., & Greiner, J. 2015, GRB Coordinates Net- arXiv:0907.3449 work, Circular Service, No. 17767, #1 (2015), 17767 Fynbo, J. P. U., Krühler, T., Leighly, K., et al. 2014, A&A, 12, 1, Krühler, T., Fynbo, J. P. U., Milvang-Jensen, B., Tanvir, N., & Jakobsson, P. arXiv:1409.4975 2012, GRB Coordinates Network, Circular Service, No. 13134, #1 (2012), Fynbo, J. P. U., Tanvir, N. R., D’Elia, V., et al. 2012b, GRB Coordinates Net- 13134 work, Circular Service, No. 14120, #1 (2012), 14120 Krühler, T., Greiner, J., Schady, P., et al. 2011, A&A, 534, A108, Fynbo, J. P. U., Tanvir, N. R., Jakobsson, P., et al. 2014, GRB Coordinates Net- arXiv:1108.0674 work, Circular Service, No. 16217, #1 (2014), 16217 Krühler, T., Ledoux, C., Fynbo, J. P. U., et al. 2013, A&A, 557, A18, Fynbo, J. P. U., Xu, D., Malesani, D., et al. 2013, GRB Coordinates Network, arXiv:1304.7003 Circular Service, No. 14286, #1 (2013), 14286 Krühler, T., Malesani, D., de Ugarte Postigo, A., Melandri, A., & Fynbo, J. P. U. Fynbo, J. U., Jensen, B. L., Gorosabel, J., et al. 2001, A&A, 369, 373 2014a, GRB Coordinates Network, Circular Service, No. 16194, #1 (2014), Gehrels, N., Ramirez-Ruiz, E., & Fox, D. 2009, ARA&A, 47, 567, 16194 arXiv:0909.1531 Krühler, T., Malesani, D., Fynbo, J. P. U., et al. 2015, A&A, 581, A125, Gehrels, N., White, N., Barthelmy, S., et al. 2004, ApJ, 611, 1005, arXiv:1505.06743 arXiv:0405233v1 Krühler, T., Malesani, D., Milvang-Jensen, B., et al. 2012, ApJ, 758, 46, Ghirlanda, G. 2007, R. Soc. London Philos. Trans. Ser. A, 365, 1385 arXiv:1205.4036 Ginsburg, A., Parikh, M., Woillez, J., et al. 2016, Astroquery V0.3.1 Krühler, T., Malesani, D., Xu, D., et al. 2016a, GRB Coordinates Network, Cir- Goldoni, P., de Ugarte Postigo, A., & Fynbo, J. P. U. 2013, GRB Coordinates cular Service, No. 19186, #1 (2016), 19186 Network, Circular Service, No. 15571, #1 (2013), 15571 Krühler, T., Malesani, D., Xu, D., et al. 2013a, GRB Coordinates Network, Cir- Goldoni, P., Flores, H., Malesani, D., et al. 2010, GRB Coordinates Network, cular Service, No. 14264, #1 (2013), 14264 Circular Service, No. 10684, #1 (2010), 10684 Krühler, T., Tanvir, N. R., de Ugarte Postigo, A., et al. 2012, GRB Coordinates Goldoni, P., Royer, F., François, P., et al. 2006, Ground-based Airborne Instrum. Network, Circular Service, No. 13930, #1 (2012), 13930 Astron. Ed. by McLean, 6269, 80 Krühler, T., Tanvir, N. R., Malesani, D., Xu, D., & Fynbo, J. P. U. 2014b, GRB Graff, P. B., Lien, A. Y., Baker, J. G., & Sakamoto, T. 2016, ApJ, 818, 55, Coordinates Network, Circular Service, No. 15900, #1 (2014), 15900 Krühler, T., Vreeswijk, P. M., & Fynbo, J. P. U. 2014c, GRB Coordinates Net- arXiv:1509.01228 work, Circular Service, No. 16401, #1 (2014), 16401 Greiner, J., Bolmer, J., Wieringa, M., et al. 2018, eprint arXiv:1802.01882, Krühler, T., Xu, D., Bolmer, J., et al. 2016b, GRB Coordinates Network, Circular arXiv:1802.01882 Service, No. 19971, #1 (2016), 19971 Greiner, J., Bornemann, W., Clemens, C., et al. 2008, PASP, 120, 405, Krühler, T., Xu, D., Fynbo, J. P. U., et al. 2015, GRB Coordinates Network, arXiv:0801.4801 Circular Service, No. 17420, #1 (2015), 17420 Greiner, J., Krühler, T., Klose, S., et al. 2011, A&A, 526, A30 Krühler, T., Xu, D., Sanchez-Ramirez, R., et al. 2013b, GRB Coordinates Net- Greiner, J., Mazzali, P. A., Kann, D. A., et al. 2015, Nature, 523, 189, work, Circular Service, No. 14390, #1 (2013), 14390 arXiv:1509.03279# Letchworth, K. L. & Benner, D. C. 2007, JQSRT, 107, 173 Greiner, J., Yu, H.-F., Krühler, T., et al. 2014, A&A, 568, A75, arXiv:1407.4210 Levan, A. J., Tanvir, N. R., Starling, R. L. C., et al. 2013, ApJ, 781, 13, Groot, P., Kaper, L., Ellerbroek, L., et al. 2010, GRB Coordinates Network, Cir- arXiv:1302.2352 cular Service, No. 10441, #1 (2010), 10441 Levesque, E. M., Kewley, L. J., Graham, J. F., & Fruchter, A. S. 2010, ApJ, 712, Groot, P. J., Galama, T. J., van Paradijs, J., et al. 1998, ApJ, 493, L27 L26 Hamuy, M., Suntzeff, N. B., Heathcote, S. R., et al. 1994, PASP, 106, 566 Lien, A., Sakamoto, T., Barthelmy, S. D., et al. 2016, ApJ, 829, 7, Hartoog, O. E., Malesani, D., Fynbo, J. P. U., et al. 2015, A&A, 580, A139, arXiv:1606.01956 arXiv:1409.4804 Lien, A., Sakamoto, T., Gehrels, N., et al. 2014, ApJ, 783, 24, arXiv:1311.4567 Hartoog, O. E., Malesani, D., Sanchez-Ramirez, R., et al. 2014, GRB Coordi- LIGO Scientific Collaboration, Virgo Collaboration, GBM, F., et al. 2017, As- nates Network, Circular Service, No. 16437, #1 (2014), 16437 trophys. J. Lett., 848, L12, arXiv:1710.05833 Hartoog, O. E., Malesani, D., Wiersema, K., et al. 2012, GRB Coordinates Net- Littlejohns, O. M., Tanvir, N. R., Willingale, R., et al. 2013, MNRAS, 436, 3640, work, Circular Service, No. 13730, #1 (2012), 13730 arXiv:1309.7045 Hartoog, O. E., Wiersema, K., Vreeswijk, P. M., et al. 2013, MNRAS, 430, 2739, López, S., D’Odorico, V., Ellison, S. L., et al. 2016, A&A, 594, A91, arXiv:1301.3912 arXiv:1607.08776 Hartoog, O. E., Xu, D., Malesani, D., et al. 2013, GRB Coordinates Network, Lyman, J. D., Levan, A. J., Tanvir, N. R., et al. 2017, MNRAS, 1817, stx220 Circular Service, No. 15494, #1 (2013), 15494 Malesani, D., Fynbo, J. P. U., Christensen, L., et al. 2009a, GRB Coordinates Heintz, K. E., Fynbo, J. P. U., Jakobsson, P., et al. 2017, A&A, 601, A83, Network, 9761 arXiv:1703.07109 Malesani, D., Goldoni, P., Fynbo, J. P. U., et al. 2009b, GRB Coordinates Net- Heintz, K. E., Malesani, D., de Ugarte Postigo, A., et al. 2016, GRB Coordinates work, 9942 Network, Circular Service, No. 20020, #1 (2016), 20020 Malesani, D., Krühler, T., Heintz, K. E., & Fynbo, J. P. U. 2016, GRB Coordi- Heintz, K. E., Malesani, D., Wiersema, K., et al. 2018, MNRAS, 474, 2738, nates Network, Circular Service, No. 20180, #1 (2016), 20180 arXiv:1711.02706v1 Malesani, D., Krühler, T., Perley, D., et al. 2013a, GRB Coordinates Network, Hjorth, J. & Bloom, J. S. 2011, in Gamma-Ray Bursts, ed. C. Kouveliotou, R. A. Circular Service, No. 14225, #1 (2013), 14225 M. J. Wijers, & S. Woosley No. 1973 (Cambridge: Cambridge University Malesani, D., Krühler, T., Xu, D., et al. 2015a, GRB Coordinates Network, Cir- Press), 169–190 cular Service, No. 17755, #1 (2015), 17755 Hjorth, J., Malesani, D., Jakobsson, P., et al. 2012, ApJ, 756, 187, Malesani, D. & Palazzi, E. 2010, GRB Coordinates Network, Circular Service, arXiv:1205.3162 No. 10631, #1 (2010), 10631 Hjorth, J., Melandri, A., Malesani, D., Krühler, T., & Xu, D. 2013, GRB Coor- Malesani, D., Schulze, S., de Ugarte Postigo, A., et al. 2012, GRB Coordinates dinates Network, Circular Service, No. 14365, #1 (2013), 14365 Network, Circular Service, No. 13649, #1 (2012), 13649 Hogg, D. W. & Fruchter, A. S. 1999, ApJ, 520, 54 Malesani, D., Tanvir, N. R., Krühler, T., et al. 2015b, GRB Coordinates Network, Horne, K. 1986, PASP, 98, 609 Circular Service, No. 18540, #1 (2015), 18540 Hunter, J. D. 2007, Comput. Sci. Eng., 9, 99, arXiv:0402594v3 Malesani, D., Xu, D., Fynbo, J. P. U., et al. 2013b, GRB Coordinates Network, Izzo, L., Thöne, C. C., Schulze, S., et al. 2017, MNRAS, 472, 4480, Circular Service, No. 14291, #1 (2013), 14291 arXiv:1704.05509 Markwardt, C. B., Barthelmy, S. D., Cummings, J. R., et al. 2016, GRB Coordi- Jakobsson, P., Frail, D. A., Fox, D. B., et al. 2005, ApJ, 629, 45 nates Network, Circular Service, No. 19974, #1 (2016), 19974

Article number, page 17 of 43 A&A proofs: manuscript no. XSGRB_sample_arxiv

Maselli, A., Melandri, A., Nava, L., et al. 2014, Sci, 343, 48, arXiv:1311.5254 Rossi, A., Piranomonte, S., Savaglio, S., et al. 2014, A&A, 47, 1, Meegan, C. A., Fishman, G. J., Wilson, R. B., et al. 1992, Nature, 355, 143 arXiv:1409.0017 Melandri, A., Bernardini, M. G., D’Avanzo, P., et al. 2015, A&A, 581, A86, Rowlinson, A., Patruno, A., & O’Brien, P. T. 2017, MNRAS, 472, 1152, arXiv:1506.03079 arXiv:1706.08538 Melandri, A., Sbarufatti, B., D’Avanzo, P., et al. 2012, MNRAS, 421, 1265, Sakamoto, T., Barthelmy, S. D., Cummings, J. R., et al. 2016, GRB Coordinates arXiv:1112.4480 Network, Circular Service, No. 19276, #1 (2016), 19276 Michałowski, M. J., Xu, D., Stevens, J., et al. 2016, eprint arXiv:1610.06928, Salvaterra, R., Campana, S., Vergani, S. D., et al. 2012, ApJ, 749, 68, arXiv:1610.06928 arXiv:1112.1700 Micol, A., Arnaboldi, M., Delmotte, N. A. R., Mascetti, L., & Retzlaff, J. 2016, Salvaterra, R., Valle, M. D., Campana, S., et al. 2009, Nature, 461, 1258, in Obs. Oper. Strateg. Process. Syst. VI, ed. A. B. Peck, R. L. Seaman, & C. R. arXiv:0906.1578 Benn, Vol. 9910 (International Society for Optics and Photonics), 991030 Sanchez-Ramirez, R., de Ugarte Postigo, A., Krühler, T., et al. 2012, GRB Co- Milvang-Jensen, B., Goldoni, P., Tanvir, N. R., et al. 2010, GRB Coordinates ordinates Network, Circular Service, No. 14035, #1 (2012), 14035 Network, Circular Service, No. 10876, #1 (2010), 10876 Sánchez-Ramírez, R., Hancock, P. J., Jóhannesson, G., et al. 2017, MNRAS, Modigliani, A., Goldoni, P., Royer, F., et al. 2010, SPIE Astron. Telesc. + In- 464, 4624, arXiv:1610.01844 strum., 7737, 773728 Savaglio, S. 2006, NJPh, 8, 195, arXiv:0609489 Moffat, A. F. J. 1969, A&A, 3, 455, arXiv:1011.1669v3 Schady, P., Krühler, T., Greiner, J., et al. 2015, A&A, 579, A126, Molinari, E., Vergani, S. D., Malesani, D., et al. 2007, A&A, 469, L13, arXiv:1505.04415v1 arXiv:0612607 Schlafly, E. F. & Finkbeiner, D. P. 2011, ApJ, 737, 103, arXiv:1012.4804 Morgan, A. N., Perley, D. A., Cenko, S. B., et al. 2014, MNRAS, 440, 1810, Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525, arXiv:1305.1928 arXiv:9710327 Noll, S., Kausch, W., Barden, M., et al. 2012, A&A, 543, A92 Schulze, S., Chapman, R., Hjorth, J., et al. 2015, ApJ, 808, 73, arXiv:1503.04246 Noterdaeme, P., Petitjean, P., Carithers, W. C., et al. 2012, A&A, 547, L1, Schulze, S., Covino, S., Flores, H., et al. 2011, GRB Coordinates Network, Cir- arXiv:1210.1213 cular Service, No. 12770, #1 (2011), 12770 Nysewander, M., Fruchter, A. S., & Pe’er, A. 2009, ApJ, 701, 824, Schulze, S., Levan, A. J., Malesani, D., et al. 2012, GRB Coordinates Network, arXiv:0806.3607 Circular Service, No. 13257, #1 (2012), 13257 Oesch, P. A., Brammer, G., van Dokkum, P. G., et al. 2016, ApJ, 819, 129, Schulze, S., Malesani, D., Cucchiara, A., et al. 2014, A&A, 102, 1 arXiv:1603.00461 Schulze, S., Wiersema, K., Xu, D., & Fynbo, J. P. U. 2014, GRB Coordinates Oliphant, T. E. 2007, CSE, 9, 10 Network, Circular Service, No. 15831, #1 (2014), 15831 Olivares E., F., Greiner, J., Schady, P., et al. 2015, A&A, 577, A44, Selsing, J., Fynbo, J. P. U., Christensen, L., & Krogager, J.-K. 2015, A&A, 87, arXiv:1502.00630 14, arXiv:1510.08058 Olivares E., F., Greiner, J., Schady, P., et al. 2012, A&A, 539, A76, Selsing, J., Heintz, K. E., Malesani, D., et al. 2016a, GRB Coordinates Network, arXiv:1110.4109 Circular Service, No. 20061, #1 (2016), 20061 Osborne, J. P., Beardmore, A. P., Evans, P. A., & Goad, M. R. 2016, GRB Coor- Selsing, J., Krühler, T., Malesani, D., et al. 2018, eprint arXiv:1707.01452, 1, dinates Network, Circular Service, No. 19969, #1 (2016), 19969 arXiv:1707.01452 Pagnini, G. & Mainardi, F. 2010, JCoAM, 233, 1590 Selsing, J., Vreeswijk, P. M., Japelj, J., et al. 2016b, GRB Coordinates Network, Palmerio, J., Krühler, T., Malesani, D., & Fynbo, J. P. U. 2017, GRB Coordinates Circular Service, No. 19274, #1 (2016), 19274 Network, Circular Service, No. 20589, #1 (2017), 20589 Sheffer, Y., Prochaska, J. X., Draine, B. T., Perley, D. A., & Bloom, J. S. 2009, Perley, D. A., Bloom, J. S., Cenko, S. B., et al. 2010, GRB Coordinates Network, ApJ, 701, L63, arXiv:0907.1285 Circular Service, No. 10399, #1 (2010), 10399 Smette, A., Sana, H., Noll, S., et al. 2015, A&A, 576, A77 Perley, D. A., Cenko, S. B., Bloom, J. S., et al. 2009, AJ, 138, 1690, Sparre, M., de Ugarte Postigo, A., Fynbo, J. P. U., et al. 2011, GRB Coordinates arXiv:0905.0001 Network, Circular Service, No. 11607, #1 (2011), 11607 Perley, D. A., Cenko, S. B., Corsi, A., et al. 2014, ApJ, 781, 37, arXiv:1307.4401 Sparre, M., Hartoog, O. E., Krühler, T., et al. 2014, ApJ, 785, 150, Perley, D. A., Krühler, T., Schady, P., et al. 2017, Mon. Not. R. Astron. Soc. Lett., arXiv:1309.2940 465, L89 Sparre, M., Sollerman, J., Fynbo, J. P. U., et al. 2011, ApJ, 735, L24 Perley, D. A., Krühler, T., Schulze, S., et al. 2016a, ApJ, 817, 7, Starling, R. L. C., Wiersema, K., Levan, A. J., et al. 2011, MNRAS, 411, 2792, arXiv:1504.02482 arXiv:1004.2919 Perley, D. A., Levan, A. J., Tanvir, N. R., et al. 2013, ApJ, 778, 128, Sudilovsky, V., Kann, D. A., Schady, P., et al. 2013, GRB Coordinates Network, arXiv:1301.5903 Circular Service, No. 15250, #1 (2013), 15250 Perley, D. A., Morgan, A. N., Updike, A., et al. 2011, AJ, 141, 36, Tanvir, N. R., Fox, D. B., Levan, A. J., et al. 2009, Nature, 461, 1254, arXiv:1009.0004 arXiv:0906.1577 Perley, D. A., Perley, R. A., Hjorth, J., et al. 2015, ApJ, 801, 102, Tanvir, N. R., Fynbo, J. P. U., Melandri, A., et al. 2012, GRB Coordinates Net- arXiv:1407.4456 work, Circular Service, No. 13890, #1 (2012), 13890 Perley, D. A., Tanvir, N. R., Hjorth, J., et al. 2016b, ApJ, 817, 8, Tanvir, N. R., Krühler, T., De Cia, A., et al. 2016a, GRB Coordinates Network, arXiv:1504.02479 Circular Service, No. 20104, #1 (2016), 20104 Pescalli, A., Ghirlanda, G., Salvaterra, R., et al. 2015, A&A, 587, A40, Tanvir, N. R., Krühler, T., Malesani, D., et al. 2015a, GRB Coordinates Network, arXiv:1506.05463 Circular Service, No. 18524, #1 (2015), 18524 Piran, T. 2005, RvMP, 76, 1143, arXiv:0405503 Tanvir, N. R., Krühler, T., Wiersema, K., et al. 2016b, GRB Coordinates Net- Piranomonte, S., Vergani, S. D., Malesani, D., et al. 2011, GRB Coordinates work, Circular Service, No. 20321, #1 (2016), 20321 Network, Circular Service, No. 12164, #1 (2011), 12164 Tanvir, N. R., Laskar, T., Levan, A. J., et al. 2017, eprint arXiv:1703.09052, Planck Collaboration, Ade, P. A. R., Aghanim, N., et al. 2016, A&A, 594, A13 arXiv:1703.09052 Pozanenko, A., Mazaeva, E., Sergeev, A., et al. 2015, GRB Coordinates Net- Tanvir, N. R., Levan, A. J., Fruchter, A. S., et al. 2013, Nature, 500, 547, work, Circular Service, No. 17903, #1 (2015), 17903 arXiv:1306.4971 Prochaska, J. X., Chen, H.-W., Dessauges-Zavadsky, M., & Bloom, J. S. 2007, Tanvir, N. R., Vergani, S., Hjorth, J., et al. 2010, GRB Coordinates Network, ApJ, 666, 267, arXiv:0703665 Circular Service, No. 11123, #1 (2010), 11123 Prochaska, J. X., Sheffer, Y., Perley, D. A., et al. 2009, ApJ, 691, L27, Tanvir, N. R., Wiersema, K., Xu, D., & Fynbo, J. P. U. 2013, GRB Coordinates arXiv:0901.0556 Network, Circular Service, No. 14882, #1 (2013), 14882 Pugliese, G., Covino, S., Krühler, T., Xu, D., & Tanvir, N. R. 2016, GRB Coor- Tanvir, N. R., Xu, D., Krühler, T., et al. 2016c, GRB Coordinates Network, Cir- dinates Network, Circular Service, No. 18982, #1 (2016), 18982 cular Service, No. 19350, #1 (2016), 19350 Pugliese, V., Xu, D., Tanvir, N. R., et al. 2015, GRB Coordinates Network, Cir- Tanvir, N. R., Xu, D., Zafar, T., Covino, S., & Schulze, S. 2015b, GRB Coordi- cular Service, No. 17672, #1 (2015), 17672 nates Network, Circular Service, No. 18080, #1 (2015), 18080 Rau, A., Krühler, T., & Greiner, J. 2013, GRB Coordinates Network, Circular Tepper García, T. 2006, MNRAS, 369, 2025, arXiv:0602124v2 Service, No. 15330, #1 (2013), 15330 The Astropy Collaboration, Robitaille, T. P., Tollerud, E. J., et al. 2013, A&A, Ricker, G. R. & Team, H. S. 2004, Am. Astron. Soc. HEAD Meet. #8, id.18.06; 558, A33, arXiv:1307.6212 Bull. Am. Astron. Soc. Vol. 36, p.942, 36, 942 Thöne, C. C., de Ugarte Postigo, A., Vreeswijk, P., et al. 2010, GRB Coordinates Rol, E., Wijers, R. A. M. J., Kouveliotou, C., Kaper, L., & Kaneko, Y. 2005, ApJ, Network, Circular Service, No. 10971, #1 (2010), 10971 624, 868, arXiv:0501375 Thöne, C. C., Fynbo, J. P. U., Goldoni, P., et al. 2013, MNRAS, 428, 3590, Roming, P. W., Kennedy, T. E., Mason, K. O., et al. 2005, SSRv, 120, 95, arXiv:1206.2337 arXiv:0507413 Thöne, C. C., Goldoni, P., Covino, S., et al. 2009, GRB Coordinates Network, Rossi, A., Klose, S., Ferrero, P., et al. 2012, A&A, 545, A77, arXiv:1202.1434 Circular Service, No. 10233, #1 (2009), 10233

Article number, page 18 of 43 Selsing et al.: The X-shooter GRB afterglow legacy sample

Thöne, C. C., Greiner, J., Savaglio, S., & Jehin, E. 2007, ApJ, 671, 628, Xu, D., Tanvir, N. R., Malesani, D., & Fynbo, J. P. U. 2015b, GRB Coordinates arXiv:0611772 Network, Circular Service, No. 18506, #1 (2015), 18506 Troja, E., Lipunov, V. M., Mundell, C. G., et al. 2017, Nature, 547, 425 Trujillo, I., Aguerri, J. A. L., Cepa, J., & Gutiérrez, C. M. 2001, MNRAS, 328, 977, arXiv:0109067v1 Ukwatta, T. N., Barthelmy, S. D., Beardmore, A. P., et al. 2016, GRB Coordi- nates Network, Circular Service, No. 19148, #1 (2016), 19148 Ukwatta, T. N. & Wózniak, P. R. 2016, MNRAS, 455, 703, arXiv:1507.07117 van der Horst, A. J., Kouveliotou, C., Gehrels, N., et al. 2009, ApJ, 699, 1087, arXiv:0905.0524 Van Der Walt, S., Colbert, S. C., & Varoquaux, G. 2011, Comput. Sci. Eng., 13, 22, arXiv:1102.1523 Xu, D., Vreeswijk, P. M., Fynbo, J. P. U., et al. 2014c, GRB Coordinates Net- van Dokkum, P. G. 2001, PASP, 113, 1420 work, Circular Service, No. 17040, #1 (2014), 17040 van Paradijs, J., Kouveliotou, C., & Wijers, R. A. M. J. 2000, ARA&A, 38, 379 Vergani, S. D., D’Avanzo, P., Levan, A. J., et al. 2010a, GRB Coordinates Net- work, Circular Service, No. 10512, #1 (2010), 10512 Vergani, S. D., D’Avanzo, P., Malesani, D., et al. 2010b, GRB Coordinates Net- work, Circular Service, No. 10495, #1 (2010), 10495 Vergani, S. D., Flores, H., Covino, S., et al. 2011, A&A, 535, A127, arXiv:1107.3841 Vergani, S. D., Levan, A. J., D’Avanzo, P., et al. 2010c, GRB Coordinates Net- work, Circular Service, No. 10513, #1 (2010), 10513 Vergani, S. D., Petitjean, P., Ledoux, C., et al. 2009, A&A, 503, 771, Zheng, W., Filippenko, A. V., Yuk, H., Zhu, Y., & Perley, D. A. 2015, GRB arXiv:0906.3269 Coordinates Network, Circular Service, No. 18273, #1 (2015), 18273 Vergani, S. D., Piranomonte, S., Hartoog, O. E., et al. 2011, GRB Coordinates Network, Circular Service, No. 12677, #1 (2011), 12677 Vergani, S. D., Salvaterra, R., Japelj, J., et al. 2015, A&A, 581, A102, arXiv:1409.7064 Vernet, J., Dekker, H., D’Odorico, S., et al. 2011, A&A, 536, A105, arXiv:1110.1944 Vernet, J., Kerber, F., Mainieri, V., et al. 2009, Proc. Int. Astron. Union, 5, 535 Vreeswijk, P., Fynbo, J., & Melandri, A. 2011, GRB Coordinates Network, Cir- cular Service, No. 12648, #1 (2011), 12648 Vreeswijk, P. M., Kaufer, A., Spyromilio, J., et al. 2010, The VLT Rapid- Zitrin, A., Labbé, I., Belli, S., et al. 2015, Astrophys. J. Lett., 810, L12, Response Mode: implementation and scientiﬁc results arXiv:1507.02679 Vreeswijk, P. M., Ledoux, C., Raassen, A. J. J., et al. 2013, A&A, 549, A22, arXiv:1209.1506 Vreeswijk, P. M., Ledoux, C., Smette, A., et al. 2007, A&A, 468, 83, arXiv:0611690 Watson, D. & Jakobsson, P. 2012, ApJ, 754, 89, arXiv:1206.0005v1 Watson, D., Zafar, T., Andersen, A. C., et al. 2013, ApJ, 768, 23, arXiv:1212.4492 Wiersema, K., Curran, P. A., Krühler, T., et al. 2012, MNRAS, 426, 2, arXiv:1203.4265 Wiersema, K., Flores, H., D’Elia, V., et al. 2011, GRB Coordinates Network, Circular Service, No. 12431, #1 (2011), 12431 Wiersema, K., Goldoni, P., Fynbo, J. P. U., et al. 2012, GRB Coordinates Net- work, Circular Service, No. 12991, #1 (2012), 12991 Wijers, R. A. M. J., Bloom, J. S., Bagla, J. S., & Natarajan, P. 1998, MNRAS, 294, L13, arXiv:9708183 Xu, D., De Ugarte Postigo, A., Leloudas, G., et al. 2013, ApJ, 776, arXiv:1305.6832v2 Xu, D., de Ugarte Postigo, A., Malesani, D., et al. 2013a, GRB Coordinates Network, Circular Service, No. 14757, #1 (2013), 14757 Xu, D., de Ugarte Postigo, A., Malesani, D., et al. 2013b, GRB Coordinates Network, Circular Service, No. 14956, #1 (2013), 14956 1 Xu, D., Fynbo, J. P. U., D’Elia, V., & Tanvir, N. R. 2012, GRB Coordinates Dark Cosmology Centre, Niels Bohr Institute, University of Copen- Network, Circular Service, No. 13460, #1 (2012), 13460 hagen, Juliane Maries Vej 30, 2100 København Ø, Denmark Xu, D., Heintz, K. E., Malesani, D., & Fynbo, J. P. U. 2017, GRB Coordinates Network, Circular Service, No. 20458, #1 (2017), 20458 Xu, D., Heintz, K. E., Malesani, D., Wiersema, K., & Fynbo, J. P. U. 2016a, GRB Coordinates Network, Circular Service, No. 19773, #1 (2016), 19773 Xu, D., Levan, A. J., de Ugarte Postigo, A., et al. 2015a, GRB Coordinates Net- 2 The Cosmic Dawn Center, Niels Bohr Institute, University of work, Circular Service, No. 17832, #1 (2015), 17832 Copenhagen, Juliane Maries Vej 30, 2100 København Ø, Denmark Xu, D., Levan, A. J., Fynbo, J. P. U., et al. 2014a, GRB Coordinates Network, Circular Service, No. 16983, #1 (2014), 16983 Xu, D., Malesani, D., Fynbo, J. P. U., et al. 2016b, GRB Coordinates Network, Circular Service, No. 19600, #1 (2016), 19600 Xu, D., Malesani, D., Krühler, T., et al. 2013c, GRB Coordinates Network, Cir- cular Service, No. 14273, #1 (2013), 14273 3 DTU Space, National Space Institute, Technical University of Den- Xu, D., Malesani, D., Schulze, S., et al. 2013d, GRB Coordinates Network, Cir- mark, Elektrovej 327, 2800 Lyngby, Denmark cular Service, No. 15451, #1 (2013), 15451 Xu, D., Malesani, D., Schulze, S., et al. 2013e, GRB Coordinates Network, Cir- cular Service, No. 14816, #1 (2013), 14816 Xu, D., Malesani, D., Tanvir, N., Krühler, T., & Fynbo, J. 2013f, GRB Coordi- nates Network, Circular Service, No. 15450, #1 (2013), 15450 Xu, D., Malesani, D., Tanvir, N. R., et al. 2014b, GRB Coordinates Network, 4 APC, Astroparticule et Cosmologie, Université Paris Diderot, Circular Service, No. 15645, #1 (2014), 15645 CNRS/IN2P3, CEA/Irfu, Observatoire de Paris, Sorbonne Paris

Article number, page 19 of 43 A&A proofs: manuscript no. XSGRB_sample_arxiv

Cité, 10, Rue Alice Domon et Léonie Duquet, 75205, Paris Cedex 13, France 5 Max-Planck-Institut für extraterrestrische Physik, Giessenbach- straße, 85748 Garching, Germany 6 INAF - Osservatorio Astronomico di Roma, Via Frascati 33, 00078, Monte Porzio Catone (Roma), Italy 7 IRFU, CEA, Université Paris-Saclay, F-91191 Gif-sur-Yvette, France 8 Université Paris Diderot, AIM, Sorbonne Paris Cité, CEA, CNRS, F-91191 Gif-sur-Yvette, France 9 European Southern Observatory, Alonso de Córdova 3107, Vitacura, Casilla 19001, Santiago 19, Chile 10 Instituto de Astrofísica de Andalucía (IAA-CSIC), Glorieta de la Astronomía s/n, E-18008, Granada, Spain. 11 INAF - Osservatorio Astronomico di Brera, via Bianchi 46, 23807, Merate (LC), Italy 12 Space Science Data Center - Agenzia Spaziale Italiana, via del Po- litecnico, s.n.c., I-00133, Roma, Italy 13 European Southern Observatory, Karl-Schwarzschild Str. 2, 85748 Garching bei München, Germany 14 GEPI, Observatoire de Paris, PSL University, CNRS, 5 place Jules Janssen 92195, Meudon, France 15 KTH Royal Institute of Technology, Department of Physics, 106 91 Stockholm, Sweden 16 The Oskar Klein Centre for Cosmoparticle Physics, AlbaNova Uni- versity Centre, 106 91 Stockholm, Sweden 17 Center for Astrophysics and Cosmology, University of Nova Gorica, Vipavska 11c, 5270 Ajdovšcina,ˇ Slovenia 18 Department of Astrophysics, IMAPP, Radboud University Ni- jmegen, PO Box 9010, 6500 GL Nijmegen, the Netherlands 19 Anton Pannekoek Institute for Astronomy, University of Amster- dam, Science Park 904, 1098 XH Amsterdam, the Netherlands 20 Centre for Astrophysics and Cosmology, Science Institute, Univer- sity of Iceland, Dunhagi 5, 107, Reykjavik, Iceland 21 Department of Physics, University of Warwick, Coventry, CV4 7AL, UK 22 IASF/INAF Bologna, via Piero Gobetti 101, 40129 Bologna, Italy 23 Sorbonne Université, CNRS, UMR7095, Institut d’Astrophysique de Paris, F-75014, Paris, France 24 Astrophysics Research Institute, Liverpool John Moores University, IC2, Liverpool Science Park, 146 Brownlow Hill, Liverpool L3 5RF, UK 25 INAF, Istituto Astroﬁsica e Planetologia Spaziali, Via Fosso del Cavaliere 100, I-00133 Roma, Italy 26 Physics Department, University of Calabria, I-87036 Arcavacata di Rende, Italy 27 Department of Particle Physics and Astrophysics, Weizmann Insti- tute of Science, 234 Herzl Street, Rehovot, 761000, Israel 28 Department of Astronomy, Stockholm University, AlbaNova, 10691 Stockholm, Sweden 29 Institut für Physik und Astronomie, Universität Potsdam, Karl- Liebknecht-Str. 24/25, 14476 Golm, Germany 30 Department of Physics and Astronomy, University of Leicester, Uni- versity Road, Leicester, LE1 7RH, United Kingdom 31 CAS Key Laboratory of Space Astronomy and Technology, National Astronomical Observatories, Chinese Academy of Sciences, Bei- jing, 100012, P.R. China 32 Australian Astronomical Observatory, PO Box 915, North Ryde, NSW 1670, Australia.

Article number, page 20 of 43 Selsing et al.: The X-shooter GRB afterglow legacy sample

Table 1. The full sample of afterglows and hosts observed in the program. We here list the burst names and details of the spectroscopic observations. The exposure times and slit widths are given in the order UVB/VIS/NIR. The column ∆t shows the time after trigger when the spectroscopic observation was started. Magacq gives the approximate magnitude (typically in the R-band) of the afterglow or the host in the acquisition image.

GRB Obs Date Exptime Slit width Airmass Seeing ∆t Magacq Redshift Notes (ks) (arcsec) (arcsec) (hr) GRB090313a 2009-03-15 6.9/6.9/6.9 1.0/0.9/0.9 1.2–1.4 1.5 45 21.6 3.374 B.1 GRB090530a 2009-05-30 4.8/4.8/4.8 1.0/1.2/1.2 1.6–2.2 1.7 20.6 22 1.266 B.2 GRB090809a 2009-08-10 7.2/7.2/7.2 1.0/0.9/0.9 1.2–1.1 1.1 10.2 21 2.737 B.3 GRB090926Aa 2009-09-27 7.2/7.2/7.2 1.0/0.9/0.9 1.4–1.5 0.7 22 17.9 2.106 B.4 GRB091018 2009-10-18 2.4/2.4/2.4 1.0/0.9/0.9 2.1–1.8 1.0 3.5 19.1 0.971 B.5 GRB091127 2009-12-02 6.0/6.0/6.0 1.0/0.9/0.9 1.1–1.2 1.0 101 21.2 0.490 B.6 GRB100205A 2010-02-08 10.8/10.8/10.8 1.0/0.9/0.9 1.9–1.8 0.9 71 >24 – B.7 GRB100219A 2010-02-20 4.8/4.8/4.8 1.0/0.9/0.9 1.3–1.1 0.8 12.5 23 4.667 B.8 GRB100316B 2010-03-16 2.4/2.4/2.4 1.0/0.9/0.9 2.0–2.4 0.6 0.7 18.2 1.180 B.9 GRB100316D-1b 2010-03-17 3.6/3.6/3.6 1.0/0.9/0.9 1.2–1.3 0.8 10 21.5 0.059 B.10 GRB100316D-2 2010-03-19 2.4/2.4/2.4 1.0/0.9/0.9 1.1–1.2 0.9 58 20.2 0.059 B.10 GRB100316D-3 2010-03-20 2.6/2.6/3.2 1.0/0.9/0.9 1.1–1.2 1.1 79 19.9 0.059 B.10 GRB100316D-4 2010-03-21 2.6/2.6/3.2 1.0/0.9/0.9 1.1–1.2 1.5 101 19.9 0.059 B.10 GRB100418A-1 2010-04-19 4.8/4.8/4.8 1.0/0.9/0.9 1.6–1.3 0.7 8.4 18.1 0.624 B.11 GRB100418A-2 2010-04-20 4.8/4.8/4.8 1.0/0.9/0.9 1.2–1.3 0.6 34 19.2 0.624 B.11 GRB100418A-3 2010-04-21 4.8/4.8/4.8 1.0/0.9/0.9 1.2–1.4 0.7 58 >24 0.624 B.11 GRB100424Ac 2013-03-11 4.8/4.8/4.8 1.0/0.9/0.9 1.1–1.2 0.9 25239 >24 2.465 B.12 GRB100425A 2010-04-25 2.4/2.4/2.4 1.0/0.9/0.9 1.5–1.3 0.7 4 20.6 1.755 B.13 GRB100615Ac 2013-03-05 4.8/4.8/4.8 1.0/0.9/0.9 1.0–1.1 0.9 23859 >24 1.398 B.14 GRB100621A 2010-06-21 2.4/2.4/2.4 1.0/0.9/0.9 1.3–1.4 1.0 7.1 22 0.542 B.15 GRB100625Ac f 2010-07-07 4.8/4.8/4.8 1.0/0.9/0.9 1.1–1.0 0.8 278.7 >24 0.452 B.16 GRB100724Aa d 2010-07-24 4.2/4.2/4.2 1.0/0.9/0.9 1.5–2.3 0.7 0.2 19.52 1.288 B.17 GRB100728Be 2010-07-29 7.2/7.2/7.2 1.0/0.9/0.9 1.5–1.1 0.6 22 23 2.106 B.18 GRB100814A-1d 2010-08-14 0.9/0.9/0.9 1.0/0.9/0.9 1.9–1.7 0.5 0.9 19 1.439 B.19 GRB100814A-2 2010-08-14 4.8/4.8/4.8 1.0/0.9/0.9 1.5–1.2 0.7 2.1 19 1.439 B.19 GRB100814A-3 2010-08-18 4.8/4.8/4.8 1.0/0.9/0.9 1.2–1.0 0.6 98 20 1.439 B.19 GRB100816Af 2010-08-17 4.8/4.8/4.8 1.0/0.9/0.9 1.8–1.6 0.8 28.4 21.6 0.805 B.20 GRB100901A 2010-09-04 2.4/2.4/2.4 1.0/0.9/0.9 1.5–1.5 1.9 66 >24 1.408 B.21 GRB101219A 2010-12-19 7.2/7.2/7.2 1.0/0.9/0.9 1.1–1.7 1.8 3.7 >24 0.718 B.22 GRB101219B-1a 2010-12-20 4.8/4.8/4.8 1.0/0.9/0.9 1.6–2.6 1.4 11.6 20 0.552 B.23 GRB101219B-2a 2011-01-05 7.2/7.2/7.2 1.0/0.9/0.9 1.2–2.0 1.0 394 22.7 0.552 B.23 GRB101219B-3a 2011-01-25 7.2/7.2/7.2 1.0/0.9/0.9 1.4–2.1 0.7 886 >24 0.552 B.23 GRB110128A 2011-01-28 7.2/7.2/7.2 1.0/0.9/0.9 2.0–1.6 0.6 5.5 22.5 2.339 B.24 GRB110407A 2011-04-08 9.6/9.6/9.6 1.0/0.9/0.9 1.4–1.3 2.1 12.4 23 – B.25 GRB110709Bc 2013-03-19 7.2/7.2/7.2 1.0/0.9/0.9 1.6–1.1 0.9 14835 >24 2.109 B.26 GRB110715Aa 2011-07-16 0.6/0.6/0.6 1.0/0.9/0.9 1.1–1.1 1.6 12.3 18.5 0.823 B.27 GRB110721Aa 2011-07-22 2.4/2.4/2.4 1.0/0.9/0.9 1.2–1.4 2.3 28.7 >24 0.382 B.28 GRB110808A 2011-08-08 2.4/2.4/2.4 1.0/0.9/0.9 1.2–1.1 1.0 3.0 21.2 1.349 B.29 GRB110818A 2011-08-19 4.8/4.8/4.8 1.0/0.9/0.9 1.3–1.3 0.9 6.2 22.3 3.36 B.30 GRB111005Aa c 2013-04-01 1.2/1.2/1.2 1.0/0.9/0.9 1.3–1.3 0.7 13052 >24 0.013 B.31 GRB111008A-1 2011-10-09 8.8/8.8/8.4 1.0/0.9/0.9 1.1–1.0 1.3 8.5 21 4.990 B.32 GRB111008A-2 2011-10-10 8.0/8.0/7.2 1.0/0.9/0.9 1.3–1.0 0.9 20.1 22 4.990 B.32 GRB111107A 2011-11-07 4.8/4.8/4.8 1.0/0.9/0.9 1.8–1.5 0.8 5.3 21.5 2.893 B.33 GRB111117Af 2011-11-19 4.8/4.8/4.8 1.0/0.9/0.9 1.5–1.4 0.7 38 >24 2.211 B.34 GRB111123A-1 2011-11-24 6.2/6.6/6.6 1.0/0.9/0.9 1.6–1.1 0.8 12.2 >24 3.152 B.35 GRB111123A-2c 2013-03-07 2.4/2.4/2.4 1.0/0.9/0.9 1.0–1.0 0.5 11266 >24 3.152 B.35 GRB111129A 2011-11-30 3.6/3.6/3.6 1.0/0.9/0.9 1.6–2.1 1.9 8.7 >24 1.080 B.36 GRB111209A-1 2011-12-10 4.8/4.8/4.8 1.0/0.9/0.9 1.1–1.2 0.8 17.7 20.1 0.677 B.37 GRB111209A-2 2011-12-29 9.6/9.6/9.6 1.0/0.9/0.9 1.2–2.0 1.0 497 23 0.677 B.37 GRB111211Aa 2011-12-13 2.4/2.4/2.4 1.0/0.9/0.9 1.4–1.6 0.6 31 19.5 0.478 B.38 GRB111228A 2011-12-29 2.4/2.4/2.4 1.0/0.9/0.9 1.4–1.4 0.7 15.9 20.1 0.716 B.39 GRB120118Bc 2013-02-13 3.6/3.6/3.6 1.0/0.9/0.9 1.1–1.0 0.7 9393 >24 2.943 B.40 GRB120119A-1 2012-01-19 2.4/2.4/2.4 1.0/0.9/0.9 1.1–1.1 0.6 1.4 17 1.728 B.41 GRB120119A-2 2012-01-19 1.2/1.2/1.2 1.0/0.9/0.9 1.8–1.9 0.5 4.5 20 1.728 B.41 GRB120119A-3c 2013-02-26 4.8/4.8/4.8 1.0/0.9/0.6JH 1.0–1.1 1.8 9694 >24 1.728 B.41 GRB120211A-1c 2013-02-17 4.8/4.8/4.8 1.0/0.9/0.9 1.1–1.4 1.3 8919 >24 2.346 B.42 GRB120211A-2c 2013-03-20 3.6/3.6/3.6 1.0/0.9/0.9 1.1–1.2 1.2 9660 >24 2.346 B.42

Article number, page 21 of 43 A&A proofs: manuscript no. XSGRB_sample_arxiv

Table 1. The full sample of afterglows or hosts observed in the program (continued).

GRB Obs Date Exptime Slit width Airmass Seeing ∆t Magacq Redshift Notes (ks) (arcsec) (arcsec) (hr) GRB120224A 2012-02-25 2.4/2.4/2.4 1.0/0.9/0.9 1.7–2.1 1.3 19.8 22.3 1.10 B.43 GRB120311Aa 2012-03-11 2.4/2.4/2.4 1.0/0.9/0.9 1.6–1.4 0.7 3.7 21.6 0.350 B.44 GRB120327A-1a 2012-03-27 2.4/2.4/2.4 1.0/0.9/0.9 1.6–1.4 0.6 2.1 18.8 2.815 B.45 GRB120327A-2a 2012-03-28 4.2/4.2/4.2 1.0/0.9/0.9 1.0–1.1 0.6 29 22.5 2.815 B.45 GRB120404A 2012-04-05 9.6/9.6/9.6 1.0/0.9/0.9JH 1.7–1.3 1.3 15.7 21.3 2.876 B.46 GRB120422A 2012-04-22 4.8/4.8/4.8 1.0/0.9/0.9 1.3–1.3 0.7 16.5 22 0.283 B.47 GRB120712A 2012-07-13 4.8/4.8/4.8 1.0/0.9/0.9 1.5–2.5 1.5 10.4 21.5 4.175 B.48 GRB120714B 2012-07-15 4.8/4.8/4.8 1.0/0.9/0.9JH 1.5–1.2 1.2 7.8 22.1 0.398 B.49 GRB120716Aa 2012-07-19 3.6/3.6/3.6 1.0/0.9/0.9JH 1.8–2.6 1.1 62 20.9 2.486 B.50 GRB120722Ab 2012-07-22 4.8/4.8/4.8 1.0/0.9/0.9 1.3–1.3 1.2 10.3 23.6 0.959 B.51 GRB120805Ab 2012-08-14 3.6/3.6/3.6 1.0/0.9/0.9JH 1.3–1.7 0.9 218 >24 3.9 B.52 GRB120815Aa 2012-08-15 2.4/2.4/2.4 1.0/0.9/0.9 1.3–1.4 0.7 1.69 18.9 2.358 B.53 GRB120909Ad 2012-09-09 1.2/1.2/1.2 1.0/0.9/0.9 1.6–1.6 1.6 1.7 21 3.929 B.54 GRB120923A 2012-09-23 9.6/9.6/9.6 1.0/0.9/0.9JH 1.2–1.4 1.0 18.5 >24 7.84 B.55 GRB121024A 2012-10-24 2.4/2.4/2.4 1.0/0.9/0.9 1.2–1.1 0.6 1.8 20 2.300 B.56 GRB121027A 2012-10-30 8.4/8.4/8.4 1.0/0.9/0.9 1.3–1.3 1.3 69.4 21.15 1.773 B.57 GRB121201A 2012-12-02 4.8/4.8/4.8 1.0/0.9/0.9JH 1.1–1.1 1.1 12.9 23 3.385 B.58 GRB121229A 2012-12-29 4.8/4.8/4.8 1.0/0.9/0.9JH 1.4–1.2 1.5 2 21.5 2.707 B.59 GRB130131Bc 2013-03-09 7.2/7.2/7.2 1.0/0.9/0.9JH 1.3–1.6 1.1 874 >24 2.539 B.60 GRB130408Aa 2013-04-08 1.2/1.2/1.2 1.0/0.9/0.9 1.0–1.0 0.9 1.9 20 3.758 B.61 GRB130418A 2013-04-18 1.2/1.2/1.2 1.0/0.9/0.9 1.4–1.3 1.2 4.6 18.5 1.222 B.62 GRB130427A 2013-04-28 1.2/1.2/1.2 1.0/0.9/0.9JH 1.8–1.8 0.8 16.5 19 0.340 B.63 GRB130427B 2013-04-28 1.2/1.2/1.2 1.0/0.9/0.9JH 1.2–1.0 1.0 20.3 22.7 2.780 B.64 GRB130603Bf 2013-06-04 2.4/2.4/2.4 1.0/0.9/0.9 1.4–1.4 1.1 8.2 21.5 0.356 B.65 GRB130606A 2013-06-07 4.2/4.2/4.2 1.0/0.9/0.9JH 1.7–1.9 0.9 7.1 19 5.91 B.66 GRB130612A 2013-06-12 1.2/1.2/1.2 1.0/0.9/0.9 1.3–1.3 1.5 1.1 21.5 2.006 B.67 GRB130615A 2013-06-15 1.2/1.2/1.2 1.0/0.9/0.9 2.1–2.2 1.0 0.8 21 2.9 B.68 GRB130701A 2013-07-01 1.2/1.2/1.2 1.0/0.9/0.9JH 2.0–2.0 1.4 5.5 19.9 1.155 B.69 GRB130925A 2013-09-25 5.88/6.0/6.9 1.0/0.9/0.9JH 1.0–1.0 0.6 3.5 >24 0.347 B.70 GRB131011Aa 2013-10-13 4.5/4.5/4.5 1.0/0.9/0.9 1.1–1.1 0.8 34.2 >24 1.874 B.71 GRB131030A 2013-10-31 3.6/3.6/3.6 1.0/0.9/0.9 1.1–1.1 1.1 3.4 18.0 1.296 B.72 GRB131103A 2013-11-05 2.4/2.4/2.4 1.0/0.9/0.9JH 1.1–1.1 1.0 5.8 20.48 0.599 B.73 GRB131105A 2013-11-05 4.8/4.8/4.8 1.0/0.9/0.9 1.3–1.4 0.8 1.3 22.4 1.686 B.74 GRB131117A 2013-11-17 4.8/4.8/4.8 1.0/0.9/0.9JH 1.3–1.2 1.7 1.1 20 4.042 B.75 GRB131231Aa 2014-01-01 2.4/2.4/2.4 1.0/0.9/0.9JH 1.4–1.3 0.9 20.2 18.5 0.642 B.76 GRB140114Ac 2014-03-28 5.4/5.4/5.4 1.0/0.9/0.9JH 1.7–1.7 1.2 1746 >24 3.0 B.77 GRB140213Aa 2014-02-14 1.2/1.2/1.2 1.0/0.9/0.9JH 1.5–1.5 0.7 5.8 19.5 1.208 B.78 GRB140301A 2014-03-02 7.2/7.2/7.2 1.0/0.9/0.9JH 1.1–1.1 0.9 9 23.1 1.416 B.79 GRB140311Aa 2014-03-13 7.6/6.3/8.4 1.0/0.9/0.9JH 1.2–1.2 0.6 32.5 >24 4.954 B.80 GRB140430Aa 2014-04-30 1.2/1.2/1.2 1.0/0.9/0.9 2.0–1.8 1.6 2.5 19 1.601 B.81 GRB140506A-1 2014-05-07 4.8/4.8/4.8 1.0/0.9/0.9 1.3–1.4 0.7 8.8 20.9 0.889 B.82 GRB140506A-2 2014-05-08 4.8/4.8/4.8 1.0/0.9/0.9 1.2–1.3 0.7 32.9 0.889 B.82 GRB140515A 2014-05-16 4.8/4.8/4.8 1.0/0.9/0.9 1.3–1.3 1.4 15.5 >24 6.327 B.83 GRB140614A 2014-06-14 2.4/2.4/2.4 1.0/0.9/0.9 1.8–1.8 0.7 3.8 21.5 4.233 B.84 GRB140622Af 2014-06-22 1.2/1.2/1.2 1.0/0.9/0.9 1.4–1.3 1.0 0.8 >24 0.959 B.85 GRB141028Aa 2014-10-29 2.4/2.4/2.4 1.0/0.9/0.9 1.5–1.4 1.0 15.4 20 2.332 B.86 GRB141031Aa c 2015-01-29 2.4/2.4/2.4 1.0/0.9/0.9 1.2–1.3 0.8 10912 >24 – B.87 GRB141109A-1 2014-11-09 2.4/2.4/2.4 1.0/0.9/0.9JH 1.5–1.7 0.8 1.9 19.2 2.993 B.88 GRB141109A-2 2014-11-10 4.3/4.3/4.5 1.0/0.9/0.9JH 1.7–2.0 0.8 25.4 2.993 B.88 GRB150206Aa 2015-02-07 2.4/2.4/2.4 1.0/0.9/0.9 2.1–1.9 0.8 10 21.9 2.087 B.89 GRB150301B 2015-03-02 3.6/3.6/3.6 1.0/0.9/0.9JH 1.2–1.2 1.1 5.1 21.0 1.517 B.90 GRB150403A 2015-04-04 2.4/2.4/2.4 1.0/0.9/0.9 1.6–1.7 0.7 10.8 19.1 2.057 B.91 GRB150423Ad f 2015-04-23 4.8/4.8/4.8 1.0/0.9/0.9 2.7–2.4 1.4 0.4 >24 1.394 B.92 GRB150428A 2015-04-28 2.4/2.4/2.4 1.0/0.9/0.9JH 1.6–1.5 0.8 3.7 >24 – B.93 GRB150514Aa 2015-05-15 2.4/2.4/2.4 1.0/0.9/0.9 2.3–2.1 0.9 28.4 19.5 0.807 B.94 GRB150518Aa 2015-05-20 2.4/2.4/2.4 1.0/0.9/0.9JH 1.3–1.3 1.7 30.7 >24 0.256 B.95 GRB150616Aa c 2015-09-12 2.4/2.4/2.4 1.0/0.9/0.9JH 1.2–1.1 1.2 2092 >24 1.188 B.96 GRB150727A 2015-07-28 3.6/3.6/3.6 1.0/0.9/0.9JH 1.2–1.2 1.4 5.0 20.5 0.313 B.97 GRB150821Ad 2015-08-21 2.4/2.4/2.4 1.0/0.9/0.9 2.0–1.8 1.3 0.2 16 0.755 B.98

Article number, page 22 of 43 Selsing et al.: The X-shooter GRB afterglow legacy sample

Table 1. The full sample of afterglows or hosts observed in the program (continued).

GRB Obs Date Exptime Slit width Airmass Seeing ∆t Magacq Redshift Notes (ks) (arcsec) (arcsec) (hr) GRB150910A 2015-09-11 1.8/1.8/1.8 1.0/0.9/0.9JH 1.9–1.9 1.3 20.1 21.2 1.359 B.99 GRB150915A 2015-09-16 4.8/4.8/4.8 1.0/0.9/0.9JH 1.1–1.1 1.6 3.3 23 1.968 B.100 GRB151021Ad 2015-10-21 4.2/4.2/4.2 1.0/0.9/0.9 1.0–1.1 1.4 0.75 18.2 2.330 B.101 GRB151027B 2015-10-28 2.4/2.4/2.4 1.0/0.9/0.9JH 1.5–1.7 1.2 5 20.5 4.063 B.102 GRB151029A 2015-10-29 1.2/1.2/1.2 1.0/0.9/0.9JH 1.9–1.7 1.1 1 20 1.423 B.103 GRB151031Ad 2015-10-31 4.2/4.2/4.2 1.0/0.9/0.9 1.1–1.1 1.1 0.3 20.4 1.167 B.104 GRB160117B 2016-01-18 4.8/4.8/4.8 1.0/0.9/0.9JH 1.1–1.2 1.1 13.5 20.8 0.870 B.105 GRB160203Ad 2016-02-03 6.6/6.6/6.6 1.0/0.9/0.9 1.0–1.8 1.0 0.3 18 3.518 B.106 GRB160228Ac 2016-03-12 4.8/4.8/4.8 1.0/0.9/0.9JH 1.7–1.7 1.0 296 >24 1.640 B.107 GRB160303Af 2016-03-04 4.8/4.8/4.8 1.0/0.9/0.9JH 1.6–1.5 0.8 19.1 >24 – B.108 GRB160314A 2016-03-15 4.8/4.8/4.8 1.0/0.9/0.9JH 1.3–1.3 0.8 13.0 21.7 0.726 B.109 GRB160410Ad f 2016-04-10 1.8/1.8/1.8 1.0/0.9/0.9 2.5–2.3 0.5 0.15 20.3 1.717 B.110 GRB160425A 2016-04-26 4.8/4.8/4.8 1.0/0.9/0.9JH 1.3–1.3 0.5 7.2 21.1 0.555 B.111 GRB160625Ba 2016-06-27 2.4/2.4/2.4 1.0/0.9/0.9JH 1.3–1.3 0.7 30 19.1 1.406 B.112 GRB160804A-1a 2016-08-04 2.4/2.4/2.4 1.0/0.9/0.9JH 1.4–1.3 0.6 22.4 21.2 0.736 B.113 GRB160804A-2a c 2016-08-27 3.6/3.6/3.6 1.0/0.9/0.9JH 1.9–1.8 0.6 574 >24 0.736 B.113 GRB161001A 2016-10-01 2.4/2.4/2.4 1.0/0.9/0.9JH 1.2–1.3 0.5 6.1 >24 0.891 B.114 GRB161007Ac 2016-10-14 2.4/2.4/2.4 1.0/0.9/0.9JH 1.6–1.6 0.7 323 >24 – B.115 GRB161014A 2016-10-15 4.8/4.8/4.8 1.0/0.9/0.9JH 1.1–1.2 0.5 11.6 21.4 2.823 B.116 GRB161023Aa 2016-10-24 1.2/1.2/1.2 1.0/0.9/0.9JH 1.2–1.2 0.9 3 17.5 2.710 B.117 GRB161117A 2016-11-17 2.4/2.4/2.4 1.0/0.9/0.9 1.8–1.6 2.6 0.73 19 1.549 B.118 GRB161219B 2016-12-21 2.4/2.4/2.4 1.0/0.9/0.9JH 1.1–1.1 0.9 35.7 19.5 0.146 B.119 GRB170113A 2017-01-14 4.8/4.8/4.8 1.0/0.9/0.9JH 1.5–1.4 0.9 15.23 21.7 1.968 B.120 GRB170202A 2017-02-03 2.4/2.4/2.4 1.0/0.9/0.9JH 1.3–1.2 0.7 9.7 20.8 3.645 B.121

Notes. (a) Not part of the statistical sample (b) Spectrum dominated by light from the host galaxy (c) Spectrum of the host galaxy taken long after the burst (d) RRM observation (e) ADC malfunction during observation (f) Short burst

Article number, page 23 of 43 A&A proofs: manuscript no. XSGRB_sample_arxiv

GRB090809A GRB090926A 1.2 1.2

1.0 1.0

0.8 0.8

0.6 0.6

0.4 0.4 Normalised Flux Normalised Flux 0.2 0.2

0.0 0.0

4200 4400 4600 4800 5000 3600 3700 3800 3900 4000 Observed Wavelength (Å) Observed Wavelength (Å) GRB 100219A GRB 100425A 1.2 1.2

1.0 1.0

0.8 0.8

0.6 0.6

0.4 0.4 Normalised Flux Normalised Flux 0.2 0.2

0.0 0.0

6600 6800 7000 7200 3250 3300 3350 3400 3450 Observed Wavelength (Å) Observed Wavelength (Å) GRB100728B GRB 110128A 1.2 1.2

1.0 1.0

0.8 0.8

0.6 0.6

0.4 0.4 Normalised Flux Normalised Flux 0.2 0.2

0.0 0.0

3400 3600 3800 4000 3800 3900 4000 4100 4200 4300 Observed Wavelength (Å) Observed Wavelength (Å) GRB110818A GRB 111008A 1.2 1.2 1.0 1.0 0.8 0.8 0.6 0.6 0.4 0.4 Normalised Flux Normalised Flux 0.2 0.2

0.0 0.0

5000 5500 6000 6500 6500 7000 7500 8000 Observed Wavelength (Å) Observed Wavelength (Å)

Article number, page 24 of 43

Fig. 10. Measurements of the hydrogen column densities for all bursts with a clear Lyman alpha absorption system. In solid black is shown the spectrum with black dotted giving the corresponding 1-σ error. The black dashed line shows zero ﬂux density. The solid red line is the absorption of column density equal to the value presented in Tab.4 with the 1- σ interval shown with dashed lines. Selsing et al.: The X-shooter GRB afterglow legacy sample

GRB111107A GRB 120119A 1.2 1.2

1.0 1.0

0.8 0.8

0.6 0.6

0.4 0.4 Normalised Flux Normalised Flux 0.2 0.2

0.0 0.0

4600 4700 4800 4900 3200 3400 3600 3800 4000 Observed Wavelength (Å) Observed Wavelength (Å) GRB 120327A GRB120404A 1.2 1.2

1.0 1.0

0.8 0.8

0.6 0.6

0.4 0.4 Normalised Flux Normalised Flux 0.2 0.2

0.0 0.0

4200 4400 4600 4800 5000 4500 4600 4700 4800 4900 Observed Wavelength (Å) Observed Wavelength (Å) GRB120712A GRB 120716A 1.2 1.2

1.0 1.0

0.8 0.8

0.6 0.6

0.4 0.4 Normalised Flux Normalised Flux 0.2 0.2

0.0 0.0

6200 6250 6300 6350 3800 4000 4200 4400 4600 Observed Wavelength (Å) Observed Wavelength (Å) GRB 120815A GRB 120909A 1.2 1.2

1.0 1.0

0.8 0.8

0.6 0.6

0.4 0.4 Normalised Flux Normalised Flux 0.2 0.2

0.0 0.0

3600 3800 4000 4200 4400 5600 5800 6000 6200 6400 Observed Wavelength (Å) Observed Wavelength (Å)

Article number, page 25 of 43

Fig. 11. * Fig. 10. continued. A&A proofs: manuscript no. XSGRB_sample_arxiv

GRB121024A GRB121027A 1.2 1.2

1.0 1.0

0.8 0.8

0.6 0.6

0.4 0.4 Normalised Flux Normalised Flux 0.2 0.2

0.0 0.0

3700 3800 3900 4000 4100 4200 4300 3500 4000 4500 Observed Wavelength (Å) Observed Wavelength (Å) GRB121201A GRB121229A 1.2 1.2 1.0 1.0 0.8 0.8 0.6 0.6

0.4 0.4 Normalised Flux Normalised Flux 0.2 0.2

0.0 0.0

4800 5000 5200 5400 5600 5800 4200 4400 4600 4800 Observed Wavelength (Å) Observed Wavelength (Å) GRB130408A GRB130427B 1.2 1.2

1.0 1.0

0.8 0.8

0.6 0.6

0.4 0.4 Normalised Flux Normalised Flux 0.2 0.2

0.0 0.0

5600 5800 6000 6200 4400 4600 4800 5000 5200 Observed Wavelength (Å) Observed Wavelength (Å) GRB130606A GRB130612A 1.2 1.2

1.0 1.0

0.8 0.8

0.6 0.6

0.4 0.4 Normalised Flux Normalised Flux 0.2 0.2

0.0 0.0

8350 8400 8450 8500 3400 3600 3800 4000 Observed Wavelength (Å) Observed Wavelength (Å)

Article number, page 26 of 43

Fig. 12. * Fig. 10. continued. Selsing et al.: The X-shooter GRB afterglow legacy sample

GRB131011A GRB131117A 1.5 1.2

1.0 1.0 0.8

0.6

0.4 0.5 Normalised Flux Normalised Flux 0.2

0.0 0.0

3400 3600 3800 4000 6000 6050 6100 6150 6200 6250 Observed Wavelength (Å) Observed Wavelength (Å) GRB140430A 1.2

1.0

0.8

0.6

0.4 Normalised Flux 0.2

0.0

3200 3300 3400 3500 Observed Wavelength (Å) GRB140515A GRB140614A 1.5 1.2

1.0 1.0 0.8

0.6

0.4 0.5 Normalised Flux Normalised Flux 0.2

0.0 0.0 8840 8860 8880 8900 8920 8940 8960 6000 6200 6400 6600 6800 Observed Wavelength (Å) Observed Wavelength (Å) GRB141028A GRB141109AA 1.2 1.2

1.0 1.0

0.8 0.8

0.6 0.6

0.4 0.4 Normalised Flux Normalised Flux 0.2 0.2

0.0 0.0

3900 3950 4000 4050 4100 4150 4200 4600 4800 5000 5200 Observed Wavelength (Å) Observed Wavelength (Å)

Article number, page 27 of 43

Fig. 13. * Fig. 10. continued. A&A proofs: manuscript no. XSGRB_sample_arxiv

GRB150206A GRB150403A 1.5 1.2

1.0 1.0 0.8

0.6

0.5 0.4 Normalised Flux Normalised Flux 0.2

0.0 0.0 3600 3800 4000 4200 3500 3600 3700 3800 3900 4000 Observed Wavelength (Å) Observed Wavelength (Å) GRB150915A GRB 151021A 1.5 1.2

1.0 1.0 0.8

0.6

0.5 0.4 Normalised Flux Normalised Flux 0.2

0.0 0.0

3400 3500 3600 3700 3800 3900 3600 3800 4000 4200 4400 Observed Wavelength (Å) Observed Wavelength (Å) GRB151027B GRB160203A 1.2 1.2 1.0 1.0 0.8 0.8 0.6 0.6

0.4 0.4 Normalised Flux Normalised Flux 0.2 0.2

0.0 0.0

6000 6100 6200 6300 5300 5400 5500 5600 5700 Observed Wavelength (Å) Observed Wavelength (Å) GRB 160410A GRB161014A 1.2 1.2

1.0 1.0

0.8 0.8

0.6 0.6

0.4 0.4 Normalised Flux Normalised Flux 0.2 0.2

0.0 0.0

3200 3250 3300 3350 3400 4400 4600 4800 5000 Observed Wavelength (Å) Observed Wavelength (Å)

Article number, page 28 of 43

Fig. 14. * Fig. 10. continued. Selsing et al.: The X-shooter GRB afterglow legacy sample

GRB161023A GRB170202A

1.0 1.2

1.0 0.8 0.8 0.6 0.6 0.4 Normalised Flux Normalised Flux 0.4

0.2 0.2

0.0 0.0 4400 4450 4500 4550 4600 5600 5700 5800 Observed Wavelength (Å) Observed Wavelength (Å)

Article number, page 29 of 43

Fig. 15. * Fig. 10. continued. A&A proofs: manuscript no. XSGRB_sample_arxiv

Appendix A: The complex error function and the 45 hours after the BAT trigger reveals a wealth of absorption Voigt profile features superposed on the afterglow continuum at a common redshift of z = 3.373. Two intervening systems at z = 1.959 When modeling the spectral PSF, we need to evaluate the Voigt- and z = 1.800 are identified based on strong Mg ii-absorption. profile. The Voigt profile, which is the convolution of the Gaus- Because this burst is observed before the instrument is science- sian and Lorentzian profiles, can, centered at zero, be written as verified, it does not enter into the statistical sample. (Pagnini & Mainardi 2010)

V(λ, σ, γ) = G(λ, σ) ⊗ L(λ, γ) Appendix B.2: GRB 090530 (z = 1.266) Z ∞ = G(ξ, σ)L(λ − ξ, γ)dξ Observed during paranalization of the instrument, these data −∞ forms the basis of GCN 15571 (Goldoni et al. 2013), but are 2 Z ∞ ξ 2 not published elsewhere. Observations began 20.6 hours after 1 − √ 1 γ = √ e 2σ dξ the BAT trigger and 4.8 ks spectroscopic integration in all three γπ λ − ξ 2 γ2 −∞ 2πσ ( ) + arms reveals the absorption signature for a host at z = 1.266 from 2 ξ Z ∞ − √ the detection of Mg ii, Mg i, Si ii, Fe ii, Al iii. Because this burst γ 1 e 2σ √ dξ. was observed before the instrument was science-veriﬁed, it does = 3/2 2 2 (A.1) 2σ π −∞ (λ − ξ) + γ not enter into the statistical sample. √ We can by making the following substitution, ξ = 2σ and √ Appendix B.3: GRB 090809 (z = 2.737) dξ = 2σdt, write it as

√ γ Observed during the first science verification period and was re- 2σ √ Z ∞ e−t2 ported in GCN 9761 (Malesani et al. 2009a) and is additionally V(λ, σ, γ) = √ 2σ √ dt π π −∞ (λ − 2σt)2 + γ2 used as the basis for the master thesis by Ása Skúladóttir (2010). γ 7.2 ks integration starting 10.2 hours after the GRB trigger yields √ Z ∞ −t2 1 2σ e a clear afterglow continuum in all arms. The simultaneous de- = √ dt. (A.2) 2 2 ii i i∗ iv 2πσ π −∞ γ tection of absorption lines identified as Lyα, Si ,O , Si , Si , √λ − t + √ 2σ 2σ C iv, Fe ii, Al ii, Al iii and Mg ii at z = 2.737 sets it as the redshift of the GRB. Because this burst is observed before the instrument This form of the convolution is closely related to the complex was science-verified, it does not enter into the statistical sample. probability function (Letchworth & Benner 2007; Abrarov & Quine 2015b), Appendix B.4: GRB 090926A (z = 2.106) i Z ∞ e−t2 W(z) = dt (A.3) Obtained during the second science verification period, this π −∞ z − t dataset forms the basis of GCN 9942 (Malesani et al. 2009b) and is additionally published in D’Elia et al.(2010). Spectroscopic for any complex argument, z = x + iy. The complex probability integration started 22 hours after the BAT trigger and from the function can be expressed as a sum of a real and an imaginary acquisition camera the optical afterglow is measured to R = 17.9 part (Benner et al. 1995; Abrarov & Quine 2015a), mag at the beginning of the observations which causes a strong continuum to be seen in all arms. An absorption trough due to W(x, y) = K(x, y) + iL(x, y) Lyα is clearly visible along with numerous metal resonance lines y Z ∞ e−t2 i Z ∞ (x − t)e−t2 C iv, Si ii, Si ii∗ Fe ii, Mg ii, all at z = 2.106. Because this burst = dt + dt, π (x − t)2 + y2 π (x − t)2 + y2 was observed before the instrument was science-verified, it does −∞ −∞ not enter into the statistical sample. (A.4) √ where the real part, Re[W(x, y)] = 2πσV(λ, σ, γ) if x = √λ 2σ Appendix B.5: GRB 091018 (z = 0.971) γ and y = √ , which can be obtained by using the complex 2σ The first burst observed during normal operation after science λ iγ argument, z = √+ , in the complex probability function. If verification was completed and the first burst that enters the sta- 2σ Im[z] ≥ 0, which is always guaranteed for the width of a spec- tistical sample. These data are the basis for GCN 10042 (de tral profile, the complex probability function equals the com- Ugarte Postigo et al. 2009b) and is published in Wiersema et al. plex error function. The complex error function has numerous, (2012). With a bright afterglow and a rapid follow-up, this spec- fast, numerical approximations where in this work we use the trum is of pristine quality. The afterglow continuum is bright scipy.special.wofz (Oliphant 2007) implementation. throughout all spectroscopic arms which allows the ready detection of Al ii, Al iii, Fe ii, Mn ii, Mg ii, Mg i, and Ca ii - all located at z = 0.971, setting it as the redshift of the GRB. Appendix B: Notes on Individual objects Appendix B.1: GRB 090313 (z = 3.373) Appendix B.6: GRB 091127 (z = 0.490) The first GRB ever observed with X-shooter, during the com- Obtained 4 days after the burst trigger, these data forms the basis missioning of the instrument, these data formed the basis of for GCN 10233 (Thöne et al. 2009) and are published in Ver- GCN 9015 (de Ugarte Postigo et al. 2009a) and are published gani et al.(2011). Due to the late follow-up and a nearby moon, in de Ugarte Postigo et al.(2010). Due to the lingering bright- the signal-to-noise of the afterglow continuum is low, especially ness of GRB 090313, 6.9 ks spectroscopic integration starting in the UVB arm. This is why no clear absorption lines are de-

Article number, page 30 of 43 Selsing et al.: The X-shooter GRB afterglow legacy sample tected against the afterglow continuum, although see Vergani 1.4, and 2.4 days after the burst, each lasting 4.8 ks. The unam- et al.(2011) which report a tentative detection of Mg ii. Emis- biguous redshift of the host, z = 0.624, is found from the si- sion lines from the underlying host are clearly visible with lines multaneous detection of emission features belonging to nebular from [O ii], Hβ, [O iii], and Hα all at z = 0.490. This bursts is ad- lines, including H i, [O ii], [O iii], [Ne iii], N ii, [S ii], [S iii], and ditionally associated with SN2009nz (Cobb et al. 2010; Berger He i as well as absorption features due to the presence of Zn ii, et al. 2011; Olivares E. et al. 2015). Cr ii, Fe ii, Mn ii, Mg ii, Mg i, Ti ii, and Ca ii, all at a consistent redshift. Temporal evolution of the ﬁne structure lines belonging to Fe ii∗ is found between the epochs. Appendix B.7: GRB 100205A (z = na) Observed 3 days after the Swift trigger. No afterglow or host detected is in 10.8 ks. The GRB is likely located at high redshift Appendix B.12: GRB 100424A (z=2.465) (GCN GCN10399; Perley et al. 2010). The spectrum has not oth- The data presented here also formed the basis of GCN 14291 erwise been published previously. (Malesani et al. 2013b). The spectrum has not otherwise been published. Observations were carried out long after the burst af- Appendix B.8: GRB 100219A (z = 4.667) terglow had faded. Emission lines from the host are detected at z = 2.465. These data presented here have also formed the basis of GCN 10441 (Groot et al. 2010) and are published in Thöne et al. (2013). Observations started 12.5 hours after the Swift trigger Appendix B.13: GRB 100425A (z=1.1755) and has a total exposure time of 4.8 ks. Absorption features, including those of Lyα and from a multitude of ions are detected The spectra presented here also formed the basis of GCN 10684 against the afterglow continuum at z = 4.667. Additionally, ab- (Goldoni et al. 2010) and are used in Skuladottir (2010), but sorption from an intervening system is found at z = 2.181. not published elsewhere. Observations started 4 hours after the Swift trigger, totaling 2.4 ks. Absorption features from Mg ii and Fe ii in the afterglow continuum are detected at z = 1.1755. Appendix B.9: GRB 100316B (z = 1.180) The data presented here also formed the basis of GCN 10495 Appendix B.14: GRB 100615A (z=1.398) (Vergani et al. 2010b). The spectrum has not otherwise been published. Observations started 44 minutes after the Swift trig- The data presented here also formed the basis of GCN 14264 ger and have a total exposure time of 2.4 ks. Absorption features (Krühler et al. 2013a), but are not published elsewhere. Host ob- from Fe ii, Al ii, Al iii, Mg ii and Mg i are well detected against servation of a dark burst (D’Elia & Stratta 2011) taken long after the afterglow continuum at z = 1.180. Additionally, strong ab- the afterglow had faded. Emission lines from the host belonging sorption lines from Fe ii and Mg ii from an intervening system to [O ii], [Ne iii], [O iii] and Hα are detected at a common redshift are found at z = 1.063. of z = 1.398.

Appendix B.10: GRB 100316D (z = 0.059) Appendix B.15: GRB 100621A (z=0.542) The data presented here also formed the basis of GCN 10512 (Vergani et al. 2010a), GCN 10513 (Vergani et al. 2010c), GCN The data presented here also formed the basis of GCN 10876 10543 (Bufano et al. 2010) and are published in Bufano et al. (Milvang-Jensen et al. 2010), but are not published elsewhere. (2012) and Starling et al.(2011). This GRB is very close by and Beginning 7.1 hours after the GRB, 2.4 ks observations reveal has an associated SN, SN2010bh, and has therefore undergone emission lines from [O ii], Hβ and [O iii] at a common redshift intense follow-up (Olivares E. et al. 2012; Cano et al. 2011; Izzo of z = 0.542 and a very weak afterglow continuum. et al. 2017). The spectra presented here consists of a subset of the entire VLT/X-shooter campaign, covering the four first observing days while the afterglow still contributes significantly to Appendix B.16: GRB 100625A (z=0.452) the total emission. The first observations started 10 hours after The data presented here is of the candidate host galaxy, taken the burst, before the SN was discovered, and targeted the star- long after the burst had faded and have not previously been pub- forming ’A’-region (Starling et al. 2011), not the GRB. A very lished. 4.8 ks of exposure reveals a weak continuum present in rich spectrum containing a multitude of emission lines puts the all arms, but an absence of emission lines. This could indicate host at z = 0.059. For three consecutive nights, 58, 79 and 101 that the host primarily contains an older stellar population. The hours after the Swift trigger, the afterglow was observed as it redshift, z = 0.452, is taken from Fong et al.(2013). transitioned into the spectrum of a high-velocity Ic-BL SN. The observations taken 79 and 101 hours after the burst are taken under programme 084.D-0265(A) (PI: Benetti), but with an identi- Appendix B.17: GRB 100724A (z = 1.288 cal setup to the first two observations. The data presented here also formed the basis of GCN 10971 Appendix B.11: GRB 100418A (z=0.624) (Thöne et al. 2010). The spectrum has not otherwise been published previously. The observations were carried out in RRM The data presented here also formed the basis of GCN 10620 starting 11 min after the GRB trigger. See section 2.3, for a (Antonelli et al. 2010) and GCN 10631 (Malesani & Palazzi description of the RRM scheme. Absorption lines from several 2010) and are published in De Ugarte Postigo et al.(2011). The ionic species are detected in the afterglow continuum at a com- burst has been followed up in three epochs of observations, 0.4, mon redshift of z = 1.288.

Article number, page 31 of 43 A&A proofs: manuscript no. XSGRB_sample_arxiv

Appendix B.18: GRB 100728B (z=2.106) burst trigger and lasting 4.8 ks, reveals absorption from Mg ii and Mg i in the host located at z = 0.552 on a featureless continuum The data presented here also formed the basis of GCN 11317 visible across the entire coverage of X-shooter. Subsequent ob- (Flores et al. 2010). The spectrum has not otherwise been pub- servations taken 16 and 37 days after the trigger shows the fading lished previously. Starting 22 hours after the burst trigger, 7.2 ks spectral signature of a SN, SN2010ma. of observations reveals a faint afterglow continuum with Lyα- and Mg ii-absorption at z = 2.106. Due to a malfunctioning ADC, the sensitivity of X-shooter is depressed with respect to Appendix B.24: GRB 110128A (z=2.339) normal operations, resulting in a lower throughout. Additionally, the position of the trace on the slit moves due to atmospheric dif- These observations form the basis of GCN 11607 (Sparre et al. ferential refraction. The presence of the DLA is confirmed in the 2011), but have not been published before. Spectroscopic inte- 2D image and despite the observational challenges that affects gration started 6.55 hours after the Swift trigger and lasted for a −2 this spectrum, we measure log(NHI/cm ) = 21.2 ± 0.15. total of 7.2 ks. The afterglow continuum is detected across the entire spectral coverage at moderate signal-to-noise. Absorption lines in the continuum are detected from Lyα,O i,C ii, Si iv,C iv, Appendix B.19: GRB 100814A (z=1.439) Si ii and Fe ii, all at a common redshift of z = 2.339. From the −2 The spectra presented here has not been published previously. broad Lyα trough, a hydrogen column density log(NHI/cm ) = The observations consist of three visits, the first beginning only 22.6 ± 0.2 is derived. An intervening system at z = 2.20 is tenta- 0.9 hours after the Swift trigger, the other two visits were 2.13 tively identified from an absorption feature, likely due to C iv. and 98.40 hours after the trigger, respectively. A bright afterglow continuum is present in all visits, allowing identification of absorption features belonging to a wide range of ions at z = 1.439. Appendix B.25: GRB 110407A (z=na) A complex velocity structure in the absorption features belonging to Mg ii, shows several components, separated by as much These observations have not been published before. Starting as 500 km/s, pointing to a likely merger scenario in the host or 12.36 hours after the BAT trigger, 4.8 ks spectroscopic integra- starburst driven outflows. tion yields a very faint trace down to ∼430 nm, only visible after binning heavily. This could indicate a redshift, z ∼ 2.5, but no emission lines or absorption lines are immediately visible to sup- Appendix B.20: GRB 100816A (z=0.805) port this. The data presented here also formed the basis of GCN 11123 (Tanvir et al. 2010). The spectrum has not otherwise been pub- Appendix B.26: GRB 110709B (z=2.109) lished previously. This short GRB was observed 28.4 hours after the GRB trigger. 4 x 1200 s of exposure reveals two distinct sets This is a late-time observation (> 1 year) and has previously been 00 of emission lines, spatially offset . 1 , very close in redshift used in Perley et al.(2016a). In this reduction of the 7.2 ks spec- space, z = 0.8034 and z = 0.8049, indicating either an interact- troscopic integration, the tentative detection of [O iii] reported in ing host or some complex velocity structure of the host. Faint un- Perley et al.(2016a) is confirmed along with low-significance derlying continua are present under both sets of emission lines. detection of Hα at the end of the spectral coverage, both at a consistent redshift, z = 2.109, securing it as the redshift of the Appendix B.21: GRB 100901A (z=1.408) GRB host. The data presented here have been published in Hartoog et al. (2013). Because of the unusual lingering brightness of this GRB, Appendix B.27: GRB 110715A (z=0.823) 2.4 ks of observations taken 65.98 hours after the GRB trigger still reveals an afterglow continuum visible across the entire These observations, starting 12.3 hours after the trigger, have spectral coverage of X-shooter. Absorption lines from a wide been published in Sánchez-Ramírez et al.(2017) and addition- range of ions sets the redshift at z = 1.408, with intervening ally formed the basis of GCN 12164 (Piranomonte et al. 2011). absorption systems at z = 1.3147 and z = 1.3179. Only a single exposure of 600 s was obtained, before strong winds interrupted the observations. A red continuum is detected across all arms and a multitude of absorption lines are super- Appendix B.22: GRB 101219A (z=0.718) posed on the afterglow continuum. We identify lines belonging to Al ii, Al iii, Zn ii, Cr ii, Fe ii, Mg ii, Mg i, Ca ii, and Ca ii, all at These data have not been published before. Starting 3.7 hours z = 0.823, marking it as the redshift of the GRB. after the GRB trigger, 7.2 ks of exposure time reveal a very faint continuum in the visual and near-infrared, only visible when heavily binning the images. No redshift estimate is available Appendix B.28: GRB 110721A (z=0.382) from these observations. Late-time Gemini-North observations reveal emission lines from the host at z = 0.718 (Chornock & This is a Fermi burst with a LAT detection and thus outside the Berger 2011). statistical sample, but nonetheless followed up due to the extremely high peak energy (Axelsson et al. 2012). Starting 28.7 Appendix B.23: GRB 101219B (z=0.552) hours after the burst trigger, 2.4 ks spectroscopic observation reveals after heavy binning, a wide, faint trace down to ∼ 580 nm, The data presented here also formed the basis of GCN 11579 offset by 2.500 relative to the centering of the slit. No good red- (de Ugarte Postigo et al. 2011b) and are published in Sparre shift measurement can be inferred from this. We have adopted et al.(2011). The first observation, taken 11.6 hours after the the redshift from GCN 12193 (Berger 2011).

Article number, page 32 of 43 Selsing et al.: The X-shooter GRB afterglow legacy sample

Appendix B.29: GRB 110808A (z=1.348) Appendix B.34: GRB 111117A (z=2.211)

This spectrum has already formed the basis of GCN 12258 These data have previously been used to form some of the basis (de Ugarte Postigo et al. 2011a), but is not published other- of Selsing et al.(2018). Starting 37.3 hours after the BAT trig- wise. Starting 3 hours after the Swift trigger, a rich spectrum ger, 4.8 ks of spectroscopic integration yields faint emission lines is obtained in 2.4 ks spectroscopic integration. The GRB after- identified at [O ii], Hβ, [O iii] and Hα, all at a common z = 2.211, glow continuum is visible across all three spectrscopic arms of marking it as the redshift of the GRB host. No afterglow contin- VLT/X-shooter with emission lines identified as [O ii], [O iii], uum is detected. This is the highest redshift short GRB detected Hα all at z = 1.348. At the same redshift, we identify absorp- to date. tion lines superposed on the afterglow continuum from Mg ii and Fe ii. Appendix B.35: GRB 111123A (z=3.151) These data formed the basis of GCN 14273 (Xu et al. 2013c), Appendix B.30: GRB 110818A (z=3.36) but are not published elsewhere. Observed twice, the first time shortly after the GRB and the second long after the burst had Starting 6.15 hours after the BAT trigger, spectroscopic integra- faded, securely sets the redshift of the host at z = 3.151 based on tion for 4.8 ks reveals a moderate signal-to-noise GRB afterglow the detection of emission lines identified as [O ii] and [O iii]. continuum, down to ∼ 500 nm. The simultaneous detection of absorption features identified as Lyα, Si ii,C iv, Al ii, Ca H, Ca K, and Mg ii, and emission from the [O iii]-doublet, securely sets Appendix B.36: GRB 111129A (z=1.080) z = 3.36 as the redshift of the GRB. These data form the ba- Starting 8.26 hours after the GRB trigger and lasting 3.6 ks, these sis of GCN 12284 (D’Avanzo et al. 2011a), but is not published observations have previously been published in Krühler et al. elsewhere. (2015). A very faint continuum is visible after severe binning and a redshift is suggested in Krühler et al.(2015), based on the Appendix B.31: GRB 111005A (z=0.013) detection of [O ii]. At this redshift, Hβ and [O ii] are located in the gap between the VIS and NIR arm and Hα is located in the The data presented here have previously been published in middle of the JH-bandgap and is therefore not detected. Michałowski et al.(2016). 2.4 ks spectroscopic integration of the host galaxy, obtained long after the burst had faded, con- Appendix B.37: GRB 111209A (z=0.677) tains bright emission lines filling the entire slit on top of a broad, underlying stellar continuum. We identify emission lines from These spectra have previously been used in Levan et al.(2013); [O ii], Hδ,Hγ,Hβ, [O iii], N ii,Hβ, [S ii], Ar iii, and [S iii], all at Greiner et al.(2015); Krühler et al.(2015); Kann et al.(2017) z = 0.013. Significant velocity structure of the lines across the and additionally formed the basis for GCN 12648 (Vreeswijk spatial direction of the slit indicates a large degree of coherent et al. 2011). The first epoch of spectroscopic observations was motion relative to the line-of-sight. initiated 17.7 hours after the BAT trigger and lasted for 4.8 ks. A very bright afterglow continuum is detected across the entire spectral coverage of X-shooter, with several absorption features Appendix B.32: GRB 111008A (z=4.989) imprinted. The absorption features are identified as Fe ii, Mg ii, Mg i, Ca H, and Ca K - all at a common redshift of z = 0.677. These data formed the basis of GCN 12431 (Wiersema et al. The second epoch, taken 20 days later, still contains a faint con- 2011) and are additionally published in Sparre et al.(2014). Ob- tinuum detected across all arms. The detection of nebular emis- servations of this GRB afterglow were initiated 7.71 hours after sion lines identified as [O ii], [O ii] and Hα at the same redshift, the BAT trigger and had a duration of 8.4 ks. A second observa- securely marks it at the redshift of this ultra-long GRB with ac- tional epoch started 20.1 hours after the GRB trigger and lasted companying GRB-SN. for 6.6 ks. The GRB afterglow continuum is well detected down to ∼760 nm, with several strong absorption features imprinted, all at a common z = 4.990. Lyα is clearly detected and we ad- Appendix B.38: GRB 111211A (z=0.478) ditionally detect lines identified as Si ii, Fe ii,C iv, Mg ii, Si ii*, These data formed the basis for GCN 12677 (Vergani et al. 2011) S ii*, O i*. An intervening DLA system is additionally detected and are also published in Krühler et al.(2015). Observations be- at z = 4.61 as seen from Lyα and Mg ii absorption. gan 31 hours after the AGILE trigger and consisted of 4 × 600 s. A bright GRB afterglow continuum is detected across the entire Appendix B.33: GRB 111107A (z=2.893) spectral coverage of X-shooter with absorption and emission features visible. We identify absorption features due to Fe ii, Mg ii, GCN 12542 (D’Avanzo et al. 2011b) is based on this spec- and Ca ii and emission lines from [O iii] and Hα, all at a common trum, but it is not published elsewhere. Spectroscopic integra- z = 0.478, which we suggest is the redshift of the GRB. Ad- tion started 5.26 hours after the Swift trigger and consists of ditionally detected in the GRB afterglow continuum are broad 4 × 1200 s integration in the UVB and VIS and 16 × 300 s in undulation, suggesting an accompanying SN. NIR, the observations ending in twilight. The GRB afterglow continuum is well detected across the arms with absorption lines Appendix B.39: GRB 111228A (z=0.716) from Lyα,C iv, Fe ii, and Mg ii, all at a consistent redshift of z = 2.893. Additionally an intervening Mg ii system is detected These data formed the basis of GCN 12770 (Schulze et al. 2011) at z = 1.998. From the Lyα absorption trough, we additionally and are also published in Krühler et al.(2015). Observations be- −2 infer log(NHI/cm ) = 21.0 ± 0.2. gan 15.9 hours after the BAT trigger and consist of 4 × 600 s.

Article number, page 33 of 43 A&A proofs: manuscript no. XSGRB_sample_arxiv

The GRB afterglow continuum is clearly detected in all the spec- 100. 4, emission lines belonging to Hβ, [O iii] and Hα are detected troscopic arms and superposed on the continuum are absorption at z = 0.350. The line belonging to Hα shows some extended features identiﬁed as due to Fe ii, Mn ii, Mg ii, Mg i, Ca H, and emission toward the afterglow continuum. The angular distance Ca K, all at z = 0.716. Supporting this redshift measurement as between the two sources correspond to a projected distance in the redshift of the GRB is the detection of nebular emission from the host plane of 6 kpc, posing a potential problem for the host [O iii]. redshift, unless the GRB occurred in a merging system. The extended emission in Hα, supports this interpretation. This burst is not a part of the statistical sample. Appendix B.40: GRB 120118B (z = 2.943)

The data presented here also formed the basis of GCN 14225 Appendix B.45: GRB 120327A (z = 2.813) (Malesani et al. 2013a), but are not published otherwise. This late-time observation of the host of GRB 120118B consists The data presented here also formed the basis of GCN 13134 of 3.6 ks exposures and contains emission lines belonging to (Krühler et al. 2012) and are published in D’Elia et al.(2014). [O ii] and [O iii] at z = 2.943, suggested to be redshift of the The observation consists of two visits, 2.13 hrs and 29.98 hrs host. after the burst, with an afterglow continuum visible in all arms for both visits. We detect absorption features from Ly-limit, Lyα, C ii/C ii∗, Si ii/Si ii∗, Al i, Fe ii and Mg ii at a consistent redshift, Appendix B.41: GRB 120119A (z = 1.728) z = 2.813. The data presented here have been examined by Japelj et al. (2015) and Zafar et al. (submitted), who both find a significant Appendix B.46: GRB 120404A (z = 2.876) amount of extinction, AV ≈ 1 mag, also supported by previous measurements (Morgan et al. 2014). Three epochs of ob- The data presented here have formed the basis of GCN 13227 servations have been obtained, the first two immediately after (D’Elia et al. 2012a), but are not published otherwise. 9.6 ks the burst, and the last one long after the afterglow had faded. integration, starting 15.7 hours after the Swift-trigger reveals a Starting 1.4 hours after the Swift trigger, the first epoch contains low-intensity afterglow continuum on which absorption from bright afterglow continuum. Rich in absorption features belong- Lyα is detected in two distinct regions at redshifts z = 2.876 ing to a multitude of ions at z = 1.728 is estimated for the host and z = 2.55. These absorption systems are confirmed by ionic with intervening systems at z = 1.476, z = 1.214, z = 0.662 absorption features at both of these redshifts. and z = 0.632. The second epoch, obtained 4.5 hours after the burst contains the fading afterglow. A third epoch is obtained > 1 year after the GRB in which emission lines from Hβ and Hα are Appendix B.47: GRB 120422A (z = 0.283) found at the redshift of the host, confirming the association of The data presented here also formed the basis of GCN 13257 the absorption line system and the host. We also detect C i in (Schulze et al. 2012) and are published in Schulze et al.(2014). absorption which indicates the presence of cold gas. Being a GRB-SN, this burst has been followed up multiple times. The data presented here only contain the first epoch in which the Appendix B.42: GRB 120211A (z = 2.346) afterglow is still visible and before the rise of SN2012bz. Start- ing 16.5 hours after the burst, 4.8 ks integration time captures The data presented here have been published in Krühler et al. both the host and the burst in emission. A blue afterglow con- (2015). Two observations of the host of GRB 120211A has been tinuum is detected at all wavelengths covered by X-shooter, on obtained, starting 2013.02.17, > 1 year after the burst has faded. which Mg ii absorption at z = 0.283 is found. Offset by 100. 75, A redshift for this object has been reported by Krühler et al. the host is clearly detected at a consistent redshift with a rich (2015) and the features seen by those authors are reproduced in emission line spectrum, the lines extending towards to burst. these reductions, confirming z = 2.346. Appendix B.48: GRB 120712A (z = 4.175) Appendix B.43: GRB 120224A (z = 1.10) The data presented here also formed the basis of GCN 13460 The data presented here have formed the basis of GCN 12991 (Xu et al. 2012), but are not published elsewhere. 4.8 ks inte- (Wiersema et al. 2012), and have also been published in Krühler gration time, starting 10.5 hours after the BAT trigger, shows a et al.(2015). Starting 19.8 hours after the GRB trigger, a total bright afterglow continuum starting at ∼ 472 nm, signifying the exposure time of 2.4 ks reveals a faint continuum, starting at ∼ onset of the Lyman alpha forest, for a GRB located at z = 4.175. 700 nm and extending all the way through 2500 nm. In the 2D- Absorption features from Lyα, Fe ii, Mg ii and Si ii are readily spectrum we detect a ∼ 2σ emission line which, if interpreted as detected at a consistent redshift. Hα, gives z = 1.10, supporting the photometric redshift (0.9 < z < 1.3) derived by Krühler et al.(2015). phot Appendix B.49: GRB 120714B (z = 0.398)

Appendix B.44: GRB 120311 (z = 0.350) The data presented here also formed the basis of GCN 13477 (Fynbo et al. 2012a) and are discussed in Klose et al. (submit- The data presented here have formed the basis of GCN 12991 ted), but are not published anywhere. Observations of this burst (Wiersema et al. 2012), but are not published otherwise. Starting started 7.8 hours after the GRB trigger, lasting 4.8 ks. A contin- just before twilight, 3.65 hours after the burst, a faint afterglow uum is visible across the entire spectral coverage of X-shooter, continuum is detected at all wavelengths. Due to the faintness of with both emission lines from [O ii], Hβ, [O iii] and Hα, as well the afterglow, no absorption features are discernible superposed as absorption from Mg ii detected at z = 0.398, securely setting on the continuum. Displaced from the afterglow continuum by it as the redshift of the GRB.

Article number, page 34 of 43 Selsing et al.: The X-shooter GRB afterglow legacy sample

Appendix B.50: GRB 120716A (z = 2.486) Lyβ, and Lyγ, visible in the Lyman alpha forest. Metal absorption lines arising from Fe ii, Ni ii, Si ii,S ii, Al ii, Al iii,C ii,O i, The data presented here also formed the basis of GCN 13494 C iv, and Zn ii are all detected along with the corresponding fine (D’Elia et al. 2012b), but are not published elsewhere. Despite structure lines (Fe ii∗, Si ii∗,O i∗,O i∗∗,C ii∗), securely anchoring observations starting 62 hours after the Swift trigger and lasting the redshift of the host. 3.6 ks, a bright afterglow is clearly seen, along with a plethora of absorption features. Absorption of Lyα-photons in the host leaves a broad trough, beyond which the Lyman alpha forest is Appendix B.55: GRB 120923A (z = 7.84) visible bluewards, all the way down to the Lyman limit. Metal absorption lines from C ii, Si ii,O i, Fe ii,C iv, Si iv, includ- This spectrum has previously been published in Tanvir et al. ing fine structure transitions identified as C ii∗, Si ii∗, Fe ii∗ and (2017). Starting 18.5 hours after the BAT trigger, the final spec- metastable Ni ii lines are all detected at z = 2.486 trum is based on 2 × 4 × 1200 s spectroscopic integration for a total of 160 minutes. Nothing is immediately visible in the 2D- spectra, however after severe binning in the dispersion direction, Appendix B.51: GRB 120722A (z = 0.959) a faint trace shows up in the NIR arm, suggesting a very high redshift. We here adopt the redshift (z = 7.84) suggested by Tan- The data presented here also formed the basis of GCN 13507 vir et al.(2017). (D’Elia et al. 2012c), but are not published elsewhere. In 4.8 ks integration time, starting 10 hours after the burst trigger, the simultaneous detection of absorption features belonging to Appendix B.56: GRB 121024A (z = 2.300) Mg ii and Fe ii superposed on a blue continuum, and emission lines from [O ii], Hγ,Hβ, [O iii] and Hα, all at z = 0.959, confi- The data presented here also formed the basis of GCN 13890 dently sets it as the redshift of the GRB. (Tanvir et al. 2012) and are published in Friis et al.(2015). Start- ing 1.8 hours after the Swift trigger, a bright afterglow continuum is visible across all arms. A broad absorption feature from Ly- Appendix B.52: GRB 120805A (z = 3.9) man alpha, along with narrow lines from C iv, Si ii, Si iv, Fe ii, ii ii ii∗ A separate reduction of this burst has been published in Krüh- S , and Al , as well as fine structure lines associated with Si ler et al.(2015), but is not used otherwise. Starting 9 days after are all detected at z = 2.300, securely setting it as the redshift of the burst trigger, this is a host observation and does not contain the GRB. any afterglow continuum. In 3.6 ks integration time, we detect a faint continuum visible from 450 nm and all the way through Appendix B.57: GRB 121027A (z = 1.773) 2100 nm, in contrast to what is found previously. The continuum from 450 - 600 nm is detected at very low significance. If The data presented here have formed the basis of GCN 13930 the drop at 450 nm is the Lyman limit, this fits with Lyman al- (Krühler et al. 2012), but are not published otherwise. Starting pha at ∼ 600 nm, giving z ∼ 3.9. The absence of nebular lines 69.6 hours after the GRB trigger, we detect the afterglow con- is due to [O ii] falling in a telluric absorption band and the rest tinuum at high significance in all arms with 8.4 ks integration, of the strong nebular lines being shifted out of the wavelength testifying to the brightness of this burst. The concurrent identi- coverage. fication of emission lines from [O iii] and absorption from C iv, Al ii, Al iii, Mg i, Mg ii, and Fe ii, tightly constrains the redshift of the burst to be z = 1.773. Appendix B.53: GRB 120815A (z = 2.358)

Not a part of the statistical sample, this burst also formed the Appendix B.58: GRB 121201A (z = 3.385) basis of GCN 13649 (Malesani et al. 2012) and is published in Krühler et al.(2013). Observations started 1.69 hours after the These data formed the basis for GCN 14035 (Sanchez-Ramirez BAT trigger and consist of 2.4 ks integration. A bright after- et al. 2012) and are additionally published in Krühler et al. glow continuum is detected across the entire spectral coverage (2015). These observations started 12.9 hours after the Swift trig- of X-shooter, with a multitude of absorption lines superposed. ger and consist of 4.8 ks spectroscopic integration under good Absorption features from the host at z = 2.358 include a DLA atmospheric conditions. The GRB afterglow continuum is well as well as metal absorption lines from N v,S ii, Si ii,O i,C iv, detected at all arms. A broad absorption trough due to Lyα is vis- Si iv, Fe ii, Al ii, Al iii, Mn ii, Mg ii, and Mg i. Also ﬁne structure ible at z = 3.385, which along with the detection of absorption lines from Ni ii and Fe ii are excited in the local environment of features identiﬁed as Si iv,C iv, Al ii, and Al iii, marks it as the the GRB. Additionally, this spectrum is one of the rare cases in redshift of the GRB. In the middle of the Lyα trough, we addi- which we detect lines from molecular hydrogen, H2. From the tionally detect Lyα emission. By modelling the Lyα absorption, −2 −2 Lyα-line we measure log(NHI/cm )=22.1 ± 0.10. Intervening we infer log(NHI/cm ) = 22.0 ± 0.3. systems are found at z = 1.539, z = 1.693, and z = 2.00. Appendix B.59: GRB 121229A (z = 2.707) Appendix B.54: GRB 120909A (z = 3.929) These data formed the basis for GCN 14120 (Fynbo et al. The data presented here have formed the basis of GCN 13730 2012b), but are not published elsewhere. Taken under poor see- (Hartoog et al. 2012), but are not published otherwise. Follow- ing conditions, a total of 4.8 ks spectroscopic integration starting up, started only 1.7 hours after the BAT trigger. This 1.2 ks ob- 2 hours after the Swift trigger yields a low signal-to-noise GRB servation captures a very bright afterglow continuum, starting at afterglow continuum is all arms. Binning the spectrum reveals 450 nm, signifying the onset of the Lyman limit for a system broad absorption troughs, which we identify as Lyβ and Lyα at at z = 3.929. Absorption from high-column density hydrogen z = 2.707. Additionally, an intervening system at z = 1.658 is leaves very prominent absorption features in the form of Lyα, detected from absorption features of Mg ii. From the absorption

Article number, page 35 of 43 A&A proofs: manuscript no. XSGRB_sample_arxiv

−2 trough due to Lyα, we infer log(NHI/cm ) = 21.7 ± 0.2. Due to the low signal-to-noise, the metal lines are weak, but the broad strong contamination in the slit, the background is slightly over absorption trough due to Lyα is detected. From the Lyα line we subtracted, causing the background to be negative in the center measure the redshift to be z = 2.780 and provide a measure of the −2 of the Lyα trough. neutral hydrogen column density, log(NHI/cm ) = 21.9 ± 0.3. The redshift is conﬁrmed by the presence of Fe ii absorption at a consistent redshift. Appendix B.60: GRB 130131B (z = 2.539)

These data formed the basis for GCN 14286 (Fynbo et al. 2013) Appendix B.65: GRB 130603B (z = 0.356) and are additionally published in Krühler et al.(2015). This is a late-time observation, taken long after the GRB afterglow had This burst is the first short GRB observed with a potential as- faded. In 7.2 ks spectroscopic integration, emission lines iden- sociated kilonova (Tanvir et al. 2013; Berger et al. 2013). GCN tified as [O ii] and [O iii] are detected at a common z = 2.539, 14757 (Xu et al. 2013a) was based on this spectrum and it is ad- which we suggest is the redshift of the GRB. ditionally published in de Ugarte Postigo et al.(2014b). Starting 8.2 hours after the Swift trigger, a total of 2.4 ks spectroscopic integration was obtained. Spectral continuum is clearly detected Appendix B.61: GRB 130408A (z = 3.758) across all arms from both host and afterglow and superposed The data presented here also formed the basis of GCN 14365 are both absorption (Ca H & K and Mg ii) and emission lines (Hjorth et al. 2013). The spectrum has not otherwise been pub- ([O ii], Hβ, [O iii], Hα, and [S ii]), all at a consistent redshift of lished previously. The observations consists of two 600 s spec- z = 0.356, which is the redshift of the GRB. tra taken 1.9 hours after the burst. We detect absorption features from a wide range of ions. We also detect intervening absorption Appendix B.66: GRB 130606A (z = 5.913) at z = 1.255 and z = 3.248. The data presented here also formed the basis of GCN 14816 (Xu et al. 2013e) and are published in Hartoog et al. Appendix B.62: GRB 130418A (z = 1.222) (2015). The observations consist of three 2 × 600 s vis- GCN 14390 (Krühler et al. 2013b) is based on this spectrum, its starting 7.1 hours after the burst at fairly high air- but it is not published elsewhere. Starting only 4.57 hours after mass. We detect absorption features from a wide range of the Swift trigger, 1.2 ks observations contain a bright GRB af- ions at z = 5.913 as well as intervening absorption at terglow continuum, visible across the entire spectral coverage of z = 2.3103, 2.5207, 3.4515, 4.4660, 4.5309, 4.5427, 4.6497 and X-shooter. Superposed on the afterglow continuum are absorp- 4.7244. tion features which we identify as C iv, Fe ii, and Mg ii, caused by an absorber at z = 1.217, and additional absorption from C iv at Appendix B.67: GRB 130612A (z = 2.006) z = 1.222. The two systems are offset by ∼ 1500 km s−1 and the proximity of the two absorption systems in velocity space, The spectral features of this spectrum have previously been re- suggests a possible association of the two systems with peculiar ported in GCN 14882 (Tanvir et al. 2013), but are not published velocity affecting the measured redshift. We adopt z = 1.222 as elsewhere. Starting only 1.1 hours after the Swift trigger, 2 × 600 the redshift of the GRB. Note, that this value is slightly different s spectroscopic integration captures a moderate signal-to-noise from the one reported in GCN 14390. afterglow continuum across the total spectral coverage of X- shooter. At a consistent redshift of z = 2.006, absorption from the metal resonance lines Fe ii, Mn ii, Mg ii, Mg i are identified. Appendix B.63: GRB 130427A (z = 0.340) Additionally, Lyα is visible as a broad absorption trough, from −2 This spectrum is also published in Xu et al.(2013) and Krüh- which we can infer log(NHI/cm ) = 22.2 ± 0.2, which is in the ler et al.(2015) and additionally has formed the basis for GCN upper end of the hydrogen column density distribution. The blue 14491 (Flores et al. 2013b). Starting 16.5 hours after the BAT part of the GRB continuum exhibits a downturn in the continuum trigger, these observations lasting 2 × 600 s contain a very bright level which could indicate the presence of a significant amount GRB afterglow continuum across the total spectral coverage of of dust along the line-of-sight. X-shooter. In absorption we identify features from the following metal resonance lines: Fe ii, Mn ii, Mg ii, Mg i, Ti ii and additional Appendix B.68: GRB 130615A (z = 2.9) line absorption from Ca ii and Na ID. Simultaneously, we find emission lines from Hα,Hβ, [O iii], [O ii]- all at common red- This spectrum has not previously been published. Starting only shift of z = 0.340, which is the redshift of the GRB. This is one 45 minutes after the BAT trigger, 2×600 s spectroscopic integra- of the most energetic GRBs observed, and its proximity along tion carried on into the beginning twilight. Observed at very high with the associated broad-lined Type Ic SN, 2013cq has caused airmass with a quickly varying background, a faint afterglow it to be one of the more well-studied GRBs (Maselli et al. 2014; trace is visible across all arms of X-shooter, down to 480 nm, Perley et al. 2014; Ackermann et al. 2014). which if interpreted as the break due to Lyα suggests z = 2.9. This supports the approximate redshift suggested in GCN 14898 (Elliott et al. 2013). Appendix B.64: GRB 130427B (z = 2.780)

This spectrum formed the basis of GCN 14493 (Flores et al. Appendix B.69: GRB 130701A (z = 1.155) 2013a), but are not published otherwise. A short, 2x600 s spectroscopic integration obtained before twilight, 20.6 hours after These data formed the basis for GCN 14956 (Xu et al. 2013b) the BAT trigger, captures a faint afterglow continuum visible and are additionally published in Krühler et al.(2015). Starting across the entire spectral coverage at low signal-to-noise. Due to 5.5 hours after the GRB trigger, 2 × 600 s reveals a bright con-

Article number, page 36 of 43 Selsing et al.: The X-shooter GRB afterglow legacy sample tinuum visible across the entire spectral coverage of X-shooter. afterglow continuum across the entire spectral coverage of X- Superposed are absorption features which we identify as due shooter. There are deviations from the continuum at both emis- to Fe ii, Mg ii, Mg i, and Ca ii - all at a consistent redshift of sion and absorption. We identify lines from Hβ, [O iii], and Hα in z = 1.155, which we take to be the redshift of the GRB. emission and Fe ii, and Mg ii in absorption. All lines are at a consistent z = 1.686, which is probably the redshift of the GRB. Absorption lines at shorter wavelengths are also detected, but at Appendix B.70: GRB 130925A (z = 0.347) low signiﬁcance due to an apparent downturn in the continuum This spectrum has already been used in GCN 15250 (Sudilovsky caused by the presence of dust local to the burst. et al. 2013) and is additionally published in Schady et al.(2015) and Krühler et al.(2015). Observations of this burst began with Appendix B.75: GRB 131117A (z = 4.042) X-shooter 3.5 hours after the Swift trigger. 6 ks spectroscopic integration captures a heavily dust obscured afterglow (AV = This spectrum has previously been used in GCN 15494 (Hartoog 5.9 ± 0.7; Greiner et al. 2014), with the spectrum primarily dom- et al. 2013), but is not published in a refereed paper. Starting only inated by host emission lines. All the nebular lines ([O ii], Hγ, 68 minutes after the BAT trigger, 4.8 ks spectroscopic integration Hβ, [O ii], Hα,N ii, [S ii]) are well detected at z = 0.347, as well secures afterglow continuum for this burst, which is measured to as those from Paδ, Paγ, and Paβ. We take this as the redshift be at z = 4.042. A moderate signal-to-noise GRB afterglow con- of the GRB. This spectrum is taken under ESO programme ID: tinuum is detected down to ∼ 610 nm, signifying the onset of the 091.A-0877(A) (PI: Schady). Lyα forest, with part of the forest also visible. Metal absorption lines from Si ii and Si iv are detected at a consistent redshift. Appendix B.71: GRB 131011A (z = 1.874) Appendix B.76: GRB 131231A (z = 0.642) These data formed the basis for GCN 15330 (Rau et al. 2013), but are not published elsewhere. Starting ∼1.5 days after the This spectrum has previously been published in Krühler et al. Fermi-GBM trigger, 4.5 ks spectroscopic integration captures a (2015) and addtionally forms the basis for GCN 15645 (Xu et al. modest signal-to-noise GRB afterglow continuum all the way 2014b). This spectrum, observed the following year (20.2 hours down to ∼320 nm. Imprinted on the continuum are absorption after the Swift trigger), consists of 4 × 600 s exposures. A high features, which we identify as due to Lyα, Fe ii, Mg ii, Mg i at signal-to-noise GRB afterglow continuum is detected all the way the same redshift, which we measure to be z = 1.874. From the through the X-shooter arms. We identify absorption features im- −2 broad absorption trough due to Lyα, we infer log(NHI/cm ) = printed in the continuum as caused by Fe ii, Mg ii, and Ca H & 22.0 ± 0.3. This spectrum is taken under ESO programme ID: K at a consistent z = 0.642. By subtracting oﬀ the bright after- 092.D-0056(A) (PI: Rau). glow continuum, we readily detect emission lines arising from [O ii], Hγ,Hβ, [O iii], and Hα in the GRB host galaxy. This GRB is the subject of a forthcoming work Kann et al. (in prep). Appendix B.72: GRB 131030A (z = 1.296)

These data has not been published before. Starting 3.4 hours af- Appendix B.77: GRB 140114A (z = 3.0) ter the Swift trigger, 6 × 600 s exposure were taken under good conditions, containing a bright GRB afterglow continuum across This spectrum has previously been published in Krühler et al. the entire spectral coverage of X-shooter. A myriad of absorp- (2015). This is a late time observation, taken long after the GRB tion features are superposed on the afterglow continuum which had faded. Despite a long integration time of 5.4 ks, no clear we identify as being caused by Si iv, Si ii,C iv, Al ii, Al iii, Zn ii, features stand out to clearly secure a redshift measurement. By Cr ii, Ni ii, Fe ii, Ni ii∗, and Fe ii∗ at z = 1.296. A very strong heavily binning the spectrum, a faint trace is visible down to 485 Mg ii-absorber is also detected, intervening the line-of-sight at nm, which if interpreted at the onset of Lyα signiﬁes z ∼ 3. We z = 1.164 with lines from Si ii,C iv, Al iii, Al iii, Fe ii, Mn ii, and adopt the redshift inferred in Krühler et al.(2015). many more. Appendix B.78: GRB 140213A (z = 1.208) Appendix B.73: GRB 131103A (z = 0.599) The data presented here also formed the basis of GCN 15831 This spectrum has already been used to form the basis for GCN (Schulze et al. 2014) and are additionally published in Krühler 15451 (Xu et al. 2013d) and is additionally published in Krühler et al.(2015). Starting 5.8 after the Swift trigger alert, 2 × 600 et al.(2015). Starting 5.8 hours after the BAT trigger, 4 × 600 s s spectroscopic integration contains a high signal-to-noise GRB exposure captures a modest signal-to-noise continuum across all afterglow continuum across the entirety of X-shooter. Imprinted arms. Imprinted on the continuum are absorption features identi- on the afterglow continuum are absorption features, which we fed as due to Fe ii and Mg ii as well as emission lines from [O ii], identify as metal resonance lines from C iv, Al ii, Al iii, Fe ii, Hδ,Hγ,Hβ, [O iii], Hα, and [N ii]. All the lines are measured at Mg ii, and Mg i. These lines are likely formed by metals in the a consistent redshift of z = 0.599, which we take as the redshift GRB host, which we measure to be at z = 1.208. of the GRB. Appendix B.79: GRB 140301A (z = 1.416) Appendix B.74: GRB 131105A (z = 1.686) These data formed the basis for GCN 15900 (Krühler et al. This spectrum has already been used in GCN 15450 (Xu et al. 2014b) and are additionally published in Krühler et al.(2015). 2013f) and is additionally published in Krühler et al.(2015). Spectroscopic follow-up began 9 hours after the BAT trigger Starting only 1.3 hours after the Swift trigger, a total of 4.8 ks and lasted for 12 × 600 s. A low signal-to-noise, spatially ex- spectroscopic integration contains a low signal-to-noise GRB tended continuum is visible across the entire spectral coverage

Article number, page 37 of 43 A&A proofs: manuscript no. XSGRB_sample_arxiv of X-shooter. The GRB afterglow continuum is visible at mod- GRB afterglow, which turns out to be at very high redshift. The erate signal-to-noise on top of the underlying host continuum. continuum is detected at moderate signal-to-noise in both the The host exhibits usual nebular emission from [O ii], Hβ, [O iii], VIS and NIR arms of X-shooter, heavily affected by absorption Hα,N ii, and [S ii] which puts it at z = 1.416. Supporting this as features. We identify the the lines belonging to Lyα, Si ii,C ii, the redshift of the GRB are absorption features from the Mg ii- C ii∗, Al ii, Al iii, Fe ii, and Mg ii at a consistent redshift of z = doublet in the GRB afterglow continuum. 4.233. From the shape of the Lyα absorption trough we measure −2 log(NHI/cm ) = 21.3 ± 0.3. Appendix B.80: GRB 140311A (z = 4.954) Appendix B.85: GRB 140622A (z = 0.959) This spectrum has not been published previously. Starting 32.5 hours after the GRB trigger onboard Swift, this observation The characteristics of this short (T90 = 0.13 ± 0.04 s; Lien et al. lasted 14 × 600 s for a total of 8.4 ks. Some loss occurred dur- 2016) GRB spectrum has been published in GCN 16437 (Har- ing the observations, which reduced integration time in the UVB toog et al. 2014), but does not appear in the refereed literature. and VIS arm slightly. The GRB afterglow continuum is clearly Spectroscopic observations began with X-shooter only 34 min- visible in the VIS and NIR arm of X-shooter. The continuum is utes after the BAT trigger and lasted for 2 × 600 s. In the spec- very rich in absorption features, with at least the following lines trum, we detect continuum across all three arms of X-shooter. In identified: Lyγ, Lyβ, Lyα, Si ii, Si iv,C iv, Al ii, Al iii, Fe ii, Mg ii, the UVB arm, the continuum is only visible after heavily bin- and Mg i. All of these lines are at z = 4.954, which we take as ning the spectrum in the dispersion direction. It is unclear how the redshift of the GRB. This spectrum is taken under ESO pro- much of the continuum is from the host galaxy and how much gramme ID: 092.D-0633(E) (PI: Greiner). is from the potential GRB afterglow. Superposed on the continuum, emission lines from [O ii], Hβ, [O iii], Hα are all detected at z = 0.959. We use this as the likely redshift of the GRB. Appendix B.81: GRB 140430A (z = 1.601) A separate reduction of this burst has been published in Krühler Appendix B.86: GRB 141028A (z = 2.332) et al.(2015), and additionally the spectrum has also been used in GCN 16194 (Krühler et al. 2014a). Observations for this burst This spectrum has already been used in GCN 16983 (Xu et al. began in twilight, 2.5 hours after the BAT trigger and lasted for 2014a), but is not published elsewhere. Starting 15.4 hours after 2×600 s. The spectrum contains a moderate signal-to-noise GRB the Swift trigger and lasting for a total of 2.4 ks, these obser- afterglow continuum all the way through the spectroscopic arms vations captures a bright GRB afterglow. Very rich in absorp- of X-shooter. We identify absorption features in the afterglow tion, the continuum is detected across the entire spectral cover- continuum from Si ii,C iv, Al ii, Fe ii, and Mg ii, and emission age of X-shooter, except the blue half of the UVB arm, where lines from [O ii] and [O iii] - all at z = 1.601 which is likely the the Lyα-forest absorbs the continuum. The redshift of the host redshift of the GRB. of the GRB is measured based in the detection of features from Lyα, Si ii,C iv,C ii, Fe ii, and Mg ii- all at a consistent z = 2.332. From the Lyα trough we infer log(N /cm−2) = 20.6±0.15. Two Appendix B.82: GRB 140506A (z = 0.889) HI intervening system at z = 1.823 and z = 2.09 are also found in The data presented here has formed the basis of GCN 16217 the spectrum based on the detection of C iv. (Fynbo et al. 2014) and is published in Fynbo et al.(2014), Krühler et al.(2015), and Heintz et al.(2017). The observa- Appendix B.87: GRB 141031A (z = na) tions consists of 4 × 600 s at 8.8 and 33 hours after the burst. We detect absorption features from a wide range of ions, to- This spectrum is an attempt at a late-time host redshift measure- gether with molecular absorption from CH+, all at z = 0.889. ment for GRB 141031A. The spectrum is taken long after the The optical/near-infrared afterglow reveals an unusual steep ex- burst has faded and consists of 4 × 600 s spectroscopic integra- tinction curve which is found to be caused by dust very close to tion. In the 2D spectrum are two sources, which are both offset the burst. from the targeted host position and are thus likely foreground objects. The spectrum does not contain anything that can be used Appendix B.83: GRB 140515A (z = 6.327) to measure a redshift from. This spectrum has previously been used in Melandri et al. (2015). Starting 15.5 hours after the Swift trigger, 8×600 s spec- Appendix B.88: GRB 141109A (z = 2.993) troscopic integration captures this very high redshift GRB af- The gross content of the spectrum has been issued in GCN 17040 terglow. The Gunn-Peterson trough is visible against the GRB (Xu et al. 2014c), but it is not used otherwise. Starting only afterglow continuum, starting at 890 nm, which along with ab- 1.9 hours after the BAT trigger, these 4 × 600 s spectra con- sorption lines from the Mg ii doublet securely sets the redshift of tain a high signal-to-noise GRB afterglow. The afterglow con- this GRB at z = 6.327. From the red wing of the Lyα-profile we tinuum is readily detected all the way down to 370 nm, although −2 measure log(NHI/cm ) = 19.0 ± 0.5, which is very low com- the bluest part is affected by the Lyα-forest. A broad absorp- pared to the measured distribution of NHI. tion trough from neutral hydrogen is clearly visible, with additional, narrower absorption features identified as due to Si ii, ∗ ∗ ∗ ∗ ∗ Appendix B.84: GRB 140614A (z = 4.233) Si ii ,C ii,C ii , Si iv,C iv, Fe ii, Fe ii ,O i , and N ii . These are all detected at a consistent z = 2.993, which due to the detec- This spectrum forms the basis for GCN 16401 (Krühler et al. tion of the locally excited fine-structure lines securely sets this 2014c), but is not published elsewhere. 4 × 600 s spectroscopic as the redshift of the GRB. Additionally, two Mg ii absorption integration, starting 3.8 hours after the BAT trigger catches this systems are detected at z = 1.67 and z = 2.5. From the broad

Article number, page 38 of 43 Selsing et al.: The X-shooter GRB afterglow legacy sample

Lyα-absorption, we measure a neutral hydrogen column density Appendix B.93: GRB 150428A (z = na) −2 of log(NHI/cm ) = 22.1 ± 0.1. This spectrum is empty, but is included here for completeness. It is not published anywhere. 4 × 600 s spectroscopic observation, Appendix B.89: GRB 150206A (z = 2.087) starting 3.7 hours after the trigger does not reveal anything con- clusive. The host association is additionally ambiguous for this, This spectrum has previously been used to form the basis for likely extincted (GCN 17767; Knust et al. 2015), GRB. GCN 17420 (Krühler et al. 2015), but it is not published in the refereed literature. Due to the short visibility of this burst, it does not enter into the statistical sample. Beginning in twi- Appendix B.94: GRB 150514A (z = 0.807) light, 4×600 s integration time contains a weak/moderate signal- to-noise GRB afterglow continuum throughout the X-shooter This spectrum has already been used in GCN 17822 (de Ugarte spectral coverage, down to ∼ 375 nm. In the continuum are ab- Postigo et al. 2015c). Spectroscopic observations began 28.4 sorption features, which we identify as due to metal lines from hours post trigger and consist of a 4 × 600 s nodding sequence. Zn ii, Fe ii, and Mg ii. Additionally, Lyα-absorption is seen at The GRB afterglow continuum is detected at moderate signal- the end of the trace. The spectral position of the lines means to-noise across the entire spectral coverage of X-shooter. Nar- that this GRB is at z = 2.087. From the Lyα trough, we infer row absorption features are imprinted in the continuum where −2 log(NHI/cm ) = 21.7 ± 0.4. The afterglow continuum appears we identify features from multiple Fe ii transitions as well as the depressed in the blue end of the spectrum, suggestion dust ex- Mg ii-doublet and the Mg iλ2852.96 resonance line. These lines tinction in the host. are all found at a position matching z = 0.807, suggesting it is the redshift of the GRB.

Appendix B.90: GRB 150301B (z = 1.517) Appendix B.95: GRB 150518A (z = 0.256) This spectrum has already been used in GCN 17523 (de Ugarte Postigo et al. 2015a), but it is not used elsewhere. Starting 5.1 This gross content of this data has previously been issued in hours after the BAT trigger, this spectrum is based on 6 × 600 s GCN 17832 (Xu et al. 2015a), but is not published. Starting > 1 spectroscopic integration. The GRB afterglow continuum is well day after the GRB trigger, 4 × 600 s spectroscopic integration detected across the entire spectral coverage of X-shooter at mod- securely allows us to measure the redshift of the host. A contin- erate signal-to-noise. Imprinted on the continuum are absorption uum is detected all the way through the spectral coverage with features which we identify as being caused by Si ii,C iv, Al ii, multiple emission lines superposed. It is not clear to what degree Fe ii, Mg ii, and Mg i- all at a similar redshift of z = 1.517, which the GRB afterglow contributes to the continuum. We identify the is likely the redshift of the GRB. emission lines as [O ii], Hβ, [O iii], Hα, [N ii], Hα, and [S ii] - all at z = 0.256. Due to the spatial proximity of this burst, it is a candidate for SN follow-up and indeed there are indications Appendix B.91: GRB 150403A (z = 2.057) of re-brightening at the burst position (GCN 17903; Pozanenko et al. 2015). This spectrum has previously been used to form the basis for GCN 17672 (Pugliese et al. 2015). 4 × 600 s spectroscopic integration, starting 10 hours after the GRB trigger, captures a bright Appendix B.96: GRB 150616A (z = 1.188) GRB afterglow. A broad absorption trough centered at ∼ 370 nm due to Lyα signals the redshift of the GRB, which is reﬁned to These host observations are taken long after the burst had faded z = 2.057 based on the additional detection of metal absorp- and have not been published before. They are included here for tion lines. We readily identify features associated with S ii, Si iv, completeness. This bursts is excluded from the sample because O i, Si ii, Si ii∗,C ii,C ii∗,C iv, Al ii, Fe ii, Fe ii∗, Mn i, Mg ii, and an observing constraint delayed the Swift slew, causing the XRT Mg i in the host of the GRB. An intervening C iv absorber is ad- observations to begin 16 minutes post trigger. No continuum is ditionally detected at z = 1.76. From the absorption trough due detected, but emission lines from [O ii], [O iii], and Hα are all to Lyα, we infer the amount of neutral hydrogen in the host along detected at z = 1.188, setting it as the redshift of the GRB host. −2 the line of sight to be log(NHI/cm ) = 21.8 ± 0.2.

Appendix B.97: GRB 150727A (z = 0.313) Appendix B.92: GRB 150423A (z = 1.394) The overall content of these observations have previously been The gross content of these observations have previously been reported in GCN 18080 (Tanvir et al. 2015b), but is not pub- presented in GCN 17755 (Malesani et al. 2015a), but is not published. Starting 5 hours after the Swift GRB trigger, these obser- lished as part of any refereed paper. This bona-ﬁde short GRB vations consist of a combined 2×1200 s and 2×600 s integration. (T90 = is 0.22 ± 0.03 s; Lien et al. 2016) was observed in RRM A blue continuum is detected across the entire spectral cover- mode and spectroscopic integration started after only 22 min- age of X-shooter, suggesting a signiﬁcant contribution from the utes. A series of stare mode observations, increasing in exposure GRB afterglow. Superposed on the continuum are emission lines time, and ending with a nodding sequence, totaling ∼ 5000 s, which we identify as Hβ, [O iii], and Hα with a measured redshift are combined to form this spectrum. A faint, almost featureless of z = 0.313. Supporting this redshift is the tentative detection of continuum is detected at low signal-to-noise all the way to the the Mg ii absorption doublet in the afterglow continuum. The rel- bluest part of the spectrum. An absorption doublet is detected in ative strength of the lines suggest that the line-forming region is the VIS arm against the GRB afterglow continuum, which we dust obscured, contrary to the story told by the blue afterglow identify as Mg ii at z = 1.394. continuum.

Article number, page 39 of 43 A&A proofs: manuscript no. XSGRB_sample_arxiv

Appendix B.98: GRB 150821A (z = 0.755) range of absorption lines which we identify as due to Si ii, Si ii∗, O i,C ii,C ii∗,C iv, Al ii, Al ii, Fe ii, and Fe ii∗. These lines are all The gross content of this spectacular spectrum has already been detected at z = 4.063 and the presence of fine-structure lines at a issued in GCN 18187 (D’Elia et al. 2015a), but the spectrum consistent redshift securely sets this as the redshift of the GRB. has not been published. Starting just 12.4 minutes after the on- GRB 151027B was the 1000th GRB detected by Swift. The af- × board trigger of Swift, a total of 4 600 s spectroscopic inte- terglow of this burst is the subject of a forthcoming publication gration was obtained. The observations began just before dawn (Greiner et al. 2018). and the last two exposures are heavily affected by the brightening sky. A bright, high signal-to-noise GRB afterglow continuum is detected all across the the spectral window and imprinted on Appendix B.103: GRB 151029A (z = 1.423) this are a myriad of absorption features. We identify individual lines from transitions in Al iii, Cr ii, Zn ii, Ni ii∗, Fe ii, Fe ii∗, Sc ii, The gross content of these observations have previously been de- Mn ii, Mg ii, Mg i, Ti ii, and Ca ii - all at a consistent redshift of scribed in GCN 18524 (Tanvir et al. 2015a). Starting only 1 hour z = 0.755. The detection of fine-structure lines, excited local to after the BAT trigger, 2 × 600 s observations, extending into the the burst, clearly marks this as the redshift of the GRB. morning twilight, were obtained of this GRB. The brightening sky affects the signal-to-noise of these spectra - especially in Appendix B.99: GRB 150910A (z = 1.359) the blue end. The GRB afterglow is well detected over the entire spectral coverage of X-shooter, with spectral features from This spectrum has not previously been published. Starting ∼ 20 Fe ii and Mg ii imprinted by an absorber at z = 1.423. The ab- hours after the BAT trigger, these observations were stopped af- sorber is likely the host of the GRB and we therefore consider ter 2 × 600s, when it became apparent that only a modest signal- z = 1.423 the redshift of the GRB. to-noise was obtainable, and the redshift had already been published in GCN 18273 (Zheng et al. 2015) and GCN 18274 (de Ugarte Postigo et al. 2015b). The continuum is detected at low Appendix B.104: GRB 151031A (z = 1.167) signal-to-noise all across X-shooter and Mg ii is detected at the suggested redshift, which we take as the redshift of the GRB. The spectral features in this spectrum has previously been reported in GCN 18540 (Malesani et al. 2015b). Observed in RRM Appendix B.100: GRB 150915A (z = 1.968) mode, these observations were initiated only 19 minutes after the BAT trigger. As per the standard RRM scheme, a series of stare The gross content of these spectra have already been issued in observation, increasing in en exposure time and ending with a GCN 18318 (D’Elia et al. 2015b). Observations began 3.3 hours regular 4 × 600 s nodding sequence, were acquired for this burst. after the Swift trigger and lasted for 4 × 1200 s. A moderate A moderate signal-to-noise continuum is detected all the way signal-to-noise afterglow continuum is detected across the entire through the X-shooter spectral coverage. In the GRB afterglow spectral coverage of X-shooter, where the bluest part is affected continuum, we identify absorption features as Fe ii, Mg ii, Mg i, by Lyα-forest absorption. Imprinted on the afterglow continuum and Ca ii at a consistent redshift of z = 1.167. Supporting this as are a wealth of both emission and absorption features, all caused the redshift of the GRB host galaxy are the simultaneous detec- in a system for which we measure a redshift of z = 1.968. We tion of nebular emission lines from [O ii], Hβ, and [O iii]. identify in emission [O ii], Hβ, [O iii], and Hα and absorption due to Lyα,C iv, Al ii, Si ii, Fe ii, and Mg ii. We additionally identify ∗ ∗ fine-structure absorption lines from Si ii and Fe ii , at a similar Appendix B.105: GRB 160117B (z = 0.870) redshift, which unequivocally sets the suggested redshift as the redshift of the GRB. The Lyα-line is affected by an atmospheric A preliminary analysis of these spectra have previously been pre- transmission drop and Lyα emission in the trough, making the sented in GCN 18886 (de Ugarte Postigo et al. 2016d). Starting measurement of the neutral hydrogen column density difficult. 13.5 hours after the Swift trigger, these spectra are formed on the −2 We measure log(NHI/cm ) 21.2 ± 0.3. basis of 4×1200 s spectroscopic integration. The GRB afterglow is clearly detected at high signal-to-noise across the entire spec- Appendix B.101: GRB 151021A (z = 2.330) tral coverage of X-shooter. Superposed on the continuum are absorption and emission features, which are identified as Fe ii, The data presented here also formed the basis of GCN 18426 Mg ii, Mg i in absorption, and [O ii], Hβ, and [O iii]. The spectral (Pugliese et al. 2016) and are not published elsewhere. The ob- positions of these lines all correspond to a redshift of z = 0.870, servation was carried out in RRM starting 44 minutes after the which most likely is the redshift of the GRB. GRB trigger. We detect absorption features from a wide range of ions at z = 2.330 as well as intervening absorption at z = 1.49. Appendix B.106: GRB 160203A (z = 3.518) Appendix B.102: GRB 151027B (z = 4.063) The data presented here also formed the basis of GCN 18982 The content of these spectra has previously been issued in GCN (Pugliese et al. 2016) and are not published elsewhere. The ob- 18506 (Xu et al. 2015b), but is not published. Beginning 5 hours servation was carried out in RRM starting 18 minutes after the after the BAT trigger, 4 spectroscopic integrations obtained in a GRB trigger. We detect absorption features from a wide range of nodding sequence, each lasting for 600s, securely allows us to ions at z = 3.518 as well as intervening absorption at z = 2.203. measure a redshift for this GRB. Beginning at ∼ 470 nm, the From the very high signal-to-noise spectrum we are able to mea- −2 Lyα-forest is clearly detected at high signal-to-noise leading up sure log(NHI/cm ) = 21.75 ± 0.10, based on the shape of the to the broad Lyα absorption trough and the onset of a bright GRB Lyα profile. This burst is the subject of a forthcoming publica- afterglow continuum. Imprinted in the afterglow continuum are a tion (Pugliese et al. in prep.).

Article number, page 40 of 43 Selsing et al.: The X-shooter GRB afterglow legacy sample

Appendix B.107: GRB 160228A (z = 1.640) the spectral coverage of X-shooter, but no clear absorption lines are discernible. It is therefore unclear how much of the contin- The gross content of these spectra have previously been pre- uum is due to the host and how much is due to the GRB. Su- sented in GCN 19186 (Krühler et al. 2016a). Taken long after the perposed on the continuum are two sets of emission lines, oﬀset × burst had faded, the target is the likely host galaxy. In 4 1200 spatially, but at very similar redshifts. We identify lines from s spectroscopic integration we detected emission lines from the [O ii], Hβ, [O iii], Hα, [N ii], and [S ii] in both systems – all at iii host. We ﬁnd [O ] and Hα at a consistent z = 1.640. We there- z = 0.555. As per GCN 19350 (Tanvir et al. 2016c), the GRB fore suggest that this is the redshift of the GRB. probably occurred in an interacting pair of galaxies.

Appendix B.108: GRB 160303A (z = na) Appendix B.112: GRB 160625B (z = 1.406) This is potentially a short burst with extended emission (GCN 19148; Ukwatta et al. 2016). The content of the spectra have GRB 160625B is the first GRB for which a signicant linear po- previously been reported in GCN 19154 (de Ugarte Postigo et al. larization was measured (Troja et al. 2017). Because this burst 2016c). Observations began 19.1 hours after the BAT trigger and was detected by Fermi, it is therefore not a part of the statis- consist of 4 × 1200 s spectroscopic integration. A very faint con- tical sample. The gross content of these spectra has previously tinuum is detected in the UVB and VIS arm after binning in the been presented in GCN 19600 (Xu et al. 2016b). Starting 30 dispersion direction, but no clear redshift can be inferred from hours after the trigger, these spectra are obtained in 4 × 600 the continuum shape. No emission lines are readily apparent, s. The GRB afterglow continuum is detected at high signal-to- suggesting little or no star formation in the host. noise throughout the spectral coverage of X-shooter. Imprinted on the afterglow continuum are a multitude of absorption features, from which we identify the responsible elements as Si ii, Appendix B.109: GRB 160314A (z = 0.726) O i, Si iv,C iv, Al ii, Al iii, Fe ii, Zn ii, Mg ii, and Mg i at a consistent redshift of z = 1.406. An additional intervening Mg ii ab- The gross content of these spectra have previously been issued sorber is detected at z = 1.319. in GCN 19192 (D’Elia et al. 2016). Spectroscopic integration began 13 hours post GRB trigger and lasted for 4 × 1200 s. In the data we see a moderate signal-to-noise continuum across the Appendix B.113: GRB 160804A (z = 0.736) entire spectral coverage of X-shooter. Clearly visible on top of the continuum are emission lines from [O ii], Hδ,Hγ,Hβ, [O iii], The data presented here also formed the basis of GCN 19773 Hα, [N ii], and [S ii] at a consistent redshift of z = 0.726. It is (Xu et al. 2016a), and are published in Heintz et al.(2018). Ob- unclear how much of the continuum is due to the host and how servations started 22.37 hours after the BAT trigger and lasted much is from the GRB afterglow. At the same redshift, we see for 2.4ks. The afterglow continuum is detected across the entire the tentative detection of Fe ii-absorption, supporting the redshift spectral coverage of X-shooter and absorption lines from Mg i, measurement and suggesting a significant afterglow contribution Mg ii, Fe ii and Al ii are found at z = 0.736. At the same redshift, to the continuum. emission lines from [O ii], [O iii], Hα,Hβ,Hγ, [N ii], [S ii], and [S iii] are found. A second epoch, lasting 3.6ks, is obtained after the afterglow had faded, confirming the emission line detections. Appendix B.110: GRB 160410A (z = 1.717) The host galaxy is found to have a roughly solar metallicity and This short GRB with extended emission (GCN 19276; Sakamoto is among the most luminous GRB hosts at z < 1. et al. 2016) was observed in RRM mode, with observations starting after only 7.7 minutes. The gross content has already been Appendix B.114: GRB 161001A (z = 0.891) presented in GCN 19274 (Selsing et al. 2016b). These observations have the shortest delay between trigger time and start of ob- The content and scope of these spectra has previously been is- servation. Spectroscopic integration began after 8.4 minutes and sued in GCN 19971 (Krühler et al. 2016b). 4 × 600 s spectro- lasted for 3 × 600 s. The last image in the nodding sequence was scopic integration was obtained of this potentially short GRB aborted because the object was setting and the telescope reached (GCN 19974; Markwardt et al. 2016), starting 6.1 hr after the the hardware limit. A blue, moderate signal-to-noise continuum GRB trigger. Close to the XRT error circle(Osborne et al. 2016) is detected across the entire spectral coverage of X-shooter with there is a bright, point-like source, also visible in archival im- multiple absorption features imprinted on it. At z = 1.717 we ages. Our spectrum reveals that this is an M-type star, based on detected the broad absorption trough due to Lyα and we addi- the presence of TiO bands, and is thus unrelated to the GRB. Two tionally find Al ii and Fe ii at a consistent redshift, setting it as more objects lie within or in the near proximity of the X-ray er- the redshift of the GRB. Additionally, two intervening C iv ab- ror circle, as visible in our images from the acquisition camera. sorption systems are detected at z = 1.444 and z = 1.581. From One is blended with the above mentioned star, and is serendip- the Lyα absorption trough, we measure a neutral hydrogen col- itously covered by the X-shooter slit, revealing emission lines −2 umn density of log(NHI/cm )= 21.2±0.2. This is the only short from [O ii], Hβ, and Hα at a consistent redshift of z = 0.891. GRB for which we have a measurement of NHI and is the subject The second source, first noted by Chen et al.(2016), is extended of an forthcoming paper (Selsing et al. in prep). and lies just outside the XRT error circle. Based on GROND imaging, Chen et al.(2016) report moderate variability, and sug- Appendix B.111: GRB 160425A (z = 0.555) gest this object as an afterglow candidate. Following our origi- nal imaging, we thus secured further late-time photometry with These spectra were already used to form the basis for GCN the X-shooter acquisition camera, starting on 2016 Oct 5.26 UT 19350 (Tanvir et al. 2016c), but are not published anywhere else. (4.22 days after the GRB), securing 3 × 120 s in each of the g Spectroscopic integration started after 7.2 hours and lasted for and r filters. No variability is apparent in either bands. In offline 4 × 1200 s. A low signal-to-noise continuum is detected across analysis of the GROND data (T. Krühler, priv. comm.), the evi-

Article number, page 41 of 43 A&A proofs: manuscript no. XSGRB_sample_arxiv

−2 dence for fading is less significant than originally reported. The infer log(NHI/cm )= 20.96±0.05. This spectrum will be further situation is thus left ambiguous, but due to the consistency of the analyzed in de Ugarte Postigo et al. (in prep). z = 0.891 object with the XRT localization, we consider this as the most likely redshift of GRB 161001A. Appendix B.118: GRB 161117A (z = 1.549) The majority of the spectral features of these spectra have pre- Appendix B.115: GRB 161007A (z = na) viously been reported in GCN 20180 (Malesani et al. 2016). ∼ These data have not been published elsewhere. Observations for Starting only 40 minutes after the BAT trigger, these observa- × GRB 161007A started 323 hr after the burst trigger and were tions consist of 4 600 s, taken under relatively poor observing centered on the potential host galaxy (see GCN 20014 (de Ugarte conditions. The GRB afterglow continuum is detected through- Postigo et al. 2016b) and GCN 20020 (Heintz et al. 2016)). out most of the X-shooter spectral coverage at low-to-moderate 4 × 600 s of observations were obtained revealing a faint con- signal-to-noise, however the continuum level drops below the ∼ tinuum rising abruptly above the noise at ∼ 685 nm and ex- noise in the bluest part of the UVB arm at 380 nm. Imprinted tending up to 2100 nm. A very low significance continuum is in the continuum are absorption features, which we identify as ii ii i detected at shorter wavelengths, down to ∼ 600 nm. No signifi- due to Fe , Mg , and Mg in an absorption system at z = 1.549. cant emission features could be identified. Two possibilities can We consider this the likely redshift of the GRB. explain the observed continuum break, but neither are fully sat- isfactory. The first option is the Lyman alpha cutoff at z ∼ 4.6; Appendix B.119: GRB 161219B (z = 0.148) at this redshift, the usually detected strong nebular lines would be shifted out of the covered wavelength range. The host would The core signatures in these spectra have already been issued in be however exceptionally bright, with an absolute magnitude at GCN 20321 (Tanvir et al. 2016b), and are additionally used in ≈ 150 nm M = −23.8 AB, that is about 2.8 mag brighter than both Ashall et al.(2017) and Cano et al.(2017). This GRB is M∗ at that redshift (e.g. Bouwens et al. 2015). The inferred X-ray very close-by and is associated with the SN Ic-BL SN2016jca hydrogen-equivalent column density would also be a record-high (GCN 20342; de Ugarte Postigo et al. 2016a). The observations 4 × 1023 cm−2 (see e.g. Campana et al. 2012). Alternatively, the consist of 4 × 600 s spectroscopic integration, starting 1.5 days continuum discontinuity could be the 400 nm break at z = 0.71; after the GRB trigger. The GRB afterglow continuum is detected such feature is however prominent only in early-type galaxies, at high signal-to-noise across the entire spectral coverage of X- which are unlikely hosts of long-duration GRBs (but see Rossi shooter. Superposed on the continuum are absorption features, et al. 2014). Given the lack of optical variability, it cannot be which we identify as being formed by the Mg ii, Mg i, and the excluded that the object within the XRT error circle is a chance Ca H & K-lines at z = 0.148. Underneath the glaring afterglow association, as the probability of a random alignment is small continuum we additionally detect nebular emission from [O ii], but non negligible (of the order of 1%). In conclusion, we can- Hβ, [O iii], Hα, and [S ii] at a consistent redshift, which is the not provide a secure redshift for GRB 161007A. likely redshift of the GRB.

Appendix B.116: GRB 161014A (z =2.823) Appendix B.120: GRB 170113A (z=1.968) These gross content of these data have previously been issued The data presented here also formed the basis of GCN 20061 in GCN 20458 (Xu et al. 2017). Starting 15.2 hours after the (Selsing et al. 2016a), but are not published elsewhere. Start- BAT trigger, these observations consist of 4 × 1200 s spectro- ing 11.6 hours after the GRB trigger, 4.8 ks of integration time scopic integration. The GRB afterglow continuum is detected at captures the afterglow continuum across all three spectroscopic moderate signal-to-noise across across the entire spectral cover- arms. A broad absorption trough due to Lyman alpha is visible, age of X-shooter. In the bluest end of the spectral coverage, the ii ii ii iv along with metal absorption features from Mg , Si ,C ,C , background noise increases, thus lowering the continuum below ii iii ii Al , Al , and Fe , all at z = 2.823. Similar to GRB 140506 detection. Superposed on the afterglow continuum are both ab- (Fynbo et al. 2014; Heintz et al. 2018), a break in the contin- sorption lines and emission lines. In absorption we identify lines uum shape is tentatively detected bluewards of 600 nm, possible from Si ii, Fe ii, and Mg ii, and in emission we detect [O ii], Hβ, signifying some anomalous form of extinction. and [O iii]. All the lines are detected at a consistent redshift of z = 1.968, which we take as the redshift of the GRB. Appendix B.117: GRB 161023A (z = 2.710) Appendix B.121: GRB 170202A (z=3.645) The gross content of these spectra have previously been reported in GCN 20104 (Tanvir et al. 2016a). This extremely bright INTE- GRB 170202A is the last GRB afterglow observed in this work. GRAL burst is not a part of the statistical sample. Starting 3 hours It concludes the follow-up eﬀort that has lasted for 8 years. The after the trigger, the spectra are based on 2 × 600 s. The GRB majority of the spectral features in this spectrum have previ- afterglow continuum is detected with the highest signal-to-noise ously been reported in GCN 20589 (Palmerio et al. 2017). Spec- of all the spectra presented here. A myriad of absorption features troscopy began 9.7 hours after the BAT trigger, and lasted for are visible against the continuum, many of which would not be 4 × 600 s. The GRB afterglow is well detected at moderate- visible against lower signal-to-noise spectra. We identify at least to-high signal-to-noise throughout the spectral coverage of X- 15 intervening absorbers between z = 1.243 and z = 2.710. In shooter, down to 430 nm. The lowest wavelength continuum the host system we clearly see absorption from Lyβ, Lyα,S ii, from 430 nm to 560 nm is absorbed by the Lyα forest, with Si ii,O i, Si iv,C iv, Al ii, Al iii, Fe ii, Mg ii, and Mg i. Several the broad Lyα absorption trough centered at 565 nm. In the ﬁne-structure transitions are additionally seen, clearly associat- continuum, we additionally identify metal absorption features ing the z = 2.710 system with the GRB. From the Lyα-line we from Si ii,C iv, Si iv, Fe ii, and Mg ii, with the additional de-

Article number, page 42 of 43 Selsing et al.: The X-shooter GRB afterglow legacy sample tection of Si ii∗ ﬁne-structure lines. With an additional detection of [O iii] and all lines detected at z = 3.645, this unequivocally marks this as the redshift of the GRB. There is an additional C iv absorber at z = 3.077. From the Lyα absorption we measure −2 log(NHI/cm )= 21.55 ± 0.10.

Article number, page 43 of 43