Issue 
A&A
Volume 608, December 2017



Article Number  A139  
Number of page(s)  10  
Section  Extragalactic astronomy  
DOI  https://doi.org/10.1051/00046361/201731791  
Published online  15 December 2017 
Modeling 237 Lymanα spectra of the MUSEWide survey^{⋆}
Institute of Theoretical Astrophysics, University of Oslo, Postboks 1029, 0315 Oslo, Norway
email: maxbg@astro.uio.no
Received: 17 August 2017
Accepted: 20 September 2017
We compare 237 Lymanα (Lyα) spectra of the MUSEWide survey to a suite of radiative transfer simulations consisting of a central luminous source within a concentric, moving shell of neutral gas, and dust. This six parameter shellmodel has been used numerously in previous studies, however, on significantly smaller datasets. We find that the shellmodel can reproduce the observed spectral shape very well – better than the also common “GaussianminusGaussian” model which we also fitted to the dataset. Specifically, we find that ~ 94% of the fits possess a goodnessoffit value of p(χ^{2}) > 0.1. The large number of spectra allows us to robustly characterize the shellmodel parameter range, and consequently, the spectral shapes typical for realistic spectra. We find that the vast majority of the Lyα spectral shapes require an outflow and only ~ 5% are wellfitted through an inflowing shell. In addition, we find ~ 46% of the spectra to be consistent with a neutral hydrogen column density < 10^{17} cm^{2} – suggestive of a nonnegligible fraction of continuum leakers in the MUSEWide sample. Furthermore, we correlate the spectral against the Lyα halo properties against each other but do not find any strong correlation.
Key words: radiative transfer / line: formation / line: profiles / galaxies: highredshift / scattering
The full Table A.1 is only available at the CDS via anonymous ftp to cdsarc.ustrasbg.fr (130.79.128.5) or via http://cdsarc.ustrasbg.fr/vizbin/qcat?J/A+A/608/A139
© ESO, 2017
1. Introduction
Not many branches of astrophysics show a similar fast transition from theory to datadriven as the study of the highz Universe through Lymanα (Lyα) emitting objects. While predicted already half a century ago (Partridge & Peebles 1967), the first detections of highredshift Lyα emitters was not until more than 25 years later (e.g., Møller & Warren 1993; Rhoads et al. 2000). They were thus preceded by groundbreaking theoretical studies – which still represent the solid foundation of our understanding of Lyα radiative transfer in astrophysical environments today (Osterbrock 1962; Adams 1972; Neufeld 1990). These theoretical studies were often accompanied by numerical radiative transfer codes (e.g., Adams 1972; Bonilha et al. 1979) – with the focus shifting towards the latter more recently (e.g., Dijkstra et al. 2006; Verhamme et al. 2006; Hansen & Oh 2006). While this numerical work together with increased computational power allowed for the study of more complex geometries (such as the output from hydrodynamical simulations) as well as the consideration of a wider parameter space, the predictive power has naturally decreased due to this multiplication of degrees of freedom and the associated degeneracies.
In addition to these difficulties, a somewhat unfortunate incident arguably slowed down the progress on theoretical side considerably, namely the apparent nonrequirement for (more) complex theoretical models in order to explain a wide range of observed Lyα spectra. Specifically, a homogeneous shell of neutral hydrogen (and dust) surrounding a Lyα and continuum emitting source (as first advocated by Ahn et al. 2003) seems to be able to reproduce many observed line shapes from z ~ 0 (Yang et al. 2017) to higher redshifts (Verhamme et al. 2008; Karman et al. 2017). This has been shown using large libraries of precomputed line shapes (Schaerer et al. 2011; Gronke et al. 2015) which can be compared to observations. In spite of this automatizing of the fitting procedure and the consequent common usage of the “shell model”, the physical meaning of its parameters is still unclear. Also, the structure of the shellmodel, that is, the existence of a large solid shell of cold, neutral gas seems somewhat contrived. Nevertheless, the shellmodel has been studied numerously and the effect of its parameters on the spectral shape are understood well (Verhamme et al. 2015). These facts make the shellmodel a common starting point of theoretical modeling of Lyα propagation – and it does so since more than a decade.
This rather slow progress on the theoretical side has not halted the planning and execution of new experiments geared towards Lyα detection and study. As a result, the pure number of detected Lyα emitting objects – even at the highest observable redshifts – has increased impressively (e.g., the latest detection of ~ 2000 Lyα emitters at z ~ 6−7, Ouchi et al. 2017; also see Barnes et al. 2014, for a list of Lyα surveys), and the quality of the detections has improved significantly (e.g., McGreer et al. 2017). This development has not only lead to the discovery of new types of astrophysical objects such as “Lymanα blobs” (e.g., Francis et al. 1996; Fynbo et al. 1999; Steidel et al. 2000; North et al. 2017), “giant Lyα halos” (e.g., Cantalupo et al. 2014; Hennawi et al. 2015; Cai et al. 2016), and “extreme equivalent width objects” (e.g., Sobral et al. 2015) but also opened up the chance of statistical usage of Lyα data, for example, to constrain the “Epoch of Reionization” (see review by Dijkstra 2014) or the leakage of ionizing photons (e.g., Dijkstra et al. 2016b; Verhamme et al. 2017).
The latest stage in this series of technological advances and consequent observational progress has been the implementation of ESO’s “Multi Unit Spectroscopic Explorer” (MUSE; Bacon et al. 2010). Equipped with 24 integral field units, this integral field spectrograph has literally added an additional dimension to Lyα observations and will thus multiply the available Lyα data even further. The unprecedented sensitivity of MUSE has already revealed the ubiquity of Lyα emission throughout the Universe (Bacon et al. 2015; Wisotzki et al. 2016; Drake et al. 2017), and the combination of imaging with spectral information opens up many new possibilities for Lyα data analysis – which will further our understanding of astrophysical Lyα emitting objects.
Recently, the data of the first installment of the MUSEWide survey has been released publicly (Herenz et al. 2017) which contains 237 Lyα spectra around z ~ 4. The sample (corresponding to ~ 20% of the final survey area) is emission lineselected which allows an unbiased spectral analysis. We take this opportunity to assemble a statistical relevant sample of shellmodel fits. Thus far, the number of modeled spectra found in the literature are often merely individual spectra (e.g. Vanzella et al. 2016; Dahle et al. 2016) or compilations up to ~ 10 s (Hashimoto et al. 2015; Karman et al. 2017; Yang et al. 2017). A (much) greater number of spectra will, thus, allow us to address the following questions:

How well can the shellmodel reproduce observed spectral lineshapes? Is there, e.g., a certain spectral characteristic that cannotbe reproduced?

What are common values for shellmodel parameters? Answering this question is useful for studies concerned with a single object and might help them to characterize how unusual their spectral shape is.

Do any correlations amongst the shellmodel or with external parameters exist?
Furthermore, this work will provide the foundational data for future work “decrypting” the shellmodel parameters. In a recent series of papers (Gronke & Dijkstra 2016; Gronke et al. 2016; Gronke et al. 2017), we’ve shown already that the shellmodel parameters do not reflect the physical properties of a simple multiphase model. However, in such a case the spectra do commonly also show more flux a linecenter and are less asymmetric (in spite of outflows) compared to observed lineprofiles. However, in the presence of smallscale structure such as tiny “cloudlets” (as theoretically predicted by McCourt et al. 2016), Lyα photons escape the – still multiphase – system as it was homogeneous which yields more realistic spectra, and makes the obtained shellmodel parameters physically meaningful. In this work, we will interpret the data also in the light of these recent results.
The paper is structured as follows: Sect. 2 describes the MUSEWidefield data set, and the fitting procedure. We present and discuss our results in Sect. 3, and conclude in Sect. 4. When required, we use the Planck Collaboration et al. (2016) parameters of a flat ΛCDM Universe (H_{0} = 67.7 km s^{1} Mpc^{1}, Ω_{m0} = 0.307, Ω_{Λ0} = 0.691), and use “log ” as logarithm with base ten.
2. Method
2.1. The MUSEWide data
The MUSEWide catalogue used in this work was described extensively in Herenz et al. (2017). Here, we summarize the main properties of the data set and explain which preprocessing steps we undertook prior to modeling.
In total, the catalogue contains 831 datasets with a spectral resolution of 2.5 Å out of which 237 contain a Lyα line. These Lyα emitting objects are in a redshift range of z_{MW} = [2.97, 6.00] with a median redshift of 3.83. The redshift z_{MW} given for these objects stems from fitting an asymmetric profile to the Lyα line (see Sect. 3.4.1 in Herenz et al. 2017) in order to obtain the peak wavelength which was converted to a redshift. The quoted redshift error is due to the fitting procedure, i.e., corresponds to an uncertainty in the peak position. This means for both, the redshift and the redshift error, no correction for radiative transfer effects was introduced. During the fitting procedure, we will leave the redshift as a free parameter. In order to use the spectra, we extract a wavelength range which allows for a shift of Δz = [−0.02, 0.01] from the given peakredshift while considering that even with a maximum shift no data point is outside the range of [± 2500 km s^{1}] which is supported by the fitting pipeline. The asymmetric Δz range is due to the fact that most Lyα peaks are redwards of nonresonant lines, i.e., we expect in most cases a z value blueward of z_{MW} (i.e., z<z_{MW}). However, we still leave enough room (Δz = 0.01 corresponds to ~ 750 km s^{1} at z = 3 – a value larger than commonly found shifts of several ~ 100 km s^{1}, Shapley et al. 2003; Steidel et al. 2010; Kulas et al. 2012) for an opposite shift. This procedure ensures that even while shifted the same number of data points is compared to the theoretical spectrum.
As some spectral data contain negative mean continuum fluxes – which are unphysical – we take a conservative approach, and remove this continuum level from all the spectra. This implies, that we will only use the Lyα line shape to fit theoretical spectra to it and ignore the surrounding continuum.
Another important property of the Herenz et al. (2017) dataset used in this work include the “confidence” category which can be either 3 for spectra which are almost certainly Lyα (~ 40% of the spectra), 2 for nearcertainty (~ 52%), and 1 for unsure detections which are, however, not spurious (~ 8%). We refer to Herenz et al. (2017) for more details about the classification.
2.2. Shellmodel fitting
The shellmodel fits were carried out via an improved version of the fitting procedure described in Gronke et al. (2015). Because the shellmodel consists of a central source emitting Lyα photons inside a moving shell of neutral hydrogen and dust, it can be parametrized through the HI column density N_{HI}, the (allabsorbing) dust optical depth τ_{d}, the outflow velocity v_{exp} (< 0 for infall), the intrinsic width of the Lyα line σ_{i}, and the effective temperature of the gas T.
In order to fit observed spectra by this model, one usually divides these parameters into sets of discrete and continuous ones. While the parameter space spanned by the discrete ones has to be covered through (computationally relatively expensive) radiative transfer simulations, the continuous ones can be modeled in postprocessing. We model the effect of σ_{i} and τ_{d} through individual weighting of photon packages (see Gronke et al. 2015, for details) with an allowed range of σ_{i} ∈ [1, 800] km s^{1} and τ_{d} ∈ [0, 5]. The remaining subspace spanned by (v_{exp}, log N_{HI}/ cm^{2}, log T/ K) we cover via a (54, 30, 8)grid which totals 12 960 discrete models. Each of the models was carried out using 20 000 “Lyα” photons with Gaussian intrinsic spectrum centered at 0 and a standard deviation of 800 km s^{1} and 150 000 “UV” photons with uniform intrinsic spectrum in the range ± 2700 km s^{1}.
While the parameters log N_{HI}/ cm^{2} and log T/ K are spaced uniformly between their respective limits [16, 21.8] and [3, 5.8], the outflow velocity v_{exp} can take the additional values of 2 km s^{1}, 5 km s^{1}, and 8 km s^{1} apart from the values with 10 km s^{1} separation between 0 and 490 km s^{1}. This is due to the strong effects small velocities can have on the resulting spectrum due to the peaked nature of the Lyα scattering crosssection. Furthermore, we use the symmetry of Lyα radiative transfer to cover inflow velocities (Dijkstra 2017).
In summary, the following steps were carried out for each spectrum:

1.
We fix a redshift prior of with meanz_{MW} − 0.005 and standard deviation of σ_{z} = 0.005. As stated above, the offset is due to the more likely observation of the red peak of the Lyα line. This prior is very wide and easily overruled by the data.

2.
Due to the multimodal nature of the likelihood landscape – with many local peaks which are, however, of dramatically different altitude – it is crucial to find the global likelihood maximum around which one can carry out the MonteCarlo optimization. In a continuous parameter space with unknown likelihood function it is impossible to know whether or not this global maximization is successful. However, we do our best by combining a basinhopping algorithm over the discrete parameters with a local maximization over the continuous parameters. This means, we perform a basinhopping algorithm with randomly chosen starting position until no better results is obtained^{1}. Each step consists of a random displacement in parameter space followed by a local maximization using the robust POWELL algorithm. In order to speed up the global maximization process by reducing the memory access, we execute the basinhopping over a local maximization using only the continuous parameters (z, τ_{d}, σ_{i}). For this process, we use the COBYLA method (Powell 1994).

