Issue |
A&A
Volume 495, Number 1, February III 2009
|
|
---|---|---|
Page(s) | 73 - 81 | |
Section | Extragalactic astronomy | |
DOI | https://doi.org/10.1051/0004-6361:200810558 | |
Published online | 15 October 2008 |
Physical properties of galaxies and their evolution in the VIMOS VLT Deep Survey![[*]](/icons/foot_motif.gif)
II. Extending the mass-metallicity relation to the range z
0.89-1.24
E. Pérez-Montero1,
- T. Contini1 - F. Lamareille1 - J. Brinchmann19 - C. J. Walcher7 - S. Charlot8,10 - M. Bolzonella3
- L. Pozzetti3
- D. Bottini2 - B. Garilli2 - V. Le Brun7 - O. Le Fèvre7 - D. Maccagni2 - R. Scaramella4,13 - M. Scodeggio2 - L. Tresse7 - G. Vettolani4 - A. Zanichelli4 - C. Adami7 - S. Arnouts23,7 - S. Bardelli3 - A. Cappi3 - P. Ciliegi3
- S. Foucaud21 - P. Franzetti2 - I. Gavignaud12 - L. Guzzo9 - O. Ilbert20 - A. Iovino9 - H. J. McCracken10,11 - B. Marano6
- C. Marinoni18 - A. Mazure7 - B. Meneux2,9 - R. Merighi3
- S. Paltani15,16 - R. Pellò1 - A. Pollo7,17 - M. Radovich5 - D. Vergani2 - G. Zamorani3
- E. Zucca3
1 - Laboratoire d'Astrophysique de Toulouse-Tarbes, Université de Toulouse, CNRS, Observatoire Midi-Pyrénées, 14 Av. E. Belin, 31400, Toulouse, France
2 -
IASF-INAF, via Bassini 15, 20133, Milano, Italy
3 -
INAF-Osservatorio Astronomico di Bologna, via Ranzani 1, 40127, Bologna, Italy
4 -
IRA-INAF, via Gobetti 101, 40129, Bologna, Italy
5 -
INAF-Osservatorio Astronomico di Capodimonte, via Moiariello 16, 80131, Napoli, Italy
6 -
Università di Bologna, Dipartimento di Astronomia, via Ranzani 1, 40127, Bologna, Italy
7 -
Laboratoire d'Astrophysique de Marseille, UMR 6110
CNRS-Université de Provence, 38 rue Frederic Joliot-Curie, 13388
Marseille Cedex 13, France
8 -
Max Planck Institut für Astrophysik, 85741, Garching, Germany
9 -
INAF-Osservatorio Astronomico di Brera, via Brera 28, 20021, Milan, Italy
10 -
Institut d'Astrophysique de Paris, UMR 7095, 98 bis Bvd Arago, 75014, Paris, France
11 -
Observatoire de Paris, LERMA, 61 Avenue de l'Observatoire, 75014, Paris, France
12 -
Astrophysical Institute Potsdam, An der Sternwarte 16, 14482, Potsdam, Germany
13 -
INAF-Osservatorio Astronomico di Roma, via di Frascati 33, 00040, Monte Porzio Catone, Italy
14 -
Universitá di Milano-Bicocca, Dipartimento di Fisica, Piazza delle Scienze 3, 20126, Milano, Italy
15 -
Integral Science Data Centre, ch. d'Écogia 16, 1290, Versoix, Switzerland
16 -
Geneva Observatory, ch. des Maillettes 51, 1290, Sauverny, Switzerland
17 -
Astronomical Observatory of the Jagiellonian University, ul Orla 171, PL-30-244, Kraków, Poland
18 -
Centre de Physique Théorique, UMR 6207 CNRS-Université de Provence, 13288, Marseille, France
19 -
Centro de Astrofísica da Universidade do Porto, Rua das Estrelas, 4150-762, Porto, Portugal
20 -
Institute for Astronomy, 2680 Woodlawn Dr., University of Hawaii, Honolulu, Hawaii, 96822, USA
21 -
School of Physics & Astronomy, University of Nottingham, University Park, Nottingham, NG72RD, UK
22 -
Max Planck Institut für Extraterrestrische Physik (MPE), Giessenbachstrasse 1,
85748 Garching bei München, Germany
23 -
Canada France Hawaii Telescope corporation, Mamalahoa Hwy,
Kamuela, 96743, USA
Received 10 July 2008 / Accepted 3 October 2008
Abstract
Aims. We present a continuation of our study about the relation between stellar mass and gas-phase metallicity in the VIMOS VLT Deep Survey (VVDS). In this work we extend the determination of metallicities up to redshift 1.24 for a sample of 42 star-forming galaxies with a mean redshift value of 0.99.
Methods. For a selected sample of emission-line galaxies, we use both diagnostic diagrams and empirical calibrations based on [O II] emission lines along with the empirical relation between the intensities of the [O III] and [Ne III] emission lines and the theoretical ratios between Balmer recombination emission lines to identify star-forming galaxies and to derive their metallicities. We derive stellar masses by fitting the whole spectral energy distribution with a set of stellar population synthesis models.
Results. These new methods allow us to extend the mass-metallicity relation to higher redshift. We show that the metallicity determinations are consistent with more established strong-line methods. Taken together this allows us to study the evolution of the mass-metallicity relation up to
with good control of systematic uncertainties. We find an evolution with redshift of the average metallicity of galaxies very similar to those reported in the literature: for a given stellar mass, galaxies at
have, on average, a metallicity
0.3 dex lower than galaxies in the local universe. However we do not see any significant metallicity evolution between redshifts
(Paper I) and
(this paper). We find also the same flattening of the mass-metallicity relation for the most massive galaxies as reported in Paper I at lower redshifts, but again no apparent evolution of the slope is seen between
and
.
Key words: galaxies: evolution - galaxies: fundamental parameters - galaxies: abundances - galaxies: starburst
1 Introduction
The formation and evolution of galaxies at different cosmological epochs are driven mainly by two processes: the star formation history and the metal enrichment. These processes, which are strongly connected in the evolution of individual galaxies, can be altered by the peculiar dynamics of the gas. The rate of star formation can be increased as a consequence of mergers or due to the inflow of gas. Regarding metallicity, it is sensitive to metal losses due to stellar winds, supernovae and active nucleus feedbacks. Therefore, the study of the dominant mechanisms of evolution of galaxies at different redshifts rely on the determination of metallicity and total stellar mass.
The correlation between stellar mass and metallicity (M-Z relation) was firstly found in the nearby universe by Lequeux et al. (1979) for dwarf irregular galaxies. Since them, numerous works have been performed in order to understand and quantify this relation in the local universe (among others
Skillman et al. 1989; Brodie & Huchra 1991; Zaritsky et al.
1994; Richer & McCall 1994;
Garnett et al. 1997; or Pilyugin & Ferrini 2000). This relation has its counterpart in
the luminosity-metallicity relation which is easier to derive when a lower number of spectrophotometric bands are
available. Besides, in the last years, this work
has been complemented with the
use of large surveys of galaxies which allow the derivation of this relation in samples
statistically more significant, (e.g. SDSS, Tremonti et al. 2004; 2dFGRS, Lamareille et al. 2004), that confirm the decrease of metal loss when galaxy stellar mass increases. Several studies of the evolution of this relation at larger redshifts have appeared due to the growing number of spectroscopic surveys of star-forming galaxies in the distant universe. This is the case of Savaglio et al. (2005) or Lamareille et al. (2006) who agree to find lower mean metallicities for the same galaxy mass bins at larger redshifts. However, they do not agree about the evolution of the slope of the M-Z relation: Lamareille et al. (2006) reported a flattening of the relation at higher redshifts whereas Savaglio et al. (2005) found a steepening of the M-Z relation.
Unfortunately, these studies are limited to z < 1 galaxies, due to the required presence of certain emission lines in the optical spectrum to derive the metallicity. Erb et al. (2006) studied this relation for a sample of star forming galaxies at
using near-infrared spectra and they find a lower value of the mean metallicity for a given galaxy stellar mass.
This evolution is apparently much more stronger at higher redshift as reported by Maiolino et al. (2008), using also IR spectra for a sample of star-forming galaxies ar z>3.
However, there is a lack of observations and available data in the redshift range 0.9 - 2, corresponding to the epoch of the universe when the maximum of star formation took place. It is therefore particularly interesting to study the
M-Z relation at this period.
The estimation of the global trend of the M-Z relation at different cosmological epochs is possible only through the analysis
of statistically significant samples with accurate determinations of the stellar mass and the gas phase oxygen abundance.
This paper is a continuation of Lamareille et al. (2008,
hereinafter Paper I), who study the M-Z relation of the star-forming galaxies detected in the VIMOS VLT Deep Survey up to redshift 0.89. In this work we extend the upper redshift limit of the study up to
,
using the relation between metallicity and emission-line strength of [O II] and [Ne III] presented in Pérez-Montero et al. (2007). In the next section we describe the sample of star-forming galaxies discussing the selection criteria and the methods to derive their stellar mass and metal content, which can go
to higher redshift than other classical methods. In Sect. 3, we describe the
M-Z relation in the redshift bin
0.89 < z <1.24, with an average value in our sample of z = 0.99 and we compare our results with those found in Paper I. Finally, we summarise our results and present our conclusions in the last section. Throughout this paper we normalize
the derived stellar masses and the absolute magnitudes with the standard
-CDM cosmology, i.e.,
h = 0.7,
and
(Spergel et al. 2007).
2 Description of the sample
The VIMOS/VLT Deep Survey (VVDS, Le Fèvre et al. 2005) is one of the widest and deepest spectrophotometric surveys of distant galaxies with a mean redshift of
.
The optical spectroscopic data belong to the first epoch sample obtained with the VIsible MultiObject Spectrograph (VIMOS) at ESO/VLT (UT3). They are divided in two deep fields (
): VVDS-02h (hereafter F02; Le Fèvre et al. 2005) and CDFS
(Le Fèvre et al. 2004) and in three wide fields (
): VVDS-10h, VVDS-14h and VVDS-22h (hereafter F10, F14 and F22 respectively, Garilli et al., in prep.).
As explained in Paper I, we divide therefore our
total sample in two: a deep one with the F02 and CDFS fields and a wide one with F10, F14 and F22 fields.
Spectroscopic data have been reduced using the VIPGI pipeline
(Scodeggio et al. 2005).
The following selection criteria have been applied on this sample. We consider only objects with a redshift known at a 75% confidence level (i.e. VVDS redshift flags 2, 3, 4 and 9) and we remove all stars, broad-line Active Galactic Nuclei (hereafter AGNs, Gavignaud et al. 2006) and duplicate observations.
The optical spectra of each galaxy have been split into two components (stellar and nebular) using the pipeline platefit_vimos as described in Paper I. This software removes the stellar continuum and absorption lines and then fits automatically all emission lines. The templates used to fit the stellar component have been derived from the libraries of Bruzual & Charlot (2003) resampled to the velocity dispersion of VVDS spectra. The fitting and measurement procedures of all emission lines (flux and equivalent width) and stellar indices in the original spectra are described in Paper I. In Fig. 1 we show an example of the fitting of the stellar continuum (top panel) made by platefit_vimos for the galaxy VVDS221485833, one of the best S/N spectra. The resultant nebular emission-lines spectrum, once this fitting is substracted, is shown in the bottom panel.
![]() |
Figure 1: VIMOS rest-frame spectrum of the galaxy VVDS220485833. The top panel illustrates the fitting of the stellar continuum obtained with the platefit_vimos tool. In the bottom panel, we show the resultant nebular emission-lines spectrum, once the fitted stellar continuum has been substracted. The labels indicate the set of emission lines used in this work to derive metallicities. |
Open with DEXTER |
2.1 Selection of star-forming galaxies
In order to correctly study the M-Z relation, we select the galaxies whose main mechanism of ionization is star formation.
We use emission-line diagnostic diagrams to identify narrow-line AGNs (Seyfert 2 and LINERS) once broad line AGNs have been
removed. In particular we use the ``blue'' diagnostic diagram based on [O III]5007 Å,
H
and [O II]3727 Å emission lines (Lamareille et al. 2004) and applied in Paper I for the redshift bin
0.5 < z < 0.9.
We extend this analysis up to redshift
1.3 using the theoretical ratios between H
and the most intense
Balmer hydrogen recombination lines at bluer wavelengths. This is the case of H
at 4340 Å and H
at 4102 Å.
We assume the following values for these ratios, typical of star-forming galaxies:
H
/H
and H
/H
.
![]() |
Figure 2: Histogram of the distribution of the ratio between the deredenned fluxes of [O III]4959+5007 Å and [Ne III]3869 Å emission lines for the five VVDS fields (VVDS-02h, VVDS-10h, VVDS-14h, VVDS-22h and CDFS). The mean value of this distribution, shown as a black solid line, is 15.23. The red solid line represents the mean value of the distribution without applying the mean reddening correction and the dashed line is mean value (=15.37) of this emission-line ratio for a representative sample of nearby galaxies studied in Pérez-Montero et al. (2007). |
Open with DEXTER |
Although these ratios slightly depend on the electron temperature of the ionized gas, this dependence is negligible in comparison with the intrinsic uncertainty of the measurement of the involved emission lines, even in the case of narrow-line AGNs, which show the lowest values of these ratios. We apply the empirical relation found between the emission-line ratio of [O III]4959+5007 Å and [Ne III]3869 Å (Pérez-Montero et al. 2007). The mean value of this ratio reported in this study for a sample of HII regions of the local universe is equal to 15.37. However, the mean value of the same ratio for the VVDS star-forming galaxies (both for the deep and wide fields) is equal to 17.37. We correct this slight difference taking into account a reddening correction corresponding to an extinction of AV = 0.5 mag. After this correction, the mean value of [O III]4959+5007 Å/[Ne III]3869 Å emission-line ratio for the VVDS galaxies is 15.23, which is almost identical to the value found for the reddening corrected HII regions of the local universe. Therefore, we decide to use the same reddening correction for the rest of the calculations in the VVDS sample. In Fig. 2, we show the distribution of the emission-line ratio [O III]4959+5007 Å/[Ne III]3869 Å for the 379 VVDS secure and candidates star-forming galaxies for which we have measured these lines in the redshift range 0.5 < z < 0.9 and with an error on the ratio lower than the ratio itself. In this figure we compare the mean value of this distribution, once extinction corrected (black solid line), the mean value found without extinction correction (red solid line) and the value found in the local universe by Pérez-Montero et al. (2007, dashed line).
![]() |
Figure 3: Spectral classification of VVDS galaxies using the diagram diagnostic based on [Ne III] emission line. In the left panel we show the diagram for 418 emission-line galaxies in the redshift range 0.5 < z < 0.9, while the right panel shows 72 emission-line galaxies in the redshift range 0.9 < z < 1.24. In these two panels, the red solid curve divides the regions of star-forming galaxies and narrow-line AGNs. The straight red solid line divides the AGN region between Seyfert 2 ( top) and LINER ( bottom). |
Open with DEXTER |
Table 1:
Distribution of star-forming galaxies and narrow-line AGNs among VVDS emission-line galaxies using the ``blue'' diagnostic diagrams
based on the relation between [O III] (blue,
0.5 < z < 0.9) or [Ne III] (blue-Ne,
0.5 < z < 1.24) and [O II] emission lines.
We show the distribution for the whole sample and for wide and deep fields. In the wide classification we include those galaxies belonging
to the deep fields having
.
![]() |
Figure 4: Stellar masses of star-forming galaxies (in solar masses) in the VVDS deep ( left panel) and wide ( right panel) fields as a function of redshift. The squares represent the star-forming galaxies selected using the diagnostic diagram based on [Ne III] emission line. The solid curve represents a logarythmic fit to the 50% mass-to-light completude levels and shows the bias imposed by the VVDS selection on apparent magnitude. |
Open with DEXTER |
The empirical limit between star-forming galaxies and narrow-line AGNs in the ``blue-Ne'' diagnostic diagram based on [Ne III]3869 Å
emission line takes the following form:
![]() |
(1) |
Following the same prescriptions as in the ``blue'' diagnostic method, the objects lying under this line are considered as star-forming galaxies, while those lying above it are considered as AGNs. This diagnostic diagram is, in fact, equivalent to that proposed by Rola et al. (1997) which is based on the same emission lines. In the AGN zone, LINERs are identified by the following criterium:
![]() |
(2) |
Adopting the emission-line ratios H






