EDP Sciences
Free Access
Issue
A&A
Volume 602, June 2017
Article Number A85
Number of page(s) 13
Section Extragalactic astronomy
DOI https://doi.org/10.1051/0004-6361/201630268
Published online 20 June 2017

© ESO, 2017

1. Introduction

Line emission from recombination of ionized hydrogen or from the decay of collisionally excited states of metal ions is a fundamental tracer of the physical conditions in H ii regions. The absolute and relative intensities of these transitions crucially depend on the ionizing source, electron density in the plasma, ionization state of the elements, and gas-phase abundances (Osterbrock 1989). This makes emission-line spectra of astronomical sources one of the most elementary diagnostics of galaxy formation and evolution (e.g., Tremonti et al. 2004; Erb et al. 2006; Förster Schreiber et al. 2009). The total intensity of the hydrogen recombination lines, for example, is proportional the number of O-type stars and thus traces the star formation rate at timescales of ~10 Myr (e.g., Kennicutt 1998). The continuum emission at the wavelength of Hα in turn originates from B- or A-type stars, which makes the Hα equivalent width (EW) a good proxy for the age of the stellar population (Leitherer et al. 1999; Levesque & Leitherer 2013).

Metal abundances have been measured through ratios of prominent emission lines from ions such as O+, O2+, N+, S+, and/or recombination lines of hydrogen (Pagel et al. 1979; Alloin et al. 1979). Given their fundamental importance in galaxy evolution and cosmology, these abundance determinations through nebular emission lines have been the focus of a large body of literature (e.g., Kobulnicky & Kewley 2004; Pilyugin & Thuan 2005; Stasińska 2006; Izotov et al. 2006; Kewley & Ellison 2008).

It is thus immediately clear that an emission-line spectrum of cosmological sources carries detailed information about the underlying stellar population and thus has been used to infer properties not only of galaxies but also of the progenitors of explosive transients. Global (e.g., Prieto et al. 2008; Li et al. 2011) or local (e.g., Anderson et al. 2010; Modjaz et al. 2011; Leloudas et al. 2011) properties of nearby supernova (SN) hosts and cosmological γ-ray bursts (GRBs; e.g., Wiersema et al. 2007; Krühler et al. 2012a; Graham & Fruchter 2013) or superluminous supernovae (e.g., Chen et al. 2013; Lunnan et al. 2014; Leloudas et al. 2015; Perley et al. 2016) have likewise been used to compare progenitor models with the expected environments.

A fundamental assumption of all these studies is the hypothesis that there is a tight relation between the two primary observables, gas-phase oxygen abundance and age of H ii regions, and the progenitor properties of metallicity and lifetime and thus initial mass. Clearly, this link is most robust when coming from an analysis of the co-spatial stellar population. Integral field spectroscopy (IFS) with high angular resolution is thus arguably the most comprehensive way of studying the environments of explosive transients. Low-redshift galaxies hosting supernovae (SNe), for example, are hence ideal targets for state-of-the-art integral field units (IFUs; e.g., Kuncarayakti et al. 2013a,b; Galbany et al. 2014).

In contrast to SNe, the vast cosmological distances of GRBs (e.g., Fynbo et al. 2009; Tanvir et al. 2009; Krühler et al. 2012b) have always posed serious limitations on using IFS for GRB-selected galaxies. Only a few GRBs are close enough such that the spatial resolution achievable with modern ground-based instrumentation yields constraints on spatial scales better than a kpc. Spatially resolved spectroscopy has therefore only been obtained for a handful of nearby GRB hosts using long-slit spectroscopy (e.g., Thöne et al. 2008; Levesque et al. 2011; Schady et al. 2015) or the previous generation of IFUs (Christensen et al. 2008; Thöne et al. 2014).

The focus of this article is on our new observations of the poster child of the connection between GRBs and SNe, GRB 980425/SN 1998bw, using the Multi-Unit Spectroscopic Explorer (MUSE; Bacon et al. 2010). GRB 980425 is the closest GRB yet discovered; the GRB (e.g. Galama et al. 1998; Kulkarni et al. 1998), SN (e.g., Iwamoto et al. 1998; Patat et al. 2001; Maeda et al. 2006), and host galaxy (e.g., Fynbo et al. 2000; Sollerman et al. 2005; Hammer et al. 2006; Michałowski et al. 2009, 2014, 2016; Le Floc’h et al. 2012; Arabsalmani et al. 2015) are extensively discussed in existing literature.

Compared to the bulk of cosmological GRBs, GRB 980425 is rather peculiar: the isotropic-equivalent release in γ-rays of GRB 980425 was ~1048 erg (Galama et al. 1998), which is a factor of ten lower than other local, low-luminosity GRBs, or around five orders of magnitudes less than conventional, higher redshift GRBs (Xu et al. 2013). No bright multiwavelength afterglow was observed for GRB 980425 despite its proximity. However, the associated SN without hydrogen or helium in its spectrum, broad metal absorption lines, and high luminosity has proven to be typical of GRB-related SNe in general (Hjorth & Bloom 2012).

The host galaxy of GRB 980425/SN 1998bw, ESO184-G82 (Lauberts & Valentijn 1989), is a barred spiral dwarf galaxy (Fynbo et al. 2000) seen nearly face on and shown in Fig. 1. It has a visible extent of approximately 67′′ by 57′′ (12 × 10 kpc) at the B = 26.5 mag isophote (Sollerman et al. 2005). Its brightness, luminosity, and stellar mass are B = 14.94 mag, MB = − 17.65 mag or L = 0.05 L, and log  (M/M) = 8.7, respectively (Sollerman et al. 2005; Michałowski et al. 2014). SN 1998bw exploded in an H ii region 12′′ distant (2 kpc projected) from its center and 860 pc to the southeast of a young star-forming region that displays signatures of Wolf-Rayet (WR) stars in its spectrum (Hammer et al. 2006), the so-called Wolf-Rayet region (Fig. 1).

thumbnail Fig. 1

False-color composite from reconstructed VRI-band images from the MUSE data cube. The image spans approximately 55′′ by 55′′, or 10 by 10 kpc at z = 0.0086.

Open with DEXTER

Despite the large set of recent literature on GRB 980425, SN 1998bw and its host mentioned above, we summarize our new data and conclusions here for three main reasons: first, the unique combination of spatial resolution and sensitivity of MUSE yields detailed maps of emission-line strength, dust reddening, and oxygen abundance, which help us to clarify some of the ambiguities around SN 1998bw and its host from previous works. Second, these maps provide the tightest and most accurate constraints on the immediate environment and underlying stellar population of SN 1998bw yet available and thus allow us to infer the progenitor properties of the GRB from its parent stellar population. And last, our new data offer an informative example of spatially resolved oxygen-abundance measurements in star-forming galaxies through strong-line diagnostics and their dependence on other physical conditions in the interstellar medium (ISM).

Throughout the paper, we adopt a flat ΛCDM cosmology with Planck parameters (H0 = 67.3 km s-1 Mpc-1, Ωm = 0.315, ΩΛ = 0.685; Planck Collaboration XVI 2014), a Chabrier (2003) initial mass function (IMF), solar abundances from Asplund et al. (2009), and report errors at the 1σ confidence level.

