arxiv: v1 [astro-ph.ga] 16 May 2017

Transcription

1 Astronomy & Astrophysics manuscript no _ap c ESO 2018 October 16, 2018 MUSE-inspired view of the quasar Q , its Lyman α blob, and its neighborhood P. L. North 1, R. A. Marino 2, C. Gorgoni 1, M. Hayes 3, D. Sluse 4, D. Chelouche 5, A. Verhamme 6, S. Cantalupo 2, and F. Courbin 1 arxiv: v1 [astro-ph.ga] 16 May Institute of Physics, Laboratory of Astrophysics, École Polytechnique Fédérale de Lausanne (EPFL), Observatoire de Sauverny, 1290 Versoix, Switzerland 2 Institute for Astronomy, ETH Zürich, Wolfgang-Pauli Strasse 27, 8093 Zürich, Switzerland 3 Department of Astronomy, Oskar Klein Centre, Stockholm University, AlbaNova University Centre, SE Stockholm, Sweden 4 STAR Institute, Quartier Agora Allée du six Août 19c, B-4000 Liège, Belgium 5 Department of Physics, University of Haifa, Mount Carmel, Haifa 31905, Israel 6 Observatoire de Genève, Université de Genève, Ch. des Maillettes 51, 1290 Versoix, Switzerland Received / Accepted ABSTRACT The radio-quiet quasar Q at redshift z = 3.08 is known to be close to a small Lyman α blob (LAB) and to be absorbed by a proximate damped Lyα (PDLA) system. Here, we present the Multi Unit Spectroscopic Explorer (MUSE) integral field spectroscopy follow-up of this quasi-stellar object (QSO). Our primary goal is to characterize this LAB in detail by mapping it both spatially and spectrally using the Lyα line, and by looking for high-ionization lines to constrain the emission mechanism. Combining the high sensitivity of the MUSE integral field spectrograph mounted on the Yepun telescope at ESO-VLT with the natural coronograph provided by the PDLA, we map the LAB down to the QSO position, after robust subtraction of QSO light in the spectral domain. In addition to confirming earlier results for the small bright component of the LAB, we unveil a faint filamentary emission protruding to the south over about 80 pkpc (physical kpc); this results in a total size of about 120 pkpc. We derive the velocity field of the LAB (assuming no transfer effects) and map the Lyα line width. Upper limits are set to the flux of the N v λ , C iv λ , He ii λ1640, and C iii] λ lines. We have discovered two probable Lyα emitters at the same redshift as the LAB and at projected distances of 265 kpc and 207 kpc from the QSO; their Lyα luminosities might well be enhanced by the QSO radiation. We also find an emission line galaxy at z = 0.33 near the line of sight to the QSO. This LAB shares the same general characteristics as the 17 others surrounding radio- quiet QSOs presented previously. However, there are indications that it may be centered on the PDLA galaxy rather than on the QSO. Key words. quasars: general - quasars: Q quasars: emission lines - galaxies: high redshift - intergalactic medium - cosmology: observations - Lyman α blobs (LABs) 1. Introduction Lyman α (Lyα) nebulae have been discovered around many quasi-stellar objects (QSOs) and radio-galaxies at high redshift through slit spectroscopy (Heckman et al. 1991; Lehnert & Becker 1998; Bunker et al. 2003; Villar-Martín et al. 2003; Weidiger et al. 2005; van Breugel et al. 2006; Willott et al. 2011; North et al. 2012), narrow-band imaging (Hu & Cowie 1987; Steidel et al. 2000; Smith & Jarvis 2007; Yang et al. 2014), and more recently, through integral field spectroscopy (hereafter IFS) (Christensen et al. 2006; Francis & McDonnell 2006; Herenz et al. 2015; Borisova et al. 2016). They easily span tens of physical kiloparsecs (pkpc) and may be larger than 100 pkpc. Those that host a radio galaxy or a radio-loud QSO are brighter than those hosting a radio-quiet QSO (RQQ), presumably because in the former, interactions between the relativistic jet and Based on observations collected at the European Organisation for Astronomical Research in the Southern Hemisphere under ESO programme 60.A-9331(A) the circum-galactic gas enhances ionization (Villar-Martín et al. 2003). The Lyα nebulae surrounding an RQQ or that host no remarkable object are also called Lyα blobs (LABs). Although they do not necessarily host a QSO or active galactic nucleus (AGN), they are found preferentially in overdense regions (Erb et al. 2011; Matsuda et al. 2004, 2009, 2011; Prescott et al. 2008; Yang et al. 2009). These nebulae are interesting because they trace the circum-galactic medium (CGM) at high redshifts, thereby probing galaxy formation, especially if their Lyα fluorescence is enhanced by the ultraviolet (UV) radiation of the QSO, as outlined by Haiman & Rees (2001) and Cantalupo et al. (2005, 2007, 2012). Although most LABs remain confined within the virial radii of their respective haloes, recent observations have uncovered remarkably large nebulae extending beyond the virial radius (Cantalupo et al. 2014; Martin et al. 2014a,b, 2015; Hennawi et al. 2015), with typical sizes of 400 pkpc. These outstanding observations, some of them made possible by new-generation IFSs like the Palomar Cosmic Web Imager (Martin et al. 2014a,b, Article number, page 1 of 16

2 A&A proofs: manuscript no _ap 2015), provide the first evidence of the cosmic web filaments predicted by large-scale ΛCDM simulations. The frequency of LABs around RQQs has remained an open question (e.g., North et al. 2012) until the new-generation IFS Multi Unit Spectroscopic Explorer (MUSE) attached to the fourth unit of the VLT at ESO was used to observe 17 RQQs at redshifts 3 < z < 4 and found an LAB around each of them (Borisova et al. 2016, hereafter B16). Therefore, in this redshift range and for a sensitivity of about erg s 1 cm 2 arcsec 2, the frequency of LABs surrounding RQQs is 100%. Moreover, the size of such nebulae always reaches or exceeds 100 pkpc, and their surface brightness (SB) profiles are consistent with a power law with a slope of about 1.8. Obviously, the ubiquity and similarity of the LABs that enshroud RQQs raise the question of their power source and emission mechanism. In addition to fluorescent emission from photoexcitation or photoionization (e.g., Geach et al. 2009), other power sources have been suggested, like cooling radiation where gravitational energy of the gas collapsing in the potential well of the dark matter halo is converted into Lyα radiation after collisional excitation and ionization (known also as the cold accretion scenario, see, e.g., Haiman et al. 2000; Fardal et al. 2001; Dijkstra & Loeb 2009; Rosdahl & Blaizot 2012). Shock heating was also proposed (Taniguchi & Shioya 2000; Mori et al. 2004). Finally, scattering of the Lyα photons emitted by a central source (like an AGN) may also occur, betrayed by a high level of linear polarization in the external parts of the nebula (Lee & Ahn 1998; Rybicki & Loeb 1999; Dijkstra & Loeb 2008, 2009; Dijkstra & Kramer 2012). This effect was first observed by Hayes et al. (2011) and may coexist with the cold accretion scenario (Trebitsch et al. 2016). Additional lines such as He iiλ1640 can potentially help distinguish between the cold accretion and the photoionization scenarios, although very deep obervations seem to be needed (Arrigoni Battaia et al. 2015). Neutral hydrogen clouds may lie on the line of sight of a QSO and produce a so-called damped Lyman alpha (DLA) absorption system, as soon as the column number density exceeds cm 2. When the velocity difference between the DLA and the background QSO is v < 3000 km s 1, the two objects are considered as being associated and the DLA is then classified as a proximate damped Lyman alpha (PDLA) system (Ellison et al. 2002, 2010). Even though the association may not be physical in the gravitational sense, it can be manifested through enhanced incidence of highly ionized species, for example, showing that the absorber is influenced by the radiation field of the QSO, as Richards (2001) phrased it. The frequency of PDLAs is close to 2% of the QSOs observed in the redshift range 2.2 < z < 5, and at redshift 3, it exceeds the frequency expected assuming that they are drawn from the same population as the intervening DLAs (Prochaska et al. 2008; Ellison et al. 2010). In addition to its intrinsic interest, a PDLA provides the advantage of acting like a natural coronagraph that extincts the QSO light in the vicinity of the Lyα line, allowing us to detect and study any emission nebula surrounding the QSO more easily (Hennawi et al. 2009). This configuration, where the PDLA lies in front of the QSO and its Lyα nebula, in principle offers the possibility to image the PDLA in silhouette against the extended screen of Lyα emission, as outlined by Hennawi et al. (2009). In this work we focus on an LAB lying along the line of sight of the QSO Q (RA 21 h 02 m 44 s.7, Dec , V= 18.6) at redshift z = 3.09 (Warren et al. 1991) 1. The LAB is emitting at 4970 Å, within the damped Lyα absorption (DLA) trough of the QSO spectrum. The DLA that we analyze in this work is at redshift z abs = 3.082, that is, only 500 km s 1 from the background QSO, if the latter is indeed at redshift z = 3.09; our revised redshift estimate is rather z = (see Sects. 3.1 and 3.3 below), implying a velocity difference of km s 1 (see Sect. 3.2 below). This would represent proper Mpc if the velocity difference were entirely due to the Hubble flow 2. An early study of this QSO and of its Lyα features has been carried out by Leibundgut & Robertson (1999). Based on long-slit spectra, they showed that the emission feature is extended both spatially and spectrally, and they provided evidence for a shift of the central wavelength as a function of the offset position from the QSO. Here, we take advantage of the IFS to provide a detailed 2D mapping of the structure. In our analysis, we assume a ΛCDM universe model with parameters H 0 = 69.6 km s 1 Mpc 1, Ω M = 0.286, and Ω Λ = (Bennett et al. 2014), and we use the cosmology calculator made available by Wright (2006) on the World Wide Web Observations and data reduction The MUSE instrument is an IFS located at the Nasmyth B focus of Yepun, the VLT UT4 telescope. It is operating in the visible wavelength range from 4800 Å to 9400 Å. Thanks to its wide field of view, 1 1 in the wide field mode, it is able to take simultaneous spectra with a spatial sampling of 0.2 per pixel (also referred to as "spaxel"). MUSE has a resolving power R = 1835 (Richard & Bacon 2014) at the redshifted Lyα wavelength, which in terms of velocity separation corresponds to a resolution v of v = c/r 163 km s 1 where c is the speed of light in the air. Or equivalently in terms of wavelength, λ = λ/r 2.7 Å, where λ 4970 Å is the wavelength of the Lyα blob emission at redshift The observations were made in June 2014 during the science verification phase (ESO program 60.A-9331(A)). The total exposure time amounts to 3 hours 20 minutes, distributed in four observing blocks, each of which included two 1497 s exposures taken at instrumental position angles differing by 90. Dithering of a few arcseconds was also applied. The data set was reduced using both the ESO MUSE pipeline (version 1.6, Weilbacher 2015) and the CubExtractor package (CubEx hereafter, version 1.6; Cantalupo, in preparation). In particular, the basic reduction steps of bias and flat-fielding corrections, wavelength solution, illumination correction, and flux calibration were made for each of the individual exposures with the ESO pipeline using the default parameters. After the calibration files were processed, the raw science data were then reduced using the scibasic and scipost recipes in order to obtain the datacubes and pixel tables for each exposure sampled in a common grid of Å. Then, we made use of 1 This is the original reference. The redshift value is mentioned as uncertain for this particular object, while the average error of most other objects is announced to be better than Article number, page 2 of 16

