A&A 423, 777-786 (2004)
DOI: 10.1051/0004-6361:20035913

Cu and Zn in the early Galaxy

G. Bihain1 - G. Israelian1 - R. Rebolo1,2 - P. Bonifacio3 - P. Molaro3


1 - Instituto de Astrofísica de Canarias, c/ Vía Láctea, s/n, 38205 La Laguna, Tenerife, Spain
2 - Consejo Superior de Investigaciones Científicas, Spain
3 - Istituto Nazionale di Astrofisica, Osservatorio Astronomico di Trieste, via G. B. Tiepolo 11, 34131 Trieste, Italy

Received 19 December 2003 / Accepted 27 April 2004

Abstract
We present Cu and Zn abundances for 38 FGK stars, mostly dwarfs, spanning a metallicity range between solar and [Fe/H] = -3. The abundances were obtained using Kurucz's local thermal equilibrium (LTE) model atmospheres and the near-UV lines of Cu I 3273.95 Å and Zn I 3302.58 Å observed at high spectral resolution. The trend of [Cu/Fe] versus [Fe/H] is almost solar for [Fe/H] > -1 and then decreases to a plateau $\langle$[Cu/Fe]$\rangle$ = -0.98 at [Fe/H] < -2.5, whereas the [Zn/Fe] trend is essentially solar for [Fe/H] > -2 and then slightly increases at lower metallicities to an average value of $\langle$[Zn/Fe]$\rangle$ = +0.18. We compare our results with previous work on these elements, and briefly discuss them in terms of nucleosynthesis processes. Predictions of halo chemical evolution fairly reproduce the trends, especially the [Cu/Fe] plateau at very low metallicities, but to a lesser extent the higher [Zn/Fe] ratios at low metallicities, indicating possibly missing yields.

Key words: stars: Population II - stars: abundances - Galaxy: evolution

1 Introduction

Massive, low- and intermediate-mass stars, type II and type Ia supernovae are likely contributors to the Galactic chemical evolution of Cu and Zn. However, the nucleosynthesis mechanisms and the relative contributions of the nucleosynthesis sites for these near- iron-peak nuclei are still uncertain. The abundances of Cu in metal-poor stars with [Fe/H] > -1.7 have been comprehensively reviewed by Peterson (1981), with no significant deviation from the solar ratio [Cu/Fe] = 0. Further studies by Luck & Bond (1985), Sneden & Crocker (1988) and Sneden et al. (1991) show that, while the [Zn/Fe] ratio remains essentially solar, there is in fact a clear decrease in [Cu/Fe] at lower metallicities, consistent with previous results for the metal-poor globular clusters M 13 and M 92 by Cohen (1978, 1979, 1980). This departure from the solar ratio is as noticeable as for the iron peak elements Cr, Mn and Co; [Co/Fe], however, is increasing at lower metallicities (e.g. Cayrel et al. 2003). More recent studies (Mishenina et al. 2002; Simmerer et al. 2003) show that the trend of Cu is non-linear in metal-poor stars and globular clusters, which presumably results from the superposition of Fe-independent ("primary'') and Fe-dependent ("secondary'') nucleosynthesis mechanisms. While there is a general consensus on the early contribution by short-lived massive stars, there are many uncertainties concerning how these stars produce Cu and Zn, and how their yields lead to the observed low-metallicity tails of the ratio distributions.

The above-mentioned studies are performed mostly with optical lines[*] (Cu I at 5105.54 Å, 5218.20 Å and 5782.12 Å; Zn I at 4722.16 Å, 4810.53 Å and 6362.35 Å). These lines become very weak in the spectra of very metal-poor stars. Therefore we decided to explore the use of the UV lines Cu I 3273.95 Å and Zn I 3302.58 Å, which are quite intense in the solar spectrum and, in the case of the Cu line, potentially detectable at very low metallicities ([Fe/H] < -3). In the following sections, we will discuss the observations, their analysis, and how these lines are indeed useful for investigating the chemical evolution of copper and zinc at low metallicities.


  \begin{figure}
