M Paalvast
email: paalvast@strw.leidenuniv.nl

A Verhamme

L A Straka

J Brinchmann

E C Herenz

D Carton

M L P Gunawardhana

L A Boogaard

S Cantalupo

T Contini

B Epinat

H Inami

R A Marino

M V Maseda

L Michel-Dansac

S Muzahid

T Nanayakkara

G Pezzulli

J Richard

J Schaye

M C Segers

T Urrutia

M Wendt

L Wisotzki

Properties and redshift evolution of star-forming galaxies with high

Keywords: Galaxies: evolution, ISM, abundances -ISM: structure, kinematics, and dynamics -Reionisation

published or not. The documents may come

Introduction

Extreme oxygen line ratios (O32 ≡ [O iii]λ5007/[O ii]λ3727> 4) were proposed recently as a potential tracer of the escape of ionising radiation from galaxies through density-bounded H ii regions (Jaskot & Oey 2013;Nakajima & Ouchi 2014). The idea is the following: if a galaxy is leaking ionising photons through a density-bounded region, the ratio of O32 can be high if the H ii region that we observe is truncated, for example if (part of) the [O ii] region is missing. We see deeper into the ionised region and the external layer of [O ii] is either nonexistent or thinner than in the classical ionisation-bounded scenario. Given their high O32 ratios, Jaskot & Oey (2013) discuss the possibility of LyC escape from "Green Pea" (GP) galaxies, a population of extremely compact, strongly star-forming galaxies in the local Universe (Cardamone et al. 2009;Izotov et al. 2011;Yang et al. 2016). Naka-jima & Ouchi (2014) and Nakajima et al. (2016) compare O32 ratios of different types of high-redshift galaxies, Lyman Break Galaxies (LBGs), and Lyman Alpha Emitters (LAEs) with GPs and Sloan Digital Sky Survey (SDSS) galaxies: high-redshift galaxies have on average higher O32 ratios than SDSS galaxies, but comparable to GPs. Furthermore, the observed O32 ratios of LAEs are larger than those of LBGs. Along the same line, GPs are also strong LAEs (Henry et al. 2015;Yang et al. 2016;Verhamme et al. 2017;Yang et al. 2017), which is very unusual for galaxies in the local Universe (Hayes et al. 2011;Wold et al. 2014).

While the O32 ratio of galaxies that are leaking ionising photons may be enhanced compared to those with a LyC escape fraction, f LyC esc , equal to zero, there are other situations that can lead to high O32 ratios. For example, the O32 ratio depends on metal-Article number, page 1 of 17 arXiv:1808.04899v1 [astro-ph.GA] 14 Aug 2018 licity: low stellar and nebular metallicities lead to higher O32 ratios (Jaskot & Oey 2013). A harder ionising spectrum will also induce higher O32 ratios, as investigated in for example Pellegrini et al. (2012), as well as a higher ionisation parameter (e.g. Stasińska et al. 2015). Furthermore, shocks could also explain these ratios, as studied in detail in Stasińska et al. (2015).

Despite intensive searches for LyC emission from galaxies, only a few LyC leakers have been identified over the last decades in the local Universe (Bergvall et al. 2006;Leitet et al. 2013;Borthakur et al. 2014;Leitherer et al. 2016), but most searches resulted in non-detections or upper limits (Siana et al. 2015;Mostardi et al. 2015;Grazian et al. 2016;Rutkowski et al. 2016Rutkowski et al. , 2017)). The discovery of the link between LyC emission and O32, however, turned the tide, as demonstrated by for example Izotov et al. (2016aIzotov et al. ( ,b, 2018a,b),b). For their studies, LyC emission was detected for all eleven galaxies at z ≈ 0.3 that were selected by their extreme O32 ratios (O32 > 4), among other criteria such as brightness, compactness, and strong Hβ equivalent widths. Furthermore, a correlation between O32 and the escape of ionising photons was found, although the scatter of f LyC esc is large (Izotov et al. 2018b). At high redshift (z ≈ 3), four galaxies with high escape fractions (> 50%) have been reported (Vanzella et al. 2015;de Barros et al. 2016;Shapley et al. 2016;Bian et al. 2017;Vanzella et al. 2018), which were selected by similar criteria. Additionally, the recent results from the Lyman Continuum Escape Survey (Fletcher et al. 2018) reveal an average escape fraction of ∼ 20% for galaxies at z ≈ 3 with strong [O iii] emission, and a weak correlation between f LyC esc and the [O iii] equivalent width for ∼ 20 galaxies with directly detected LyC emission. Although the combination of these selection criteria has resulted in relatively few galaxies with confirmed LyC emission yet, the detection of extreme O32 emission from a local low-mass GP analogue (Micheva et al. 2017) might, however, suggest that low-mass extreme O32 emitters, and thus possible low-mass LyC emitters, are more common than the bright GP samples suggest. A statistical study of the O32 ratios of emission-line selected galaxies over a broad range of stellar masses has, however, not been performed so far.

The unique capabilities of the Multi-Unit Spectroscopic Explorer (MUSE) [START_REF] Bacon | Proc. SPIE[END_REF]) allow us to study the properties of galaxies with extreme O32 ratios and how common they are in emission-line selected samples. For this study we combine four MUSE Guaranteed Time Observing (GTO) surveys and collect a sample of mainly emission-line detected galaxies with a high specific star formation rate and stellar masses between ∼ 10 6 and ∼ 10 10 , from which we compute the distribution of O32 ratios in a blind survey of star-forming galaxies. We will here present the properties and occurrences of extreme oxygen emitters spanning the redshift range 0.28 < z < 0.85, where both lines are in the MUSE spectral range, in the largest statistical sample of emission-line detected galaxies in three-dimensional spectral data. This article is organised as follows: in Sect. 2 we describe the data from different programmes that we used for this study; in Sect. 2.6 we describe the sample selection; in Sect. 3 we investigate the occurrence of high O32 ratios and study potential correlations with stellar mass ( M * ), star formation rate (SFR), and the metallicity indicator R23 line ratio; in Sect. 4 we study the incidence rate of galaxies with high O32 ratios as a function of M * and z, and we also discuss how our results compare to nebular models with no escape of ionising photons. We end with a discussion on the most extreme oxygen emitters. Throughout this paper we adopt a cosmology with H 0 = 70 km s -1 Mpc -1 , Ω m = 0.3 and Ω Λ = 0.7.

Data

Different MUSE surveys

For this study we used data taken with MUSE, as part of GTO observations, covering a wavelength range of 4800-9300 Å. We selected data from four surveys that together span an area of more than 55 arcmin 2 . Below follows a description of the different surveys.

Hubble Ultra Deep Field survey