2. Observations

We observed ESO184-G82 (z = 0.0086, or DL = 37 Mpc) using MUSE (Bacon et al. 2010) at the ESO Very Large Telescope (VLT) during the two clear nights starting on 2015-05-14 and 2015-05-15 in a classical observing run from Paranal. On each night, we obtained four dithered exposures of 450 s integration each, totaling 3600 s on source. The on-target frames were supplemented by an offset pointing to blank sky for 200 s. For absolute flux calibration, the spectrophotometric standard LTT3218 was observed at the beginning of each night. The full width half maximum of the stellar point spread function, which defines our spatial resolution, is between 0.′′9 (at 9000 Å) and 1.′′1 (at 5000 Å) in the MUSE data.

The MUSE instrument is a state-of-the-art integral field spectrograph that splits the light into 24 identical subunits. In the wide-field mode, each of these sub-IFUs disperses a 60′′ × 2.5′′ patch of the sky onto a single CCD. In this way, MUSE covers a continuous sky region of 60′′ × 60′′ in the wavelength range between 4750 Å and 9300 Å when operated in its nominal configuration. With its excellent total throughput, small spaxel size (0.′′2 × 0.′′2), and decent resolving power (1800 <R< 3600 increasing from blue to red wavelengths), MUSE offers an unprecedented combination of sensitivity, spatial resolution, and field of view for IFUs (Bacon et al. 2010).

3. Data reduction

We reduced the MUSE data with the pipeline supplied through ESO1 in its version 1.2.1 (Weilbacher et al. 2014), which applies corrections for bias level, flat-fields, illumination level, and geometric distortions. The pipeline also performs the wavelength calibration using daytime arc-lamp frames, which is subsequently refined by skylines in the science data. The sky background was subtracted using an offset pointing to blank sky and making use of algorithms from the Zurich Atmospheric Package (Soto et al. 2016). The exposures from the two different nights were then corrected for slight pointing offsets between night one and two, stacked using variance weighting, and dereddened based on the Galactic foreground EBV = 0.05 mag (Schlafly & Finkbeiner 2011) assuming an average Milky Way extinction law (Pei 1992) and RV = 3.08.

The spectrum of the star at , Decl(J2000) = − 52°50′21.′′1 was used to correct for telluric absorption via molecfit (Smette et al. 2015). By fitting the three prominent telluric absorption bands within the MUSE wavelength coverage (centered around 6870 Å, 7600 Å, and 7630 Å) with a physical model of the atmospheric molecular oxygen and water vapor content, we derived the telluric absorption for this single star, which we then subsequently applied to all spaxels.

The final data cube has slight astrometric offsets, which we corrected by tying the position of stars in the field of MUSE to coordinates from a reference image taken on 2000-10-25 with the SOFI imager on the ESO New Technology Telescope. We then measured the position of the SN in the reference frame, mapping it onto the MUSE cube with an accuracy of around 50 mas. Figure 1 shows a false-color image reconstructed from the MUSE cube where the position of SN 1998bw derived from the reference image is indicated.

Similarly, we use existing photometry to corroborate our flux calibration through the V-, RC- and IC-band magnitudes of star 1 of Clocchiatti et al. (2011) and synthetic photometry from the MUSE data cube, yielding differences of ΔV = 0.05 ± 0.03 mag, ΔRC = 0.05 ± 0.06 mag, and ΔIC = 0.00 ± 0.05 mag. After applying a linear fit in wavelength to these correction terms, we can accurately reproduce the optical colors of the host galaxy (Sollerman et al. 2005) to better than 0.02 mag.

4. Analysis and discussion

4.1. Separating gas-phase and stellar component

thumbnail Fig. 2

Separating stellar and gas-phase components at the wavelength around Hβ and [O iii](λλ4959, 5007). Top: original spectrum (black) and fitted stellar component (red). Middle: zoom-in to the continuum. Bottom: resulting spectrum of the gas-phase contribution (blue) together with the error spectrum (grey).

Open with DEXTER

As we are primarily interested in the absolute and relative strengths of the nebular lines and thus the ionized gas component of the galaxy, we needed to remove the stellar Balmer absorption for accurate line flux measurements, in particular for the Hβ line (Fig. 2). The strength of the stellar absorption is primarily a function of the age of the underlying stellar population. It thus depends on the position within a galaxy and needs to be accurately modeled for reliable constraints on the Balmer decrement, which we use to measure dust reddening maps.

We separated the stellar and gas-phase components of the galaxy by fitting a linear superposition of template spectra, based on the Bruzual & Charlot (2003) models, to the MUSE data. We divided the full field of view into regions with a size of 0.′′4 × 0.′′4 (or 2 × 2 spaxels) and extracted spectra for each of the ~22 000 subregions. These spectra were then fit with stellar population models using starlight (Cid Fernandes et al. 2005, 2009) following methods that we previously applied to MUSE data (Galbany et al. 2016a; Kuncarayakti et al. 2016; Prieto et al. 2016). The 2 × 2 co-adding effectively increases the signal-to-noise ratio (S/N) at the expense of spatial resolution for the stellar properties, but is necessary to robustly perform an automated fit in particular in the fainter regions of the galaxy. We then linearly scale the best-fit stellar template to the intensity in single spaxels. Subtracting this stellar component from the original data leaves us with the spectral contribution of the gas-phase only (Fig. 2).

After having isolated the stellar component from the observed spectrum, it is trivial to produce maps of line flux (integral over the nebular emission line), continuum (average of the stellar component), and equivalent width (flux divided by continuum).

Plotting the characteristic emission-line flux ratios of [O iii]/Hβ versus [N ii]/Hα for individual spaxels in the BPT diagram (Baldwin et al. 1981) allows us to immediately discard active galactic nuclei or shocks as the ionization source and to confirm that the radiation from massive stars is the origin of the nebular lines in ESO184-G82 (Fig. 3).

thumbnail Fig. 3

Spaxel BPT diagram for ESO184-G82. The solid line denotes the ridge line of SDSS galaxies (Brinchmann et al. 2008), whereas the dotted line represents the z = 0 classification line between star formation and AGN ionization (Kewley et al. 2013). The explosion site, a galaxy average (including all spaxels), and an H ii-region average (including only spaxels with EWHα> 10 Å) are indicated.

Open with DEXTER

4.2. Equivalent width maps and stellar population ages

The Hα equivalent width map, which is a rather direct proxy of stellar population age, in particular for young stellar populations, is shown in Fig. 4. The spaxel closest to the SN position has an Hα equivalent width of EWHα = 98 ± 4 Å. The spaxels within a radius of 70 pc yield EWHα = 92 ± 15 Å.

thumbnail Fig. 4

Reconstructed images from the MUSE data cube. Each panel shows the host of GRB 980425 in a different nebular line. Top left: Hα. Top right: [O iii](λ5007). Bottom left: [S ii](λ6731). Bottom right: [S iii](λ9070). All panels are approximately 55′′ by 55′′, or 10 kpc by 10 kpc at the GRB redshift. The effective spatial resolution is given by the point spread function (PSF) indicated in the lower left corner of each image with a full width half maximum (FWHM) of approximately 0.′′9 or 160 pc.