In Fig. 3 we show the [O II]/H
vs. [Ne III]/H
relation
for the sample of VVDS galaxies with the appropriate lines, in the redshift ranges
0.5 < z < 0.9(left panel) and
0.9 < z < 1.24 (right panel).
The dashed lines delineate the uncertainty domain (
0.15 dex around the empirical limit) in the spectral classification: objects lying in this region are considered as ``candidates'' to star-forming galaxies.
The results of this analysis are summarised in Table 1. We make also a comparison with the classification
derived from the ``blue'' diagnostic diagram based on [O II] and
[O III] emission lines for galaxies
in the redshift range
0.5 < z < 0.9.
We see that the fraction of galaxies classified as candidates to star-forming from the diagram based on
[Ne III] is much higher (at all redshift bins) that the one based on [O III] emission line. Taking into account
that the rate of coincidence between blue and blue-Ne diagnostics for the galaxies for which it is possible to apply at the
same time the two methods reaches 98%, we think that the sample of galaxies with accurate [Ne III] measurements (i.e. with a high enough S/N)
is somehow biased towards a population with enhanced [Ne III] emission.
This is indicated by the range of log ([Ne III]/H)
values
in both panels of Fig. 3 which are between -1 and 0, while in
Fig. 3 of Paper I we see that
the ratio log ([O III]/H
), which is approximately an order of magnitude larger
than log ([Ne III]/H
), varies between -0.5 and 1.
This means that in our sample we are certainly affected by some selection effect that
preserves us to properly detect all the star forming galaxies with a value of
log ([Ne III]/H
)
between -1.5 and -1. Therefore, the relative number of candidates is larger
in this subsample.
This same effect is seen when looking at the relative number of AGNs in the redshift bin 0.5 < z < 0.9, especially in the wide VVDS sample. On the other hand, when we compare the relative numbers of AGNs in the redshift regime z > 0.9, we do not see a significant increase compared with the blue diagnostic diagram at lower redshift. Note, finally, that we did not identify any LINER in our VVDS sample.
For the following analysis of stellar masses and metallicities, we consider both the samples of candidates and secure star-forming galaxies. Unfortunately, there are no deep enough X-ray or far-IR observations of the VVDS fields to check the presence of hidden AGNs.
2.2 Estimation of stellar masses
The procedure to derive stellar masses is based on the comparison between observations, basically the
photometric Spectral Energy Distribution (hereafter SED) of galaxies, combined in some cases with their
stellar absorption index H
and D4000 break measured in VIMOS spectra, and population synthesis models.
As described in Paper I and introduced by SDSS collaborators (Kauffmann et al. 2003; Tremonti et al. 2004; Brinchmann et al. 2005),
we are using the Bayesian approch.
To summarise, the population synthesis models are based on Bruzual & Charlot (2003), with the addition of secondary bursts to
the standard declining exponential star formation history (stochastic library, Salim et al. 2005; Gallazzi et al. 2005).
The SED of the models are compared with the available photometry for VVDS galaxies:
CFH12k observations in
bands for F02, F22 and F10 fields, and in
I band for F14 fields
(McCracken et al. 2003;
Le Fèvre et al. 2004), completed by U band observations in the
F02 field
(Radovich et al. 2004), CFHTLS
(Canada-France-Hawaii Telescope Legacy Survey) in
u*g'r'i'z' bands
for F02 and F22 fields and JKs photometry available for a substantial fraction of the objects
in the F02 field (Iovino et al. 2005;
Temporin et al. 2006). For the CDFS field, we use CFH12k
(Arnouts et al. 2001) and HST bviz (Giavalisco et al. 2004)
observations. Finally, objects in F02 and F22 fields have been cross-matched with the UKIDSS public
catolog (Warren et al. 2007), providing additional observations in
JK bands.
As in Paper I, we correct for the fact that both minimum and maximal stellar detected masses depend on redshift. We thus define a limiting mass for each redshift bin, which is the minimum mass that would be detected at a given redshift for a given limiting apparent magnitude. For our M-Z relation study we do not consider galaxies with a stellar mass lower that this minimum mass derived in the corresponding redshift bin.
The limiting mass for the redshift range
0.89 <z< 1.24 is derived empirically using the same procedure as in Paper I.
We derive the limiting mass for each galaxy following the expression:
![]() |
(3) |
where