\par\includegraphics[width=12cm,clip]{0913fig1.eps}\end{figure} Figure 1: Spectra of UV copper and zinc lines in metal-poor dwarfs.
Open with DEXTER

2 Observations

The studied sample contains 38 FGK metal-poor stars, principally dwarfs. Most of the high-resolution spectra were obtained at the 4.2 m WHT/UES (La Palma) during 1992-1993. Spectra of the stars HD 3795, LHS 540, HD 211998, HD 166913, HD 218502 and HD 128279 were obtained at the 3.6 m/CASPEC (La Silla) during 1993; and of the stars HD 22879 and HD 132475 at the 3.9 m AAT/UCLES (Siding Spring) during 1990-1991. The spectra of the very metal-poor stars LP 815-43, G275-4 and G64-12 were obtained at the 8.2 m VLT Kueyen/UVES (Paranal) during 1999-2000 (see details in Israelian et al. 2001). Finally, the spectra of the very metal-poor giants HD 2796 and BD -18 $\hbox {$^\circ $ }$5550 were obtained at the WHT/UES during 2000, as described in Israelian & Rebolo (2001). The échelle spectra covered the spectral range 3025-3825 Å, with spectral resolution $R \sim$ 40 000-60 000 and signal-to-noise ratio $S/N \sim$ 80-150 at 3300 Å. All the data were reduced using standard IRAF[*] routines for bias, scattered light, flat field correction and wavelength calibration. We present in Fig. 1 the spectra of five stars with different metallicities. The normalization of the spectra was achieved by fitting a low-order polynomial to continuum regions previously identified in the solar spectrum (Kurucz et al. 1984), in the neighborhood of the Cu and Zn lines.


  \begin{figure}
\par\includegraphics[width=14cm,clip]{0913fig2.ps}\end{figure} Figure 2: Expanded plots of the solar spectral regions ( dotted lines) fitted with the synthetic spectrum ( solid lines).
Open with DEXTER

3 Spectral synthesis and stellar parameters

Model atmospheres with "overshooting'' were chosen from the ATLAS9 grid (Kurucz 1993) by standard interpolation through [Fe/H], $T_{\rm eff}$ and log g. To analyze the blended lines of copper and zinc, 5 Å wide synthetic spectra were computed with the LTE code MOOG (Sneden 1973). The atomic parameters of the Cu I resonance lines at 3247.54 Å (log gf = -0.062, $\chi$ = 0.000) and at 3273.95 Å (log gf = -0.359, $\chi$ = 0.000), and those of the Zn I line at 3302.58 Å (log gf = -0.057, $\chi$ = 4.030) were taken from VALD-2 (Kupka et al. 1999). The wavelengths and the relative strengths of the hyperfine components of the copper isotopes 63Cu and 65Cu were obtained from Kurucz (2003, private communication), and the solar isotopic ratio was adopted from Anders & Grevesse (1989) for all the stars. The oscillator strengths of the lines were calibrated with the solar spectrum, using a model atmosphere with overshooting from Kurucz (1992, private communication) for the parameters $T_{\rm eff,\odot}$ = 5777 K, log $g_{\odot}$ = 4.438, $V_{\rm t,\odot}$ = 1,0 km s-1 and the Unsöld damping approximation without enhancement. The calibration provided log gf = -0.167 for the Zn I line (see Fig. 2). In the case of the Cu I lines log gf was not modified, because its modification strongly affected the suitable fit to the wings (as obtained from the cog driver of MOOG, the lines were in the linear damping part of the curve of growth) and did not improve the fit to the cores. In fact, the cores of such lines are usually formed close to the chromospheric temperature minimum, which exceeds that of the upper photosphere. Since no chromosphere is included in our model, these cores could not be reproduced as in the observed spectrum.