Open with DEXTER

Assuming a single stellar population from an instantaneous starburst, this EWHα corresponds to stellar population ages between 5 Myr and 8 Myr2 from various models at metallicities of Z = 0.2Z (see, e.g., Levesque & Leitherer 2013; Kuncarayakti et al. 2016, and references therein). The relatively large range in age is almost entirely due to the spread from different stellar evolution models or initial mass functions. These ages corresponds to lifetimes of stars with zero-age main sequence (ZAMS) masses between approximately 25 and 40 M(Fagotto et al. 1994; Meynet & Maeder 2005). This is consistent with the ejected oxygen mass of the SN (MO ~ 5 − 6 M evolved from a MZAMS ~ 30 − 40Mstar) derived through modeling the SN 1998bw nebular spectra (Mazzali et al. 2001; Maeda et al. 2006).

These considerations are only valid, of course, if the progenitor was born where it exploded, and was not ejected from the nearby WR region (Hammer et al. 2006). However, this region is so young that timing arguments make an ejection scenario quite contrived: Very high EW values of Hα and nebular transitions are observed in the center of the WR region (EWHα> 900 Å and EW[ O iii ] (λ5007)> 850 Å)3. Together with the detection of strong He i(λ4922) with an EWHeI(λ4922) = 3.5 ± 0.2 Å, they ascertain population ages younger than 3 Myr (see also Appendix A) in instantaneous starburst models, or MZAMS ≳ 60M(see, e.g., Thöne et al. 2015, and references therein for a similar case). The discrepancy with respect to the progenitor mass from SN modeling (e.g., Mazzali et al. 2001; Maeda et al. 2006) then suggests that the GRB progenitor was indeed not born in the WR region, but rather formed in situ.

To travel to the explosion site in less time than the age of the WR region (<3 Myr), the progenitor would require very high peculiar velocities v of v> 260 km s-1. This is again a strict lower limit, as projection effects further increase the necessary velocities. Scenarios that are believed to give rise to these kinds of massive runaway stars are dynamical few-body encounters or binary supernovae, but both seem unfeasible here. A dynamical ejection produces hyper-velocity stars in only very rare and extreme cases (Hoogerwerf et al. 2001; Perets & Šubr 2012), and the probability of potential GRB progenitors (MZAMS ≳ 20 M) receiving velocity kicks of >200 km s-1 from a companion SN is also practically zero (Eldridge et al. 2011). This makes a binary supernova origin highly implausible as there simply would not be enough time to evolve and explode the primary and eject the secondary to a distance 860 pc. Also the fraction of stars ejected by dynamical encounters is of course a function of the elapsed time after starburst, and reaches only 0.01 or 0.03 at 1 Myr or 3 Myr at MZAMS ~ 35 M (Banerjee et al. 2012), again leaving little time for the ejected star to travel as far as 860 pc (or further).

Given the presence of massive stars in the vicinity of the SN position (Fynbo et al. 2000), the substantial level of recent star formation, as evidenced by high EWs of nebular lines at the SN position (Fig. 4), and the consistency between MZAMS derived from the age of the stellar population and the SN 1998bw nebular spectra, we see no compelling reason to invoke an artificial ejection from the nearby H ii region to explain the GRB location within its host.

4.3. Dust distribution

thumbnail Fig. 5

Dust-reddening distribution in ESO184-G82 as measured through the Balmer decrement. We only show spaxels in which Hβ is detected with a significance of at least 3σ. The image spans 34′′ by 38′′, or 6.1 kpc by 6.8 kpc. The circle denoting the position of SN 1998bw has a diameter of 160 pc.

Open with DEXTER

Because SN 1998bw has an exquisite photometric and spectroscopic data coverage, it is widely used as a comparison object. Deriving its intrinsic properties such as luminosity or color, however, requires an accurate knowledge of the absorbing dust column in its host, which we infer here through Balmer lines. We convert line fluxes of Hα and Hβ into a map of color excess EBV as shown in Fig. 5 using Eq. (5) of Krühler et al. (2015). This procedure assumes standard ratios of the Balmer lines from Osterbrock (1989) at a temperature or 104 K and electron density ne ~ 100 cm-3. These assumptions are broadly consistent with the values measured for ESO184-G82 (Table 1). Our results depend only marginally on our choice of an average Milky Way extinction law (Pei 1992) with RV = 3.08, as the difference between reddening laws in the Local Group is small in the wavelength range probed by Hα and Hβ.

Figure 5 shows very little dust in general in ESO184-G82, and also only minor evidence for dust at the spaxel closest to the actual GRB or SN position (EBV = 0.04 ± 0.03 mag or AV = 0.13 ± 0.09 mag) and its immediate environment ( in the nine spaxels closest to the GRB position). The only locations where we observe significant dust reddening are the centers of H ii regions as exemplified by the WR region (Fig. 5). These areas show a centrally symmetric substructure in dust extinction that is decreasing from the inside out and peaking at EBV ~ 0.25 mag or AV ~ 0.8 mag. A galaxy-integrated spectrum yields EBV = 0.05 ± 0.02 mag or AV = 0.15 ± 0.06 mag, which is in remarkable agreement with the average optical depth τV derived from modeling the UV-to-radio spectral energy distribution (Michałowski et al. 2014).

These low dust reddening values are somewhat in tension with results from optical spectroscopy in previous works (Hammer et al. 2006; Christensen et al. 2008). While we can reproduce our MUSE values with a re-reduction of the archival data used in Hammer et al. (2006), as shown in Appendix A, the origin of the mismatch to the spatially resolved data of Christensen et al. (2008) remains unclear4.

The reliable spectrophotometric calibration of the MUSE data (Sect. 3), our accurate stellar continuum modeling (Sect. 4.1) and the consistency with the spatially resolved dust properties from long-wavelength data (Michałowski et al. 2014) leads us to believe that our new data now resolve the apparent conflict between an unexpectedly large reddening at the SN position derived in previous works and the SN itself, which did not show strong evidence of dust obscuration (e.g., Iwamoto et al. 1998; Patat et al. 2001). These data provide further confidence in using SN 1998bw as an only marginally reddened SN template for comparison to other events (e.g., Zeh et al. 2004; Schulze et al. 2014; Kann et al. 2017, and references therein).

4.4. Metallicity diagnostics

thumbnail Fig. 6

Face-to-face comparison between the ESO184-G82 maps of O3N2, often used for abundance determinations and [S iii]/[S ii], a tracer of the ionization state of the hot gas. Only spaxels with S/N> 3 are shown. Image dimensions are similar to Fig. 5.

Open with DEXTER

4.4.1. Initial considerations