3 P. L. North et al.: MUSE-inspired view of the quasar Q , its Lyman α blob, and its neighborhood 3.1. Redshift of the QSO We compared the QSO spectrum extracted from the MUSE datacube with the template QSO spectrum built by Vanden Berk et al. (2001) on the basis of SDSS spectra. The latter have about the same resolution as MUSE, making them appropriate for the task. The most useful emission line for redshift determination is C ivλ1549, although other features such as Si ivλ1396/o iv]λ1402 and C iii]λ1908/fe iii can also help. After enhancing the line amplitude with respect to the continuum in the template spectrum (by a factor of 1.8 for C iv and 1.2 for C iii]), the best visual match is obtained for z ± 0.003, (1) as shown in Fig. 2 This is slightly larger than z = 3.09 estimated by Warren et al. (1991), but within their 0.01 error bar. Our redshift estimate is largely based on the C iv line, which is blueshifted on average by 810 km s 1 with respect to the systemic velocity for radio-quiet quasars (Richards et al. 2011); however, this should also be the case of the template spectrum, so this is not likely to cause any severe bias. Fig. 1. White-light image obtained by integrating the datacube from 4980 Å to 9350 Å. The positions of the four objects discussed in this paper are indicated. The field of view is one arcminute on a side. the CubEx package to compute a more accurate astrometry, to improve the quality of the flat-field correction (CubeFix tool), and to perform a flux-conserving sky subtraction on the individual flatfield-corrected cubes (CubeSharp). We refer to B16 for a detailed description of the CubEx package. Finally, the eight individual exposures were combined into a final cube using a 3σ clipping and mean statistics. The resulting point-spread function (PSF) measured at 7000 Å on the final combined cube has a fullwidth at half-maximum (FWHM) of At the wavelength of the emitter (4968 Å), it is The white-light image obtained by collapsing the datacube along the wavelength axis from 4980 Å to 9350 Å is shown in Fig. 1. The position of the QSO is indicated, as well as that of the LAEs found at same redshift, and of an emission line galaxy at z = We examine these four objects in turn below. 3. Analysis of the datacube in the QSO vicinity First of all, we have redetermined the redshift of the QSO (Sect. 3.1), and the redhift of the PDLA using metallic absorption lines (Sect. 3.2). Then, we have performed two independent analyses of the immediate neighborhood of the QSO, aimed first at characterizing the LAB. The analysis is complicated by the velocity offset between the PDLA absorption trough and the nebular emission, causing the latter to fall not in the middle of the former, but on its red wing. Thus, the QSO contribution has to be subtracted, even though it remains relatively small, at least in its immediate vicinity, namely within of it. We used the method proposed by Leibundgut & Robertson (1999), which consists of defining the QSO intrinsic spectrum and subtract it in the wavelength domain for each spaxel; for the sake of consistency, we also performed an analysis similar to that presented by B16 (subtraction of the QSO PSF in the spatial plane for each relevant wavelength bin). We describe these two methods in more detail below (Sects. 3.3 and 4.1) Redshift of the PDLA according to metallic absorption lines Ten metallic absorption lines can be easily identified in the PDLA system within the MUSE spectral range. These are Si ii 1193, Si iii 1206, N v 1238, Si ii 1260, O i 1302, Si ii 1304, C ii 1334, Si ii 1526, Fe ii 1608, and Al ii The N v 1238 line is the faintest of all and has the largest residual from the regression of observed versus laboratory wavelengths, so we discarded it. We also discarded the C iv lines because they fall in the middle of telluric lines. The difference between the observed wavelength corrected for the air refraction index (nair ) and the redshifted vacuum wavelength is shown in Fig. 3. A least-squares fit weighted by the observed equivalent widths (EW) of the nine lines gives a redshift zabs = ± (2) with an rms scatter of Å, much smaller than the spectral resolution. Our estimate agrees with the value of given by Leibundgut & Robertson (1999) (based on the Fe ii 1144, Si ii 1193, and Si iii 1206 lines) to within 1.5 σ. There is no sign of any multiple absorption components. The redshift difference of the PDLA relative to the QSO can be translated into a velocity difference using the formulae given by Davis & Scrimgeour (2014): An object with a peculiar redshift z p with respect to a QSO at redshift zem will appear to the observer to have z such that 1 + z = (1 + zem )(1 + z p ) = z p = z zem zem (3) The peculiar velocity v p is linked to z p in the following way for relativistic and non-relativistic velocities, respectively: vp = c (1 + z p )2 1 ' c z p = 930 ± 220 km s 1, (1 + z p )2 + 1 (4) where the error reflects the uncertainty on the QSO redshift. Thus we have v p < 3000 km s 1, justifying the classification of the absorption system as a PDLA. Article number, page 3 of 16