Stellar parameters for the model atmospheres were obtained as follows. The effective temperatures ( $T_{\rm eff}$) were estimated using the Alonso et al. (1996, 1999b) calibrations versus V-K colours for dwarfs and giants, respectively. The V magnitudes were taken from the Hipparcos Catalog (Perryman et al. 1997), Alonso et al. (1994) or else from the SIMBAD database. The K magnitudes were taken from Alonso et al. (1994, 1998), or else from Carney & Aaronson (1979; HD 170153, HD 165908, HD 225239 and HD 157214) and Laird et al. (1988; G170-47). For HD 166913, a $K_{\rm s}$ magnitude was taken from the 2MASS Catalogue, and then converted to K using a relation determined from a comparison of the magnitudes in these bands for similar stars in the sample. Since no accurate K or $K_{\rm s}$ values were available for the stars HD 211998 and HD 3795, the $T_{\rm eff}$-(B-V) calibrations from Alonso et al. (1996, 1999b) were used, with (B-V) from the Hipparcos Catalogue or the SIMBAD database. LTE metallicities were adopted from the literature, mainly from Fulbright (2000), Mishenina et al. (2002) and Thévenin (1998). The differences in the values provided by these three sources were typically not greater than 0.12 dex for the stars in common, therefore justifying the use of an average metallicity (see Table 1). Finally, the initial surface gravities were taken from the same references, with the exception of Thévenin & Idiart (1999) which was used instead of Thévenin (1998), because a higher priority was given to their NLTE log g values (see comments below). Stars with surface gravities log g > 3.5 were classified as dwarfs, while stars with log g < 2.5 were classified as giants, so that the corresponding $T_{\rm eff}$ colour calibration was used for each group. When stars had an intermediate log g value, $T_{\rm eff}$ was adopted as the average of the temperatures obtained with the calibrations for dwarfs and for giants, respectively.

We found a very good agreement between our $T_{\rm eff}$ values and the spectroscopic $T_{\rm eff}$ values from Mishenina & Kovtyukh (2001) for the 23 stars in common, $\langle$ $T_{\rm eff}-T_{\rm eff, MK}$$\rangle$ = +26 K, with a standard deviation $\sigma$ = 75 K (see Fig. 3). A comparison with the spectroscopic $T_{\rm eff}$ values from Fulbright (2000), however, showed a systematic difference of $\langle$ $T_{\rm eff}-T_{\rm eff, F}$$\rangle$ = +99 K with $\sigma$ = 82 K (25 stars). Since the 21 stars for this latter comparison were also included in the comparison with Mishenina & Kovtyukh's sample, this difference should be attributed to the different methods used. Fulbright determines an excitation temperature for the dwarfs (by requiring the same iron abundance for Fe I lines with high and low excitation potential), whereas Mishenina & Kovtyukh determine a temperature from the fits to the wings of the H$_{\alpha}$ lines. Comparing with Thévenin's (1998) $T_{\rm eff}$ values, we found $\langle$ $T_{\rm
eff}-T_{\rm eff,T}$$\rangle$ = +79 K with $\sigma$ = 97 K (23 stars). These are selected from the catalogue of Cayrel de Strobel et al. (1992), and redetermined from colours in the rare case that they are initially determined from Fe I lines (Thévenin 2003, private communication). Comparing for the remaining stars, we found similar values to those in the references used for the metallicities, with a difference of less than 80 K, apart for the star HD 128279, for which we obtained 5130 K, and Gratton et al. (2000) obtain 5394 K. This latter value is derived from dereddened (B-V)0 and (b-v)0 colours using colour- $T_{\rm eff}$ transformations (Kurucz 1995), and is higher, as for nearly all the stars in common with the present sample.


  \begin{figure}
\par\includegraphics[width=8.65cm,clip]{0913fig3.eps}\end{figure} Figure 3: Comparison of the $T_{\rm eff}$ obtained in this study with spectroscopic $T_{\rm eff}$ from Fulbright ( triangles) and Mishenina & Kovtyukh ( squares).
Open with DEXTER