Metal abundances of H ii regions are a central observable to study cosmochemical evolution, and a large set of literature is devoted to the various possibilities, their advantages, and perils to infer abundances from H ii-region spectra (e.g., Pagel et al. 1979; McGaugh 1991; Pilyugin & Thuan 2005; Kewley & Ellison 2008). Very briefly, the most common methods to deduce chemical abundances, and from those, the abundance of oxygen (traditionally expressed in 12 + log (O / H)) make use of either photoionization models (e.g., Evans & Dopita 1985; Dopita et al. 2000; Kewley & Dopita 2002) or empirical correlations between certain strong-line ratios and oxygen abundances derived through electron temperatures Te from collisionally excited lines (CELs; e.g., Pettini & Pagel 2004; Marino et al. 2013). Commonly used ratios are, for example, [N ii]/[O ii], [O iii]/[N ii], [N ii]/Hα, or R23 = ([O ii] + [O iii])/Hβ, which have been (re)calibrated numerous times against different samples of Te-based abundances or photoionization models, yielding a large set of different calibrators in the literature (e.g., Kewley & Dopita 2002; Kobulnicky & Kewley 2004; Pilyugin & Thuan 2005; Nagao et al. 2006; Maiolino et al. 2008).

The most fundamental problem in using and interpreting the oxygen abundances derived in this way is that different methods are only very rarely consistent (e.g., Kewley & Ellison 2008), which gives rise to the abundance determination problem (Peimbert 1967). Methods based on temperature-sensitive CELs typically show abundances that are lower by 0.20.4 dex with respect to photoionization-based methods or abundances derived with temperatures from recombination lines (e.g., López-Sánchez et al. 2012, and references therein), in particular in the high-metallicity region. A possible solution to the abundance determination problem are small-scale temperature fluctuations (e.g., Peimbert 2003; Esteban et al. 2004) or an electron population that is distributed somewhat differently than in thermal Maxwell-Boltzmann equilibrium (Nicholls et al. 2012; López-Sánchez et al. 2012). But until these discrepancies are fully resolved, element abundances from emission lines remain the subject of large controversy.

A second, independent problem relates to the observational difficulties in robustly measuring emission-line fluxes for lines in different wavelength ranges for faint, high-redshift galaxies. As a consequence of various observational constraints, the available data are typically limited to a handful of strong lines. This is similarly true for our observations here, as the MUSE data do not cover the strong [O ii](λλ3726, 3729) doublet nor [O iii](λ4363), which is one of the most commonly used, temperature-sensitive CEL. A very popular emission-line diagnostic in the literature has thus been the logarithm of the ratio of [O iii] /Hβ to [N ii]/Hαor short O3N2 (e.g., Pettini & Pagel 2004; Marino et al. 2013) because of its independence on dust reddening and the relative observational ease with which it can be measured even at z ~ 2.

4.4.2. Specific problems of empirical metallicity diagnostics

From the very different ionization potentials (IPs) of N and O+ (14.5 eV versus 35.1 eV), it is directly clear that O3N2 should also carry a strong dependence on the ionization parameter in addition to its inverse proportionality to oxygen abundance (e.g., Alloin et al. 1979; Ho et al. 2015). In Fig. 6, we plot the O3N2 map of ESO184-G82, which immediately translates into a map of oxygen abundance (12 + log (O / H) would be lowest where O3N2 is highest) in common strong-line diagnostics (Pettini & Pagel 2004).

When inspecting the left part of Fig. 6, however, it becomes apparent that O3N2-based oxygen abundances (and similarly for those from N2) produce abundance maps that are hard to understand in a physical context: O3N2 varies significantly on kpc scales leading to an unexpected5 chemically inhomogeneous structure within individual H ii regions with their central abundances up to 0.3 dex lower than their outer edges. However, we believe that the significant gradients in O3N2/N2 observed in most of the H ii regions are unlikely to be due to a genuine variation in oxygen abundance, but are more likely to result from a changing ionization parameter. We explore this hypothesis further in the following sections.

4.4.3. Ionization map

Our MUSE data are of sufficient depth and quality to test how strongly O3N2 is affected by ionization empirically through the ratio of [S iii] (IP = 23.3 eV) to [S ii] (IP = 10.3 eV); this ratio is widely considered as one of the best tracers of the ionization parameter (Diaz et al. 1991) as it shows, in contrast to [O iii]/[O ii], only very little dependence on abundance itself (Kewley & Dopita 2002; Dors et al. 2011). The resulting map6 is shown in Fig. 6 and clearly highlights that the center of the H ii regions have the highest values of [S iii]/[S ii], and thus ionization parameter.

This adds further support to our initial conjecture and ascertains that the radially symmetric structure of individual H ii regions in the O3N2 map is due to an increase of the ionization parameter toward the center and not a decrease in chemical abundance. Simple O3N2 (or N2-based) diagnostic line ratios are thus inadequate to produce accurate maps of oxygen abundance at the level of detail of our MUSE data (see also Appendix B). Comparing the specific measurements of oxygen abundance in Table 1, it becomes clear that the O3N2-based diagnostic systematically underestimates 12 + log (O / H) for high ionization parameters as observed in the WR region, and overestimates 12 + log (O / H) for regions of lower ionization such as the SN explosion site.

4.4.4. Metallicity map

thumbnail Fig. 7

Map of 12 + log (O / H) as obtained through the [S ii] and [N ii]-based method of Dopita et al. (2016). Only spaxels with S/N> 3 are shown. Image dimensions are 34′′ by 38′′ or 6.1 kpc by 6.8 kpc, similar to Fig. 5. Methods to estimate 12 + log (O / H) based on electron temperatures return values that are 0.1–0.2 dex higher.

Open with DEXTER

After rejecting common empirical methods using O3N2 or N2 as accurate metallicity tracers due to their ionization dependence, we turn to diagnostics based on photoionization models. Unfortunately, most of the previous strong-line methods rely in one way or another on the strong [O ii](λλ3726,3729) doublet (Kewley & Dopita 2002), which is not directly available to us here. Also [S iii]/[S ii] is a good ionization tracer, but [S iii] is relatively faint and not detected in most of our spaxels (see Fig. 6).

A recently published method based on photoionization modeling (Dopita et al. 2016) seems to perfectly fit to our data; it relies solely on Hα, [N ii], and [S ii], which are all strong and well within the wavelength range of MUSE. The method, shortened as D16 in the following, is introduced as “effectively independent of both ionization parameter and ISM pressure” (Dopita et al. 2016); Fig. 7 shows the respective map of 12 + log (O / H) in regions where all necessary lines are detected at a S/N of at least 3.

Clearly, the strong abundance gradient over individual H ii regions, as would have been deduced from O3N2, is not observed in this diagnostic. Instead, the oxygen abundance map shows a relatively smooth behavior with a decreasing overall metallicity from the center of the galaxy toward the outside (see also Sect. 4.5). The spaxel abundance at the SN position is 12 + log (O / H) = 8.00 or 0.20 Z. The immediate environment is consistent with this value and homogeneous: spaxels within a radius of 70 pc to the SN position yield 12 + log (O / H) = 8.04 ± 0.06. The WR region displays a somewhat higher oxygen abundance with 12 + log (O / H) = 8.22 ± 0.03 or 0.34 ± 0.03 Z around the peak of the Hα emission.