The MUSE Hubble Ultra Deep (HUDF) survey (Bacon et al. 2017) consists of two fields of different size and depth in the original HUDF region. The Medium Deep Mosaic Field (e.g. UDF-mosaic) is a mosaic of nine pointings (≈3 x 3 ) with a depth of approximately ten hours. The UDF-10 Ultra Deep Field (e.g. UDF-10) is a deeper observation of a single pointing within the UDF-mosaic region, with a depth of approximately 31 hours, covering ≈1.15 arcmin 2 . The data reduction is described by Bacon et al. (2017) and is based on the MUSE standard pipeline version 1.7dev (Weilbacher et al. 2014). For the construction of the redshift catalogue (Inami et al. 2017), Hubble Space Telescope (HST) priors are used, as well as a blind search for emission lines in the datacube using the software ORIGIN (Mary et al. in prep).

MUSE-Wide survey

The MUSE-Wide survey aims to cover a larger field of view than the HUDF survey with relatively short exposures of one hour per pointing. The data release of the full survey will be presented in Urrutia et al (in prep). Here we focus on the first 24 pointings of MUSE-Wide, which together cover an area of 22.2 arcmin 2 . We use the source catalogue that is presented in Herenz et al. (2017), which is created using the emission-line detection software LSDCat (Herenz & Wisotzki 2017), but we do not use their supplied spectra since we extract the spectra of all surveys consistently.

MUSE QuBES survey

The MUSE Quasar Blind Emitter Survey (QuBES) consists of 24 individual fields centred on quasars. All datacubes have a minimum total exposure time of two hours, with selected fields observed to a depth of ten hours based on the availability of high-quality archival auxiliary data such as quasi-stellar object (QSO) spectra. The standard data reduction was carried out with the MUSE Data Reduction System (DRS; Weilbacher et al. 2014). Post-processing procedures for additional integral field unit (IFU) normalisation and sky subtraction were carried out using the CubEx package (Cantalupo et al. in prep;see Marino et al. 2018 andBorisova et al. 2016 for details). The survey will be fully described in Straka et al. (in prep) and Segers et al. (in prep). A subset of 21 fields is used for this study based on the availability of galaxy catalogues. For this study, the presence of the QSO in the field is not important. The galaxy catalogues are compiled as follows. First, white light images are created from the MUSE datacubes by summing the entire cube along the wavelength axis. Then, SExtractor (Bertin & Arnouts 1996) Paalvast et al.: O32 ratios of galaxies in MUSE data is used to detect any objects down to 1σ. The spectra for each object are extracted by selecting associated pixels in the segmentation maps produced by SExtractor. The spectra of the resulting detections are then inspected by eye in order to determine the galaxy redshifts based on nebular emission lines and stellar absorption features.

Galaxy groups survey

The last dataset added to our sample is that of the Galaxy group survey (Epinat et al. in prep), which targets galaxy groups at intermediate redshift (z ≈ 0.5-0.7) from the zCOSMOS 20k group catalogue (Knobel et al. 2012). We selected data from three galaxy groups, namely COSMOS-Gr30, 34, and 84, with the deepest MUSE data of 9.75, 5.25, and 5.25 hours respectively, ~1 x 1 each. The data reduction followed the same approach as the Hubble Ultra Deep Field and is described in Epinat et al. (2018) for the galaxy group field COSMOS-Gr30. For the construction of the redshift catalogues, galaxies were selected from the COSMOS photometric catalogue by Laigle et al. (2016), complemented by emission-line detection using ORIGIN for the deepest field COSMOS-Gr30 (see also Epinat et al. 2018).

Spectrum extraction and emission-line flux measurements

The spectra of all sources are extracted from the datacubes using the same method in order to make the line flux measurements comparable. We followed the approach of Inami et al. (2017) and extracted the spectra using a mask region, which is the HST segmentation map convolved by the MUSE point spread function for all surveys except the MUSE QuBES survey, for which there is no HST coverage. Spectra in this survey are extracted using a mask region that is constructed from the MUSE white light image. We then used the simple unweighted sum of the flux in the mask region and used this as the spectrum of each galaxy. We measured the line fluxes of the galaxies in the catalogues using the software platefit (Tremonti et al. 2004;Brinchmann et al. 2008). Since this method is the same as that used by Inami et al. (2017) to construct the HUDF emission-line catalogue, the line flux measurements that we use here are identical to theirs.

Deriving SFRs and dust extinction

The SFRs of our galaxies are calculated using the method described in Brinchmann et al. (2013). In short, we simultaneously fit the Charlot & Longhetti (2001) models to the brightest (signal-to-noise ratio (S/N)>3) optical emission lines. From this we estimate the SFR marginalising over: metallicity, ionisation parameter, dust-to-metal ratio, and the optical depth of the dust attenuation (τ V ). The advantage of using a multi-emission-line approach over using a single Balmer line to calculate the SFR is that it is less affected by sky line contamination of a single line and should therefore provide a more robust SFR. Also, this method provides an estimate of τ V , which we adopt to correct the emission-line fluxes for dust extinction.

Calculating stellar masses

We obtained stellar masses for the galaxies by performing spectral energy distribution (SED) fitting using the fast (Fitting and Assessment of Synthetic Templates) algorithm (Kriek et al. 2009), where we used the Bruzual & Charlot (2003) library and assumed exponentially declining star formation histories (SFHs) with a Calzetti et al. (2000) extinction law and a Chabrier (2003) initial mass function (IMF). To test the influence of our SFH assumption, we compare our stellar masses with those derived by the magphys code (da Cunha et al. 2008). For these M * estimations, the photometry is fitted to stellar population synthesis models assuming random bursts of star formation in addition to an exponentially declining SFH. We find consistent stellar masses, for example the median difference between M * derived from fast and M * from magphys equals log M * /M = 0.08 with a standard deviation of log M * /M = 0.3, from which we conclude that the influence of the assumed SFH is small for this sample. For the UDF, we used HST Advanced Camera for Surveys (ACS) and Wide Field Camera 3 (WFC3) photometry from the catalogue of Rafelski et al. (2015). The same approach is applied to the MUSE-Wide (photometry from Skelton et al. 2014) and the Galaxy groups (photometry from Laigle et al. 2016). Unfortunately there is no deep photometry available for the MUSE-Qubes survey. We therefore used a set of 11 400 Å -wide boxcar filters, where emission lines are masked, to compute a pseudophotometric SED, as will be described in more detail in Segers et al. (in prep). Since not all the bright emission lines lie within the MUSE spectral range for our redshifts, we decided to leave the photometry uncorrected for bright emission lines. To test the possible effect of strong emission lines on M * , we compared our stellar masses to those based on the photometry that is corrected for emission lines in the MUSE spectral range. We find that for the bulk of the galaxies in our sample, this effect is negligible; for example, the maximum difference between the emission-line corrected and the non-corrected M * for the extreme O32 emitters, which will be introduced in Sect. 3, corresponds to log M * /M = 0.03.