The surface gravities (log g) were estimated using Hipparcos parallaxes (Perryman et al. 1997) through the relation proposed by Nissen et al. (1997), [g] = [M] + 4 [ $T_{\rm eff}$] + 0.4 ( $M_{\rm bol}-M_{\rm bol,\odot}$), where [x] = log $x/{x_{\odot}}$. The solar bolometric magnitude $M_{\rm bol,\odot}$ = 4.75 mag was taken from Allen (1976), and the bolometric magnitudes of the stars were determined using interpolated bolometric corrections from the grids given by Alonso et al. (1999a): $M_{\rm bol}$ = K + BC(K) + 5 log $\pi$ + 5, where $\pi$ is the trigonometric parallax. For the two stars without accurate Kand $K_{\rm s}$ magnitudes BC(V) was interpolated from the grids of Alonso et al. (1995, 1999b). The relative standard error in parallax, $\sigma(\pi)/\pi$, was typically of order 0.10 and smaller than 0.30. The stellar masses (M) were derived from their position in the (log  $T_{\rm eff}$, log  $L/L_{\odot}$) plane, using the isochrones from Bergbusch & VandenBerg (2001). Most of the stars were found to be in the main sequence and the subgiant domain, while the remainder were in the red giant branch. The average stellar mass for the sample was 0.78 $M_{\odot}$. The error in log g arose mainly from the parallax uncertainty: an error of $\sigma(\pi)/\pi$ = 0.10 implied (via the bolometric magnitude) $\Delta$log g = 0.10 dex, while typical errors in the mass $\sigma(M)$ = 0.045 $M_{\odot}$ and in the effective temperature $\sigma$( $T_{\rm eff}$) = 75 K implied errors of only $\Delta$log g = 0.025 dex.


  \begin{figure}
\par\includegraphics[width=8.8cm,clip]{0913fig4.eps}\end{figure} Figure 4: Comparison of the log g obtained in this study with spectroscopic log g from Fulbright ( triangles) and Mishenina & Kovtyukh ( squares).
Open with DEXTER

The surface gravities of the dwarf stars agreed in general with the LTE spectroscopic values from Fulbright (2000) and Mishenina & Kovtyukh (2001), see Fig. 4. However, those of the very metal-poor subgiants HD 87140 and G170-47 were different, with $\Delta$log g > +0.40 dex. This discrepancy between trigonometric and LTE spectroscopic surface gravities in late-type stars is discussed by Nissen et al. (1997) and Allende Prieto et al. (1999). Moreover, trigonometric surface gravities agree better with NLTE spectroscopic determinations (Thévenin & Idiart 1999). Since no NLTE value is published for HD 87140, we could only compare with a trigonometric value, for instance that given by Gratton et al. (2000) and we found a difference of 0.05 dex. For G170-47, we could compare with the NLTE "best choice'' determination from Israelian et al. (2001) and found a difference of 0.07 dex. As for the remaining stars, we found log g values very similar to those in the references we used for metallicities, with differences of less than 0.17 dex.

Table 1: Stellar atmospheric parameters and Cu and Zn abundances. Some parameters are from 1) Castro et al. (1999); 2) Chen et al. (2000); 3) Fulbright (2000); 4) Spite & Spite (1982); 5) Gratton et al. (2000); 6) Mishenina et al. (2002); 7) Thévenin (1998); 8) Israelian et al. (2001); 9) Israelian & Rebolo (2001). Solar Cu and Zn abundances are from Anders & Grevesse (1989).

For the three very metal-poor dwarfs G64-12, G275-4 and LP815-43 observed with VLT/UVES, the atmospheric parameters were taken from Israelian et al. (2001). For the two very metal-poor giants BD -18 $\hbox {$^\circ $ }$5550 and HD 2796, the spectroscopic log g and  $T_{\rm eff}$ were taken from Israelian & Rebolo (2001). Finally, the microturbulence was fixed at $V_{\rm t}$ = 1.0 km s-1 for the dwarfs and 1.5 km s-1 for the giants.


  \begin{figure}
\par\includegraphics[width=14cm,clip]{0913fig5.ps}\end{figure} Figure 5: Synthetic spectral fits to the copper and zinc lines in the spectra of several metal-poor stars of our sample.
Open with DEXTER

   
3.1 Cu and Zn abundances

