Open Access
Issue
A&A
Volume 688, August 2024
Article Number A28
Number of page(s) 18
Section Catalogs and data
DOI https://doi.org/10.1051/0004-6361/202348185
Published online 30 July 2024

© The Authors 2024

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.

1 Introduction

Galaxies hosting strong events of star formation can provide keys to our understanding of early galaxy evolution, since such extreme star-forming objects are expected to be more frequent toward the first epochs of the Universe. In this sense, these highly star-forming galaxies provide useful windows to the high-z Universe, allowing us to study the photon budget required for its reionization (e.g. Stark 2016; Dayal & Ferrara 2018). The nature, properties and prevalence of these galaxies at the Epoch of Reionization (EoR) remain to be fully understood.

All these galaxies have been named, as a class, as extreme-emission-line galaxies (EELGs) since their (optical) spectra show very strong emission lines with extreme equivalent width (EW). EELGs host a considerable number of massive, young, star-forming systems producing a substantial amount of pho-toionizing radiation. This radiation is absorbed by the surrounding metal-poor gas (e.g., Ravindranath et al. 2020). Since the gas is expected to have a perforated irregular distribution, hard ionizing radiation could escape (e.g., Bergvall et al. 2006; Izotov et al. 2016; Fletcher et al. 2019; Pérez-Montero et al. 2020). It turns these systems into perfect laboratories for the sutdy of LyC escape. Therefore, the discovery and characterization of local EELG analogs give us a way to gain more insight into the first moments of galaxy evolution when leakage of LyC photons from EELGs provided a major contribution to the EoR photon budget (e.g., Erb et al. 2016; Yang et al. 2017a; Naidu et al. 2022; Matthee et al. 2022).

Extreme-emission-line galaxies present spectra with very high EWs of [O III] and Hβ emission lines, and it has been found that these galaxies also used to be strong Lyα emitters (Tang et al. 2021). On the other hand, most typical LyC leakers presenting a very high EWs (Lyα) are frequently found among Green Pea galaxies, an outstanding family among the prototypes of EELGs (Izotov et al. 2016, 2018; Pérez-Montero et al. 2021).

As of today, observational studies have discovered thousands of EELGs, but most searches have produced a mixed bag of denominations, typically inspired by their compact optical appearance and/or color in plates; examples of these are Green Peas, Blue Berries, ELDots/HαDots, or SDSS extreme [H II] galaxies and their cohorts, to cite a few (e.g., Terlevich et al. 1991; Cardamone et al. 2009; Amorín et al. 2010, 2012; Bekki 2015; Yang et al. 2017b; Salzer et al. 2020), all illustrating how EELG selection can depend on the redshift range explored and the exact criteria and technique applied. The search for EELGs is somewhat complementary to the historical searches for blue galaxies of very compact morphology that were identified and analyzed in pioneering works since the middle of the last century (e.g., Haro 1956; Zwicky 1966; Markarian 1967). However, it is important to emphasize here that many of EELG searches performed so far have been based on broad-band photometry data, thus avoiding direct spectroscopic or narrow-band-based selection based on emission lines, which remain scarce (see, e.g., Indahl et al. 2021; Lumbreras-Calle et al. 2022; Iglesias-Páramo et al. 2022). This situation may have somehow supported the persistence of systematics and/or obvious bias in current samples of EELG, given the typical very faint continua and the strong emission lines of these galaxies1. A clear example are the differences in the number and occurrence of these objects. Paalvast et al. (2018) explored galaxies with extreme [O III]/[O II] ratios, combining data from several MUSE Guaranteed Time Observing (GTO) programs (0.28 < ɀ < 0.85) and founding 15 of these systems in their sample (3.7% of their total sample). Cardamone et al. (2009), imposing conditions on rSDSS, redshift, optical colours, and morphology, reported a spatial density of Green Peas of approximately two per square degree, a critical benchmark in EELG studies. In the local Universe, Yang et al. (2017a) found a density of 0.003 EELGs per square degree for “blueberries” in the local Universe (ɀ ≤ 0.05), while Lumbreras-Calle et al. (2022) identified 466 EELGs within 2000 square degrees, equivalent to a density of about 0.23 per square degree, in the J-PLUS survey (ɀ ≤ 0.06). Furthermore, Amorín et al. (2015) reported 165 EELGs in the 1.7 square degree area covered by the zCoSMOS-bright survey, indicating a density of approximately 97 per square degree in the 0.11 ≤ ɀ ≤ 0.93 redshift range. As Iglesias-Páramo et al. (2022) noted, precise comparisons among different works are prevented by the different selection criteria and observational limits of each sample.

There is not a clear limit for the rest-frame equivalent width (EWo) to define a galaxy as an EELG. In this work, we used a similar definition for EELGs to that used by Iglesias-Páramo et al. (2022), that is, objects that show at least one emission line ([O II]λλ3726, 3729, [O III]λλ4959, 5007 or Hα) with EWo >300 Å. Nevertheless, we are aware that previous works have assumed different selection limits for EELGs (e.g., Amorín et al. 2015; Salzer et al. 2020; Iglesias-Páramo et al. 2022 and references therein); therefore, for the sake of consistency with previous works we enlarged our search producing an extended sample including all sources presenting EWo ≥200 Å. In this way, we gain consistency in order to compare with literature work, and we were able to examine the behavior of the main properties of EELGs around the adopted EWo limit; this an especially relevant exercise when looking for any possible evolution through redshift.

In this work, we took advantage of the exceptional depth and quality of the recent MUSE Hubble Ultra-Deep Field Survey data release (HUDFS; Bacon et al. 2023) to carry out a purely spectroscopic search for EELG with MUSE. EELGs usually show a compact appearance, although some low-surface-brightness structure can be seen for some of them (see Appendix A of Iglesias-Páramo et al. 2022 for examples). Here, we were able to spatially resolve the ionization structure and the gaseous kinematics of this type of object. This paper is organized as follows. In the next section, we summarize the observations and criteria to select the sample. Section 3 shows the basic properties of our EELG candidates. The results and discussion are presented in Sect. 4. The main conclusions are summarized in Sect. 5.

2 Data and sample

2.1 MUSE Hubble Ultra-Deep Field