4 A&A proofs: manuscript no _ap Fig. 2. Top left: Observed QSO spectrum (black curve) compared to the normalized SDSS template of Vanden Berk et al. (2001) redshifted to z = (red curve), and to the same template after enhancement of the carbon emission lines (cyan curve, see text). Bottom left: Ratio of the observed QSO spectrum to the redshifted and normalized template with enhanced C lines. Right panels: Enlarged portions of the QSO spectrum showing the metallic lines used to determine the redshift of the PDLA. Fig. 3. Difference between the observed wavelength corrected for the air refraction index and the redshifted laboratory (vacuum) wavelength versus laboratory wavelength for the nine lines mentioned in the text. The red horizontal solid line corresponds to the EW-weighted regression, while the red dashed lines show the ±1σ rms scatter Fit of the Lyman alpha absorption trough and subtraction of the QSO contribution From the whole datacube given by the pipeline, we wish to extract the spectrum of the LAB alone in the vicinity of the QSO. The goal is then to subtract the absorption profile from the observed spectrum to obtain a clear signal of the emission feature without the contribution of the QSO. We first used the ds9 tool to determine the position of the LAB and roughly estimate its extent. This allowed us to select a region encompassing most of the LAB luminosity, and of course the whole QSO contribution. The region we selected is elliptical, with semi-axes of 2.3 and 2.8, the position angle (orientation with respect to the north direction, counted positively to the east) of the semi-major axis being 30. The procedure to isolate the spectrum of the LAB is then the following: We use the software VPFit (Carswell & Webb 2014) to fit a Voigt profile to the absorption trough. The needed error vector was extracted from the variance datacube, and we dropped all irrelevant data, namely (i) absorption lines of the Lyα forest, on the blue side of the DLA trough, (ii) some faint absorption features to the red side of the DLA trough (5000 Å < λ < 5100 Å), (iii) some narrow absorption lines on the blue wing of the DLA trough ( 4930 Å < λ < 4940 Å), and (iv) the Lyα blob emission. These regions are shown in Fig. 4. In order to obtain a meaningful fit, we need to reconstruct the original energy distribution of the QSO (we call it the continuum somewhat abusively, because here lies the broad Lyα line). This is done by iteration over the following steps: We start by performing a first fit, feeding VPFit with our measured spectrum (blue curve in Fig. 4) and a preliminary continuum (close to the green curve in Fig. 4). VPFit will then provide a Voigt profile fit of the DLA. We divide the original spectrum by the DLA fit, obtaining the function f r shown in Fig. 5 (blue line). The upper envelope of f r should tend toward 1 outside the emission region of the LAB (and outside the bottom of the absorption trough) because Article number, page 4 of 16

