A&A 442, 961-992 (2005)
DOI: 10.1051/0004-6361:20042491
C. Charbonnel1,2
-
F. Primas3,2,
1 -
Geneva Observatory, 51 chemin des Maillettes, 1290 Sauverny,
Switzerland
2 -
Laboratoire d'Astrophysique de Toulouse et de Tarbes, CNRS UMR 5572, 14 Av. E.
Belin, 31400 Toulouse, France
3 -
European Southern Observatory, Karl-Schwarzschild Str. 2,
85748 Garching b. München, Germany
Received 7 December 2004 / Accepted 8 May 2005
Abstract
Thanks to the accurate determination of the baryon density of the universe
by the recent cosmic microwave background experiments, updated predictions
of the standard model of Big Bang nucleosynthesis now yield the initial
abundance of the primordial light elements with unprecedented precision.
In the case of 7Li, the CMB+SBBN value is significantly higher than the
generally reported abundances for Pop II stars along the so-called Spite
plateau. In view of the crucial importance of this disagreement, which has
cosmological, galactic and stellar implications, we decided to tackle the
most critical issues of the problem by revisiting a large sample of
literature Li data in halo stars that we assembled following some strict
selection criteria on the quality of the original analyses.
In the first part of the paper we focus on the systematic uncertainties
affecting the determination of the Li abundances, one of our main goal
being to look for the "highest observational accuracy achievable" for
one of the largest sets of Li abundances ever assembled. We explore in
great detail the temperature scale issue with a special emphasis on reddening.
We derive four sets of effective temperatures by applying the same
colour-
calibration but making four different assumptions about
reddening and determine the LTE lithium values for each of them. We compute
the NLTE corrections and apply them to the LTE lithium abundances.
We then focus on our "best'' (i.e. most consistent) set of temperatures
in order to discuss the inferred mean Li value and dispersion in several
and metallicity intervals. The resulting mean Li values along the
plateau for [Fe/H]
-1.5 are
and
when the lowest
effective temperature considered is taken equal to 5700 K and 6000 K
respectively. This is a factor of
2.48 to 2.81
(depending on the adopted SBBN model and on the
effective temperature range chosen to delimit the plateau)
lower than the CMB+SBBN determination. We find no evidence of intrinsic
dispersion.
Assuming the correctness of the CMB+SBBN prediction, we are
then left with the conclusion that the Li abundance along the plateau is
not the pristine one, but that halo stars have undergone surface depletion
during their evolution.
In the second part of the paper we further dissect our sample in search
of new constraints on Li depletion in halo stars. By means of the Hipparcos
parallaxes, we derive the evolutionary status of each of our sample stars,
and re-discuss our derived Li abundances. A very surprising result emerges
for the first time from this examination.
Namely, the mean Li value as well as the dispersion
appear to be lower (although fully compatible within the errors) for the
dwarfs than for the turnoff and subgiant stars.
For our most homogeneous dwarfs-only sample with [Fe/H] -1.5,
the mean Li abundances are
and
when the lowest
effective temperature considered is taken equal to 5700 K
and 6000 K respectively. This is a factor of 2.52 to 3.06
(depending on the selected range in
for the plateau and on the SBBN
predictions we compare to) lower than the CMB+SBBN primordial value.
Instead, for the post-main sequence stars the corresponding
values are
and
,
which correspond to a
depletion factor of 2.28 to 2.52.
These results, together with the finding that all the stars with Li abnormalities (strong deficiency or high content) lie on or originate from the hot side of the plateau, lead us to suggest that the most massive of the halo stars have had a slightly different Li history than their less massive contemporaries. In turn, this puts strong new constraints on the possible depletion mechanisms and reinforces Li as a stellar tomographer.
Key words: stars: abundances - stars: Population II - stars: evolution - Galaxy: abundances - Galaxy: halo - cosmology: early Universe
From the time of its discovery in the stellar atmospheres of very metal-poor Population II stars (Spite & Spite 1982a), lithium has been considered a key diagnostic to test and constrain our understanding and description of the primordial Universe, of stellar interiors and evolution, and of spallation physics (for the latter two especially so if combined with abundances of beryllium and boron).
Table 1:
WMAP+SBBN primordial abundances predicted when
(or
0.0009;
Spergel et al. 2003) is adopted.
Lithium-7 is one of the four primordial isotopes that have been formed
in observable quantities by nuclear reactions during the first minutes
of the Universe (e.g. Olive et al. 2000, and references therein).
Together with deuterium, helium-3 and helium-4, knowledge of its primordial
abundance provides one of the main observational constraints on the
baryon-to-photon ratio
which is the only
free parameter of the standard Big Bang Nucleosynthesis (SBBN).
Among these light elements, lithium is one of the easiest to observe (its resonant doublet falls at 670 nm, i.e. easily accessible from ground telescopes even of small sizes), which explains the wealth of data available in the literature.
Another reason for such a large literature database is connected to the
important finding of a remarkably flat and constant Li abundance among
Galactic halo dwarf stars spanning a wide range of effective temperatures
and metallicities (the so-called Spite plateau, cf. Spite & Spite 1982a,b).
This result came as a surprise. At that time indeed it was generally believed
that the primordial A(Li) abundance
was in the range 3.0-3.3, which corresponds to the value measured in
meteorites and also to the maximum value detected in Population I stars,
both in the field and in open clusters (see the review by Boesgaard &
Steigman 1985). If this were the case, then the "constant'' but lower value
of lithium along the Spite plateau would have required that all the oldest
stars of our Galaxy had suffered a uniform depletion by about a factor of 10. Although several mechanisms could be conjectured to modify the surface
lithium abundance (proto-stellar destruction, microscopic diffusion,
turbulent diffusion, mass loss), they were also suspected to depend on the
stellar mass (i.e., on the effective temperature). Consequently the
constancy of the Li plateau was used (and is actually still used very
often) as an argument to say that these processes were in fact not efficient
in Pop II stars.
The other interpretation was then that the plateau value represents the
amount of Li
produced during the Big Bang, and that the Galaxy had been enriched in its
Li content by a factor of at least 10 since its birth.
The fact that lithium is produced in several other nucleosynthetic sites (i.e.
fusion, spallation reactions, late stellar evolutionary
stages, like AGB stars, novae, etc.; see Romano et al. 2003, and references
therein), none of which has been quantitatively and accurately estimated
nor strongly constrained by observations, complicates the final
interpretation of its Galactic evolution.
Recent results on cosmic microwave background anisotropies, most particularly
from the Wilkinson Microwave Anisotropy Probe (WMAP) experiment
(Bennet et al. 2003; Spergel et al. 2003) allowed an unprecedented precision
on the determination of the baryon-to-photon ratio
and revealed that
Li seems to lie between the two extreme
solutions discussed above. The WMAP data alone lead to
,
or
.
When combined with additional CMB experiments (CBI, Pearson et al. 2003;
ACBAR, Kuo et al. 2002) and with measurements of the power spectrum (2dF
Galaxy Redshift Survey, Percival et al. 2001; Ly
forest, Croft
et al. 2002; Gnedin & Hamilton 2002), the resulting
values are
or
.
With this value of
,
updated SBBN predictions now allow a precise
determination of the primordial abundances of the light elements D, 3He,
4He and 7Li that we can compare with observations in low-metallicity
environments. The WMAP+SBBN determinations of these abundances in the
most two recent studies (Coc et al. 2004; Cyburt 2004; Serpico et al. 2004) are summarised in
Table 1.
A very good agreement is achieved between the primordial abundance of
deuterium derived from WMAP+SBBN and the average value of D/H
observations in cosmological clouds along the line of sight of quasars
(Kirkman et al. 2003). On the other hand the observational data of 3He
in galactic HII regions are scarce and must be corrected for contamination
of the observed gas by ejecta from earlier generations of stars (e.g. Tosi
1998; Charbonnel 2002). The upper limit to the primordial abundance
recommended by Bania et al. (2002) is however quite consistent with the
CMB-derived value.
Finally the CMB-predicted primordial 4He abundance is higher than the
values derived from the determinations in complex low-metallicity HII regions
(both galactic and extra-galactic) and the extrapolation to zero oxygen
abundance (Izotov & Thuan 2004; Olive & Skillman 2004 and references
therein). However the difference is relatively modest (2-3)
and it may
simply call for further exploration of the systematic effects in the
abundance analysis.
The most critical case concerns 7Li, the CMB-derived primordial abundance
of which is clearly higher (by about a factor of 3) than the current
determinations in low-metallicity halo stars.
This result seems to be very robust with respect to the nuclear uncertainties
on the SBBN reactions although Coc et al. (2004) show that the discrepancy
could be resolved by an increase of a factor of 100 of the
7Be(d,p)24He reaction rate. Although this is not supported by the
data currently available, this issue has to be further investigated
experimentally. Should this nuclear solution be excluded, we would then be
left with the astrophysical solution.
Namely, with the conclusion that the Li abundance that we see at the surface
of halo stars is not the pristine one, but that these stars have undergone
surface lithium depletion at some point during their evolution. This
possibility has been discussed many times in the literature. Several physical
mechanisms have been invoked, but all the current models encounter
considerable difficulties to reconcile a non negligible
depletion of lithium with both the flatness and the small dispersion along
the so-called Li plateau (see the review by Pinsonneault et al. 2000; and
Talon & Charbonnel 2004, for more recent references).
The challenge thus still remains to identify the process (or processes) by
which a reduction by a factor of
3 could occur so uniformly in stars
over a large range in effective temperature and metallicity.
With the CMB constraint, we are now entering a golden age for Li as both a
baryometer and a stellar tomographer.
In this quest however one has still to pay special attention to the
observational analysis and determination of the lithium abundances in
the most metal-poor, thus the oldest stars of our Galaxy. As a matter of
fact and despite the large number of spectroscopic data that has become
available in the last two decades, there are still on-going debates on
the patterns of the plateau, like its thinness, the possible existence
of a spread and of a dependence of A(Li) with metallicity and effective
temperature. These characteristics must be precisely determined in order
to constrain the physical processes which lead to Li depletion in Pop II
stars as well as those of Galactic production.
Deliyannis et al. (1993) were among the first to present
evidence for the existence of dispersion (of the order of 20%
about the mean, derived from the "equivalent width-colour'' plane), followed
by Thorburn (1994) and Norris et al. (1994), the latter being the first to
also have found a dependence of A(Li) on both
and metallicity.
Molaro et al. (1995) counter-argued these findings showing that when a fully
consistently determined temperature scale is used (in their case, the Fuhrmann
et al. 1994 scale, derived from Balmer lines fitting), no dispersion nor
tilt is found (Li abundances are mostly sensitive to
), but they were
shortly followed by Ryan et al. (1996) who once again confirmed the slopes.
They argued that the Molaro et al. sample was plagued by the inclusion of
subgiants that may have affected their final outcome. On the intrinsic scatter
issue, they noted that there were some stars, characterized by very similar
parameters (colour,
and metallicity), but that turned out from
multiple measurements to have very different Li abundances. The debate kept
being very alive: Spite et al. (1996) explored further the
scale
issue (comparing different
determinations) and found that the rms
scatter of the Li abundance was between 0.06 and 0.08 dex, hence very small
if real. Bonifacio & Molaro (1997) re-selected their sample, this time
excluding possible outliers (like the abovementioned subgiants) and once
again came to the conclusion of no intrinsic dispersion nor dependence of
A(Li) on metallicity, but of a tiny trend with the temperature. They concluded
that the finding or not of a trend with effective temperature may well depend
on the adopted
scale.
It is only towards the end of the 90s when a
full agreement on the absence of intrinsic dispersion was reached: Ryan et al. (1999), by analysing a new sample of 23 stars covering a narrow range
in
(6050-6350 K) and in metallicity (-3.5--2.5), claimed
that the intrinsic spread is effectively zero, i.e. 0.031 dex, at the
1
level (to be compared to their formal errors of 0.033 dex).
However, they still recovered the dependence on metallicity, at the level
of dA(Li)/d[Fe/H] = 0.118
0.023 (1
)
dex per dex, i.e. very
similar to the slopes previously found. The trend with
, if any, is
likely to be meaningless because of the very narrow
range there
explored.
This is when we started developing our project. By comparing the data samples
analysed by different authors (starting with the Bonifacio & Molaro 1997;
and Ryan et al. 1996, 1999, because of their final, opposite claims)
we noticed that for some stars, common to several analyses, very discrepant
Li abundances were reported, which could have clearly influenced some of the
early claims for dispersion, and they could still play a role in the current
debate about the existence of a slope between Li and
and [Fe/H].
In order to further tackle these issues, we decided to re-analyse the large
sample of Li abundances available in the literature (Sects. 3 and 4) from a
different perspective. One of our main goals is to focus on the systematic
uncertainties affecting the determination of Li abundances. First, because
the Li abundance is strongly dependent on the assumed temperature, we explore
further the temperature scale issue (Sect. 5) with the aim of deriving what
the best achievable accuracy may be for a temperature scale derived in a
consistent manner. We put special emphasis on reddening, usually an
underestimated source of error.
We derive four sets of effective temperatures by applying the same
colour-
calibration but making four different assumptions about
reddening.
We then derive the LTE lithium abundances for each of these sets and compute
the NLTE corrections. Then we select our "best'' (in terms of consistency)
set of temperatures in order to
determine the mean Li abundance and the dispersion for one of the largest
sample of halo stars ever studied in a consistent way
(Sects. 6 and 7). This allows us to derive preliminary results on the mean
lithium abundance and dispersion which can be compared to previous analyses
(Sect. 8).
Secondly we look afresh at our Li abundances together with the evolutionary status of each target (Sect. 9) in order to get clues on the internal processes that may have been involved in modifying the Li abundances along the Plateau. This is why we did not restrict a priori our sample to any specific evolutionary status. We then discuss the lithium abundance along the plateau for the dwarf stars only (Sect. 10) and look at its behaviour in subgiants (Sect. 11). We test whether our results on the dispersion and trends can be altered by the presence of binary stars in the sample (Sect. 12), and we finally inspect the cases of stars with extreme lithium abundances (Sect. 13). Then we discuss the current status of Pop II stellar models in view of our observational results (Sect. 14). Finally we summarise our results and conclude on some remaining open questions (Sect. 15).
The database of Li abundances measured in Galactic stars and available in the literature is huge. During this work, we restricted our search to the main observational analyses from the early 90 onwards.
Table 2: Literature sources.
Table 3: The data sample and its main characteristics, as found in the literature. The whole table is available on-line.
We assembled our final data sample from 13 literature sources (cf. Table 2).
The main ones were the works by Ryan et al. (1996, 1999, hereafter
R96, R99), Bonifacio & Molaro (1997, hereafter BM97), and Pilachowski et al. (1993, hereafter PSB93), the last one in order to include a large set
of subgiant stars. R96, R99, and PSB93 analysed newly observed spectra and
derived new temperatures, whereas BM97 collected a large sample of stars for
which the equivalent width of the Li I line had already been measured and
re-computed the Li abundance based on a new temperature scale. This first
list of targets was then
complemented by objects taken from Hobbs & Thorburn (1991), Thorburn (1994,
for those very few stars which had not been already re-analysed by Ryan et al. 1996), Molaro et al. (1995), Spite et al. (1996), Ryan & Deliyannis
(1998), Gutierrez et al. (1999), Fulbright (2000), Ryan et al. (2001a,b),
and Ford et al. (2002). We gave our preference to original works, i.e.
works that added new observations or that re-analysed literature values based
on a new temperature scale. This exercise left us with a sample of 146 stars,
covering the metallicity range between [Fe/H] = -1.0 and -3.5, the
temperature interval
= 4500-6500 K, and the surface gravity range
between log g = 3.0-5.0. In other words, we are sampling the main sequence,
subgiant and giant evolutionary stages. Table 3
presents the data sample and its main characteristics
in terms of nomenclature (including cross-identifications), stellar parameters
and the equivalent widths of the Li I line as found in the literature. For
those objects, for which multiple determinations are available, the minimum
and maximum values are listed. For a more detailed comparison, Fig. 1
shows how the equivalent widths measured (used) by some of the literature
sources listed in Table 1 compare to each other. For the purpose of this test,
we selected those works that had the largest number of stars in common.
For completion, we note that there have been three other recent works that have also made use of the large database of Li measurements available from the literature. Two of them had different scientific goals and they both used (after a critical selection) the Li abundances as found in the literature: Romano et al. (1999) re-assessed the Galactic evolution of lithium, whereas Pinsonneault et al. (2002) compared the most recent Li abundances to theoretical predictions of models including rotational mixing and examined them for trends with metallicity. This is why they do not appear in our list of literature sources. The third, most recent and most similar work to ours is the one from Meléndez & Ramírez (2004), who studied the behavior of the A(Li) plateau and its trends in a sample of 41 dwarf stars. An improved InfraRed Flux Method-based temperature scale was derived (Ramírez & Meléndez 2005a,b, hereafter RM05a,b) and used to compute the Li abundances (from equivalent widths taken from the literature). Because RM05a,b became public at a late stage in our refereeing process, we did not update our input targets list, but instead we decided to discuss and compare their and our results when relevant.
Analysing lithium is not very difficult. Lithium appears in a stellar
absorption spectrum with few transitions, namely the resonant line at
670.7 nm, and a much weaker signature at 610.4 nm, only recently
explored in the most metal-poor stars (cf. Bonifacio & Molaro 1998;
Ford et al. 2002). The 670.7 nm line falls in a clean spectral region,
especially in metal-deficient stars. From its equivalent width it is
easy to derive an abundance, once the stellar parameters of the object
under investigation have been determined. The sensitivity of the final
Li abundance to surface gravity, metallicity, and microturbulence is
not very significant (see below), whereas an uncertainty of 70 K
in
(commonly quoted as a reasonable uncertainty on this parameter,
for solar-type stars) translates into
0.056 dex on the
final lithium abundance. Because the effective temperature
is clearly the most critical parameter, it will be discussed separately, in
the next section, where we provide a more detailed description of what we
have learned from its derivation. Here, we will briefly comment on the other
input stellar parameters and how they were determined.
The surface gravity was first determined from an inspection of the (b-y) vs. c1 diagramme (cf. Fig. 2), which allowed us to assign a preliminary log g value to each of our stars. These first-guess values were first checked versus those quoted in the literature sources used to assemble our sample. Then we finally attributed to each star the log g value deduced from its position in the Herzsprung-Russel diagram (see Sect. 9.2).
![]() |
Figure 1: The comparison of EWs as measured/adopted by different authors for stars in common. The literature sources used for these comparisons are reported in the upper left and lower right corners for the y- and x-axis respectively. |
Open with DEXTER |
![]() |
Figure 2: The entire sample plotted in the (b-y) vs. c1 diagram from which we assigned the first set of gravities. Different symbols identify different values of log g, as summarised in the lower left corner of the figure. |
Open with DEXTER |
In order to evaluate the sensitivity of the derived A(Li) abundances on the
stellar gravity, we performed our tests at
= 5250 K and 6000 K, for
log g = 3.0 and 4.0, and for two metallicities, [Fe/H] = -1.0 and -2.5. On
the average, we found that a change of 1 dex in log g affects the lithium
abundance by
0.018 dex only, with very little dependence on the
effective temperature, or on the equivalent width of the lithium line. Our
findings agree very well with the common statements that the dependence of
A(Li) on the stellar gravity is negligible in the error budget (e.g.
Meléndez & Ramírez 2004).
Table 4: Photometry and reddening excesses. The whole table is available on-line.
The metallicity was first selected from the same literature sources
from which we assembled our data sample. As one can see from inspecting
Table 3 the agreement between different literature sources (when available)
is in general quite satisfactory. Only in few cases, a large discrepancy
is present but was fortunately solved because observed spectra were in hand.
One such example is the star HIP 88827, for which BM97 reported
[Fe/H] = -0.91 (taken from Alonso et al. 1996) and R96 measured -2.10:
from a high resolution, high S/N UVES spectrum taken at the VLT, its
metallicity has been recently derived from a reliable set of Fe II lines
and found to be -2.4 (Nissen et al. 2002). Because of the general good
agreement, our final metallicities are simply the weighted mean of all the
values available for a given star (with the exception, of course, of the very
few discrepant cases mentioned above). All the [Fe/H] values determined
spectroscopically from high resolution, high S/N spectra were double weighted.
The uncertainty on the final [Fe/H] values is the
of the
weighted mean (cf. Table 6).
We cross-checked our metallicities also photometrically, via the calibration
of Schuster & Nissen (1989, hereafter SN89 - cf. Eq. #(3)), and taking
advantage of the availability of ubvy-
photometry for the entire
sample. The satisfactory agreement we found was judged more than sufficient
for our test-purposes, therefore we did not explore any further the possible
systematic differences between metallicities derived spectroscopically and
photometrically.
Although the Li abundance is only slightly dependent on the adopted
metallicity (0.2-0.3 dex uncertainties in [Fe/H] affect the Li abundance
less than 0.01 dex, cf. R96), one has to remember that a possible error in
the metallicity may affect the A(Li) vs. [Fe/H] trend, and more importantly
the
derived for that star and consequently its lithium abundance.
Our calculations confirm what already found by R99: if
[Fe/H] is off by 0.15 dex (a reasonable uncertainty for this parameter,
especially since we have assembled our sample from a variety of analyses),
the effect on
is almost negligible (
20 K) which implies
an uncertainty on the Li abundance of less than 0.02 dex. This
sensitivity applies, of course, to the parameters space spanned by our
stars, with a tendency of finding larger dependences as the metallicity
increases.
A similar, almost negligible dependence, is found also between Li and microturbulence, for which 0.5 km s-1 in
correspond to
0.005 dex in A(Li). Because of this negligible dependence and after some
checks of previous literature works, that included both dwarf and (sub)giant
stars, we decided to run all our calculations assuming
km s-1. Because
of the very small dependence of A(Li) on microturbulence this choice gives
identical results to what was implemented by PSB93, who let
varying
smoothly between 1.0 and 2.0 km s-1 going from the hotter to the cooler stars
of their sample.
The main goal of any lithium analysis is to determine a fully consistent temperature scale for all the targets under examination, as the lithium abundance is strongly dependent on this stellar parameter. This approach is usually considered a guarantee of the absence of spurious differences possibly arising by having applied different criteria to the derivation of the effective temperature. Ideally, one would like to determine this parameter from first principles, i.e. to derive direct temperatures for metal-poor dwarfs. In practice, this has been achieved sofar for very few and very bright targets only (cf. RM05a, for a summary of what is currently available).
Stellar temperatures can be determined spectroscopically (e.g. via profile
fitting of the wings of some of the Balmer lines, or from minimising the
slope between the iron abundance - as derived from Fe I lines - and their
excitation potential), or from photometry. Since we have assembled our data
sample from the literature (i.e. no newly observed spectra), photometry is
the only choice we have. Furthermore, because of the size of the sample only
Strömgren uvby-
photometry (among the photometric indices most
sensitive to stellar temperatures) is available for all our stars, thanks
to the extensive photometry by Schuster & Nissen (1988), supplemented
by unpublished photometry by Schuster (private communication) for
approximately 10% of our stars. Table 4 (available in its entirety only
on-line) summarises all the photometry we have used, together with the
different colour excesses we have derived and that we will now discuss.
For reference and test purposes, it also includes (B-V) values
taken from Hipparcos (ESA 1997) and from the literature, but we remind the
reader that (B-V) is not a good temperature indicator.
We derived
from the (b-y)0 colour index using the IRFM calibrations
of Alonso et al. (1996, 1999, plus the erratum from 2001) for dwarf and giant
stars (cf. their Eqs. #6 and #14 respectively). The evolutionary status
assigned to each object for the determination of the gravity (see
previous section) was used to decide if a star was a dwarf or a post-main
sequence object.
In order to overcome the known problem of Alonso's calibration (i.e.
diverging towards high values at the lowest metallicities, cf.
R99, Nissen et al. 2002), we adopted a lower limit of
in the
equation. Two comments are mandatory here. Firstly, although it is
important to keep in mind this a posteriori solution when discussing
the findings on the most metal-poor stars of our sample, we note that Nissen
et al. (2002) found no worrysome behavior (due to the assumption of this lower
limit on [Fe/H]) when comparing effective temperatures derived via Alonso's
colour-
calibrations from (b-y) and (V-K). Secondly, we note that RM05b
have argued for the first time that this characteristic of the Alonso's
calibration is not a problem, in the sense that it is not a numerical
artifact due to the quadratic dependence on [Fe/H] but it is intrinsic to
the IRFM. In their recent work, they see a similar effect not only in the
vs. (b-y) plane, but also for other colour indices. If confirmed,
this would imply that adopting a lower limit on the metallicity is not
justified, thus different effective temperatures would be derived. How
different can be seen from Fig. 3, where we plot the effective
temperatures of all our stars with
derived from applying
or not this lower limit at -2.1. Knowing what the dependence of A(Li) on
is, one can already have an idea of what the effect will be on our
final Li abundances. We will come back to this when discussing our results.
![]() |
Figure 3:
A direct comparison of effective temperatures derived from the
Alonso's calibration, but adopting a lower limit of
![]() ![]() |
Open with DEXTER |
The interstellar reddening excess was estimated from the
calibration of SN89, including a zero-point correction of +0.005 mag (Nissen
1994). The sample of stars from which Schuster & Nissen (1988) derived the
abovementioned relation span the metallicity range
,
and specific ranges of Strömgren colour indices, namely
,
,
,
and
.
However, despite (b-y), c1, m1 indices are available for all our stars,
some of them fall outside the ranges of validity for the calibrating
equations we have planned to use, and in few cases the
index is
missing. For consistency, one is then left with two
possibilities: a) to reduce the sample to only those objects for
which the complete set of valid Strömgren photometric indices is
available (our sample would then be reduced to approximately 90 stars);
b) to replace the missing information with other methods and carefully
investigate the effects of such "pollution'' (unavoidable in order
to keep the number statistics high) on the final output(s). Option a)
clearly represents the simplest path, and it will be used as our benchmark.
Here, instead, we want to describe in some detail the series of compromises
(i.e. sample pollutions) we had to introduce in order to keep working with
our whole data sample. At the end, we will compare a) to b), and discuss
how these two approaches affect the final Li abundances.
Estimating the reddening correction is probably the weakest point of any
photometrically-based
scale. Interstellar reddening excesses are
rarely quoted and discussed in spectroscopic analyses, despite they can
affect significantly the determination of any stellar parameter, the
effective temperatures in particular. In order to be consistent within
the Strömgren photometry framework, we decided to use the
relation to evaluate the interstellar reddening excess. However, for 24 stars (
20% of the sample), the
index was not found. A
common solution shared by several analyses has been to assume zero
reddening, based on the fact that the stars likely belong to the solar
neighborhood. Alternatively, one could derive E(b-y) from other colour
excesses, if available. Either way, this is an approximation, that we
consider as the first compromise on our data-sample. When referring to
this sub-sample of stars, we will call it the
-sample.
We decided to derive E(b-y) from E(B-V) via the formula
(Crawford 1975), which is based on a 1/
reddening law and
on the central wavelengths of the bandpasses. For the E(B-V) colour
excesses we simply averaged the E(B-V) values derived from IR Dust Maps
(Schlegel et al. 1998, hereafter S98) and the models of large scale visual
interstellar extinction by Hakkila et al. (1997, hereafter H97). If these
two sources (which will be extensively discussed in Sect. 5.4) were found to
diverge significantly (e.g. in the case of HIP 49616, for which we find 0.177
from S98 maps and 0.024 from H97 models), that object was dropped from our
list. The
values were given a much lower weight, being their
original source not always available. However, when found in agreement with
the other two E(B-V) sources, they were included in the straight average.
Because this solution includes a mixture of E(B-V) sources, from now on we
will refer to these values as
.
Out of the 24 stars we have
without the
index, 14 were rescued. Figure 4 shows a direct
comparison between (b-y)0 values which have been corrected for reddening
derived respectively from
(on the x-axis) and from
values and Crawford's formula (on the y-axis). This plot includes all the
stars of our sample for which both methods could be safely applied. There
is clearly some scatter around the 1:1 relation (on the order of 0.015 mag),
and a systematic tendency of deriving larger reddenings from the
indirect formula.
![]() |
Figure 4:
A direct comparison of the (b-y)0 index, de-reddened respectively
from the ![]() ![]() |
Open with DEXTER |
We mentioned above that the equation to be used in the derivation of the
interstellar reddening is valid (i.e. has been tested) only for specific
ranges of the Strömgren indices. Taking the photometry at face value,
our sample includes a total of 35 stars, for which at least one of the
Strömgren indices does not fulfill these criteria: 7 stars have the
(b-y) index outside the 0.254-0.550 interval, 5 stars have the c1 index
outside 0.116-0.540, another 8 have the m1 index falling outside
0.033-0.470, and 21 have the index outside 2.55-2.68 (six of
which also have some of the other indices off).
For the 15 (i.e., 21-6) stars with
(only) outside the allowed
range, we were able to apply the same solution described in Sect. 5.1, i.e.
we derived E(b-y) from
values, to 6 of them. If the remaining
9 stars are dropped, together with those 20 (7+5+8) that have (b-y) or m1 or
c1 outside the allowed ranges, we are then left with a total of 111 stars, of
which 20 are "polluted'' because their E(b-y) was derived from
(of the initial 39 = 24+15 stars that belonged to this sample, 5 were dropped
because one or more of their Strömgren indices fell outside the allowed
ranges). We note, however, that so far we took all the photometric indices
at face value, whilst each of them has its own associated uncertainty. If
one were to take this into account, then the application of the validity
ranges would allow some flexibility. For some stars (7 out of 20) such
approach seems very reasonable: HIP 103337, for instance, has the (b-y) index
off by 0.002 mag, which translates to a 10 K effect; all "m1'' drop-outs
(except two) have their photometric index off by only few thousandths of a
magnitude (at most by 0.006 mag), which affect their final effective
temperatures between 3.5 K and 20 K. Because of these considerations, we
decided to keep these 7 objects in a separate (polluted) sample, which we
call the ubvy-sample. In summary, we are then left with a total of 118 stars.
Should we have considered also the
values with a 3-digits precision,
six more targets would now belong to the ubvy-sample. Since this could
influence our final discussion of the A(Li) plateau, we will take a closer
look at them once their Li abundances have been derived (cf. Sect. 7).
Knude (1979) showed that interstellar reddening is caused primarily by
small dust clouds with a typical reddening of
.
If
true, this would make any correction for reddening values smaller than
0.03 almost meaningless. This is why in the past it has been common
practice to correct only those stars for which the derived excess was
comparable to or larger than this value. SN89, for
instance, chose
E(b-y)=0.025 as their reference value and performed
the corrections
(b-y)0 = (b-y) - E(b-y),
c0 = c1 - 0.2E(b-y), and
m0 = m1 + 0.3E(b-y) for all stars with
.
However, the commonly accepted picture of the nature and appearence of
the interstellar reddening has significantly evolved since Knude's work
and seems to suggest a patchy distribution of interstellar dust (implying
that values smaller than
E(b-y)=0.03 may be real) with a void of about
70-75 pc around the Sun (e.g. Lallement et al. 2003). In principle, one
should then apply the interstellar reddening excesses derived from the
(b-y)0 relation to all the stars. In practice, one has to face the
problem of correcting also those stars for which negative E(b-y) values
are derived (28 objects in our case). This means switching to the index as the main
indicator, thus losing precision for those stars
that are close enough to be in the void around the Sun (Nissen, priv.
comm.). This explains why it is still common procedure to choose a minimum
E(b-y) threshold below which the colour indices are not corrected for
reddening excesses. For instance, Nissen et al. (2002) have arbitrarily
chosen
E(b-y)=0.015, which corresponds to twice the sigma of E(b-y).
Being unsure of what is the best approach to follow, we tested the effects
of the abovementioned solutions by deriving different sets of
temperatures, under the following assumptions: 1. all stars were de-reddened
[
(1)]; 2. all stars were dereddened except those with negative E(b-y)values [
(2)]; 3. all stars were dereddened except those with
[
(3)]; 4. all stars were dereddened except those with
negative E(b-y) values and those with
and d<70 pc
[
(4)]. Our four derived sets of effective temperatures are compared
in Fig. 5 (
(1) always on the x-axis) and show very little differences,
the smallest in the case of
(1) compared to
(2) (
27 K around
the mean, cf. top panel). Looking at these two sets of
in
more detail (the crosses in the top panel of Fig. 5 identify the 28 differing
stars), one notices that for one third of these stars there is practically
no difference (these are the stars for which very small reddening excesses,
in absolute value, were found to be negative). Except for 3 targets, for
which the difference in effective temperature is larger than 100 K, all
the others are within
50 K, with a systematic tendency of the
(2)
values to be slightly higher than
(1) ones.
![]() |
Figure 5:
How the sets of temperatures
![]() ![]() ![]() ![]() |
Open with DEXTER |
The fact that a very small difference in E(b-y) (such as 0.005 mag,
which accounts approximately for half of the common uncertainty) translates
already into 35 K in effective temperature is a strong indication that
accounting for interstellar reddening excesses plays an important role in
deriving an accurate photometrically-based
scale - especially when
the abundance of the element(s) under investigation is very sensitive to
like in the case of lithium. Therefore, it is important to comment
also on the other sources of interstellar reddenings available in the
literature.
In order to complete our comparison tests, we decided to derive reddening values also from the most recent tabulated values, i.e. those derived from infrared mapping of the dust emission distribution (Schlegel et al. 1998, S98 for short), and from the models of large scale visual interstellar extinction (Hakkila et al. 1997, hereafter H97). This was done not only to thoroughly test our photometrically-based reddening values, but especially to evaluate a posteriori the effect of mixing different sources of reddening on the derived temperature scale (a common approach to many spectroscopic abundance analyses).
Schlegel et al. (1998) estimated the dust column densities from the COBE/DIRBE
and IRAS/SISSA infrared maps of dust emission over the entire sky, and
transformed them to reddenings by using colours of elliptical galaxies. In
other words, these maps give reddening values as if the objects lie outside
the Galaxy, hence they may overestimate the real reddening, especially for
relatively nearby objects. Also, at low latitudes ()
the removal
of IR point sources is not optimal, hence the derived reddening values may
be strongly affected. The quoted errors are of the order of 16%.
Hakkila et al. (1997) instead, developed a numerical algorithm to model the
large scale Galactic clumpy distribution of obscured interstellar gas and
dust by using published results of large-scale visual interstellar extinction.
It is concentrated towards the Galactic plane and it varies as a function of
Galactic longitude and latitude. These estimates depend on the assumed
distance, which is one of the input parameters to the algorithm. A word of
caution concerns its inability in identifying small scale (less than 1 )
extinction variations, and the fact that reddening estimates for mid-Galactic
latitudes (
)
and for distances between 1 and 5 kpc
in the Galactic plane are more unsecure. The quoted errors are not very
meaningful since they represent the mean of the errors as reported in the
original studies, hence they are likely overestimated.
Table 4 presents an overview on how reddening excesses derived from different
methods compare to each other. Columns 2 to 7 report Strömgren photometry
and
values as derived from the (b-y)0-
calibration
(see previous sub-sections), while Cols. 8 and 9 list the Hipparcos-
and literature-based (B-V) values. Column 10 reports the E(B-V) values as
derived from the S98 maps, while Col. 12 lists the values as derived from
the H97 algorithm. The corresponding E(b-y), derived via the Crawford (1975)
relation are reported in Cols. 11 and 13 respectively.
Two remarks are important. First, some very high values of E(B-V) are derived from the S98 maps (by using the IDL code made available by the same authors), which look unrealistic when compared to all the other E(B-V) values. Since the main purpose of our tests is to have some feeling on the possible scatter introduced by mixing reddening values taken from different sources, without being biased by outliers, we did not investigate these high values any further. Hence, they have been discarded from all comparison figures and tests. However, one can easily identify them in Table 4, since no corresponding E(b-y) value was derived (same applies also to Hakkila-based E(B-V) values - though for significantly fewer stars). One should also note that in Table 4 there remain some suspiciously high E(B-V) values.
![]() |
Figure 6:
Comparison between E(B-V) values derived from different methods:
filled dots show how E(B-V) colour excess compare when derived from the
InfraRed Dust Maps of Schlegel et al. (1998, on the x-axis) and of
Burstein & Heiles (1982, on the y-axis); crosses represent a similar
comparison between
![]() ![]() |
Open with DEXTER |
![]() |
Figure 7:
Comparison between
![]() |
Open with DEXTER |
Secondly, in order to survey as many choices of reddening as possible,
Table 4 includes also E(B-V) values as found in the literature sources from
which we assembled our data sample
(Col. 14 labeled
)
and as
derived from the neutral hydrogen H I column density distribution
of Burstein & Heiles (1982, BH for short -
Col. 15 labeled
)
in correlation with deep galaxy counts. A partial summary of Table 4 is
provided in Fig. 6, where two comparisons are plotted simultaneously: the
filled circles show the relation between E(B-V) values derived from the S98
(on the x-axis) and from the BH maps (on the y-axis). The crosses represent
the comparison between E(B-V) values derived from the S98 IR dust maps (on
the x-axis) and from the H97 model (on the y-axis). The 1:1 relation is
plotted for comparison. From this figure, one immediately notices the lack
of a tight correlation between the E(B-V) values derived from the maps of
Schlegel et al. and Burstein & Heiles. The comparison with H97 clearly
shows that Hakkila's model of the Galactic interstellar dust distribution
tends to give reddening values much lower than the S98 IRDM. Another way to
look at these comparisons is by inspecting Fig. 7, where the difference
between the different reddening estimates (always with respect to S98 values)
are plotted versus
.
From this figure it becomes clear that
there are systematic differences between S98 reddening values and the other
two estimates, although there does not seem to be any dependence on the
distance, except for a slightly larger dispersion at small distances (cf.
Fig. 8).
![]() |
Figure 8: Comparison between E(B-V) values derived from various sources (identified from the acronym given in the right bottom corner of each pan) and the distance of the object, as derived from Hipparcos parallaxes. |
Open with DEXTER |
In summary, based on the abovementioned arguments, it is very hard to defend
the position that we know reddening better than 0.007-0.010 mag (2)
which correspond already to affecting the temperature determination by
50-70 K (and in turn the lithium abundance by
0.05 dex).
![]() |
Figure 9:
How our final set of temperatures
![]() |
Open with DEXTER |
Table 5: Final temperatures as derived from different reddening corrections. The whole table is available on-line.
Table 5 summarises
all sets of temperatures that have resulted from our several
digressions: Cols. 2 to 5 report the effective temperatures derived by applying the reddening
relation and making different assumptions about it (cf.
Sect. 5.3), whereas Cols. 6 and 7 represent sets of temperatures as
derived by using different sources of reddening. The latter two columns
are useful comparisons to check how much scatter and slope on the
A(Li)-plateau could originate just by mixing different reddening sources
in the same analysis. Figure 10 clearly shows how different the derived
temperatures can be when the reddenings are derived from different methods:
if they are taken from the IRDM of S98, for instance, the derived temperatures
are significantly higher (+178 K on the average, with a dispersion around
the mean of
236 K). Of course, depending on for which targets one may
need to select a different source of reddening excess, there may well be some
artificial scatter emerging among their Li abundances.
![]() |
Figure 10:
Comparison between our final set of effective temperatures
![]() ![]() ![]() |
Open with DEXTER |
Another relevant comparison to make would be the one between our
scale and the
values used in the original literature sources
from which our list of targets was assembled. We present these comparisons
in Fig. 11, in the form of
(2) -
(X) vs.
(2), where
(2)
is our preferred and finally selected set of temperatures and
(X)
represents the temperature scales used in the original works from which our
list of targets was assembled (cf. Table 2). Here, we selected to plot those
literature analyses from which we took the largest numbers of stars, except
for the bottom panel where the comparison with the very recent RM05a
scale is presented. We remind the reader that our list of targets
does not include objects from this work because the analyses were carried
out almost simultaneously. As one can see, our temperature scale (
(2))
tend to be always higher, the largest difference being with PSB93 (on average
+158 K, with a dispersion around the mean of
136 K). The smallest
differences, on average, are between us and BM97 (+6 K) and RM05a
(practically zero), but the dispersions around these means remain on the
order of
100-150 K.
Also, in Table 5 all our targets are grouped in three
different sub-samples. As stated at the beginning of Sect. 5, in order to be
fully consistent with the analytical method one chooses to follow, one is
usually forced to work with a much smaller sample of stars compared to the
initial data-set: in our case, 91 stars compared to the original 146. In
order to avoid this drastic reduction, we explored alternative solutions
which allowed us to retain a larger number of objects
(118), but at the price of contaminating part of the sample as follows:
Sample #1 is the clean sample: it includes 91 stars, for which the complete set of Strömgren photometric indices are available, and for which the SN89 calibrations can be successfully applied.
Sample #2 is the
sample: it includes 20 stars, for
which the reddening value E(b-y) was derived from averaging different
sources of E(B-V) values (S98, H97, and literature - cf. Sect. 5.1 for details
on how this average was performed), via Crawford's formula (1975). We
note that no correction was applied to these stars to compensate for the
offset seen in Fig. 4.
Sample #3 is the ubvy sample: it includes 7 stars, for which one of the ubvy photometric indices (b-y, c1, m1) falls just slightly outside (see last paragraph of Sect. 5.2) the allowed intervals for the application of the Nissen & Schuster (1989) calibrations.
The lithium abundance for all the stars of our sample was determined from
the equivalent widths (EWs) of the 670.7 nm line as reported in the
literature works from which we assembled the data sample. Table 6 lists the mean equivalent width and its
1
uncertainty that were used in our computations of the Li abundance:
we opted for the mean value because of a satisfactory overall agreement found
in the literature (cf. Table 3, for the corresponding references, listed in
the last column of the table).
The lithium abundance was derived under the assumption of Local Thermodynamic Equilibrium (LTE) using Kurucz (1993) WIDTH9 and model atmospheres with the overshooting option switched off (cf. Castelli et al. 1997 for the models, and Molaro et al. 1995 for comparisons between different versions of Kurucz model atmospheres). The gf-value we used for the Li I is 0.171, for which the VALD database reports an accuracy of 3%.
As Carlsson et al. (1994) have shown, the LTE approximation when deriving Li
abundances for cool stars is not a realistic representation of the physics
present in the atmospheric layers where the 670.7 nm line forms. Since
non-LTE corrections vary in sign and size when spanning a large range of
stellar parameters (being larger for cooler stars), ignoring these
corrections may clearly affect any interpretation of the A(Li)-plateau and
its possible slope with effective temperature and metallicity. Therefore,
NLTE corrections were computed with the interpolation code made available
by Carlsson et al. (1994) and applied to our LTE A(Li) abundances. We note
that no NLTE correction could be derived for those stars with A(Li)
abundances smaller than 0.6, because the interpolation code works on a
given range of input parameters (for instance, in the case of the Li
abundance, the range is A(Li) = 0.6-4.2). Furthermore, despite these ranges
of input parameters, not all the combinations are covered in the table which
contains the tabulated coefficients from which the NLTE corrections are
computed. In our sample, we had only few of these cases, namely two stars
had a metallicity lower than the minimum threshold (-3.0), and another
couple of stars had their effective temperatures slightly higher of the
maximum threshold (by 4 and 23 K respectively). For these objects, we
rounded off their parameters to the nearest allowed value, and computed
the NLTE correction.
![]() |
Figure 11:
Comparison between our final set of effective temperatures
![]() |
Open with DEXTER |
Table 6:
Final Li abundances as derived by using different sets of
.
The whole table is available on-line.
The resulting LTE Li abundances (both LTE and NLTE are given only for our
final set of
) are presented in Table 6, where we give also other
relevant parameters like the metallicity and its 1
uncertainty as
derived from a critical analysis of the literature.
Three different sets of lithium abundances are reported for each target,
depending on which set of
was used (cf. Table 5). For the A(Li)
values (which were derived from our final set of temperatures, i.e.
(2))
we also give the associated
1
error
(Cols. 9, 10). Similar uncertainties apply also to the lithium abundances
listed under the last two columns of the table.
Our sample was assembled from the literature, following some selection criteria on the quality of the analyses, i.e. high resolution and high S/N. Since the Li I line falls in a very clean spectral region, with very few neighboring absorption lines, the placement of the continuum is usually quite accurate (on the order of 1-2%) if the data quality is high. This uncertainty is usually included in the uncertainty associated to the equivalent width measurement.
For the latter, because of the generally quite satisfactory agreement
between different literature sources (on a given target, cf. Fig. 1) we
decided to use as our final EW the arithmetic mean of all the measurements,
and take the dispersion around the mean as the uncertainty on each
measurement. When only one measurement was available, the associated
uncertainty is the error quoted in the original work. Table 6 reports both
the error on each EW and the corresponding 1
uncertainty on the Li
abundance. Except for few cases, the latter are well below 0.1 dex. The
uncertainty due to a
3% error in the log gf value is
0.013 dex.
Lithium abundances are known to be very sensitive to the chosen effective
temperature, but their dependence on the other parameters, i.e. gravity,
metallicity, and microturbulence is negligible. Common uncertainties on
log g, [Fe/H] and
(
0.25 dex,
0.15 dex, and
0.3 km s-1
respectively) affect the final Li abundances by at most 0.005 dex, 0.015 dex,
and 0.003 dex. When summed under quadrature, the resulting uncertainty is
around 0.017 dex only.
On the contrary, the dependence of Li abundances on the effective temperature
is much stronger. An uncertainty of 70 K in
(commonly quoted as
a reasonable uncertainty on this parameter) translates into a
0.05 dex
on the lithium abundance. For this work, we have considered only the
uncertainties associated to the photometric indices (b-y) (generally quoted
to be around 0.008 mag, cf. Nissen et al. 2002) and
from which we have
derived our reddening estimates (generally quoted to be around 0.011 mag).
When summed under quadrature, this gives us an average uncertainty on our
effective temperatures of
75 K, which corresponds to
0.054 dex in
A(Li).
Combining all the uncertainties together, we find that depending on the
(EW) error on A(Li) our best achievable accuracy is 0.06 dex.
In the worst cases it could be as high as 0.15 dex, but one should notice
that for all the stars for which a Li abundance uncertainty larger than
0.1 dex has been derived, the equivalent width of the Li I line is always
quite small (in the 5-15 mÅ range) with a very significant 1
EW
quoted error. Although we do not have the observed spectra available for
further checking, this indicates that S/N ratios on the order of 100 are
probably too low for accurate measurements of weak Li lines. If one were
to exclude those stars with very small equivalent widths (and large quoted
uncertainties) then our final (individual) accuracies range between 0.06
and 0.1 dex.
Last but not least, one should not forget that our abundances were
derived based on Kurucz non-overshooting model atmospheres and that most
of the current Li analyses are carried out under LTE assumptions (with NLTE
corrections applied to them) and with one-dimensional model atmospheres.
Choosing a different treatment of the convective motions in the model
atmospheres (i.e. choosing the so-called Kurucz overshooting models) has
the effect of deriving slightly higher Li abundances (by 0.08 dex)
but both sets of models carry similar uncertainties which are difficult to
quantitatively assess. NLTE corrections also carry their own uncertainty,
but this is small and well within the average abundance errors, according
to Carlsson et al. (1994). Finally, although Li abundances determined using
3D-hydrodynamical model atmospheres and corrected for NLTE effects differ
from 1D NLTE Li abundances by less than 0.1 dex (0.05 dex for the few stars
that have been investigated so far, cf. Asplund et al. 2003; and Barklem et al. 2003), these same authors warn about possible dependences on temperature
and metallicity, that could clearly affect any discussion on the existence
or lack thereof of a slope in the A(Li)-plateau with
and/or
metallicity. Hence, no conclusion can be final, until 3D NLTE effects on
Li abundances are mapped on a larger stellar parameters space.
Since one of the main focuses of this work is a critical assessment of how
accurately Li abundances in halo stars can be determined and the Li plateau
can be characterised (via its width, spread and slopes with
and
metallicity), comparisons between our derived Li abundances and previous
analyses are not very significant, especially since our analysis is not
based on newly observed spectra. Our work does not aim at showing that,
with our consistently determined temperature scale, we can now better
describe the properties of the A(Li)-plateau. On the contrary, our analysis
has so far pointed out that although this is clearly a must, even with such
a careful determination of the temperature scale, many uncertainties remain
especially for large samples.
However, since the opposite findings by some earlier analyses (e.g. R96, R99, and BM97) have indeed been among the initial triggers of this work, we think it is useful to further comment on few points.
First of all, we note that most of the discrepancies originally present among
some stars common to R96, R99, and BM97 (which some of the earlier claims for
dispersion and/or slopes may have originated from) can be fully explained
by differences in the stellar parameters adopted during the various analyses.
Table 7 summarises these comparisons for some objects common to these works,
with Cols. 2 to 5 reporting the differences in
, metallicity, EW (if
any), and the A(Li) abundance as reported by the original investigators
(always given as "R96-BM97''). When, for a given target, a double entry
is present, this second row of values corresponds to "R99-BM97''. The
last column of the table reports the remaining difference in the final A(Li)
value, after having taken into account the differences in
, [Fe/H], and EW.
Table 7:
Discrepant cases in the original works of R96, R99, and BM97.
Columns identified by a symbol represent the difference on that
given quantity (
, [Fe/H], EW, and lithium abundance) between R96 and
BM97 (when a target has a second entry, this refers to R99 - BM97).
This is indeed very positive, but in general we find a little worrysome the
revisions made by R99 to some of their measured equivalent widths, compared
to their own previous analysis from 1996. Although both analyses are based
on high quality, high S/N spectra, EW measurements of the Li I line for the
same star differ by as much as 6 mÅ which for the
range covered by
the R99 sample correspond to 0.06 dex in lithium abundance.
Also, the effective temperatures of the 10 stars common to the analyses of R99, BM97, and this work span respectively
220 K, 358 K, 520 K. In other words, the data-sample that Ryan et al.
carefully chose to span a very narrow range of effective temperatures (and
that indeed did so on their
scale), it is found to cover a much larger
interval when the temperature is derived following different prescriptions.
And all three analyses used self-consistent methods!
Another example stressing the weakness of the temperature scale issue
comes from the comparison of our results with the recent work by Meléndez &
Ramírez (2004). The latter have carried out a study similar to ours, in
which 62 halo dwarfs (of which, in the end, 41 were used to discuss the mean
Li abundance of the plateau) were analysed based on a newly derived and
improved IRFM-based
scale (Ramírez & Meléndez 2005a,b) and
using Li equivalent width measurements available in the literature.
A quick check between the effective temperatures reported in their Table 1
(Meléndez & Ramírez 2004) and our values, for the 32 stars we have
in common, shows both an offset (their temperature scale is hotter) and a
larger temperature interval spanned (737 K versus 583 K, the latter from
our analysis). The offset implies that their mean Li-plateau value will be
higher than ours. The fact that their temperature scale
is hotter than ours is not in contradiction with what is shown in the bottom
panel of Fig. 11, where a much larger sample of stars is plotted.
Still related to the determination of an accurate temperature scale, and as
already seen in Sect. 5, reddening plays an important role. For instance,
R99 noted that their reddening estimates based on two different methods
(Strömgren photometry and reddening maps) showed a clear discrepancy of
about 0.02 mag, with the Strömgren-based E(b-y) values being higher despite
the expected relationship
.
The solution chosen by
these authors was to give higher weight to the maps-based reddening values
based on the consideration that their targets were bright, hence a low
intrinsic reddening might be expected. Therefore, they systematically lowered
all the Strömgren-based E(b-y) values by 0.02 mag before averaging the two
methods. As we do not know the final answer either, we cannot say if this is
a good solution or not. At least for those objects we have in common these stars
fall well outside the inner 50-70 pc of the solar neighborhood
(they span distances up to 1 kpc), for which a low intrinsic reddening could be questionable.
Even more instructive is to evaluate the realistic situation in which one
may refold on selecting reddening values from different sources, such as
InfraRed dust maps and/or from the literature, because one method alone
cannot be applied homogeneously to the entire sample under investigation.
Keeping in mind that the exact details of such situation are very difficult
to foresee, hence to reproduce, we tried to mimick such case by plotting a
mix of lithium abundances:
Fig. 12 shows A(Li)
versus
, where
the lithium abundances were derived by selecting
(2) and
(S98) for
the clean sample (half-half), and
(H97) and
(mix) for the
and ubvy samples respectively. This likely represents an
extreme case, but it certainly gives an idea of what effect could be expected.
Also, please note that the A(Li) values corresponding to
(mix) were
computed only for those stars for which this solution was applied (cf. Sects. 5.1
and 5.2).
In summary, despite the seeming convergence at least on the
absence of dispersion, the finding of discordant results is not surprising
if some (or all) of the above-mentioned points are kept
in mind. At the moment the only claim for a tilted A(Li)
-plateau
is with metallicity, but most "metallicities''
quoted in the literature are still derived from neutral iron lines, the
formation of which is subject to NLTE conditions. Furthermore, R99 (as well
as our work) have used metallicity values extracted from a careful inspection
of the literature, which does not guarantee the homogeneity required for
discussing the A(Li) vs. [Fe/H] trend. Since the whole discussion of
a possible slope of the lithium plateau with metallicity is centred on a
small dependence, we encourage future analyses of Li abundances to determine
the metallicities spectroscopically, possibly from Fe II lines (insensitive
to NLTE effects).
![]() |
Figure 12:
Lithium abundances as determined by using the sets of temperatures
![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
A first visual comparison of the different A(Li)
abundances we
have derived for each star is presented in Fig. 13, where the first
three (from the top) panels show how differently the
plateau may appear when the lithium abundance has been derived from a
different set of temperatures. The bottom pan is shown only for comparison
purposes, and represents the A(Li)
abundances as derived from
photometry that has been dereddened using the BH H I maps.
Figure 13 shows how claims of dispersion or lack thereof
from the same data-sample are perfectly plausible
(see also our discussion in Sect. 7).
![]() |
Figure 13: Lithium abundances computed with four different sets of temperatures as derived by using different reddening estimates (as reported in the upper right corner of each pan). Upside-down open triangles represent abundance upper limits. |
Open with DEXTER |
![]() |
Figure 14:
Final lithium abundances with associated 1![]() ![]() |
Open with DEXTER |
In addition, since our complete sample is a non-homogeneous sample
(because of the compromises made on the derivation of the effective temperature
for some of the stars), any discussion of the width and slope of the
A(Li)-plateau requires the separation of the complete sample into
the clean, ,
and ubvy sub-samples. In
Figs. 14 and 15,
we plot the A(Li)
values vs. our
scale (i.e.
(2)) and [Fe/H] respectively, with
the clean sample always plotted in the top panel and the
and ubvy samples in the lower panel. Moreover, open
symbols refer to objects with
in
Fig. 14 and
5700 K in
Fig. 15, and filled symbols
represent respectively stars with [Fe/H] >-1.5 and
< 5700 K.
Upside-down triangles always identify abundance upper limits.
In the rest of the paper, we give our results
regarding the characteristics of the plateau for the clean
sample on one hand, and for the complete (i.e., clean +
+ ubvy) sample on the other hand.
Additionally, we will consider as "plateau stars'' those with
K and
.
The metallicity limit is taken in order to avoid any contamination
by lithium production from the various possible stellar sources
(e.g., Travaglio et al. 2001).
The cutoff in
is chosen for comparison reasons with
previous analyses in the literature. However for the purposes of
constraining stellar evolution models we will also discuss the cases
where the lower limit in
is increased to 6000 K
in order to avoid proto-stellar lithium destruction (in the case
of the dwarfs, see Sect. 10) or dilution at the very beginning of
the first dredge-up phase (in the case of the post-turnoff stars,
see Sect. 11).
Our results will be given considering the stars
with lithium abundance determinations only, the case of stars
with upper limits being considered separately in
Sect. 13
.
![]() |
Figure 15:
Final lithium abundances with associated 1![]() ![]() |
Open with DEXTER |
Considering the stars with [Fe/H]
and in the case of the most relaxed
lower limit (5700 K, cf.
Figs. 14 and 15) we obtain
For the stars with
6000 K, we
find a mean lithium abundance of
Similar conclusions on dispersion can be drawn if no
lower limit on [Fe/H] (in our specific case this was set to -2.1) is
assumed in the
- (b-y)0 calibration (see Sect. 5 for more details).
The higher effective temperatures thus derived and shown in Fig. 3 would have
the only effect of increasing our mean A(Li) plateau values by
0.025-0.035 dex.
We would like to note that the 6 stars that could be
moved to the ubvy sample (if we were to consider a 3 digit precision
on ,
cf. end of Sect. 5.2) are not relevant in the final
discussion of the Li plateau spread as they are cool stars that have
their lithium already depleted (5 of them) or a normal plateau lithium
abundance (1 star with A(Li)
).
The absence of intrinsic dispersion that we get is in agreement with the results of Molaro et al. (1995), Spite et al. (1996), Bonifacio & Molaro (1997), Ryan et al. (1999) and Meléndez & Ramírez (2004).
Furthermore, because all these works (ours included) find quite consistent
A(Li) plateau values (note that the apparently higher plateau value found
by Meléndez & Ramírez (2004, A(Li)
)
should be corrected
downwards by 0.08 dex because of the different Kurucz model atmospheres
employed in their and our study, making it then in closer agreement with our
findings), it seems that in general the relatively low lithium abundance
(when compared to the CMB+SBBN result) seen in metal poor halo stars is a
very robust result. Assuming the correctness of
the CMB constraint on the value of the baryon-to-photon ratio we are
then left with the conclusion that the Li abundance seen at the surface of
halo stars is not the pristine one, but that these stars have undergone
surface lithium depletion at some point during their evolution.
Let us try now to look for some constraints on the depletion mechanism(s).
Intrinsic absolute magnitudes MV are derived from the mV and the
parallaxes given in the Hipparcos catalogue.
We determine the bolometric corrections BC by using the relations between
BC and V-I (these quantities being also taken from the Hipparcos catalogue)
given by Lejeune et al. (1998) and which
are metallicity-dependent. We use the values of [Fe/H] derived in our final
analysis (Table 6, Col. 2). We first iterate using the log g values
available in the literature, and finally attribute to our stars the log g
value derived from their position in the HR diagram
(Table 9, Col. 10).
Finally, we compute the stellar luminosity and the associated error from
the one sigma error on the parallax. All the relevant quantities are given
in Tables 8 and 9.
Table 8: Quantities extracted from the HIPPARCOS catalogue for our sample stars ((B-V) and E(B-V) are given in Table 4). The whole table is available on-line.
Table 9: Characteristics and evolutionary status of the sample stars. The whole table is available on-line.
The resulting HR diagrams are shown in Fig. 16
for our complete sample split in four metallicity bins.
The evolutionary status of each star has been determined on the basis
of these HR diagrams, and is given in Table 9.
Each star has been assigned to one of the following sub-classes:
5 identifies the dwarfs (i.e., main-sequence stars), 4.5 stands for the
turnoff stars, 4 for the subgiants (i.e., stars crossing the Hertzsprung gap),
3.5 for the stars at the base of the RGB, and 3 for the stars on the RGB.
"Post-main sequence stars'' can be easily located on the HRD (see
Fig. 16) as those stars with Log(
)
higher
than
0.4. In this luminosity range, we identify as turnoff and subgiant
stars those with
5600 K, whereas cooler stars are
classified as "base RGB'' or RGB stars depending on their luminosity. This
limit in effective temperature is chosen because it corresponds to the
approximate value where lithium dilution is expected to occur at the very
beginning of the first dredge-up (i.e., Deliyannis et al. 1990; Charbonnel
1995). The distribution in these sub-classes is given in Table 10. For most
of our objects, the classification was performed unambiguously. Some stars
(21 in total) were first classified as "unsure'', the uncertainty being due
to the large errorbar on the derived luminosity inherited from the error on
the Hipparcos parallax (see Fig. 16). However, by
cross-checking the luminosity and the corresponding gravity obtained from
Hipparcos data with the MV and gravity values quoted in the
literature, we were able to attribute a relatively certain evolutionary status
(given in the last row of
Table 10) to all these objects, thus keeping them
in our statistical analysis.
As a result, the clean sample (respectively the complete sample)
contains 49 (59) dwarfs, 5 (5) turnoff stars, 29 (31) subgiants,
5 (18) base RGB stars and 2 (5) RGB stars.
In Fig. 17
we plot A(Li)
vs.
for four metallicity bins
and indicate the evolutionary status of the stars by different symbols.
It is interesting to note that for [Fe/H]
-2, all the coolest stars
(i.e., with
5500 K) are actually evolved stars
(see also Fig. 18).
This fact must be taken into account when comparing stellar evolution
models and observations.
In the [Fe/H] range between -2 and -1.5, the decrease of lithium
relative to the plateau dwarf stars appears at
5600-5700 K.
![]() |
Figure 16:
HR diagram for our complete sample stars for separate metallicity bins
(using
![]() |
Open with DEXTER |
Table 10: Distribution of the evolutionary status in our different sub-classes. The last line indicates the status that we could attribute to the stars which were first labelled as "unsure'' (see the text for more details).
In their careful analysis of the classical mechanisms that could alter the
surface Li abundance of halo stars at different phases of their evolution,
Deliyannis et al. (1990, hereafter D90) focused on the
separation of the halo population into pre- and post-turnoff groups as
an essential prerequisite to understand the Li observations.
At that time they considered a limited sample of halo stars. They used the
trigonometric parallaxes and V-colors from the Yale Parallax Catalogue
(van Altena et al. 1989) when available to determine MV; otherwise MV
was taken from the original observational papers.
We have 36 objects in common with D90's sample.
For 19 of them the same evolutionary status has been attributed in both D90
and our study. However 13 of the stars which were identified as "pre-turnoff''
objects by D90 appear to be more evolved stars on the basis of their more
precise and reliable Hipparcos parallaxes.
Conversely 4 out of the 7 stars which were claimed to be "post-turnoff''
objects in D90's study are dwarf stars.
Ryan & Deliyannis (1998) took advantage of the Hipparcos parallaxes (when available) in order to separate their sample of halo stars cooler than the plateau into dwarf and subgiant classes. The 14 stars that we have in common with their sample have been attributed the same evolutionary status in both analyses (theirs and ours).
Note that our definition of a subgiant is different and more strict, as far
as the effective temperature range is concerned, than the one adopted by
PSB93,
who studied the lithium abundances for 79 so-called halo subgiants.
PSB93 constructed their sample on the basis of ubvy photometry (more
precisely from their location in the c1 vs. (b-y) plane) from
several catalogs of metal-poor stars. Hipparcos parallaxes are now available
for almost all their objects, allowing a more precise determination of their
location in the HRD. In the present study we did not consider the most evolved
stars of PSB93's sample, i.e., those with
below 4800 K. These
objects indeed experience some extra-mixing beyond the first dredge-up and
lose then all the information about their initial lithium content (see
Charbonnel 1995; Weiss & Charbonnel 2004). They are thus of no help in
discussing the initial lithium content.
![]() |
Figure 17:
A(Li)
![]() ![]() |
Open with DEXTER |
![]() |
Figure 18:
A(Li)
![]() ![]() |
Open with DEXTER |
We now concentrate our analysis on the dwarf stars. Their lithium abundance is plotted in Fig. 18.
Table 11: Number of dwarfs in each subsample.
Table 12:
Average value of A(Li)
,
average and standard deviation, and rms for the dwarf stars :
mean
(rms). We separate the clean and the complete sample
and give the results for 2 limits in effective temperature.
For the reasons given in Sect. 8.1, we identify as "plateau stars'' those with
and with
higher than 5700 or 6000 K. The
corresponding numbers of stars in each bin is given in
Table 11.
For the time being, we keep in our analysis
the possible binary stars (but see Sect. 12), but we eliminate the stars with
lithium upper limits, that will be discussed separetely in Sect. 13. We note
that in the metallicity and effective temperature ranges we have chosen to
define the plateau only two dwarf stars of the clean sample have Li
upper limits, namely HIP 72561 and 100682. With A(Li)
and
1.088 respectively (with
and -2.83,
= 6388 and 6362 K), both stars fall well below the plateau (see
Fig. 18).
In order to investigate the existence of trends in the
[A(Li)
,
] and [A(Li)
,
[Fe/H]] planes
we performed univariate fits by means of four estimators:
(1) The least-squares fit with errors in the independant variable only
(routine FIT of Press et al. 1982);
(2) The least-squares fit with errors in both variables
(routine FITEXY of Press et al. 1982);
(3) The BCES (bivariate correlated errors and intrinsic scatter) of
Akritas & Bershady (1996);
(4) BCES simulations bootstrap based on 10 000 samples.
All univariate fits in the A(Li)
versus
plane
were computed considering an error on
of both 75 and 100 K.
The average value of the lithium abundance is given in Table 12, together with the average, standard deviation and root mean square of the estimated observational errors for several subsets of data (in terms of metallicity and effective temperature intervals, as well as in terms of clean vs. complete samples) of the dwarf sample.
If we consider the entire metallicity range [Fe/H] 1.5,
the straight average value of the lithium abundance for the plateau dwarf stars
with
5700 K is
When we restrict our analysis to the dwarf stars with
6000 K, the mean value
increases to
We thus find no evidence of an intrinsic dispersion in the Li abundances along the plateau. As expected, this result which was already presented in Sect. 8.2, is now reinforced having eliminated the "pollution'' by evolved stars and focussed on the dwarf stars only.
As we already mentioned in Sect. 8.2, the assumption of a lower limit
on [Fe/H] in our
- (b-y)0 calibration does not affect our
conclusions on dispersion. The only effect when one does not consider this limit
is to increase the mean A(Li) value of the dwarf sample by
less than 0.02 dex for the entire metallicity range. If we consider
only our most metal-poor subsamples the increase in the mean A(Li) value
is of course slightly higher (from 0.02 up 0.08 dex depending on the
lowest limit on
used for the plateau definition).
The mean lithium values for the dwarf stars are slightly lower, although fully compatible within the quoted errors, than the mean values given in Sect. 8.2 for the entire sample of stars (i.e., that in which we did not discriminate the stars according to their evolutionary status). This point will be discussed at length in Sect. 11.
When we consider the dwarfs with
5700 K and
[Fe/H]
-1.5, we find a small slope in the A(Li)
-
plane of
This is very similar to the slope found by BM97 of
0.02/100 K, and slightly lower than those found by Thorburn (1994) of
0.034/100 K and by R96 of 0.0408/100 K. On the contrary,
Meléndez & Ramírez (2004) do not find any dependence of A(Li) on either
or [Fe/H].
In order to focus on the physical processes that may affect the surface
lithium abundance only when the stars are on the main sequence (in other
words, to avoid any pre-main sequence depletion; see e.g. D90, Richard et al.
2004) we made the same computations for the stars with
6000 K. In this case and for the whole [Fe/H]
-1.5
range, we get a slightly smaller slope, i.e.
![]() |
Figure 19:
A(Li)
![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
We note that if we split the [Fe/H] interval as done for example in
Fig. 18, we find a negative slope of -0.025 for the stars with [Fe/H]
-2.5, both for the clean and complete samples for the stars with
6000 K. Small but
positive slopes are detected when we look at the -2.5 < [Fe/H]
-2.0
and -2.0 < [Fe/H]
-1.5 intervals.
In Fig. 19 we plot A(Li)
vs. [Fe/H]
for the plateau dwarfs of our entire sample but we look for trends
up to [Fe/H] = -1.5 only.
The highest slope (
0.05 up to [Fe/H] = -1.5)
is obtained for the most extended sample,
i.e., for the dwarf stars of the complete sample with
5700 K, whereas when we consider only the dwarf stars of the clean
sample (still with
5700 K) the correlation becomes
slightly flatter (the slope is
0.02).
When we focus on the dwarf stars with
6000 K, the
dependence of A(Li) with metallicity
remains small (with a slope of
0.03 and 0.02 dex up to [Fe/H] = -1.5
for the complete and clean samples respectively).
In Fig. 19 we also show the values when
no lower limit on [Fe/H] is assumed in the
- (b-y)0 calibration.
A slightly lower slope would have been derived under this assumption.
Our results are intermediate between those of BM97 (-0.05 to 0.00, i.e., no slope) and those of R99 (+0.118) and Thorburn (1994; +0.13). We note that the preliminary results of Asplund et al. (2005) indicate a dependence of A(Li) on metallicity (the latter computed both as the abundance of iron and oxygen), characterised by a slope of 0.10.
We remind the reader that lithium synthesis by galactic cosmic rays (GCR) leads
indeed to an increase of its abundance with metallicity, but this contribution
safely can be considered small up to
(Molaro et al. 1997)
in agreement with our finding.
This can be estimated using the observed correlation between 9Be and metallicity
together with the Li/Be ratio predicted by GCR theory (Walker et al. 1985).
As discussed already in Sect. 7.2, the trends with metallicity here derived must be taken with caution, as we have gathered our metallicities from the literature. This does not guarantee the homogeneity required for an accurate study of the metallicity dependence. Additionally such an analysis would require homogeneous spectroscopic determinations from Fe II lines insensitive to NLTE effects rather than from neutral iron lines as it is the case for most of the currently available Li analyses (hence also for our values). For these reasons we will not continue the discussion nor derive any conclusions on the lithium-metallicity relation.
As already mentioned in Sect. 9.2, the current sample does not contain
cool (i.e., with
5500 K) dwarf stars
at low-metallicity (i.e., [Fe/H]
2).
This is the first surprise of our check on the evolutionary status
of our sample stars. It thus appears that for the cooler and more
metal-deficient dwarfs we have no direct clue to the lithium behavior.
We expect these objects to exhibit the same pattern as that seen in cool open
cluster and Pop I field stars as well as in our sample stars at higher
metallicity (bottom panels of Fig. 18),
although this remains to be confirmed observationally.
In the [Fe/H] range between -2 and -1.5, substantial lithium depletion
relative to the plateau value sets in at
5600-5700 K.
Between 5700 and 5100 K, the lithium depletion slope is
0.21 dex per 100 K.
![]() |
Figure 20:
A(Li)
![]() ![]() |
Open with DEXTER |
Table 13:
Straight average value of A(Li)
,
average and standard
deviation for the turnoff and subgiant stars:
mean
(rms). We separate the clean
and the complete samples and present the results for 2 limits in effective temperature.
As already mentioned in Sect. 10, we put a lower limit on the effective
temperature (5700 and 6000 K) for our analysis, this time in order to avoid
any lithium depletion due to the first dredge-up. We consider stars both at
the turnoff and on the subgiant branch as defined in Sect. 9.2.
Two stars are excluded from the analysis, namely HIP 81276 ([Fe/H] = -1.50,
= 6523 K) of the clean sample and HIP 83320
([Fe/H] = -2.56,
= 5984 K) of the ubvy sample,
because of their A(Li)
upper limits (of 1.705 and 1.287 respectively,
see Fig. 20).
The average value of A(Li)
,
the average and standard
deviation and the root mean square of the observational errors are given
in Table 13 for various metallicity and effective temperature intervals and
for the clean and complete samples.
A very surprising result emerges here for the first time:
Whatever the subsample we consider (i.e., clean or complete
sample,
5700 or 6000 K), the mean value of A(Li)
always appears to be higher for the subgiant stars than for the dwarfs,
except for the most metal-deficient objects (i.e., with [Fe/H]
2.5)
where the mean lithium abundance is very similar in both evolutionary statuses.
In the case of the evolved stars with
5700 K and
,
When we restrict our inspection to the stars with
6000 K and
[Fe/H]
-1.5, the mean value of A(Li)
for the evolved stars
is
Once again, we note that our conclusions are not affected by having assumed -2.1
as the lower limit on the metallicity
to be used in our colour-
calibration (see Sect. 5).
The only effect of not considering this limit would be a slight increase
of the Li mean value by
0.04-0.05 dex.
Although the numbers for the dwarf and evolved stars are fully compatible
within the quoted errors, the difference between the mean lithium values
of both populations is stricking.
As already discussed in Sect. 7.1, the dependence of the lithium abundance on
gravity is weaker than that on effective temperature. At
= 6000 K
and [Fe/H] = -1.5, the typical effect is at most +/-0.01 dex in A(Li)
for -/+1 dex in log g. We can thus see that even an error of 1 dex on the
attributed gravity (which is very unlikely in view of the good precision of
the Hipparcos parallax for most of the stars) cannot explain the difference
on the mean lithium abundance that we obtain between the dwarf and subgiant
stars. Neither can the dependence of our colour-
calibration on
log g. There is indeed a dependence for giant stars because (b-y) measures
the slope of the Paschen continuum, which in turn is affected by a change
downwards of the gravity (cf. RM05b). However, what we call "evolved'' stars
are objects that have just passed the turn-off or that are located on the
subgiant branch: they are not real giants, as defined in RM05b.
Our analysis reaveals another remarkable feature: Post-main sequence Population II stars appear to exhibit a non negligible Li dispersion. As can be seen for example in Fig. 20, this is already true at the turnoff and all along the evolution traced by the effective temperature.
For the subgiants near 6000 K, PSB93 found a small spread in the lithium abundance around a mean value of 2.1. We see from the previously quoted numbers that the dispersion is actually not negligible in our sample of slightly evolved stars, independently of the adopted limit in effective temperature (5700 or 6000 K).
![]() |
Figure 21:
Lithium abundance versus
![]() ![]() |
Open with DEXTER |
For the more evolved (i.e., cooler) stars,
the lithium abundance is expected to decrease due to the first dredge-up
at the effective temperature of 5700 K, this dilution episode being
completed around 5200 K (i.e., Deliyannis et al. 1990; Charbonnel 1995).
This is what we observe, even though this region is not very well sampled.
However if all the stars had left the main sequence with the same lithium
content, they should share a common lithium abundance after the first
dredge-up (remember that these stars have approximately the same initial
mass).
Instead, the evolved stars with
lower than
5500 K which have already undergone the first dredge-up dilution
exhibit a relatively large lithium dispersion.
The question of the lithium dispersion among metal-poor stars has
already been discussed in the context of globular cluster studies.
In Fig. 21 we show the data in the only two
globular clusters for which Li abundances have been reported in stars
down to the turnoff, namely M 92 and NGC 6397 (middle and lower panel
respectively). These are NLTE values, i.e. we corrected the LTE values
reported in the literature for NLTE corrections which were computed
following the same prescriptions we used for our field stars sample.
Unfortunately, the unavailability of Strömgren photometry for all
the globular cluster stars analysed in M 92 and NGC 6397 did not allow
us to re-derive their lithium abundances based on our temperature
scale.
For M 92 we rely on two analyses: i) Boesgaard et al. (1998, hereafter B98),
from which we report the Li values for both temperature scales (Carney 1983,
C83; and King 1993, K93) discussed in their original paper; ii) Bonifacio
(2002, hereafter B02) which is based on B98 equivalent widths and on the
temperature calibration by Bonifacio et al. (2002, hereafter BSP02).
For NGC 6397, the Li abundances are taken from Pasquini & Molaro (1996),
Castilho et al. (2000), Thévenin et al. (2001) and Bonifacio et al.
(2002, BSP02).
In Fig. 21 the globular cluster data are compared
with the NLTE Li abundances of our evolved field stars in the metallicity
range around that of the clusters (-2.7 [Fe/H]
-1.7,
black symbols in the upper panel).
B98 investigated the lithium abundance in seven stars near the turnoff of
the old and metal-poor cluster M 92 ([Fe/H] = -2.52, King et al. 1998).
They reported a dispersion of a factor of 2.6 for the subgiant stars
in a region around
5800 K or 5950 K (depending on the adopted
temperature scale, C83 or K93 respectively) as can be seen in the middle panel
of Fig. 21. Contrary to B98, B02 concluded that there
is no strong evidence for intrinsic dispersion in Li abundances among the M 92
stars, although a dispersion as large as 0.18 dex is possible. B02 actually
warned the reader that no definitive conclusion can be drawn on the intrinsic
dispersion in this stellar cluster on the basis of the currently available
spectra. Better observations with higher S/N ratios are still needed
for these very faint M 92 stars.
Although the data points are relatively scarce in this region, one
sees from Fig. 21 that some dispersion indeed exists
among field stars in the
range of B98's data.
The stars HIP 36430, 92775 and 102718 which have respectively
[Fe/H] = -2.10, -2.18 and -1.80,
= 5985,
5762 and 5832 K, and very similar
Log(
) =
,
and
,
present significant
differences in their NLTE-Li abundance:
,
and
.
This corresponds to the dispersion claimed by B98 in M 92.
In addition at approximatively the same effective temperature the slightly
more metal-poor star HIP 83320 ([Fe/H] = -2.56,
= 5984 K,
Log(
) = 0.77
1.17) shows only an upper limit for Li
(A(Li)
1.287)
.
B98 discussed the case of the halo subgiant BD + 23
3912 (HIP 99423) which
was found by King et al. (1996) to have a remarkably high lithium abundance
of
.
This star is also in our sample, with the following
parameters: [Fe/H] = -1.54,
= 5806 K, and
Log(
) =
.
We find a relatively high NLTE Li abundance of 2.620 for this object. Although its relatively high metallicity
makes this star not very relevant when compared to M 92 stars, we think that
it reinforces the case for dispersion.
![]() |
Figure 22: Duplicity and variability among our complete sample. Asterisks: confirmed single- or double-lined binaries (Latham et al. 2002; Carney et al. 2003, 1994). Open squares: stars marked as spectroscopic binaries in Bonifacio & Molaro (1997). Open triangles: suspected binaries (Latham et al. 2002). Crosses: binaries or stars in double/multiple system as reported in SIMBAD. Filled circles: single stars. Open circles: variable stars as reported in SIMBAD. |
Open with DEXTER |
The lithium behavior in NGC 6397 ([Fe/H] = -2.02, Thévenin et al. 2001)
is instead completely different: stars with
6000 K
share the same lithium abundance with essentially no intrinsic
scatter. No lithium abundance has been reported up
to now in this cluster for stars over the
5600 to 6000 K range.
Observations in this region would be useful.
Table 14:
Average value of A(Li)
and average deviation
for [Fe/H]
-1.5:
When all the stars are considered regardless of their possible
duplicity or variability (Col. a), when all the suspected and confirmed
binaries are omitted (Col. b), and when only the confirmed binaries are
excluded (Col. c). In cases (b) and (c) the stars reported as variable
in SIMBAD are not considered. The first and second lines contain the values
for the clean and the complete samples respectively.
After a careful inspection of the literature in order to identify single- and
double-lined binaries (we mainly used the extensive surveys and listings by
Latham et al. 2002; Carney et al. 1994, 2003) we then tested if any of our above
conclusions on the mean lithium abundance and on the dispersion may have been
affected by the inclusion of binary stars. This is especially important because
the literature-based EW values may have been derived at a time when the binary
nature of the star was still unknown. In order to get the most conservative
answer, we also checked the SIMBAD database for binary and variable stars, and
we included also the stars marked as "binaries'' in BM97. All the stars thus
identified are marked with different symbols in Table 6 (2nd column, see the
legend at the bottom of the version of the table provided on-line) and
the situation can be visualized in Fig. 22.
Over the whole metallicity and effective temperature ranges sampled in this
paper, 27 (10, 3) stars from the clean sample are confirmed binaries
(suspected binaries, variable stars). The corresponding numbers for the ubvy and
samples are 1 (0, 0) and 4 (1, 3).
We compute the average value of A(Li)
and the average
deviation of the observational errors for two cases: We first omit all
the suspected and confirmed binaries, and then we exclude only the confirmed
binaries. In both cases, the stars reported as variable in SIMBAD are not
considered.
The results are reported in Table 14 for the dwarf stars on one hand and for the turnoff and subgiant stars together on the other hand. In each case the first and second lines refer respectively to the values of the clean and complete samples. We also recall the values discussed previously, which were obtained regardless of the duplicity or variability of the stars (see Col. a).
First, we see that in all the considered cases, the absolute numbers for the mean lithium abundance are always slightly higher (although fully compatible within the errors) when binary and variable stars are excluded. Secondly, the dispersion increases but only slightly. This confirms the early finding by Molaro (1991) that known binaries do not exhibit lithium abundances significantly different than the other stars and that they do not introduce significant scatter into the plateau. Finally, the finding that plateau dwarf stars exhibit a lower lithium mean value (although fully compatible within the quoted errors) than their evolved counterparts is robust and resists the duplicity check.
Table 15:
Sample stars with Li upper limits. The stars with a * lie in the
ranges in effective temperature and metallicity chosen to delimitate the Li plateau.
The others are either cooler or more metal-rich. The
indicates the stars
for which a v sin i value higher than 4.5 km s-1 could
be derived by Ryan et al. (2002) or Ryan & Elliott (2004). See the text for more
details.
In all the above discussions we have quoted the values for the mean lithium abundance and dispersion as well as for the trends with effective temperature and metallicity obtained without taking into account the stars with Li upper limits. The main characteristics of these Li-depleted stars are summarized in Table 15. In the metallicity and effective temperature ranges we have chosen to delimit the plateau, we have two dwarfs (HIP 72561 and 100682) and two more evolved stars (HIP 81276 and 83320). The first 3 of these objects have relatively high effective temperatures and are either dwarfs or subgiants lying very close to the turnoff, while the last one is clearly crossing the Hertzsprung gap.
Several studies have been devoted to these so-called ultra-lithium-deficient
(hereafter ULDs) halo stars.
Understanding their nature is crucial.
They should indeed be excluded from the investigations on the
primordial Li abundance if they belong to a very special class of objects.
However if they appear to be the extreme representatives of a process that
has affected all the stars along the plateau, they would be precious clues
on Li depletion in halo stars.
Up to now, no mechanism has been unambiguously identified as responsible
for the ULD phenomenon.
Norris et al. (1997) looked at the abundances of many elements in the
spectra of HIP 72561, 83320 and 100682 and found no common abundance
anomaly that could be associated with the Li deficiency (see also
Ryan et al. 1998; and Ryan & Elliot 2004). HIP 100682
presents an overabundance of the heavy-neutron capture elements (a feature
which is shared by 25
of halo objects with similar [Fe/H])
whereas none of the other two ULDs shows such an enrichment.
Hobbs et al. (1991) and Thorburn (1992) suggested that HIP 83320 and 100682 might be the progeny of blue stragglers that had depleted their Li to undetectable levels earlier on and are now evolving redwards, through the temperature range of the plateau. However Thorburn noticed that there are too many ULDs compared to the number of known halo blue stragglers.
The possible role of duplicity in producing the observed Li depletion in the ULDs via mass transfer has been discussed several times in the literature. Ryan et al. (2001a) proposed that the same mass transfer processes responsible for the field blue stragglers should also produce sub-turnoff-mass objects which would be indistinguishable from normal stars except for their Li abundance. This could in principle lead to Li depletion and Ba enhancement as observed in HIP 100682 if the donor companion was an AGB star. On the other hand, there also exist the possibility that the mass transfer from an RGB or an AGB star is not always accompanied by other abundance anomalies. HIP 72561 and 81276 appear to be multiple systems and HIP 100682 is a suspected binary (Carney et al. 1994, 2003; Norris et al. 1997). HIP 83320 however is a single star.
Ryan et al. (2002, hereafter R02) looked in detail at a sample of
18 halo main-sequence turnoff stars including 4 ULDs.
They discovered that 3 out of these 4 Li-depleted stars (but none of the
Li-normal stars) exhibited unusually broad absorbtion lines that could be
attributed to rotational broadening. Among these objects are HIP 72561
and 81276 for which values of
and
km s-1were respectively inferred for v sini (Note that the "Li-normal''
stars have undetectable rotation, generally below 3 km s-1).
The 3 broadened objects of R02's sample all have relatively high effective
temperature and they lie close to the turnoff. Furthermore, they are all
confirmed binaries with orbits that are not tidally synchronised. These
complementary features lead R02
to draw a connection between Li depletion, rapid rotation
and mass and angular momentum transfer from a companion, now the white dwarf
remnant of a star initially more massive than the present plateau stars, as
confirmed by the projected companion masses inferred by Carney et al. (2001).
Ryan & Elliott (2004) looked again at line-broadening in ULDs and showed
that 5 out of 8 have rotation velocities in excess of 4 km s-1 and
that 4 out of 5 are confirmed binaries. The invoked mass transfer thus has
not resulted in the merger of the components.
In the scenario proposed by Ryan and collaborators, the Li depletion could be due either to mixing triggered during the accretion event, or to Li deficiency of the donor only or of both companions prior to the mass transfer. If additional evidence could be found in favor of this scenario for the formation of ULDs, it would become clear that these objects are not useful to infer the primordial abundance nor to constrain the classical depletion mechanism(s). However, although this explanation is appealing and may work for the line-broadened and hot ULD stars, the unbroadened and/or single ULDs still await a plausible interpretation. As we shall discuss in more detail later on, the fact that all these stars lie on or originate from the hottest side of the plateau may reveal an alternative explanation.
In our sample some other stars that are cooler and/or more metal-rich than the plateau also have Li upper limits, namely the dwarfs HIP 67655 and 67863 and the post-turnoff stars HIP 55022 and 60719 (see Table 15).
HIP 55022 is a relatively hot and metal-rich subgiant that lies very close
to the turnoff. It is a confirmed binary, and could be the high-metallicity
counterpart of the ULDs previously described.
Ryan & Elliott (2004) give a v sini value of
km s-1 for this star.
All the other three objects are single stars. HIP 67655 and 67863 are cool dwarf stars that lie in the effective temperature region where the lithium depletion relative to the plateau is strong (see Sect. 10.5). Therefore, the non-detection of Li in these two stars does not appear to be a real abnormalty.
HIP 60719 is a single star at the base of the RGB that is undergoing Li dilution. Although this could explain its Li upper limit, this star clearly lies below the other evolved stars of our sample indicating that it has undergone a more severe lithium depletion. It is slightly more evolved than the single subgiant HIP 83320 which was discussed previously. For none of these two stars can mass transfer be advocated to explain the Li behavior. Note that HIP 60719 and 83320 additionally contribute to the Li dispersion observed in evolved stars and discussed in Sect. 11.2.
On the other extreme, some stars of our sample present relatively high lithium abundances, namely HIP 86694 and 99423. Both of them are post-turnoff stars. We have tried to identify in our sample other stars that share the same characteristics as these objects, but differ in their Li abundance.
No real pair star could be attributed to HIP 86694. This star
appears in several independent studies which all confirm its high Li content.
A(Li)
values between 2.481 and 2.554 are obtained on the basis
of the extreme literature EW determinations.
HIP 99423 has the same characteristics as HIP 73385 and 114962
(Teff of 5806, 5831, 5883 K; [Fe/H] of -1.54, -1.62, -1.52;
Log
of 0.482, 0.552, 0.554) but a significantly higher
Li abundance (2.620, 2.383 and 2.245). Although HIP 99423 seems to be
a single star, its "pairs'' are suspected or confirmed binaries.
We summarise our findings concerning the stars that present Li abnormalities, i.e., Li deficiencies or high Li values:
Although we do not aim to homogeneize the 6Li/7Li data available in
the literature, we need to comment on the very few halo stars in which 6Li
has been detected. This isotope is extremely difficult to observe, being only
a weak component, blending with the much stronger 7Li doublet at 670.7 nm.
The isotopic separation is 0.16 Å only. So far it has been
detected in only a few stars which all belong to our sample, namely
HIP 8572 (G271-162), HIP 48152 (HD 84937) and HIP 96115 (BD + 263578).
For all the other halo stars that have been looked at, only upper limits could
be derived (see Smith et al. 1993, 1998, 2001; Hobbs & Thorburn 1994, 1997;
Cayrel et al. 1999a; Hobbs et al. 1999; Nissen et al. 2000).
In Table 16 we recall the main characteristics of our sample stars for which
positive detections of 6Li have been reported and published. Note that the
3 stars under scrutiny are at the turnoff (i.e., with relatively high
)
and have relatively low [Fe/H] values
.
They have no additional peculiarities except for 6Li when compared to other
stars in which this isotope has been looked for but not detected (see e.g.
Smith et al. 1998).
Table 16:
Sample stars with published 6Li detection.
The data are from Smith et al. (1998, sln98), Cayrel et al. (1999a, c99) and
Nissen et al. (2000, n00). Columns 2 and 3 give the
(2) and [Fe/H]
values that we attributed to the stars while Cols. 4 and 5 display the
values quoted in the original papers. For completeness, we note that a
re-analysis by Asplund et al. (2005) of HIP 8572 and HIP 96115 from very
high resolution and S/N UVES (Dekker et al. 1999) spectra finds a slightly
less than 2
detection for HIP 8572 and no 6Li detection for HIP
96115.
Additionally, recent results from Asplund et al. (2005) indicate
nine new 6Li detections, in stars spanning a range in metallicity from
down to -2.7. Four out of these nine objects belong to our
sample and again turn out to be stars with relatively high effective
temperatures, close to the turnoff, and with relatively low metallicity.
6Li, together with 9Be, and 10B (the exact contribution from
-spallation in supernovae to 11B being still under debate), is
believed to originate primarily (plus some extra contribution by
reactions and other stellar sources, the latter important
at higher metallicites) from spallation reactions in the interstellar
medium between cosmic-ray (CR)
-particles and protons and heavy
nuclei like CNO (Reeves et al. 1973) . Basic considerations about the
production rate and environment of these light nuclides predict a quadratic
slope in the logarithmic plane ((light nuclide, e.g. 6Li),metallicity),
whereas the almost linear slope derived from spectroscopic analyses of Be
and B abundances in stars of the Galactic halo (e.g. Boesgaard et al. 1999
for Be, Primas et al. 1999 for B) seems instead
to indicate a primary (instead of secondary) origin, which in turn requires
the need for a production mechanism independent of the metallicity in the
ISM. One can generally refer to
the two abovementioned scenarios as to classical vs. reverse
spallation reactions, with the latter having several implementations based on
different assumptions.
Despite the small existing number of 6Li detections, this light nuclide seems to have a different history, showing a much flatter evolution (i.e. most of GCR-based models underproduce the amount of 6Li observed in halo stars).
Because 6Li is the most fragile isotope to proton destruction, its presence in the atmosphere of Pop II stars is usually considered as a very severe limit on the amount of 7Li depletion. This argument is very often used in favor of the Li plateau abundance being the primordial value (Brown & Schramm 1988), a conclusion that is challenged now by the CMB constraint. In this context several questions remain open. First, what is the pre-stellar value of 6Li/7Li (Probably not the one we observe now, because of the different sensitivity of both isotopes to nuclear destruction)? How can we explain the presence of 6Li in some of the plateau stars? Why do some turnoff stars exhibit some 6Li in their spectra while the majority does not?
One possibility calls for some 6Li production at the stellar surface by either stellar flares (Deliyannis & Malaney 1995, but see the criticism by Lemoine et al. 1997) or Galactic cosmic rays (Lambert 1995). The other possibility is more related to the questions about the origin of 6Li in metal-poor stars and about a possible scatter in the 6Li/7Li ratio in the ISM at a given metallicity. The type of cosmic ray sources and the production mechanims operating in the early, forming galaxy are still very controversial (see e.g. Vangioni-Flam et al. 1999). Some models for the formation of light elements by cosmic ray processes predict a scatter of one order of magnitude in the abundances of 6Li, Be and B relative to Fe. This is the case of the bimodal superbubble model by Parizot & Drury (1999) and of the supernovae-driven chemical evolution model for the Galactic halo by Suzuki et al. (1999). The superbubble scenario is actually challenged by several observations, for example the fact that there seems to exist no association of core collapse SN with superbubbles but rather to HII regions which reflect the metallicity of the ambiant interstellar medium rather than that of the SN (see the discussion in Prantzos 2004).
Suzuki & Inoue (2004) recently presented an additional mechanism for cosmic ray production of 6Li by virialisation shocks during hierarchical structure formation of the Galaxy. In their model, cosmic rays accelerated by this source dominate the production of 6Li (compared to SNe as the main source of acceleration), without co-producing Be and B. This seems to better account for the observed 6Li data in halo stars, and in particular for the "plateau'' of 6Li/H reported by Asplund et al. (2005), than the SN-driven cosmic ray scenarios. Large 6Li scatter should be a natural consequence of this process, a higher initial content being expected in stars belonging to the inner halo compared to those of the outer halo. This is an attractive scenario, although very difficult to estimate quantitatively due to the many uncertainties on the physics of structure formation and on the energetics of the implied cosmic rays.
Clearly the detection of 6Li in low metallicity halo stars constitutes a challenge to cosmic ray spallation and/or stellar mixing along the plateau. Alternatively it may be the case that the lithium isotope ratio probes Physics beyond the Standard Model. Indeed Jedamzik (2004) found that the decay of a supersymetric particle like the gravitino or neutralino around 103 s after the Big Bang may at the same time yield to a reduction of the primordial 7Li/H by a factor 2-3 and produce 6Li to the magnitude observed in halo stars. However before we turn to such possibilities, the astrophysical uncertainties and solutions must be critically assessed.
Since the discovery of lithium in Pop II stars, not only the observers have been very active on the subject. Theoreticians, too, have taken advantage of the lithium diagnostic to probe the interior and the evolution of Pop II stars. Although these low-mass objects seem to be very simple at first sight, their Li behavior is paradoxal and has still not been fully understood. Part of the difficulty comes from the fact that nothing is known about the initial conditions of an important quantity like the angular momentum distribution. We do not aim here to discuss all the literature published on the subject (see the review by Pinsonneault et al. 2000 and Talon & Charbonnel 2004, for more recent references). We rather use the specific findings of our analysis to propose a synthesis, i.e., to extract from the theoretical debate some of the most adapted and promising directions that require further investigation.
The (incorrectly) so-called "standard'' stellar models have long been
very popular in the plateau debate, mainly because they predict
no variation of the surface Li abundance with time and consequently
no Li dispersion from star to star along the plateau.
The corresponding Li isochrones thus naturally support the conclusion that
the Spite plateau reflects the cosmological Li abundance.
This belief has been shaken by the recent CMB measurements.
In addition to and maybe even more fundamental
(from the stellar physicist point of view)
than the CMB result, the validity of these stellar models can be refuted
on the basis of some of their assumptions, the most critical of which being
the absence of transport of chemicals except in the convection zones.
This hypothesis denies the fundamental nature of a star which
is a gaseous mixture of elements with various atomic masses.
As a result the stellar gas cannot be in equilibrium as a whole until
each component reaches its own equilibrium via gravitational
separation and thermal diffusion (Eddington 1929).
Thanks to helioseismology it is now fully recognised that atomic diffusion
must be an integral part of stellar
evolution computations and in particular of the standard
model
of the Sun (i.e., Richard et al. 1996; Bahcall et al. 1997) and of low-mass
stars.
The relative unpopularity of the Pop II standard models rests mainly on the fact that pure atomic diffusion leads to a degree of surface Li depletion which increases with effective temperature along the plateau (Michaud et al. 1984; Deliyannis et al. 1990; Proffitt & Michaud 1991; Chaboyer & Demarque 1994; Vauclair & Charbonnel 1995; Salaris & Weiss 2001; Richard et al. 2002). As confirmed again in the present analysis, this feature is not observed, although some of the hottest dwarf stars with Li deficiency may exhibit the effects of atomic diffusion (see Sect. 13.1, and discussion below). This difficulty can be simply read as the signature of some macroscopic processes that minimize the effects of atomic diffusion which is always present in stellar interiors and cannot be arbitrarily turned off. Such a suggestion was already made in a more general context by Eddington (1929) who pointed out that some mixing must occur in the stellar radiative zone in order to prevent the gravitational settling and the thermal diffusion of heavy elements, the effects of which were not always observed. Consistently, such a process is also required in Pop II stars to counteract the settling of heavier elements in order to explain the close similarity of iron abundances in near turnoff, sub-giant and lower RGB stars in globular clusters, and to reproduce the observed morphologies of globular cluster color-magnitude diagrams (see VandenBerg et al. 2002).
Several processes have been invoked to counteract atomic diffusion in stellar interiors. In the case of Pop II stars the possible candidates that have been studied (in models including or not atomic diffusion) are : Rotation (Vauclair 1988; Chaboyer & Demarque 1994; Pinsonneault et al. 1991, 1999, 2002; Vauclair & Charbonnel 1995), stellar wind (Vauclair & Charbonnel 1995), interaction between meridional circulation and helium settling (Théado & Vauclair 2001). All these models have real difficulties, unless some adhoc assumptions are made, in reconciling a non-negligible Li destruction with both the flatness and the extremely small dispersion on the plateau as definitively needed in view of the difference between the WMAP constraint and the plateau value discussed in this paper. This does not mean of course that the physical processes invoked are not at work in Pop II stars but it indicates at least that their theoretical description is still incorrect or incomplete.
Richard et al. (2002, hereafter Ri02, 2004) re-investigated the case
of the Li plateau with a new generation of Pop II models that include
self-consistently all the effects of atomic diffusion in the presence of
weak turbulence whose nature is not postulated a priori.
They discuss in great detail the parametrization of turbulence that has
to be included in order to
reproduce the constancy of the Li plateau (see also Proffitt & Michaud 1991
and Vauclair & Charbonnel 1998 and references therein).
The mean Li value we derive in the present paper for the plateau dwarf stars
favors such a model with an intermediate efficiency of turbulence (actually
between the so-called T6.0 and T6.25 models of Ri02 ).
In addition Ri02 show that acceptable variations in the
turbulence (i.e., from no macroscopic motion to that needed to fit the
plateau) can lead to variations of the Li abundance as high as 0.5-0.6 dex
at the turnoff. As can be seen in their Figs. 14 and 16 the Li abundance at
that evolutionary stage is indeed a very sensitive function of the exact
position in the HR diagram and of the adopted turbulence.
The abundance variations and in particular the Li deficiencies that we find
in our data at or just after the turnoff as well as those found by B98
in M 92 are similar to those expected at this evolutionary stage by Ri02 in
the case of variations of turbulence from star to star.
However in Ri02's models these abundance variations are
theoretically erased by dilution in the subgiants when they reach
6000 K,
although the dispersion persists in the data at and below this effective
temperature.
This difficulty may be alleviated by assuming that some stars
have undergone even stronger turbulence that lead to stronger Li destruction.
As previously mentioned, Richard and collaborators do not postulate
the physical mechanism that causes the turbulence required by their models.
Our present results bring a very important piece to the puzzle.
We could show that the turnoff and more evolved
Pop II stars present a slightly higher Li mean value as well as
a larger Li dispersion than the
less massive dwarfs.
This points towards a mechanism, the efficiency of which changes as one
reaches the extreme blue edge of the Li plateau.
Such behavior corresponds to that of the generation and filtering of
internal gravity waves in Pop II stars. As shown by Talon & Charbonnel (2004)
gravity waves are indeed very efficient in dwarf stars along the
plateau up to an effective temperature of 6300 K.
There they dominate the transport of angular momentum
and should lead to a quasi-solid rotation state of the stellar radiative zones
on very short time-scales.
As a result the surface Li depletion is expected to be independant of the
initial angular momentum distribution in this range of effective temperatures.
This should alleviate the difficulty encountered by the classical
rotating models which predict that a range of initial angular momenta
generates a range of Li depletion and that the scatter increases
with the average Li depletion.
In more massive stars however the efficiency of the gravity waves
strongly decreases and
internal differential rotation is expected to be
maintained under the effect of meridional circulation and turbulence
induced by rotation
.
Consequently variations of the original angular momentum from star to star
would lead to more Li dispersion and to more frequent abnormalities
in the case of the most massive stars where gravity waves are not fully efficient.
In other words internal gravity waves are expected to dominate completely
the transport of angular momentum and should lead to higher Li homogeneity
in the less massive, rigid rotators than in the stars lying on or originating
from the hot side of the plateau.
The proper treatment of the effects of gravity waves together with those
of atomic diffusion, meridional circulation and shear turbulence has
now to be undertaken in stellar evolution models in order to test the real
and interactive consequences of all these complex mechanisms
(Talon & Charbonnel 2005).
In the first part we focused on the systematic uncertainties affecting the
determination of Li abundances. We explored in detail the temperature
scale issue and put special emphasis on reddening with the aim of deriving a
tool as consistent as possible for all our sample stars. In order to do so,
we chose to derive photometric temperatures using Strömgren uvby-
photometry which was the only one available for our entire sample. During
these steps, we identified one of the major drawbacks of such determinations,
namely an accurate estimate of the reddening affecting each of our stars.
We derived four sets of effective temperatures based on different assumptions
for the interstellar reddening excess values, E(b-y).
We tried to evaluate the effect of using reddening values taken from different
sources on the derived temperature scales and in turn on the derived A(Li)
abundance, showing that an unpredictable mix of different reddening sources could
be held responsible for opposite findings, on the same dataset, about the presence
of dispersion and/or slope on the Li plateau.
Finally, we selected as our best and final
scale the one (
(2)) that has been
derived from de-reddening all the stars, except those with negative E(b-y)values.
In order to keep as many stars as possible in our analysis we had to make some
compromises on the derivation of the effective temperature of some objects.
This contaminated our complete sample of
118 stars which was finally subdivided as follows:
1) the clean sample which contains 91 stars for which
the complete set of Strömgren photometric indices are available
and for which the Schuster & Nissen (1989) calibration
for the interstellar reddening excess is applicable;
2) the
sample which includes
20 stars for which the reddening E(b-y) value was derived from
averaging different sources of E(B-V) via the Crawford's formula;
3) the ubvy sample which contains 7 stars for which one of the
ubvy photometric indices was found to fall just slightly outside
the allowed intervals for the application of the Schuster & Nissen
calibrations.
We made several tests to quantify the influence of these compromises
on the statistical analysis and on our final assessments.
In order to guarantee the absence of spurious differences and
conclusions due to the use of different criteria in the determination
of the effective temperature,
all our results (summarised hereafter) regarding the characteristics
of the plateau were given for the clean sample on one hand and
for the complete (i.e., clean +
+ ubvy) sample on
the other hand.
This approach should provide as high accuracy and reliability as possible
for one of the largest sample yet studied.
We then derived the lithium abundances for the various subsamples using
the different sets of temperature. This was done using the arithmetic
mean of the equivalent widths and the 1
uncertainty of the 670.7 nm
line as reported in the literature from which we had assembled the sample.
The lithium abundance was first derived under LTE assumptions for
all the
scales studied here and then NLTE corrections were applied. With these
NLTE Li abundances we determined the mean Li value and dispersion
along the plateau for our sample as a whole.
In order to avoid any contamination by lithium production from various
stellar sources we then restricted our discussion only to those stars
with [Fe/H] lower than -1.5. In this metallicity range, we considered as
plateau stars those with an effective temperature higher than 5700 or 6000 K.
Stars with Li upper limits were excluded from the analysis and are discussed
separately later.
We found sligthly different results for the mean lithium abundance
and for the dispersion depending on the lowest limit on
and on the sample under consideration, i.e. clean or complete.
We note however that the pollution due to our compromises on the
derivation of the effective temperature of some of our objects has a
negligible impact on our conclusions. The effect of having assumed
a lower limit on [Fe/H] in the colour-
calibration
used here to derive our sets of effective temperatures is also negligible.
This assumption has only the effect of raising the A(Li) plateau
abundances by 0.03 dex, on average, making our conclusions quite robust.
In the case of the stars of the clean sample with
6000 K
we obtain
In the second part of the present paper we further pushed the constraints on lithium depletion in halo stars. Using our homogeneized data we looked at the Li plateau by considering the evolutionary status of each star. This could be done using the Hipparcos parallaxes which were available for almost all our initial sample stars. This step of the analysis proved to be crucial since a contamination exists from post-main sequence stars, which has to be removed in order to precisely determine the depletion factor along the plateau. Several conclusions could then be drawn.
Again the mean lithium abundance for the dwarf stars depends on
the lowest effective temperature chosen to delimit the plateau and
slightly varies when one considers the clean sample only or
the complete one.
The mean lithium plateau value for the dwarf stars of
the clean sample with
6000 K is
A very surprising result was found for the first time:
Whatever the subsample we considered, the mean value of A(Li)
always appears to be higher (although compatible within the errors)
for the subgiant stars than for the dwarfs, except for the most metal-deficient
objects (i.e., with
)
where the mean lithium abundance is very
similar in both evolutionary status.
The mean lithium value for the post-main sequence stars of
the clean sample with
6000 K is
We checked that the above conclusions do not change when we exclude the stars that belong to multiple systems or show variability. We could confirm that binary stars do not exhibit lithium abundances significantly different to their single counterparts.
We finished our close examination of our sample by looking at the stars that
present deviant Li abundances, i.e., the stars with strong Li deficiency
(the so-called ULDs) and those with an abnormally high Li content.
We found that all of them lie on or originate from the hot side of the plateau.
This agrees with our finding that the turnoff and subgiant stars
present a slightly higher Li mean value and dispersion than the dwarfs.
These results indicate that the post-main sequence halo stars experienced
a Li history slightly different from that of the less massive plateau dwarfs.
We suggest that such a behavior may be the signature of a transport process of the chemicals
and of angular momentum whose efficiency changes on the blue edge of the plateau.
This in agreement with the fact that most of the ULDs are presently rotating
faster than the Li-normal stars.
Our analysis provided thus some crucial clues to the internal processes
that may be involved in modifying the surface Li abundances in halo stars.
Since internal gravity waves coupled with rotation-induced mixing are
expected to lead to higher Li homogeneity with
in the plateau
stars than in the more massive stars lying on or originating from the
hot side of the plateau, such a model is favoured.
Although we excluded the ULDs stars from the analysis and focused on
the "normal-Li'' stars to derive the Li mean values and the trends
with
and metallicity, these objects should not be
excluded from hydrodynamical studies of the Li depletion mechanisms
that affect the Pop II stars.
Acknowledgements
We dedicate this paper to the two new little lithium boys Angel David and Martin Lou. F.P. is indebted to CNRS and the University Paul Sabatier in Toulouse for the Visiting Researcher contracts awarded. F.P. warmly thanks the Laboratoire d'Astrophysique de Toulouse et Tarbes, France for the warm hospitality received during her extended and frequent visits. C.C. thanks the french Programme National de Physique Stellaire and Programme National Galaxies for financial support. Part of this work has been carried out at the University of Washington (Seattle, USA) where both F.P. and C.C. participated in the INT Programme on Nucleosynthesis (April 2002). We used the SIMBAD data base operated at the CDS (Strasbourg, France).
Table 3: The data sample and its main characteristics, as found in the literature.
Table 4: Photometry and reddening excesses.
Table 5: Final temperatures as derived from different reddening corrections.
Table 6:
Final Li abundances as derived by using different sets of
.
Table 8: Quantities extracted from the HIPPARCOS catalogue for our sample stars ((B-V) and E(B-V) are given in Table 4).
Table 9: Characteristics and evolutionary status of the sample stars.