The MUSE instrument is a powerful tool for optical (4750– 9350 Å) spectroscopic studies thanks to a combination of a wide field of view (FoV) of 1′× 1′ with a spatial sampling of 0.2″ × 0.2″ per spaxel and high resolving power (ranging from 4R4 ~ 1800 in the blue to R ~ 4000 Å in the red). Our study is based in the second data release (DR2) of the MUSE Hubble Ultra-Deep Field Surveys (Bacon et al. 2023), which includes the deepest spectroscopic survey ever performed. With an achieved depth of 141 h and a field of view of 1 arcmin of diameter, it is a key addition to the previous MUSE spectroscopic survey data release (Bacon et al. 2017; Inami et al. 2017; DR1, a 3 × 3 arcmin2 mosaic of nine MUSE fields at a 10-h depth and a single 1 × 1 arcmin2 at 31-h depth) in the Hubble Ultra-Deep Field (Beck-with et al. 2006; HUDF) area. This second data release (DR2) incorporates the previous ones reprocessed with the same tools and methodology, providing a homogeneous data set of deep spectroscopic observations in the HUDF region.

In addition to the calibrated data, the data cubes have an extension with the propagated variance estimate, which is used in the line fitting process (see Sect. 3.1). Reduced and calibrated data cubes, the source catalogs, and an interactive browser are available to the public through the AMUSED web interface2.

2.2 Sample selection

Our selection process begins by analyzing the information offered by the AMUSED interface. The main catalog gives a summary of source properties, and the final data cubes and associated images (e.g., Hubble images) are public for each object. We follow the criteria in Iglesias-Páramo et al. (2022) and define our primary sample of EELG as those objects that show at least one emission line with EWo >300 Å. We use the EWo measured by Bacon et al. (2023) using the python code pyPlatefit3. Briefly, pyPlatefit is an improved Python-based version of the PLATEFIT IDL software (Brinchmann et al. 2004) developed to fit MUSE data. It performs a stellar continuum fitting using a simple population model (Bruzual & Charlot 2003) and then fits the emission lines with Gaussian profiles after subtracting the continuum. To characterize the local continuum, pyPlatefit masks the emission lines and uses windows of 50 Å on the left and right.

Using these criteria, of the 2221 cataloged objects, we found a total of two, five, and three galaxies presenting rest-frame EWo ≥300 Å in [O II], [O III], and Hα, respectively. One of these objects presents EWo ≥300Å in two lines ([O III] and Hα). Thus, we selected nine candidates (ID 91, 891, 895, 1795, 2478, 2532, 6465, 7373, and 7601) fulfilling the criteria for our primary EELG sample. A detailed closer inspection of the spectra of these nine candidates was performed on a one-by-one basis. In doing so, we found that the spectra of two of the candidates (ID 1795 and 7601) show low signal-to-noise ratios (SN) over the entire wavelength range, giving nonrealistic line fitting in the AMUSED database. Therefore, these two candidates were discarded from our EELG primary sample. Table 1 shows the basic properties of the sample; these include quantities extracted from AMUSED (MUSE ID, coordinates, redshift, magnitude in the HST F775W band, stellar mass, star formation rate, line of detection, line flux, corresponding luminosity,4 and rest-frame EWo), and measurements of this work for line-flux and rest-frame EWo. All these relevant spectrophotometric properties are considered in the discussion.

Furthermore, as mentioned before, we also produced an extended EELG sample with sources presenting 200 A ≤ EWo ≤ 300 A in at least one of the lines considered above ([O II]λλ3726, 3729, [O III]λλ4959, 5007 or Hα). Using these extended criteria, we found a total of three, eight, and one galaxies showing 200 Å≤ EWo ≤300 Å in [O II], [O III], and Hα, respectively. One of the objects detected in [O III] presents a EWo ≥300 Å in Hα (ID 7373) and one object shows 200 Å ≤ EWo ≤ 300 Å in both [O III] and Hα (ID 6474). Therefore, this gives ten candidates (ID 1093, 1426, 1561, 1699, 6474, 6865, 1863, 7105, 7985, 8000) for the extended EELG sample. We performed a detailed closer inspection of the spectra of these ten candidates, as was done for the primary sample. Following the same criteria, we end up discarding four galaxies (ID 1863, 7105, 7985 and 8000). The galaxies selected for our extended EELG sample are listed in Table 2.

Figure 1 shows the distribution of redshift, magF775w, EWo[O III]5 and L[O III]. These candidates expand the following ranges: 0.1 < ɀ < 0.9, 22 < magF775w < 28, 200 Å < EWo < 1200 A and 41.9log(L[O III]/ergs)44.2$41.9 \le \log \left( {{L_{[{\rm{O <sc>III</sc>}}]}}/{{{\rm{erg}}} \over {\rm{s}}}} \right) \le 44.2$, respectively. We can see thatboth primary andextendedsamples share the samerange of basic properties. We note that the low number of objects in our sample is totally expected based on reported densities of similar objects in the literature (e.g., Iglesias-Páramo et al. 2022). Furthermore, we do not find systems with a redshift ≤0.1 to be consistent with previous works based in the local Universe (e.g., Yang et al. 2017b; Lumbreras-Calle et al. 2022).

Table 1

Basic properties of EELG primary sample candidates.

Table 2

Same as Table 1, but for EELG extended sample candidates.

thumbnail Fig. 1

Histograms of distribution of properties of EELG sample (see text for details).

3 Properties of the EELG sample

3.1 EELG spectra and line fluxes

Figure 2 shows the spectra of both EELG primary and extended samples. All the selected objects present characteristic spectra of strong-emission-line galaxies, consisting of a faint continuum and very prominent emission lines.

As we indicate in Sect. 2.2, our sample selection is based in the information given by the AMUSED web interface and the fittings performed by Bacon et al. (2023). The equivalent-width values depend on the continuum level, and this can be difficult to characterise. Therefore, we decide to re-evaluate the EWo values of our samples to confirm the EELG classification. Emission-line fitting was performed with pyPlatefit, including improvements with respect the previous fittings. Firstly, we use the parameter ziter, which allows a second iteration improving the fitting. Secondly, we perform 100 Monte Carlo simulations to compute uncertainties following standard procedures in full spectra-fitting techniques (e.g., del Moral-Castro et al. 2019, 2020; Bittner et al. 2019; Domínguez-Gómez et al. 2023, and references therein). Briefly, we derive these uncertainties creating multiple realisations of the original spectrum adding random noise consistent with the quality of the observed spectrum to the best fit and running pyPlatefit again. The considered errors are the standard deviation of the 100 recovered values. Tables C.1 and D.1 present the fluxes measured for the emission lines for the primary and extended samples.

