Search

A dedicated gamma-ray burst (GRB) afterglow observing program was performed between 2007 and 2016 with GROND, a seven-channel optical and near-infrared imager at the 2.2m telescope of the Max-Planck Society at ESO/La Silla. In this first of a series of papers, we describe the GRB observing plan, providing first readings of all so far unpublished GRB afterglow measurements and some observing statistics. In total, we observed 514 GRBs with GROND, including 434 Swift-detected GRBs, representing 81\% of the observable Swift sample. For GROND-observations within 30 min of the GRB trigger, the optical/NIR afterglow detection rate is 81\% for long- and 57\% for short-duration GRBs. We report the discovery of ten new GRB afterglows plus one candidate, along with redshift estimates (partly improved) for four GRBs and new host detections for seven GRBs. We identify the (already known) afterglow of GRB 140209A as the sixth GRB exhibiting a 2175 Angstroem dust feature. As a side result, we identified two blazars, with one at a redshift of z=3.8 (in the GRB 131209A field)., Comment: 20 pages, 31 figures, accepted for publication in Astronomy & Astrophysics

Core-collapse supernovae are explosions of massive stars at the end of their evolution. They are responsible for metal production and for halting star formation, having a significant impact on galaxy evolution. The details of these processes depend on the nature of supernova progenitors, but it is unclear if Type Ic supernovae (without hydrogen or helium lines in their spectra) originate from core-collapses of very massive stars (> 30 Msun) or from less massive stars in binary systems. Here we show that Type II (with hydrogen lines) and Ic supernovae are located in environments with similar molecular gas densities, therefore their progenitors have comparable lifetimes and initial masses. This supports a binary interaction for most Type Ic supernova progenitors, which explains the lack of hydrogen and helium lines. This finding can be implemented in sub-grid prescriptions in numerical cosmological simulations to improve the feedback and chemical mixing., Comment: 9 pages, 3 figures, published in Nature Communications

Long GRB hosts at z<1 are usually low-mass, low metallicity star-forming galaxies. Here we present the until now most detailed, spatially resolved study of the host of GRB 171205A, a grand-design barred spiral galaxy at z=0.036. Our analysis includes MUSE integral field spectroscopy, complemented by high spatial resolution UV/VIS HST imaging and CO(1-0) and HI 21cm data. The GRB is located in a small star-forming region in a spiral arm of the galaxy at a deprojected distance of ~ 8 kpc from the center. The galaxy shows a smooth negative metallicity gradient and the metallicity at the GRB site is half solar, slightly below the mean metallicity at the corresponding distance from the center. Star formation in this galaxy is concentrated in a few HII regions between 5-7 kpc from the center and at the end of the bar, inwards of the GRB region, however, the HII region hosting the GRB is in the top 10% of regions with highest specific star-formation rate. The stellar population at the GRB site has a very young component (< 5 Myr) contributing a significant part of the light. Ionized and molecular gas show only minor deviations at the end of the bar. A parallel study found an asymmetric HI distribution and some additional gas near the position of the GRB, which might explain the star-forming region of the GRB site. Our study shows that long GRBs can occur in many types of star-forming galaxies, however, the actual GRB sites consistently have low metallicity, high star formation and a young population. Furthermore, gas inflow or interactions triggering the star formation producing the GRB progenitor might not be evident in ionized or even molecular gas but only in HI., Comment: 22 pages, 18 figures, Appendix with additional figures, A&A under revision

The relatively red wavelength range (4800-9300{\AA}) of the VLT Multi Unit Spectroscopic Explorer (MUSE) limits which metallicity diagnostics can be used; in particular excluding those requiring the [O ii]{\lambda}{\lambda}3726,29 doublet. We assess various strong line diagnostics by comparing to sulphur Te-based metallicity measurements for a sample of 671 HII regions from 36 nearby galaxies from the MUSE Atlas of Disks (MAD) survey. We find that the O3N2 and N2 diagnostics return a narrower range of metallicities which lie up to ~0.3 dex below Te-based measurements, with a clear dependence on both metallicity and ionisation parameter. The N2S2H{\alpha} diagnostic shows a near-linear relation with the Te-based metallicities, although with a systematic downward offset of ~0.2 dex, but no clear dependence on ionisation parameter. These results imply that the N2S2H{\alpha} diagnostic produces the most reliable results when studying the distribution of metals within galaxies with MUSE. On sub-HII region scales, the O3N2 and N2 diagnostics measure metallicity decreasing towards the centres of HII regions, contrary to expectations. The S-calibration and N2S2H{\alpha} diagnostics show no evidence of this, and show a positive relationship between ionisation parameter and metallicity at 12 + log(O/H)> 8.4, implying the relationship between ionisation parameter and metallicity differs on local and global scales. We also present HIIdentify, a python tool developed to identify HII regions within galaxies from H{\alpha} emission maps. All segmentation maps and measured emission line strengths for the 4408 HII regions identified within the MAD sample are available to download., Comment: Accepted for publication in MNRAS. 21 pages, 13 figures

Much of what is known of the chemical composition of the universe is based on emission line spectra from star forming galaxies. Emission-based inferences are, nevertheless, model-dependent and they are dominated by light from luminous star forming regions. An alternative and sensitive probe of the metallicity of galaxies is through absorption lines imprinted on the luminous afterglow spectra of long gamma ray bursts (GRBs) from neutral material within their host galaxy. We present results from a JWST/NIRSpec programme to investigate for the first time the relation between the metallicity of neutral gas probed in absorption by GRB afterglows and the metallicity of the star forming regions for the same host galaxy sample. Using an initial sample of eight GRB host galaxies at z=2.1-4.7, we find a tight relation between absorption and emission line metallicities when using the recently proposed $\hat{R}$ metallicity diagnostic (+/-0.2dex). This agreement implies a relatively chemically-homogeneous multi-phase interstellar medium, and indicates that absorption and emission line probes can be directly compared. However, the relation is less clear when using other diagnostics, such as R23 and R3. We also find possible evidence of an elevated N/O ratio in the host galaxy of GRB090323 at z=3.58, consistent with what has been seen in other $z>4$ galaxies. Ultimate confirmation of an enhanced N/O ratio and of the relation between absorption and emission line metallicities will require a more direct determination of the emission line metallicity via the detection of temperature-sensitive auroral lines in our GRB host galaxy sample., Comment: Accepted for publication in MNRAS; 24 pages, 15 figures

The chemical enrichment of dust and metals in the interstellar medium (ISM) of galaxies throughout cosmic time is one of the key driving processes of galaxy evolution. Here we study the evolution of the gas-phase metallicities, dust-to-gas (DTG), and dust-to-metal (DTM) ratios of 36 star-forming galaxies at $1.7 < z < 6.3$ probed by gamma-ray bursts (GRBs). We compile all GRB-selected galaxies with intermediate (R=7000) to high (R>40,000) resolution spectroscopic data for which at least one refractory (e.g. Fe) and one volatile (e.g. S or Zn) element have been detected at S/N>3. This is to ensure that accurate abundances and dust depletion patterns can be obtained. We first derive the redshift evolution of the dust-corrected, absorption-line based gas-phase metallicity [M/H]$_{\rm tot}$ in these galaxies, for which we determine a linear relation with redshift ${\rm [M/H]_{tot}}(z) = (-0.21\pm 0.04)z -(0.47\pm 0.14)$. We then examine the DTG and DTM ratios as a function of redshift and through three orders of magnitude in metallicity, quantifying the relative dust abundance both through the direct line-of-sight visual extinction $A_V$ and the derived depletion level. We use a novel method to derive the DTG and DTM mass ratios for each GRB sightline, summing up the mass of all the depleted elements in the dust-phase. We find that the DTG and DTM mass ratios are both strongly correlated with the gas-phase metallicity and show a mild evolution with redshift as well. While these results are subject to a variety of caveats related to the physical environments and the narrow pencil-beam sightlines through the ISM probed by the GRBs, they provide strong implications for studies of dust masses to infer the gas and metal content of high-redshift galaxies, and particularly demonstrate the large offset from the average Galactic value in the low-metallicity, high-redshift regime., Comment: Accepted in A&A