It is of course reasonable to ask now whether the new D16 diagnostic provides more reliable constraints on oxygen abundance than previous methods given the significant differences that exists between all of them (Bian et al. 2017). In addition, for low-mass galaxies as is the case here, this diagnostic seems to return lower oxygen abundances than previous methods (Kashino et al. 2016). To elaborate further on the [S ii]-based diagnostic, we reproduce the combined Eqs. (1) and (2) from Dopita et al. (2016), i.e., (1)where [N ii] is the flux in the [N ii](λ6484) line and [S ii] the flux in the [S ii](λλ6717,6731) doublet. The primary observable is thus the nitrogen to sulfur ratio, which is a tracer of the nitrogen to oxygen ratio7. Because nitrogen is also produced in intermediate-mass stars, N/O starts to depend on 12 + log (O / H) above 12 + log (O / H) ~ 7.8 (e.g., Izotov & Thuan 1999; Pérez-Montero et al. 2013, 2016), and as expected, a N/O map via Amorín et al. (2010) is very similar to the map of oxygen abundance in the respective strong-line diagnostic.

The map of oxygen abundance then fundamentally relies on the N/O-to-O/H calibration, and it is in principle not impossible that the applied calibration is slightly offset, in particular in the low-metallicity region. However, the calibration sample for N/O-to-O/H in the metallicity range of interest is based on low-metallicity blue compact dwarf galaxies (Izotov & Thuan 1999), which are not dissimilar in physical properties to our galaxy. In addition, there could also be internal variations in N/O at a given 12 + log (O / H) or differences between the host of SN 1998bw to the calibration sample. For example, infall of primordial gas (e.g., Michałowski et al. 2015) would decrease 12 + log (O / H), but leave N/O unaffected (Kashino et al. 2016). Another point of concern would be an anomalously high N/O ratio for the SN region as claimed by Hammer et al. (2006). Both of these effects would lead us to overpredict the actual oxygen abundance in the SN region via the D16 diagnostic. The oxygen abundances derived through Eq. (1), however, are already lower than the values from temperature-sensitive collisionally excited lines (Sect. 4.4.5), which are sensitive to the actual oxygen abundance. We hence find no evidence for a systematic overestimation of the metallicity in the D16 scale at either the explosion site or the WR region.

4.4.5. Metallicities based on electron temperatures

Given the strong constraints that these measurements of oxygen abundance could imply for the SN 1998bw progenitor metallicity, we further seek to corroborate our earlier abundances through those from temperature-sensitive lines. Unfortunately, the brightest and most commonly used auroral line [O iii](λ4363) is not covered by the MUSE wavelength response and the other temperature-sensitive lines are too faint to be detected in most of the individual spaxels. However, we clearly detect [S iii](λ6312) both in the WR region (Fig. 8) and in integrated spectra around the SN position.

Following Nicholls et al. (2013), we derive an electron temperature Te using [S iii] in the central part of the WR region of around Te( [ S iii ] ) = (0.94 ± 0.02) × 104 K. Figure 8 contains maps of the crucial emission lines, i.e., nebular [S iii](λ9069) and auroral [S iii](λ6312); the resulting temperatures; and sulfur abundances. These values correspond to somewhat larger temperatures of [O iii] of around Te( [ O iii ] ) ~ 1.07 × 104 K (Izotov et al. 2006; Binette et al. 2012). The flux in the doublets of [O ii](λλ7320, 7330), [O iii](λλ4959, 5007) and Hβ then yields a central abundance of the WR region around 12 + log (O / H) ~ 8.3.

Similar values are obtained when using solar abundances to convert the measured sulfur to an oxygen abundance. These are slightly (0.10.2 dex) higher than those implied by the strong-line diagnostic from Sect. 4.4.4, but critically depend on the sulfur-to-oxygen temperature conversion or assumed abundance. They are thus subject to some systematic uncertainties. However, it is clear that the strong structure in O3N2 or N2 over the H ii region are not observed in electron temperatures.

Emission lines at the explosion site are substantially fainter and [S iii](λ6312) is only detected in a stacked spectrum extracted from three by three spaxels around the SN position, which yields Te( [ S iii ] ) = (1.24 ± 0.18) × 104 K. This is slightly higher than in the WR region, but with large uncertainties. The corresponding oxygen abundance derived in a similar manner as above is 12 + log (O / H) = 7.92 ± 0.27.

We further use here our reduction of the public archival VLT long-slit spectroscopy (Appendix A) as it covers both the WR region and SN position and extends below 4000 Å. From the well-detected [O iii](λ4363) line, we measure temperatures of Te( [ O iii ] ) = (1.05 ± 0.05) × 104 K and Te( [ O iii ] ) = (1.29 ± 0.18) × 104 K for the WR and SN region, which are broadly consistent with the estimates from MUSE through [S iii](λ6312). They imply oxygen abundances of 12 + log (O / H) = 8.36 ± 0.07 and 8.09 ± 0.15 for WR and SN region, respectively.

Our relatively low oxygen abundances at the location of explosion site, WR region, as well as the galaxy-integrated value (12 + log (O / H) ~ 8.3; Sect. 4.6) correspond to 0.2 to 0.5 times the solar value, and thus explain many of the long-wavelength properties of ESO184-G82 as typical for metal-poor dwarf galaxies without invoking a deficiency in molecular gas (Michałowski et al. 2016).

thumbnail Fig. 8

Electron temperatures (Te) from [S iii] and corresponding sulfur abundances in the WR region. Top left: flux map of nebular [S iii](λ9069). Top right: flux map of auroral [S iii](λ6312). Both flux maps are in units of 10-17erg cm-2 s-1. Bottom left: electron temperatures from [S iii]. Bottom right: 12 + log (S / H) as derived from Te. The solar abundance [O/S] is 1.57 (Asplund et al. 2009), so the 12 + log (S / H) scale corresponds to 12 + log (O / H) = 8.2 to 8.4. All panels are approximately 6′′ by 6′′, or 1 kpc by 1 kpc.

Open with DEXTER

4.5. Metallicity gradient

The metallicity of galaxies is often observed to decrease with the distance from their centers (e.g., Zaritsky et al. 1994; Sánchez et al. 2014), which is also seen in the hosts of SNe (Galbany et al. 2016b) and here (Fig. 7). These gradients are important to understand for spatially unresolved studies at high redshift, where positional offsets can be measured, but abundances are only derived in a galaxy-integrated manner.

Using a linear regression on the data in Fig. 9, the metallicity gradient is best fit with a slope of − 0.25 ± 0.02 dex/R25 in relative, or − 0.06 dex kpc-1 in physical scales, which is well in the range that was previously reported for galaxies of comparable stellar mass (Ho et al. 2015). Here, R25 is the radius at the B = 25 mag arcsec2 isophote. The O3N2-based diagnostics return slightly steeper, but generally compatible values within errors. The linear fit to the oxygen abundance data is a reasonable description of the data, except for the very center (deprojected distance <1 kpc), where the metallicity is seen to decrease more steeply. Limiting the fit range to this region, the slope in the central kpc is − 0.18 dex kpc-1 (see dashed line in Fig. 9), but it depends strongly on the metallicity diagnostic8.

Table 1

Physical properties inferred from MUSE IFU spectroscopy for the GRB 980425 host galaxy ESO184-G82.

