A&A 465, 587-591 (2007)
DOI: 10.1051/0004-6361:20066709

Lithium abundances in metal-poor stars[*]
(Research Note)

J. R. Shi1,2 - T. Gehren 2 - H. W. Zhang2,3 - J. L. Zeng1,2 - G. Zhao1


1 - National Astronomical Observatories, Chinese Academy of Sciences, Beijing 100012, PR China
2 - Institut für Astronomie und Astrophysik der Universität München, Scheinerstrasse 1, 81679 München, Germany
3 - Peking University, Beijing 100871, PR China

Received 7 November 2006 / Accepted 2 January 2007

Abstract
Aims. Lithium abundances for 19 metal-poor stars are determined using high-resolution spectroscopy. The abundances of stars on the lithium plateau are discussed.
Methods. All abundance results are derived from NLTE statistical equilibrium calculations and spectrum synthesis methods.
Results. In agreement with previous analyses it is found that excitation and de-excitation due to hydrogen collisions are negligible for the lithium line formation process, while charge transfer reactions are an important source of thermalization. However, the resulting NLTE effects on the determination of lithium abundances for metal-poor stars are negligible (<0.06 dex).
Conclusions. The mean lithium abundance for stars on the lithium plateau determined from NLTE analyses is A(Li) $\sim$ 2.26, while it is 2.21 dex when charge transfer reactions are included. The latter result enhances the discrepancy between the observed lithium abundances and the primordial lithium abundance as inferred by the WMAP analysis of the cosmic microwave background. This discrepancy may be explained by metal diffusion.

Key words: line: formation - line: profile - stars: abundance - stars: late-type - galaxies: evolution

1 Introduction

Lithium abundances in metal-poor stars are of importance for cosmology, the evolution of the Galaxy, Big Bang nucleosynthesis and stellar structure. Spite & Spite (1982) first detected that the 7Li abundance in halo dwarfs is independent of temperature and metallicity, suggesting that they represent the primordial abundance from Big Bang nucleosynthesis (BBN). However, based on the WMAP results (Wilkinson Microwave Anisotropy Probe, see Spergel et al. 2003), the BBN prediction for primordial 7Li is between 2 and 3 times higher than observed in the atmospheres of halo dwarfs (Coc et al. 2004; Cuoco et al. 2004; Cyburt 2004; Cyburt et al. 2005), and it remains to be resolved what type of systematic effects are responsible for this discrepancy.

Norris et al. (1994) and Thorburn (1994) have found a dependence of A(Li) on both temperature and metallicity, which was confirmed by later work, such as that of Ryan et al. (1996, 1999, 2001), Zhang & Zhao (2003), and Boesgaard et al. (2005). However, Molaro et al. (1995) and Bonifacio & Molaro (1997) questioned this trend. Their result is also supported by Meléndez & Ramírez (2004) and Charbonnel & Primas (2005). The latter authors noted that for some stars, common to several analyses, very discrepant Li abundances were reported, which could have influenced some of the early claims for a significant abundance dispersion. They suggested that this could also play a role in the current debate about the existence of a slope of lithium abundances with temperature and metallicity.

As Li I has a low ionization potential, non-local thermodynamic equilibrium (NLTE) effects can be expected to be important (Carlsson et al. 1994; Asplund et al. 2003). Carlsson et al. (1994) have shown that the LTE approximation - when deriving lithium abundances for cool stars - is not representative of the physical properties of the atmospheric layers in which the 6707 Å line forms. Since NLTE corrections vary in sign and size when spanning a large range of stellar parameters, ignoring these corrections may affect any interpretation of the A(Li)-plateau and its possible slope with effective temperature and metallicity. Recently, Barklem et al. (2003) found that the collisional excitation and de-excitation (Li$^\ast+$ H  $ \rightleftharpoons$ Li$^\ast+$ H) reactions are negligible for the line formation process, while the charge transfer reactions (Li$^\ast$ + H  $ \rightleftharpoons$ Li+ + H-) are important in thermalizing lithium. Asplund (2005) suggested that new 1D calculation including the charge transfer data and improved line-blocking treatment would be worthwhile.

In this work we present new results of NLTE studies for metal-poor stars. In Sect. 2 we comment on the observational techniques. The models are discussed in Sect. 3, while in Sect. 4 we describe the methods of analysis and uncertainty for our results. The results and discussions are given in Sect. 5, and the summary is presented in the last section.