GRB 220627A is a rare burst with two distinct gamma-ray emission episodes separated by almost 1000 s that triggered the Fermi Gamma-ray Burst Monitor twice. High-energy GeV emission was detected by the Fermi Large Area Telescope coincident with the first emission episode but not the second. The discovery of the optical afterglow with MeerLICHT led to MUSE observations which secured the burst redshift to z=3.08, making this the most distant ultra-long gamma-ray burst (GRB) detected to date. The progenitors of some ultra-long GRBs have been suggested in the literature to be different to those of normal long GRBs. Our aim is to determine whether the afterglow and host properties of GRB 220627A agree with this interpretation. We performed empirical and theoretical modelling of the afterglow data within the external forward shock framework, and determined the metallicity of the GRB environment through modelling the absorption lines in the MUSE spectrum. Our optical data show evidence for a jet break in the light curve at ~1.2 days, while our theoretical modelling shows a preference for a homogeneous circumburst medium. Our forward shock parameters are typical for the wider GRB population, and we find that the environment of the burst is characterised by a sub-solar metallicity. Our observations and modelling of GRB 220627A do not suggest that a different progenitor compared to the progenitor of normal long GRBs is required. We find that more observations of ultra-long GRBs are needed to determine if they form a separate population with distinct prompt and afterglow features, and possibly distinct progenitors., Comment: 17 pages, accepted for publication in A&A

We analyse the joint distribution of dust attenuation and projected axis ratios, together with galaxy size and surface brightness profile information, to infer lessons on the dust content and star/dust geometry within star-forming galaxies at 0 < z <2.5. To do so, we make use of large observational datasets from KiDS+VIKING+HSC-SSP and extend the analysis out to redshift z = 2.5 using the HST surveys CANDELS and 3D-DASH. We construct suites of SKIRT radiative transfer models for idealized galaxies observed under random viewing angles with the aim of reproducing the aforementioned distributions, including the level and inclination dependence of dust attenuation. We find that attenuation-based dust mass estimates are at odds with constraints from far-infrared observations, especially at higher redshifts, when assuming smooth star and dust geometries of equal extent. We demonstrate that UV-to-near-IR and far-infrared constraints can be reconciled by invoking clumpier dust geometries for galaxies at higher redshifts and/or very compact dust cores. We discuss implications for the significant wavelength- and redshift-dependent differences between half-light and half-mass radii that result from spatially varying dust columns within -- especially massive -- star-forming galaxies., Comment: Accepted for publication in MNRAS

We report the discovery of Swift J221951-484240 (hereafter: J221951), a luminous slow-evolving blue transient that was detected by the Neil Gehrels Swift Observatory Ultra-violet/Optical Telescope (Swift/UVOT) during the follow-up of Gravitational Wave alert S190930t, to which it is unrelated. Swift/UVOT photometry shows the UV spectral energy distribution of the transient to be well modelled by a slowly shrinking black body with an approximately constant temperature of T~2.5x10^4 K. At a redshift z=0.5205, J221951 had a peak absolute magnitude of M_u,AB = -23 mag, peak bolometric luminosity L_max=1.1x10^45 erg s^-1 and a total radiated energy of E>2.6x10^52 erg. The archival WISE IR photometry shows a slow rise prior to a peak near the discovery date. Spectroscopic UV observations display broad absorption lines in N V and O VI, pointing toward an outflow at coronal temperatures. The lack of emission in the higher H~Lyman lines, N I and other neutral lines is consistent with a viewing angle close to the plane of the accretion or debris disc. The origin of J221951 can not be determined with certainty but has properties consistent with a tidal disruption event and the turn-on of an active galactic nucleus., Comment: 37 pages (25 main + 12 supplementary), submitted to MNRAS

Context: The extreme luminosity of gamma-ray bursts (GRBs) makes them powerful beacons for studies of the distant Universe. The most luminous bursts are typically detected at moderate/high redshift, where the volume for seeing such rare events is maximized and the star-formation activity is greater than at z = 0. For distant events, not all observations are feasible, such as at TeV energies. Aims: Here we present a spectroscopic redshift measurement for the exceptional GRB 221009A, the brightest GRB observed to date with emission extending well into the TeV regime. Methods: We used the X-shooter spectrograph at the ESO Very Large Telescope (VLT) to obtain simultaneous optical to near-IR spectroscopy of the burst afterglow 0.5 days after the explosion. Results: The spectra exhibit both absorption and emission lines from material in a host galaxy at z = 0.151. Thus GRB 221009A was a relatively nearby burst with a luminosity distance of 745 Mpc. Its host galaxy properties (star-formation rate and metallicity) are consistent with those of LGRB hosts at low redshift. This redshift measurement yields information on the energy of the burst. The inferred isotropic energy release, $E_{\rm iso} > 5 \times 10^{54}$ erg, lies at the high end of the distribution, making GRB 221009A one of the nearest and also most energetic GRBs observed to date. We estimate that such a combination (nearby as well as intrinsically bright) occurs between once every few decades to once per millennium., Comment: 9 pages, 4 figures, submitted to Astronomy & Astrophysics

The study of the properties of galaxies in the first billion years after the Big Bang is one of the major topic of current astrophysics. Optical/near-infrared spectroscopy of the afterglows of long Gamma-ray bursts (GRBs) provide a powerful diagnostic tool to probe the interstellar medium (ISM) of their host galaxies and foreground absorbers, even up to the highest redshifts. We analyze the VLT/X-shooter afterglow spectrum of GRB 210905A, triggered by the Swift Neil Gehrels Observatory, and detect neutral-hydrogen, low-ionization, high-ionization, and fine-structure absorption lines from a complex system at z=6.3118, that we associate with the GRB host galaxy. We study the ISM properties of the host system, revealing the metallicity, kinematics and chemical abundance pattern. The total metallicity of the z~6.3 system is [M/H]=-1.72+/-0.13, after correcting for dust-depletion and taking into account alpha-element enhancement. In addition, we determine the overall amount of dust and dust-to-metal mass ratio (DTM) ([Zn/Fe]_fit=0.33+/-0.09, DTM=0.18+/-0.03). We find indications of nucleosynthesis due to massive stars and evidence of peculiar over-abundance of aluminium. From the analysis of fine-structure lines, we determine distances of several kpc for the low-ionization gas clouds closest to the GRB. Those farther distances are possibly due to the high number of ionizing photons. Using the HST/F140W image of the GRB field, we show the GRB host galaxy as well as multiple objects within 2" from the GRB. We discuss the galaxy structure and kinematics that could explain our observations, also taking into account a tentative detection of Lyman-alpha emission. Deep spectroscopic observations with VLT/MUSE and JWST will offer the unique possibility of combining our results with the ionized-gas properties, with the goal of better understanding how galaxies in the reionization era form and evolve., Comment: Accepted Publication (In Press on A&A) - 22 pages, 10 figures, 6 tables - Appendix: 6 figures, 3 tables