Thirdly, three objects (MUSE ID 91, 7373, and 1093) present a cut in the spectrum (from ~5800 to ~5970Å). Then, to avoid possible problems with the fitting to the continuum, we fit the wavelength range from 6000 to 9300 Å. As for the rest of objects, the whole wavelength range was considered (4750–9300 Å). We notice that the object 2478 presents a very high EWo([O III]λ5007). This value seems to be too high, since the F([O III]λ4959) / ([O III]λ5007) is ~1/3. Therefore, EWo([O III]λ5007) should be ~3·EWo([O III]λ4959). However, it is ~30. Through an inspection of the line catalog given by AMUSED, we found that the continuum associated with the [O III]λ5007 emission line is almost an order of magnitude lower than the continuum level of the [O III]λ4959 emission line. However, Fcont([O III]λ4959)/Fcont([O III]λ5007) ~1 for the rest of the objects. Therefore, we decided to use the continuum level of the [O III]λ4959 emission line to compute the EWo([O III]λ5007) of MUSE ID 2478.

Finally, we studied the effect of the window width used by pyPlatefit to characterize the local continuum6 on the EWo values. We used pyPlatefit to change the internal configuration to use windows of 100 and 200 Å confirming the EELG classification except for the object MUSE ID 7373, where the EWo values for windows of 100 and 200 Å are 102±11 and 106±7, respectively. Therefore, we added a flag of caution for this galaxy in Table 1.

thumbnail Fig. 2

Spectra for both our primary (top panel) and extended (bottom panel) EELG samples. MUSE ID, data set, redshift, and emission line used for the classification are indicated for each spectra.

3.2 Spectroscopic and spatially resolved properties

Figure 3 presents an example of photometric and spectroscopic properties derived for the galaxy ID 891. The upper row shows the HST 775W and MUSE white-light images, as well as the MUSE spectrum with an inset zoomed-in view around the Hγ and [O III]4363 Å emission lines. In the bottom row, we present the spatially resolved spectroscopic and kinematic information of the galaxy, which we discuss in Sect. 4.2. Emission-line fitting was performed using pyPlatefit spaxel by spaxel. The velocity field shows a clear pattern of rotation similar to those of nearby disk galaxies (i.e., Barrera-Ballesteros et al. 2014; del Moral-Castro et al. 2019). Only spaxels with significant line detections (S/N > 6) are shown. We checked that the results are not contingent upon the choice of the minimum S/N by repeating the analysis with different thresholds (3, 6, and 9). We prefer to maintain S/N = 6 to preserve the best possible spatial resolution for this analysis, combining it with a robust selection. Appendices A and B show similar images for the EELG primary and extended samples, respectively.

4 Results and discussion

4.1 Physical properties and chemical abundances

In order to compute the physical properties of the galaxies of our EELG sample,7 we followed the expressions in Pérez-Montero (2017). The Balmer decrement was used for the derivation of the reddening coefficient C(Hβ) assuming the Case B approximation for electron temperature 15 000 K and density 100 cm−3 (Storey & Hummer 1995), and the extinction law by Cardelli et al. (1989) with RV = 3.1. In any case, no relevant change is expected if the exact electron temperature and density of each object are used. Whenever possible, the ratio of the bright Balmer lines Hα/Hβ was used. When the flux measurement of Hα (ID 1426) or Hβ (ID 6465) is uncertain, or Hα is not observed (ID 2478). The C(Hβ) values for all the EELGs are presented in Table 3. After inspection of the errors of the flux of the Balmer lines considered above (see Appendices C and D), a formal 20% relative error (face value) was adopted for C(Hβ).

For two galaxies (ID 2478, ID 6465), the computed C(Hβ) was zero or slightly negative, though consistent to within the errors with the adopted value of C(Hβ) = 0. A C(Hβ) = 0 was adopted also for ID 2532 since only one Balmer line, Hβ, was measured for this galaxy. In what follows, all the fluxes used for the derivation of the physical properties and chemical abundances of the selected EELGs have been reddening-corrected accordingly.

Overall, the extinction values derived for our sample of EELGs are low, typically C(Hβ) ≈ 0.1 (Av ≈ 0.22), and the maximum is around Av ≈ 0.45, whereas a significant fraction of the sample shows values consistent with no extinction. This result, given the negligible Galactic extinction (Schlafly & Finkbeiner 2011), would translate into a low (or very modest) dust component associated with our EELGs, which is somewhat expected for low-metallicity (dwarf) galaxies.

For seven galaxies, the flux in the auroral [O III]λ4363 temperature-sensitive line was measured with a good signal-to-noise level, illustrating the high quality of the data. Therefore, we performed a direct calculation of their electron temperature, with significant implications for the chemical abundance derivation. We followed Pérez-Montero (2017) for the computation of the T[O III] electron temperature using the ([O III]λ5007 +[O III]λ4959)/[O III]λ4363 line ratio. The values of T[O III] obtained, along with their corresponding errors derived via error propagation from the measured flux error, are presented in Table 3. The T[O III] electron temperatures are high, ranging from 1.3 × 104 K to over 20 × 104 K; these temperature values are typically found in low-metallicity star-forming galaxies (e.g., Amorín et al. 2015; Kehrig et al. 2016; Izotov et al. 2019; Pérez-Montero et al. 2021).