A comparison between this metallicity gradient to the typical values of (projected) GRB distances from the galaxy centers9 of 1.3 kpc (Blanchard et al. 2016) does not provide strong reason to suggest that the average measurement of GRB host metallicities from spatially unresolved data is significantly skewed when compared to the GRB site metallicity. There is, however, a non-negligible fraction of GRBs at substantial distances to their hosts (10% of GRBs are located at offsets >3 kpc) where metallicity gradients might lead to overestimates of the GRB site abundance from unresolved spectra for single objects.

thumbnail Fig. 9

Metallicity gradient for ESO184-G82 in the D16 scale. Each data point is a spaxel-based measurement of oxygen abundance, where individual spaxels were binned with their direct neighbors to enhance the S/N if necessary.

Open with DEXTER

4.6. ESO184-G82 if seen at high redshift

GRB 980425 is the closest GRB since the discovery of GRB afterglows and provides an exceptional opportunity to measure the physical parameters of its host with high precision and high spatial resolution. Typically, information on the environments of GRBs or similar kinds of transients at higher redshift has only been obtained through galaxy-integrated measurements (Krühler et al. 2015; Japelj et al. 2016; Vergani et al. 2017), where it is not necessarily obvious how well the measured parameters actually correspond to properties at the GRB location. We thus compare the physical parameters derived from the explosion site spectrum and the WR region, extracted from the nine spaxels closest to the SN position and the WR region center, to the parameters obtained from a galaxy-integrated spectrum of ESO184-G82. The latter is a simple sum of all spaxels within a radius of 25′′ around the center of the galaxy, which should mimic how the galaxy would appear if it were unresolved and at high redshift.

First, we measure the total star formation rate of ESO184-G82 by adding up the Hα flux of all individual spaxels after the respective spaxel-based reddening correction, which yields SFR = 0.22 ± 0.02M yr-1. If we instead use the Hα flux from the galaxy-integrated spectrum and then apply a single reddening correction, we derive SFR = 0.17 ± 0.02M yr-1. Both values broadly agree with far-infrared, [O i], and [C ii]-based SFRs (Michałowski et al. 2014, 2016), and the narrowband Hα image from Sollerman et al. (2005) once their Salpeter IMF is taken into account (see also Savaglio et al. 2009). The difference from other Hα-based SFR measurements (Hammer et al. 2006; Christensen et al. 2008) is entirely due to either an overestimated dust correction or a different IMF.

We then obtain further physical parameters from the extracted spectra of the explosion site, the central part of the Wolf-Rayet region, and the integrated-galaxy spectrum as summarized in Table 1. Generally, there is decent agreement between many of the physical properties of the galaxy and the SN site spectrum. Not unexpectedly, resolved measurements of EWs are significantly higher, whereas both dust reddening (by 0.03 mag) and oxygen abundance (by 0.10.2 dex) are marginally lower at the explosion site than for a galaxy-integrated spectrum. However, ESO184-G82 does not provide strong evidence that GRB position spectra are markedly different from galaxy-integrated values except for the obvious mismatch in equivalent width measurements. Whether this observation remains valid for a larger sample of GRB hosts remains to be seen, of course.

Comparing Table 1 with Figs. 7 and 9 also illustrates that the oxygen abundance derived from a galaxy-integrated spectrum is not the central abundance of the galaxy (see also, e.g., Galbany et al. 2016b). As the metallicity determination in an unresolved case is a SFR-weighted measurement, it is dominated by H ii regions that are primarily located somewhat offset from the center, as illustrated in Fig. 4, and thus corresponds to a measurement at a distance of around 2 kpc to 3 kpc from the center of ESO184-G82.

In general, GRB positions tend to be highly concentrated on the UV-brightest regions of their hosts (e.g., Svensson et al. 2010; Blanchard et al. 2016; Lyman et al. 2017). These are the same regions that dominate the nebular-line emission in an unresolved, SFR-weighted spectrum from which physical parameters such as metallicity or dust content are typically inferred at higher redshift.

5. Conclusions and summary

The first solid observational evidence that linked GRBs with core-collapse supernovae was provided by SN 1998bw, and since then this object has become prototypical for SNe following GRBs in terms of luminosity, temporal evolution and spectral properties. In particular the lack of hydrogen and helium and the presence of broad metal absorption lines, characteristic of a photosphere expanding at high velocities, have been observed in most GRB-SNe (Hjorth & Bloom 2012; Cano et al. 2016).

In this article, we use spatially resolved spectroscopy obtained with the novel integral field unit MUSE at the VLT to study the properties of the hot gas phase and the stellar population in the SN 1998bw host galaxy with a particular emphasis on the physical parameters at the SN explosion site. Our data cover the largest part of the host galaxy with individual spaxels of size 35 × 35 pc. The effective spatial resolution is limited by the atmospheric conditions and characterized by the point spread function with a full width half maximum of 160 pc. The main results derived from our analysis can be summarized as follows:

  • (i)

    GRB 980425/SN 1998bwexploded in a young (58 Myr) and dust-poor () environment. The age of the stellar population corresponds to lifetimes of stars with MZAMS between approximately 25 M and 40 M. Our measurement of dust reddening for the explosion site is significantly lower than previous estimates and resolves the mismatch between an apparently dust-rich environment and a non-extinguished SN light curve.

  • (ii)

    The Wolf-Rayet region is located at a projected distance of 860 pc to the GRB explosion site and is extremely young (<3 Myr). This young age would imply implausibly high peculiar velocities for runaway stars and makes scenarios where the GRB progenitor is ejected from the WR region very contrived. The progenitor mass obtained from modeling the nebular spectra of SN 1998bw (Maeda et al. 2006) is consistent with that derived from the explosion environment and is inconsistent with an origin in the WR region. This strongly suggests that the GRB formed in situ.

  • (iii)

    The total, dust-corrected SFR of ESO184-G82 is SFR = 0.22 ± 0.02 M yr-1 inferred from integrating the dust-corrected Hα flux in the respective spaxels. The galaxy has a metallicity gradient of − 0.06 dex kpc-1, suggesting that the typical offsets of at most few kpc between GRBs and their host center have a relatively small impact on the abundance determination for higher redshift GRB sites.

  • (iv)

    Most physical parameters that are derived from emission lines are dominated by the most star-forming regions. An integrated spectrum of the full galaxy and can thus lead to an adequate representation of the actual dust extinction or oxygen abundance at the GRB location. A noteworthy exception are equivalent widths, where an integrated spectrum significantly underestimates the values from the explosion site.

Despite considerable systematic uncertainty stemming from the validity and accuracy of specific strong-line metallicity diagnostics, we reach here the following robust conclusions with respect to abundances:

  • (v)

    Empirical strong-line methods using [O iii]/Hβ and/or [N ii]/Hα fail to produce accurate maps of metallicity at the level of detail probed by our MUSE observation. They significantly under- or overestimate 12 + log (O / H) in regions of high or low ionization parameter, respectively, and therefore return an unphysical radial gradient in individual H ii regions with their centers appearing artificially deprived in metals.

  • (vi)

    A recent method based on photoionization models and [S ii] from Dopita et al. (2016) does not show an obvious dependence of inferred abundance on ionization and returns values that are broadly similar, but somewhat lower (by 0.1 dex to 0.2 dex) than those from electron temperatures via the observed ratios of auroral-to-nebular [S iii] or [O iii].

  • (vii)

    Taking the above considerations into account, we consider 12 + log (O / H) ~ 8.2 at the SN position, 12 + log (O / H) ~ 8.3 for a galaxy-integrated spectrum, and 12 + log (O / H) ~ 8.4 for the nearby WR region as our best estimates of the respective oxygen abundances10.