Long ($>2$ s) gamma ray bursts (GRBs) are associated with explosions of massive stars, although in three instances, supernovae (SNe) have not been detected, despite deep observations. With new HI line and archival optical integral field spectroscopy data, we characterize the interstellar medium (ISM) of the host galaxy of one of these events, GRB 111005A, in order to shed light on the unclear nature of these peculiar objects. We found that the atomic gas, radio continuum, and rotational patterns are in general very smooth throughout the galaxy, which does not indicate a recent gas inflow or outflow. There is also no gas concentration around the GRB position. The ISM in this galaxy differs from that in hosts of other GRBs and SNe, which may suggest that the progenitor of GRB 111005A was not an explosion of a very massive star (e.g. a compact object merger). However, there are subtle irregularities of the GRB 111005A host (most at a $2\sigma$ level), which may point to a weak gas inflow or interaction. Since in the SE part of the host there is 15% more atomic gas and twice less molecular gas than in NW part, the molecular gas fraction is low. In the SE part there is also a region with very high H$\alpha$ equivalent width. There is more continuum 1.4 GHz emission to the SE and an S-shaped warp in the UV. Finally, there is also a low-metallicity region 3.5" (1 kpc) from the GRB position. Two galaxies within 300 kpc or a past merger can be responsible for these irregularities., Comment: Accepted for publication in ApJS, 18 pages, 10 figures

Radial metallicity trends provide a key indicator of physical processes such as star formation and radial gas migration within a galaxy. Large IFU surveys allow for detailed studies of these radial variations, with recent observations detecting central dips in the metallicity, which may trace the impact of various evolutionary processes. However, the origin of these dips has not been conclusively determined, with suggestions that they may be diagnostic dependent. In this paper, we use the SDSS-IV MaNGA survey to investigate whether the observed dips represent genuine decreases in the central metallicity, or if they could be an artefact of the diagnostic used. Using a sub-sample of 758 local star-forming galaxies at low inclinations, we investigate in detail the impact of using different strong line diagnostics on the shapes of the returned profiles, and the prevalence of dips. We find no clear evidence of the dips being caused by changing values of the ionisation parameter within galaxies. To investigate physical causes, we explore both global and spatially-resolved parameters, finding that galaxies exhibiting central dips in the O3N2 metallicity profile have on average lower H$\alpha$EW values out to $R/R_\rm{e} \sim 1.5$, and higher values of D$_N$(4000) in the central regions. We additionally find a higher prevalence of dips in galaxies with high stellar mass, and lower values of global specific star formation rate, suggesting a possible link to central quenching. Nevertheless, these results are dependent on the diagnostic used, suggesting caution should be taken when interpreting observed features in galaxy metallicity gradients., Comment: Accepted for publication in MNRAS. 23 pages; 21 figures

Short-GRB progenitors could come in various flavors, depending on the nature of the merging compact stellar objects (including a stellar-mass black hole or not) or depending on their ages (millions or billions of years). At a redshift of z=0.122, the nearly face-on spiral host of the short GRB 080905A is one of the closest short-GRB host galaxies identified so far. This made it a preferred target to explore spatially resolved star-formation and to investigate the afterglow position in the context of its star formation structures. We used VLT/MUSE integral-field unit observations, supplemented by ATCA 5.5/9.0 GHz radio-continuum measurements and publicly available HST data, to study the star-formation activity in the GRB 080905A host galaxy. The MUSE observations reveal that the entire host is characterized by strong line emission. Using the Halpha line flux, we measure for the entire galaxy an SFR of about 1.6 Msun/yr, consistent with its non-detection by ATCA. Several individual star-forming regions are scattered across the host. The most luminous region has a Halpha luminosity that is nearly four times as high as the luminosity of the Tarantula nebula in the Large Magellanic Cloud. Even though star-forming activity can be traced as close to about 3 kpc (in projection) distance to the GRB explosion site, stellar population synthesis calculations show that none of the Halpha-bright star-forming regions is a likely birthplace of the short-GRB progenitor., Comment: The Astrophysical Journal Supplement Series, in press

Short Gamma-Ray Bursts (SGRBs) are produced by the coalescence of compact binary systems which are remnants of massive stars. GRB 160410A is classified as a short-duration GRB with extended emission and is currently the farthest SGRB with a redshift determined from an afterglow spectrum and also one of the brightest SGRBs to date. The fast reaction to the Neil Gehrels Swift Observatory alert allowed us to obtain a spectrum of the afterglow using the X-shooter spectrograph at the Very Large Telescope (VLT). The spectrum shows several absorption features at a redshift of z=1.7177, in addition, we detect two intervening systems at z=1.581 and z=1.444. The spectrum shows ly-alpha in absorption with a column density of log N(HI)=21.2+/-0.2 cm$^{-2}$ which, together with FeII, CII, SiII, AlII and OI, allow us to perform the first study of chemical abundances in a SGRB host galaxy. We determine a metallicity of [X/H]=-2.3+/-0.2 for FeII and -2.5+/-0.2 for SiII and no dust depletion. We also find no evidence for extinction in the afterglow spectral energy distribution (SED) modeling. The environment has a low degree of ionisation and the CIV and SiIV lines are completely absent. We do not detect an underlying host galaxy down to deep limits. Additionally, we compare GRB 160410A to GRB 201221D, another high-z short GRB that shows absorption lines at z=1.045 and an underlying massive host galaxy., Comment: 25 pages, 22 figures. Accepted for publication in Monthly Notices of the Royal Astronomical Society (MNRAS)

We investigate the thermal emission and extinction from dust associated with the nearby superluminous supernova (SLSN) 2018bsz. Our dataset has daily cadence and simultaneous optical and near-infrared coverage up to ~ 100 days, together with late time (+ 1.7 yr) MIR observations. At 230 days after light curve peak the SN is not detected in the optical, but shows a surprisingly strong near-infrared excess, with r - J > 3 mag and r - Ks > 5 mag. The time evolution of the infrared light curve enables us to investigate if the mid-infrared emission is from newly formed dust inside the SN ejecta, from a pre-existing circumstellar envelope, or interstellar material heated by the radiation from the SN. We find the latter two scenarios can be ruled out, and a scenario where new dust is forming in the SN ejecta at epochs > 200 days can self-consistently reproduce the evolution of the SN flux. We can fit the spectral energy distribution well at +230 d with 5 x 10^-4 solar mass of carbon dust, increasing over the following several hundred days to 10^-2 solar mass by +535 d. SN 2018bsz is the first SLSN showing evidence for dust formation within the SN ejecta, and appears to form ten times more dust than normal core-collapse SNe at similar epochs. Together with their preference for low mass, low metallicity host galaxies, we suggest that SLSNe may be a significant contributor to dust formation in the early Universe., Comment: 14 pages and 7 figures in main text, 12 pages and 6 figures in appendix. The observational data will be updated once the paper is accepted