3.
From this initial position, we start the affine invariant MonteCarlo sampler emcee (Goodman & Weare 2010; ForemanMackey et al. 2013) with 200 walkers and 1240 steps. Specifically, we place the walkers on several positions identified through the optimization process weighted by their likelihood. In case a betterfitting spectrum (Δlnp ≥ 0.5) than the starting one is found, the MC sampling is restarted from this position (with 1200 steps).

4.
For further analysis, we discard the first “burnin” steps for which the mean (over the walkers) of each parameter does not disagree by more than 15% off the final value. From the remaining chains we compute the median values of each shellmodel parameter and its uncertainties.
Fig. 1 Signaltonoise ratio (S/N) of the Lyα line versus the quality of the fit (p(χ^{2}) computed as described in Sect. 3.1). The color coding denotes the “confidence” category as defined in Herenz et al. (2017). 
Fig. 2 The four spectra with the lowest p(χ^{2})values of (from left to rightmost panel) { < 0.1, 0.5,6,14 } × 10^{3}. The data is shown in darkblue, and the best fit shellmodel spectrum with a solid red line (the dotted red line is the intrinsic spectrum). In addition, the semitransparent black lines show 10 randomly drawn spectra from the burnedin chains. Apart from the spectrum in the third panel which has been graded with “confidence” level 2 the other three spectra have been classified as “3”. 
Fig. 3 Uncertainties on the obtained fitting parameters as a function of the S/N. Specifically, we show half of the 84th minus the 16th percentile of the corresponding posterior distributions. The color corresponds to the “confidence” category assigned to each spectrum by Herenz et al. (2017) with red for the highest confidence (3), followed by green and blue for the lowest confidence class. 
3. Results
3.1. Fit quality
In order to assess the quality of the shellmodel fits, we computed the p(χ^{2}) values of the bestfit spectra. Following Pearson (1900) we compute (where σ_{i} is the flux error in bin i and the sum is taken over the N data points of the spectrum), and quote as p(χ^{2}) the cumulative probability of a χ^{2}distribution with N degrees of freedom having values greater than χ^{2}.
Figure 1 shows the distribution of the resulting p(χ^{2}) values as a function of the S/N of the Lyα line. Also shown in Fig. 1 is the “confidence” class of the spectrum (as discussed in Sect. 2.1). One can notice that the p(χ^{2}) values are roughly uniformly distributed between the extremes of 3 × 10^{13} and 0.994 with a median of 0.773. Overall, the fit quality is very good with 94% [73%] of the spectra with a p(χ^{2}) > 0.1 [> 0.5].
As visible from Fig. 1, there seems to be no strong overall correlation neither versus the S/N (Spearman’s rank correlation coefficient of r_{S} ≈ −0.08, with a twosided pvalue of a hypothesis test to check whether the two variables are correlated p_{S} ≈ 0.24) nor with any of the (median) shellmodel parameters (p_{S}> 10^{3}) – with the exception of σ_{i} where r_{S} ≈ −0.26 (p_{S} ≈ 4 × 10^{5}). This tentative anticorrelation between σ_{i} and the shellmodel fit quality might be suggestive of a stronger influence of radiativetransfer for spectra with lower σ_{i}.
Figure 2 shows the four worst fits, i.e., with the lowest p(χ^{2}) values. It is noteworthy that even for those spectra the fit does not “look that bad”, i.e., the model does not seem to fail catastrophically. The large χ^{2} values in the examples shown in Fig. 2 can be attributed to slightly lower (negative) continuum on the blue side of the emission (two leftmost panels), or some (again with negative intensity) outliers (two rightmost panels). When increasing the error bars only by 10% (20%), the overall p(χ^{2}) increases from to ()^{2}. This would not be the case if the theoretical model cannot reproduce certain features which are present in the observed spectra – such as an additional peak – and lets us conclude that overall, the shellmodel seems to provide (surprisingly) good fits to all the 237 Lyα spectra.
Fig. 4 Correlation matrix between the median obtained shellmodel parameters. The color coding denotes the “confidence” category the spectrum was assigned (blue = 1, green = 2, red = 3). While the lower triangle shows scatter plots with each spectrum marked individually, the same data is plotted as contours of kerneldensity estimators in the upper right triangle. The panels along the diagonal show stacked histograms of the (projected) distribution. 
3.2. Fitting results
Through the MonteCarlo sampling, we obtain an estimate of the full posterior distribution which we quantify through one and twodimensional projections, and the 16th, 50th and 84th percentiles. Due to the large number of data points, we cannot show all the three percentiles (i.e., the median and the uncertainty on the parameters) in the same plot. Therefore, Fig. 3 shows first the uncertainty on the fitting parameters which we quantified through (where q_{X} it the Xth percentile) as a function of the S/N. Also, the color coding shows the “confidence” group with red the highest and blue the lowest confidence spectra. Clearly, with improved spectral quality, the fitting uncertainty reduces. For instance, in order to obtain a ≲ 0.5 dex uncertainty on the column density, the S/N should be ≳ 20. However, as is also apparent by comparing the various panels of Fig. 3 the improvement is not equally strong for the six parameters. For instance, the dustoptical depth shows mostly (with a few exceptions) a very large uncertainty (≳ 1), and the effective temperature can only be estimated within an order of magnitude (see also Gronke et al. 2015, for similar results using mock spectra).
Note that a while a greater S/N will improve the constraints on the shellmodel parameters, the increased spectral resolution will – most likely – not help significantly. The MUSE spectral resolution of 2.5 Å corresponds to a resolving power of R ~ [1900, 3400] for the Lyα line in the considered redshiftrange. However, above R ~ 1000 the uncertainty does not depend significantly on the spectral resolution (see Appendix A in Gronke et al. 2015). Better resolution spectra provide, nevertheless, important additional information as potentially features arise which cannot be captured by the shellmodel (as, e.g., the extended redwing in McGreer et al. 2017; or the reversed absorption draught in Yang et al. 2016) let us peek beyond its simplicity. They might, thus, give insights into the physical meaning of the shellmodel parameters^{3}.
Figure 4 shows the one and twodimensional projections of the distribution of the obtained median shellmodel parameters. Here, we again split the dataset according to the “confidence” group provided. No strong correlation is visible between the parameters. A slight tendency for more shift (which we show in velocity space as Δv; see below) for larger HI column densities or larger σ_{i} are the only (weak) correlations we could identify.
The distribution of each individual parameter might help to understand what the shellmodel parameters mean physically. Here we discuss them onebyone:

The neutral hydrogen column density is near lognormallydistributed (d’Agostino & Pearson 1973, test ofnormalization yieldsp_{AP} ≈ 0.14) with mean 18.78 and standard deviation 0.88. This means the limits on this parameter are sufficiently far away, i.e., a wider covering of the parameter space would most probable not help to obtain better fits. Interestingly, ~ 46% of the spectral shapes [~ 36% of the “category 3” spectra] are consistent (within 1  σ) with a column density of N_{HI}< 10^{17} cm^{2}. This could indicate a potentially large number of Lymancontinuum leakers amongst the MUSEWide galaxies. However, better quality spectra are required (in particular, for the fainter galaxies) in order to draw any conclusions. With the current data set only one object clearly prefers such a low column density, that is, the 84th percentile is lower than the threshold. We would also like to caution that even if this suspicion substantiates with better quality data, i.e., also then the Lyα spectrum points towards N_{HI} < 10^{17} cm^{2}, it is not clear if these objects would be detected as continuum leakers due to a potential (non)directional dependence of the Lyα spectrum (Dijkstra et al. 2016a; Eide et al., in prep.; but see Verhamme et al. 2012).

The “effective temperature” of the shell obtains most likely takes values of T ~ 10^{4}K (median log (T/ K) ~ 4.38) which is close to what is expected of neutral gas. However, also larger (and smaller) values are found – and, in particular temperatures of T ~ 10^{5}K which are thermally unstable. If the shellmodel temperature had some physical meaning, these values could be explained with subgrid turbulence which enters as T_{eff} = T_{0} + ⟨ v^{2} ⟩ m/k_{B}, i.e., dispersion velocities of would correspond to an effective temperature of T ~ 10^{5}K. Note also that mostly ~ 1 dex uncertainty is associated with the effective temperature (cf. Fig. 3).

The width of the intrinsic spectrum σ_{i} occupies the range . The distribution is furthermore faster declining towards smaller values of σ_{i} which makes the canonical value of σ_{i} ~ 13 km s^{1} – corresponding to thermal motion of a gas with T ~ 10^{4}K – unlikely. However, both turbulent motion of the emitting gas (in the case of, e.g., Lyα production through fluorescence, Hogan & Weymann 1987; MasRibas & Dijkstra 2016) and of cold gas relatively nearby the emitting source will lead to a broadening of the Lyα line prior to additional radiative transfer effects, and these larger σ_{i} values are relatively easy to justify. Interestingly, although even broader intrinsic lines of σ_{i} ≳ 400 km s^{1} are allowed by the fitting pipeline (up to σ_{i} ≤ 800 km s^{1}), they are not favored by the data.

The resulting (absorbing) dust optical depth of the shellmodel is centered around . Within the geometry of the homogeneous shell, these values would correspond to extremely low escape fractions of which seem highly unlikely for the considered objects. Furthermore, if converted to a metallicity (using the SMC conversion factor of Z/Z_{⊙} = 9.3 × 10^{20} cm^{2}/N_{HI}τ_{d} which assumes an albedo of A = 0.32, Pei 1992; Li & Draine 2001; Laursen et al. 2009) the resulting metallicities of are unrealistically high. Therefore, the literal interpretation of the dust optical depth is not very probable. However, as the uncertainty on the τ_{d} is mostly greater than unity (compare Fig. 3), lower dust contents are consistent with the spectral shape. The reason for these preferred large τ_{d} values can be understood when considering of how an increased dust optical depth changes the emergent Lyα spectrum. A common misconception – probably due to the (near)frequency independence of the dust crosssection around λ ≈ 1216 Å (Pei 1992; Laursen et al. 2009) – is that the effect of dust does not change the Lyα spectrum (but merely the Lyα escape fraction). While this assumption holds for many cases, it does not universally. This is because dust attenuates the photons proportional to their path length through cold, dusty medium, and thus, if there is a correlation between emergent frequency and trajectory length, the shape of the spectrum will be changed. For instance, in “usual” doublepeaked Lyα profiles through outflowing media (with an enhanced red side), the Lyα photons forming the blue peak have travelled further through neutral hydrogen. Therefore, an increased dust content will lower the blue peak and increase the spectral asymmetry. This might be the reason for the high τ_{d} preference of the data. However, for the observed spectra the blue side might be attenuated not by dust but by neutral gas in the CGM and/or IGM which scatters Lyα photons out of the lineofsight, and thus, effectively absorb the blue side of the spectrum (Dijkstra et al. 2007; Laursen et al. 2011).

The shell outflow velocities (i.e., v_{exp}> 0) are approximately normally distributed with mean ~ 211 km s^{1} and standard deviation ~ 94 km s^{1}. For lower quality spectra (categories 1 and 2), we found also inflow velocities with similar magnitude. The values of v_{exp} thus agree with previous measured outflow velocities in the nearby (Kunth et al. 1998; RiveraThorsen et al. 2015) and distant Universe (Shapley et al. 2003; Kulas et al. 2012). Overall, we find (within 1  σ) only ~ 5% of the spectra to be better fit by inflow velocities [~ 3% of the confidence 3 spectra] and ~ 25% [~ 9%] to be consistent with inflows. These limits could potentially constrain the covering fraction of cold, inflowing gas and thus yield insights into the gas budget of highz galaxies. We caution, however, that the directional dependence of Lyα spectra (in such scenarios) is yet unclear (see Fig. D.1 in Trebitsch et al. 2016, for an angleaveraged spectrum of an filamentary accreting object which shows signs of inflows).
Fig. 5 Extent of the Lyα halo and it’s brightness versus the six obtained shellmodel parameters. The upper and central rows show the Kron radius (in kpc) and the line flux (with ) as a function of the fitting parameters (see Sect. 3.3 for more details). The lower row shows the flux ratio between 3 × and 1 × R_{Kron}. The color coding corresponds again to the “confidence” category (blue, green, red in increasing order). Also shown in each panel are the Spearman correlation coefficients.

