A&A 438, 251-256 (2005)
DOI: 10.1051/0004-6361:20052750

Abundances of Na, Mg and Al in stars with giant planets[*],[*]

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

   
1 Introduction

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].

   
2 Data and spectral analysis

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 $\log{gf}$ 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 $T_{{\rm eff}}=5777$ K, $\log{g}=4.44$ dex, $\xi_t=1.00$ km s-1 and $\log{\epsilon}\rm (Fe)=7.47$ 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.

   
2.1 Errors

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, $\log{g}$, 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 $T_{\rm eff}$, 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].

   
3 Na, Mg and Al in planet-host stars

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).


  \begin{figure}
\par\psfig{width=8.5cm,file=2750fig1.eps}\end{figure} Figure 1: [X/Fe] as a function of $T_{\rm eff}$. The dotted line represents the solar [X/Fe] value.
Open with DEXTER


  \begin{figure}
\par\psfig{width=8.5cm,file=2750fig2.eps}\end{figure} Figure 2: Same as Fig. 1 but for $\log{g}$.
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]$~\sim~$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.


  \begin{figure}
\par\psfig{width=8.8cm,file=2750fig3.eps}\end{figure} 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 $\langle$[X/H]$\rangle$ 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.

   
3.1 [X/Fe] vs. [Fe/H]

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).

   
4 Galactic chemical evolution trends

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.


  \begin{figure}
\par\psfig{width=8.5cm,file=2750fig4.eps}\end{figure} 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


  \begin{figure}
\par\psfig{width=8.5cm,file=2750fig5.eps}\end{figure} Figure 5: Abundance ratios [X/H] for Na, Mg and Al as functions of [Fe/H]. Symbols as in Fig. 4.
Open with DEXTER

   
4.1 Sodium

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 $-0.7<\mbox{[Fe/H]}<0.0$. For solar metallicity the average [Na/Fe] abundances are clearly below solar. For $\rm [Fe/H]>0$, 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 $\sim$0.25 dex. The average [Na/Fe] of the stars with $\rm [Fe/H]>0.0$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 $\sim$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.

   
4.2 Magnesium

As for Na, the [Mg/Fe] distribution decreases with increasing metallicity until $\rm [Fe/H]\sim0.0$, although the decrease is stronger than for Na (the [Na/Fe] values always remain close to solar for $\rm [Fe/H]<0.0$). For higher metallicities, however, the [Mg/Fe] values become constant at an average of $0.06\pm0.06$ 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 $\rm [Fe/H]>0.0$. 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).

   
4.3 Aluminum

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 $0.07\pm0.07$ 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.

4.4 Anomalies and planets

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, $\gamma$)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 $\rightarrow$N or C $\rightarrow$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.

   
5 Concluding remarks

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 $\rm [Fe/H]>0.2$ 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.

References

 

  
Online Material

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 $\sigma $ 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.



Copyright ESO 2005