5 P. L. North et al.: MUSE-inspired view of the quasar Q , its Lyman α blob, and its neighborhood Fig. 4. QSO spectrum with the Lyα DLA feature. The blue line is the measured spectrum, with the emission feature at the bottom of the damped Lyα absorption trough. The Lyα forest lies on the blue side of the DLA Lyα absorption. The red curve shows the Voigt fit of the DLA. The green curve is the reconstructed "continuum", actually the broad Lyα emission line of the QSO (see text). Note that the reconstructed QSO line is only an interpolated guess between 4950 Å and 4958 Å. The thin dashed magenta line is the SDSS QSO template of Vanden Berk et al. (2001) redshifted to z = and normalized to match the green curve as well as possible. one has to ensure that the wings of the Lyα absorption are well fitted, regardless of the other absorption lines. The closer f r is to 1, the better is the continuum. We then manually define the upper envelope of the f r function, neglecting both the bottom of the PDLA absorption (where f r is not defined) and the LAB emission (see Fig. 5, red line). This provides a vector by which we multiply the previous estimate of the continuum to obtain an improved estimae. We note that by construction, this vector remains uncertain in the wavelength range 4948 Å 4977 Å, where we are forced to interpolate the continuum. If the latter were to present any discontinuity or peak within that range, neither this method nor any other would be able to unveil it. Four such iterations allowed us to converge to a satisfactory fit of the QSO+PDLA spectrum, shown in Fig. 4. We note that the shape of the resulting broad emission Lyα line is remarkably close to that of the template spectrum, after the latter has been properly normalized (i.e., divided by the linear function y = λ [Å] ). We can now scale the fitted spectrum and subtract it from the observed spectrum in each spaxel concerned by the QSO contribution. This is done in the following way. Let I re f be the spectrum extracted in the first part (the blue curve in Fig. 4), and s i j the spectrum in a particular spaxel located at position (i, j) in the image. The scaling factor h i j for a specific spectrum is given by the ratio between the sum of s i j, selected in a small wavelength range D λ = [ Å, Å] where the spectrum is nearly Fig. 5. Illustration of the second and third step of the procedure for estimating the QSO continuum. The highest values of f r (blue curve) must tend to 1 (outside the bottom of the PDLA absorption and the LAB emission), when the profile for the DLA converges to a good solution; the blue curve diverges in the wavelength range corresponding to the bottom of the PDLA absorption and to the LAB emission, where it has to be interpolated. The red curve shows the spline interpolation performed over the highest points of the blue curve. The spline function has then to be multiplied by the previous estimate of the continuum to create a better estimate. flat, and the sum of I re f in the same range: h i j = s i j (λ) λ λ D λ, (5) I re f (λ) λ λ D λ where λ = 1.25 Å is the spectral sampling of our image. Thanks to the scaling factor, we can estimate a fit for the PDLA profile in each spaxel. If F re f is the QSO+PDLA spectrum computed above (red curve in Fig. 4), the QSO+PDLA spectrum f i j in each spaxel is given by f i j = h i j F re f. (6) Therefore, the LAB spectrum s i j in each spaxel, without the QSO contribution, is given by subtracting the estimated fit f i j to the measured spectrum s i j : s i j = s i j f i j. (7) We now have a new subcube s i j from which the signal of the QSO has been subtracted. 4. Exploration of the datacube at larger scales We have explored the datacube at scales larger than just a few arcseconds from the QSO, both visually and in an automated way. Visual exploration was made by scanning the cube in increasing wavelength bins and by displaying the spectrum of the few emission line objects that were detected. It also allowed us to discover a faint extension of the LAB protruding to the southsouthwest for some 10. Automated exploration aimed at defining the extension of the LAB in a more objective and sensitive way was performed with Article number, page 5 of 16

6 A&A proofs: manuscript no _ap CubEx, which in addition to its post-processing capabilities, is also able to exploit simultaneously the 3D information of the datacube to automatically extract extended emission above some predefined signal to noise ratio (S/N) threshold. More details are given hereafter and by B16. This method allowed us to deeply explore the filamentary structure we had assumed to be there based on visual exploration, which is oriented to the south over more than 10 ( 100 pkpc), and to better define its shape and surface brightness PSF quasar subtraction from each monochromatic image We also explored an alternative empirical method that is implemented in CubEx to minimize the contribution of the QSO PSF and obtain the pure extended Lyman alpha emission. The principle of this method has been widely used in imaging astronomy (Husemann et al. 2013) and is very well suited to IFS data (Herenz et al. 2015). In summary, an empirical PSF is scaled and subtracted from the QSO at each wavelength layer. We used the CubePSFSub tool, which is also part of the CubEx package, to create the pseudo-nb images (150 spectral pixels or 187 Å) at the position of the QSO, then the flux measured in the 1 1 central region was rescaled in each empirical PSF image and finally subtracted from the corresponding wavelength layer in the datacube (see Sect. 3.1 in B16 for further details). Although this method produces excellent results far away from the quasar, it does not provide robust results in the central 1 region. Ultimately, we used the latter CubEx approach to detect the extended LAB and the approach described in Sect. 3.3 to study the PDLA properties. 5. Results 5.1. Properties of the PDLA and of the bright component of the Lyman alpha blob We first focus on the properties of the PDLA obtained through the VPFit code. Together with the fitted absorption profile, both the column density and the redshift of the PDLA are derived: log(n) = ± ; z abs = ± (8) The redshift value obtained here differs by no more than from the value obtained based on metallic lines and is thus entirely consistent with it (see Sect. 3.2 above). We took the refractive index n of air into account that causes a shift in the observed wavelengths: λ air = λ vac /n. We used the Ciddor equation (Ciddor 1996), taking the atmospheric conditions at the time of the measurement into account (air temperature 12 C, pressure hpa, and relative humidity 3%, as reported in the header of the original FITS file). The errors are the formal errors provided by VPFit and are certainly underestimated, as they do not include uncertainties linked to the continuum determination. In their pioneer study, Leibundgut & Robertson (1999) obtained a column density and redshift log(n) = ± 0.03 ; z abs = Our estimate of the column density is higher than theirs by 2.6σ (referring to their error estimate), and we find a slightly lower redshift (the difference amounts to 7.6σ according to our own error estimate). The latter discrepancy might be due to our lower spectral resolution (R 1840 instead of 3300), which might cause the metallic absorption lines in the blue wing of the Lyα Flux F [erg s 1 cm 2 Å 1 ] Bimodal Gaussian fit 1 = ± 0.2 Å 2 = ± 0.3 Å Emission wavelength em [Å] Fig. 6. Spectrum of the emission feature averaged over an elliptic surface of , centered on the QSO, but free of its flux. The heavy continuous curve is a fit performed with the sum of two Gaussians (dash-dotted lines). The wavelengths indicated are corrected to vacuum. Table 1. Central position λ (corrected to vacuum) and FWHM of the two Gaussian components in the LAB spectrum. The third line gives the same quantities for the total emission profile. λ FWHM FWHM [Å] [km s 1 ] [Å] first peak ± ± ± 0.7 second peak ± ± ± 0.4 first moment ± ± ± 0.20 absorption trough to be less easily recognized and be partly included in the Voigt profile fit. This would result in a slightly overestimated column densitiy and in a slight blue shift. However, our VPFit estimate of the redshift, based on the Lyα line alone, agrees within only 0.5σ with the PDLA redshift estimate based on metallic lines. We can determine the average spectral shape of the LAB emission feature over the elliptic surface defined previously by subtracting the DLA fit to the measured spectrum. The result, shown in Fig. 6, is a clearly asymmetric line that can formally be fitted by a pair of Gaussian functions (Fig. 6). The parameters of the two Gaussians are listed in Table 1. We note that these two functions are only used as a convenient way of representing the average emission line and do not bear straightforward physical meaning. The total Lyα flux F Lyα of the emission feature is computed by integrating the flux over the double Gaussian profile: F Lyα = (2.53 ± 0.05) erg s 1 cm 2. (9) The mean redshift of the LAB is defined as z = λ em λ 0 λ 0 = ± , (10) where λ 0 = Å is the Lyα wavelength at rest (in the vacuum) and λ em is the observed wavelength of the emission corrected to vacuum. The (vacuum) wavelength λ em = ± 0.1 Å is determined by the first-order moment of the distribution. Article number, page 6 of 16