While not technically being a shellmodel parameter, the redshift of the source is crucial for the fit. We display the redshift in Fig. 4 as shift between the Lyα peak z_{MW}(as provided by Herenz et al. 2017) and the modeled intrinsic redshift z. Furthermore, we convert the difference to velocity space resulting in (1)We find Δv values of ~ (~ for only the highest “confidence” spectra), which implies an emission blueward of the main spectral peak for the majority of analyzed spectra – as commonly found for other Lyα emitting objects (e.g., Erb et al. 2014; Trainor et al. 2015). Only ~ 3.8% (~ 4.2% of the “confidence” 3 spectra) of the spectra show Δv> 0 (above 1  σ). On the other hand, ~ 49% of the spectra (~ 29% of the “confidence” three spectra) are consistent with Δz ~ 0 (within 1  σ). Another potentially interesting point is the correlation between N_{HI} and Δz, i.e., larger shifts for higher column densities which follows directly from the frequency diffusion in the models. Observationally, it has been established that this shift is anticorrelated with the Lyα equivalent width (Erb et al. 2014) which invites to the interpretation of facilitated Lyα escape in systems with less frequency diffusion.
Note that we did not find an evolution of any of the above discussed parameters with redshift. We found, however, that only spectra with high S/N (≳ 20) occupy the parameter space with Δv ≲ 200 km s^{1} or N_{HI} ≳ 10^{19.5} cm^{2} – similarly to what we discussed above for v_{exp}. This is likely due to the fact that the spectral features required for, e.g., a large column density (for instance, the extended wing) is only apparent for higher quality spectra.
All the above described shellmodel parameters, that is, the 16th, 50th, and 84th percentiles of the posterior distributions for each object are publicly available^{4}.
3.3. Correlations with halo properties
The combined modeling of both the spectrum and surface brightness distribution is a crucial next step to understand how Lyα photons escape galaxies and what information they obtain on their trajectory as using both constraints breaks degeneracies which otherwise exist. Integral field spectrograph such as MUSE are the perfect tool to study these two angles simultaneously. As a first step towards the full joined modeling, we want to correlate the spectral with spatial information.
Included in the MUSEWide catalogue is the Kron (1980) radius R_{Kron} of each object as calculated by LSDcat (Herenz & Wisotzki 2017) as well as the flux within { 1,2,3 } × R_{Kron}. LSDcat calculates R_{Kron} following the twodimensional definition of Bertin & Arnouts (1996) which reads (2)where the summations are carried out over the 2D pseudonarrowband images created from the 3D MUSE datacubes.
Figure 5 shows the R_{Kron} (upper row), the flux within R_{Kron} (central row), and the ratio of the flux within 3R_{Kron} and R_{Kron} (lower row) versus the six shellmodel parameters. The latter quantity is a measure how steep the surface brightness profile is falling off. In each panel, we also display the Spearman correlation coefficient. We cannot identify a clear correlation between the fitting parameters and the halo properties. However, some tentative correlations might be apparent. For instance, σ_{i} and the halo extent as well as the brightness of the source – in particular the high quality spectra (shown as red points) are somewhat correlated. Also Δv and the flux within R_{Kron} show some level of anticorrelation. In order to confirm or rule out these tentative correlations (as well as the ones discussed in Sect. 3.2), the uncertainty of the shellmodel parameters should be lower than the correlation range of interest. This mean a similar number of spectra with better signaltonoise (S/N ≳ 20; cf. Fig. 3) would be most useful – a goal which can be achieved in the near future given that this first release only covers about one fifth of the total survey area, and the “MUSEWide” scanning strategy is very shallow (one hour per field Herenz et al. 2017).
3.4. Contrasting the “shell” with the “GaussianminusGaussian” model
As an alternative to the radiativetransfer based shellmodel, we fit the 237 Lyα spectra using the a “GaussianminusGaussian” model (which has been previously used in the literature modeling Lyα spectra; e.g., Jensen et al. 2013). This model also consists of six free parameters, namely the location (i.e., the mean μ) and width (the standard deviation σ) of two normal distributions which have been scaled by an arbitrary amplitude (the total area under the curve A). Please note that in concordance with the name, we do not allow negative amplitudes (which could correspond to a “GaussianplusGaussian” model), and also don’t allow for negative fluxes, that is, we set the data point which would be negative to zero.
Fig. 6 Quality of the fits of the “shell” versus the “GaussianminusGaussian” models. The color stands for the three “confidence” categories. Also included (as a “guide to the eye”) are linear fits to the data, and the identity function (as black dashed line). 
Fig. 7 Example spectral fits of the “shell” and “GaussianminusGaussian” models with all combinations of fitting failures/successes. The blue lines show the data, the red solid line the shell model fit (with the intrinsic spectrum as dotted red line), and the green dashed line shows the “GaussianminusGaussian” best fit. The numbers in the panel show the p(χ^{2}) values of the best fits in the corresponding colors. 
Figure 6 shows a comparison between the fit qualities of the two models using the p(χ^{2}) values of the bestfitting parameters. The color coding stands again for the “confidence” category assigned to each spectrum. Clearly, the shellmodel can reproduce most spectra better than the (even simpler – but with the same number of free parameters) “GaussianminusGaussian” model leading to higher p(χ^{2}) values for ~ 9% of the spectra^{5}. In particular, in some cases the “GaussianminusGaussian” model fails to reproduce the observed spectral shape completely (i.e., p(χ^{2}) ≲ 0.1) while the shellmodel can fit same spectrum successfully. Overall, this improvement of the fit quality seems slightly larger for higher quality spectra (that is, higher “confidence” values); however, this trend is not statistically significant with the current data set. On the other hand, spectra that cannot be modeled using the shellmodel can neither be reproduced using the “GaussianminusGaussian” fit. Instead, both models yield comparable low p(χ^{2}) values.
Figure 7 illustrates the reason for superiority of the “shell” over the “GaussianminusGaussian”model. Each panel shows an observed Lyα spectrum (in blue), the best shellmodel fit (in solid red with dotted red the intrinsic spectrum), and (in dashed green) the best “GaussianminusGaussian” fit. The two numbers in the top left and right corner of each panel show the p(χ^{2}) values of the bestfit shellmodel and “GaussianminusGaussian” fit, respectively. This means, the leftmost panel of Fig. 7 shows an example where the shellmodel yields a much better quality fit. The reason here is clearly the double peaked nature which cannot be captured by the “GaussianminusGaussian” model. The second panel from the left shows the same spectrum as in Fig. 2, i.e., a case where the shellmodel cannot reproduce the observed spectrum well. However, fitting two Gaussian curves yields a similar fit and, thus, also p(χ^{2}) < 0.01. As mentioned in Sect. 3.1 the reason for these low p(χ^{2}) values in this case are the small error bars on the data combined with the slight negative flux bluewards of the line. It is noteworthy that such wide, relatively symmetric spectra can naturally be well reproduced using a simple Gaussian curve. This means per construction the “GaussianminusGaussian” model suits this spectral type well. However, also the shellmodel can feature such a spectrum – given a low optical depth of the shell (through low column density and high outflow velocity) and a wide intrinsic spectrum. This will lead to hardly any radiative transfer effects, and thus, to a Gaussian emergent spectrum.
The third panel (from the left) in Fig. 7 displays an singlepeaked, asymmetric spectrum. This is an example where both models yield good fits to the data. This is also the case in the right panel of Fig. 7 (which we chose because it is the spectrum where the p(χ^{2}) difference is greatest in favour for the “GaussianminusGaussian” model). Interestingly, here the two normal distributions can also reproduce the two peaks of the data. This is due to the lower blue peak (compared to the first panel) and the larger error within it.
A potential physical interpretation of the “GaussianminusGaussian” model is some intrinsic spectrum (e.g., shaped by within the interstellar medium) which is then processed by the circum and/or intergalactic medium (see, e.g., Gronke & Dijkstra 2016; and Dijkstra et al. 2007, respectively). The latter step would lead to a “Lyα halo” surrounding the galaxy and – as (with increasing distance from the galaxy) most scatterings are out of the lineofsight – can be modeled as an effective absorption. This would mean that an increased level of absorption should correlated positively with the halo size. However, we do not find any correlation of the “absorbed flux”^{6} (r_{S} ≈ −0.07) with R_{Kron}. Neither we find any correlation with any other halo property discussed in Sect. 3.3. As both the intrinsic spectrum as well as the absorption feature may differ from a normal distribution – which makes the “GaussianminusGaussian” model a simplification – this result is maybe not too surprising.
One can surely extend the “GaussianminusGaussian” model with additional free parameters (by, e.g., allowing for negative amplitudes and nonzero skewness) in order to provide better fits to the observed Lyα spectra. However, one should keep in mind that the advantage of a radiativetransfer based model (such as the shellmodel) is that its spectral shapes are at least a subset of all the possible Lyα spectra. This is not necessarily so for “more artificial” models. Also, these models to not directly alter the spectral shape by changing the frequency diffusion process, and, thus it physical content of the fitting parameters is even more questionable than in the shellmodel.
4. Conclusions
We fitted the publicly available Lyα spectra of the “MUSEWide survey” published in Herenz et al. (2017) using the shellmodel. This dataset contains 237 spectra around redshift z ~ 3.8. We carried out the fitting procedure using the pipeline described in (Gronke et al. 2015) which yields robustly the bestfit parameters as well as the respective confidence intervals. Our main findings can be described as follows:

Overall, the shellmodel provides excellent fits to a majority of thespectra (with > 90% possessing a p(χ^{2}) > 0.1; see Sect. 3.1) with a better fits than the ‘GaussianminusGaussian’ model (Sect. 3.4). Even the worst shellmodel fits can reproduce the spectral shape reasonably well (cf. Fig. 2). This implies that this simple model can practically fit all 237 spectra.

We identified common parameter ranges which are typical for observed spectra. This allows future studies to generate more realistic Lyα spectra, and to compare newly obtained shellmodel parameters to this distribution. For this purpose, we make all the fits publicly available^{7}. With the exception of the shell’s dust optical depth we find none of the parameters completely unphysical (see Sect. 3.2), i.e., within the range what has been derived by other means. For instance, we recovered from only the spectral shape outflows in the range of ~ 200 km s^{1} (in particular, the higher quality spectra) which has been found using metalabsorption lines (e.g., Steidel et al. 2010).

Two recovered parameter distributions (which can have physical meaning in a multiphase medium with many clumps per sightline; see Gronke et al. 2017) are of particular interest: (i) the majority of spectra (~ 70%) are better fit through an outflowing shell and only ~ 5% are prefer an inflowing shell. Given firmer underlying data this could yield insights into the gas budget of highredshift galaxies; (ii) half of the spectra can be fit with column densities of N_{HI} ≲ 10^{17} cm^{2}. This suggests a possible large fraction of Lymancontinuum leakers in the MUSEWide data set and supports the importance of Lyα emitting galaxies on the ionizing budget (Dressler et al. 2015; Verhamme et al. 2015; Dijkstra et al. 2016b).