Redshift distribution

For our main sample, we only select galaxies with spectra of sufficient quality for our study, that is, spectra with S/N > 3 in the [O iii]λ5007 line at 0.28 < z < 0.85, since for this redshift interval the [O ii] and [O iii] emission lines fall within the MUSE wavelength range. For galaxies that meet this criterion, but have S/N([O ii]) < 3, we use the 3-σ lower limit for the O32 ratio. We show two examples of spectra that meet these criteria in Fig. 1. This leads to a total sample of 815 galaxies, of which the redshift distribution for each survey is shown in Fig. 2. This histogram shows the number of galaxies in the redshift bins, separated based on the survey from which the data originates. Because the depths of the different surveys that we combine here are not uniform, care is required for selecting galaxies for our final sample.

Sample selection

Observations of star-forming galaxies have shown that the star formation rate (SFR) and stellar mass (M * ) of these galaxies are tightly related (e.g. Brinchmann et al. 2004;Noeske et al. 2007). This relation, also called the star-forming main sequence (SFMS), has been studied intensively in the last decade (e.g. Whitaker et al. 2012Whitaker et al. , 2014;;Lee et al. 2015;Santini et al. 2017). Boogaard et al. (2018) have constrained the SFMS for galaxies that are detected in deep MUSE data, modelling them as a Gaussian distribution around a three-dimensional (3-D) plane, taking into account and obtaining the redshift evolution of the SFR - Redshift distribution of 815 galaxies in total, in the range 0.28 < z < 0.85, with high confidence, spectroscopically determined MUSE redshift, in the four surveys that we are using for this project.

M * relation, resulting in log SFR = 0.83 +0.07 -0.06 log

M * M 0 -0.83 +0.05 -0.05 + 1.74 +0.66 -0.68 log 1 + z 1 + z 0 ± 0.44 +0.05 -0.04 , (1)
with M 0 = 10 8.5 M * and z 0 = 0.55. Because this relation is derived for a sample of galaxies with deep photometry and high S/N Balmer lines from MUSE spectra of galaxies with stellar masses down to log M * /M ≈ 7, and comparable to the mass range of the galaxies in our sample, we use this z-dependent M * -SFR relation to select galaxies for our final sample. We calculate how much the SFR of a galaxy is offset from the redshiftcorrected SFMS, which we will herein refer to as the 'distance to the SFMS' (∆ SFMS) given in dex with a sign such that objects above the SFMS have a positive distance. The distribution of the distances to the SFMS is shown in Fig. 3 for each survey that we use for this study separately. The dashed line represents the SFMS, with, at the left side, galaxies below and, at the right side, galaxies above the SFMS. For all four surveys, this distribution peaks within 1-σ from the SFMS. Since our sample is mostly emission-line selected, we expect our sample to be most complete at high SFRs. However, below the SFMS it is likely that a fraction of the galaxies are below the detection threshold. We therefore conservatively select galaxies with a distance > 0 dex from the SFMS, resulting in a final sample of 406 galaxies. Applying such a selection based on a fixed distance from the z-dependent SFMS relation also ensures that we select the same fraction of star-forming galaxies at each redshift.

For an overview of the stellar mass and the SFR distribution of the four surveys, we refer the reader to Figs. A.1 and A.2 in the Appendix.

Based on the stellar mass distribution, we assume that all the galaxies in the sample are pure star-forming systems, since studies show that the fraction of active galactic nucleus (AGN) host galaxies is low at these masses. For example almost all AGN hosts at 0 < z < 1 have log M * /M 10.2 (Vitale et al. 2013) and the fraction of AGNs over all galaxies with log M * /M ≈ 10 at z 0.3 is ∼1 % (Kauffmann et al. 2003;Haggard et al. 2010).

Results

In this section we explore whether there is a correlation between galactic properties and O32. We also study the location of these emitters with respect to the SFMS, since Chisholm et al. (2017) report that LyC leakers lie further above the SFMS.

Redshift distribution of extreme O32 emitters

In Fig. 4 we show the redshift distribution of our sample after the selection that is described in the previous section. In light grey we show the number density of all the galaxies, while in dark 1) in dex of the four different surveys used for this study. We show the 1-σ variation around the mean µ (grey area) calculated from the fitted normal distribution (black line). We included all galaxies with a distance > 0 dex (dashed line), so all galaxies that lie above the main sequence. The number in the upper right corner of the diagram shows the number of galaxies above this threshold. In our final sample we include 406 galaxies in total. In dark grey we show the redshift distribution of the galaxies with O32 > 1 and in light grey the total sample.

grey we show that of galaxies with O32 > 1. For galaxies with S/N([O ii]) < 3, we use the 3-σ upper limits on [O ii], resulting in 3-σ lower limits on O32. We determined the ratio between the number of galaxies with O32 > 1 to the total number of galaxies in a bin. A χ 2 statistical test against the null hypothesis that this fraction in each bin is independent of redshift gives χ 2 = 1.9 a P-value of ∼0.99, which indicates that the fraction of galaxies with oxygen ratios above unity does not evolve as a function of redshift. We adopted the same approach for the extreme oxygen emitters with O32 > 4, as shown in red, which yielded a comparable result ( χ 2 = 4.4, P-value ≈ 0.88). This latter result is not particularly robust, because we only have 15 extreme emitters in our sample.

O32 emitters on the star-formation main sequence

The redshift-corrected distribution in the log SFR -log M * plane is shown in Fig. 5. We normalised our results to z = 0, and show the SFMS of Eq. 1 corrected to z = 0 (dashed line). The dashed line shows this relation for z = 0, above which we selected galaxies in our sample, as described in Sect. 2.6. The points are coloured by the logarithm of the dust-corrected O32 ratio. The symbols show galaxies in the final selection (circles), galaxies below the selection threshold (triangles), and galaxies with a lower limit on the O32 ratio (squares).

gen ratios larger than unity (blue points) are overall more abundant above the SFMS (104 out of 406 galaxies, 26%) than below this relation (46 out of 324 galaxies, 14%). The number of extreme O32 emitters deviates even more; 15 versus one galaxies in O32 regime above and below the SFMS, respectively. Regarding only galaxies in the final sample, there is no clear correlation between distance from the SFMS and O32 ratio. Moreover, we also find that galaxies with O32 > 1 are more common at low masses (logM * /M < 9). A two-dimensional plot of O32 versus the distance to the SFMS is shown in Fig. 6. We will now turn to a more quantitative discussion of the relation between the oxygen ratio and SFR and M * .

O32 as a function of stellar mass