The electron density, Ne, was derived using the line ratios from the [O II]λλ3726, 3729 and [S II]λλ6717, 6731 doublets whenever observed (redshift permitting; see Fig. 5 for an example), with both line doublets being resolved in the MUSE data. For galaxies ID 91, 6465, 1093, 1426, 1561, 1699, 6865, we have measurements of the [O II]λλ3726, 3729 doublet, whereas the [S II]λλ6717, 6731 doublet is measured for ID 891, 895, 6474. For three galaxies, ID 2478, 2532, 7373, we have no information to derive the electron density. We assume an electron temperature of 104 K (Sanders et al. 2016; Harshan et al. 2020). After a close inspection of the measured flux errors, in this work we adopted a formal relative error for Ne of 20% (face value) and a value of the low density limit for Ne ≤ 50 cm−3. The electron densities derived for our EELG galaxies are shown in Table 3. These values of electron density indicate that most galaxies of the sample have a low density, and Ne in four of them is consistent with the low-density limit. However, for three galaxies, IDs 1699, 6474, and 6865, substantially high electron densities are derived, with Ne ≈ 500 to 600 cm−3. Interestingly, two of them, ID 6474 and ID 6865, show the highest electron temperatures derived.

After obtaining the physical properties of the galaxies above, we derived their metallicity using the spectroscopic information gathered. We calculated the oxygen chemical abundance of the ionized gas of these EELGs, as a direct measurement of their gas metallicity. To do so, we proceeded as follows.

1. For galaxies (ID 91, 10938, 6865) for which we have measurements of the [O II]λλ3726, 3729 and [O III]λλ5007, 4959 flux, together with the electron temperature T[O III], we can perform a direct derivation of the ionic O+/H+ and O++/H+ abundances, adding them up to derive the total oxygen abundance, expressed as 12+log(O/H). For this calculation, a two-ionization-zone scheme was adopted following Pérez-Montero (2017), using T[O III] for the O++ zone and estimating the O+ zone temperature T[O II]. The 12+log(O/H) of these galaxies, along with their corresponding errors estimated from error propagation of line fluxes and temperature, are quoted in Table 3.

2. For galaxies (ID 891, 895, 2478, 6474) for which we derived the electron temperature, T[O III], but no [O II]λλ3726, 3729 lines are measured, we cannot proceed with the direct method as above. In these cases, we rely on the anticorrelation between the total oxygen abundance of the ionized gas and its electron temperature, as observed in Amorín et al. (2015); this is a direct, theoretically expected consequence of the cooling of the ionized gas via oxygen lines. The values of 12+log(O/H) for these galaxies, along with the corresponding statistical error ±0.08 dex, were derived following the expression in Amorín et al. (2015) and are presented in Table 3. For the sake of comparison, we derived 12+log(O/H) also for the three galaxies treated with the direct method above, giving a very good agreement within the errors.

3. For the rest of the galaxies, we only have flux measurements of the [O II]λλ3726, 3729 and [O III]λλ5007, 4959 lines, for ID 6465, 1426, 1561, 1699, or even only [O III]λλ5007, 4959 for ID 2532, 7373. In this third case, in order to determine the oxygen abundance we have to rely on an empirical calibration of abundance based on bright lines. Among the empirical calibrations, the one based on the R23 parameter, originally defined by Pagel et al. (1979) as R23 = ([O II]λλ3726, 3729 +[O III]λλ5007, 4959)/Hβ, stands out as a useful and robust tool against possible ionization structure effects (e.g., Relaño et al. 2010; Kehrig et al. 2016). In this work, we made use of the R23 calibration (Izotov et al. 2019) to derive the oxygen abundance 12+log(O/H). However, the R23 calibration is bi-valuated, and additional information is required to choose the appropriate branch; thus, we applied it as in Williams et al. (2023), selecting the low-metallicity regime taking into account the high [O III]λλ5007, 4959/[O II]λλ3726, 3729 measured for the EELG galaxies (e.g., Sanders et al. 2016). The derived 12+log(O/H) abundances and their corresponding statistical errors of ±0.1 dex are presented in Table 3. We also show the values of R23 when available, or logO3=[O III]λλ5007, 4959)/Hβ when [O II]λλ3726, 3729 lines are not measured.

It has been reported that the ionized interstellar medium in some extreme-emission-line galaxies could be experiencing density-bounded conditions (e.g., Jaskot & Oey 2013). These objects would show extreme [O III]λλ5007, 4959/[O II]λλ3737, 3729 flux ratios (typically ≳40), suggestive of the presence of nonstandard conditions, which may complicate the applicability of abundance calibrations. The abundance calibration applied in this work is selected from Izotov et al. (2019), and it was specifically derived for the class of extreme emission line low-metallicity galaxies studied here. All sample galaxies for which the empirical calibration has been applied present moderate to low values of the [O III]λλ5007, 4959/[O II]λλ3737, 3729 ratio (see Tables C.1, D.1); hence, any potential effect associated with density bounded conditions is expected to be minor. The analysis of the oxygen abundances derived for our EELG sample indicates that these galaxies are all metal poor, with 7.35 ≤ 12 + log(O/H) ≤ 8.05, with four of them showing a metallicity below 0.1 solar. The very low metallicities obtained in our search for EELGs have significant implications for their characterization as local analogs of the earliest star-forming galaxies.

Table 3

Physical properties and chemical abundances.

thumbnail Fig. 3

MUSE ID 891: (a) HST F775W image with the spatial coverage of the MUSE data cube; (b) MUSE white light image of the galaxy; (c) spectrum of the object with an inset window showing a zoomed-in view around the Hγ and [O III]λ4363 emission lines; and maps of (d) S/N, (e) normalized flux, (f) relative radial velocity, and (g) velocity dispersion (σ), for the emission line indicated in the figure header.

thumbnail Fig. 4

Correlations between luminosity, stellar mass, star formation rate, and σ[O III]. We assume typical uncertainties of 0.3dex for mass and star formation rate (see Bacon et al. 2023).

thumbnail Fig. 5

Electron density sensitive lines. Left panel: fit to [O II]λλ3726, 3729 doublet in spectrum of galaxy ID 91. Right panel: fit to [S II]λλ6717, 6731 doublet in spectrum of galaxy ID 891. In each plot, the observed spectrum (blue), its fit (orange), and the residuals (green) are shown.

4.2 Kinematical properties, structure, and velocity maps

By inspecting the kinematics maps presented in Sect. 3.2, we can separate our sample in three different classes. These are listed below.