We did not find any strong correlation neither amongst the shellmodel parameters nor with the Lyα halo properties. However, a weak correlation between Δv and N_{HI}, and a tentative correlation between σ_{i} and Δv exist.
This work serves also as a “proofofconcept”, i.e., that the modeling of a large set of observed Lyα spectra – as will be common in the near future – is possible. A similar analysis can be improved by increasing the S/N of the spectra (cf. Fig. 3), and by providing priors on the systemic redshift (thus, confining Δv to values close to zero).
Naturally, if such spectral features are present they would lower the fitquality of the shellmodel discussed in Sect. 3.1. However, this has to be tested with higherresolution data as the current dataset seems to be well fit by the shellmodel.
The full list is available at the CDS (see Appendix A). It is also available under http://bit.ly/aspectraofMW
The individual parameters are shown in Table A.1. It is also available at http://bit.ly/aspectraofMW.
Acknowledgments
I thank the anonymous referee for constructive and helpful feedback. I am grateful to my PhD supervisor Mark Dijkstra for his guidance and comments on the draft. Furthermore, I thank Edmund Christian Herenz for a careful read of the draft. Also, I would like to thank him and the entire MUSE team for making the data used in this work publicly available. Finally, I wish to thank Lutz Wisotzki for inviting me to the Leibnitz Institute for Astrophysics in Potsdam. This research made use of a number of open source software such as the Python programming language and its packages astropy (Astropy Collaboration et al. 2013), IPython (Pérez & Granger 2007), matplotlib (Hunter 2007), seaborn (Waskom et al. 2014), and SciPy (Jones et al. 2001).
References
 Adams, T. F. 1972, ApJ, 174, 439 [NASA ADS] [CrossRef] [Google Scholar]
 Ahn, S.H., Lee, H.W., & Lee, H. M. 2003, MNRAS, 340, 863 [NASA ADS] [CrossRef] [Google Scholar]
 Astropy Collaboration, Robitaille, T. P., Tollerud, E. J., et al. 2013, A&A, 558, A33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bacon, R., Accardo, M., Adjali, L., et al. 2010, in Groundbased and Airborne Instrumentation for Astronomy III, Proc. SPIE, 7735, 773508 [Google Scholar]
 Bacon, R., Brinchmann, J., Richard, J., et al. 2015, A&A, 575, A75 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Barnes, L. A., Garel, T., & Kacprzak, G. G. 2014, PASP, 126, 969 [NASA ADS] [CrossRef] [Google Scholar]
 Bertin, E., & Arnouts, S. 1996, A&AS, 117, 393 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bonilha, J. R. M., Ferch, R., Salpeter, E. E., Slater, G., & Noerdlinger, P. D. 1979, ApJ, 233, 649 [NASA ADS] [CrossRef] [Google Scholar]
 Cai, Z., Fan, X., Yang, Y., et al. 2016, ApJ, 837, 1 [Google Scholar]
 Cantalupo, S., ArrigoniBattaia, F., Prochaska, J. X., Hennawi, J. F., & Madau, P. 2014, Nature, 506, 63 [NASA ADS] [CrossRef] [Google Scholar]
 d’Agostino, R., & Pearson, E. S. 1973, Biometrika, 60, 613 [Google Scholar]
 Dahle, H., Aghanim, N., Guennou, L., et al. 2016, A&A, 590, L4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Dijkstra, M. 2014, PASA, 31, e040 [Google Scholar]
 Dijkstra, M. 2017, ArXiv eprints [arXiv:1704.03416] [Google Scholar]
 Dijkstra, M., Haiman, Z., & Spaans, M. 2006, ApJ, 649, 14 [NASA ADS] [CrossRef] [Google Scholar]
 Dijkstra, M., Lidz, A., & Wyithe, J. S. B. 2007, MNRAS, 377, 1175 [NASA ADS] [CrossRef] [Google Scholar]
 Dijkstra, M., Gronke, M., & Sobral, D. 2016a, ApJ, 823, 74 [NASA ADS] [CrossRef] [Google Scholar]
 Dijkstra, M., Gronke, M., & Venkatesan, A. 2016b, ApJ, 828, 71 [NASA ADS] [CrossRef] [Google Scholar]
 Drake, A. B., Guiderdoni, B., Blaizot, J., et al. 2017, MNRAS, 471, 267 [NASA ADS] [CrossRef] [Google Scholar]
 Dressler, A., Henry, A., Martin, C. L., et al. 2015, ApJ, 806, 19 [NASA ADS] [CrossRef] [Google Scholar]
 Erb, D. K., Steidel, C. C., Trainor, R. F., et al. 2014, ApJ, 795, 33 [NASA ADS] [CrossRef] [Google Scholar]
 ForemanMackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 125, 306 [CrossRef] [Google Scholar]
 Francis, P. J., Woodgate, B. E., Warren, S. J., et al. 1996, ApJ, 457, 490 [NASA ADS] [CrossRef] [Google Scholar]
 Fynbo, J. U., Møller, P., & Warren, S. J. 1999, MNRAS, 305, 849 [NASA ADS] [CrossRef] [Google Scholar]
 Goodman, J., & Weare, J. 2010, Commun. Appl. Math. Comput. Sci., 5, 65 [Google Scholar]
 Gronke, M., & Dijkstra, M. 2016, ApJ, 826, 14 [NASA ADS] [CrossRef] [Google Scholar]
 Gronke, M., Bull, P., & Dijkstra, M. 2015, ApJ, 812, 123 [NASA ADS] [CrossRef] [Google Scholar]
 Gronke, M., Dijkstra, M., Mccourt, M., & Peng Oh, S. 2016, ApJ, 833, 1 [Google Scholar]
 Gronke, M., Dijkstra, M., McCourt, M., & Oh, S. P. 2017, A&A, 607, A71 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Hansen, M., & Oh, S. P. 2006, MNRAS, 367, 979 [NASA ADS] [CrossRef] [Google Scholar]
 Hashimoto, T., Verhamme, A., Ouchi, M., et al. 2015, ApJ, 812, 157 [NASA ADS] [CrossRef] [Google Scholar]
 Hennawi, J. F., Prochaska, J. X., Cantalupo, S., & ArrigoniBattaia, F. 2015, Science, 348, 779 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Herenz, E. C., & Wisotzki, L. 2017, A&A, 602, A111 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Herenz, E. C., Urrutia, T., Wisotzki, L., et al. 2017, A&A, 606, A12 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Hogan, C. J., & Weymann, R. J. 1987, MNRAS, 225, 1P [NASA ADS] [CrossRef] [Google Scholar]
 Hunter, J. D. 2007, Comput. Sci. Engin., 9, 90 [Google Scholar]
 Jensen, H., Laursen, P., Mellema, G., et al. 2013, MNRAS, 428, 1366 [NASA ADS] [CrossRef] [Google Scholar]
 Jones, E., Oliphant, T., Peterson, P., et al. 2001, SciPy: Open source scientific tools for Python, http://www.scipy.org/ [Google Scholar]
 Karman, W., Caputi, K. I., Caminha, G. B., et al. 2017, A&A, 599, A28 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Kron, R. G. 1980, ApJS, 43, 305 [NASA ADS] [CrossRef] [Google Scholar]
 Kulas, K. R., Shapley, A. E., Kollmeier, J. A., et al. 2012, ApJ, 745, 33 [NASA ADS] [CrossRef] [Google Scholar]
 Kunth, D., MasHesse, J. M., Terlevich, E., et al. 1998, A&A, 334, 11 [NASA ADS] [Google Scholar]
 Laursen, P., SommerLarsen, J., & Andersen, A. C. 2009, ApJ, 704, 1640 [NASA ADS] [CrossRef] [Google Scholar]
 Laursen, P., SommerLarsen, J., & Razoumov, A. O. 2011, ApJ, 728, 52 [NASA ADS] [CrossRef] [Google Scholar]
 Li, A., & Draine, B. T. 2001, ApJ, 554, 778 [Google Scholar]
 MasRibas, L., & Dijkstra, M. 2016, ApJ, 822, 1 [Google Scholar]
 McCourt, M., Oh, S. P., O’Leary, R. M., & Madigan, A.M. 2016, ArXiv eprints [arXiv:1610.01164] [Google Scholar]
 McGreer, I. D., Clément, B., Mainali, R., et al. 2017, MNRAS, submitted [arXiv:1706.09428] [Google Scholar]
 Møller, P., & Warren, S. J. 1993, A&A, 270, 43 [NASA ADS] [CrossRef] [Google Scholar]
 Neufeld, D. A. 1990, ApJ, 350, 216 [NASA ADS] [CrossRef] [Google Scholar]
 North, P. L., Marino, R. A., Gorgoni, C., et al. 2017, A&A, 604, A23 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Osterbrock, D. E. 1962, ApJ, 135, 195 [NASA ADS] [CrossRef] [Google Scholar]
 Ouchi, M., Harikane, Y., Shibuya, T., et al. 2017, PASJ, submitted [arXiv:1704.07455] [Google Scholar]
 Partridge, R. B., & Peebles, P. J. E. 1967, ApJ, 147, 868 [NASA ADS] [CrossRef] [Google Scholar]
 Pearson, K. 1900, The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science, 50, 157 [Google Scholar]
 Pei, Y. C. 1992, ApJ, 395, 130 [NASA ADS] [CrossRef] [Google Scholar]
 Pérez, F., & Granger, B. E. 2007, Comput. Sci. Engin., 9, 21 [Google Scholar]
 Planck Collaboration, Ade, P. A. R., Aghanim, N., et al. 2016, A&A, 594, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Powell, M. J. 1994, in Advances in Optimization and Numerical Analysis (Springer), 51 [Google Scholar]
 Rhoads, J. E., Malhotra, S., Dey, A., et al. 2000, ApJ, 545, L85 [NASA ADS] [CrossRef] [Google Scholar]
 RiveraThorsen, T. E., Hayes, M., Östlin, G., et al. 2015, ApJ, 805, 14 [NASA ADS] [CrossRef] [Google Scholar]
 Schaerer, D., Hayes, M., Verhamme, A., & Teyssier, R. 2011, A&A, 531, A12 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Shapley, A. E., Steidel, C. C., Pettini, M., & Adelberger, K. L. 2003, ApJ, 58, 65 [NASA ADS] [CrossRef] [Google Scholar]
 Sobral, D., Matthee, J., Darvish, B., et al. 2015, ApJ, 808, 139 [NASA ADS] [CrossRef] [Google Scholar]
 Steidel, C. C., Adelberger, K. L., Shapley, A. E., et al. 2000, ApJ, 532, 170 [NASA ADS] [CrossRef] [Google Scholar]
 Steidel, C. C., Erb, D. K., Shapley, A. E., et al. 2010, ApJ, 717, 289 [NASA ADS] [CrossRef] [Google Scholar]
 Trainor, R. F., Steidel, C. C., Strom, A. L., & Rudie, G. C. 2015, ApJ, 809, 89 [NASA ADS] [CrossRef] [Google Scholar]
 Trebitsch, M., Verhamme, A., Blaizot, J., & Rosdahl, J. 2016, A&A, 593, A122 [Google Scholar]
 Vanzella, E., De Barros, S., Cupani, G., et al. 2016, ApJ, 821, L27 [NASA ADS] [CrossRef] [Google Scholar]
 Verhamme, A., Schaerer, D., & Maselli, A. 2006, A&A, 460, 397 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Verhamme, A., Schaerer, D., Atek, H., & Tapken, C. 2008, A&A, 491, 26 [Google Scholar]
 Verhamme, A., Dubois, Y., Blaizot, J., et al. 2012, A&A, 546, A111 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Verhamme, A., Orlitová, I., Schaerer, D., & Hayes, M. 2015, A&A, 578, A7 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Verhamme, A., Orlitová, I., Schaerer, D., et al. 2017, A&A, 597, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Waskom, M., Botvinnik, O., Hobson, P., et al. 2014, seaborn: v0.5.0 (November 2014), https://doi.org/10.5281/zenodo.12710 [Google Scholar]
 Wisotzki, L., Bacon, R., Blaizot, J., et al. 2016, A&A, 587, A98 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Yang, H., Malhotra, S., Gronke, M., et al. 2016, ApJ, 820, 130 [NASA ADS] [CrossRef] [Google Scholar]
 Yang, H., Malhotra, S., Gronke, M., et al. 2017, ApJ, 844, 171 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A: Table of shellmodel parameters