In Fig. 7 we plot the log O32 line ratio versus the stellar mass (black dots). The red squares denote the median values in both the x-and y-directions, in the intervals of log O32 between -1.0 and 1.0 with a step size of 0.5. Additionally we show 3-σ lower limits of the O32 ratio for galaxies without an The median values in the O32 bins show a clear anticorrelation between the oxygen ratio and the stellar mass. However, the Spearman's rank correlation coefficient for individual galaxies equals -0.30 (P-value ≈ 0), indicating no clear trend between log O32 and log M * . The stellar masses of galaxies with O32 > 4 are lower than the average stellar mass of the entire population, for example all the extreme oxygen emitters in our sample have stellar masses below 10 9 . Although the O32 ratios of the extreme emitters in our sample are similar to those of the confirmed LyC leakers, their stellar masses are smaller than those of most of the leaking galaxies from Izotov et al. (2016aIzotov et al. ( ,b, 2018a)), de Barros et al. (2016) and Vanzella et al. (2016). Because their galaxies were, besides their extreme O32 ratios, selected by their brightness to increase the possibility to directly detect LyC photons at z ≈ 0.3, their mass is not necessarily a reflection of the typical mass of an LyC emitter. For example, Izotov et al. (2017) show that galaxies at z < 0.1, which are only selected by extreme O32 ratios, all have masses between 10 6 -10 7 M . In addition, Mrk 71, a near green pea analogue and a LyC emitter candidate, has a stellar mass around 10 5 (Micheva et al. 2017), suggesting that the mass of the bulk of LyC emitters might be lower than what is derived from confirmed LyC leakers. The approximately 20 recently discovered galaxies with LyC emission at z ≈ 3 (Fletcher et al. 2018) show an anti-correlation between the LyC escape fraction and the stellar mass, similar to our results.

O32 as a function of SFR