The abundances of Cu and Zn are listed in Table 1. The wings of the intense Cu I lines were fitted with the synthetic spectrum by differential analysis with respect to the Sun. Because these lines were in the damping or saturation parts of the curve of growth for stars with solar temperature and metallicities -2 < [Fe/H] < 0, the effects of changing the abundances were carefully checked. In general an uncertainty of 0.1 dex was adopted for the derived abundances. For metallicities lower than -2 the lines were unsaturated, and the adjustments could be achieved more precisely.

Table 2: Sensitivities of the copper and zinc abundances to variations in the atmospheric parameters of the stars HD 94028 ( $T_{\rm eff}$ = 6058 K, log g = 4.28, $V_{\rm t}$ = 1.0 km s-1), HD 140283 ( $T_{\rm eff}$ = 5723 K, log g = 3.68, $V_{\rm t}$ = 1.0 km s-1) and G170-47 ( $T_{\rm eff}$ = 5154 K, log g = 2.93, $V_{\rm t}$ = 1.25 km s-1), and to the overshooting.

The uncertainties in $T_{\rm eff}$ and log g led to average errors of respectively 0.11 and 0.02 dex in the Cu abundance, and 0.035 and 0.02 dex in the Zn abundance, while the $V_{\rm t}$ uncertainty led to an average error of 0.07 in the Cu abundance and 0.05 dex in the Zn abundance. The main contribution to the error came from the uncertainty in the localization of the continuum, especially in the spectra of the more metal-rich stars and in the noisy spectra; this uncertainty implied an error of 0.1-0.2 dex in the abundances. Additionally the noise in the copper and zinc lines led to an error of 0.05-0.1 dex. In the case of the Cu I $\lambda$3273.95 Å line, an extra error of 0.01-0.04 dex was taken into account for the uncertainty in the contribution of the weak Fe I $\lambda$3273.95 Å line.

A few spectra also permitted the analysis of the other Cu I line at 3247.54 Å. In the case of the very metal-poor dwarfs G64-12, G275-4, LP 815-43, an agreement to within 0.15 dex was found between the Cu abundances given by the two lines. Those in the spectrum of the giant BD -18 $\hbox {$^\circ $ }$5550, however, provided a difference of 0.25 dex, suggesting that the differential analysis with respect to the Sun was less reliable for this star (note that for the other giant, HD 2796, the difference was still greater, reason why its Cu abundance was not listed in Table 1). As for the zinc line, the spectra of only six stars out of the thirteen with [Fe/H] < -2 permitted a clear detection. The spectra of the remaining stars presented a very weak feature at 3302.58 Å, from which no reliable abundance could be derived; thus upper limits to the abundances were established. Figure 5 shows the fits for three metal-poor stars.


  \begin{figure}
\par\includegraphics[width=12cm,clip]{0913fig6.eps}\end{figure} Figure 6: [Cu/Fe] and [Zn/Fe] ratios vs. [Fe/H]. The upper limits of Zn abundances are represented by the arrows, the Cu abundance in the giant BD -18 $\hbox {$^\circ $ }$5550 by the square and the fit to the Cu data obtained by Mishenina et al. by the dotted line. The solid lines represent the predictions of the Prantzos model and the horizontal lines represent the elapsed time since the beginning of the halo formation.
Open with DEXTER

Since Kurucz's models atmospheres without overshooting give a better reproduction of colour indices, Balmer profiles and Procyon data than Kurucz's models with overshooting (Castelli et al. 1997), we checked the changes produced on the abundances using the former models. For the stars HD 94028, HD 140283 and G147-70 (Table 2), these models implied Cu and Zn abundances systematically lower by about 0.15 and 0.10 dex, respectively, differences small enough to preserve the main properties of our results discussed in Sect. 4.