We present a pilot search of CO emission in three H$_2$-absorbing, long-duration gamma-ray burst (GRB) host galaxies at z~2-3. We used the Atacama Large Millimeter/sub-millimeter Array (ALMA) to target the CO(3-2) emission line and report non-detections for all three hosts. These are used to place limits on the host molecular gas masses, assuming a metallicity-dependent CO-to-H$_2$ conversion factor ($\alpha_{\rm CO}$). We find, $M_{\rm mol} < 3.5\times 10^{10}\,M_{\odot}$ (GRB\,080607), $M_{\rm mol} < 4.7\times 10^{11}\,M_{\odot}$ (GRB\,120815A), and $M_{\rm mol} < 8.9\times 10^{11}\,M_{\odot}$ (GRB\,181020A). The high limits on the molecular gas mass for the latter two cases are a consequence of their low stellar masses $M_\star$ ($M_\star \lesssim 10^{8}\,M_{\odot}$) and low gas-phase metallicities ($Z\sim 0.03\,Z_{\odot}$). The limit on the $M_{\rm mol}/M_\star$ ratio derived for GRB\,080607, however, is consistent with the average population of star-forming galaxies at similar redshifts and stellar masses. We discuss the broader implications for a metallicity-dependent CO-to-H$_2$ conversion factor, and demonstrate that the canonical Galactic $\alpha_{\rm CO}$, will severely underestimate the actual molecular gas mass for all galaxies at $z>1$ with $M_\star < 10^{10}\,M_\odot$. To better quantify this we develop a simple approach to estimate the relevant $\alpha_{\rm CO}$ factor based only on the redshift and stellar mass of individual galaxies. The elevated conversion factors will make these galaxies appear CO-"dark" and difficult to detect in emission, as is the case for the majority of GRB hosts. GRB spectroscopy thus offers a complementary approach to identify low-metallicity, star-forming galaxies with abundant molecular gas reservoirs at high redshifts that are otherwise missed by current ALMA surveys., Comment: 7 pages, 5 Figs, accepted for publication in MNRAS

Large-scale surveys over the last years have revealed about 300 QSOs at redshift above 6. Follow-up observations identified surprising properties, such as the very high black hole (BH) masses, spatial correlations with surrounding cold gas of the host galaxy, or high CIV-MgII velocity shifts. In particular, the discovery of luminous high-redshift quasars suggests that at least some black holes likely have large masses at birth and grow efficiently. We aim at quantifying quasar pairs at high redshift for a large sample of objects. This provides a new key constraint on a combination of parameters related to the origin and assembly for the most massive black holes: BH formation efficiency and clustering, growth efficiency and relative contribution of BH mergers. We observed 116 spectroscopically confirmed QSOs around redshift 6 with the simultaneous 7-channel imager GROND in order to search for companions. Applying identical colour-colour cuts as for those which led to the spectroscopically confirmed QSO, we perform LePHARE fits to the 26 best QSO pair candidates, and obtained spectroscopic observations for 11 of those. e do not find any QSO pair with a companion brighter than M1450(AB) < -26 mag within our 0.1-3.3 h^-1 cMpc search radius, in contrast to the serendipitous findings in the redshift range 4--5. However, a low fraction of such pairs at this luminosity and redshift is consistent with indications from present-day cosmological-scale galaxy evolution models. In turn, the incidence of L- and T-type brown dwarfs which occupy a similar colour space as z ~ 6 QSOs, is higher than expected, by a factor of 5 and 20, respectively., Comment: 20 pages, 18 figures, A&A (accepted)

We present a modified version of the L-GALAXIES 2020 semi-analytic model of galaxy evolution, which includes significantly increased direct metal enrichment of the circumgalactic medium (CGM) by supernovae (SNe). These more metal-rich outflows do not require increased mass-loading factors, in contrast to some other galaxy evolution models. This modified L-GALAXIES 2020 model is able to simultaneously reproduce the gas-phase metallicity $(Z_{\rm g})$ and stellar metallicity $(Z_{*})$ radial profiles observed in nearby disc galaxies by MaNGA and MUSE, as well as the observed mass - metallicity relations for gas and stars at $z=0$ and their evolution back to $z\sim{}2-3$. A direct CGM enrichment fraction of $\sim{}90\%$ for SNe-II is preferred. We find that massive disc galaxies have slightly flatter $Z_{\rm g}$ profiles than their lower-mass counterparts in L-GALAXIES 2020, due to more efficient enrichment of their outskirts via inside-out growth and metal-rich accretion. Such a weak, positive correlation between stellar mass and $Z_{\rm g}$ profile slope is also seen in our MaNGA-DR15 sample of 571 star-forming disc galaxies. Although, below ${\rm log}(M_{*}/{\rm M}_{\odot})\sim{}10.0$ this observational result is strongly dependent on the metallicity diagnostic and morphological selection chosen. In addition, a lowered maximum SN-II progenitor mass of $25{\rm M}_{\odot}$, reflecting recent theoretical and observational estimates, can also provide a good match to observed metallicity profiles at $z=0$ in L-GALAXIES 2020. However, this model version fails to reproduce an evolution in $Z_{\rm g}$ at fixed mass over cosmic time, or the magnesium abundances observed in the intracluster medium (ICM)., Comment: 22 pages, 19 figures. Accepted by MNRAS. (Missing citation hyperlinks corrected)

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

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

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

We present comprehensive multiwavelength radio to X-ray observations of GRB 181201A spanning from $\approx150$ s to $\approx163$ days after the burst, comprising the first joint ALMA-VLA-GMRT observations of a gamma-ray burst (GRB) afterglow. The radio and mm-band data reveal a distinct signature at $\approx3.9$ days, which we interpret as reverse shock (RS) emission. Our observations present the first time that a single radio-frequency spectral energy distribution can be decomposed directly into RS and forward shock (FS) components. We perform detailed modeling of the full multiwavelength data set, using Markov Chain Monte Carlo sampling to construct the joint posterior density function of the underlying physical parameters describing the RS and FS synchrotron emission. We uncover and account for all degeneracies in the model parameters. The joint RS-FS modeling reveals a weakly magnetized ($\sigma\approx3\times10^{-3}$), mildly relativistic RS, from which we derive an initial bulk Lorentz factor of $\Gamma_0\approx103$ for the GRB jet. Our results support the hypothesis that low-density environments are conducive to the observability of RS emission. We compare our observations to other events with strong RS detections, and find a likely observational bias selecting for longer lasting, non-relativistic reverse shocks. We present and begin to address new challenges in modeling posed by the present generation of comprehensive, multi-frequency data sets., Comment: Accepted version

There is observational evidence of a dearth in core-collapse supernova (ccSN) explosions from stars with zero age main sequence (ZAMS) mass M_0~17-30 Msol, referred to as the 'red supergiant problem'. However, simulations now predict that above 20Msol we should indeed only expect stars within certain pockets of M_0 to produce a visible SN explosion. Validating these predictions requires large numbers of ccSNe of different types with measured M_0, which is challenging. In this paper we explore the reliability of using host galaxy emission lines and the Halpha equivalent width to constrain the age, and thus the M_0 of ccSNe progenitors. We use Binary Population and Spectral Synthesis models to infer a stellar population age from MUSE observations of the ionised gas properties and Halpha EW at the location of eleven ccSNe with reliable M_0 measurements. Comparing our results to published M_0 values, we find that models that do not consider binary systems yield stellar ages that are systematically too young (thus M_0 too large), whereas accounting for binary system interactions typically overpredict the stellar age (thus underpredict M_0). Taking into account the effects of photon leakage bring our M_0 estimates in much closer agreement with expectations. These results highlight the need for careful modelling of diffuse environments, such as are present in the vicinity of type II SNe, before ionised emission line spectra can be used as reliable tracers of progenitor stellar age., Comment: 17 pages and 5 figures (excluding appendix). Replaced to match published version in MNRAS

