A&A 437, 1127-1133 (2005)
DOI: 10.1051/0004-6361:20052895
N. C. Santos1,2 - G. Israelian3 - M. Mayor2 - J. P. Bento1 - P. C. Almeida1 - S. G. Sousa1 - A. Ecuvillon3
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 17 February 2005 / Accepted 23 March 2005
Abstract
We present stellar parameters and metallicities for
29 planet-host stars, as well as for a large volume-limited sample of
53 stars not known to be orbited by any planetary-mass companion.
These stars add to the results presented in our previous series
of papers, providing two large and uniform samples of
119 planet-hosts and 94 "single'' stars with accurate stellar parameters
and [Fe/H] estimates. The analysis of the results further confirms that stars
with planets are metal-rich when compared with average field dwarfs.
Important biases that may compromise future studies are also discussed.
Finally, we compare the metallicity distributions for single planet-hosts
and planet-hosts in multiple stellar systems. The results show that
a small difference cannot be excluded, in the sense that the latter sample
is slighly overmetallic. However, more data are needed to confirm this correlation.
Key words: stars: abundances - stars: fundamental parameters - stars: chemically peculiar - stars: planetary systems - stars: planetary systems: formation
The discovery that stars with giant planet companions are particularly metal-rich when compared to single field dwarfs is providing crucial clues to our understanding of the processes of planetary formation and evolution. This fact, first discussed by Gonzalez (1997) and Gonzalez (1998), was later proved by the uniform analysis of large samples of planet-host and field "single'' stars (Santos et al. 2001).
Further refinements to these studies have shown that the probability of finding (and presumably forming) a planet, at least of the kind now discovered by radial-velocity surveys, is a strong function of the stellar metal content (Reid 2002; Santos et al. 2003; Fischer & Valenti 2005; Santos et al. 2004b,2001). Although the exact dependence is still not completely settled (see Santos et al. 2004b, and references therein), it was shown that about 25% of solar type stars with twice the solar metallicity are orbited by a planet, while less than 5% of stars with the metal content of the Sun have an orbiting planetary-mass companion.
Interestingly, these observational results are now being corroborated by theoretical models of planet formation, and give support to the core-accretion model to explain the origin of the now discovered giant planets (Alibert et al. 2004; Pollack et al. 1996). Current models can even reasonably explain the observed [Fe/H] distribution of planet-host stars (Ida & Lin 2004a,b). However, and as discussed in Santos et al. (2004b), current data still cannot exclude that giant planets may also be formed by the disk instability model (Mayer et al. 2002; Boss 1997,2002), in particular planets orbiting low metallicity stars.
In the last few years we have published a series of papers regarding the study of the stellar metallicity-giant planet connection (Santos et al. 2003,2004b,2001,2000) (hereafter Papers I, II, III and IV, respectively). In the last three of these papers we have presented a spectroscopic analysis of a sample of 41 "single'' field stars, used as a reference to study the excess metallicity observed in planet hosts. In Paper IV we presented a uniform spectroscopic analysis of 98 extra-solar planet-host stars, including most of the known objects. The continuation of this work is extremely important, since the search for new possible correlations between the presence (and properties) of exoplanets and the chemical abundances of their host stars needs increased samples.
In this paper we add 53 new stars to our "comparison'' sample, as well as precise and uniform metallicity estimates for a further 29 known planet-hosts, most of which have no previous spectroscopic metallicity estimate. In Sects. 2 and 3 we present the observed samples, and in Sect. 4 we review the metallicity distribution of planet-host stars, as well as of the whole volume-limited comparison sample of 94 solar-type F-G-K dwarfs. The analysis confirms, as expected, that stars with planets are clearly metal-rich when compared to average field dwarfs. In Sect. 4.1 we review a few possible biases that should be considered when studying the metallicity of planet-host stars. In Sect. 5 we compare the metallicity distributions for "single'' planet-host stars with that found for planet-hosts in multiple stellar systems. We conclude in Sect. 6.
In Papers II through IV we compared the metallicity distribution for stars
with planets with a volume-limited sample of stars
having right-ascensions (RA) between 20
and 9
,
and a
distance below 20 pc as derived from Hipparcos parallaxes (ESA 1997).
All these stars belong to the CORALIE southern planet search sample (Udry et al. 2000),
have spectral types between F8 and M1, and were not found to harbor any planetary companion.
In order to extend the number of comparison stars in our samples, we have now
observed a complementary sample of CORALIE stars with RA between 9
and 20
(and d< 20 pc). These stars were observed using the FEROS spectrograph at the 2.2-m ESO/MPI telescope (ESO, La Silla) in March 2004 (observing run ID 072.C-0033).
Data reduction was done using the FEROS pipeline, and the wavelength calibration
was done using the spectrum of a ThAr lamp. The final spectra have a resolution
R =
= 48 000, and the S/N between
200 and 400.
A few stars from the former comparison sample for which we could not obtain a spectrum before were also observed, both using FEROS and the CORALIE spectrograph at the 1.2-m Euler Swiss telescope, also at La Silla observatory (R = 50 000).
The comparison sample presented in this paper is listed in Table 2. Together with Table 5 of Paper IV, this constitutes a large volume-limited sample of 94 "single'' stars (with no known planetary-mass companions). Amid the stars in this new volume-limited comparison sample there is HD 147513 and the very metal rich planet-host star HD 160691 (Santos et al. 2004a; Butler et al. 2001). For the list of other planet-hosts that belong to the whole volume-limited comparison sample we refer to Santos et al. (2004b).
In addition to the new comparison sample, spectra of some planet hosts was also obtained,
both during the March run and in two complementary runs, one in October 2004 (using FEROS -
observing run ID 074.C-0135) and another in December 2004, using the UVES spectrograph
at the VLT 8.2-m Kueyen telescope (run ID 074.C-0134). These latter spectra were reduced using
the UVES pipeline and have a resolution R
75 000.
Table 1: Planet-host stars studied in this paper. See text for more details.
Table 2: Comparison sample stars presented in this paper. See text for more details.
For HD 104985 and HD 219449 we further obtained spectra
with a resolution of 57 000 using the SARG spectrograph, at the TNG telescope
(La Palma Observatory, Spain). For a few other recently-announced planet-hosts,
stellar parameters were derived using the HARPS spectrograph (R
100 000 - for
details see e.g. Pepe et al. 2004). These parameters were also presented in
the papers announcing the discovery of the exoplanets (see Table 1).
The list of planet-host stars presented in this paper
is shown in Table 1. For planet-discovery references see the
Geneva Planet-Search web pages.
Adding the 97 planet-host stars whose [Fe/H] values were listed in Paper IV to the
results listed in Table 1, we now have
a uniform sample of metallicities and stellar parameter for 119 stars
with orbiting planets. From the whole known planet-host star sample, the
only 4 objects for which we could not obtain metallicity estimates are GJ 876
and GJ 436 (both M-dwarfs - see Bonfils et al. 2005), HD 45350 and HD 13189 (a giant star). To these we add the transiting planets detected by the OGLE survey
(Pont et al. 2004; Konacki et al. 2003; Bouchy et al. 2004), three currently with poor estimates given
the low quality of the now available spectra, as well as
TrEs-1 (Alonso et al. 2004).
The spectroscopic analysis was done in LTE using the
2002 version of the code MOOG (Sneden 1973) and a grid of Kurucz Atlas plane-parallel
model atmospheres (Kurucz 1993).
Stellar parameters and metallicities were derived as in Santos et al. (2004b),
based on the Equivalent Widths (EW) of 39 Fe I and 12 Fe II week lines, and by imposing
excitation and ionization equilibrium. The line EW were measured using IRAF
"splot'' and "bplot''
routines within the echelle package. A Fortran 77 code that uses a Downhill Simplex Method
(Press et al. 1992) was then used to find the best solution in
the (stellar) parameter space once the line EWs are measured.
The results of this analysis are presented in Tables 1 and 2. As shown in Paper IV, this methodology gives excellent results for the derived stellar parameters, compatible with other spectroscopic results. A comparison between the EW and stellar parameters derived from these spectra and data obtained using other instruments also shows that no significant differences exist (e.g. Santos et al. 2004b) - see also Table 1.
The errors in
,
,
and [Fe/H] were derived as
in Paper IV, and are typically of the order of 50 K, 0.2 dex, 0.1 km s-1and 0.05 dex, respectively. Higher uncertainties are found for
lower temperature stars.
In the tables we also present the stellar masses. These were
derived from an interpolation of the theoretical isochrones of Schaller et al. (1992),
Schaerer et al. (1993a) and Schaerer et al. (1993b), using MV computed using Hipparcos
parallaxes (ESA 1997), a Bolometric Correction (BC) from Flower (1996)
and the
obtained from the spectroscopy. Although this method
to estimate individual stellar masses may suffer from important
uncertainties (Fernandes & Santos 2004), from a statistical point of view it
should provide good results.
We adopt a typical relative error of 0.05
for the masses. In some cases,
no mass estimates are presented, since these involved large extrapolations of the
isochrones.
From the derived stellar mass, the spectroscopic temperature and BC, we could also obtain "trigonometric'' surface gravities for our stars using the basic formula presented in Paper IV. The results (see Tables 1 and 2) show that the two gravity estimates are perfectly compatible within the errors, as expected from our previous studies (see Paper IV). Considering also all planet-hosts and comparison sample stars presented in Paper IV, the average difference between the two estimates is nearly zero (0.004 dex), with a stardard deviation of 0.14 dex.
The stellar parameters for HD 41004 A may suffer from the fact that the spectrum of this star is blended with the spectrum of the fainter M-dwarf companion HD 41004 B (see e.g. Santos et al. 2002). The separation between the two components of the binary is only 0.5'', smaller than the diameter of the FEROS fiber. However, given the magnitude difference between the two components, the flux of HD 41004 B should represent a mere 3% of the total flux of the primary; the impact on the determination of the line EW for HD 41004 A is thus probably limited.
In a recent paper, Ramírez & Meléndez (2004) derived temperatures and radii for
extra-solar planet-host stars using the IRFM technique. In their work
they compare the temperatures derived in Paper IV, as well
as by Ribas et al. (2003) (who used infra-red photometry), with the ones obtained by them,
and conclude that our values are systematically shifted by 100 K.
One argument used by Ramírez & Meléndez (2004) in favor of their temperature scale is the fact that for 5 stars that have accurate stellar radii measured (and thus accurate temperatures Ramírez & Meléndez 2004), the results from the IRFM differ by only 14 K with respect to the measured "direct'' temperatures (see their Table 1).
Four of the stars in their Table 1 were also analyzed by us.
These are HD 128620, HD 128621 (see Table 2),
HD 10700 and HD 209458 (see Paper IV). Interestingly, a comparison of
the obtained
shows that the difference between our estimates and the "direct''
temperatures listed by Ramírez & Meléndez (2004) is only of 7 K. In other words,
we do not find any strong reason to support that the IRFM scale of temperatures
is more accurate than the one derived from the method used in our studies
(Paper IV and this paper).
The major goal of the current study is to compare the [Fe/H] values for planet-host stars and "single'' field dwarfs in an uniform way. The crucial point is thus to use uniform samples, analyzed using the same line-lists, model atmospheres and techniques to derive stellar parameters and metallicities. This goal is already achieved, and thus our conclusions are not dependent on the temperature scale used. The definition of a correct scale, which is indeed very important when an absolute study is needed, is thus not a limitation in our case.
In Fig. 1 we present the metallicity distribution for the 119 planet-host stars with metallicity estimates (Paper IV and Table 1), as well as for the whole volume-limited comparison sample of 94 "single'' dwarfs (Paper IV and Table 2).
As it is well seen from the plot, and confirming previous studies, planet hosts have, on average, a [Fe/H] that is 0.24 dex higher than that found for stars with no detected planetary mass companions. The Kolmogorov-Smirnov probability that the two samples belong to the same population is of the order of 10-12.
Although it is not our goal here to study in detail the metallicity distribution of planet-hosts (see next section), or to study the frequency of planets as a function of stellar metallicity (see Paper IV), the addition of new stars to our comparison sample reinforces the results presented in previous works (e.g. Papers II and IV). For instance, the average [Fe/H] of the stars in Table 2 is -0.11, similar to the one found for the 41 comparison stars presented in Paper IV.
The metallicity distribution for our 94 comparison sample stars though symetrical, seems to have a small depression for stars with [Fe/H] values between -0.1 and -0.2 dex. A similar scarcity of planet-host stars seems to exist at the same [Fe/H] interval. The cause for this feature, also tentatively observed in our previous samples (for iron and for other elements - Bodaghee et al. 2003; Santos et al. 2001) as well as in the volume-limited sample of Fischer & Valenti (2005), is not clear, but may be simply related to small number statistics. It is interesting to note that this metallicity regime is at the frontier of the Galactic thin and thick disk populations (e.g. Fuhrmann 2004). Further studies of the UVW velocities of the stars in our sample as well as of other chemical elements are necessary to better tackle this issue.
![]() |
Figure 1: Metallicity distribution for 119 planet-host stars (dashed line) and for a volume-limited "comparison'' sample of 94 stars with no known planetary-mass companions (continuous line). The metallicities were taken from Table 1, Table 2 and from Tables 2-5 of Santos et al. (2004b). The difference between the average metallicity of two samples is 0.24 dex. Inset: cumulative functions for both distributions. A Kolmogorov-Smirnov test shows that the probability that both distributions belong to the same population is of the order of 10-12. |
Open with DEXTER |
In Paper III we have shown that the whole metallicity "excess'' observed for stars with giant planets cannot be explained by any observational or sampling bias. However, the striking fact that stars with planets are particularly metal-rich has led some planet-search programs to base their samples on metal-rich stars (e.g. Fischer et al. 2004; Tinney et al. 2002). Such surveys, now giving their first results, are adding to the lists of known planet hosts some objects that should be considered with care if a consistent statistical study about the stellar metallicity-giant planet connection is pursued.
The impact of these surveys may have consequences on any future statistical studies of the properties of exoplanets orbiting stars with different chemical composition. For example, the first result of these metallicity-biased surveys will be the discovery of short period giant planets (hot-Jupiters). Any study of the relation between the metallicity of the stars and the orbital period of the orbiting planet (Sozzetti 2004; Gonzalez 1998; Queloz et al. 2000; Santos et al. 2003) will thus be seriously compromised, as such systematic errors are very difficult to take into account.
Such was not taken in consideration when plotting stars in Fig. 1, as the goal of the present paper is simply to present stellar metallicities for an extended comparison sample, and for a few planet-host stars whose stellar parameters and [Fe/H] were not presented in Paper IV. A detailed and unbiased comparison of the [Fe/H] distributions for planet-hosts and "single'' stars implies the use of samples for which we are sure not to have any important bias. This was the case of our previous studies (e.g. Paper IV), where we have compared planet-hosts found amid the stars in the CORALIE planet-search sample (Udry et al. 2000) with the [Fe/H] distribution of the whole sample itself; with CORALIE, a target was never chosen based on its high metal content.
In Fig. 2 we present a comparison between the [Fe/H] distributions of 22 stars with planets that are currently known to belong to multiple stellar
systems - mostly binary systems, hereafter PIB -
with the same distribution for "single'' planet-host stars. See Eggenberger et al. (2004) and Eggenberger et al. (in prep.) for references.
![]() |
Figure 2:
Metallicity distribution for 22 planet-host stars in binary stellar systems (PIB - filled line) when compared with the same distribution for single planet-hosts (dashed line).
The metallicities were taken from Table 1 and from Santos et al. (2004b).
The difference between the average metallicity of two samples is ![]() |
Open with DEXTER |
The comparison shows that no major difference seems to exist between the two samples. Even if PIB stars seem to be slighly over-metallic when compared to single planet-hosts, the Kolmogorov-Smirnov probability that the two samples belong to the same population is 0.08. If we exclude HD 114762 from the PIB sample (this star has a companion with a minimum mass above 10 times the mass of Jupiter, likely a brown dwarf Latham et al. 1989), this probability decreases to 0.04, only marginally significant.
It has been proposed that planet migration processes may act differently in single and multiple stellar systems, and that migration rates may be increased if the planet-host has a stellar companion (for a discussion see Eggenberger et al. 2004). Increasing the migration efficiency could then increase the scattering within the disc of planetesimals (e.g. Murray et al. 1998). In such a case, we could expect that PIB stars, or at least those with "close'' stellar companions, could have been polluted by a larger amount of planetesimal (metal-rich) material.
The fact that we do not find a clear difference between PIB stars and single planet-hosts suggests that strong planetesimal accretion is not usual both for single planet-hosts (Cody & Sasselov 2004; Pinsonneault et al. 2001; Santos et al. 2003) nor for planet-host stars in multiple stellar systems. However, if the slight correlation discussed above is real it may imply, for example, that stellar pollution played some role (although probably small) in planet-hosts in multiple stellar systems. A follow-up of these results as more planets are discovered is necessary to resolve this issue.
We also checked if within the PIB star sample there was any correlation between stellar metallicity and the separation of the stellar components. No major correlation was found.
In this paper we have presented accurate stellar metallicities for a volume-limited sample of 53 comparison "single'' stars. This sample was observed to complement the already existing list of 41 comparison dwarfs presented in Papers II and IV.
We further present metallicities and stellar parameters for 29 planet-host stars, most of which were not studied in our previous series of papers (Papers I-IV). Together with the 97 planet host stars studied in Paper IV, we now have a large sample of 119 stars with planets (almost all the known targets) with precise metallicity and stellar parameter estimates.
These two large samples (of planet-hosts and "single'' stars) were analyzed using the same line-lists, model atmospheres and techniques to derive stellar parameters and [Fe/H]. Therefore, they constitute very uniform samples, unique and appropriate for future studies of the stellar metallicity-giant planet connection.
We have used these samples to investigate whether stars with planets belonging to multiple stellar systems have a different [Fe/H] distribution when compared with single planet-hosts. The results show that a slight tendency, although not statistically significant, may exist. The enlargement of the samples is needed to resolve this issue.
Acknowledgements
We would like to thank the anonymous referee for a prompt report. 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.