2 Observations

The spectra of our samples were obtained from 1997 to 2003 with the fiber-coupled Cassegrain échelle spectrograph FOCES (Pfeiffer et al. 1998) mounted on the 2.2 m telescope of the Calar Alto Observatory. The standard configuration was a 20482 24 $\mu$ CCD that covered the wavelength region 3900-9000 Å with $\lambda/ \Delta\lambda \sim 40~000$. Some stars were observed with the 15 $\mu$ CCD resulting in $R \sim65~000$. All stars were observed at least twice with a signal-to-noise ratio ranging from $S/N \sim 100$ up to $S/N \sim 400$. Our sample contains main sequence/turnoff stars presented in Gehren et al. (2004, 2006).

Data extraction followed the standard automatic IDL program environment designed for the FOCES spectrograph (Pfeiffer et al. 1998). All echelle images including flatfield and ThAr were corrected for bias and scattered light background. Objects and ThAr exposures were extracted and corrected for flatfield response. Bad pixels were detected and as far as possible removed by comparison of at least 2 single exposures (see Gehren et al. 2004, for details). Figure 1 shows a portion of the spectrum of HD 19445 around 6707 Å.


  \begin{figure}
\par\includegraphics[width=8cm,clip]{6709fig1.eps}\end{figure} Figure 1: The Li I line near $\lambda $ 6707 Å in the program star HD 19445.
Open with DEXTER

3 NLTE line formation

3.1 Model atom

The lithium model atom employed in this paper is illustrated by its Grotrian diagram in Fig. 2. We extend the number of levels to n=10. The model is closed with the Li+ ground state containing 49 levels and all transitions between them. Fine structure splitting is neglected with the exception of ${2{\rm p}}~^2{\rm P} {1/2}^{\rm o}$ and ${2{\rm p}}~^2{\rm P} {3/2}^{\rm o}$, and a balance between these two levels is assumed.


  \begin{figure}
\par\includegraphics[angle=90,width=8.5cm,clip]{6709fig2.eps}\end{figure} Figure 2: Grotrian diagram for Li I.
Open with DEXTER

The level energies and oscillator strengths are taken from the TOPBASE data when available. All other energy levels and transition probabilities are assumed to follow a hydrogen-like approximation. We adopt the photoionization cross-sections provided by the Opacity Project via the TOPBASE system (Cunto & Mendoza 1992) for the 25 terms available; for the remaining terms hydrogen-like approximations are used. As usual, background opacities are calculated with an opacity sampling code based on the line lists made available by Kurucz (1992).

As in previous calculations, all bound-bound electron collision strengths for allowed transitions are approximated by the simple formula of van Regemorter (1962). Forbidden collisions in our representation follow Allen's (1973) formula for which no experimental verification is available. Ionization by electronic collisions is calculated from Seaton's (1962) classical path approximation. The collision rate is therefore proportional to the photoionization cross-section at the threshold.

The cross-sections for excitation and de-excitation due to collisions with neutral hydrogen are taken from Drawin (1968, 1969); here we adopt a multiplication factor of ${\it SH} = 0.05$. As already noted by Asplund et al. (2003) and Barklem et al. (2003), hydrogen collisions are negligible for the Li line formation process. The cross-sections for charge transfer reactions are taken from Barklem et al. (2003).

3.2 Model atmospheres

Our analyses are all based on the ODF line-blanketed LTE model atmospheres generated and discussed by Fuhrmann et al. (1997). Their main characteristics are:

3.3 Stellar parameters

For most of the stars we adopt the stellar parameters determined by Korn et al. (2003) and Gehren et al. (2006), where the effective temperatures are derived from the wings of the Balmer lines; surface gravities are based on the HIPPARCOS parallaxes. Iron abundances are based on Fe II lines, and the microturbulence velocities are estimated by requiring that the iron abundances derived from Fe II lines should not depend on equivalent width. The uncertainties for the temperature, surface gravity, metal abundance and microturbulence velocities are generally assumed to be $\pm$80 K, 0.05 dex, 0.05 dex and 0.1 km s-1 respectively.

Following previous work, we employ the macroturbulence parameters $\zeta$ in the radial-tangential form and adopt values as described in Gray (1984). The projected rotational velocity is obtained as the residual to the observed line profiles, after having removed the known instrumental profile obtained from the Moon spectra. It is, however, not used in this analysis.

