A&A 459, 85-101 (2006)
DOI: 10.1051/0004-6361:20065216

Gas metallicity diagnostics in star-forming galaxies[*]

T. Nagao1,2 - R. Maiolino1,3 - A. Marconi1


1 - INAF - Osservatorio Astrofisico di Arcetri, Largo Enrico Fermi 5, 50125 Firenze, Italy
2 - National Astronomical Observatory of Japan, 2-21-1 Osawa, Mitaka, Tokyo 181-8588, Japan
3 - INAF - Osservatorio Astrofisico di Roma, Via di Frascati 33, 00040 Monte Porzio Catone, Italy

Received 15 March 2006 / Accepted 10 August 2006

Abstract
Generally the gas metallicity in distant galaxies can only be inferred by using a few prominent emission lines. Various theoretical models have been used to predict the relationship between emission line fluxes and metallicity, suggesting that some line ratios can be used as diagnostics of the gas metallicity in galaxies. However, accurate empirical calibrations of these emission line flux ratios from real galaxy spectra spanning a wide metallicity range are still lacking. In this paper we provide such empirical calibrations by using the combination of two sets of spectroscopic data: one consisting of low-metallicity galaxies with a measurement of [O III]$\lambda $4363 taken from the literature, including spectra from the Sloan Digital Sky Survey (SDSS), and the other one consisting of galaxies in the SDSS database whose gas metallicity has been determined from various strong emission lines in their spectra. This combined data set constitutes the largest sample of galaxies with information on the gas metallicity available so far and spanning the widest metallicity range. By using these data we obtain accurate empirical relations between gas metallicity and several emission line diagnostics, including the R23 parameter, the [N II]$\lambda $6584/H$\alpha $ and [O III]$\lambda $5007/[N II]$\lambda $6584 ratios. Our empirical diagrams show that the line ratio [O III]$\lambda $5007/[O II]$\lambda $3727 is a useful tool to break the degeneracy in the R23 parameter when no information on the [N II]$\lambda $6584 line is available. The line ratio [Ne III]$\lambda $3869/[O II]$\lambda $3727 also results to be a useful metallicity indicator for high-z galaxies, especially when the R23 parameter or other diagnostics involving [O III]$\lambda $5007 or [N II]$\lambda $6584 are not available. Additional, useful diagnostics newly proposed in this paper are the line ratios of (H$\alpha $+[N II] $\lambda \lambda $6548,6584)/[S II]$\lambda $6720, [O III]$\lambda $5007/H$\beta $, and [O II]$\lambda $3727/H$\beta $. Finally, we compare these empirical relations with photoionization models. We find that the empirical R23-metallicity sequence is strongly discrepant with respect to the trend expected by models with constant ionization parameter. Such a discrepancy is also found for other line ratios. These discrepancies provide evidence for a strong metallicity dependence of the average ionization parameter in galaxies. In particular, we find that the average ionization parameter in galaxies increases by $\sim$0.7 dex as the metallicity decreases from 2 $Z_{\odot}$ to 0.05 $Z_{\odot}$, with a small dispersion. This result should warn about the use of theoretical models with constant ionization parameter to infer metallicities from observed line ratios.

Key words: galaxies: abundances - galaxies: evolution - galaxies: general - galaxies: ISM - H II regions

1 Introduction

The gas metallicity is one of the most important tools to investigate the evolutionary history of galaxies. This is because the gas metallicity of galaxies is basically determined by their star-formation history. Recent observational studies allowed the investigation of the gas metallicity even in high-z galaxies beyond z = 1, such as Lyman-break galaxies (e.g., Teplitz et al. 2000a,b; Pettini et al. 2001), submillimeter-selected high-z galaxies (Swinbank et al. 2004), and so on (see also, e.g., Kobulnicky & Kewley 2004; Savaglio et al. 2005; Maier et al. 2006; Liang et al. 2006; Erb et al. 2006). Such observational insights on the metallicity evolution of galaxies are now giving constraints on the theoretical understandings of the formation and the evolution of galaxies (e.g., Bicker et al. 2004).

However, metallicity measurements in distant galaxies are not straightforward. Information on the gas temperature is required for a precise determination of the gas metallicity, but the gas temperature can be accurately inferred only when the fluxes of auroral emission lines such as [O III]$\lambda $4363 and [N II]$\lambda $5755 are known, and these are generally too weak to be measured in faint distant galaxies. The measurement of the auroral emission lines is difficult even for galaxies in the local universe especially when the gas metallicity is high, because the collisional excitation of the auroral transitions is suppressed due to efficient cooling through far-infrared fine-structure emission lines (e.g., Ferland et al. 1984; Nagao et al. 2006a). Therefore, in general we have to rely on some relations between gas metallicity and flux ratios of strong emission-lines to estimate the gas metallicity in most galaxies. Extensive studies have been performed to calibrate such metallicity diagnostics by using only strong emission lines. One of the most frequently used metallicity diagnostics is the R23 parameter, defined as