The immediate environment of GRB 980425 thus indicates a progenitor with Z ~ 0.3 Z and MZAMS between ~25 M and ~40 M. Despite the preeminent efficiency, sensitivity, and optical quality of MUSE, similar studies for a larger sample of GRB hosts remain observationally challenging as a result of the limitations in image quality implied by atmospheric turbulence. However, MUSE will soon be equipped with an adaptive optics module, with which the achievable spatial resolution could decrease to 0.′′3 over the 1′ × 1′ field of view. This brings similar studies (e.g., Izzo et al. 2017) into the realm of possibility for at least a handful of close GRBs and will allow us to reassess our constraints on GRB progenitors on a statistically more significant sample.


2

Constraints from other, age-sensitive lines, such as Hβ or He i (González Delgado et al. 1999), and the stellar continuum (Appendix A) generally agree with this age range.

3

These are strictly lower limits. The fact that the compact WR region is convolved with the seeing-introduced spatial scale of ~0.′′9 leads to a smoothed EW distribution. In fact, archival VLT spectroscopy discussed in the Appendix were taken under significantly better atmospheric conditions and yield EW[ O iii ] (λ5007) ~ 2000 Å (Appendix A).

4

Potential reasons include uncertainties in the flux calibration or the stellar Balmer absorption correction; the disagreement in AV seems strongest where the Hβ fluxes are lowest.

5

Despite the filamentary structure of nearby giant H ii regions such as 30 Doradus, they are usually adequately described with abundances that are homogeneous throughout the region (e.g., Pellegrini et al. 2011, and references therein).

6

As MUSE does not cover the wavelength range of [S iii](λ9532), we use a theoretical value of [ S iii ] (λ9532) = 2.44 × [ S iii ] (λ9069) (Mendoza & Zeippen 1982) here.

7

Sulfur and oxygen are both α-process elements produced in massive stars and observed to track each other well in different environments (see, e.g., Fig. 6 in Izotov et al. 2006).

8

A O3N2-based method returns a decreasing abundance toward the center, similar to what is seen also in a third of SN hosts in this diagnostic (Galbany et al. 2016b).

9

GRBs are commonly found in dwarf galaxies with an irregular morphology (e.g., Fruchter et al. 2006; Lyman et al. 2017). Galaxy centers and thus positional offsets are hence sometimes difficult to determine robustly.

10

These values are based on the Dopita et al. (2016) diagnostic, but scaled upward by 0.15 dex as indicated by the corresponding Te abundances.

11

We would measure significant dust reddening if we did not correct our Hβ, Hγ, or Hδ fluxes for stellar Balmer absorption.

Acknowledgments

We are grateful to the referee as well as S. Schulze, M. Michałowski, G. Leloudas, and D. A. Kann for very helpful and constructive comments that increased the quality and strength of this paper. It is a pleasure to thank D. Malesani for providing broadband imaging of the SN field and L. Christensen, J. Greiner, R. Yates, T. W. Chen, J. Graham, and P. Wiseman for helpful discussions. T.K. and P.S. acknowledge support through the Sofja Kovalevskaja Award to Patricia Schady from the Alexander von Humboldt Foundation of Germany. L.G. was supported in part by the US National Science Foundation under Grant AST-1311862. We acknowledge the use of NumPy and SciPy (Walt et al. 2011) for computing and matplotlib (Hunter 2007) for creating all plots in this manuscript.

References

Appendix A: Analysis of archival long-slit spectra

The ESO archive contains a number of public, high S/N, but low-resolution spectra obtained with the FOcal Reducer/low dispersion Spectrograph 2 (FORS2; Appenzeller et al. 1998), which are particularly interesting in the context of this work. These are the same spectra used by Hammer et al. (2006) and were taken with a 1.′′0 slit on 2004-07-15 with the grism 600B (3450 Å to 6050 Å at a resolving power R of R ~ 850) and on 2004-07-16 with grism 600RI (5300 Å to 8450 Å and R ~ 1050). The position angle during both observations was 28° from east such that both the SN and WR region are covered by the slit (Fig. 1). The total exposure time was 1350 s (each three individual frames with 300 s and 150 s) in the 600B setup and 2250 s in the 600RI (each five single frames with 300 s and 150 s) setup. We reduced and analyzed the archival FORS2 data using standard procedures, specifically the ESO FORS2 pipeline in its version 5.3.8 and self-written methods and algorithms in python (Krühler et al. 2015).

thumbnail Fig. A.1

FORS2 600B spectrum of the brightest pixel of the WR region. Black is the original spectrum and red represents the fitted stellar component. The upper panel shows the full flux range, while the lower panel is a zoom into the stellar continuum.

Open with DEXTER

thumbnail Fig. A.2

FORS2 600B spectrum of the explosion site, i.e., the SN region. Black is the original spectrum and red represents the fitted stellar component. The upper panel shows the full flux range, while the lower panel is a zoom into the stellar continuum.

Open with DEXTER

Table A.1

Emission-line fluxes and EW measurements from archival FORS2 long-slit spectra for the GRB 980425 host galaxy.

In particular the 600B spectra (Figs. A.1 and A.2) were taken under excellent atmospheric conditions, leading to a width of the spectral line spread function of FWHM = 0.′′6 at 5000 Å as evidenced by the trace of a bright stellar source serendipitously on the slit at a distance of 22.′′0 to the WR region (seen at RA(J2000) = 19h35m02s.00, Dec(J2000) = − 52°50′21.′′1 in Fig. 1). The 600RI data have somewhat worse spatial resolution (FWHM = 1.′′0 at 7000 Å). This mismatch leads to the complication that the two different spectra are convolved with very different spatial scales. Because of the small angular size of the WR region core (0.′′1 from HST imaging) and the gradient in physical properties, ratios of lines observed in the two different setups (e.g., for a determination of the dust reddening) are thus clearly nontrivial to interpret.

Table A.1 contains line fluxes and equivalent widths from our analysis of the FORS2 spectra that are sometimes significantly different from the original Hammer et al. (2006) values. Their actual errors on the line fluxes remain unfortunately unclear, but a substantial uncertainty of at least 25% must be present; this value was estimated from the obvious discrepancy between their line-flux ratios of [O iii](λ5007)/[O iii](λ4959) and [N ii](λ6584)/[N ii](λ6548) to the theoretical value of 2.98 set by transition probabilities and observed in high-quality SDSS spectra (e.g., Storey & Zeippen 2000; Osterbrock & Ferland 2006; Wesson et al. 2016).

From the Hβ/Hγ and Hβ/Hδ ratio, we measured EBV = 0.22 ± 0.03 mag for the WR and mag for the SN region, which are both perfectly consistent with the MUSE data11. Also other properties are very similar to our IFU-based values as given in Table 1.