We present ALMA 97.5 GHz total intensity and linear polarization observations of the mm-band afterglow of GRB 190114C spanning 2.2 to 5.2 hours after the burst. We detect linear polarization at the $\approx 5\,\sigma$ level, decreasing from $\Pi=(0.87\pm0.13)\%$ to $(0.60\pm0.19)\%$, and evolving in polarization position angle from $(10\pm5)^\circ$ to $(-44\pm12)^\circ$ during the course of the observations. This represents the first detection of polarized millimeter emission in a $\gamma$-ray burst. We show that the optical and X-ray observations between $0.03$ days and $\sim0.3$ days are consistent with a fast cooling forward shock expanding into a wind environment. However, the optical observations at $\lesssim0.03$ days, as well as the radio and millimeter observations arise from a separate component, which we interpret as emission from the reverse-shocked ejecta. Using the measured linear polarization, we constrain the coherence scale of tangled magnetic fields in the ejecta to an angular size of $\theta_{\rm B} \approx10^{-3}$ radian, while the rotation of the polarization angle rules out the presence of large scale, ordered axisymmetric magnetic fields, and in particular a large scale toroidal field, in the jet., Comment: Published in ApJ Letters

The number of long-duration Gamma Ray Burst (LGRB) host galaxies with measured metallicities and host masses has now grown to over one hundred, allowing us to investigate how the distributions of both these properties change with redshift. Using the combined host galaxy metallicity sample from Graham & Fruchter (2013) and Kr\"uhler et al. (2015), we find a surprising lack of evolution in the LGRB metallicity distribution across different redshifts. In particular, the fraction of LGRB hosts with relatively high metallicity (12+log(O/H) $\geq$ 8.4) remains essentially constant out to z = 2.5. This result is at odds with the evolution in the mass-metallicity relation of typical galaxies, which become progressively more metal poor with increasing redshift. A similar result is found when converting the LGRB host galaxy mass distribution taken from the SHOALS (Swift GRB Host Galaxy Legacy Survey) sample to a corresponding metallicity distribution by applying a redshift-dependent mass-metallicity relation. The SHOALS sample is compiled using an unbiased selection function implying that the observed lack of evolution in the host galaxy high metallicity distribution is not caused by selection effects. However, the LGRB host galaxy metallicities estimated from the stellar mass are typically a quarter dex higher at all redshifts than the metallicity we measure spectroscopically. This implies that using mass-metallicity relationships to estimate host metallicities will thus produce a substantial systematic bias., Comment: 23 pages, 13 figures, 12 tables, 3 appendices, ApJ in press

Unusual stellar explosions represent an opportunity to learn about both stellar and galaxy evolution. Mapping the atomic gas in host galaxies of such transients can lead to an understanding of the conditions triggering them. We provide resolved atomic gas observations of the host galaxy, CGCG137-068, of the unusual, poorly-understood transient AT2018cow searching for clues to understand its nature. We test whether it is consistent with a recent inflow of atomic gas from the intergalactic medium, as suggested for host galaxies of gamma-ray bursts (GRBs) and some supernovae (SNe). We observed the HI hyperfine structure line of the AT2018cow host with the Giant Metrewave Radio Telescope. There is no unusual atomic gas concentration near the position of AT2018cow. The gas distribution is much more regular than those of GRB/SN hosts. The AT2018cow host has an atomic gas mass lower by 0.24 dex than predicted from its star formation rate (SFR) and is at the lower edge of the galaxy main sequence. In the continuum we detected the emission of AT2018cow and of a star-forming region in the north-eastern part of the bar (away from AT2018cow). This region hosts a third of the galaxy's SFR. The absence of atomic gas concentration close to AT2018cow, along with a normal SFR and regular HI velocity field, sets CGCG137-068 apart from GRB/SN hosts studied in HI. The environment of AT2018cow therefore suggests that its progenitor may not have been a massive star. Our findings are consistent with an origin of the transient that does not require a connection between its progenitor and gas concentration or inflow: an exploding low-mass star, a tidal disruption event, a merger of white dwarfs, or a merger between a neutron star and a giant star. We interpret the recently reported atomic gas ring in CGCG137-068 as a result of internal processes connected with gravitational resonances caused by the bar., Comment: Accepted to A&A. 9 pages, 4 figures, 2 tables. Main changes: 1) inclusion of the HI data from August 2018; 2) the atomic gas ring discovered by Roychowdhury et al. (2019) is visible, but we present a different interpretation of its origin (internal mechanisms); 3) the faint low-S/N companions were spurious

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

We investigate electron temperature (Te) and gas-phase oxygen abundance (Z_Te) measurements for galaxies in the local Universe (z < 0.25). Our sample comprises spectra from a total of 264 emission-line systems, ranging from individual HII regions to whole galaxies, including 23 composite HII regions from "star-forming main sequence" galaxies in the MaNGA survey. We utilise 130 of these systems with directly measurable Te(OII) to calibrate a new metallicity-dependent Te(OIII) - Te(OII) relation that provides a better representation of our varied dataset than existing relations from the literature. We also provide an alternative Te(OIII) - Te(NII) calibration. This new Te method is then used to obtain accurate Z_Te estimates and form the mass - metallicity relation (MZR) for a sample of 118 local galaxies. We find that all the Te(OIII) - Te(OII) relations considered here systematically under-estimate Z_Te for low-ionisation systems by up to 0.6 dex. We determine that this is due to such systems having an intrinsically higher O+ abundance than O++ abundance, rendering Z_Te estimates based only on [OIII] lines inaccurate. We therefore provide an empirical correction based on strong emission lines to account for this bias when using our new Te(OIII) - Te(OII) and Te(OIII) - Te(NII) relations. This allows for accurate metallicities (1sigma = 0.08 dex) to be derived for any low-redshift system with an [OIII]4363 detection, regardless of its physical size or ionisation state. The MZR formed from our dataset is in very good agreement with those formed from direct measurements of metal recombination lines and blue supergiant absorption lines, in contrast to most other Te-based and strong-line-based MZRs. Our new Te method therefore provides an accurate and precise way of obtaining Z_Te for a large and diverse range of star-forming systems in the local Universe., Comment: 15 pages (excluding appendices), 18 figures, 4 tables. Accepted for publication in A&A

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

We present observations of supernova (SN) 2017ens, discovered by the ATLAS survey and identified as a hot blue object through the GREAT program. The redshift z=0.1086 implies a peak brightness of M_g=-21.1 mag, placing the object within the regime of superluminous supernovae. We observe a dramatic spectral evolution, from initially being blue and featureless, to later developing features similar to those of the broadlined Type Ic SN 1998bw, and finally showing ~2000 km s^-1 wide H-alpha and H-beta emission. Relatively narrow Balmer emission (reminiscent of a SN IIn) is present at all times. We also detect coronal lines, indicative of a dense circumstellar medium. We constrain the progenitor wind velocity to ~50-60 km s^-1 based on P-Cygni profiles, which is far slower than those present in Wolf-Rayet stars. This may suggest that the progenitor passed through a luminous blue variable phase, or that the wind is instead from a binary companion red supergiant star. At late times we see the ~2000 km s^-1 wide H-alpha emission persisting at high luminosity (~3x10^40 erg s^-1) for at least 100 day, perhaps indicative of additional mass loss at high velocities that could have been ejected by a pulsational pair instability., Comment: Accepted version in ApJL