As found by Cayrel et al. (2003), the abundances derived from the spectra of metal-poor stars may be overestimated if the continuum scattering is approximated as an additional opacity source in the spectral synthesis code. The overestimation is especially high for $\lambda$ < 400 nm, where continuum scattering becomes important relative to continuous absorption. Since the spectral synthesis code MOOG approximates scattering by absorption, we evaluated the influence on the abundances by using another LTE spectral synthesis code, TurboSpectrum (Alvarez & Plez 1998), which takes proper account of the continuum scattering. For the only two giants, HD 2796 and BD -18 $\hbox {$^\circ $ }$5550, the Cu abundances derived using TurboSpectrum were lower by more than 0.2 dex than those derived using MOOG, while the Zn abundances were the same. As for the dwarfs G275-4 and G64-12, the abundances for both elements were higher by less than +0.1 dex. We therefore concluded that, for the dwarfs, the Cu and Zn abundances obtained from the UV lines were not significantly affected by this source of errors, while for the giants only the Zn abundances were. We thus removed the contribution of the giants from our conclusion on the Cu trend.

We compared our Cu and Zn abundances with those of Mishenina et al. (2002) obtained from optical lines, and found a good agreement: $\Delta$[Cu/H] = $\langle[{\rm Cu/H}]-[{\rm Cu/H}]_{M}\rangle$ = -0.04 dex, with a standard deviation $\sigma$ = 0.15 dex (16 stars; no giant data to compare), and $\Delta$[Zn/H] = 0.03 dex, with $\sigma$ = 0.13 dex (16 stars), respectively. Converting our values with our abundance sensitivities to their adopted stellar atmospheric parameters provided a similar agreement: $\Delta$[Cu/H] = -0.09 dex, with $\sigma$ = 0.14 dex, and $\Delta$[Zn/H] = -0.02 dex, with $\sigma$ = 0.13 dex, respectively. The slight abundance decrease was due principally to their lower effective temperatures and microturbulence. The dwarf HD 19445 (see Fig. 5) was the most metal-poor star for Cu abundance comparison with Mishenina et al. (2002), with [Fe/H] = -2.04. We found [Cu/H] = -2.69, a value 0.37 dex lower than that of Mishenina et al. (2002), and lower by $\sim$0.1 dex than the upper limit estimate of Sneden et al. (1991) (for their adopted atmospheric parameters). We also compared our results with the abundances obtained from optical lines by Sneden et al. (1991), and we found a rough agreement for Cu, $\Delta$[Cu/H] = 0.13 dex, with $\sigma$ = 0.36 dex (5 stars only[*]; no giant data to compare), and for Zn, $\Delta$[Zn/H] = 0.13 dex, with $\sigma$ = 0.13 dex (8 stars), respectively. Converting our values with our abundance sensitivities to their adopted stellar atmospheric parameters also revealed a rough agreement for Cu, $\Delta$[Cu/H] = 0.17 dex, with $\sigma$ = 0.24 dex, and a good agreement for Zn, $\Delta$[Zn/H] = 0.04 dex, with $\sigma$ = 0.14 dex, respectively.

4 Trends of [Cu/Fe] and [Zn/Fe]

4.1 Results

Our derived [Cu/Fe] and [Zn/Fe] ratios are shown in Fig. 6. The [Cu/Fe] ratio is almost solar down to [Fe/H] $\sim$ -1 and then decreases, until it reaches a plateau $\langle$[Cu/Fe]$\rangle$ = -0.98 at [Fe/H] < -2.5 (the ratio is averaged over the five stars that are not giants). This slanted s-shape trend in dwarfs is similar to that outlined by Mishenina et al. (2002) (see Fig. 6). However, the [Cu/Fe] plateau is $\sim$0.3 dex lower than the plateau obtained in the range -2.7 < [Fe/H] < -1.7 in their halo giants and in the globular clusters studied by Simmerer et al. (2003). It could be noticed that underabundances near [Cu/Fe] = -1 are also found in Mishenina et al. (2002) (the two dwarfs BD +41 $\hbox {$^\circ $ }$3931 and HD 140283[*], Fig. 9, with [Cu/Fe] = -0.97 and $\sim$-0.95, respectively) and in Sneden et al. (1991) (the giants HD 2665 and HD 122563, Fig. 7, with [Cu/Fe] = -0.87 and -0.93). In addition, Sneden et al. (1991) provide an upper limit [Cu/H] < -3.93 for the giant BD -18 $\hbox {$^\circ $ }$5550, and we find [Cu/H] = -3.94 ([Cu/Fe] = -0.93), a value that would be even lower if we were to take into account a better treatment of the UV spectral synthesis (see Sect. 3.1).