1. Resolved with kinematic pattern9: Spatially resolved galaxies that show a clear pattern of rotation (i.e., presence of a disk, with the velocity increasing outward). ~43% of EELGs in the primary sample are in this class. They include the following: (i) iD 891: The HST F775W continuum image reflects a nonnuclear shape, with several high-surface-brightness knots connected by a diffuse low-surface-brightness component. This structure is not shown in the MUSE white light image (given the lower spatial resolution), which presents a head-tail shape (i.e., a main, bright, star-forming clump located at the center of the image and a low-surface-brightness tail). The S/N map shows its highest values in the inner galaxy, though these peak slightly south of the flux maximum, which is likely a consequence of the different line profile shapes and the higher velocity dispersion measured. The ionized-gas-velocity field presents a clear rotating disk structure, showing higher values in the region corresponding to the knots identified in the HST image. We note that these knots could be reminiscent of the structures seen in tidal dwarf and/or cometary galaxies (e.g., Lagos et al. 2016; Roche et al. 2023). These galaxies are usually considered the result of a tidal process or interactions of disk galaxies. No apparent companion is envisaged in this figure; however, further investigation on its true nature is beyond the scope of this paper; (ii) ID 895: The HST image shows an elongated shape, which is also reflected by the MUSE white light image. The [O III] S/N map shows a nuclear shape, while the Hα map also presents an elongated shape, with the highest values peaking slightly south of the flux maximum. Both [O III] and Hα flux maps present a nuclear shape peaking at the center of the galaxy. The [O III] and Hα velocity field show clear rotation disk structures; (iii) ID 91: The HST F775W image presents a nuclear shape. This shape is also shown in the MUSE white-light, S/N, and flux maps. The [O III] velocity field shows a rotating pattern.

2. Resolved without kinematic pattern: These are spatially resolved galaxies, but, unlike the previous category, they do not show a rotation pattern. This class represents ~14% and ~83% of our EELGs (ID: 2478) and extended samples, respectively (ID: 1093, 1426, 1561, 6474, and 6865). The HST F775W images reflect a nuclear shape, except for the ID 6474 object, which presents an elongated shape with a low-surface-brightness tail. The nuclear shape is also presented in the MUSE white light, S/N, and flux maps for the six galaxies.

3. Not resolved: These are galaxies whose structure cannot be resolved. ~43% and ~17% of both, EELGs primary (ID: 2532, 6465 and 7373) and extended (ID: 1699) samples, respectively, are included in this class.

We note that we applied a similar approach to Guérou et al. (2017), considering a galaxy spatially resolved if the galaxy has an area of at least 16 MUSE spaxels (i.e., about 1.5 times the PSF FWHM size of the data cube). We are aware that mergers may mimic the appearance of rotation patterns in observations of finite spatial resolution. Therefore, observations with better spatial resolution (e.g., JWST or MUSE-NFM observations) are necessary to reliably determine the internal structure of these systems.

4.3 Global properties and evolution

Figure 4 displays the relations among the derived [O III] luminosity and velocity dispersion, stellar mass and SFR (both parameters taken from the AMUSED database, as derived using the Prospector code10), for both the primary and extended samples. inspection of this figure first reveals that, when comparing the behavior of both samples, no clear differences can be found between them, i.e. the EWo limits assumed to select our EELGs do not seem to introduce a bias in their observed properties. Overall, all four parameters show good correlations, with the SFR versus stellar mass plot delineating a very clear main sequence for our sample of EELGs. From the ionized gas, the derived [O III] line luminosity and σ correlate, both parameters showing some degree of correlation versus stellar mass and SFR, most notably between stellar mass and [O III] luminosity. This is relevant since we are comparing properties, e.g. [O III] luminosity and stellar mass, which were independently derived, corresponding to gaseous and stellar components of the EELG galaxies, respectively.

In Fig. 6, we present the mass-metallicity relation showing average loci corresponding to nearby (SDSS) and high-redshift galaxies. Our EELG galaxies11, sampling the low-stellar-mass (masspro ≤9.0) and low-metallicity (12+log(O/H)≤8.0) regions of the plot, are consistent with the mass-metallicity relation derived for ɀ~2.2 (Sanders et al. 2016) galaxies. They are also close to the locus defined by very-high-redshift objects (ɀ = 7.7– 8–5 and ɀ = 6.1–6.3), as recently studied by Schaerer et al. (2022) and Sun et al. (2023) using JWST data. This result illustrates how the EELGs in our sample can be considered true local analogs of these high-redshift systems. The next generation of ground-based 30-m telescopes, together with additional observations from space, will allow us to fill this figure with low metallicity and low massive EELG galaxies at high-redshift and their local analogs.

thumbnail Fig. 6

Mass-metallicity relation. Oxygen chemical abundance and stellar mass for galaxies in the primary (blue circles) and extended (orange circles) samples in this work, compared to mean relations observed at ɀ~0 (Andrews & Martini 2013), ɀ<0.22 (Duarte Puertas et al. 2022), and ɀ = 2.2 (Sanders et al. 2016). Black stars and red squares show recent results for high-redshift galaxies from JWST data (Schaerer et al. 2022; Sun et al. 2023).

5 Summary

For decades, EELGs searches have been performed on broadband photometry data. Recently, thanks to the development of new instruments and recent surveys such as SDSS or J-PAS, direct spectroscopy or narrow-band searches have been carried out (e.g., Indahl et al. 2021; Iglesias-Páramo et al. 2022). However, the number of EELG searches and analyses performed to date based on IFU data is (to our knowledge) very scarce. In fact, this work could be considered among the first EELG searches performed based on deep IFU data surveys, since previous works have mostly studied specific objects (e.g., Lofthouse et al. 2017; Bosch et al. 2019; Arroyo-Polonio et al. 2023).

Taking advantage of the deepest IFU survey to date, we present spatially resolved kinematic maps for a sample of intermediate redshift EELGs observed from the ground. Furthermore, we confirm that these systems are true local analogs of high-redshift systems (ɀ = 6–8) studied by the JWST. While these cases already reveal the fantastic potential of the MUSE Ultra Deep Field data set, our small sample prevents a more detailed study of the basic properties of EELGs. However, this work can be regarded as a successful pilot study, demonstrating the power of the presented methodology to analyze strong emission-line galaxies using IFU data. Future deep observations with MUSE or ongoing instruments such as HARMONI will result in a larger sample to which to apply the methodology presented here.

Acknowledgements

