A&A 459, 85-101 (2006)
DOI: 10.1051/0004-6361:20065216
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]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]
6584/H
and
[O III]
5007/[N II]
6584 ratios.
Our empirical diagrams show that the line ratio
[O III]
5007/[O II]
3727 is a
useful tool to break the degeneracy in the R23 parameter
when no information on the [N II]
6584 line is
available. The line ratio
[Ne III]
3869/[O II]
3727 also
results to be a useful metallicity indicator for high-z
galaxies, especially when the R23 parameter or other
diagnostics involving [O III]
5007 or
[N II]
6584 are not available.
Additional, useful diagnostics newly proposed in this
paper are the line ratios of
(H
+[N II]
6548,6584)/[S II]
6720,
[O III]
5007/H
,
and
[O II]
3727/H
.
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
0.7 dex
as the metallicity decreases from 2
to 0.05
,
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
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]4363 and [N II]
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
![]() |
(1) |
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]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
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 (
,
where
is the surface flux of hydrogen-ionizing
photon and
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
(
O/H
), whose metallicity is accurately
determined through the [O III]
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) (
O/H
;
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.
The gas-phase oxygen abundance is well determined when
the flux of [O III]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]
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
[
(log
]
(146 spectra).
Here we adopt the uncertainty
[
(log
)] 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
(log
is 139. Hereafter we call this sample "sample A''.
![]() |
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 O/H) < 7.6. Therefore the
statistical reliability
of the empirical calibration of metallicity diagnostics would
be extremely poor at
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
O/H) < 7.6are at z < 0.02. This means that the [O II]
3727
flux cannot be measured for the latter galaxies due to the limited
wavelength coverage of the SDSS spectroscopy
(
). Therefore, if using
only sample A, we could not calibrate the diagnostics involving
[O II]
3727 [R23 and F([N II]
6584)/F([O II]
3727)]
in the metallicity range
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]
4363-detected galaxies,
which are described in the following
subsection.
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]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 [
(S+)], gas
temperature of [O III]-emitting region
[
]
and oxygen abundance [
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]
4959) but calculate the ratio of
[F([O II]
3727
([O III]
5007)]
/F(H
)
since only the flux of
[O III]
5007 (without that of
[O III]
4959) is given in some reference papers.
![]() |
Figure 2:
Emission-line flux ratios of
[S II]![]() ![]() |
Open with DEXTER |
![]() |
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, [O/H
O/H
], 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
(log
.
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,
[
O/H
O/H)
], 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.
For high-metallicity galaxies, we referred to the oxygen
abundance derived by Tremonti et al. (2004), who derived the
metallicities of 53 000 galaxies in the SDSS database.
They used the fluxes of many strong emission lines
([O II]
3727, H
4861,
[O III]
5007, H
,
[N II]
6584,
[S II]
6717, and [S II]
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]
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:
Note that in sample C we do not put constraints on
(log
.
This is to prevent
sample C to be devoid of galaxies at metallicities
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.
![]() |
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 |
![]() |
Figure 5: Frequency distribution of the oxygen abundance of galaxies in sample C. |
Open with DEXTER |
![]() |
Figure 6:
Emission-line flux ratios of R23
(=[F([O II]![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 7:
Same as Fig. 6 but for the emission-line flux ratios of
F([N II]![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
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
(cataloged value) for all the related
emission lines (e.g., H
,
[O II]
3726,
[O II]
3729 and [O III]
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]
5007)/F([O II]
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:
![]() |
(2) |
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]6584)/F([S II]
6720)
versus the oxygen abundance shows an apparent
discontinuity between sample A+B and sample C, where
F([S II]
6720) denotes the sum of
F([S II]
6717) and
F([S II]
6731). We will discuss the issue
of this discontinuity in Sect. 4.1.
The diagram of R23 versus the oxygen abundance shows
a -shaped distribution with a peak at
O/H
.
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
O/H) < 8.55 with a bin width of
[log(O/H)] = 0.1, except at lowest and highest
oxygen abundances, where the bin width is wider (i.e.,
O/H) < 7.45 and
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
O/H) < 9.25 with a
bin width of
[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.
![]() |
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 [![]() |
Open with DEXTER |
![]() |
Figure 9:
Same as Fig. 8 but for the emission-line flux ratios of
F([N II]![]() ![]() ![]() ![]() ![]() ![]() |
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
O/H) < 9.25 (or
), 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
![]() |
(3) |
![]() |
(4) |
Table 5:
Coefficients of the
best-fit polynomials for the observed relations between
the emission-line flux ratios and the oxygen abundance,
where
[
log(Z/
O/H)-8.69].
Finally we recall that these relations are valid only
in the range
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.
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]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]
5007)/F([O II]
3727)
and
F([N II]
6584)/F([S II]
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]
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]4363
emission and thus the oxygen abundance may tend to be
underestimated accordingly. This is because the strength of
the [O III]
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
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]
6584)/F([S II]
6720),
it goes in the opposite direction to account for the
discrepancy seen in
F([O III]
5007)/F([N II]
6584)
and
F([O III]
5007)/ F([O II]
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)
.
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]
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]5007)/F([N II]
6584)
and
F([O III]
5007)/F([O II]
3727),
although the difference in the ionization parameter
should not cause a significant difference in the ratio of
F([N II]
6584)/F([S II]
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.
![]() |
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
[![]() |
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]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
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 [
O/H) > 8.5] is much larger
than the low metallicity sub-sample [
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
[
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
for all of the lines H
,
[O II]
3726,
[O II]
3729, and [O III]
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
,
H
and [N II]
6584,
not for [O II]
3726,
[O II]
3729, and [O III]
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
O/H) > 9, it is suggested that
our calibration for the R23 may overestimate the gas
metallicity at
O/H) > 9 by a factor of
dex at
O/H
.
The calibration of the diagnostic flux ratio
F([N II]6584)/F(H
)
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
(e.g., Erb et al. 2006).
In Fig. 11, we compare the empirical calibrations of
F([N II]
6584)/F(H
)
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
O/H) < 8.5.
However, the difference is significant (
dex)
only at metallicities
O/H) < 7.5 and
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]
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
O/H) = 8.5corresponds to
,
still in the sub-solar
metallicity domain.
![]() |
Figure 11:
Comparison of our results with the previous
empirical metallicity calibrations for
F([N II]![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 12:
Comparison of our results with the previous
empirical metallicity calibrations for
F([O III]![]() ![]() ![]() |
Open with DEXTER |
In Fig. 12, we compare the empirical calibrations of
F([O III]5007)/F([N II]
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,
O/H) < 8. However Pettini & Pagel (2004)
correctly mentioned that the flux ratio
F([O III]
5007)/F([N II]
6584
is of little use when
F([O III]
5007)/F([N II]
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]
6584 flux.
![]() |
Figure 13:
The averaged flux ratios and the best-fit polynomial
functions of the metallicity dependence of R23,
F([N II]![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 14:
Same as Fig. 13 but for the emission-line flux ratios of
F([N II]![]() ![]() ![]() ![]() ![]() ![]() |
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]6584)/F(H
)
that is
taken from Kobulnicky & Kewley (2004). Since the explicit
analytic expression for the metallicity dependence of
F([O III]
5007)/F([O II]
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
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
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
,
or equivalently,
(where
).
See Kewley & Dopita (2002) for details on the calculations.
Note that they adopted
O/H)
(Anders & Grevesse 1989) and expressed the metallicity in
units of
(
O/H)
). However,
since we adopt a more recent value for the solar abundance,
O/H)
(Allende Prieto et al. 2001),
the metallicity notation is different when the
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 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]
4363 line
(
dex), especially in the low metallicity
range (i.e.,
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
,
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 [
O/H
;
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]
5007)/F([N II]
6584)
(Fig. 13) and
F([O III]
5007)/F([O II]
3727)
(Fig. 14). By focusing on these two U-sensitive flux
ratios, we can see that the ionization parameter increases
by
0.7 dex with decreasing oxygen abundance from
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]5007)/F([N II]
6584)
and
F([O III]
5007)/F([O II]
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]
5007)/F([N II]
6584)
and
F([O III]
5007)/F([O II]
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]
5007)/F([O II]
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]5007)/F([O II]
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]5007)/F([N II]
6584)
and
F([O III]
5007)/ F([O II]
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
0.5 (in logarithm) at
O/H
.
This corresponds to an RMS of the ionization parameter of
0.5 dex. This is the reason why the very
U-sensitive flux ratio,
F([O III]
5007)/F([O II]
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]
5007)/F([N II]
6584).
This flux ratio is predicted to decrease with
the oxygen abundance below
O/H
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]
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]
5007)/F([N II]
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
[
O/H)].
As for the U-insensitive diagnostic flux ratios,
F([N II]6584)/F(H
),
F([N II]
6584)/F([O II]
3727)
and
F([N II]
6584)/F([S II]
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]
6584)/F([O II]
3727)
and
F([N II]
6584)/F([S II]
6720)
in the low-metallicity range,
O/H
,
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 O/H
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]
4363 emission is
extremely weak in higher metallicity galaxies. This means
that we can measure the [O III]
4363 flux of
galaxies with
O/H
(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.
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]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]
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]
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
m can detect the
[N II]
6584 emission of galaxies only at
,
and the K-band atmospheric window
limits the highest redshift to
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]
6584 emission of the objects at
,
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]
6584 line is that it becomes very weak and difficult to measure at
low metallicities: [N II]
6584/H
at
O/H) < 8.5.
Our results on the empirical metallicity dependences
suggest that one does not need [N II]6584
any more to distinguish the upper- and lower-branches
of the R23 sequence. This is because the flux
ratio of
F([O III]
5007)/F([O II]
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
O/H
,
where the
empirically determined flux ratio of
F([O III]
5007)/F([O II]
3727)
is
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]
5007)/F([O II]
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]
5007)/F([O II]
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]
5007)/F([O II]
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]3869)/F([O II]
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.
![]() |
Figure 15:
( Upper) Emission-line flux ratios of
F([Ne III]![]() ![]() |
Open with DEXTER |
This flux ratio can be measured for galaxies up to
with optical instruments, up to
with
near-infrared instruments on the ground-based facilities, and
up to
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]
3869
becomes comparable to [O II]
3727 and therefore
easier to detect
[F([Ne III]
3869)/F([O II]
3727)
> 0.2 at
].
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]3869)/F([O II]
3727)
[and F([O III]
5007)/F([O II]
3727),
too] may be the effect of AGNs. Since AGNs also tend to
show higher ratios of
F([Ne III]
3869)/F([O II]
3727)
and
F([O III]
5007)/ F([O II]
3727),
galaxies harboring an AGN may be misidentified as
low-metallicity galaxies. However, we can identify AGNs
through the detection of
He II
4686 and/or [Ne V]
3426.
Nagao et al. (2001) reported that typical type-2 AGNs show
F([Ne V]
3426)/F([O II]
3727
,
and typical type-1 AGNs show even higher ratio (
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]
5007)/F(H
)
versus
F([O II]
3727)/F(H
)
(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
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
)
and F([N II]). In
particular, we have examined the metallicity dependence of
the flux ratio
F(H
+[N II]
6548,6584)/F([S II]
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
-shaped distribution
with a minimum at
O/H
(i.e.,
). 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
)
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
-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
[
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
[
O/H)].
![]() |
Figure 16:
Same as Fig. 15 but for the emission-line flux ratio of
F(H![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 17:
Same as Fig. 15 but for the emission-line flux ratio of
F([O III]![]() ![]() |
Open with DEXTER |
![]() |
Figure 18:
Same as Fig. 15 but for the emission-line flux ratio of
F([O II]![]() ![]() |
Open with DEXTER |
![]() |
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 |
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]3869)/F([O II]
3727)
extends the diagnostic capability of any observing
facility to significantly higher redshift.
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 O/H)>8.3 (sample C), and of 328 spectra of low metallicity galaxies with a measurement of the
[O III]
4363 line (sample A+B), which probe
metallicities
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]6584)/F(H
),
F([O III]
5007)/F([N II]
6584),
F([N II]
6584)/F([O II]
3727),
F([N II]
6584)/ F([S II]
6720),
F([O III]
5007)/F([O II]
3727), and
F([Ne III]
3869)/F([O II]
3727).
The calibrations are performed within the metallicity
range
O/H
.
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]
5007)/F([O II]
3727)
can be used to break the degeneracy of the R23 parameter when F([N II]
6584)/F(H
)
is not available. The
F([O III]
5007)/F([O II]
3727)
ratio is particularly useful at high redshift, where
H
and [N II]
6584 are shifted outside
the observed band. Another promising metallicity tracer at
high-z is the ratio
F([Ne III]
3869)/F([O II]
3727),
which is found to anti-correlate with metallicity. The
F([Ne III]
3869)/F([O II]
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
(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(O+) and
(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).
Table 1: Re-calculated properties of the compiled low-metallicity galaxies.