In Fig. 8 we show the oxygen ratio as a function of SFR with the same colour code as used in Fig. 7. The median values of the SFR decrease with increasing O32. For individual galaxies there is again no clear correlation (Spearman's rank correlation coefficient ≈ -0.35, P-value ≈ 0). The SFR of the confirmed LyC emitters, visualised by the green stars and square, is about two orders of magnitude larger than the median of our galaxies with comparable O32 emission.

When comparing the O32 ratio with the specific star formation rate (sSFR = SFR/M * ), we find that these values are also not correlated and that this also holds for the median values in the O32 bins. The sSFR of the confirmed LyC emitters is on average one order of magnitude larger than the sSFR of our sample.

O32 as a function of R23

R23 is a diagnostic to estimate the gas-phase oxygen abundances (later referred to as metallicity, Z) of galaxies and is based on the ratio of [O ii] + [O iii] over Hβ (Alloin et al. 1979;[START_REF] Osterbrock | Astrophysics of gaseous nebulae and active galactic nuclei[END_REF]Pagel et al. , 1980;;McGaugh 1991), given by (McGaugh 1991

) R23 = [O ii]λ3727 + [O iii]λ4959 + [O iii]λ5007 Hβ . (2)
Since it only relies on the blue rest-frame spectrum, this diagnostic is often used when the red emission lines are out of the spectrum. However, the relation between R23 and Z is degenerate and therefore additional lines are still necessary to constrain the metallicity. In Fig. 9 we plot the logarithm of the oxygen line ratio against the logarithm of R23. At high O32 (log O32 -0.2), we visually determine a trend between O32 and R23, which is followed by the LyC leakers of Izotov et al. (2016aIzotov et al. ( ,b, 2018a)). However, below log O32 ≈ -0.2 the data points scatter in the log R23 direction, as a result of uncertain Hβ measurements as we will return to in the discussion. There we will also discuss how the stellar and gas-phase metallicity influences the oxygen ratio.

Discussion

In the previous section we gave an overview of the galaxy properties of our sample. Here we aim to determine what processes and properties are responsible for a high oxygen line ratio and under which circumstances it is likely to find extreme O32 ratios.

The occurrence of high O32

Although in previous studies 'extreme' oxygen emitters are defined as having O32 > 4, we use O32 > 1 as the threshold for 'high' O32 emitters to study their occurrence, which will be justified later in this section. This results in a significant sample of 104 high O32 emitters, in contrast to applying the O32 > 4 threshold, which would only leave 15 galaxies in the extreme regime. In order to study how our selection criterion influences our results and to study the redshift dependence of our results, we constructed a comparison sample from SDSS data.

Creating a comparison sample from SDSS data

We created a comparison sample of star-forming galaxies in SDSS DR7, from which the derivation of the measurements is detailed in Brinchmann et al. (2004) and Tremonti et al. (2004).

For each galaxy in the MUSE sample we selected a local analogue in SDSS that lies at the same position in the redshiftcorrected SFR -M * plane, thus both samples contain the same number of galaxies. However, because the spectra and consequently the emission-line flux measurement of SDSS galaxies are, unlike the MUSE galaxies, biased by aperture effects due to a finite fiber size, we first corrected for this. Assuming that the slope of the SFMS from Eq. 1 is unaffected, we refitted the relation to the SDSS data, resulting in a shift of +0.2 in the log SFR-2.93 log(1+z) direction. Each SDSS galaxy in the same redshift-corrected SFR -M * position is a potential analogue of a MUSE galaxy. We therefore selected all SDSS galaxies within the 1 -σ error bars of the position of the MUSE galaxy. We used the median value and the 16% and 84% percentiles of the [O iii] and [O ii] emission-line fluxes of all selected SDSS galaxies galaxy as the fluxes and errors of the analogue.

Figure 10 shows the distribution of O32 of our sample (left) and the SDSS comparison sample (middle). We compared the O32 of each galaxy in our sample with its counterpart in the SDSS sample. The result of this is shown in the right panel of Fig. 10. At ∆ log O32 = 0 (black dashed line), the oxygen ratio of the galaxy in our sample is the same as its SDSS analogue. The median value of the ∆ log O32 is at 0.13 dex (blue solid line), which means that the O32 ratio of the MUSE galaxies on average exceeds the O32 of galaxies in the comparison sample. The dashed blue lines, however, indicate that this difference is within the 1-σ error bars.

Incidence rate of high O32 as a function of mass

We divided the sample into three subsets based on stellar mass, with log M * /M in the ranges [7.0, 8.0], [8.0, 9.0], and [9.0, 10.0]. We then derived the fraction of galaxies with O32 > 1 in each mass bin, shown as the black line in Fig. 11, where the points indicate the centres of the mass bins. One-σ errors are derived by bootstrapping the sample 10,000 times (grey area). We then applied the same approach to the SDSS comparison sample (see the red points, line, and shaded area).

The fraction of galaxies in the MUSE sample with O32 > 1 decreases with increasing M * (black line); we find ∼30 % for galaxies with stellar masses between 10 7 and 10 9 M , but ∼10 % in the highest mass bin. This trend is comparable to that of the median bins between O32 and M * , which we described in Sect.

3.3 and Fig. 7. The SDSS comparison sample follows a comparable trend, but the fractions are offset by 0.05-0.2 towards lower fractions. Below log M * /M ≈ 8 the SDSS sample is, however, incomplete (see Appendix B), so we caution that the SDSS results in the lowest mass bin are likely to be biased as a result.

Adopting a metallicity-dependent threshold for the O32 ratio

High oxygen ratios can also be driven by low metallicity systems (see for example the nebular models of Gutkin et al. 2016) due to the harder ionising spectrum of low metallicity stars and less efficient cooling. Here we perform the same analysis as in the previous section, but instead of the fixed threshold at O32 = 1, we adopt a metallicity dependent threshold on O32 that we derive as follows. For each galaxy we derive the metallicity Z using the redshift dependent M * -Z relation of Zahid et al. (2014).

We then set the threshold for this galaxy equal to the O32 ratio from photo-ionisation models from Gutkin et al. (2016) with this metallicity and the ionisation parameter set to log U = -3. We show the relation between the metallicity dependent O32 threshold and stellar mass in Fig. 12, for the minimum and maximum redshifts of our sample (z = 0.28 and z = 0.85). Our assumption that galaxies are highly ionised above log U = -3 results in an O32 threshold for O32 between 0.5 and 1.7. This is the reason for setting the fixed threshold to O32 > 1 in the previous section.

The incidence rate of galaxies with O32 above this metallicitydependent threshold in each mass bin is shown in Fig. 13 for the MUSE sample (black/grey) and the comparison sample (red). We see that with the Z-dependent threshold there is no longer a strong trend between M * and the fraction of high O32 emitters. For the bin with most galaxies (8 < log M * /M < 9, see Table 1), the SDSS and MUSE results are in agreement. This indicates that the relation between M * and the incidence rate of high O32 emitters is most likely the result of the relation between metallicity and oxygen ratio.

Evolution of the fraction of galaxies with high O32

We split the sample into three redshift bins of equal number, with median redshifts of z = 0.42, z = 0.63, and z = 0.74. We then calculated the fractions of galaxies above the metallicitydependent threshold and show the result in Fig. 14. For comparison we also show the incidence rates of the entire SDSS comparison sample, which has a median redshift of z = 0.03.

In the lowest mass bin there seems to be a weak trend between the incidence rate and the redshift, for example the fraction in the highest z bin is significantly higher than the fraction of the SDSS comparison sample. However, the lowest redshift subsample in this mass bin is larger (44) than the intermediate and high-redshift bin (26 and 28 respectively; see also Table 1), which may indicate that we only include the most extreme starforming systems in the high-redshift sample and this can explain the results that we observe in the lowest mass bin. In the two highest mass bins, there is no significant difference between the fraction of O32 > 1 at different redshifts. In Fig. 15 the results for the galaxies with stellar masses between log M * /M = 8 and log M * /M = 9 are presented against the look-back time and redshift, which we calculated using the median redshifts of the redshift bins.

The incidence rate of high oxygen emitters when controlled for metallicity is thus independent of redshift for our sample. Since we selected the galaxies based on their distance to the redshift-dependent main SFMS (Eq. 1), we selected the same fraction of star-forming galaxies at each redshift, rather than selecting the same kind (same SFR and M * and thus sSFR) of galaxies over cosmic time. The incidence rates that we here derived may therefore be extrapolated to the entire population of star-forming galaxies, and suggest that the fraction of high O32 emitters are constant over time between z = 0.28 and z = 0.85.

In the current paradigm, the O32 ratio of high-redshift galaxies are believed to be more extreme than those of low-redshift galaxies. Results from the MOSFIRE Deep Evolution Field (MOSDEF) and the Keck Baryonic Structure Survey (KBSS)-MOSFIRE surveys of galaxies at z ∼ 2.3 show that their O32 ratios are offset towards significantly larger values in comparison to those of local galaxies (Steidel et al. 2014;Sanders et al. 2016;Strom et al. 2017). This may be interpreted by a harder stellar radiation field at fixed mass at higher redshift (Steidel et al. 2014;Strom et al. 2017), or as the result of lower metallicities of high-redshift galaxies at fixed mass (Sanders et al. 2016). These results are supported by cosmological simulations of massive galaxies (Hirschmann et al. 2017) that show that the ionisation parameter and the [O iii]/Hβ ratio increases with redshift at fixed M * at 0 < z < 4. Our results at 0.28 < z < 0.85 , however, support another scenario where the O32 ratio is constant over cosmic time. This difference may be the result of the different M * regime that is probed in the high-redshift surveys (9 log M * /M 11) and in the cosmological simulations (9.5 log M * /M 11.5) with respect to the stellar mass of the galaxies in this work (7 log M * /M 10). Our data sample consists of a combination of several surveys with depths between one and 30 hours. The catalogues for most fields are a mixture of emission lines and continuum-detected galaxies, resulting in samples of different completeness limits in M * and emission-line flux. In Appendix B we study how such a possible incompleteness of our data sample alters our results by simulating different completenesses in flux and M * of SDSS data, and find that the effect is negligible.

We studied the impact of our SFR calibration method on the results and compared them with Hβ-derived SFRs that are dereddened by the Hβ/Hγ fraction. We also re-analysed the data by adopting different stellar libraries for M * calculations and derived similar results. However, for one of the surveys that is used for this study, the MUSE QuBES, we used the MUSE spectrum for the SED fitting instead of deep photometry as we used for the data of the other surveys. We are aware that this induces uncertainty on the mass estimates and therefore re-analysed the results in Sect. 4.1.3 without the data of the MUSE QuBES survey and acquired comparable results.

Can nebular models with no escape of ionising photons

predict the observed O32?

In Sect. 3.5 we discussed the behaviour of our galaxies in the log O32 versus R23 diagram (see also Fig. 9). Here we compare these results with nebular models to study if they are consistent with each other and whether we can derive the nebular metallicity of our galaxies with the R23 method (see Fig. 16). We added the grids of line ratios from nebular models that are calculated by Gutkin et al. (2016) by the coloured squares, where each colour represents a model of a fixed metallicity Z (see the figure legend; metallicity here refers to the combination of nebular and stellar oxygen abundances, since these are kept constant for these models). We connected the grids of models with constant metallicity by dashed lines, where the ionisation parameter U increases towards the upper right from log U = -4 to log U = -1. For the calculation of these models the ionising photon escape fraction was assumed to be zero. The lower limits on O32, which are also upper limits on R23, are not shown here since in this plot because it is difficult to compare these galaxies with models.

Deriving the nebular metallicity of our galaxies by comparing it to the model results is not straightforward, due to the degeneracy of the R23-metallicity relation. However, from Fig. 16 it is clear that the R23 ratio indicates that the metallicity of the majority of the galaxies in our sample is sub-solar and around 0.006 (≈ 1/3 Z ). Hence the R23 of many of the galaxies exceeds the maximum R23 predicted by the models, which can partly be explained by an uncertain Hβ flux measurement (galaxies with an Hβ signal to noise lower than three are shown by open squares).

We have differentiated between models with ionisation parameters between log U = -4 and log U = -2 (solid line) and those of higher values, since observations show that the bulk of star-forming galaxies have log U < -2 (e.g. Shirazi et al. 2014). The ionisation parameter of galaxies with extreme O32 ratios might however exceed those of normal galaxies, as pointed out by Stasińska et al. (2015). However, if the O32 ratio of our galaxies exceeds the predicted ratio of models with log U > -2, the escape of LyC photons is also a likely scenario. For this reason the galaxy with the highest O32 ratio is a promising LyC escape candidate. Although the O32 ratios of the other extreme emitters in our sample (with O32 > 4) in this diagram are similar to those of the confirmed LyC leakers (green stars), the logaritm of R23 of our galaxies scatters within 0.4 of the value of the LyC leakers. Most of them however imply either low stellar and nebular metallicities, high ionisation parameters, the escape of ionising photons, or a combination of these factors. However, comparing the data that lie inside the model grid to these nebular models with no escape of ionising photons is thus not sufficient to determine if LyC escape is responsible for the extreme oxygen emission.

Extreme O32 emitters

In the previous sections we discussed how the properties of galaxies with extreme oxygen ratios are related to the entire sample. Here we study the robustness of the O32 measurement and investigate the electron temperatures of the [O iii] regime in extreme galaxies.

To confirm that the O32 ratios are well measured, we compare them to the [Ne iii]λ3869/[O ii] ratio. This ratio is an alternative diagnostic of the ionisation parameter because of its tight relation to the O32 ratio (Levesque & Richardson 2014). Even though [Ne iii]λ3869 is more than an order of magnitude fainter than [O iii]λ5007, the [Ne iii]/[O ii] ratio is less affected by reddening than the O32 ratio. The relationship between [O iii]λ5007/[O ii]λ3727 and [Ne iii]λ3869/[O ii]λ3727 for the extreme emitters with a significant [Ne iii]λ3869 detection (σ > 3) is shown in Fig. 17. The solid grey line shows the predictions of the Starburst99/Mappings III photo-ionisation models of Levesque & Richardson (2014) with Z = 0.2 Z . We offset the models by a factor of +0.6 in the log [O iii]λ5007/[O ii]λ3727 direction to take into account the discrepancy between these models and the observations of Izotov et al. (2006) and Jaskot & Oey (2013), as discussed by Levesque & Richardson (2014). We use these results to conclude that the [Ne iii]λ3869/[O ii] ratios are consistent with extreme O32 ratios. There are, however, offsets between our observations and the corrected Levesque & Richardson (2014) relation of up to 0.2 dex, which are most likely caused by the lower significance of the [Ne iii]λ3869 line and/or by an offset in the dust attenuation estimate. We note however, that if we use these offsets to correct the [O iii] line fluxes, the O32 ratios will still be in the regime of the extreme oxygen emitters.

For most of the extreme emitters, we have high S/N measurements (S/N > 3) of the auroral 0.0005 0.001 0.002 0.006 0.01 0.014 0.017 0.02 Fig. 16. Logarithm of O32 versus the logarithm of R23. The black points are galaxies from our sample, open symbols reflect galaxies with S/N (Hβ) < 3), while the green stars are the confirmed LyC leakers from Izotov et al. (2016aIzotov et al. ( ,b, 2018a)). The squares that are connected by solid (log U < -2) and dashed coloured lines show the position in this diagram of the nebular model of Gutkin et al. (2016) for constant gas-phase metallicity (see legend) and with increasing log U towards the upper right, with -4 < log U < -1 and step size 0.5. gion (see [START_REF] Osterbrock | Astrophysics of gaseous nebulae and active galactic nuclei[END_REF], using the nebular.ionic routine in Pyneb (Luridiana et al. 2015), assuming an electron density of n e = 100 cm -3 . The [O iii]λ4363/[O iii]λ5007 ratios and the corresponding values of T e are shown in Fig. 18. A complete discussion on this is outside the scope of the paper. However, there is presumably a positive correlation between O32 and T e since the electron temperatures of the galaxies with the most extreme O32 ratios exceed those of the other extreme galaxies. Furthermore, comparing the [O iii]λ4363/[O iii]λ5007 and O32 ratios to the predictions from the Gutkin et al. (2016) nebular models, the metallicity of the galaxies can be constrained, varying from Z≈ 0.0001 (Z/Z ≈0.005) for the galaxy with the highest O32 ratio (diamond) to Z≈ 0.002 -0.004 (Z/Z ≈0.15) for the galaxies with O32≈ 4. Again, for the calculation of these nebular models, the fraction of escaping ionising photons is assumed to equal zero. Although such a scenario would boost the O32 ratio, there is no obvious reason to expect an effect on the [O iii]λ4363/[O iii]λ5007 ratio, and thus the real metallicity of a leaking system would be even lower than what is predicted by the models. Other scenarios, such as a top-heavy initial mass function, might increase the hardness of the ionising spectrum and could, however, also affect the position of a galaxy in this plot.

UDF object 6865

In Fig. 19 indicates that the metallicity of this galaxy is extremely low (Z≈ 0.0001), which deviates somewhat from what is predicted from the R23 ratio (Z≈ 0.0005 -0.001), and the ionisation parameter is high (log U ≈ -1.5). However, all models assume that 2016) models of different metallicity as indicated by the colours, with increasing log U towards the upper right, and log U < -2 (solid line) and log U < -2 (dashed line), and step size of log U = 0.5 between the squares. The colours of the circles are used as in Fig. 17.