We thank the referee for carefully reading our manuscript and providing constructive comments improving the paper. This study uses data provided by the Muse Hubble Ultra Deep survey, publicly available at the AMUSED database (https://amused.univ-lyon1.fr/). I.M.C. acknowledges funding from U. La Laguna through the Margarita Salas Fellowship from the Spanish Ministry of Universities ref. UNI/551/2021-May 26 and under the EU Next Generation funds, and from ANID programme FONDE-CYT Postdoctorado 3230653. J.M.V., J.I.P., A.A.P. acknowledge financial support from grants CEX2021-001131-S and PID2019-107408GB-C44 funded by MCIN/AEI/ 10.13039/501100011033. This paper makes use of python (http://www.python.org); Matplotlib (Hunter 2007): a suite of open-source python modules that provide a framework for creating scientific plots; Astropy (http://www.astropy.org): a community-developed core Python package and an ecosystem of tools and resources for astronomy (Astropy Collaboration et al. 2013, 2018, 2022); MPDAF (https://mpdaf.readthedocs.io): an open-source Python package that provides tools to work with MUSE-specific data; and pyplatefit (https://github.com/musevlt/pyplatefit) (Bacon et al. 2023): a python routine for emission and absorption lines fit of MUSE spectra. and AMUSED (https://amused.univ-lyon1.fr/): a public web interface for inspection and retrieval of extragalactic MUSE data products. IMC thank Alejandro S. Borlaff for helpful discussions on the HST broad-band images.

Appendix A Spectra and images of the primary EELG sample candidates

In this appendix, we show both the spectra and the images of the selected EELG candidates.

thumbnail Fig. A.1

Same as Fig. 3, but for MUSE ID galaxy 895.

thumbnail Fig. A.2

Same as Fig. 3, but for MUSE ID galaxy 91.

thumbnail Fig. A.3

Same as Fig. 3, but for MUSE ID galaxy 2478.

thumbnail Fig. A.4

Same as Fig. 3, but for MUSE ID galaxy 2532.

thumbnail Fig. A.5

Same as Fig. 3, but for MUSE ID galaxy 6465.

thumbnail Fig. A.6

Same as Fig. 3, but for MUSE ID galaxy 7373.

Appendix B Spectra and images of the EELG extended sample

Here, we present both the spectra and the images of the extended sample.

thumbnail Fig. B.1

Same as Fig. 3, but for MUSE ID galaxy 1093.

thumbnail Fig. B.2

Same as Fig. 3, but for MUSE ID galaxy 1426.

thumbnail Fig. B.3

Same as Fig. 3, but for MUSE ID galaxy 1561.

thumbnail Fig. B.4

Same as Fig. 3, but for MUSE ID galaxy 1699.

thumbnail Fig. B.5

Same as Fig. 3, but for MUSE ID galaxy 6474.

thumbnail Fig. B.6

Same as Fig. 3, but for MUSE ID galaxy 6865.

Appendix C Line fluxes of the EELG primary sample candidates

In this appendix, we show the fluxes and equivalent width measured with pyPlatefit for our EELG sample. Only lines with an S/N higher than three are included in the tables.

Table C.1

Observed line fluxes for the primary sample. (1) Emission line; (2) rest-frame wavelength; (3) observed line fluxes for ID 91; (4) observed line fluxes for ID 891; (5) observed line fluxes for ID 895; (6) observed line fluxes for ID 2478; (7) observed line fluxes for ID 2532; (8) observed line fluxes for ID 6465; (9) observed line fluxes for ID 7373. a Uncertain value.

Appendix D Line fluxes of the EELG extended sample candidates

Here, we present the fluxes and equivalent width measured with pyPlatefit for our extended sample. Only lines with an S/N higher than three are included in the tables.

Table D.1

Same as Table C.1, but for the extended sample. a Uncertain value.

References

  1. Amorín, R. O., Pérez-Montero, E., & Vflchez, J. M. 2010, ApJ, 715, L128 [CrossRef] [Google Scholar]
  2. Amorín, R., Vflchez, J. M., Hägele, G. F., et al. 2012, ApJ, 754, L22 [CrossRef] [Google Scholar]
  3. Amorín, R., Pérez-Montero, E., Contini, T., et al. 2015, A&A, 578, A105 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  4. Andrews, B. H., & Martini, P. 2013, ApJ, 765, 140 [NASA ADS] [CrossRef] [Google Scholar]
  5. Arroyo-Polonio, A., Iglesias-Páramo, J., Kehrig, C., et al. 2023, A&A, 677, A114 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Astropy Collaboration (Robitaille, T. P., et al.) 2013, A&A, 558, A33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Astropy Collaboration (Price-Whelan, A. M., et al.) 2018, AJ, 156, 123 [Google Scholar]
  8. Astropy Collaboration (Price-Whelan, A. M., et al.) 2022, ApJ, 935, 167 [NASA ADS] [CrossRef] [Google Scholar]
  9. Bacon, R., Conseil, S., Mary, D., et al. 2017, A&A, 608, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Bacon, R., Brinchmann, J., Conseil, S., et al. 2023, A&A, 670, A4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Barrera-Ballesteros, J. K., Falcón-Barroso, J., García-Lorenzo, B., et al. 2014, A&A, 568, A70 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Beckwith, S. V. W., Stiavelli, M., Koekemoer, A. M., et al. 2006, AJ, 132, 1729 [Google Scholar]
  13. Bekki, K. 2015, MNRAS, 454, L41 [NASA ADS] [CrossRef] [Google Scholar]
  14. Bergvall, N., Zackrisson, E., Andersson, B. G., et al. 2006, A&A, 448, 513 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Bittner, A., Falcón-Barroso, J., Nedelchev, B., et al. 2019, A&A, 628, A117 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Bosch, G., Hägele, G. F., Amorín, R., et al. 2019, MNRAS, 489, 1787 [NASA ADS] [CrossRef] [Google Scholar]
  17. Brinchmann, J., Charlot, S., White, S. D. M., et al. 2004, MNRAS, 351, 1151 [Google Scholar]
  18. Bruzual, G., & Charlot, S. 2003, MNRAS, 344, 1000 [NASA ADS] [CrossRef] [Google Scholar]
  19. Cardamone, C., Schawinski, K., Sarzi, M., et al. 2009, MNRAS, 399, 1191 [NASA ADS] [CrossRef] [Google Scholar]
  20. Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 345, 245 [Google Scholar]
  21. Dayal, P., & Ferrara, A. 2018, Phys. Rep., 780, 1 [Google Scholar]
  22. del Moral-Castro, I., García-Lorenzo, B., Ramos Almeida, C., et al. 2019, MNRAS, 485, 3794 [NASA ADS] [CrossRef] [Google Scholar]
  23. del Moral-Castro, I., García-Lorenzo, B., Ramos Almeida, C., et al. 2020, A&A, 639, A9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  24. Domínguez-Gómez, J., Pérez, I., Ruiz-Lara, T., et al. 2023, Nature, 619, 269 [CrossRef] [Google Scholar]
  25. Duarte Puertas, S., Vilchez, J. M., Iglesias-Páramo, J., et al. 2022, A&A, 666, A186 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  26. Erb, D. K., Pettini, M., Steidel, C. C., et al. 2016, ApJ, 830, 52 [NASA ADS] [CrossRef] [Google Scholar]
  27. Fletcher, T. J., Tang, M., Robertson, B. E., et al. 2019, ApJ, 878, 87 [Google Scholar]
  28. Guérou, A., Krajnovic, D., Epinat, B., et al. 2017, A&A, 608, A5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  29. Haro, G. 1956, Bol. Observatorios Tonantzintla Tacubaya, 2, 8 [NASA ADS] [Google Scholar]
  30. Harshan, A., Gupta, A., Tran, K.-V., et al. 2020, ApJ, 892, 77 [NASA ADS] [CrossRef] [Google Scholar]
  31. Hunter, J. D. 2007, Comput. Sci. Eng., 9, 90 [NASA ADS] [CrossRef] [Google Scholar]
  32. Iglesias-Páramo, J., Arroyo, A., Kehrig, C., et al. 2022, A&A, 665, A95 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  33. Inami, H., Bacon, R., Brinchmann, J., et al. 2017, A&A, 608, A2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. Indahl, B., Zeimann, G., Hill, G. J., et al. 2021, ApJ, 916, 11 [NASA ADS] [CrossRef] [Google Scholar]
  35. Izotov, Y. I., Schaerer, D., Thuan, T. X., et al. 2016, MNRAS, 461, 3683 [Google Scholar]
  36. Izotov, Y. I., Worseck, G., Schaerer, D., et al. 2018, MNRAS, 478, 4851 [Google Scholar]
  37. Izotov, Y. I., Guseva, N. G., Fricke, K. J., & Henkel, C. 2019, A&A, 623, A40 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  38. Jaskot, A. E., & Oey, M. S. 2013, ApJ, 766, 91 [Google Scholar]
  39. Kehrig, C., Vflchez, J. M., Pérez-Montero, E., et al. 2016, MNRAS, 459, 2992 [NASA ADS] [CrossRef] [Google Scholar]
  40. Lagos, P., Demarco, R., Papaderos, P., et al. 2016, MNRAS, 456, 1549 [NASA ADS] [CrossRef] [Google Scholar]
  41. Lofthouse, E. K., Houghton, R. C. W., & Kaviraj, S. 2017, MNRAS, 471, 2311 [NASA ADS] [CrossRef] [Google Scholar]
  42. Lumbreras-Calle, A., López-Sanjuan, C., Sobral, D., et al. 2022, A&A, 668, A60 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. Markarian, B. E. 1967, Astrofizika, 3, 24 [NASA ADS] [Google Scholar]
  44. Matthee, J., Naidu, R. P., Pezzulli, G., et al. 2022, MNRAS, 512, 5960 [NASA ADS] [CrossRef] [Google Scholar]
  45. Naidu, R. P., Matthee, J., Oesch, P. A., et al. 2022, MNRAS, 510, 4582 [CrossRef] [Google Scholar]
  46. Paalvast, M., Verhamme, A., Straka, L. A., et al. 2018, A&A, 618, A40 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  47. Pagel, B. E. J., Edmunds, M. G., Blackwell, D. E., Chun, M. S., & Smith, G. 1979, MNRAS, 189, 95 [NASA ADS] [CrossRef] [Google Scholar]
  48. Pérez-Montero, E. 2017, PASP, 129, 043001 [CrossRef] [Google Scholar]
  49. Pérez-Montero, E., Kehrig, C., Vflchez, J. M., et al. 2020, A&A, 643, A80 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  50. Pérez-Montero, E., Amorín, R., Sánchez Almeida, J., et al. 2021, MNRAS, 504, 1237 [CrossRef] [Google Scholar]
  51. Ravindranath, S., Monroe, T., Jaskot, A., Ferguson, H. C., & Tumlinson, J. 2020, ApJ, 896, 170 [NASA ADS] [CrossRef] [Google Scholar]
  52. Relaño, M., Monreal-Ibero, A., Vílchez, J. M., & Kennicutt, R. C. 2010, MNRAS, 402, 1635 [CrossRef] [Google Scholar]
  53. Roche, N., Vílchez, J. M., Iglesias-Páramo, J., et al. 2023, MNRAS, 523, 270 [NASA ADS] [CrossRef] [Google Scholar]
  54. Salzer, J. J., Gronwall, C., Lipovetsky, V. A., et al. 2000, AJ, 120, 80 [NASA ADS] [CrossRef] [Google Scholar]
  55. Salzer, J. J., Feddersen, J. R., Derloshon, K., et al. 2020, AJ, 160, 242 [NASA ADS] [CrossRef] [Google Scholar]
  56. Sanders, R. L., Shapley, A. E., Kriek, M., et al. 2016, ApJ, 816, 23 [Google Scholar]
  57. Schaerer, D., Marques-Chaves, R., Barrufet, L., et al. 2022, A&A, 665, L4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  58. Schlafly, E. F., & Finkbeiner, D. P. 2011, ApJ, 737, 103 [Google Scholar]
  59. Stark, D. P. 2016, ARA&A, 54, 761 [Google Scholar]
  60. Storey, P. J., & Hummer, D. G. 1995, MNRAS, 272, 41 [NASA ADS] [CrossRef] [Google Scholar]
  61. Sun, F., Egami, E., Pirzkal, N., et al. 2023, ApJ, 953, 53 [NASA ADS] [CrossRef] [Google Scholar]
  62. Tang, M., Stark, D. P., Chevallard, J., et al. 2021, MNRAS, 503, 4105 [NASA ADS] [CrossRef] [Google Scholar]
  63. Terlevich, R., Melnick, J., Masegosa, J., Moles, M., & Copetti, M. V. F. 1991, A&AS, 91, 285 [NASA ADS] [Google Scholar]
  64. Williams, H., Kelly, P. L., Chen, W., et al. 2023, Science, 380, 416 [NASA ADS] [CrossRef] [Google Scholar]
  65. Yang, H., Malhotra, S., Gronke, M., et al. 2017a, ApJ, 844, 171 [NASA ADS] [CrossRef] [Google Scholar]
  66. Yang, H., Malhotra, S., Rhoads, J. E., & Wang, J. 2017b, ApJ, 847, 38 [NASA ADS] [CrossRef] [Google Scholar]
  67. Zwicky, F. 1966, ApJ, 143, 192 [NASA ADS] [CrossRef] [Google Scholar]

1

Pioneering spectroscopic surveys were carried out using wide field slitless spectroscopy techniques (see Salzer et al. 2000 and references there).

4

We note that this luminosity is not corrected by the extinction of each galaxy. See Sect. 4.1.

5

We note that two galaxies (MUSE iD: 891 and 7373) are included in the primary sample performing its condition in EWoHα, though their EWo[O III] are lower than 300 A.

6

Although, the broad-band Hubble Space Telescope (HST) observations are deeper than those of MUSE, the filter widths (F606W: 4632.15–7179.43 Å and F775W: 6801.41–8630.74 Å; see Fig. 1 of Bacon et al. 2023) are too large to ignore the presence of nonuniform spectral-energy-distribution (SED) components and other emission features in the characterization of the local continuum.

7

A closer inspection of all individual line profiles confirmed the presence of narrow lines in all the sources, discarding AGN contamination. Furthermore, the BPT diagram was created for three of them, locating the galaxies in the star formation region.

8

A second fitting in the rest-frame 4000–5800 Å was perfomed to fit the [O II]λλ3726, 3729 doublet.

9

We note that this class presents significantly larger velocity amplitudes.

10

All the details of the derivation of stellar mass and SFR can be found in Bacon et al. (2023).

11

Only eight of the galaxies in Table 3 have stellar mass values in the AMUSED database.

All Tables

Table 1

Basic properties of EELG primary sample candidates.

Table 2

Same as Table 1, but for EELG extended sample candidates.

Table 3

Physical properties and chemical abundances.

Table C.1

Observed line fluxes for the primary sample. (1) Emission line; (2) rest-frame wavelength; (3) observed line fluxes for ID 91; (4) observed line fluxes for ID 891; (5) observed line fluxes for ID 895; (6) observed line fluxes for ID 2478; (7) observed line fluxes for ID 2532; (8) observed line fluxes for ID 6465; (9) observed line fluxes for ID 7373. a Uncertain value.

Table D.1

Same as Table C.1, but for the extended sample. a Uncertain value.

All Figures

thumbnail Fig. 1

Histograms of distribution of properties of EELG sample (see text for details).

In the text
thumbnail Fig. 2

Spectra for both our primary (top panel) and extended (bottom panel) EELG samples. MUSE ID, data set, redshift, and emission line used for the classification are indicated for each spectra.

In the text
thumbnail Fig. 3

MUSE ID 891: (a) HST F775W image with the spatial coverage of the MUSE data cube; (b) MUSE white light image of the galaxy; (c) spectrum of the object with an inset window showing a zoomed-in view around the Hγ and [O III]λ4363 emission lines; and maps of (d) S/N, (e) normalized flux, (f) relative radial velocity, and (g) velocity dispersion (σ), for the emission line indicated in the figure header.

In the text
thumbnail Fig. 4

Correlations between luminosity, stellar mass, star formation rate, and σ[O III]. We assume typical uncertainties of 0.3dex for mass and star formation rate (see Bacon et al. 2023).

In the text
thumbnail Fig. 5

Electron density sensitive lines. Left panel: fit to [O II]λλ3726, 3729 doublet in spectrum of galaxy ID 91. Right panel: fit to [S II]λλ6717, 6731 doublet in spectrum of galaxy ID 891. In each plot, the observed spectrum (blue), its fit (orange), and the residuals (green) are shown.

In the text
thumbnail Fig. 6

Mass-metallicity relation. Oxygen chemical abundance and stellar mass for galaxies in the primary (blue circles) and extended (orange circles) samples in this work, compared to mean relations observed at ɀ~0 (Andrews & Martini 2013), ɀ<0.22 (Duarte Puertas et al. 2022), and ɀ = 2.2 (Sanders et al. 2016). Black stars and red squares show recent results for high-redshift galaxies from JWST data (Schaerer et al. 2022; Sun et al. 2023).

In the text
thumbnail Fig. A.1

Same as Fig. 3, but for MUSE ID galaxy 895.

In the text
thumbnail Fig. A.2

Same as Fig. 3, but for MUSE ID galaxy 91.

In the text
thumbnail Fig. A.3

Same as Fig. 3, but for MUSE ID galaxy 2478.

In the text
thumbnail Fig. A.4

Same as Fig. 3, but for MUSE ID galaxy 2532.

In the text
thumbnail Fig. A.5

Same as Fig. 3, but for MUSE ID galaxy 6465.

In the text
thumbnail Fig. A.6

Same as Fig. 3, but for MUSE ID galaxy 7373.

In the text
thumbnail Fig. B.1

Same as Fig. 3, but for MUSE ID galaxy 1093.

In the text
thumbnail Fig. B.2

Same as Fig. 3, but for MUSE ID galaxy 1426.

In the text
thumbnail Fig. B.3

Same as Fig. 3, but for MUSE ID galaxy 1561.

In the text
thumbnail Fig. B.4

Same as Fig. 3, but for MUSE ID galaxy 1699.

In the text
thumbnail Fig. B.5

Same as Fig. 3, but for MUSE ID galaxy 6474.

In the text
thumbnail Fig. B.6

Same as Fig. 3, but for MUSE ID galaxy 6865.

In the text

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

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

Initial download of the metrics may take a while.