7 P. L. North et al.: MUSE-inspired view of the quasar Q , its Lyman α blob, and its neighborhood This is a reasonable choice because it takes the skewness in the shape of the spectrum into account. The redshift of this bright central part of the LAB is therefore very close to but different from that of the QSO: the difference amounts to z p = or 462 km s 1, about twice the uncertainty on the QSO redshift (assuming z QSO = ± 0.003). The velocity difference between the PDLA and the LAB is then v = v(pdla) v(lab) = 473 ± 24 km s 1, in perfect agreement with the 490 km s 1 given by Leibundgut & Robertson (1999). From the redshift, we determine the luminosity distance d L and the total Lyα luminosity of the blob: δ L Lyα = 4πd 2 L F Lyα = (2.18 ± 0.05) erg s 1 (11) 2 α with d L = Mpc. The error given here is a formal error provided by the fit; a more realistic error would be closer to 10%, taking into account the uncertainty of the subtracted QSO spectrum. The FWHM of the average emission spectrum (not corrected for the instrumental width) is: FWHM = 9.27 ± 0.20 Å or FWHM = 559 ± 12 km s 1 (12) if we naively admit that the line width reflects the velocity field alone, without any transfer effect. This rather modest value already indicates that the gas is not subject to violent motions, especially as part of the line width may be accounted to transfer effects, making our FWHM estimate an upper limit rather than a true estimate of the kinematics. Leibundgut & Robertson (1999) obtain FWHM values varying from 300 to 540 km s 1 (or even 650 km s 1, but for a line affected by a cosmic ray), which is clearly smaller. However, their values refer to individual observations of various parts of the LAB, so they do not include the large scale velocity field (or transfer effects). They also use a higher spectral resolution, which tends to decrease the measured FWHM, although this accounts for no more than a difference of 25 km s Spatial shape of the central LAB component Summing the successive layers in the wavelength range of the LAB, we obtain an image of the blob and can determine its extension and shape. In each spaxel at position (i, j) of the resulting stacked image we have the sum of all pixel values I i j (λ), located at the same position (i, j), over the interval D λ ranging from 4950 Å to 4990 Å encompassing the LAB emission: Ī i j = I i j (λ) λ. (13) λ D λ The result is shown in Fig. 7. The coordinates of the blob center are RA = 21 h 02 m s, Dec = (14) We see that this bright LAB component has a roughly circular shape with an extension of about 3 4. This corresponds to a diameter of kpc at the LAB redshift. Fig. 7. Image of the LAB obtained by addition of the successive frames from 4950 Å to 4990 Å. The signal of the QSO has been removed, as explained in the text, and the position of the QSO is shown by the green dot. The red curves are isophotes spanning SBs from erg s 1 cm 2 spaxel 1 to erg s 1 cm 2 spaxel 1, or erg s 1 cm 2 arcsec 1 to erg s 1 cm 2 arcsec 1. The outermost isophote corresponds to S/N 1.8, while the inner corresponds to S/N 12 per spaxel. There is a slight offset of the QSO position by about 0.4 with respect to the isophotal center of the blob. The relative coordinates are defined with respect to the field center Velocity map of the LAB central part To create the velocity map of the blob 4, we fit a Gaussian profile on each spectrum s i j of the subcube. From the mean position value of the Gaussian fit, we obtain the redshift z. We wish to obtain the peculiar redshift z p at each location of the blob, however, which has a global redshift z em = By applying relations (3) and (4) to each spectrum of the subcube, we obtain the velocity map of the blob, shown in the left panel of Fig. 8. We can clearly distinguish the presence of two regions, each with a relatively uniform velocity, and separated by a steep north-south gradient from about 80 km s 1 at the LAB center, 0.4 south of the QSO, to about 100 km s 1 2 south of the QSO. Leibundgut & Robertson (1999) have shown some evidence for this feature (see their Fig. 4c), but their data could not provide a complete map such as ours, which also shows an east-west velocity gradient 1 to the east of the LAB center. We also note that our peculiar velocities overall range from about 120 km s 1 to +80 km s 1. This 200 km s 1 difference is not much larger than our spectral resolution, which is 163 km s 1, but it is quite significant. It agrees well with the 240 km s 1 difference found by Leibundgut & Robertson (1999, Table 2 and Fig. 4). The right panel of Fig. 8 is a map of the FWHM (obtained through the Gaussian fit) of the Lyα emission line. When we consider both panels, a rough overall trend of an increasing FWHM with increasing velocity seems to appear: the FWHM is around 300 km s 1 2 south of the QSO, but closer to 600 km s 1 at the LAB center. This is also consistent with Leibundgut & Robertson (1999, Fig. 4d). 4 The velocity word used here is a proxy for line position and is strictly justified only when transfer effects can be neglected. Article number, page 7 of 16

8 A&A proofs: manuscript no _ap Fig. 8. Left: Velocity map of the Lyα emission feature. The color code highlights the velocity along the line of sight in each spaxel of the image, the zero-point being set at z = or λ vac = Å. The curves are the isophotes labeled in units of erg s 1 cm 2 spaxel 1. Below erg s 1 cm 2 spaxel 1, the S/N is very low and everything outside this isophote is strongly affected by background noise. The green point near the center shows the position of the QSO. The dashed rectangle represents the region of the image used in Fig. 9; we have enhanced it by lowering the contrast in the area that surrounds it. Right: Map of the FWHMs of the Lyα emission feature. The FWHM in each spaxel is computed from the Gaussian fit of the emission spectrum. For comparison, the average FWHM measured from the total emission spectrum is 560 km s Fig. 9. Map of spectra recorded in individual spaxels. The total area corresponds to the dashed rectangle in Fig. 8; it covers spaxels ( ). North is up and east to the left. The blue lines show the measured spectra, and the red curves the Gaussian fits. The spectral range is 4955 Å 4980 Å. The black number inside each subplot (top left corner) is the peculiar velocity in km s 1 computed from the Gaussian fit; the red number just below is the FWHM of the emission line in km s 1 ; the purple number (top right) is the amplitude of the Gaussian in units of erg s 1 cm 2 Å 1 spaxel 1. The dashed vertical line marks the average emission wavelength ( Å). The quasar lies in the upper left quadrant, in the spaxel [7,5] just next to spaxel [7,4] (where [1,1] and [7,8] define the lower left and upper right spaxels). Article number, page 8 of 16

