A&A 438, 251-256 (2005)
DOI: 10.1051/0004-6361:20052750
P. Beirão1 - N. C. Santos1,2 - G. Israelian3 - M. Mayor2
1 - Centro de Astronomia e Astrofísica da Universidade de Lisboa,
Observatório Astronómico de Lisboa, Tapada da Ajuda, 1349-018
Lisboa, Portugal
2 -
Observatoire de Genève, 51 ch. des
Maillettes, 1290 Sauverny, Switzerland
3 -
Instituto de Astrofísica de Canarias, 38200
La Laguna, Tenerife, Spain
Received 24 January 2005 / Accepted 31 March 2005
Abstract
We present Na, Mg and Al abundances in a set of 98 stars with known giant
planets, and in a comparison sample of 41 "single'' stars. The results
show that the [X/H] abundances (with X = Na, Mg and Al) are, on average, higher in
stars with giant planets, a result similar to the one found for iron.
However, we did not find any strong difference in the [X/Fe] ratios, for
a fixed [Fe/H], between the two samples of stars in the region where
the samples overlap. The data was used to study the Galactic
chemical evolution trends for Na, Mg and Al and to discuss
the possible influence of planets on this evolution. The results,
similar to those obtained by other authors, show that the [X/Fe]
ratios all decrease as a function of metallicity
up to solar values. While for Mg and Al this trend then becomes relatively
constant, for Na we find indications of an upturn up to [Fe/H] values
close to 0.25 dex. For metallicities above this value the [Na/Fe] becomes
constant.
Key words: stars: abundances - stars: fundamental parameters - stars: planetary systems: formation - galaxy: abundances - solar neighbourhood
The study of the chemical abundances of planet-host stars is providing important clues about the processes of planet formation and evolution. Soon after the discovery of the first exoplanets, it had been suggested that their host stars were particularly metal-rich compared to average field dwarfs (Gonzalez 1998). As new planets were added to the lists, this fact became more and more clear (e.g. Reid 2002; Santos et al. 2003; Gonzalez et al. 2001; Santos et al. 2004c,2001). Today, it is well accepted that the probability of finding a planet is a strongly increasing function of the stellar metallicity, at least for [Fe/H] values above solar (Santos et al. 2004a).
Besides the analysis of the global metal content of the planet-host stars, some efforts have also been made to analyze the chemical abundances of specific elements. These include light element abundances (e.g. Garcia Lopez & Perez de Taoro 1998; Israelian et al. 2003,2004; Santos et al. 2002; Deliyannis et al. 2000; Ryan 2000; Santos et al. 2004b; Gonzalez & Laws 2000), as well as other metals (e.g. Ecuvillon et al. 2004a; Smith et al. 2001; Gonzalez & Laws 2000; Santos et al. 2000; Sadakane et al. 2002; Takeda et al. 2001; Bodaghee et al. 2003; Ecuvillon et al. 2004b; Gonzalez et al. 2001). In general, these studies agree that planet host stars are globally metal-rich compared to single field dwarfs. With a few exceptions (e.g. Israelian et al. 2003; Gonzalez et al. 2001) no major unexpected trends were suggested.
In their study of Na, Mg and Al abundances in stars harboring giant planets, Gonzalez et al. (2001) have presented evidence that the ratios [Na/Fe], [Mg/Fe] and [Al/Fe] may be smaller, for a given [Fe/H], in stars with planets when compared to "single'' dwarfs. Unfortunately, Gonzalez et al. (2001) had to compare their Na, Mg and Al abundances for planet-hosts with those derived by other authors for field stars. The cause of the observed difference, not seen by other authors (Sadakane et al. 2002), may thus be simply due to the use of non-uniform samples.
To address this issue, we present in this paper a uniform analysis of the abundances of Na, Mg and Al in a sample of 98 stars with planets and a comparison sample of 41 "single'' field stars. Stellar parameters for both samples have been obtained from a very uniform spectroscopic analysis (Santos et al. 2004c). Furthermore, the abundances of the three elements studied were derived from the analysis of high-resolution spectra using the same line-lists for every star. The results show that planet host stars have, on average, [X/H] values (with X = Na, Mg and Al) above those found in field dwarfs. However, no major differences were found to exist between the [X/Fe] ratios in planet hosts and those found in our comparison sample stars for a given value of [Fe/H].
With the exception of HD 330075 (Pepe et al. 2004), the spectra used in this paper were those taken by Santos et al. (2004c) in their study of 98 planet-host stars and 41 comparison sample "single'' dwarfs. For a more thorough presentation of the samples see Santos et al. (2004c, and references therein). HD 219542 B was excluded from the planet-host star list as the presence of a planet around this star has been disproved (Desidera et al. 2004).
The stellar parameters used are also the same has listed by Santos et al. (2004c)
and Pepe et al. (2004)
. These were
derived in a very uniform way, using high resolution spectra and the same line-lists and
model atmospheres for all the stars.
Table 1: List of Na, Mg and Al lines and atomic parameters.
Abundances for the elements studied here were derived from the analysis of a list of spectral lines initially taken from the literature (Sadakane et al. 2002; Edvardsson et al. 1993; Chen et al. 2003; Smith et al. 2001; Feltzing & Gustafsson 1998; Gonzalez & Laws 2000; Gonzalez et al. 2001). A careful choice of these was done based on the analysis of the Kurucz Solar Flux Atlas (Kurucz et al. 1984) to include only those lines that are not blended in the solar spectrum. All the lines used have wavelengths between 5000 and 9000 Å.
To derive semi-empirical atomic
values for the lines, we used their Equivalent Widths
(EW) measured in the solar spectrum, and performed an
inverted solar analysis. The solar abundances for each element were taken
from Anders & Grevesse (1989), and the solar model was considered to
have
K,
dex,
km s-1 and
dex. The final list of lines with their atomic parameters
is presented in Table 1.
For each star, line EW were then measured in the stellar spectra using the IRAF "splot'' task within the echelle package. The analysis was done in Local Thermodynamic Equilibrium (LTE) using a revised version of the code MOOG (Sneden 1973) (with the abfind driver) and a grid of (Kurucz 1993) ATLAS9 atmospheres. In a few particular cases, some lines were not used in the analysis due to the poor quality of the spectrum in the region of interest.
In Tables 4 through 6 we summarize the derived abundances for all stars with and without planetary-mass companions.
Errors and uncertainties can affect abundance measurements in various ways. Individual lines can be affected by errors in the EW, caused by unnoticed blends or a poor location of the continuum. Given that we have usually measured more than one line of each element in a given star, the dispersion around the average abundance gives us an idea about the errors due to these effects. Tables 4 through 6 show that the observed dispersions are usually below 0.05 dex, attesting to the quality of the line-list used.
Table 2 shows the abundance sensitivity to changes in
effective temperature,
,
microturbulence and [Fe/H] in three of our targets.
A perturbation of 50 K in effective temperature, 0.15 dex in log g, 0.10 km s-1
in microturbulence and 0.05 dex in [Fe/H]
(typical uncertainties on our parameters - Santos et al. 2004c)
leads to a total typical uncertainty below 0.06 dex in the derived abundances.
Systematic errors in the stellar parameters and/or NLTE effects (the latter not taken
into account in the current analysis) may also affect the derived
abundances (see e.g. Shi et al. 2004).
In Figs. 1 and 2 we present the abundance ratios
[Na/Fe], [Mg/Fe] and [Al/Fe] as functions of the effective temperature and surface
gravity. The results show that a very small dependence seems to exist for the abundances
as a function of
,
in particular for [Na/Fe]. A fit to the
points in Fig. 1
shows, however, that the dependence is always below 0.10 dex/1000 K. Shi et al. (2004)
have shown that the Na lines used in the current
paper are particularly insensitive to NLTE effects.
Table 2: Abundance sensitivity to changes of 50 K in effective temperature, 0.15 dex in log g, 0.10 km s-1 in microturbulence and 0.10 dex in [Fe/H].
In Fig. 3 we present the distributions of [X/H] (with X = Na, Mg and Al) for stars with planets and comparison stars, as well as their cumulative functions. These histograms are similar to those presented for [Fe/H] in Santos et al. (2004c), and indicate that the excess metallicity for stars with planets is, as expected, not unique to iron. In Table 3, we list the average abundances for each element in the two samples, as well as its mean standard deviation and the differences in abundances between stars with planets and comparison stars. These differences vary between 0.20 dex (Mg) and 0.28 dex (Na), and are comparable to those found for other metals (Bodaghee et al. 2003; Ecuvillon et al. 2004a,b; Santos et al. 2004c).
![]() |
Figure 1:
[X/Fe] as a function of
|
| Open with DEXTER | |
![]() |
Figure 2:
Same as Fig. 1 but for |
| Open with DEXTER | |
The [Na/H] and [Al/H] distributions for planet-host stars also reveal an interesting feature
already discussed in the literature (e.g. Bodaghee et al. 2003; Santos et al. 2001)
for other elements: the distributions are not symmetrical, increasing as a function
of increasing [X/H]. Possible interpretations
for this are discussed in e.g. Santos et al. (2004c,2001). The [Mg/H] distribution also
looks slightly bimodal. The reason for the lack of stars with [Mg/H]
0.3 is not
clear, and may be simply due to small number statistics.
We have investigated possible dependences of the Na, Mg and Al abundances with the parameters of the orbiting planets (period, mass, eccentricity). No statistically significant differences were found.
![]() |
Figure 3: Upper panels: distributions of [Na/H], [Mg/H] and [Al/H] for stars with planets (solid histogram) and for the comparison sample (dashed histogram). Lower panels: cumulative functions for both samples. In all cases, Kolmogorov-Smirnov probabilities are of the order of 10-8 that the two samples belong to the same population. |
| Open with DEXTER | |
Table 3:
Average abundances
[X/H]
for stars with planets and for comparison
sample stars. The rms around the mean and the abundance difference between the two
samples for each element are also listed.
In Gonzalez & Laws (2000) and Gonzalez et al. (2001) some possible anomalies concerning Na, Mg and Al were discussed. These authors have found that the [X/Fe] abundance ratios (with X = Na, Mg and Al) for stars with planets seemed to be slightly lower than those found for field dwarfs. Such a difference was, however, not found by other authors (Sadakane et al. 2002).
In Figs. 4 and 5 we present plots of [X/Fe] and [X/H] vs. [Fe/H] for the three elements. Filled circles represent planet-host stars, while open circles denote "single'' field dwarfs. As we can see from the plots, there is no clear and distinctive difference between the two samples. We do not confirm that planet-host stars have lower [X/Fe] ratios for a given metallicity ([Fe/H]). The abundance distributions of stars with planets are high-[Fe/H] extensions to the curves traced by the field dwarfs, and no discontinuity seems to exist.
We should caution, however, that the number of comparison stars plotted is quite small, in particular for the metal-rich domain. An extension of the samples may thus be needed to support this discussion. Such a study is currently in progress.
Aluminum and magnesium were shown to be good tracers of the thin- and thick-disk populations, the latter presenting higher [X/Fe] values for a given [Fe/H] (e.g. Fuhrmann 2004; Bensby et al. 2003). The fact that no clear difference exists between planet-hosts and single field dwarfs suggests that both groups of objects belong (statistically) to the same population (the thin disk).
Although the main goal of this work is to compare the abundances of Na, Mg and Al for planet-hosts and stars without giant planets, the results presented also give us a chance to study the chemical evolution trends of the Galaxy, and in particular in the metal-rich domain. Given the absence of significant differences between stars with planets and comparison stars, we will consider the samples as a whole for the rest of the paper.
![]() |
Figure 4: Abundance ratios [X/Fe] for Na, Mg and Al as functions of metallicity. Filled circles are planet-host stars and the empty triangles are comparison stars. Dotted lines represent solar values. |
| Open with DEXTER | |
![]() |
Figure 5: Abundance ratios [X/H] for Na, Mg and Al as functions of [Fe/H]. Symbols as in Fig. 4. |
| Open with DEXTER | |
From Fig. 4, we can see that the distribution of [Na/Fe] presents a shallow
decrease as a function of [Fe/H] in the domain
.
For solar metallicity the average [Na/Fe] abundances are clearly below solar.
For
,
we notice an uptrend. Interestingly, there
is some hint that the [Na/Fe] values may become constant (and near the solar value)
for [Fe/H] above
0.25 dex. The average [Na/Fe] of the stars with
is solar (0.00 dex, with a dispersion of 0.10 dex).
Unfortunately, the abundance distribution is affected
by the dispersion of the points. It is not clear whether
this dispersion is real or an artifact of the analysis (see Sect. 2.1).
The global trend observed in this paper confirms the results presented
by other authors (e.g. Shi et al. 2004; Edvardsson et al. 1993; Reddy et al. 2003; Feltzing & Gustafsson 1998; Bensby et al. 2003), with the exception of Chen et al. (2000).
The possible constance of the [Na/Fe] values
for [Fe/H] above
0.25 dex is also tentatively seen
in Fig. 13 of Bensby et al. (2003), although it is not present
in the NLTE study of Shi et al. (2004).
Sodium, as well as aluminum, is though to be mostly a product of Ne and C burning in massive stars, the same objects that will later give rise to SN II, the major producers of Mg (Arnett 1996). As we can see from Fig. 4, however, Na and Mg do not behave similarly as a function of stellar metallicity. This difference, also discussed in Bensby et al. (2003), may indicate the existence of different sources for one or both of these elements (SN Ia and AGB stars - Mowlavi 1999; Karakas & Lattanzio 2003; Tsujimoto et al. 1995), a dependence of the SNe yields of both elements with stellar metallicity, or a variation in the spectrum of the progenitor masses inside the metallicity regime studied here. Interestingly, Na abundances increase (globaly) almost exactly as iron (see Fig. 5), which may indicate a similar origin for both elements.
As for Na, the [Mg/Fe] distribution decreases with increasing metallicity
until
,
although the decrease is stronger than for Na
(the [Na/Fe] values always remain close to solar for
).
For higher metallicities, however, the [Mg/Fe] values become constant at an average of
dex (where the error bar represents the rms), and no clear
upturn in observed, although a slight increase of [Mg/Fe] as a function
of [Fe/H] cannot be excluded. The average [Mg/Fe] (for a given [Fe/H]) remains above
solar for the whole metallicity regime studied.
The Mg trend observed in Fig. 4 is very similar to the one presented by other authors (e.g. Chen et al. 2000; Reddy et al. 2003; Feltzing & Gustafsson 1998; Bensby et al. 2003).
Comparison of the Mg abundances thought to be produced by massive stars (Arnett 1996)
with the ones found for other alpha-elements like Si, Ca and Ti, also
produced by SN Ia events (e.g. Thielemann et al. 2002), should provide
evidence about the origin of these elements. If we compare the Mg abundances
derived in the current paper with those for other alpha elements
presented in Bodaghee et al. (2003) (see their Fig. 2), we see a global agreement
with silicon, although for metallicities below solar, the [Si/Fe] abundances
tend to be lower than those found for [Mg/Fe]. For
calcium (presenting a constant downward trend as a function of [Fe/H]), and to
a lesser extent for Ti, however, we do not find the same trends for
.
This may indicate that the yields for these two latter elements may be
metallicity-dependent (e.g. due to a progenitor mass-metallicity dependence),
or that different sources for these different elements give important contributions
(e.g. SN Ia and AGB stars - Karakas & Lattanzio 2003; Thielemann et al. 2002).
Our result regarding this issue is similar to the one found by Bensby et al. (2003).
The Al distribution is very similar to the one observed for Mg,
with the [Al/Fe] values decreasing with increasing [Fe/H] up to solar values,
then becoming relatively constant for metallicities above solar at an average
value of
dex.
Contrary to Mg, however, for solar metallicities
the average [Al/Fe] of our stars is close to 0.0 dex.
Again, we cannot exclude a very shallow uptrend for above-solar metallicities.
Globally, the trend observed for Al is similar to the one already presented by other authors (e.g. Edvardsson et al. 1993; Reddy et al. 2003; Feltzing & Gustafsson 1998), although the slight upturn seen in Fig. 13 of Bensby et al. (2003) for [Fe/H] above solar is not clear in our data. Also, as for Na, our results are considerably different from those obtained by Chen et al. (2000).
The very similar trend observed for Al and Mg suggests a similar origin for both elements.
In Sect. 3.1 we have seen that no large differences in the [X/Fe] ratios were found between planet-hosts and single stars. However, it is still possible that some specific trends found for the metal-rich tail of the Galactic chemical evolution are due to the presence of planets. Indeed, planets are more frequent around metal-rich stars (e.g. Santos et al. 2004c). It should thus be no surprise if their presence influences the yields of certain elements for metal-rich stars in the AGB phase responsible for the chemical enrichment of the Galaxy.
It is known that in addition to the CNO elements, Na, Mg and Al abundances are
correlated in a large number of globular cluster stars (Kraft 1994) showing
the general effects of mixing. It also became clear that Na and Al
enhancements observed in many globular-cluster giants (Kraft 1994) can be
produced by the 22Ne(p,
)23Na reaction of the NeNa-cycle
(Langer et al. 1993, and references therein). For this to occur, one would need
material that had undergone significant O
N or C
N processing,
while in the O-depleted layers we could expect a transformation of the 26Mg and
25Mg isotopes into Al in the MgAl-cycle (Langer et al. 1993).
Langer et al. (1997) have also shown that the use of
updated rates of Ne-Na cycle reactions can explain the Na-O abundance anticorrelation observed in
globular-cluster giants, and the Na-N correlation
observed in field halo giants.
It is thus apparent that giants are not reliable
for studying the primordial abundances of CNO-cycle elements. Clearly, it is
not safe to use CNO, Na, Al and Mg elements in giants as probes of Galactic
chemical evolution since rotation and mixing can alter surface abundances of these
elements.
Furthermore, it is not clear how the yields of these elements will be altered in
metal-rich giants due to the presence of planets.
As a tentative example, if a star during its AGB evolution engulfs
(short period) planets (many such planets have been found by current radial-velocity
surveys
),
it will change its angular momentum and rotation. As a consequence, its internal mixing
history may be changed, possibly bringing different quantities of produced elements
to the stellar outer layers, thus changing the chemical yields. Planet
engulfment, eventually followed by extra-mixing, has already been suggested
to explain the existence of Li-rich giant stars (e.g. Alexander 1967; Siess & Livio 1999).
Overall, these processes could produce important changes in the chemical evolution trends
of the Galaxy, in particular for the metal-rich domain.
In this paper we have derived abundances of Na, Mg and Al in a set of 98 planet-host stars and 41 comparison stars without known planets. The results were derived through a detailed and uniform spectroscopic analysis and gave us the possibility to compare the two samples. The main goal of this comparison is to look for differences connected to the presence of giant planets, but it was also used to explore chemical evolution trends in the Galaxy.
Through a comparison of Na, Mg and Al abundances of planet-hosts and stars without planets, we found an overabundance similar to that already found in iron. However, the abundance ratios [X/Fe] (with X = Na, Mg and Al) found for planet-hosts were not found to be significantly different from those obtained for our comparison stars for a given [Fe/H]. Finally, the trends in the [X/Fe] vs. [Fe/H] plots derived in this paper are similar to the ones presented in the literature, although some new interesting details seem to appear in the metal-rich domain. A discussion of possible implications for the chemical evolution of the Galaxy is presented.
Given the lack of "single'' comparison stars with
dex, the
conclusions presented above must be seen as preliminary.
Indeed, future studies of elemental abundances in planet-host stars
should include a comparison with a more significant number of "single''
dwarfs having metallicities above solar. Such work is currently in progress.
Acknowledgements
Support from Fundação para a Ciência e Tecnologia (Portugal) to N.C.S. in the form of a scholarship is gratefully acknowledged. We wish to thank the Swiss National Science Foundation (Swiss NSF) for the continuous support for this project.
Table 4:
Atmospheric parameters (from Santos et al. 2004c) and abundances of
Na, Mg and Al for stars with planets from HD 142 to HD 95128. The number of spectral lines used is
given by n while
denotes the rms around the average.
Table 5: Same as Table 4 for stars with planets from HD 106252 to HD 330075.
Table 6: Same as Table 4 for the comparison sample of stars without known planets.