\begin{displaymath}R_{23} =
\frac{ F({\rm [OII]}\lambda3727) +
F({\rm [OIII]...
...F({\rm [OIII]}\lambda5007) }
{ F({\rm H}\beta\lambda4861) },
\end{displaymath} (1)

where F([O II]$\lambda $3727), F([O III]$\lambda $4959) and so on denote the emission-line fluxes of [O II]$\lambda $3727, [O III]$\lambda $4959 and so on, respectively. The R23 was proposed by Pagel et al. (1979), and its calibration to the oxygen abundance has been improved by various photoionization model calculations (e.g., McGaugh 1991; Kewley & Dopita 2002).

Table 2: Statistical properties of the samples.

One serious problem of this indicator is that a certain value of R23 has two different solutions, a low-metallicity solution and a high metallicity one. Therefore additional, or alternative, diagnostics aimed at removing the R23 degeneracy have been proposed (e.g., Alloin et al. 1979; Denicoló et al. 2002; Kewley & Dopita 2002; Pettini & Pagel 2004). However, most of these methods exploit the [N II]$\lambda $6584 line, which has the problem of being very weak at sub-solar metallicities (hence difficult to measure) and the problem of being rapidly shifted outside the spectral band of many surveys at high redshift (e.g., unusable beyond $z \sim 0.5$ in optical spectra). On the theoretical side various models have been presented, which provide the ratios among the most prominent emission lines as a function of metallicity (e.g. Kewley & Dopita 2002). However, model predictions strongly depend on the assumed physical parameters of the ionized gas, and in particular on the ionization parameter ( $U \equiv \Phi_{\rm H} / (c n_{\rm H})$, where $\Phi_{\rm H}$ is the surface flux of hydrogen-ionizing photon and $n_{\rm H}$ is the gas density). As a consequence, an accurate correspondence between individual diagnostics (line ratios) and metallicity cannot be established, because of the lack of information on the physical conditions of the gas. Summarizing, many gas metallicity diagnostics proposed so far are either ambiguous or unusable when applied to the spectra of distant galaxies.

The goal of this paper is to obtain accurate, empirical calibrations between metallicity and individual diagnostics involving a few strong emission lines, which can be applied to the spectra of distant galaxies. In particular, we re-calibrate diagnostics already proposed in the past, but we also propose new diagnostics which appear particularly suited for distant galaxies. This work is obtained by combining two large data sets. The first one is composed of recent spectroscopic observations of low-metallicity galaxies ( $7.0 \la 12+\log($O/H$)\la8.5$), whose metallicity is accurately determined through the [O III]$\lambda $4363 line (Sect. 2.1). This dataset consists of two subsamples; one is taken from the database of the Sloan Digital Sky Survey (SDSS; York et al. 2000; Strauss et al. 2002) (Sect. 2.1.1) and the other is taken from the literature (Sect. 2.1.2). The second data set is a subsample of galaxies in the SDSS database, whose metallicity has been derived by Tremonti et al. (2004) ( $8.2\la 12+\log($O/H$)\la9.2$; see Sect. 2.2). These combined data sets provide the largest sample of galaxies with information on the gas metallicities and spanning more than 2 dex in metallicity.

2 Data

2.1 Spectroscopic data of low-metallicity galaxies with a [OIII]$\lambda $4363 measurement

2.1.1 SDSS data (sample A)

The gas-phase oxygen abundance is well determined when the flux of [O III]$\lambda $4363 is measured (e.g., Osterbrock 1989). Although such measurements have been performed for more than a hundred low-metallicity galaxies, simple compilation of those earlier results may introduce some unexpected biases and uncertainties. This is because the data were collected by various (heterogeneous) observations with different properties (aperture size, wavelength resolution, and so on) and because the method of calculating the oxygen abundance is different for different authors. Recently, Izotov et al. (2006b) reported their systematic measurements of the oxygen abundance for low-metallicity galaxies in the SDSS Data Release 3 (DR3; Abazajian et al. 2005) by using the [O III]$\lambda $4363 emission-line flux. The extinction-corrected emission-line fluxes of galaxies with a measurement of the oxygen abundance provided by Izotov et al. (2006b) are the ideal data for the empirical calibration of metallicity diagnostics, because the data were obtained and measured in a homogeneous way and because the oxygen abundance is also calculated with a common method. The number of spectra analyzed by Izotov et al. (2006b) is 309. Among them, we use the data with a relatively small error in the oxygen abundance [$\Delta$(log $\frac{\rm O}{\rm H}) \leq 0.05$] (146 spectra). Here we adopt the uncertainty [$\Delta$(log $\frac{\rm O}{\rm H}$)] given in Table 2 of Izotov et al. (2006). Since some spectra in the database of Izotov et al. (2006b) are duplicated for the same objects, the number of galaxies with $\Delta$(log $\frac{\rm O}{\rm H}) \leq 0.05$ is 139. Hereafter we call this sample "sample A''.


  \begin{figure}
\par\includegraphics[angle=-90,width=6cm,clip]{5216f01.eps}
\end{figure} Figure 1: Oxygen abundance of the galaxies in sample A, derived by Izotov et al. (2006b), as a function of redshift.
Open with DEXTER

However, this sample has two problems when used to accurately calibrate metallicity diagnostics. First, most of the galaxies in sample A have a relatively high oxygen abundance, and only 6 of them have $12+\log($O/H) < 7.6. Therefore the statistical reliability of the empirical calibration of metallicity diagnostics would be extremely poor at $12+\log($O/H) < 7.6if using only this sample. Second, in sample A, there is the remarkable tendency for lower-metallicity galaxies to have lower redshift. In Fig. 1, the oxygen abundance of galaxies in sample A is shown as a function of redshift. The origin of this apparent tendency is likely due to the fact that SDSS is not a volume-limited survey; that is, galaxies at higher redshift have preferentially higher luminosity, and thus higher metallicity. In particular, all of the galaxies with $12+\log($O/H) < 7.6are at z < 0.02. This means that the [O II]$\lambda $3727 flux cannot be measured for the latter galaxies due to the limited wavelength coverage of the SDSS spectroscopy ( $\lambda_{\rm obs} \ga 3800~{\rm\AA}$). Therefore, if using only sample A, we could not calibrate the diagnostics involving [O II]$\lambda $3727 [R23 and F([N II]$\lambda $6584)/F([O II]$\lambda $3727)] in the metallicity range $12+\log($O/H) < 7.6. This is a serious problem, because R23 is one of the most frequently used metallicity diagnostics and thus should be calibrated in a wide metallicity range. In conclusion, the accurate calibration of various metallicity diagnostics in a wide metallicity range cannot be achieved by using only sample A. We therefore collected additional data of [O III]$\lambda $4363-detected galaxies, which are described in the following subsection.

2.1.2 Other data from literature (sample B)

In order to increase the number of low-metallicity galaxies with a measurement of the oxygen abundance, we compiled the reddening-corrected emission-line flux data of galaxies with a [O III]$\lambda $4363 measurement from the literature. The sample of compiled galaxies is given in Table 1. For the objects whose spectroscopic properties have been reported by more than one paper independently, we chose the one with higher signal-to-noise ratio. When both the spectroscopic properties of the whole galaxy and of parts of it have been reported, we compiled both of them (e.g., Mrk 116). Consequently, the number of the compiled objects is 157. To minimize possible systematic errors owing to the different methods on the calculation of the oxygen abundance, we re-calculate their oxygen abundance by adopting the same method used for sample A (Izotov et al. 2006b). The re-calculated R23 parameter, gas density of [S II]-emitting region [$n_{\rm H}$(S+)], gas temperature of [O III]-emitting region [ $t{\rm (O^{2+})}$] and oxygen abundance [ $12+\log($O/H)] are given in Table 1, along with the reference to the data of the emission-line flux ratios. To calculate R23, we did not use F([O III]$\lambda $4959) but calculate the ratio of [F([O II]$\lambda $3727 $) + 1.327\times F$([O III]$\lambda $5007)] /F(H$\beta $) since only the flux of [O III]$\lambda $5007 (without that of [O III]$\lambda $4959) is given in some reference papers.


  \begin{figure}
\par\includegraphics[angle=-90,width=6cm,clip]{5216f02.eps}
\end{figure} Figure 2: Emission-line flux ratios of [S II]$\lambda $6717/[S II]$\lambda $6731 for the compiled low-metallicity galaxies, as a function of the oxygen abundance derived by us (see Sect. 2.1.2). The upper horizontal dotted line denotes the theoretical low-density limit of this flux ratio and the lower dotted line denotes the high-density limit.
Open with DEXTER


  \begin{figure}
\par\includegraphics[angle=-90,width=6cm,clip]{5216f03.eps}
\end{figure} Figure 3: Oxygen abundances of the compiled low-metallicity galaxies re-calculated by us (see Sect. 2.1.2) are plotted as a function of the oxygen abundances given in the original references. Dotted line is not the best-fit line but a reference line for the case when the two quantities are the same.
Open with DEXTER

To check whether our adopted method causes possible systematic difference in the oxygen abundance from the values given in the original references, we compare the oxygen abundances re-calculated by us and those given in the original papers in Fig. 3. Apparently, there is no systematic difference between our results and the results given in the literature. The mean and the RMS of the difference, [$12+\log($O/H $)_{\rm this \ work}]
- [12+\log($O/H $)_{\rm literature}$], are +0.001 and 0.041, respectively. This mean value of the difference is smaller than the typical error of the re-calculated oxygen abundance. Among the 157 objects given in Table 1 and plotted in Fig. 3, we use only data with $\Delta$(log $\frac{\rm O}{\rm H}) \leq 0.05$. This constraint results in a sample of 120 objects, that is hereafter called "sample B''. Note that the mean and the RMS of the difference, [$12+\log($O/H $)_{\rm this \ work}]
- [12+\log($O/H) $_{\rm literature}$], for sample B are -0.008 and 0.038. Again the mean value of the difference is smaller than the typical error on the oxygen abundance.

2.2 Spectroscopic data of high-metallicity galaxies (sample C)

For high-metallicity galaxies, we referred to the oxygen abundance derived by Tremonti et al. (2004), who derived the metallicities of $\sim$53 000 galaxies in the SDSS database. They used the fluxes of many strong emission lines ([O II]$\lambda $3727, H $\beta \lambda$4861, [O III]$\lambda $5007, H$\alpha $, [N II]$\lambda $6584, [S II]$\lambda $6717, and [S II]$\lambda $6731) and comparing them with photoionization models (Ferland et al. 1998). Although they presented the results of their analysis on the spectra of the SDSS Data Release 2 (DR2; Abazajian et al. 2004), they also provide the results of their recent analysis on the spectra from Data Release 4 (DR4; Adelman-McCarthy et al. 2006) on their web site[*]. Their estimate of the oxygen abundance does not rely only on a single metallicity diagnostic flux ratio, but uses all the optical prominent emission lines (see also, Charlot & Longhetti 2001; Brinchmann et al. 2004). Therefore, among galaxies without [O III]$\lambda $4363 flux, their sample is currently the best one in terms of both sample size and reliability.

The oxygen-abundance catalog of the SDSS DR4 galaxies contains 567 486 objects. The objects in the catalog are classified into the five classes; "star-forming galaxies'', "low S/N star-forming galaxies'', "composite'', "active galactic nuclei (AGNs)'', and "unclassificable''. We referred only to the galaxies belonging the first class (141 317 objects). The emission-line fluxes of these galaxies are obtained from the emission-line flux catalog of the SDSS DR4 galaxies provided on the same web site as the oxygen-abundance catalog. The emission-line fluxes given in this catalog were measured from the stellar-continuum subtracted spectra with the latest high spectral resolution population synthesis models by Bruzual & Charlot (2003), and thus more reliable than the flux data provided on the SDSS Data Archive Server. Since the emission-line fluxes given in their catalog are not corrected for dust extinction, we corrected them by using the Balmer decrement method with the reddening curve of Cardelli et al. (1989). We then removed the duplicated objects and the objects observed in some problematic plates (see the SDSS web page[*]) from the cross sample  of the oxygen-abundance catalog and the emission-line flux catalog. Then we select only objects satisfying all of the following five criteria:

1.
The redshift is higher than 0.028;
2.
Both H$\alpha $ and H$\beta $ emission lines have S/N $\geq$ 10;
3.
log $\frac{F({\rm [OIII]\lambda5007})}{F({\rm H\beta})}
> \frac{0.61}{{\rm log}\frac{F({\rm [NII]\lambda6584})}
{F({\rm H\alpha})}-0.05} + 1.3$;
4.
The fiber aperture covers at least 20% of the total $g^\prime$-band photons;
5.
The uncertainty on the estimated stellar mass is less than 0.5 dex (i.e., log $M_{97.5} -\log~M_{2.5} <
0.5$, where M97.5 and M2.5 are the 97.5th and 2.5th percentile of probability distribution of the estimated stellar mass; see Kauffmann et al. 2003b for more details).
The first criterion ensures the accurate measurement of F([O II]$\lambda $3727). The second criterion selects emission-line galaxies. The S/N values are taken from the emission-line flux catalog; however, the errors in this catalog are underestimated with a factor of $\sim$2 (see the web site for more details). The third criterion is required to reject AGNs from the sample, following Kauffmann et al. (2003a), although the AGN removal had been already examined in the process of checking the classification in the oxygen-abundance catalog mentioned above. The fourth criterion is required to avoid significant aperture effects on the flux ratios. Kewley et al. (2005) reported that the introduced systematic error of the metallicity determination reaches up to $\sim$40% when the fiber covers only 20% of the total $g^\prime$-band photons (see also Tremonti et al. 2004). The fifth criterion is not directly relevant to this study, but is required in our companion paper (Nagao et al. 2006b).

Note that in sample C we do not put constraints on $\Delta$(log $\frac{\rm O}{\rm H})$. This is to prevent sample C to be devoid of galaxies at metallicities $12+\log($O/H) < 7.6 so that sample C and sample A+B have some overlap in terms of metallicities. Indeed, due to the luminosity-metallicity relation, low metallicity galaxies in sample C are on average fainter and therefore tend to have larger errors. However, the averaged gas metallicities in sample C are most likely reliable even at low metallicities, thanks to the large number of objects in the sample. Our final sample of emission-line galaxies consists of 48 497 objects, which is hereafter called "sample C''. The redshift distribution of galaxies in sample C is shown in Fig. 4. Its median value is 0.085, while the mean and the RMS are 0.092 and 0.040, respectively. The distribution of the oxygen abundance of galaxies in sample C is shown in Fig. 5. Its median value is 9.016, while the mean and the RMS are 8.976 and 0.166. The means and the RMSs of the oxygen abundance of samples A, B, and C are summarized in Table 2.


  \begin{figure}
\par\includegraphics[angle=-90,width=6cm,clip]{5216f04.eps}
\end{figure} Figure 4: Frequency distribution of the redshift of the SDSS DR4 galaxies after our sample selection (sample C) described in Sect. 2.2. Galaxies at z<0.028 are not included (see text).
Open with DEXTER


  \begin{figure}
\par\includegraphics[angle=-90,width=6cm,clip]{5216f05.eps}
\end{figure} Figure 5: Frequency distribution of the oxygen abundance of galaxies in sample C.
Open with DEXTER


  \begin{figure}
\par\includegraphics[angle=-90,width=6.5cm,clip]{5216f06.eps}
\end{figure} Figure 6: Emission-line flux ratios of R23 (=[F([O II]$\lambda $3727 $) + 1.327\times F$([O III]$\lambda $5007)]/F(H$\beta $)), F([N II]$\lambda $6584)/ F(H$\alpha $), and F([O III]$\lambda $5007)/F([N II]$\lambda $6584) for galaxies in sample A (red triangles), in sample B (green circles) and in sample C (black dots), as a function of the oxygen abundance. The compiled low-metallicity galaxies with an error of the oxygen abundance larger than 0.05 dex are not plotted. Dotted lines denote the solar metallicity [$12+\log($O/H) = 8.69].
Open with DEXTER


  \begin{figure}
\par\includegraphics[angle=-90,width=6cm,clip]{5216f07.eps}
\end{figure} Figure 7: Same as Fig. 6 but for the emission-line flux ratios of F([N II]$\lambda $6584)/F([O II]$\lambda $3727), F([N II]$\lambda $6584)/F([S II]$\lambda $6720), and F([O III]$\lambda $5007)/F([O II]$\lambda $3727).
Open with DEXTER

3 Results

In Figs. 6 and 7, we plot emission-line flux ratios for the galaxies in samples A, B, and C. To avoid noisy objects in sample C, we consider only those with $S/N \geq 10$ (cataloged value) for all the related emission lines (e.g., H$\beta $, [O II]$\lambda $3726, [O II]$\lambda $3729 and [O III]$\lambda $5007 for the case of R23). In addition to R23, all the other flux ratios investigated here are metallicity-sensitive flux ratios and sometimes regarded as metallicity diagnostics (see, e.g., Kewley & Dopita 2002; Pettini & Pagel 2004; Kobulnicky & Kewley 2004). Among them, F([O III]$\lambda $5007)/F([O II]$\lambda $3727) is sensitive also to the ionization parameter and thus it has not been regarded as a good metallicity diagnostic flux ratio (see Kewley & Dopita 2002). Instead, this flux ratio has been used to investigate the ionization parameter, and has been sometimes used in the following form:

\begin{displaymath}O_{32} =
\frac{ F({\rm [OIII]}\lambda4959) +
F({\rm [OIII]}\lambda5007) }
{ F({\rm [OII]}\lambda3727) }
\end{displaymath} (2)

(e.g., Kobulnicky & Kewley 2004).

The data sequences in the diagnostics-metallicity diagrams are mostly continuous for different samples, and accordingly the whole sample shows clear relations between various metallicity diagnostics and the oxygen abundance. Since there are no apparent systematic differences in the diagnostics-metallicity sequences between sample A and sample B, we combined these two samples and identified as "sample A+B'' hereafter in order to improve the statistics at low metallicities. The statistical properties of the sample A+B are given in Table 2. The diagram of F([N II]$\lambda $6584)/F([S II]$\lambda $6720) versus the oxygen abundance shows an apparent discontinuity between sample A+B and sample C, where F([S II]$\lambda $6720) denotes the sum of F([S II]$\lambda $6717) and F([S II]$\lambda $6731). We will discuss the issue of this discontinuity in Sect. 4.1.

The diagram of R23 versus the oxygen abundance shows a $\cap$-shaped distribution with a peak at $12+\log($O/H $)
\sim 8.0$. This is consistent with the previous studies on the empirical relation between R23 and the oxygen abundance based on smaller samples of observational data (e.g., Edmunds & Pagel 1984; McGaugh 1991; Miller & Hodge 1996; Castellanos et al. 2002; Lee et al. 2003a; Bresolin et al. 2004, 2005; Pilyugin & Thuan 2005). As discussed in Sect. 4.3, however, this appears to be systematically different from previous predictions of photoionization models.

To investigate the relation between the flux ratios and the oxygen abundances quantitatively, we calculate the means and the RMSs of the flux ratios of galaxies within bins of oxygen abundance. For sample A+B, we calculate them in the range  $7.05 < 12+\log($O/H) < 8.55 with a bin width of $\Delta$[log(O/H)] = 0.1, except at lowest and highest oxygen abundances, where the bin width is wider (i.e., $7.05 < 12+\log($O/H) < 7.45 and $8.35 < 12+\log($O/H) < 8.55) due to the small number of sources in these ranges. All of the metallicity bins contain at least 6 galaxies. The results are given in Table 3. We also calculate the mean and the RMS of flux ratios for sample C in the range $8.15 < 12+\log($O/H) < 9.25 with a bin width of $\Delta$[log(O/H)] = 0.1 dex. The results are given in Table 4.

Table 3: Means and RMSs of emission-line flux ratios of the galaxies in the sample A+Ba.

Table 4: Means and RMSs of emission-line flux ratios of the galaxies in sample Ca.


  \begin{figure}
\par\includegraphics[angle=-90,width=6.5cm,clip]{5216f08.eps}
\end{figure} Figure 8: Same as Fig. 6 but means and the RMS values are shown in each bin of oxygen abundance, instead of individual data. Filled squares and filled stars denote the mean flux ratios for galaxies in sample C and those for galaxies in sample A+B, respectively. The errorbar denotes the RMS. The dashed line denotes the best-fit polynomial function, as described in the text. Dashed lines denote the solar metallicity [$12+\log($O/H) = 8.69].
Open with DEXTER


  \begin{figure}
\par\includegraphics[angle=-90,width=6.5cm,clip]{5216f09.eps}
\end{figure} Figure 9: Same as Fig. 8 but for the emission-line flux ratios of F([N II]$\lambda $6584)/F([O II]$\lambda $3727), F([N II]$\lambda $6584)/F([S II]$\lambda $6720), and F([O III]$\lambda $5007)/F([O II]$\lambda $3727).
Open with DEXTER

The calculated mean and the RMS of the flux ratios for each metallicity bin are shown in Figs. 8 and 9. We then fit the observed sequences between flux ratios and oxygen abundance with polynomial functions in the range $7.05 < 12+\log($O/H) < 9.25 (or $0.02 < Z_{\rm gas}/Z_\odot < 4$), and the results of the fits are also shown in Figs. 8 and 9. We decided to fit 3rd-order polynomial functions for the binned data, not for the individual data, in order to avoid giving too much weights to the high metallicity range (where most of the data are). In Table 5, the coefficients of the best-fit polynomial functions are provided, according to the formula

\begin{displaymath}{\rm log} \ (F_1 / F_2) =
\sum_{N}^{}{a_N [{\rm log}(Z/Z_\odot)]^N}
\end{displaymath} (3)

where F1 / F2 is the line flux ratio (or the R23parameter) and N= (0, 1, 2, 3). Here we adopt $12+\log($O/H) $_\odot = 8.69$ for $Z_{\odot}$ (Allende Prieto et al. 2001). For the convenience of the reader, in Table 6 we also give the coefficients of the best-fit polynomial functions in the following form:

\begin{displaymath}{\rm log} \ (F_1 / F_2) =
\sum_{N}^{}{b_N [12 + {\rm log (O/H)}]^N}.
\end{displaymath} (4)

The expected uncertainties on the derived metallicities from the diagnostic flux ratios calibrated here can be estimated using the RMS values of the flux ratios given for each mass bin (Tables 3 and 4). By looking at the RMS plotted in Figs. 8 and 9, we can recognize that, for instance, the diagnostic flux ratios give highly uncertain metallicities when F([N II]$\lambda $6584)/F([O II]$\lambda $3727 $)
\la 0.05$, F([O III]$\lambda $5007)/F([O II]$\lambda $3727$)
\ga 2$, and F([N II]$\lambda $6584)/ F([S II]$\lambda $6720$)
\la 0.3$. On the contrary, the metallicity is well determined ( $\Delta Z \la 0.2$ dex) when F([O III]$\lambda $5007)/F([N II]$\lambda $6584$)
\la 10$ and F([N II]$\lambda $6584)/F([O II]$\lambda $3727 $)
\ga 0.05$.

Table 5: Coefficients of the best-fit polynomials for the observed relations between the emission-line flux ratios and the oxygen abundance, where $\log ~R = a_0 + a_1 x + a_2 x^2 + a_3 x^3$ [$x \equiv $ log(Z/ $Z_\odot) \equiv 12+\log($O/H)-8.69].

Finally we recall that these relations are valid only in the range $7.05 < 12+\log($O/H) < 9.25(which is however much wider than in any previous work). We warn on the use of these relations outside such metallicity range, since it would rely only on their extrapolation.

4 Discussion

4.1 Consistency between sample A+B and sample C

Before interpreting the results, we discuss on the consistency of the two main samples, i.e., galaxies with (A+B) and without (C) [O III]$\lambda $4363 measurements. As mentioned in Sect. 3, the relation between some emission-line flux ratios and the oxygen abundance is not smoothly connected between the two samples (A+B and C), and this is especially significant for the flux ratios of F([O III]$\lambda $5007)/F([O II]$\lambda $3727) and F([N II]$\lambda $6584)/F([S II]$\lambda $6720), but is also seen in other cases (Figs. 8 and 9). One of the possible reasons for this discrepancy is a systematic error in the estimate of the oxygen abundance for one (or both) of the two different methods, which in one case consists in using the gas temperature inferred through [O III]$\lambda $4363 emission (see Sect. 2.1) and in the other case is using all of optical strong emission lines (Tremonti et al. 2004).

Kobulnicky et al. (1999) investigated a possible systematic error in the former method, that is, the gas temperature may be overestimated through the [O III]$\lambda $4363 emission and thus the oxygen abundance may tend to be underestimated accordingly. This is because the strength of the [O III]$\lambda $4363 emission significantly depends on the gas temperature and thus spectra obtained by a global aperture toward a galaxy are biased towards higher gas-temperature H II regions (see also Peimbert 1967). According to their analysis, the overestimation of the gas temperature could be more serious in low-metallicity systems and could reach up to $\Delta T_{\rm e} = 1000{-}3000$ K, which results in the systematic underestimation of the oxygen abundance of 0.05-0.2 dex. However, although this effect may partly account for the discrepancy of the metallicity dependence of F([N II]$\lambda $6584)/F([S II]$\lambda $6720), it goes in the opposite direction to account for the discrepancy seen in F([O III]$\lambda $5007)/F([N II]$\lambda $6584) and F([O III]$\lambda $5007)/ F([O II]$\lambda $3727). Therefore, the effect of the biased temperature measurement is not the dominant origin of the discontinuities seen in Figs. 8 and 9.

A systematic error in the oxygen abundance may exist in the method of Tremonti et al. (2004). They estimated the oxygen abundance by comparing photoionization models with some optical emission-line fluxes, which were measured on the spectra after subtraction of the stellar component. Although their method of the stellar-component subtraction is a sophisticated one and uses the most recent population synthesis models of Bruzual & Charlot (2003), it is not clear whether the measurement of emission lines lying on the deep and complex stellar absorption features is completely free from some possible systematic errors. A possible improper subtraction of the stellar absorption features may lead to systematic errors on the fluxes of Balmer lines, which might result in a systematic error in the estimation of the gas metallicity of galaxies in sample C. The subtraction of stellar absorption features may be inaccurate also in sample A+B. For instance, in some earlier works the stellar subtraction was performed by simply assuming EW(H$\beta $) $_{\rm abs} = 1~{\rm\AA}$. This over-simplified assumption may introduce systematic errors in the derived gas metallicity and the emission-line flux ratios given in Table 1. Another possible source of uncertainty in the method of Tremonti et al. (2004) is the use of the [N II]$\lambda $6584 flux and its comparison to models. Most photoionization models assume that the relative nitrogen abundance scales with the metallicity linearly when the primary nitrogen creation dominates, and scales quadratically when the secondary nitrogen creation is dominant. However, the transition metallicity between the two modes is uncertain. An inaccurate value of the transition metallicity (which is indeed uncertain) may lead to systematic errors in the estimation of metallicity especially at low metallicities, which could be one of the possible origin of the discrepancy seen in Figs. 8 and 9.

The discrepancy in the metallicity dependences of emission-line flux ratios may also be a consequence of the selection of spectroscopic targets. While galaxies in sample C are basically selected in terms of their apparent magnitude and thus not largely biased toward any specific population, galaxies in sample B could be biased toward very strong emission-line galaxies (galaxies in sample A are in a composite situation; see Izotov et al. 2006). This is because the motivation behind most of the original observations, such as the studies on the primordial helium abundance (see the original references given in Table 1), required very accurate measurements of emission-line flux ratios. For a given metallicity, galaxies with stronger emission lines tend to be characterized by a higher ionization parameter, which may result into larger flux ratios of F([O III]$\lambda $5007)/F([N II]$\lambda $6584) and F([O III]$\lambda $5007)/F([O II]$\lambda $3727), although the difference in the ionization parameter should not cause a significant difference in the ratio of F([N II]$\lambda $6584)/F([S II]$\lambda $6720). We will discuss the effect of the ionization parameter on the discrepancy further in Sect. 4.3.

Actually some or all of the above matters could contribute to the discontinuity in the metallicity dependences of emission-line flux ratios, and their discrimination or their accurate correction are not feasible. We thus simply adopt the results of the fit described in Sect. 3 and not take the effects of the possible systematic errors into account in the following discussion. However, it should be noted that this rather complex situation is caused by relying on two independent methods to measure the oxygen abundance. This problem will be solved if a large sample of galaxies with a wide range of the oxygen abundance is investigated by using a unique method throughout the concerned metallicity range.

4.2 Comparison with previous empirical calibrations


  \begin{figure}
\par\includegraphics[angle=-90,width=8cm,clip]{5216f10.eps}
\end{figure} Figure 10: Comparison of our results with the previous empirical metallicity calibrations for the R23 parameter. Solid red line denotes our calibration. Blue, green, and magenta lines denote the calibration given by Tremonti et al. (2004), Edmunds & Pagel (1984), and Zaritsky et al. (1994), respectively. Symbols and errorbars are the same as those in Fig. 8. Vertical dotted line denotes the solar metallicity [$12+\log($O/H) = 8.69].
Open with DEXTER

We compare the results of our calibrations with previous empirical calibrations. In particular, in Fig. 10, we compare the empirical calibrations of R23 derived by us with those obtained by Tremonti et al. (2004), Edmunds & Pagel (1984), and Zaritsky et al. (1994). While there is a reasonable agreement between our result and the result from previous calibration for the lower branch (Edmunds & Pagel 1984), there are some systematic discrepancies for the upper branch. We should in particular discuss the difference between our calibration and that of Tremonti et al. (2004), since our calibration in the high-metallicity range is based on the metallicity of the SDSS galaxies (in sample C) derived by Tremonti et al. (2004). The calibration by Tremonti et al. (2004), which is provided only for the upper branch, is clearly flatter than ours. This discrepancy may be ascribed to the combination of various possible factors. Our calibration also includes the fit of the new sample of [O III]$\lambda $4363-detected galaxies, which are not included in Tremonti et al. (2004), and this is certainly one of the reasons for the discrepancy. However, the latter issue cannot completely account for the discrepancy, since the Tremonti et al. (2004) calibration fails to reproduce the SDSS data at $12+\log($O/H) < 8.5 (as shown in Fig. 10). It is likely that an additional source of the discrepancy is the different strategy of fitting the analytical function to the data. While we fit the third polynomial function to the binned data, Tremonti et al. (2004) fit the function to the whole sample of individual SDSS galaxies. Since the number of high metallicity galaxies [$12+\log($O/H) > 8.5] is much larger than the low metallicity sub-sample [$12+\log($O/H) < 8.5] as shown in Fig. 5, the analytical fit of Tremonti et al. (2004) is dominated by the high-metallicity part of the R23 diagram. Finally, the modest discrepancy at high metallicities [$12+\log($O/H) > 9] may be partly attributed to the difference in the sample selection criteria. As described in Sect. 3, we select the galaxies in the sample C with $S/N \geq 10$ for all of the lines H$\beta $, [O II]$\lambda $3726, [O II]$\lambda $3729, and [O III]$\lambda $5007 (note that the [O II] doublet lines are measured separately in the original catalog we used), while Tremonti et al. (2004) adopted the S/N criteria only for H$\beta $, H$\alpha $ and [N II]$\lambda $6584, not for [O II]$\lambda $3726, [O II]$\lambda $3729, and [O III]$\lambda $5007. Our selection criteria may preferentially reject objects at high metallicities with respect to those of Tremonti et al. (2004), because forbidden lines such as [O II] and [O III] become weak when gas metallicity is high due to the suppressed collisional excitation mechanism (e.g., Ferland et al. 1984; Nagao et al. 2006a). This effect may result in our selective choice of objects with strong [O II] and [O III] emission in a given metallicity bin, which could make our calibration to be steeper at high metallicities. Since the difference in the calibration between ours and that of Tremonti et al. (2004) is significant at $12+\log($O/H) > 9, it is suggested that our calibration for the R23 may overestimate the gas metallicity at $12+\log($O/H) > 9 by a factor of  $\Delta Z \sim 0.1$ dex at $12+\log($O/H $) \sim 9.1$.

The calibration of the diagnostic flux ratio F([N II]$\lambda $6584)/F(H$\alpha $) is especially important, because wavelength separation of the two lines is small (i.e., not sensitive to dust reddening and requiring only small wavelength coverage) and thus it is used as a diagnostic of the gas metallicity of galaxies at $z \la 2$ (e.g., Erb et al. 2006). In Fig. 11, we compare the empirical calibrations of F([N II]$\lambda $6584)/F(H$\alpha $) derived by us, with those derived by Pettini & Pagel (2004) and Denicoló et al. (2002). Our result agree reasonably well with Denicoló et al. (2002) only at sub solar metallicity, and there is a systematic difference in the slope between our result and the result reported by Pettini & Pagel (2004). The latter discrepancy may be due to the reduced metallicity range of the sample of Pettini & Pagel (2004), indeed most of their objects are distributed within $7.7 < 12+\log($O/H) < 8.5. However, the difference is significant ( $\Delta Z \ga 0.2$ dex) only at metallicities $12+\log($O/H) < 7.5 and $12+\log($O/H) > 8.5. Although the difference in the lowest-metallicity range is not a serious problem (because in this metallicity range the expected [N II]$\lambda $6584 flux is extremely weak and thus its measurement would be very challenging and probably inaccurate), it is important to pay attention to the difference in the high-metallicity range. Note that such "high-metallicity'' range (where the discrepancy with Pettini & Pagel 2004 occurs) is not so metal rich - the metallicity $12+\log($O/H) = 8.5corresponds to $Z = 0.65~Z_\odot$, still in the sub-solar metallicity domain.


  \begin{figure}
\par\includegraphics[angle=-90,width=8cm,clip]{5216f11.eps}
\end{figure} Figure 11: Comparison of our results with the previous empirical metallicity calibrations for F([N II]$\lambda $6584)/F(H$\alpha $). Solid red line denotes our calibration. Blue and magenta lines denote the calibration given by Pettini & Pagel (2004) and Denicoló et al. (2002), respectively. Symbols and errorbars are the same as those in Fig. 8. Vertical dotted line denotes the solar metallicity [$12+\log($O/H) = 8.69].
Open with DEXTER


  \begin{figure}
\par\includegraphics[angle=-90,width=8cm,clip]{5216f12.eps}
\end{figure} Figure 12: Comparison of our results with the previous empirical metallicity calibrations for F([O III]$\lambda $5007)/F([N II]$\lambda $6584). Solid red line denotes our calibration, and blue line denotes the calibration given by Pettini & Pagel (2004). Symbols and errorbars are the same as those in Fig. 8. Vertical dotted line denotes the solar metallicity [$12+\log($O/H) = 8.69].
Open with DEXTER

In Fig. 12, we compare the empirical calibrations of F([O III]$\lambda $5007)/F([N II]$\lambda $6584 derived by us with the one derived by Pettini & Pagel (2004). The difference between the two calibrations is more serious than that seen in Fig. 11 in the low metallicity range, $12+\log($O/H) < 8. However Pettini & Pagel (2004) correctly mentioned that the flux ratio F([O III]$\lambda $5007)/F([N II]$\lambda $6584 is of little use when F([O III]$\lambda $5007)/F([N II] $\lambda 6584
\ga 100$ because of the saturation of this diagnostic. The behavior of this diagnostic flux ratio in the low-metallicity range would be important to derive the upper limits on the metallicity from an upper limit of the [N II]$\lambda $6584 flux.

4.3 Comparison with photoionization models


  \begin{figure}
\par\includegraphics[angle=-90,width=6.5cm,clip]{5216f13.eps}
\end{figure} Figure 13: The averaged flux ratios and the best-fit polynomial functions of the metallicity dependence of R23, F([N II]$\lambda $6584)/F(H$\alpha $) and F([O III]$\lambda $5007)/F([N II]$\lambda $6584) (dashed line) are compared with the predictions of photoionization models (solid lines: Kewley & Dopita 2002; Kobulnicky & Kewley 2004). The lines with a digit 1, 2, 3, 4, 5, 6, and 7 denote the model predictions with the ionization parameter of log U = -3.8, -3.5, -3.2, -2.9, -2.6, -2.3, and -2.0, respectively. Dotted line denotes the solar metallicity [$12+\log($O/H) = 8.69].
Open with DEXTER


  \begin{figure}
\par\includegraphics[angle=-90,width=6.5cm,clip]{5216f14.eps}
\end{figure} Figure 14: Same as Fig. 13 but for the emission-line flux ratios of F([N II]$\lambda $6584)/F([O II]$\lambda $3727), F([N II]$\lambda $6584)/F([S II]$\lambda $6720), and F([O III]$\lambda $5007)/F([O II]$\lambda $3727).
Open with DEXTER

To interpret the metallicity dependences of the emission-line flux ratios, we compare observational data with the predictions of photoionization models. In Figs. 13 and 14, we show the empirical metallicity dependences and the theoretical metallicity dependences of some metallicity diagnostics, where the latter are taken from Kewley & Dopita (2002) except for F([N II]$\lambda $6584)/F(H$\alpha $) that is taken from Kobulnicky & Kewley (2004). Since the explicit analytic expression for the metallicity dependence of F([O III]$\lambda $5007)/F([O II]$\lambda $3727) is not given by Kewley & Dopita (2002), we derive the polynomial expression of the theoretical metallicity dependence by fitting the results given in Table 2 of Kewley & Dopita (2002). The photoionization models presented by Kewley & Dopita (2002) and Kobulnicky & Kewley (2004) were calculated by the photoionization code MAPPINGS III (Sutherland & Dopita 1993) combined with the stellar population synthesis codes PEGASE (Fioc & Rocca-Volmerange 1997) and STARBURST99 (Leitherer et al. 1999), for the range $7.6 < 12+\log($O/H) < 9.4. They assume that stars and gas have the same metallicity, which is a reasonable assumption given that photoionization is due to hot, young stars, presumably recently formed from the same gas that they are photoionizing. In their calculations, nitrogen is assumed to be a secondary nucleosynthesis element at $12+\log($O/H) > 8.3, and a primary nucleosynthesis element at lower metallicity. Effects of dust grains on the depletion of gas-phase heavy elements and on the radiative transfer are consistently taken into account. Their calculations cover the range of ionization parameters $-3.8 \leq {\rm log} \ U \leq -2.0$, or equivalently, $5\times10^6~{\rm cm ~s}^{-1} \leq q
\leq 3\times10^8~{\rm cm~s}^{-1}$ (where $U \equiv q/c$). See Kewley & Dopita (2002) for details on the calculations. Note that they adopted $12+\log($O/H) $_\odot = 8.93$ (Anders & Grevesse 1989) and expressed the metallicity in units of $Z_{\odot}$ ($12+\log($O/H) $_\odot = 8.93$). However, since we adopt a more recent value for the solar abundance, $12+\log($O/H) $_\odot = 8.69$ (Allende Prieto et al. 2001), the metallicity notation is different when the $Z_{\odot}$ unit is used, which should be kept in mind to compare our results with their predictions.

The most remarkable matter in the comparison between the empirical and theoretical metallicity dependences of emission-line flux ratios is the significant discrepancy in the theoretically-expected R23-sequence with respect to the observed trend. This is especially significant at low metallicity range $12+\log($O/H) < 8. Shi et al. (2006) also recently reported that a previous theoretical calibration of R23 (see McGaugh 1991; Kobulnicky et al. 1999) overpredicts the gas metallicity with respect to the metallicity measured through the gas temperature determined with [O III]$\lambda $4363 line ( $\Delta Z \sim 0.2$ dex), especially in the low metallicity range (i.e., $12+\log($O/H) < 8). This discrepancy is not due to an improper compilation in our data, because it has been reported also in the earlier works that the empirical peak of R23 is seen around 12+(O/H $)
\sim 8.0$, as mentioned already in Sect. 3. The discrepancy cannot be ascribed to problems to the model results of Kewley & Dopita (2002) either, because other theoretical works also predict higher peak metallicity of R23 independently [$12+\log($O/H$) \ga 8.3$; e.g., Kobulnicky et al. 1999]. One possible idea to reconcile this discrepancy is that the ionization parameter of the gas is higher than the parameter range which Kewley & Dopita (2002) covers, especially in low-metallicity objects. If the ionization parameter correlates negatively with the gas metallicity and it reaches up to log U > -2 at the lowest metallicities, photoionization models would predict larger values of R23 in the lower-metallicity range with respect to constant-U models. This idea appears to be consistent with the behaviors of the empirical sequences in the U-sensitive flux ratios, F([O III]$\lambda $5007)/F([N II]$\lambda $6584) (Fig. 13) and F([O III]$\lambda $5007)/F([O II]$\lambda $3727) (Fig. 14). By focusing on these two U-sensitive flux ratios, we can see that the ionization parameter increases by $\sim$0.7 dex with decreasing oxygen abundance from $12+\log($O/H) = 9.0 to 7.5, supporting the above interpretation. Although the absolute value of the required ionization parameter appear to be inconsistent between R23 and the latter two U-sensitive flux ratios, the inferred absolute U values depends also on some model assumptions such as the spectral energy distribution (SED) of ionizing photons or the relative elemental abundance ratios, which also change as a function of metallicity. We thus conclude that the metallicity dependence of the ionization parameter (hereafter "U-Z relation'') causes the discrepancy between the empirical R23 distribution and the model predictions with a constant ionization parameter.

Note that the F([O III]$\lambda $5007)/F([N II]$\lambda $6584) and F([O III]$\lambda $5007)/F([O II]$\lambda $3727) ratios are also sensitive to the hardness of the ionizing radiation, which is a strong function of the stellar metallicity. This effect can in principle also contribute to the dependence of F([O III]$\lambda $5007)/F([N II]$\lambda $6584) and F([O III]$\lambda $5007)/F([O II]$\lambda $3727) ratios on metallicity. However, the models by Kewley & Dopita (2002) plotted in Figs. 13 and 14 already take into account the hardening of the stellar spectra as a function of metallicity. Therefore, the discrepancy between constant-U models and the data indicates that the hardening of the ionizing spectra must be associated with a variation of U with metallicity. In particular, the dependence of the F([O III]$\lambda $5007)/F([O II]$\lambda $3727) ratio on metallicity cannot entirely be ascribed only to the hardening of the ionizing radiation, but also to a U-Z relation.

The inferred U-Z relation is a very interesting result. Maier et al. (2006) also recently reported that the lower-metallicity galaxies tend to be characterized by a higher ionization parameter (see also Maier et al. 2004, who reported the correlation between the absolute B magnitude and the flux ratio of F([O III]$\lambda $5007)/F([O II]$\lambda $3727) among galaxies in the local universe). Although a detailed theoretical interpretation of this empirical relation goes beyond the scope of this paper, in the following we discuss two possible qualitative interpretations. One possible origin of this effect may be associated with the mass-metallicity relation and with the mass-age relation in local galaxies. According to these relations, higher metallicity galaxies are associated with more massive and older systems. H II regions ionized by later stellar populations are expected to be characterized by lower ionization parameters, due to the lower luminosity of the ionizing stars. Another possible explanation may be the (plausible) relation between gas metallicity and stellar metallicity, and in particular that lower metallicity gas is ionized by lower metallicity stars. For a given stellar mass, lower metallicity stars emit a harder and stronger radiation field, therefore giving a higher ionization parameter. The latter effect would naturally yield a U-Z relationship. The former are just qualitative interpretations. However, a thorough investigation of this phenomenon will requite detailed observational studies of stellar population in star forming galaxies.

The comparison of the empirical and the theoretical sequences of the two U-sensitive diagnostic flux ratios, F([O III]$\lambda $5007)/F([N II]$\lambda $6584) and F([O III]$\lambda $5007)/ F([O II]$\lambda $3727), also suggests the fact that the dispersion of the ionization parameter for a given metallicity should be relatively small. The typical RMS of the two flux ratios are $\sim$0.5 (in logarithm) at $12+\log($O/H $) \sim 7.5$. This corresponds to an RMS of the ionization parameter of $\sim$0.5 dex. This is the reason why the very U-sensitive flux ratio, F([O III]$\lambda $5007)/F([O II]$\lambda $3727), shows a clear metallicity dependence as seen in Fig. 13. The U-metallicity relationship is also important to understand the behavior of the empirical metallicity dependence of the flux ratio F([O III]$\lambda $5007)/F([N II]$\lambda $6584). This flux ratio is predicted to decrease with the oxygen abundance below $12+\log($O/H $) \sim 7.6$ by photoionization models with a constant ionization parameter. Owing to the metallicity dependence of the ionization parameter, this flux ratio does not show the "turnover'' seen in R23 and thus it is very useful to investigate the gas metallicity of galaxies without the measurement of F([O III]$\lambda $4363). Another implication of these results is that one should not use constant-U photoionization models to derive the oxygen abundance from the observed flux ratios, not only from F([O III]$\lambda $5007)/F([N II]$\lambda $6584 but also from any other metallicity diagnostics, which introduce systematic errors in the calibration.

Table 6: Coefficients of the best-fit polynomials for the observed relations between the emission-line flux ratios and the oxygen abundance, where $\log~R = b_0 + b_1 y + b_2 y^2 + b_3 y^3$ [ $y \equiv 12+\log($O/H)].

The empirical relations provided in this paper (Tables 5 and 6) are very useful to avoid such systematic errors to derive the gas metallicity by using only strong emission lines.

As for the U-insensitive diagnostic flux ratios, F([N II]$\lambda $6584)/F(H$\alpha $), F([N II]$\lambda $6584)/F([O II]$\lambda $3727) and F([N II]$\lambda $6584)/F([S II]$\lambda $6720), there are no significant discrepancies between the empirical sequence and the theoretical sequence (with a constant ionization parameter). This indirectly supports the above interpretation that the apparent discrepancy in R23 between the empirical sequence and the results of photoionization model is caused by the effect of the ionization parameter. Note that there is little or no metallicity dependence of the flux ratios of F([N II]$\lambda $6584)/F([O II]$\lambda $3727) and F([N II]$\lambda $6584)/F([S II]$\lambda $6720) in the low-metallicity range, $12+\log($O/H$) \la 8.0$, in terms both of empirical and theoretical dependences. Therefore these diagnostic flux ratios are useful only for the high metallicity galaxies.

The photoionization models presented in Figs. 13 and 14 suggest an additional interpretation of the discrepancy in some diagnostics between the two samples discussed in Sect. 4.1 (i.e., the discontinuity between sample A+B and sample C). Focusing on the metallicity range of $12+\log($O/H $)
\sim 8.3$ where the two datasets of sample A+B and sample C overlap, we note that the trend of the discrepancy suggests that the galaxies in sample A+B have higher ionization parameter than the galaxies in sample C. This supports the interpretation that the discrepancy is at least partly caused by the selection effect, i.e., galaxies with higher ionization parameter are selectively picked up in sample A+B. Then, what causes this selection effect? This may be related with the fact that the [O III]$\lambda $4363 emission is extremely weak in higher metallicity galaxies. This means that we can measure the [O III]$\lambda $4363 flux of galaxies with $12+\log($O/H $)
\sim 8.3$ (the highest metallicity in the galaxies in sample A+B) only when the [O III] emission is very strong, which corresponds to a very high ionization parameter.

4.4 Implications for studies of high-redshift galaxies and new diagnostics

Although the R23 method is thought to be a good metallicity diagnostic, various other diagnostics (some of which are investigated in this paper) have been proposed up to now. Indeed one of the main problems of the R23 method is that there are two solutions for a given R23 value and thus one cannot obtain a unique metallicity solution. Most of the newly proposed diagnostics use the [N II]$\lambda $6584 line to remove the degeneracy in R23, because the secondary nucleosynthesis of nitrogen makes this line emission very sensitive to the gas metallicity. However, there are two non-negligible problems with the use of the [N II]$\lambda $6584 line. First, especially for low-metallicity systems, the contribution of the primary nucleosynthesis and the secondary nucleosynthesis in the nitrogen abundance is not well understood, which leads to an uncertainty in the relative nitrogen abundance as a function of the metallicity. Second, the [N II]$\lambda $6584 emission is in the red part of the rest-frame optical spectrum of galaxies, which prevents its application to the observational investigations of high-z systems. For example, the optical detectors with a sensitivity up to $\lambda \sim 1~\mu$m can detect the [N II]$\lambda $6584 emission of galaxies only at $z \la 0.52$, and the K-band atmospheric window limits the highest redshift to $z \sim 2.7$ for ground-based facilities. Although one of the undoubtfully interesting targets for the JWST is the population related to the cosmic reionization, the sensitivity of NIRSpec (Posselt et al. 2004) boarded on JWST can examine the [N II]$\lambda $6584 emission of the objects at $z \la 6.6$, where the cosmic reionization has already nearly ended (e.g., Kashikawa et al. 2006; Fan et al. 2006). Another problem associated with the [N II]$\lambda $6584 line is that it becomes very weak and difficult to measure at low metallicities: [N II]$\lambda $6584/H $\alpha < 0.1$at $12+\log($O/H) < 8.5.

Our results on the empirical metallicity dependences suggest that one does not need [N II]$\lambda $6584 any more to distinguish the upper- and lower-branches of the R23 sequence. This is because the flux ratio of F([O III]$\lambda $5007)/F([O II]$\lambda $3727) is also a good metallicity diagnostics, thanks to the small dispersion of the ionization parameter at a given metallicity. The empirical R23 sequence peaks at $12+\log($O/H $)
\sim 8.0$, where the empirically determined flux ratio of F([O III]$\lambda $5007)/F([O II]$\lambda $3727) is $\sim$2. Therefore one can recognize whether the observed R23 belongs to the upper-branch of the R23 sequence or not, depending on whether F([O III]$\lambda $5007)/F([O II]$\lambda $3727) < 2 or not. Note that this result is consistent with an earlier remark by Maier et al. (2004) that the flux ratio of F([O III]$\lambda $5007)/F([O II]$\lambda $3727) can be used to distinguish the upper- and lower-branches of the R23 sequence. Our work gives the physical explanation for this idea (the U-Z relation) and a criterion to distinguish the degeneracy [F([O III]$\lambda $5007)/F([O II]$\lambda $3727) < 2] on the remark by Maier et al. (2004).

The above result is due to the fact that the ionization parameter has a strong metallicity dependence, and it thus implies that the ionization parameter itself is a sort of metallicity diagnostic. Motivated by this, we examine the metallicity dependence of the flux ratio F([Ne III]$\lambda $3869)/F([O II]$\lambda $3727), in Fig. 15. The reasons for focusing on this flux ratio are: (a) the two emission lines have different ionization degrees, their ratio should have a strong dependence on the ionization parameter and therefore is a possible good metallicity diagnostics; (b) their wavelength separation is very small and thus their flux ratio is not significantly affected by dust reddening; and (c) the two lines are located at a blue part in the rest-frame optical spectrum and thus their flux ratio could be a powerful diagnostic even for high-z galaxies. As expected, this flux ratio shows a clear metallicity dependence, which is apparently seen in Fig. 14. In Tables 7 and 8, the mean and the RMS of this flux ratio for within each bins of oxygen abundance are given, just similar to Tables 3 and 4 (Sect. 3). To obtain the analytic expression of this relation, we fit the observed sequence with a second-order polynomial function. The coefficients of the fit are given in Tables 9 and 10.


  \begin{figure}
\par\includegraphics[angle=-90,width=6.5cm,clip]{5216f15.eps}
\end{figure} Figure 15: ( Upper) Emission-line flux ratios of F([Ne III]$\lambda $3869)/ F([O II]$\lambda $3727) of the galaxies in sample A (red triangles), those in sample B (green circles), and those in sample C (black dots), as a function of the oxygen abundance. As for sample C, only the objects with S/N([Ne III]) > 10 and S/N([O II]) > 10 (cataloged values) are plotted. ( Lower) Same as the upper panel but the mean and the RMS values are shown for each bin of the oxygen abundance, instead of the individual data. Filled stars and filled squares denote the mean flux ratios for galaxies in sample A+B and those for galaxies in sample C, respectively. The errorbar denotes the RMS for each metallicity bin. The dashed line denotes the best-fit polynomial (second-order) function.
Open with DEXTER

This flux ratio can be measured for galaxies up to $z \sim 1.6$ with optical instruments, up to $z \sim 5.2$ with near-infrared instruments on the ground-based facilities, and up to $z \sim 12$ with JWST/NIRSpec, therefore this flux ratio is a promising tool for metallicity studies at high redshift. In particular, it is useful for low metallicity galaxies, for which the intensity of [Ne III]$\lambda $3869 becomes comparable to [O II]$\lambda $3727 and therefore easier to detect [F([Ne III]$\lambda $3869)/F([O II]$\lambda $3727) > 0.2 at $12+\log{\rm O/H} <8$]. Detailed theoretical calibrations on this flux ratio are required, taking the metallicity dependence of the ionization parameter into account, which go beyond the scope of this paper.

One possible caveat for the use of the diagnostic flux ratios of F([Ne III]$\lambda $3869)/F([O II]$\lambda $3727) [and F([O III]$\lambda $5007)/F([O II]$\lambda $3727), too] may be the effect of AGNs. Since AGNs also tend to show higher ratios of F([Ne III]$\lambda $3869)/F([O II]$\lambda $3727) and F([O III]$\lambda $5007)/ F([O II]$\lambda $3727), galaxies harboring an AGN may be misidentified as low-metallicity galaxies. However, we can identify AGNs through the detection of He II$\lambda $4686 and/or [Ne V]$\lambda $3426. Nagao et al. (2001) reported that typical type-2 AGNs show F([Ne V]$\lambda $3426)/F([O II]$\lambda $3727 $)
\sim 0.4$, and typical type-1 AGNs show even higher ratio ($\ga$1). Lamareille et al. (2004) also reported that AGNs and star-forming galaxies can be distinguished by using diagnostic diagrams using only the blue part of the spectrum, i.e., O32 versus R23 and F([O III]$\lambda $5007)/F(H$\beta $) versus F([O II]$\lambda $3727)/F(H$\beta $) (see also Rola et al. 1997). These suggest that we can easily distinguish AGNs from low-metallicity galaxies by using only diagnostics available in the blue part of the spectrum, even with moderate quality spectroscopic data. Another caveat for the use of some diagnostic flux ratios calibrated in this paper especially for high-z galaxies is that several of the empirical relations rely on the U-Z relation. It is not obvious that the U-Z relation found in the local galaxies also holds for high-z galaxies. If the U-Z relation is a consequence of the relation between gas and stellar metallicity, as discussed in the previous section, then the relation is not expected to evolve and should remain valid at any redshift. Instead, if the U-Z relation is a consequence of the mass-metallicity relation which evolves with redshift (Savaglio et al. 2005; Erb et al. 2006; see also Maier et al. 2004), then also the U-Z relation may evolve with redshift and may require a re-calibration of our empirical relations at high redshift. The latter case would be a serious problem for several studies at high redshift. Indeed, most of the gas metallicity diagnostics discussed in this paper, including the ones most widely used (e.g. R23), are significantly affected by the dependence on the ionization parameter.

Another difficulty to measure the gas metallicities of high-z galaxies is the faintness of targets, which sometimes prevents from measuring accurate emission-line fluxes. The use of low-resolution grating to improve the signal-to-noise ratio may yield to a blending of the H$\alpha $ and [N II] emission lines, which results in poor determinations of the gas metallicity. Therefore It may be useful to investigate metallicity diagnostics which use the sum of F(H$\alpha $) and F([N II]). In particular, we have examined the metallicity dependence of the flux ratio F(H$\alpha $+[N II] $\lambda \lambda $6548,6584)/F([S II]$\lambda $6720) in Fig. 16. This group of lines can be measured even in low-resolution spectra and even in spectra covering a relatively narrow wavelength range, and therefore may be particularly useful in high-z studies. There is a clear dependence of this flux ratio on the oxygen abundance, seen as a $\cup$-shaped distribution with a minimum at $12+\log($O/H $) \sim 8.7$ (i.e., $Z_{\rm gas}
\sim Z_\odot$). The mean and the RMS of this flux ratio for each bin of oxygen abundance are given in Tables 7 and 8, and the coefficients of the fit are given in Tables 9 and 10. The observed distribution of this flux ratio is naturally expected, since the behavior of the nitrogen emission as a secondary element should dominate at the super-solar metallicity range, while F([N II]) and F([S II]) should become weak with respect to F(H$\alpha $) at the low-metallicities due to the decrease of the corresponding ions. Although this diagnostic like R23 has two solutions when the ratio is below 10, this ratio seems useful for low-metallicity galaxies where it is larger than 10, in which case it is possible to state that the object belongs to the lower branch of the $\cup$-shaped distribution. This diagnostic is also useful when the [S II] emission is not detected (this is frequently the case when high-zfaint galaxies are concerned). In this case, we can calculate a lower limit for this flux ratio, and we can derive accordingly an upper limit to the gas metallicity if the lower limit is larger than 10. Note that this diagnostic is essentially independent of dust reddening.

Table 7: Means and RMSs of additional emission-line flux ratios of the galaxies in the sample A+Ba.

Table 8: Means and RMSs of additional emission-line flux ratios of the galaxies in sample Ca.

Table 9: Coefficients of the best-fit polynomials for the observed relations between the additional emission-line flux ratios and the oxygen abundance, where $\log ~R = a_0 + a_1 x + a_2 x^2 + a_3 x^3$ [ $x \equiv \log(Z/Z_\odot) \equiv 12+\log($O/H)-8.69].

Table 10: Coefficients of the best-fit polynomials for the observed relations between the additional emission-line flux ratios and the oxygen abundance, where $\log~R = b_0 + b_1 y + b_2 y^2 + b_3 y^3$ [ $y \equiv 12+\log($O/H)].


  \begin{figure}
\par\includegraphics[angle=-90,width=6.5cm,clip]{5216f16.eps}
\end{figure} Figure 16: Same as Fig. 15 but for the emission-line flux ratio of F(H$\alpha $+[N II] $\lambda \lambda $6548,6584)/F([S II]$\lambda $6720).
Open with DEXTER


  \begin{figure}
\par\includegraphics[angle=-90,width=6.5cm,clip]{5216f17.eps}
\end{figure} Figure 17: Same as Fig. 15 but for the emission-line flux ratio of F([O III]$\lambda $5007)/F(H$\beta $).
Open with DEXTER


  \begin{figure}
\par\includegraphics[angle=-90,width=6.5cm,clip]{5216f18.eps}
\end{figure} Figure 18: Same as Fig. 15 but for the emission-line flux ratio of F([O II]$\lambda $3727)/F(H$\beta $).
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=8cm,clip]{5216f19.eps}
\end{figure} Figure 19: Schematic view of the availability of various metallicity diagnostics for each redshift. The black solid curves indicate the effect of redshift for some of the diagnostic lines discussed in this paper. The colored boxes indicate the wavelength coverages of optical spectrometers (blue), of ground-based near-IR spectrometers (magenta), and of NIRSpec/MIRI on board of JWST (red). The marks on the right of the diagram indicate the maximum redshift at which some of the metallicity diagnostics can be used with the various facilities.
Open with DEXTER

Finally, We have also investigated the flux ratios of F([O III]$\lambda $5007)/F(H$\beta $) and F([O II]$\lambda $3727)/F(H$\beta $) (Figs. 17 and 18). The empirical calibrations for these two flux ratios may be useful when either [O II] or [O III] are not available, because out of the wavelength range or on a strong OH airglow emission line, or in a region of bad atmospheric transmission. The means and the RMSs of these two flux ratios for each bin of oxygen abundance are given in Tables 7 and 8, and the coefficients of the fit are given in Tables 9 and 10. As expected, both of the two flux ratios again show $\cap$-shaped distributions just similar to the R23 parameter. The F([O III]$\lambda $5007)/F(H$\beta $) ratio is useful especially for high metallicity galaxies, because the targets should belong to the upper-branch of the distribution when F([O III]$\lambda $5007)/F(H $\beta) < 1$(although this might be wrong when extremely metal-poor galaxies [$12+\log($O/H) < 7.0] are concerned). Note that this flux ratio is very sensitive to the oxygen abundance and the dispersion of the data is small at F([O III]$\lambda $5007)/F(H $\beta) < 1$.

Figure 19 summarizes the use of some of the metallicity diagnostics discussed in this paper as a function of redshift and for various observing facilities, and in particular optical spectrometers, ground-based near-IR spectrometers and NIRSpec on board of JWST. In principle (i.e., sensitivity permitting), MIRI on board of JWST will be able to observe the same diagnostics at even higher redshifts. Note that the ratio F([Ne III]$\lambda $3869)/F([O II]$\lambda $3727) extends the diagnostic capability of any observing facility to significantly higher redshift.

5 Summary

We have combined two large spectroscopic datasets to derive empirical calibrations for gas metallicity diagnostics involving strong emission lines. The two datasets consist of about 50 000 spectra from the SDSS DR4, which probe metallicities $12+\log($O/H)>8.3 (sample C), and of 328 spectra of low metallicity galaxies with a measurement of the [O III]$\lambda $4363 line (sample A+B), which probe metallicities $12+\log($O/H)<8.4. Together, these two samples provide the largest dataset of galaxies with known metallicity currently available, and spanning more than 2 dex in metallicity.

We have provided empirical calibrations both for metallicity diagnostics already proposed in the past and for new metallicity indicators proposed in this paper. We have given an analytical description for the metallicity dependence of the following diagnostics and line ratios: R23, F([N II]$\lambda $6584)/F(H$\alpha $), F([O III]$\lambda $5007)/F([N II]$\lambda $6584), F([N II]$\lambda $6584)/F([O II]$\lambda $3727), F([N II]$\lambda $6584)/ F([S II]$\lambda $6720), F([O III]$\lambda $5007)/F([O II]$\lambda $3727), and F([Ne III]$\lambda $3869)/F([O II]$\lambda $3727). The calibrations are performed within the metallicity range $7.0 \leq 12+\log($O/H $) \leq 9.2$. All of the investigated flux ratios show strong dependences on metallicity, at least in some metallicity ranges. We have shown that the monotonic metallicity dependence of the ratio F([O III]$\lambda $5007)/F([O II]$\lambda $3727) can be used to break the degeneracy of the R23 parameter when F([N II]$\lambda $6584)/F(H$\alpha $) is not available. The F([O III]$\lambda $5007)/F([O II]$\lambda $3727) ratio is particularly useful at high redshift, where H$\alpha $ and [N II]$\lambda $6584 are shifted outside the observed band. Another promising metallicity tracer at high-z is the ratio F([Ne III]$\lambda $3869)/F([O II]$\lambda $3727), which is found to anti-correlate with metallicity. The F([Ne III]$\lambda $3869)/F([O II]$\lambda $3727), ratio is particularly useful at high redshift, where most of the other diagnostic lines are shifted outside the observed band.

We have also investigated the observed relationships through a comparison with photoionization models. Some of the diagnostics investigated in this paper are strongly dependent on the ionization parameter U. The observed trends of these diagnostics highlight a clear, inverse relationship between ionization parameter and metallicity in galaxies. Such a strong U-Z relationship is also required to explain the trend observed for the R23 parameter. The U-Z relationship is relatively tight and, indeed, we have found that at any given metallicity the ionization parameter has a small dispersion ($\sim$0.5 dex). The strong relationship between ionization parameter and metallicity in galaxies should warn about the use of simple models, which assume constant ionization parameter, to infer gas metallicities from line ratios.

Acknowledgements
We thank J. Lee for comments on the flux data of the KISS galaxies, Y. I. Izotov for comments on the relation between $T_{\rm e}$(O+) and $T_{\rm e}$(O2+), L. Kewley for providing us their model results, and C. Maier, M. Onodera, and the anonymous referee for useful comments on this work. This work is based on the SDSS data catalogs released from Max Planck Institute for Astrophysics (MPA) and John Hopkins University (JHU), and produced by S. Charlot, G. Kauffmann, S. White, T. Heckman, C. Tremonti, and J. Brinchmann. T.N. acknowledges financial support from the Japan Society for the Promotion of Science (JSPS) through JSPS Research Fellowship for Young Scientists. R.M. and A.M. acknowledge financial support from the Italian Space Agency (ASI) and the Italian Institute for Astrophysics (INAF).

References

 

  
Online Material

Table 1: Re-calculated properties of the compiled low-metallicity galaxies.



Copyright ESO 2006