9 P. L. North et al.: MUSE-inspired view of the quasar Q , its Lyman α blob, and its neighborhood δ " 20" 10" E 100 kpc N Q " 0" 10" 20" 30" α 2000 Fig. 10. Optimally extracted map of the LAB, obtained using the CubEx software on the PSF and continuum subtracted MUSE datacube, and by collapsing along the wavelength axis the datacube voxels inside the 3D segmentation map. The spatial projection of the 3D mask is indicated by the thick contour that corresponds to an SB of about erg s 1 cm 2 arcsec 2 and an S/N of 2.3. The thin contours indicate the propagated S/N in the images, the outermost one corresponding to S/N= 2 and the other ones to S/Ns increasing by steps of 6 units (see B16 for a detailed explanation). The QSO position is shown by the black dot. Note the long tail extending southward of the QSO; this was unknown to date. To verify that the different spectra in each spaxel are well fitted by the Gaussian functions, we provide Fig. 9. Each subplot of this figure shows the total flux recorded in each spaxel. The total area of the figure corresponds to the dashed rectangle in Fig. 8. The Gaussian fit appears satisfactory in each subplot where the noise is low enough, which occurs when the fitted amplitude is larger than erg s 1 cm 2 Å Properties of the Lyman alpha blob at larger distance from the QSO and other objects SBLy α (10 18 ergs 1 cm 2 arcsec 2 ) Using the CubExtractor software, we were able to map the whole LAB, especially at large distance from the QSO, as seen in Fig. 10. The LAB extends over 100 pkpc at least, like the 17 other LABs observed around radio-quiet QSOs by B16. The large extension is mainly due to a filament protruding southward and the Lya flux measured within the 3D mask (thick contour) is erg s 1 cm 2. This is lower than what we found for the bright part of the LAB, probably because of the poor subtraction of the QSO light at the QSO position, which artificially reduces the LAB intensity there. The filament is relatively thin but appears resolved, being wider than 1 or 8 pkpc. The morphology of this LAB is quite similar to that of several LABs discovered by Matsuda et al. (2011) through the Subaru Lyα blob survey in the SSA22 area at the same redshift. In particular, SSA22-Sb6-LAB4 has almost exactly the same shape and extent, with a thin filament protruding to the southwest. The only difference is the absence of any obvious QSO. The corresponding velocity map is shown in the left panel of Fig. 11, the zero velocity corresponding to λ vac = Å (maximum of the emission intensity), that is, to 196 km s 1 in Fig. 8 (left panel). It is based on the first moment of the emission line, while the velocity map of the immediate vicinity of the QSO (Fig. 8) is obtained through Gaussian fits. This might partly explain the larger velocity contrast ( 300 km s 1 instead of 200 km s 1 ) near the QSO. Another cause for the discrepancy may be the approximate subtraction of the QSO PSF. Nevertheless, the two velocity maps agree regarding the north south velocity gradient, which appears very steep about 1 south of the QSO. At larger scale, the southern component of the LAB overall has a velocity close to 0 km s 1, while in the close QSO neighborhood, the velocity is higher than 100 km s 1. The map of velocity dispersion (as a proxy for the line width) is shown in the right panel of Fig. 11. We note the very small dispersion in the southern component of the LAB, which is essentially equal to the MUSE resolution limit. The same is true of the southern part of the main LAB component. Near the QSO, the Lyα line is broader on an oblique band stretching from the SW of the QSO to the SE of it. Our previous map (Fig. 8, right panel) does not show this feature so clearly, but does display widths between 240 and 600 km s 1, which are similar. There is rough agreement between the two maps in the sense that the Lyα line tends to be broader near the QSO than far away from it Surface brightness profile and size luminosity relation We show the mean radial variation of the azimuthally averaged surface brightness (SB) in Fig. 12, as measured on the SB map of Fig. 10. The average is taken over circular annuli. Although the LAB is strongly asymmetric, it is worth comparing this radial variation with that of other blobs. As mentioned above, the innermost part (radius 1 ) is not reliable, but the rest is, and the SB curve shown here was obtained in exactly the same way as the 17 others published by B16. It is essentially free of any pollution from the QSO light because it shows only the result beyond 10 kpc of it, that is, beyond 1.3. The red line shows a power-law fit (see B16, Appendix B) with an index α = 1.51, slightly lower than the average shown by the purple line and corresponding to α = 1.8. Only four of the 19 objects (including two nebulae associated with radio loud QSOs) of B16 have a shallower profile. The exceptionally elongated and filamentary shape of our LAB might be responsible for this unusual profile. Figure 12 also shows the SB profile of the Slug Nebula (Cantalupo et al. 2014). We see that in the range of radii considered, the profile of our LAB lies almost an order of magnitude below the average of the B16 sample. Overall, our LAB compares well with the sample of B16, even though it lies on the faint side. This is also confirmed by the relation between luminosity and size shown in Fig. 14. The size is defined here as the major axis of the nebula. The objects observed by B16 are represented as blue stars, while the full magenta dots are the LABs discovered by Matsuda et al. (2011). Interestingly, the latter largely overlap the former, but extend to lower luminosities and smaller sizes. This is probably an observational bias: B16 have selected the brightest QSOs, so one may expect that the luminosity of the associated LABs is enhanced by fluorescence due to the QSO UV radiation. On the other hand, Matsuda s LABs do not harbor any obvious QSO. Thus, our object indeed lies on the faint end of the LABs surrounding bright QSOs, but appears to be a rather average one when compared to the general LAB population, as far as the sample of Matsuda et al. (2011) is representative. Article number, page 9 of 16