Finally, we exploit the bluer response of the FORS2 600B grism to obtain an age of the H ii region through a fit using composite stellar templates in starlight in a similar manner as in the main text. The stellar population in the WR region is extremely young and the Bruzual & Charlot (2003) templates with an ages between 1 Myr and 3 Myr dominate the best fit by contributing ~60% to 90% to the total observed starlight in various fits using different spectral templates and reddening laws. This age is consistent with the extremely high EW of [O iii] and Hα (Table A.1) and the constraints derived in the main text.

The SN region has prominent stellar components with ages of 5 Myr and 40 Myr each contributing around 30% to the best-fit composite template; the younger of these components is again consistent with the age estimate from the Hα EW in the main text.

Appendix B: Dependence of O3N2-based oxygen abundance on the ionization parameter

thumbnail Fig. B.1

Dependence of the inferred oxygen abundance in two strong-line diagnostic ratios (upper panel: O3N2 from Pettini & Pagel 2004, lower panel: D16) on the ionization parameter U (from [S iii]/[S ii]) for the WR region. Each data point corresponds to a single spaxel and the grey region indicates the constraints from the temperature-sensitive [O iii](λ4363) emission line.

Open with DEXTER

To better illustrate how a changing ionization parameter U affects the metallicity measurement in the different strong-line diagnostics, we plot U (defined as ionizing photons per hydrogen atoms) versus 12 + log (O / H) in the Pettini & Pagel (2004) or Dopita et al. (2016) scale of the brightest H ii region in ESO184-G82 in Fig. B.1. Here, we use the [S ii]/[S iii] ratio (Fig. 6) to calculate U via photoionization models (Dors et al. 2011). Different parameterizations in U (e.g., Morisset et al. 2016) do not change Fig. B.1 significantly. The [S ii]/[S iii] ratio has the advantage that it is nearly insensitive to metallicity, so we should not measure a strong correlation between both quantities in accurate metallicity diagnostics. It is again clear, however, that the O3N2 metallicity scale only reproduces the oxygen abundance at lower ionization parameters and systematically underpredicts it at higher U. In contrast, the Dopita et al. (2016) diagnostic is mostly independent of ionization parameter, but seems offset by an average ~0.15 dex toward lower 12 + log (O / H) (Fig. B.1, lower panel).

All Tables

Table 1

Physical properties inferred from MUSE IFU spectroscopy for the GRB 980425 host galaxy ESO184-G82.

Table A.1

Emission-line fluxes and EW measurements from archival FORS2 long-slit spectra for the GRB 980425 host galaxy.

All Figures

thumbnail Fig. 1

False-color composite from reconstructed VRI-band images from the MUSE data cube. The image spans approximately 55′′ by 55′′, or 10 by 10 kpc at z = 0.0086.

Open with DEXTER
In the text
thumbnail Fig. 2

Separating stellar and gas-phase components at the wavelength around Hβ and [O iii](λλ4959, 5007). Top: original spectrum (black) and fitted stellar component (red). Middle: zoom-in to the continuum. Bottom: resulting spectrum of the gas-phase contribution (blue) together with the error spectrum (grey).

Open with DEXTER
In the text
thumbnail Fig. 3

Spaxel BPT diagram for ESO184-G82. The solid line denotes the ridge line of SDSS galaxies (Brinchmann et al. 2008), whereas the dotted line represents the z = 0 classification line between star formation and AGN ionization (Kewley et al. 2013). The explosion site, a galaxy average (including all spaxels), and an H ii-region average (including only spaxels with EWHα> 10 Å) are indicated.

Open with DEXTER
In the text
thumbnail Fig. 4

Reconstructed images from the MUSE data cube. Each panel shows the host of GRB 980425 in a different nebular line. Top left: Hα. Top right: [O iii](λ5007). Bottom left: [S ii](λ6731). Bottom right: [S iii](λ9070). All panels are approximately 55′′ by 55′′, or 10 kpc by 10 kpc at the GRB redshift. The effective spatial resolution is given by the point spread function (PSF) indicated in the lower left corner of each image with a full width half maximum (FWHM) of approximately 0.′′9 or 160 pc.

Open with DEXTER
In the text
thumbnail Fig. 5

Dust-reddening distribution in ESO184-G82 as measured through the Balmer decrement. We only show spaxels in which Hβ is detected with a significance of at least 3σ. The image spans 34′′ by 38′′, or 6.1 kpc by 6.8 kpc. The circle denoting the position of SN 1998bw has a diameter of 160 pc.

Open with DEXTER
In the text
thumbnail Fig. 6

Face-to-face comparison between the ESO184-G82 maps of O3N2, often used for abundance determinations and [S iii]/[S ii], a tracer of the ionization state of the hot gas. Only spaxels with S/N> 3 are shown. Image dimensions are similar to Fig. 5.

Open with DEXTER
In the text
thumbnail Fig. 7

Map of 12 + log (O / H) as obtained through the [S ii] and [N ii]-based method of Dopita et al. (2016). Only spaxels with S/N> 3 are shown. Image dimensions are 34′′ by 38′′ or 6.1 kpc by 6.8 kpc, similar to Fig. 5. Methods to estimate 12 + log (O / H) based on electron temperatures return values that are 0.1–0.2 dex higher.

Open with DEXTER
In the text
thumbnail Fig. 8

Electron temperatures (Te) from [S iii] and corresponding sulfur abundances in the WR region. Top left: flux map of nebular [S iii](λ9069). Top right: flux map of auroral [S iii](λ6312). Both flux maps are in units of 10-17erg cm-2 s-1. Bottom left: electron temperatures from [S iii]. Bottom right: 12 + log (S / H) as derived from Te. The solar abundance [O/S] is 1.57 (Asplund et al. 2009), so the 12 + log (S / H) scale corresponds to 12 + log (O / H) = 8.2 to 8.4. All panels are approximately 6′′ by 6′′, or 1 kpc by 1 kpc.

Open with DEXTER
In the text
thumbnail Fig. 9

Metallicity gradient for ESO184-G82 in the D16 scale. Each data point is a spaxel-based measurement of oxygen abundance, where individual spaxels were binned with their direct neighbors to enhance the S/N if necessary.

Open with DEXTER
In the text
thumbnail Fig. A.1

FORS2 600B spectrum of the brightest pixel of the WR region. Black is the original spectrum and red represents the fitted stellar component. The upper panel shows the full flux range, while the lower panel is a zoom into the stellar continuum.

Open with DEXTER
In the text
thumbnail Fig. A.2

FORS2 600B spectrum of the explosion site, i.e., the SN region. Black is the original spectrum and red represents the fitted stellar component. The upper panel shows the full flux range, while the lower panel is a zoom into the stellar continuum.

Open with DEXTER
In the text
thumbnail Fig. B.1

Dependence of the inferred oxygen abundance in two strong-line diagnostic ratios (upper panel: O3N2 from Pettini & Pagel 2004, lower panel: D16) on the ionization parameter U (from [S iii]/[S ii]) for the WR region. Each data point corresponds to a single spaxel and the grey region indicates the constraints from the temperature-sensitive [O iii](λ4363) emission line.

Open with DEXTER
In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.