there is no escape of ionising photons, which may affect the O32 and R23 ratios. Due to the combination of the relatively low stellar mass, log(M * /M )= 8.16 +0.16 -0.06 , and the redshift, we are not able to precisely derive the size of the object, because the apparent size is comparable to the PSF of the HST image and the source is not resolved in the MUSE data. This, however, indicates that the galaxy is compact, like the objects of Izotov et al. (2016a,b, ). Together with the comparison of its oxygen ratio to those of nebular models, this suggests that this galaxy may be an LyC emitter candidate.

Conclusions

We constructed a sample of emission-line galaxies in the redshift range 0.28 < z < 0.85 that are detected in data from four MUSE GTO surveys. The galaxies are selected based on their position in the SFR -M * plane in a way that we only included galaxies that are above the redshift-dependent SFMS from Boogaard et al. (2018). In this regime we expect the sample to be independent of selection effects. Our final sample consists of 406 galaxies, of which 104 (26%) have a high O32 ratio (O32 > 1) and 15 galaxies are extreme emitters with O32> 4 (3.7%). We studied the O32 ratio as a function of the position in the (redshift-corrected) SFR versus M * diagram, as a function of stellar mass M * , SFR, and metallicity indicator R23. We then studied the incidence rate of galaxies with high oxygen ratios, which is defined by either a fixed threshold, O32 > 1, or by a metallicity-dependent threshold as a function of M * and redshift. The main conclusions of this study are:

-Galaxies with a high oxygen ratio are more common at lower masses (M * < 9) and above the SFMS. There is no clear correlation between distance from the SFMS and the O32 ratio for galaxies in our final sample that are above the SFMS (Fig. 5).

-We find no correlation between O32 ratio and M * , although the median values in O32 bins seems to be anti-correlated (Fig. 7).

-We observe the same trend between the median values of O32 and SFR, but again no significant correlation for individual galaxies. The SFR of most of our extreme emitters is two to three orders of magnitude smaller than those of confirmed leakers (Fig. 8).

-The fraction of galaxies with high O32 ratios is independent of stellar mass when we use a metallicity-dependent O32 threshold (Fig. 13).

-We find no significant correlation between the fraction of high O32 emitters and redshift, suggesting that there is no redshift evolution of the number of high O32 in the redshift range 0.28 < z < 0.85 (Figs. 14 and 15).

-Comparing O32 and R23 of our galaxies with those of nebular models with no escape of ionising photons, we find that some of the high oxygen emitters can be reproduced Paalvast et al.: O32 ratios of galaxies in MUSE data by models with a high ionisation parameter (log U ≈ -2), a very low stellar and nebular metallicity (smaller than ∼ 1/3 Z ), or a combination of both. However, our extreme emitters are in the same regime as the confirmed leakers from Izotov et al. (2016aIzotov et al. ( ,b, 2018a) ) and we therefore cannot exclude the escape of ionising photons from these galaxies. The O32 ratio of our most extreme oxygen emitter can only be explained by models with very high ionisation parameter (log U > -2), from which we conclude that this galaxy may be a LyC leaker candidate (Fig. 16).

-For galaxies with a significant [O iii]λ4363 detection, we derived the [O iii]λ4363/[O iii]λ5007 ratio and the electron temperature and find that these values are similar to or larger than those predicted by nebular models with extremely low metallicity, high ionisation parameters, and constant SFR at t = 3 × 10 8 years. From this we conclude that a part of the extreme O32 emitters may have light-weighted ages of t < 3 × 10 8 years (Fig. 18).

Fig. 1 .

 1 Fig. 1. Example spectra of O32 emitters in our sample. The upper panel shows the spectrum of the galaxy with the most extreme oxygen ratio (dust-corrected O32 = 23) from the UDF-mosaic catalogue with id = 6865 and z = 0.83. The orange shaded lines show the Hγ, Hβ, and [O iii]λλ4959, 5007 of neighbouring sources at z ≈ 0.62. An example spectrum of a galaxy with a lower O32 ratio, in this case with dust-corrected O32 = 0.25, is shown in the lower panel (UDF-mosaic, id = 892, z = 0.74).

Fig. 3 .

 3 Fig. 3. Distribution of the distance to the redshift-dependent SFMS from Boogaard et al. (2018) (Equation1) in dex of the four different surveys used for this study. We show the 1-σ variation around the mean µ (grey area) calculated from the fitted normal distribution (black line). We included all galaxies with a distance > 0 dex (dashed line), so all galaxies that lie above the main sequence. The number in the upper right corner of the diagram shows the number of galaxies above this threshold. In our final sample we include 406 galaxies in total.

Fig. 4 .

 4 Fig. 4. Redshift frequency of extreme O32 emitters with O32 > 4 (red).In dark grey we show the redshift distribution of the galaxies with O32 > 1 and in light grey the total sample.

Fig. 5 .

 5 Fig.5. SFR -M * diagram where we have corrected the SFRs to z = 0 using the redshift evolution from Eq. 1. The dashed line shows this relation for z = 0, above which we selected galaxies in our sample, as described in Sect. 2.6. The points are coloured by the logarithm of the dust-corrected O32 ratio. The symbols show galaxies in the final selection (circles), galaxies below the selection threshold (triangles), and galaxies with a lower limit on the O32 ratio (squares).

Fig. 6 .

 6 Fig.6. O32 ratio versus the distance to the SFMS. Only galaxies above ∆ SFMS = 0 (dashed line) are included in the final sample. Lower limits on O32 are shown by orange triangles.

Fig. 7 .

 7 Fig.7. Stellar mass as a function of log O32 (black dots) and lower limits on O32 (orange triangles). The median per bin with 1-σ errors is shown by the red squares and error bars. We defined the bins by intervals of logO32 between -1.0 and 1.0 with a step size of 0.5. Three-σ lower limits on the O32 ratio for sources with a detection of [O ii] below our S/N cut are shown by the orange triangles. The green stars indicate the positions of the confirmed Lyman continuum leakers fromIzotov et al. (2016aIzotov et al. ( ,b, 2018a) ) and the green square is a LyC emitter at z=3.2 (deBarros et al. 2016;Vanzella et al. 2016). Galaxies above the grey dashed line have extreme O32 ratios (O32 > 4).

Fig. 8 .Fig. 9 .

 89 Fig.8. O32 ratio versus SFR. The colours are described in the caption of Fig.7. All SFRs that are shown here are derived from emission lines, as described in Sect. 2.6.

Fig. 10 .Fig. 11 .

 1011 Fig.10. Number of galaxies per log O32 bin for our sample of 406 galaxies (left panel). The dashed lines correspond to O32 = 1 (left) and O32 = 4 (right). In the middle panel we show a similar plot for our SDSS comparison sample 1 (see the text for details). The distribution of the difference in O32 of each galaxy in our sample with its counterpart in the SDSS comparison sample is shown in the right panel. The median and the 1-σ spread of the distribution are shown by the solid and dashed blue lines respectively.

Fig. 12 .

 12 Fig.12. Metallicity-dependent O32 threshold as a function of stellar mass, for the minimum and maximum redshift of our sample, derived using the redshift-dependent M * -Z relation ofZahid et al. (2014) and the O32 ratio from photo-ionisation models fromGutkin et al. (2016) (see the text for details).

Fig. 13 .

 13 Fig. 13. Fraction of galaxies with O32 greater than the metallicitydependent threshold for our MUSE sample (black/grey) and the SDSS comparison sample (red) (see the text for details).