10 A&A proofs: manuscript no _ap δ " 10" E N Q kpc Velocity (kms 1 ) δ " 10" E N Q kpc FWHM (kms 1 ) 0" 0" 10" 20" α " 0" 10" 20" α Fig. 11. Left: Large-scale velocity map of the LAB given by the CubEx software on the basis of the first moment of the Lyα emission line. The QSO position is marked by the black dot. The zero velocity corresponds to λ vac = Å, or to 196 km s 1 in Fig. 8 (left panel). Right: Large-scale map of the velocity dispersion in the LAB given by the CubEx software on the basis of the second moment of the Lyα emission line. SB (ergs 1 cm 2 arcsec 2 ) Borisova+16 (QSO) Cantalupo+14 (Slug Nebula) Q (This work) α = C p = distance from QSO (kpc) Fig. 12. Radial variation of the SB of the LAB (black dots; distance in physical kpc). The red line is a least-squares fit of a power law to the data. The average profile of the 17 nebulae associated with radio-quiet QSOs (B16) is shown as a purple line. The yellow line represents the Slug Nebula (Cantalupo et al. 2014), the dashed part of it is polluted by the QSO light. Correcting the observed SB for redshift dimming would bring the yellow curve very close to the purple curve, but would slightly increase the distance between the red and purple curves. The upper limit of the gray area corresponds to the 2σ noise on the profile Upper limits to high-ionization lines We examined the QSO-subtracted datacube for any presence of the N vλ , C ivλ , He iiλ1640.4, and C iii]λ We defined 2σ upper limits to the flux of these lines, using the same method as in B16, and compared them with Table 2. Upper limits to the SB of lines due to highly ionized species. Line, λ (vac.) Flux (2σ limit) [Å] [10 18 erg/s/cm 2 /arcsec 2 ] Flux(line) F(Lyα) N v < C iv < He ii < C iii] < the flux of the Lyα line defined in the same datacube. The results are given in Table 2, and Fig. 13 shows the ratio of the He ii and Lyα fluxes versus the ratio of the C iv and Lyα fluxes for our object and for the B16 sample. Here again, our object appears quite similar to the other nebulae found around radio quiet QSOs. Arrigoni Battaia et al. (2015) have shown that in order to rule out the photoionization scenario, one would need to reach upper limits as low as 0.05 and 0.07 on the He ii/lyα and C iv/lyα ratios, respectively. This is indicated by the green area in Fig. 13: photoionization as modeled by Arrigoni Battaia et al. (2015) is possible within that area, but would imply unrealistic hydrogen column densities (> cm 2 ) beyond it (the authors assume spherical clouds of cool gas). None of the high-ionization lines available in the MUSE spectral range can be seen, and our upper limit on the C iv/lyα flux ratio seems compelling enough to exclude photoionization by the QSO as the main driver of the LAB Lyα luminosity. The upper limit on the He ii/lyα ratio is just below the border mentioned, so that based on this criterion alone, our object would be on the verge of being compatible with photoionization by the central source. Interestingly, only one among the B16 radio quiet targets lies in the photoionization region in Fig. 13, but since its position corresponds to upper limits only, it remains impossible to tell how far its luminosity can really be explained by photoionization. All other objects have a too low line ratio along either one or both axes. Article number, page 10 of 16