3.4 Atomic line data

We adopt the oscillator strengths from the NIST data base, namely $\log gf = 0.002$ and -0.299 for the 6707.76 and 6707.91 Å lines, respectively. While collisional broadening parameters describing van der Waals interaction with hydrogen atoms are taken from Barklem et al. (1998), we note that the lines are not sensitive to the damping constant either.

Table 1: Atmospheric parameters and lithium abundances.

4 Lithium abundances

4.1 Abundance determination

Abundance determinations are carried out using a spectral synthesis method based on the IDL/Fortran SIU software package of Reetz (1993). The synthetic spectra are convolved with macroturbulence, rotational and instrumental broadening profiles in order to match the observed spectral lines. Lithium abundances are obtained with a best fit to the observed line profile. They are listed in Table 1, together with the stellar parameters. Three different sets of lithium abundances are reported for each star, namely LTE and NLTE, the latter without and with charge reaction processes. The internal errors in the derived lithium abandances are less than 0.05 dex. An example of comparison of synthetic profiles with the observed ones for lithium lines is shown in Fig. 3.


  \begin{figure}
\par\includegraphics[width=8cm,clip]{6709fig3.eps}\end{figure} Figure 3: Synthetic NLTE flux profiles of Li lines compared with the observed spectra. The solid lines are the observed spectra, while the dashed lines are NLTE profiles.
Open with DEXTER

Without the charge reactions, our results show that the NLTE corrections are negligible (about 0.02 dex) for stars on the plateau, while they are about 0.1 dex for low temperature metal-poor main sequence and subgiant stars, i.e., HD 45282, and they change their sign with the stellar parameters. Due to the overpopulation of the ground state, the NLTE effects are all negative, and they decrease the lithium abundances of the plateau stars by about -0.05 dex when charge reactions are included.

4.2 Uncertainty in fundamental parameters

It is known that lithium abundances are very sensitive to the effective temperature chosen, whereas their dependence on the other parameters, i.e. gravity, metallicity and microturbulence is negligible. Common uncertainties on $\log g$, [Fe/H] and $\xi$ ($\pm$0.07 dex, $\pm$0.05 dex and $\pm$0.1 km s-1, respectively) will not affect the final lithium abundances. Fields et al. (2005) discussed the influence of the temperature scale for halo dwarfs on lithium abundances and chemical evolution in detail. An uncertainty of $\pm$80 K in $T_{\rm eff}$ for stars of our sample translates into about $\pm$0.06 dex in the lithium abundance.

5 Results and discussion

Lithium abundances of metal-poor stars have been determined by several groups. We have 16 stars in common with Charbonnel & Primas (2005) and five stars with Meléndez & Ramírez (2004). The average differences between our [Li/H] and theirs are $0.10 \pm 0.10$ and - $0.08 \pm 0.07$ for NLTE without the charge reactions. We note that the temperatures adopted for our sample stars are slightly higher than those of Charbonnel & Primas (+ $33 \pm 110$ K), but lower than determined by Meléndez & Ramírez (- $51 \pm 112.52$ K); part of the differences can be explained by the different temperatures adopted. We have four stars in common with Asplund et al. (2006) and two stars with Boesgaard et al. (2005); the average differences between our [Li/H] and theirs are $0.07 \pm 0.04$ and $0.03 \pm 0.03$, respectively. Though not perfect, the agreement is satisfactory.

Figure 4 shows the overall distribution of A(Li) in our sample as a function of temperature. It is obvious that the lithium abundances depend on temperatures for stars with $T_{\rm eff} < 6000$ K, whereas no such trend can be found for stars with $T_{\rm eff} > 6000$ K.


  \begin{figure}
\par\includegraphics[width=7.5cm,clip]{6709fig4.eps}\end{figure} Figure 4: Li abundances as a function of temperature. Pluses (+) are the LTE lithium abundances; filled circles ($\bullet $) are the NLTE lithium abundances without charge transfer and open circles ($\circ $) are the lithium abundances with charge transfer.
Open with DEXTER