2.3 Estimation of metallicities
We estimate metallicities as the oxygen abundances of the ionized gas phase from
the O
parameter (Pérez-Montero et al. 2007). This parameter is equivalent
to O23 (also known as R23), but it is based on the empirical relation existing
between [O III] and [Ne III] emission lines:
![]() |
(4) |
This calibration allows to derive metallicities using the theoretical ratios between H





![]() |
(5) |
We keep the factor 15.37 obtained in Pérez-Montero et al. (2007) for a sample of HII regions of the local universe once the observed fluxes have been dereddened. As shown in Fig. 2, we obtain a very similar factor to correct the flux of the involved emission lines, with a reddening equivalent to an extinction of AV=0.5 mag.
![]() |
Figure 5:
Comparison between O/H abundances derived for the star-forming galaxies of the SDSS DR4 catalogue
using the O
|
Open with DEXTER |




![]() |
Figure 6:
Relation between the O
|
Open with DEXTER |
![]() |
Figure 7:
Comparison between the oxygen abundances derived with the O23 and O
|
Open with DEXTER |
As described in Pérez-Montero et al. (2007), the
O
parameter has the
same problems than O23, the most important being its degeneracy with metallicity.
This means that for a given value of the O
parameter we can find two possible values of the metallicity.
Among the numerous calibrations of the O23 parameter there is always a calibration for high
metallicities (upper branch) and another one for low metallicities (lower branch) (e.g. Skillman 1989).
It is thus mandatory to be able to distinguish between these two metallicity regimes. The ratio [Ne III]/[O II]
has been proposed by Nagao et al. (2006) as an empirical calibrator of the metallicity but, due to its
high dependence on ionization parameter, it does not allow to retrieve accurate determinations of O/H.
However, it can be used to choose the appropriate branch (low or high metallicity regime) for the O
parameter.
The results obtained by Pérez-Montero et al. (2007) indicate that galaxies with log ([Ne III]/[O II]) < -1.0are in the high metallicity regime, although the contrary is uncertain. In
Fig. 6, we show the
relation between O
parameter and log ([Ne III]/[O II]) for the VVDS galaxies harbouring
the involved emission lines. We see that only a fraction of the objects lie
in the high metallicity regime using this criterion.
For the rest of galaxies with log ([Ne III]/[O II]) > -1.0, we constrain the metallicity regime by comparing the derived metallicities in low- and high-metallicity branches to the value derived with the Bayesian method described in Paper I or with the O23 parameter in the case of the redshift bin
0.5 < z < 0.9. After this comparison we conclude that 95% of the sample of star-forming galaxies lie in the high metallicity regime.
The selected objects have values of O
in the range
0.0-1.5, corresponding to oxygen abundances in the
range
in the upper branch.
In order to be consistent with metallicity estimates performed in Paper I, we use the calibrations of the upper branch of O23 parameter based on the models of Charlot & Longhetti (2001). As it has been shown by Kewley & Ellison (2008), it is very important to keep the consistency between the different empirical calibrators used to derive the metallicity of galaxies. If not, systematic differences as high as 0.7 dex can be observed for the M-Z relation in the local universe as derived from the SDSS.
We check the consistency of this new method at higher redshift by comparing the metallicities
derived using the O
and O23 parameters in the VVDS star-forming galaxies at z < 0.89.
Up to this redshift, all the emission lines required to compute these two parameters are observed
in the VIMOS optical spectra. For this purpose we take into account only the oxygen abundances derived
with an accuracy better than 0.3 dex. The results of this comparison are shown in Fig. 7.
We distinguish three cases where O
parameter is estimated using H
(268 galaxies), H
(217 galaxies) and H
(187 galaxies) respectively. Dashed lines represent the 1:1 relation while the
solid lines are the best linear fit to the data.
The best correlation is found in the first case, where the average of the residuals is 0.01 dex with a standard
deviation of 0.17 dex. For the second case (middle panel) the mean value of the residuals is also 0.01 dex, but the standard
deviation increases to 0.25 dex. Finally, in the third case (right panel), we obtain a mean value of the residuals of 0.04 dex
and a standard deviation of 0.24 dex. The same kind of analysis comparing the metallicities derived from the O
parameter
on one side and the Bayesian Method on the other side gives a mean value of the residuals of 0.03 dex only with a standard deviation
of 0.31 dex.
Table 2: Sample of 42 star-forming VVDS galaxies in the redshift regime 0.89 < z < 1.24 with a determination of both the stellar mass and the oxygen abundance. We report also their spectroscopic redshift and I-band (AB system) apparent magnitude.
![]() |
Figure 8:
Relation between the gas-phase metallicity and the stellar mass (in solar masses) for VVDS star-forming galaxies in the redshift range
0.89 < z < 1.24. In left panel, we show the 42 individual VVDS star-forming galaxies (black circles: deep fields; white circles: wide fields).
Vertical solid lines represent the corresponding limiting mass to the deep and wide samples.
The solid red line represents the Tremonti et al. (2004)
M-Z relation for the SDSS sample.
We show as well the typical uncertainty of the O
|
Open with DEXTER |
3 The mass-metallicity relation
We study the M-Z relation for the redshift bin
0.89 < z < 1.24, using the stellar masses and the metallicities
derived in previous sections. We take only galaxies with a stellar mass higher than the
limiting mass imposed in Fig. 4 for a mean redshift z=1, corresponding to a logarithmic stellar mass in solar masses
of 8.9 for the deep fields and 9.6 for the wide fields.
We take only galaxies with an uncertainty in the calculation of the oxygen abundance lower than 0.3 dex, so
they remain 42 star-forming VVDS galaxies: 27 in the deep fields and 15 in the wide fields.
The galaxies are listed in Table 2 along with their
apparent magnitude, spectroscopic redshift,
and derived metallicities and stellar masses. In the case of metallicity, all the errors due to
the uncertainties on the emission-line fluxes are lower than the intrinsic uncertainty of the empirical
parameter based on [Ne III] emission lines, which is 0.2 dex for the calibration of the upper branch
of the O
parameter (Pérez-Montero et al. 2007). The spectra of all these galaxies are available
via the CENCOS database
. We have to notice that, after applying the selection criteria
described above, no galaxies from the 10h and 14h VVDS fields remain in the final sample.
In left panel of Fig. 8 we show the relation between stellar mass and oxygen abundance for these 42 galaxies. For comparison, we show also the M-Z relation derived in the local universe (z < 0.3) by Tremonti et al. (2004) with SDSS data (solid line). Note that the comparison can be made directly as the same method/calibration (namely Charlot & Longhetti 2001, see Sect. 2.3) has been used to derive metallicities in the local and high-redshift galaxies. The respective limiting mass for both deep and wide fields are shown as vertical solid lines.
In order to test the global evolution in metallicity of star-forming galaxies we make the same analysis carried out
in Paper I, consisting in the calculation of the mean shift on metallicity by fitting the same curve deduced for
the SDSS sample to our data points and, hence, the derivation of a different zero-point.
These fits are shown in left panel of Fig. 8 as a blue dashed-point line for the galaxies in the deep fields and
a red dashed-point line for the wide field galaxies. The mean difference in metallicity found with respect to the local M-Z relation
are -0.29 0.21 and -0.12
0.15 for the wide and deep field galaxies respectively. As it was shown in Paper I, we see
a stronger evolution of the mean metallicity of galaxies in the wide sample. As the limiting
mass of the wide sample is much higher than in the deep sample, this could be due to a a more
significant evolution in metallicity for the most massive galaxies. The
M-Z relation derived in the redshift range
0.89 < z < 1.24 is rather similar to the one derived in Paper I for the redshift bin
0.7 < z < 0.9.
However, this difference between the wide and deep samples could be due to different selection effects. First the exposure time was higher for the deep fields which could lead to a higher dispersion in the metallicity distribution of these objects, as the number of high-metallicity galaxies, i.e. with fainter emission lines, is larger. There is also a larger number of low-mass galaxies in the deep sample, for which the metallicity is also lower compared with the average value.
We assess also the shape of the M-Z relation calculating a running median of the metallicity for the whole sample (both deep and
wide fields) along the stellar mass axis in order to compare it with the results from Paper I and the relation deduced
for the local universe using the SDSS sample (Tremonti et al. 2004).
The method to calculate these median is identical to that described in Paper I: we have taken four bins
in stellar mass, with a similar number of objects in each bin (11) and then we calculate the median
of 1000 bootstraps estimates.
These medians are shown in the right panel of Fig. 8 as black squares.
The error bars associated to the median metallicities of galaxies with z>0.9 are much higher than the ones derived
in Paper I at lower redshifts, due to the lower number of objects in each mass bin.
The M-Z relation found for
0.9<z<1.24 is fully compatible with the results of Paper I for z>0.5 in absolute value
and shape, pointing to a significant global chemical evolution for redshifts between 0 and 0.5.
However, for a redshift z>0.5, we do not find a significant global
evolution within the error bars up to z = 1.
Regarding the shape of the M-Z relation, we see a flattening as compared the with the M-Z relation in the local universe.
This could be a consequence of more effective yields for the more massive galaxies,
where the evolution is stronger. However, again, we do not appreciate any significant evolution of the slope of the
M-Z relation between redshifts 0.5 and 1.
It is expected somehow that the sample of star-forming galaxies considered in this work can be affected by some selection effects which lead to a higher fraction of objects with low metallicities. Firstly, the selection criterion on
[Ne III] emission lines biases the sample towards lower metallicities as pointed by Nagao et al. (2006) and Pérez-Montero et al. (2007). Secondly, the criterion on hydrogen recombination lines (H
and H
)
implies a sample of galaxies with higher star-formation efficiencies. According to Ellison et al. (2008), this leads also to galaxies with lower metallicities. Nevertheless, the consistency between the metallicities derived using the O
parameter in this redshift range and those at lower redshifts indicates that the probable selection effects inherent to our sample do not have a significant weight in the overall deduced metallicity.
4 Summary and conclusions
We have carried out a continuation of our study of the mass-metallicity relation of star-forming VVDS galaxies
extending our analysis up to redshift
.
We selected galaxies with a reliable measurement of [O II], [Ne III] and
H
or H
emission lines and we identified star-forming galaxies with the diagnostic diagram based
on these lines and presented in Pérez-Montero et al. (2007). Then, we derived oxygen gas-phase
abundances using the empirical parameter O
.
We have shown that these new methods
are consistent with those based on [O III] if we use the same calibrations.
As it is based on weaker emission lines that, generally, are more affected by absorption effects,
the only problem appears as a selection effect, biaising the sample towards galaxies with bright
emission lines and, hence, towards lower metallicities.
The relative number of narrow-line AGNs at higher redshifts is thus to be taken with caution.
The final sample used to derive the relation between metallicity and stellar mass consists of 42 VVDS star-forming galaxies with uncertainties lower than 0.3 dex in both oxygen abundance and stellar mass. These have been selected to be larger than the previously calculated limiting mass for this redshift range. Since the number of selected objects is sensibly lower than in our previous study (Paper I), the sample is not as statistically significant but we have found fully consistent results.
The zero-point and the slope of the M-Z relation in the 0.9 < z < 1.24 redshift range for both the deep and wide fields are rather similar to the values reported in Paper I for z>0.5, showing a significant global cosmic chemical evolution in relation with the local universe. We detect also a flattening of the M-Z relation in comparison with the SDSS-based relation in the local universe (Tremonti et al. 2004). This evolution is stronger for the galaxies of the wide sample, consistently with the fact that most massive galaxies have undergone a stronger chemical evolution. The evolution of massive galaxies compared with the local sample can be understood as a consequence of lower effective yields for the most massive galaxies at larger redshifts, according to an hierarchical scenario of formation and evolution of galaxies.
Nevertheless, despite the probable selection effects of the final sample,
we do not appreciate any significant evolution neither in the absolute value of the metallicity
nor in the shape of the M-Z relation for redshifts larger than 0.5. This can be a consequence of the smaller
time interval in the redshift range 0.7<z<1 combined with a larger uncertainty due to the lower number of
objects with an appropriate signal-to-noise ratio for the used emission-lines. In any case, our results
are fully compatible with those published for
in the DEEP2 sample (Shapley et al.
2005; Liu et al. 2008) or in the GOODS sample
(Cowie & Barger 2008),
who find 0.2-0.3 dex lower metallicities in relation with the local universe.
Acknowledgements
We thank C. Tremonti for having made the original platefit code available to the VVDS collaboration. We also thank the referee, Roberto Maiolino, for many valuable suggestions and comments, which have helped us to improve this paper. This work has been partially supported by the CNRS-INSU and its Programmes Nationaux de Galaxies et de Cosmologie (France) and by project AYA-2004-08260-C03-03 of the Spanish National Plan for Astronomy and Astrophysics. The VLT-VIMOS observations have been carried out on guarateed time (GTO) allocated by the European Southern Observatory (ESO) to the VIRMOS consortium, under a contractual agreement between the Centre National de la Recherche Scientifique of France, heading a consortium of French and Italian institutes, and ESO, to design, manufacture and test the VIMOS instrument. Based on observations obtained with MegaPrime/MegaCam, a joint project of CFHT and CEA/DAPNIA, at the Canada-France- Hawaii Telescope (CFHT) which is operated by the National Research Council (NRC) of Canada, the Institut National des Science de l'Univers of the Centre National de la Recherche Scientifique (CNRS) of France and the University of Hawaii. This work is based in part on data products produced at TERAPIX and the Canadian Astronomy Data Centre as part of the Canada-France-Hawaii Telescope Legacy Survey, a collaboration project of NRC and CNRS.
References
- Arnouts, S., Vandame, B., Benoist, C., et al. 2001, A&A, 379, 740 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Brinchmann, J., Charlot, S., White, S. D. M., et al. 2004, MNRAS, 351, 1151 [NASA ADS] [CrossRef]
- Brodie, J. P., & Huchra, J. P. 1991, ApJ, 379, 157 [NASA ADS] [CrossRef] (In the text)
- Bruzual, G., & Charlot, S. 2003, MNRAS, 344, 1000 [NASA ADS] [CrossRef] (In the text)
- Charlot, S., & Longhetti, M. 2001, MNRAS, 323, 887 [NASA ADS] [CrossRef] (In the text)
- Cowie, L. L., & Barger, A. J. 2008 [arXiv:0806.3457]
- Ellison, S. L., York, B. A., Murphy, M. T., et al. 2008, MNRAS, 383, 30 [NASA ADS]
- Erb, D. K., Shapley, A. E., Pettini, M., et al. 2006, ApJ, 644, 813 [NASA ADS] [CrossRef]
- Gallazzi, A., Charlot, S., Brinchmann, J., White, S. D. M., & Tremonti, C. A. 2005, MNRAS, 362, 41 [NASA ADS] [CrossRef] (In the text)
- Garnett, D. R., Shields, G. A., Skillman, E. D., Sagan, S. P., & Dufour, R. J. 1997, ApJ, 489, 63 [NASA ADS] [CrossRef] (In the text)
- Gavignaud, L., Bongiorno, A., Paltani, S., et al. 2006, A&A, 457, 79 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Giavalisco, M., Ferguson, H. C., Koekemoer, A. M., et al., 2004, ApJ, 600, L93 [NASA ADS] [CrossRef] (In the text)
- Iovino, A., McCracken, M. J., Garilli, B., et al. 2005, A&A, 442, 423 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Kauffmann, G., Heckman, T. M., Tremonti, C., et al. 2003, MNRAS, 341, 33 [NASA ADS] [CrossRef] (In the text)
- Kewley, L. J., & Ellison, S. L. 2008, AJ, 801, accepted [arXiv:0801.1849] (In the text)
- Kobulnicky, H. A., & Phillips, A. C. 2003, ApJ, 599, 1006 [NASA ADS] [CrossRef] (In the text)
- Lamareille, F., Mouhcine, M., Contini, T., Lewis, I., & Maddox, S. 2004, MNRAS, 350, 396 [NASA ADS] [CrossRef] (In the text)
- Lamareille, F., Contini, T., Le Borgne, J.-F., et al. 2006, A&A, 417, 839 (In the text)
- Lamareille, F., Brinchmann, J., Contini, et al. 2009, A&A, 495, 53 [CrossRef] [EDP Sciences] (Paper I)
- Le Fèvre, O., Meilleur, Y., McCracken, M. J., et al. 2004, A&A, 417, 839 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Le Fèvre, O., Vettolani, G., Garilli, B., et al. 2005, A&A, 439, 845 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Le Fèvre, O., Vettolani, G., Garilli, B., et al. 2006, ApJ, 647, 790
- Lequeux, J., Peimbert, M., Rayo, J. F., Serrano, A., & Torres-Peimbert, S. 1979, A&A, 80, 155 [NASA ADS] (In the text)
- Liu, X., Shapley, A. E., Coil, A. L., Brinchmann, J., & Ma, C.-P. 2008, ApJ, 678, 758 [NASA ADS] [CrossRef] (In the text)
- Maiolino, R., Nagao, T., Grazian, A., et al. 2008, A&A, 488, 463 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- McCracken, H. J., Radovich, M., Bertin, E., et al. 2003, A&A, 410, 17 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Nagao, T., Maiolino, R., & Marconi, A., 2006, A&A, 459, 85 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Pérez-Montero, E., Hägele, G. F., Contini, T., & Díaz, A. I. 2007, MNRAS, 381, 125 [NASA ADS] [CrossRef] (In the text)
- Pilyugin, L. S., & Ferrini, F. 2000, A&A, 358, 72 [NASA ADS] (In the text)
- Radovich, M., Arnaboldi, M., Ripepi, V., et al. 2004, A&A, 417, 51 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Richer, M. G., & McCall, M. L. 1995, ApJ, 445, 642 [NASA ADS] [CrossRef]
- Rola, C. S., Terlevich, E., & Terlevich, R. J. 1997, MNRAS, 289, 419 [NASA ADS] (In the text)
- Salim, S., Charlot, S., Rich, R. M., et al. 2005, ApJ, 619, L39 [NASA ADS] [CrossRef] (In the text)
- Savaglio, S., Glazebrook, L., Le Borgne, D., et al. 2005, ApJ, 635, 260 [NASA ADS] [CrossRef] (In the text)
- Scodeggio, M., Franzetti, P., Garilli, B., et al. 2005, PASP, 117, 1284 [NASA ADS] [CrossRef] (In the text)
- Shapley, A. E., Coil, A. L., Ma, C.-P., & Bundy, K. 2005, ApJ, 635, 1006 [NASA ADS] [CrossRef] (In the text)
- Skillman, E. D. 1989, ApJ, 347, 883 [NASA ADS] [CrossRef] (In the text)
- Skillman, E. D., Kennicutt, R. C., & Hodge, P. W. 1989, ApJ, 347, 875 [NASA ADS] [CrossRef]
- Spergel, D. N., Verde, L., Peiris, H. V., et al. 2003, ApJS, 148, 175 [NASA ADS] [CrossRef]
- Temporin, S., Iovino, A., Bolzonella, M., et al. 2008, A&A, 482, 81 [NASA ADS] [CrossRef] [EDP Sciences]
- Tremonti, C. A., Heckman, T. M., Kauffmann, G., et al. 2004, ApJ, 613, 898 [NASA ADS] [CrossRef] (In the text)
- Warren, S. J., Cross, N. J. G., Dye, S., et al. 2007, (In the text)
- Zaritsky, D., Kennicutt, R. C., & Huchra, J. P. 1994, ApJ, 420, 87 [NASA ADS] [CrossRef] (In the text)
Footnotes
- ... Survey
- Based on data obtained with the European Southern Observatory Very Large Telescope, Paranal, Chile, program 070.A-9007, and on data obtained at the Canada-France-Hawaii Telescope, operated by the CNRS in France, CNRC in Canada and the University of Hawaii.
- ...
- Postdoctoral fellow of the Spanish Fundación Española para la Ciencia y la Tecnología.
- ... database
- http://cencosw.oamp.fr
All Tables
Table 1:
Distribution of star-forming galaxies and narrow-line AGNs among VVDS emission-line galaxies using the ``blue'' diagnostic diagrams
based on the relation between [O III] (blue,
0.5 < z < 0.9) or [Ne III] (blue-Ne,
0.5 < z < 1.24) and [O II] emission lines.
We show the distribution for the whole sample and for wide and deep fields. In the wide classification we include those galaxies belonging
to the deep fields having
.
Table 2: Sample of 42 star-forming VVDS galaxies in the redshift regime 0.89 < z < 1.24 with a determination of both the stellar mass and the oxygen abundance. We report also their spectroscopic redshift and I-band (AB system) apparent magnitude.
All Figures
![]() |
Figure 1: VIMOS rest-frame spectrum of the galaxy VVDS220485833. The top panel illustrates the fitting of the stellar continuum obtained with the platefit_vimos tool. In the bottom panel, we show the resultant nebular emission-lines spectrum, once the fitted stellar continuum has been substracted. The labels indicate the set of emission lines used in this work to derive metallicities. |
Open with DEXTER | |
In the text |
![]() |
Figure 2: Histogram of the distribution of the ratio between the deredenned fluxes of [O III]4959+5007 Å and [Ne III]3869 Å emission lines for the five VVDS fields (VVDS-02h, VVDS-10h, VVDS-14h, VVDS-22h and CDFS). The mean value of this distribution, shown as a black solid line, is 15.23. The red solid line represents the mean value of the distribution without applying the mean reddening correction and the dashed line is mean value (=15.37) of this emission-line ratio for a representative sample of nearby galaxies studied in Pérez-Montero et al. (2007). |
Open with DEXTER | |
In the text |
![]() |
Figure 3: Spectral classification of VVDS galaxies using the diagram diagnostic based on [Ne III] emission line. In the left panel we show the diagram for 418 emission-line galaxies in the redshift range 0.5 < z < 0.9, while the right panel shows 72 emission-line galaxies in the redshift range 0.9 < z < 1.24. In these two panels, the red solid curve divides the regions of star-forming galaxies and narrow-line AGNs. The straight red solid line divides the AGN region between Seyfert 2 ( top) and LINER ( bottom). |
Open with DEXTER | |
In the text |
![]() |
Figure 4: Stellar masses of star-forming galaxies (in solar masses) in the VVDS deep ( left panel) and wide ( right panel) fields as a function of redshift. The squares represent the star-forming galaxies selected using the diagnostic diagram based on [Ne III] emission line. The solid curve represents a logarythmic fit to the 50% mass-to-light completude levels and shows the bias imposed by the VVDS selection on apparent magnitude. |
Open with DEXTER | |
In the text |
![]() |
Figure 5:
Comparison between O/H abundances derived for the star-forming galaxies of the SDSS DR4 catalogue
using the O
|
Open with DEXTER | |
In the text |
![]() |
Figure 6:
Relation between the O
|
Open with DEXTER | |
In the text |
![]() |
Figure 7:
Comparison between the oxygen abundances derived with the O23 and O
|
Open with DEXTER | |
In the text |
![]() |
Figure 8:
Relation between the gas-phase metallicity and the stellar mass (in solar masses) for VVDS star-forming galaxies in the redshift range
0.89 < z < 1.24. In left panel, we show the 42 individual VVDS star-forming galaxies (black circles: deep fields; white circles: wide fields).
Vertical solid lines represent the corresponding limiting mass to the deep and wide samples.
The solid red line represents the Tremonti et al. (2004)
M-Z relation for the SDSS sample.
We show as well the typical uncertainty of the O
|
Open with DEXTER | |
In the text |
Copyright ESO 2009
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.