11 P. L. North et al.: MUSE-inspired view of the quasar Q , its Lyman α blob, and its neighborhood Fig. 13. Ratios of high-ionization line fluxes for our object (red) and for those of B16 (Table 2) surrounding radio-quiet QSOs (blue). The He ii line is detected in only one object of the latter sample, while the C iv line is detected in three objects. The green hatched area represents the range of line ratios compatible with photoionization models (Arrigoni Battaia et al. 2015). Fig. 14. Relation size luminosity for our object (red) and for those of B16 that are radio quiet (blue). The magenta full dots represent the 14 giant LABs discovered by Matsuda et al. (2011) Two other probable Lyman alpha emitters at the same redshift We have detected two other possible Lyα emission features in the neighborhood of the QSO; their coordinates are given in Table 3. The brighter feature (LAE1) lies 34 away from the QSO, which corresponds to a projected distance of about 265 pkpc, and it emits within the range 4960 Å to 4975 Å. The vicinity to the QSO and the emission at wavelengths identical to the central LAB suggest that this object is emitting in the Lyα line. Figure 15 (left panel) shows the extension of this emitter as found with CubEx, which is about 4.3 or 34 pkpc in diameter; if the eastern extension is real and physically belongs to the LAE, then the major axis of the LAE would reach 6.4 or 50 pkpc. This appears very large, but is consistent with the finding of Wisotzki et al. (2016), who found many examples of LAEs with Lyα emission extending to 2 radius at similar redshift. The emission feature of this nebula consists of two peaks, as shown in the right panel of Fig. 15. This is very typical of some Lyα emitters (LAEs). LAE , for instance, shows a very similar profile (Yamada et al. 2012, Fig. 2). The latter authors concluded from their observation of 91 LAEs that more than 40% of them show a double peak, and that in most cases the red component is stronger than the blue one, which is naturally explained by transfer effects in expanding media, even inhomogeneous ones (Verhamme et al. 2006, 2008, 2012; Gronke & Dijkstra 2016). Its total Lyα flux is F Lyα = (7.6 ± 0.4) erg s 1 cm 2 and the positions and widths of the two peaks are summarized in Table 3. The emission wavelength of the red peak coincides closely (i.e., within 0.6 Å or 36 km s 1 ) with the first moment of the emission line of the LAB surrounding the QSO. Thus, assuming that the double-peak appearance is due to transfer in an expanding medium, the unabsorbed line would be centered at about 230 km s 1 with respect to the LAB emission. Such a small velocity difference suggests that this LAE is indeed close to the QSO+LAB system. At this point, the question arises whether LAE1 might owe at least part of its luminosity to photoionization by the QSO (fluorescence scenario). In this case, even in the most extreme case of a purely external photoionizing source, a double peak similar to the observed peak would be produced (Cantalupo et al. 2005; Verhamme et al. 2015), so the answer is positive. Yamada et al. (2012) observed 91 LAEs with a spectral resolution R 1700, very similar to that of MUSE. They show in their Fig. 5 the empirical relation between the velocity separation between the red and blue peaks and the FWHM of the stronger peak. Our LAE1 object has a rather small velocity separation (392±36 km s 1 ) for the width of the stronger (red) peak, sending its representative point onto the lower envelope defined by the objects of Yamada et al. (2012). It is also reminiscent of the Lyman continuum emitters (LCEs) examined by Verhamme et al. (2017). Using Eq. 11 and the luminosity distance corresponding to the average redshift z = , we obtain the luminosity L Lyα = (6.5 ± 0.3) erg s 1 We do not show any velocity or FWHM map for this emitter because of its small spatial extension. In a white-light image (obtained by integrating the datacube over all wavelengths), we see a very faint object at the position of the Lyα emission, which confirms that we are not in the extreme case of the dark galaxies seen by Cantalupo et al. (2012). The rest-frame equivalent width (EW 0 ) obtained using the same fivepixels radius aperture for both the Lyα line and the continuum (broad band 2000 Å image reconstructed from the MUSE Article number, page 11 of 16

12 A&A proofs: manuscript no _ap δ " 15" 12" 9" 6" 3" E N SBLy α (10 18 ergs 1 cm 2 arcsec 2 ) Flux Fλ [erg s 1 cm 2 Å 1 ] Bimodal Gaussian fit λ 1 = ± 0.6 Å λ 2 = ± 0.2 Å 0" 0" 3" 6" 9" 12" 15" 18" α Emission wavelength λ em [Å] Fig. 15. Left: CubEx optimally extracted map of the probable Lyα emitter LAE1 located 34 away from the QSO. The key to the curves is the same as for Fig. 10. Right: Spectrum of the LAE1 object, summed on a disk with 2.0 radius. The broad continuous curve shows a fit performed with a bimodal Gaussian functions. The continuum has been subtracted and the wavelengths are corrected to vacuum. Table 3. Equatorial coordinates of the two LAEs detected at same redshift as the LAB. The central positions λ (corrected to vacuum) and FWHM of the two Gaussian fits to each LAE Lyα emission line are given as well. The v quantity is the velocity difference between the peak position and the first moment of the average Lyα profile of the LAB surrounding the QSO. Object α (J2000) δ peak λ z FWHM v [h:m:s] [ : :"] [Å] [km s 1 ] [Å] [km s 1 ] LAE1 21 : 02 : : 52 : 42.7 blue ± ± ± ± ± 36 red ± ± ± ± ± 12 LAE2 21 : 02 : : 52 : 42.6 blue ± ± ± ± ± 48 red ± ± ± ± ± 18 datacube and assuming a β slope of 2) is EW 0 < 130 Å, while a dark galaxy is defined by EW 0 > 240 Å. Therefore, LAE1 must be powered from inside, for instance, by some stellar formation, although we cannot exclude some contribution of fluorescence due to the QSO, as mentioned above. The second object, LAE2, appears very similar to LAE1 although it is much fainter and less extended (Fig. 16, left panel). The blue component of the emission line seems brighter compared to the red component than it is in LAE1 (see Fig. 16, right panel), although it is also more affected by the noise 5. LAE2 lies 26.5 away from the QSO, representing a projected distance of 207 pkpc and also shares the same redshift as the LAB one. Its coordinates and characteristics are given in Table 3. Its total Lyα flux within a 1.0 radius is F Lyα = (12.0 ± 1.2) erg s 1 cm 2 For the luminosity distance corresponding to the average redshift z = , we obtain the luminosity L Lyα = (1.0 ± 0.1) erg s 1. 5 The spectrum is drawn here from a datacube from which continuum sources were subtracted, thereby avoiding contamination by galaxies close to the line of sight. The rest-frame equivalent width of LAE2 has been estimated in a r = 6 pixels aperture at EW 0 17 ± 4 Å. The blue peaks of LAE1 and LAE2 occur at exactly the same redshift, but this should be considered a mere coincidence, especially because for LAE1, the fit of the blue peak is very sensitive to noise because the two peaks are not clearly separated. The presence of these two LAEs in the field probably reflects the overdensity that is generally related with a QSO Starburst galaxy near the line of sight of the QSO An emission line galaxy is detected about 5.5 to the NW of the QSO. The H α, H β, [O ii]λ3727, [O iii]λ4959, 5007 lines are clearly identified, while the [N ii]λ6584 line marginally so. The five stronger lines allow us to determine the redshift of this galaxy at z = ± We searched the QSO spectrum for a possible sodium absorption linked to this galaxy in the strong Na i D lines at 5889 Å and 5895 Å (rest frame). As one arcsecond corresponds to 4.8 pkpc at z = 0.33, the impact parameter is about 26 pkpc, and one may expect a non-negligible absorption if the galaxy is not too compact. We found a faint (2% depth) unresolved line at Å that might correspond to Na i λ5889 (which has the larger oscillator strength), but it is offset by a velocity Article number, page 12 of 16

13 P. L. North et al.: MUSE-inspired view of the quasar Q , its Lyman α blob, and its neighborhood δ " 12" 9" 6" 3" E N SBLy α (10 18 ergs 1 cm 2 arcsec 2 ) Flux Fλ [erg s 1 cm 2 Å 1 ] Bimodal Gaussian fit λ 1 = ± 0.8 Å λ 2 = ± 0.3 Å 0" 0" 3" 6" 9" 12" 15" 18" α Emission wavelength λ em [Å] Fig. 16. Left: Same as Fig. 15, but for LAE2. Right: Same as Fig. 15, but for LAE2. The spectrum is summed over a disk with 1.0 radius. Table 4. Equatorial coordinates and line intensities of the starburst galaxy at z = 0.33 α J2000 [h:mn:s] δ J2000 [ : : ] 21 : 02 : : 53 : Line Observed flux [10 18 erg s 1 cm 2 ] [O ii] H β 14.9 [O iii] [O iii] H α 40.3 [N ii] v = 280 km s 1 with respect to the expected redshifted line. If this feature really were the Na i λ5889 line, then it would correspond to a Na i column density on the order of cm 2. If all Na were neutral and none of it were locked in dust, this would correspond to an H column density of about cm 2 for solar metallicity. These numbers appear reasonable in view of the Na column densities ( cm 2 to cm 2 ) found by Schwartz & Martin (2004) for local dwarf starburst galaxies; however, the non-detection of the Na iλ5895 line makes the above estimates a speculation. We now focus on other properties of this galaxy. The Balmer decrement can be estimated through the ratio H α /H β 2.71: it is slightly lower than the usually adopted value of 2.86, suggesting no significant dust reddening. The [N ii]λ6584/h α and [O iii]λ5007/h β ratios point to a rather typical star-forming galaxy (see e.g. Lamareille et al. 2004, Fig. 1). It has a relatively high excitation: with a ratio [O iii]λ5007/[o ii]λ3727 = 3.46, it resembles the galaxies of the DR7 sample of Stasińska et al. (2015), although the [O iii]λ4363 line is not clearly detected. The [O iii]/[o ii] ratio, although rather high, is lower than that of the LCEs of Verhamme et al. (2017). Its properties are listed in Table 4 and its spectrum is shown in Fig. 17. The oxygen abundance can be estimated from the line strengths given in Table 4 and from the O3N2 log ( [O iii]λ5007/h β [N ii]λ6584/h α ) and N2 log([n ii]λ6584/hα ) criteria (Marino Fig. 17. MUSE spectrum of the starburst galaxy (top) and spectrum of a nearby sky background (bottom). The brightest lines are identified. The very bright [O i] λ telluric line appears, by chance, in the galaxy spectrum but not in the sky spectrum, because it is not subtracted correctly in most parts of the datacube. The inset is the monochromatic MUSE image at Å (redshifted [O iii] λ5007) centered on the QSO, showing the position of the emission line galaxy; the image is 30 across; the bright object near the top is a Galactic red giant star. et al. 2013, Eqs. (2) and (4), respectively). The main uncertainty, about a factor of two, lies in the [N ii]λ6584 line intensity. We obtain 12 + log(o/h) = 8.13 ± 0.07 (O3N2) (15) = 8.17 ± 0.14 (N2), (16) assuming a ±0.3 dex error on both criteria. This corresponds to [O/H] 0.6 assuming a solar oxygen abundance of 8.69 (Allende Prieto et al. 2001), and agrees well with local starburst Article number, page 13 of 16