Considering zinc, the [Zn/Fe] ratio is essentially solar down to [Fe/H] $\sim$ -2.0, with a possible excess of $\langle$[Zn/Fe]$\rangle$ = +0.18 at lower metallicities. It is interesting to compare our results with existing high quality surveys which include observations of Zn: those of Mishenina et al. (2002) and Gratton et al. (2003) in a metallicity range comparable to that of this study, and that of Cayrel et al. (2003), which is complementary and covers considerably lower metallicities. The comparison of our data with those of these three studies is shown in Fig. 7. The ensemble of all the data gives a clear picture of the evolution of the Zn/Fe ratio with metallicity. A solar ratio of Zn/Fe is supported by the data of both Mishenina et al. (2002) and Gratton et al. (2003) down to a metallicity of about [Fe/H] $\sim$ -2.5. The rise, which is hinted at by the lower metallicity tail of our sample, is consistent with the robust rise from the Cayrel et al. (2003) sample. The regression line for [Zn/Fe] provided in Table 7 of Cayrel et al. (2003) intercepts the solar ratio for [Fe/H] $\sim$ -2.1. It is worth noting that our data are consistent with the data of both Mishenina et al. (2002) and Gratton et al. (2003) over the common metallicity range, and that our Zn abundances are derived from Zn I 3302.58 Å, while the other groups use optical lines. Two stars, HD 166913 and HD 132475, however, seem to defy agreement, by presenting [Zn/Fe] greater than 0.45. In the case of HD 166913, the apparent high ratio [Zn/Fe] = 0.45 stems not from the Zn abundance itself, which is 0.11 dex above that derived by Gratton et al. (2003) (once the different atmospheric parameters are taken into account, using the abundance sensitivities given in Table 2 for HD 94028), but rather from our much lower Fe abundance, adopted from Thévenin et al. (1998) and 0.23 dex lower than Gratton et al.'s (2003). For the latter comparison, we use the abundance sensitivities from Gratton et al. (2003; Table 9), since none are given in Thévenin et al. (1998). In the case of HD 132475, the high ratio [Zn/Fe] = 0.50 stems from its reddening, E(b-y) = 0.046 (Schuster & Nissen 1989), because this reddening implies an appreciably higher effective temperature and thus a lower ratio [Zn/Fe] (if referred to Table 9, Gratton et al. 2003).


  \begin{figure}
\par\includegraphics[width=8.1cm,clip]{0913fig7.ps}\end{figure} Figure 7: [Zn/Fe] ratios vs. [Fe/H] from our measurements ( filled dots), from Cayrel et al. ( filled triangles), from Mishenina et al. ( asterisks) and Gratton et al. ( crossed squares). The [Zn/Fe] values of Gratton et al. are decreased by 0.09 dex for normalization to meteoritic rather than to solar abundances.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=7.9cm,clip]{0913fig8.eps}\end{figure} Figure 8: [Cu/Zn] ratios vs. [Zn/H]. The solid line represents the predictions of the Prantzos model.
Open with DEXTER

4.2 Cu and Zn nucleosynthesis

In spite of considerable theoretical effort, the production sites of Cu and Zn are not clear (e.g. the discussion in ishenina et al. 2002, and references therein).