The distribution of A(Li) with [Fe/H] for all stars in the sample is shown in Fig. 5, while Fig. 6 is restricted only to stars with $T_{\rm eff}$ greater than 6000 K. From these two figures, it is clearly seen that A(Li) is independent of metallicity for our sample. The abundances obtained in this work for plateau stars show a very small star-to-star scatter. This is a consequence of employing internally very accurate temperatures based on the Balmer lines. For stars with $T_{\rm eff} > 6000$ K, the mean A(Li) is $2.27 \pm 0.04$, $2.26 \pm 0.04$ and $2.21 \pm 0.04$ for LTE, and for NLTE without and with charge reactions, respectively. Our Li abundance is higher than that of Ryan et al. (1999), but lower than that of Meléndez & Ramírez (2004).


  \begin{figure}
\par\includegraphics[width=7.5cm,clip]{6709fig5.eps}\end{figure} Figure 5: A(Li) abundances analysis as a function of [Fe/H]. Symbols are the same as in Fig. 4.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=7.5cm,clip]{6709fig6.eps}\end{figure} Figure 6: Same as Fig. 3, but only with the stars with $T_{\rm eff} > 6000$ K. Symbols are the same as in Fig. 4.
Open with DEXTER

When charge reactions are included, our NLTE results are in good agreement with those of Barklem et al. (2003), namely the derived lithium abundances are lowered by 0.05 ...0.10 dex. This trend enhances the discrepancy between the observed lithium abundances in metal-poor stars and the primordial lithium abundance as inferred from the WMAP analysis of the cosmic microwave background. As often discussed, a possible reason may be the depletion of atmospheric 7Li. Due to turbulent mixing, mass loss and/or atomic diffusion, uniform depletion factors of order 0.2 dex (Vauclair & Charbonnel 1998) and as high as 0.3 dex (Salaris & Weiss 2001; Richard et al. 2005) have been discussed. This possibility faces the strong constraint that the observed lithium abundances show very little dispersion. It is unlikely that the stellar processes that depend on temperature, mass, and rotation velocity of the star would destroy 7Li always by the same amount. Coc et al. (2004) suggested that the 7Be(d,p)$2\alpha$ reaction would solve the lithium problem if its cross-section were very much higher than assumed. However, recent experimental results show this cross-section is ten times overestimated (Angulo et al. 2005).

Recently, Korn et al. (2006) suggested that metal diffusion may be responsible for the low apparent stellar lithium abundances in the atmospheres of old stars. Their atomic diffusion model, added to standard stellar evolution models, leads to a decrease of the original lithium abundance by about 0.25 dex. At the time of formation of these very old stars, therefore, the cosmic matter would have had a lithium abundance of A(Li) $\simeq$  $ 2.21 + 0.25 = 2.46 \pm
0.04$, a value in good agreement with the primordial lithium abundance based on the cosmic microwave background and Big Bang nucleosynthesis, A(Li) WMAP $\simeq$  $ 2.54 \pm 0.10$.

6 Conclusions

We have determined the lithium abundances of 19 metal-poor stars spanning the range -2.5 < [Fe/H] < -1.0. All abundances are derived from NLTE statistical equilibrium calculations. Based on our results, we come to the following conclusions.

1.
For plateau stars our lithium abundances are independent of the temperatures and metallicities.

2.
The lithium abundances obtained in this work for plateau stars show a very small star-to-star scatter. The mean lithium abundance obtained for the plateau stars are A(Li) = 2.26 and 2.21 dex without and with the charge reactions included. Our result is higher than the primordial lithium abundance proposed by Ryan et al. (2001) (A(Li) $\approx$ 2.0 $\ldots 2.1$ dex), while it is lower than that obtained by Meléndez & Ramírez (2004).

3.
NLTE corrections are relatively small (<0.1 dex). NLTE abundances are decreased compared to LTE by about -0.05 for all sample stars when charge reactions are included.

4.
While metal diffusion seems to remove most of the problems related to the cosmic lithium abundance plateau, it is important to note the remaining influence of 3D hydrodynamic models on the final value of the lithium abundance.

Acknowledgements
S.J.R. thanks the Institute of Astronomy and Astrophysics of Munich University for warm hospitality during a productive stay in 2004. We thank the referee for helpful comments. This research was supported by the National Natural Science Foundation of China under grants Nos. 10433010 and 10521001, and by the Deutsche Forschungsgemeinschaft under contracts Ge 490/26-1 and 446 CHV112.

References

 

Copyright ESO 2007