Twenty years ago, GRB 980425/SN 1998bw revealed that long Gamma-Ray Bursts (GRBs) are physically associated with broad-lined type Ic supernovae. Since then more than 1000 long GRBs have been localized to high angular precision, but only in about 50 cases the underlying supernova (SN) component was identified. Using the multi-channel imager GROND (Gamma-Ray Burst Optical Near-Infrared Detector) at ESO/La Silla, during the last ten years we have devoted a substantial amount of observing time to reveal and to study SN components in long-GRB afterglows. Here we report on four more GRB-SNe (associated with GRBs 071112C, 111228A, 120714B, and 130831A) which were discovered and/or followed-up with GROND and whose redshifts lie between z=0.4 and 0.8. We study their afterglow light curves, follow the associated SN bumps over several weeks, and characterize their host galaxies. Using SN 1998bw as a template, the derived SN explosion parameters are fully consistent with the corresponding properties of the so-far known GRB-SN ensemble, with no evidence for an evolution of their properties as a function of redshift. In two cases (GRB 120714B/SN 2012eb at z=0.398 and GRB 130831A/SN 2013fu at z=0.479) additional Very Large Telescope (VLT) spectroscopy of the associated SNe revealed a photospheric expansion velocity at maximum light of about 40 000 and 20 000 km/s, respectively. For GRB 120714B, which was an intermediate-luminosity burst, we find additional evidence for a blackbody component in the light of the optical transient at early times, similar to what has been detected in some GRB-SNe at lower redshifts., Comment: submitted to Astron. Astroph., revised version

We present the ATLAS discovery and initial analysis of the first 18 days of the unusual transient event, ATLAS18qqn/AT2018cow. It is characterized by a high peak luminosity ($\sim$1.7 $\times$ 10$^{44}$ erg s$^{-1}$), rapidly evolving light curves ($>$5 mag rise to peak in $\sim$3.5 days), and hot blackbody spectra, peaking at $\sim$27000 K that are relatively featureless and unchanging over the first two weeks. The bolometric light curve cannot be powered by radioactive decay under realistic assumptions. The detection of high-energy emission may suggest a central engine as the powering source. Using a magnetar model, we estimated an ejected mass of $0.1-0.4$ \msol, which lies between that of low-energy core-collapse events and the kilonova, AT2017gfo. The spectra cooled rapidly from 27000 to 15000 K in just over 2 weeks but remained smooth and featureless. Broad and shallow emission lines appear after about 20 days, and we tentatively identify them as He I although they would be redshifted from their rest wavelengths. We rule out that there are any features in the spectra due to intermediate mass elements up to and including the Fe-group. The presence of r-process elements cannot be ruled out. If these lines are due to He, then we suggest a low-mass star with residual He as a potential progenitor. Alternatively, models of magnetars formed in neutron-star mergers give plausible matches to the data., Comment: Accepted for publication in ApJL, updated to reflect the published version

Long gamma-ray bursts give us the chance to study both their extreme physics and the star-forming galaxies in which they form. GRB 100418A, at a z = 0.6239, had a bright optical and radio afterglow, and a luminous star-forming host galaxy. This allowed us to study the radiation of the explosion as well as the interstellar medium of the host both in absorption and emission. We collected photometric data from radio to X-ray wavelengths to study the evolution of the afterglow and the contribution of a possible supernova and three X-shooter spectra obtained during the first 60 hr. The light curve shows a very fast optical rebrightening, with an amplitude of 3 magnitudes, starting 2.4 hr after the GRB onset. This cannot be explained by a standard external shock model and requires other contributions, such as late central-engine activity. Two weeks after the burst we detect an excess in the light curve consistent with a SN with peak absolute magnitude M_V = -18.5 mag, among the faintest GRB-SNe detected to date. The host galaxy shows two components in emission, with velocities differing by 130 km s^-1, but otherwise having similar properties. While some absorption and emission components coincide, the absorbing gas spans much higher velocities, indicating the presence of gas beyond the star-forming regions. The host has a star-formation rate of 12.2 M_sol yr^-1, a metallicity of 12 + log(O/H) = 8.55 and a mass of 1.6x10^9 M_sol. GRB 100418A is a member of a class of afterglow light curves which show a steep rebrightening in the optical during the first day, which cannot be explained by traditional models. Its very faint associated SN shows that GRB-SNe can have a larger dispersion in luminosities than previously seen. Furthermore, we have obtained a complete view of the host of GRB 100418A owing to its spectrum, which contains a remarkable number of both emission and absorption lines., Comment: 23 pages, 14 figures, accepted for publication in A&A