Fig. 14 .Fig. 15 .

 1415 Fig.14. Incidence rate of galaxies with O32 above the metallicitydependent threshold for three equally sized redshift-selected subsets, with median redshifts of z = 0.42 (green), z = 0.63 (purple), and z = 0.74 (orange) and the entire SDSS comparison sample (red) with median redshift z = 0.03. The vertical lines reflect the 1-σ errors in each bin.

  [O iii]λ4363 emission line. As such we can use the [O iii]λ4363/[O iii]λ5007 ratio to diagnose the electron temperatures (T e ) in the [O iii] regime of the H ii re-Paalvast et al.: O32 ratios of galaxies in MUSE data

  we show an HST image of the F775W filter and MUSE images of the [O ii] and [O iii] lines of the most extreme oxygen emitter in our sample, which has O32 = 23. This galaxy is observed in the UDF-mosaic at z = 0.83 (identified in the MUSE catalogue of Inami et al. 2017 as id = 6865). The spectrum of this object is shown in the upper panel of Fig. 1, from which we derived an extinction τ V = 0.49. The field is very crowded as can be seen in the HST image. The MUSE narrow band images of the [O ii]λ3727 and [O iii]λ5007 lines, however, confirm that the measured line fluxes originate from this source in the centre of the images. The log R23, log [Ne iii]/[O ii], and log [O iii]λ4363/[O iii]λ5007 of this object are shown by the diamond in Figs. 16, 17, and 18. The log [O iii]λ4363/[O iii]λ5007

Fig. 19 .

 19 Fig. 19. HST F77W image (left), a MUSE narrow band image of the [O ii]λ3727 line (middle), and a MUSE narrow band image of the [O iii]λ5007 line (right) of the most extreme oxygen emitter with O32 = 23.

  2018a

  Fig. A.3. SFR -M * diagram, where the SFR is corrected to z = 0 using the redshift evolution from Eq. 1, similar to Fig. 5. The symbols are coloured according to the survey from which the galaxies originate.

Fig. B. 2 .

 2 Fig. B.2. Fraction of O32 above the metallicity-dependent threshold for high O32 emitters versus the distance to the SFMS. The colours are used in the same way as in the previous figure. The upper panel shows the entire SDSS sample, the other panels for a M * range are specified in green. The median of the fraction of the entire sample from the upper left panel is presented by the black dashed line.

Table 1 .

 1 Number of galaxies in mass and redshift bins.7< log M * /M <8 8 < log M * /M < 9 9 < log M

	/M <10 total

*

  Fig. 17. Log O32 versus the log [Ne iii]λ3869/[O ii] ratio for our extreme emitters (O32 >4) with a [Ne iii]λ3869 detection of at least 3 σ. The galaxies of which we do not have a significant [O ii] detection are presented by red circles, which are the 3-σ upper limits on both ratios. The rest of the extreme emitters sample is shown by black circles and the most extreme O32 emitter by the black diamond. The grey line corresponds to the Levesque & Richardson (2014) relation for Z = 0.2 Z between log O32 and log [Ne iii]λ3869/[O ii] , which is offset by +0.6 in the y-direction.Fig. 18. Log O32 versus the [O iii]λ4363/[O iii]λ5007 ratio (lower xaxis) and the electron temperature in the [O iii] regime (lower x-axis) calculated using the nebular.ionic routine in Pyneb (Luridiana et al. 2015), assuming n e = 100 cm 3 . The coloured lines indicate the predictions from the Gutkin et al. (

		1.4	LR14 relation + 0.6
		1.2
	log O32	1.0
		0.8
		0.6
		0.8	0.6	0.4 log [NeIII]/[OII] 0.2 0.0	0.2
		1.2		1.4	1.6 10 -4 T e	1.8 2.0 2.2 2.4
		1.2 1.4	0.0001 0.0002 0.0005 0.001 0.002 0.004
	log O32	1.0
		0.8
		0.6
		2.0	1.9	1.8 log [OIII] 4363 /[OIII] 5007 1.7 1.6	1.5	1.4

Acknowledgements. We thank the referee for a constructive report that helped improve the paper. AV is supported by a Marie Heim-Vögtlin fellowship of the Swiss National Foundation. JB acknowledges support from the Fundação para a Ciência e a Technologia (FCT) through national funds (UID/FIS/04434/2013) and Investigador FCT contract IF/01654/2014/CP1215/CT0003., and by FEDER through COMPETE2020 (POCI-01-0145-FEDER-007672). TC acknowledges support from the ANR FOGHAR (ANR-13-BS05-0010-02), the OCEVU Labex (ANR-11-LABX-0060), and the A*MIDEX project (ANR-11-IDEX-0001-02) funded by the "Investissements d'avenir" French government programme. JS and SM acknowledge support from The Netherlands Organisation for Scientific Research (NWO), VICI grant 639.043.409. SC gratefully acknowledges support from Swiss National Science Foundation grant PP00P2_163824.

Appendix A: Distribution of galaxy properties

In addition to the histograms of the distances to the SFMS of the four different surveys (Fig. 3), we show here the M * (Fig. A.1) and the SFR (Fig. A.2) histograms. Because deep HST photometry is used a priori for the detection of galaxies in the UDF survey, the median M * of the UDF is at lower stellar mass than that of the other surveys. However, the median stellar mass of the Galaxy Groups survey is higher, at log M * /M ≈ 9.2, because this survey is aimed to observe galaxy groups that contain higher mass galaxies. The SFR distributions reflect the differences of the depth of the MUSE data of the different surveys. The median SFR of the shallowest survey, MUSE-Wide, is therefore approximately one order of magnitude larger than that of the other surveys.

In Fig.

A.3 we show originating survey of individual galaxies in a redshift-corrected SFR -M * diagram. Below this relation, galaxies from the MUSE QuBES and the MUSE-Wide surveys are more densely populated, which means that most of the galaxies that are not selected for the final sample are from these surveys.

A&A proofs: manuscript no. MUSE_O32 6 7 8 9 10 11 12

Appendix B: Stellar mass and flux completeness in SDSS data

Here we discuss the influence of flux incompleteness and M * incompleteness on the O32 ratio and the high O32 fractions.

To study how potential incompletenesses affect our results, we use star-forming galaxies selected from the SDSS DR7 sample, which is introduced in Sect. 4.1.1. We assume this sample is complete down to log M * /M = 9 (see Fig. 4 in Baldry et al. 2008). In Here we show the fraction of galaxies with O32 above the metallicity-dependent threshold, defined in Sect. 4.1.3, versus the distance to the SFMS. The top left panel shows the relations for all galaxies and in the other three panels we show the relations in different mass intervals. The upper right panel with 8.5 < log M * /M < 9.0 is mass-incomplete, however, if we compare the two lowest panels, we see that the trends are comparable to the median of the entire sample (dashed black line) and independent of the Hβ flux cut. This suggest that the trend that we derived between the fraction of high oxygen emitters and mass is not very sensitive to a possible mass incompleteness.