The very different behavior of Cu and Zn implies that n-capture during He-burning cannot be the main mode of production for both. Concerning the production of Zn, there are reasons to believe that neither the s- nor the r-process can be very relevant (Cayrel et al. 2003). If the s-process were important, the [Zn/Fe] ratio ought to decrease with decreasing metallicity, as is the case for typical s-process products such as Y, Sr or Ba. From Fig. 7, it is instead clear that this ratio is initially constant and then increases at lower metallicities. This observational argument supports the conclusions drawn from the theoretical computations of Gallino et al. (1998). Concerning the relevance of the r-process, we may note that in the very metal-poor star CS 31082-001 ([Fe/H] = -2.90), in which all the r-process elements are greatly enhanced, [Zn/Fe] = +0.18 (Hill et al. 2002), remarkably similar to the behaviour of other stars of comparable metallicity.

We may also extend the above argument to Cu; in fact, no Cu is detected in CS 31082-001, although no upper limit is published. A greatly enhanced Cu would surely be detected in the high quality data of Hill et al. (2002). Therefore we may conclude that strong Cu production in the r-process is unlikely.

We compare our results with new chemical evolution predictions from Prantzos (2003, private communication), which follow those from Goswami & Prantzos (2000). The model is based on the stellar initial mass function (IMF) from Kroupa et al. (1993) and assume an evolution of the halo for a gas outflow rate which is equal to 8 times the star formation rate. We find that the Prantzos model (2003) agrees most reasonably with the observations - especially for Cu, except for some points (see Fig. 6). Model and observations agree on the existence of a plateau for Cu at very low metallicities ([Fe/H] < -2.5). This supports the adopted SNII yields (Woosley & Weaver 1995) averaged over the stellar initial mass function from Kroupa et al. (1993). Nevertheless, we stress that since higher ratios than those predicted are also observed at low metallicity (Sneden et al. 1991; Mishenina et al. 2002), we can only conclude on a possibly lower limit for the relative production of copper. As for Zn, we observe local agreements and a systematic disagreement appearing at [Fe/H] < -2, indicating, to a first approximation, that some Zn sources are missing.

5 Conclusions

In summary, we present Cu and Zn abundances for a sample of 38 stars, principally dwarfs, with metallicities ranging from solar to -3. The abundances were obtained using the intense UV lines Cu I 3273.95 Å and Zn I 3302.58 Å, and were generally in agreement with the results from previous studies. We showed that the Cu abundance indicator is adequate for the study of abundances in very metal-poor stars, permitting us to detect [Cu/Fe] $\sim$ -1 in dwarfs at [Fe/H] = -3. This research can be extended to dwarf stars of even lower metallicities and with high resolution spectra. Contrary to the slanted s-shape trend found for [Cu/Fe], the [Zn/Fe] trend is approximately solar in most stars with [Fe/H] > -2; this is a common result, which supports the use of Zn as a metallicity tracer in damped Ly$\alpha$ systems. However, at [Fe/H] < -2, Zn appears to be slightly overabundant ($\langle$[Zn/Fe]$\rangle$ = +0.18) and it is in fact at these metallicities that the rise also becomes noticeable in the sample of stars studied by Cayrel et al. (2003).

The trends found can be reasonably reproduced by the halo chemical evolution predictions from Prantzos (2003, private communication). The ratio [Cu/Fe] $\sim$ -1 for -3 < [Fe/H] < -2.5 is predicted in the Prantzos model for SNII yields from Woosley & Weaver (1995). It probably represents, at this evolutionary stage, a lower limit of Cu production, since greater ratios are also observed in previous studies. The [Zn/Fe] trend is in good general agreement with that predicted, except in the low-metallicity tail, where the predicted ratios should be systematically increased.

Acknowledgements
We would like to thank Peter A. Bergbusch and Don A. VandenBerg for providing the software producing the isochrones, and François Thévenin for information about the Catalogue III/193. We would also like to thank Robert L. Kurucz for information about the intense copper lines and their hyperfine structure, Nicolas Prantzos for providing the unpublished trends of copper and zinc, and Bertrand Plez for the TurboSpectrum code. This research has made use of the SIMBAD database, operated at CDS, Strasbourg, France, and also of data results from the Two Micron All Sky Survey, which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis Center, funded by the National Aeronautics and Space Administration and the National Science Foundation.

References

 

Copyright ESO 2004