Table A.1 shows the retrieved shellmodel parameters, the MUSE UNIQUE ID (Herenz et al. 2017), and the p(χ^{2}) values of the bestfit. In this table, we use (as in the rest of this work) the notation where A denotes the median and b (c) the difference to the 16th and 84th percentile.
Note that the full version of Table A.1 is available at the CDS and under http://bit.ly/aspectraofMW
Retrieved shellmodel parameters alongside with the MUSE UNIQUE ID and the p(χ^{2}) value of the bestfit.
All Tables
Retrieved shellmodel parameters alongside with the MUSE UNIQUE ID and the p(χ^{2}) value of the bestfit.
All Figures
Fig. 1 Signaltonoise ratio (S/N) of the Lyα line versus the quality of the fit (p(χ^{2}) computed as described in Sect. 3.1). The color coding denotes the “confidence” category as defined in Herenz et al. (2017). 

In the text 
Fig. 2 The four spectra with the lowest p(χ^{2})values of (from left to rightmost panel) { < 0.1, 0.5,6,14 } × 10^{3}. The data is shown in darkblue, and the best fit shellmodel spectrum with a solid red line (the dotted red line is the intrinsic spectrum). In addition, the semitransparent black lines show 10 randomly drawn spectra from the burnedin chains. Apart from the spectrum in the third panel which has been graded with “confidence” level 2 the other three spectra have been classified as “3”. 

In the text 
Fig. 3 Uncertainties on the obtained fitting parameters as a function of the S/N. Specifically, we show half of the 84th minus the 16th percentile of the corresponding posterior distributions. The color corresponds to the “confidence” category assigned to each spectrum by Herenz et al. (2017) with red for the highest confidence (3), followed by green and blue for the lowest confidence class. 

In the text 
Fig. 4 Correlation matrix between the median obtained shellmodel parameters. The color coding denotes the “confidence” category the spectrum was assigned (blue = 1, green = 2, red = 3). While the lower triangle shows scatter plots with each spectrum marked individually, the same data is plotted as contours of kerneldensity estimators in the upper right triangle. The panels along the diagonal show stacked histograms of the (projected) distribution. 

In the text 
Fig. 5 Extent of the Lyα halo and it’s brightness versus the six obtained shellmodel parameters. The upper and central rows show the Kron radius (in kpc) and the line flux (with ) as a function of the fitting parameters (see Sect. 3.3 for more details). The lower row shows the flux ratio between 3 × and 1 × R_{Kron}. The color coding corresponds again to the “confidence” category (blue, green, red in increasing order). Also shown in each panel are the Spearman correlation coefficients. 

In the text 
Fig. 6 Quality of the fits of the “shell” versus the “GaussianminusGaussian” models. The color stands for the three “confidence” categories. Also included (as a “guide to the eye”) are linear fits to the data, and the identity function (as black dashed line). 

In the text 
Fig. 7 Example spectral fits of the “shell” and “GaussianminusGaussian” models with all combinations of fitting failures/successes. The blue lines show the data, the red solid line the shell model fit (with the intrinsic spectrum as dotted red line), and the green dashed line shows the “GaussianminusGaussian” best fit. The numbers in the panel show the p(χ^{2}) values of the best fits in the corresponding colors. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext 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 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.