The extinction profiles in Gamma-Ray Burst (GRB) afterglow spectral energy distributions (SEDs) are usually described by the Small Magellanic Cloud (SMC)-type extinction curve. In different empirical extinction laws, the total-to-selective extinction, RV, is an important quantity because of its relation to dust grain sizes and compositions. We here analyse a sample of 17 GRBs (0.34

We obtained CO(2-1) observations of seven GRB hosts with the APEX and IRAM 30m telescopes. We analysed these data together with all other hosts with previous CO observations. We obtained detections for 3 GRB hosts (980425, 080207, and 111005A) and upper limits for the remaining 4 (031203, 060505, 060814, and 100316D). In our entire sample of 12 CO-observed GRB hosts, 3 are clearly deficient in molecular gas, even taking into account their metallicity (980425, 060814, and 080517). Four others are close to the best-fit line for other star-forming galaxies on the SFR-MH2 plot (051022, 060505, 080207, and 100316D). One host is clearly molecule rich (111005A). Finally, the data for 4 GRB hosts are not deep enough to judge whether they are molecule deficient (000418, 030329, 031203, and 090423). The median value of the molecular gas depletion time, MH2/SFR, of GRB hosts is ~0.3 dex below that of other star-forming galaxies, but this result has low statistical significance. A Kolmogorov-Smirnov test performed on MH2/SFR shows an only ~2sigma difference between GRB hosts and other galaxies. This difference can partly be explained by metallicity effects, since the significance decreases to ~1sigma for MH2/SFR versus~metallicity. We found that any molecular gas deficiency of GRB hosts has low statistical significance and that it can be attributed to their lower metallicities; and thus the sample of GRB hosts has molecular properties that are consistent with those of other galaxies, and they can be treated as representative star-forming galaxies. Given the concentration of atomic gas recently found close to GRB and supernova sites, indicating recent gas inflow, our results about the weak molecular deficiency imply that such an inflow does not enhance the SFRs significantly, or that atomic gas converts efficiently into the molecular phase, which fuels star formation., Comment: Published in A&A. 12 pages, 5 figures, 4 tables

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

Deriving physical parameters from gamma-ray burst afterglow observations remains a challenge, even now, 20 years after the discovery of afterglows. The main reason for the lack of progress is that the peak of the synchrotron emission is in the sub-mm range, thus requiring radio observations in conjunction with X-ray/optical/near-infrared data in order to measure the corresponding spectral slopes and consequently remove the ambiguity wrt. slow vs. fast cooling and the ordering of the characteristic frequencies. We observed GRB 151027B, the 1000th Swift-detected GRB, with GROND in the optical-NIR, ALMA in the sub-millimeter, ATCA in the radio band, and combine this with public Swift-XRT X-ray data. While some observations at crucial times only return upper limits or surprising features, the fireball model is narrowly constrained by our data set, and allows us to draw a consistent picture with a fully-determined parameter set. Surprisingly, we find rapid, large-amplitude flux density variations in the radio band which are extreme not only for GRBs, but generally for any radio source. We interpret these as scintillation effects, though the extreme nature requires either the scattering screen to be at much smaller distance than usually assumed, multiple screens, or a combination of the two., Comment: 10 pages, 8 figures, accepted for publ. in A&A

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

THESEUS is a space mission concept aimed at exploiting Gamma-Ray Bursts for investigating the early Universe and at providing a substantial advancement of multi-messenger and time-domain astrophysics. These goals will be achieved through a unique combination of instruments allowing GRB and X-ray transient detection over a broad field of view (more than 1sr) with 0.5-1 arcmin localization, an energy band extending from several MeV down to 0.3 keV and high sensitivity to transient sources in the soft X-ray domain, as well as on-board prompt (few minutes) follow-up with a 0.7 m class IR telescope with both imaging and spectroscopic capabilities. THESEUS will be perfectly suited for addressing the main open issues in cosmology such as, e.g., star formation rate and metallicity evolution of the inter-stellar and intra-galactic medium up to redshift $\sim$10, signatures of Pop III stars, sources and physics of re-ionization, and the faint end of the galaxy luminosity function. In addition, it will provide unprecedented capability to monitor the X-ray variable sky, thus detecting, localizing, and identifying the electromagnetic counterparts to sources of gravitational radiation, which may be routinely detected in the late '20s / early '30s by next generation facilities like aLIGO/ aVirgo, eLISA, KAGRA, and Einstein Telescope. THESEUS will also provide powerful synergies with the next generation of multi-wavelength observatories (e.g., LSST, ELT, SKA, CTA, ATHENA)., Comment: Accepted for publication in Advances in Space Research. Partly based on the proposal submitted on October 2016 in response to the ESA Call for next M5 mission, with expanded and updated science sections

Dust is known to be produced in the envelopes of AGB stars, the expanded shells of supernova (SN) remnants, and in situ grain growth in the ISM, although the corresponding efficiency of each of these dust formation mechanisms at different redshifts remains a topic of debate. During the first Gyr after the Big Bang, it is widely believed that there was not enough time to form AGB stars in high numbers, so that the dust at this epoch is expected to be purely from SNe, or subsequent grain growth in the ISM. The time period corresponding to z ~5-6 is thus expected to display the transition from SN-only dust to a mixture of both formation channels as we know it today. Here we aim to use afterglow observations of GRBs at redshifts larger than $z > 4$ in order to derive host galaxy dust column densities along their line-of-sight and to test if a SN-type dust extinction curve is required for some of the bursts. GRB afterglow observations were performed with the 7-channel GROND Detector at the 2.2m MPI telescope in La Silla, Chile and combined with data gathered with XRT. We increase the number of measured $A_V$ values for GRBs at z > 4 by a factor of ~2-3 and find that, in contrast to samples at mostly lower redshift, all of the GRB afterglows have a visual extinction of $A_V$ < 0.5 mag. Analysis of the GROND detection thresholds and results from a Monte-Carlo simulation show that, although we partly suffer from an observational bias against highly extinguished sight-lines, GRB host galaxies at 4 < z < 6 seem to contain on average less dust than at z ~ 2. Additionally, we find that all of the GRBs can be modeled with locally measured extinction curves and that the SN-like dust extinction curve provides a better fit for only two of the afterglow SEDs. For the first time we also report a photometric redshift of $z = 7.88$ for GRB 100905A, making it one of the most distant GRBs known to date., Comment: 26 pages, 37 figures

BL Lacertae (BL Lac) objects are the prominent members of the third {\it Fermi} Large Area Telescope catalog of $\gamma$-ray sources. Half of the BL Lac population ($\sim$ 300) lack redshift measurements, which is due to the absence of lines in their optical spectrum, thereby making it difficult to utilize spectroscopic methods. Our photometric drop-out technique can be used to establish the redshift for a fraction of these sources. This work employed 6 filters mounted on the $Swift$-UVOT and 4 optical filters on two telescopes, the 0.65 m SARA-CTIO in Chile and 1.0 m SARA-ORM in the Canary Islands, Spain. A sample of 15 sources was extracted from the $Swift$ archival data for which 6 filter UVOT observations were conducted. By complementing the {\it Swift} observations with the SARA ones, we were able to discover two high redshift sources: 3FGL J1155.4-3417 and 3FGL J1156.7-2250 at $z=1.83^{+0.10}_{-0.13}$ and $z=1.73^{+0.11}_{-0.19}$ , respectively, resulting from the dropouts in the powerlaw template fits to these data. The discoveries add to the important (26 total) sample of high-redshift BL Lacs. While the sample of high-z BL Lacs is still rather small, these objects do not seem to fit well within known schemes of the blazar population and represent the best probes of the extragalactic background light., Comment: 10 pages, 6 figures (accepted by ApJ)

The collapsar model has proved highly successful in explaining the properties of long gamma-ray bursts (GRBs), with the most direct confirmation being the detection of a supernova (SN) coincident with the majority of nearby long GRBs. Within this model, a long GRB is produced by the core-collapse of a metal-poor, rapidly rotating, massive star. The detection of some long GRBs in metal-rich environments, and more fundamentally the three examples of long GRBs (GRB 060505, GRB 060614 and GRB 111005A) with no coincident SN detection down to very deep limits is in strong contention with theoretical expectations. In this paper we present MUSE observations of the host galaxy of GRB 111005A, which is the most recent and compelling example yet of a SN-less, long GRB. At z=0.01326, GRB 111005A is the third closest GRB ever detected, and second closest long duration GRB, enabling the nearby environment to be studied at a resolution of 270 pc. From the analysis of the MUSE data cube, we find GRB 111005A to have occurred within a metal-rich environment with little signs of ongoing star formation. Spectral analysis at the position of the GRB indicates the presence of an old stellar population (tau > 10 Myr), which limits the mass of the GRB progenitor to M_ZAMS<15 Msolar, in direct conflict with the collapsar model. Our deep limits on the presence of any SN emission combined with the environmental conditions at the position of GRB 111005A necessitate the exploration of a novel long GRB formation mechanism that is unrelated to massive stars., Comment: Now accepted by A&A. Manuscript replaced to match accepted version. Some additional discussion added, and velocity map of the host galaxy now included

Superluminous supernovae (SLSNe) are found predominantly in dwarf galaxies, indicating that their progenitors have a low metallicity. However, the most nearby SLSN to date, SN 2017egm, occurred in the spiral galaxy NGC 3191, which has a relatively high stellar mass and correspondingly high metallicity. In this paper, we present detailed analysis of the nearby environment of SN 2017egm using MaNGA IFU data, which provides spectral data on kiloparsec scales. From the velocity map we find no evidence that SN 2017egm occurred within some intervening satellite galaxy, and at the SN position most metallicity diagnostics yield a solar and above solar metallicity (12 + log (O/H) = 8.8-9.1). Additionally we measure a small H-alpha equivalent width (EW) at the SN position of just 34 Angs, which is one of the lowest EWs measured at any SLSN or Gamma-Ray Burst position, and indicative of the progenitor star being comparatively old. We also compare the observed properties of NGC 3191 with other SLSN host galaxies. The solar-metallicity environment at the position of SN 2017egm presents a challenge to our theoretical understanding, and our spatially resolved spectral analysis provides further constraints on the progenitors of SLSNe., Comment: Accepted version in ApJ Letter. Thank you for useful comments

Since the launch of the highly successful and ongoing Swift mission, the field of gamma-ray bursts (GRBs) has undergone a revolution. The arcsecond GRB localizations available within just a few minutes of the GRB alert has signified the continual sampling of the GRB evolution through the prompt to afterglow phases revealing unexpected flaring and plateau phases, the first detection of a kilonova coincident with a short GRB, and the identification of samples of low-luminosity, ultra-long, and highly dust extinguished GRBs. The increased numbers of GRB afterglows, GRB-supernova detections, redshifts, and host galaxy associations has greatly improved our understanding of what produces and powers these immense, cosmological explosions. Nevertheless, more high quality data often also reveal greater complexity. In this review, I summarize some of the milestones made in GRB research during the Swift era, and how previous widely accepted theoretical models have had to adapt to accommodate the new wealth of observational data., Comment: Article replaced to match published version

Short-duration GRBs with extended emission form a subclass of short GRBs, comprising about 15% of the short-duration sample. Afterglow detections of short GRBs are also rare (about 30%) due to their smaller luminosity. We present a multi-band data set of the short burst with extended emission GRB 150424A, comprising of GROND observations, complemented with data from Swift/UVOT, Swift/XRT, HST, Keck/LRIS and data points from the literature. The GRB 150424A afterglow shows an extended plateau phase, lasting about 8hrs. The analysis of this unique GRB afterglow might shed light on the understanding of afterglow plateau emission, the nature of which is still under debate. We present a phenomenological analysis by applying fireball closure relations, and interpret the findings in the context of the fireball model. We discuss the plausibility of a magnetar as a central engine, being responsible for additional and prolonged energy injection into the fireball. We find convincing evidence for energy injection into the afterglow of GRB 150424A. We find that a magnetar spin down as source for a prolonged energy injection requires that at least 4% of the spin-down energy is converted to radiation., Comment: 13 pages, 8 figures

Afterglows of Gamma-Ray Bursts (GRBs) are simple in the most basic model, but can show many complex features. The ultra-long duration GRB 111209A, one of the longest GRBs ever detected, also has the best-monitored afterglow in this rare class of GRBs. We want to address the question whether GRB 111209A was a special event beyond its extreme duration alone, and whether it is a classical GRB or another kind of high-energy transient. The afterglow may yield significant clues. We present afterglow photometry obtained in seven bands with the GROND imager as well as in further seven bands with the UVOT telescope on-board the Neil Gehrels Swift Observatory. The light curve is analysed by multi-band modelling and joint fitting with power-laws and broken power-laws, and we use the contemporaneous GROND data to study the evolution of the spectral energy distribution. We compare the optical afterglow to a large ensemble we have analysed in earlier works, and especially to that of another ultra-long event, GRB 130925A. We furthermore undertake a photometric study of the host galaxy. We find a strong, chromatic rebrightening event at approx 0.8 days after the GRB, during which the spectral slope becomes redder. After this, the light curve decays achromatically, with evidence for a break at about 9 days after the trigger. The afterglow luminosity is found to not be exceptional. We find that a double-jet model is able to explain the chromatic rebrightening. The afterglow features have been detected in other events and are not unique. The duration aside, the GRB prompt emission and afterglow parameters of GRB 111209A are in agreement with the known distributions for these parameters. While the central engine of this event may differ from that of classical GRBs, there are multiple lines of evidence pointing to GRB 111209A resulting from the core-collapse of a massive star with a stripped envelope., Comment: A&A, in press, expanded discussion on energetics, and updated otherwise. 32 pages, 17 pages main paper, 2 pages Appendix, 11 pages data tables. arXiv admin note: substantial text overlap with arXiv:1606.06791

We reveal multiple components of an interacting galaxy system at $z\approx3.35$ through a detailed analysis of the exquisite high-resolution Keck/HIRES spectrum of the afterglow of a gamma-ray burst (GRB). Through Voigt-profile fitting of absorption lines from the Lyman-series, we constrain the neutral hydrogen column density to $N_{\mathrm{HI}} \leq 10^{18.35}$ cm$^{-2}$ for the densest of four distinct systems at the host redshift of GRB~080810, among the lowest $N_{\mathrm{HI}}$ ever observed in a GRB host, despite the line of sight passing within a projected 5 kpc of the galaxy centres. By detailed analysis of the corresponding metal absorption lines, we derive chemical, ionic and kinematic properties of the individual absorbing systems, and thus build a picture of the host as a whole. Striking differences between the systems imply that the line of sight passes through several phases of gas: the star-forming regions of the GRB host; enriched material in the form of a galactic outflow; the hot and ionised halo of a second, interacting galaxy falling towards the host at a line-of-sight velocity of 700 km s$^{-1}$; and a cool, metal-poor cloud which may represent one of the best candidates yet for the inflow of metal-poor gas from the intergalactic medium., Comment: 18 pages, 15 figures. A&A, Accepted. Updated to accepted version

The most powerful sources among the blazar family are MeV blazars. Often detected at $z>2$, they usually display high X- and \gm-ray luminosities, larger-than-average jet powers and black hole masses $\gtrsim 10^9 M_{\odot}$. In the present work we perform a multiwavelength study of three high redshift blazars: 3FGL J0325.5+2223 ($z=2.06$), 3FGL J0449.0+1121 ($z= 2.15$), and 3FGL J0453.2$-$2808 ($z=2.56$), analysing quasi simultaneous data from GROND, \swift-UVOT and XRT, \nustar, and \fermi-LAT. Our main focus is on the hard X-ray band recently unveiled by \nustar~(3$-$79 keV) where these objects show a hard spectrum which enables us to constrain the inverse Compton peak and the jet power. We found that all three targets resemble the most powerful blazars, with the synchrotron peak located in the sub-millimeter range and the inverse Compton peak in the MeV range, and therefore belong to the MeV blazar class. Using a simple one zone leptonic emission model to reproduce the spectral energy distributions, we conclude that a simple combination of synchrotron and accretion disk emission reproduces the infrared-optical spectra while the X-ray to \gm-ray part is well reproduced by the inverse Compton scattering of low energy photons supplied by the broad line region. The black hole masses for each of the three sources are calculated to be $\gtrsim 4 \times 10^{8} M_{\odot}$. The three studied sources have jet power at the level of, or beyond, the accretion luminosity., Comment: 4 figures, 3 tables, accepted for publication in ApJ

Spatially-resolved spectroscopy of the environments of explosive transients carries detailed information about the physical properties of the stellar population that gave rise to the explosion, and thus the progenitor itself. Here, we present new observations of ESO184-G82, the galaxy hosting the archetype of the $\gamma$-ray burst/supernova connection, GRB 980425/SN 1998bw, obtained with the integral-field spectrograph MUSE mounted at the Very Large Telescope. These observations have yielded detailed maps of emission-line strength for various nebular lines, as well as physical parameters such as dust extinction, stellar age, and oxygen abundance on spatial scales of 160 pc. The immediate environment of GRB 980425 is young (5-8 Myr) and consistent with a mildly-extinguished ($A_V\sim0.1\ \mathrm{mag}$) progenitor of zero-age main-sequence mass between 25 $M_{\odot}$ and 40 $M_{\odot}$ and oxygen abundance 12+log(O/H)~8.2 ($Z\sim0.3\ {Z}_\odot$), which is slightly lower than the one of an integrated measurement of the galaxy (12+log(O/H)~8.3) and a prominent nearby HII region (12+log(O/H)~8.4). This region is significantly younger than the explosion site, and we argue that a scenario in which the GRB progenitor formed in this environment and was subsequently ejected appears very unlikely. We show that empirical strong-line methods based on [OIII] and/or [NII] are inadequate to produce accurate maps of oxygen abundance at the level of detail of our MUSE observation as these methods strongly depend on the ionization state of the gas. The metallicity gradient in ESO184-G82 is -0.06 dex kpc$^{-1}$, indicating that the typical offsets of at most few kpc for cosmological GRBs on average have a small impact on oxygen abundance measurements at higher redshift., Comment: 13 pages, 12 figures, A&A in press