Issue |
A&A
Volume 566, June 2014
|
|
---|---|---|
Article Number | A83 | |
Number of page(s) | 16 | |
Section | Planets and planetary systems | |
DOI | https://doi.org/10.1051/0004-6361/201423803 | |
Published online | 18 June 2014 |
Chemical abundances of stars with brown-dwarf companions⋆
1 Instituto de Astrofísica de Canarias, 38200 La Laguna, Tenerife, Spain
e-mail: dmata@iac.es; jonay@iac.es
2 Departamento de Astrofísica, Universidad de La Laguna, 38206 La Laguna, Tenerife, Spain
3 Centro de Astrofísica, Universidade do Porto, Rua das Estrelas, 4150-762 Porto, Portugal
4 Departamento de Física e Astronomia, Faculdade de Cîencias, Universidade do Porto, Portugal
5 Observatoire Astronomique de l’Université de Genève, 51 Ch. des Maillettes, 1290 Sauverny, Versoix, Switzerland
6 European Space Agency, European Space Astronomy Centre, PO Box 78, Villanueva de la Cañada, 28691 Madrid, Spain
Received: 13 March 2014
Accepted: 2 April 2014
Context. It is well known that stars with giant planets are, on average, more metal-rich than stars without giant planets, whereas stars with detected low-mass planets do not require to be metal-rich.
Aims. With the aim of studying the weak boundary that separates giant planets and brown dwarfs (BDs) and their formation mechanism, we analyze the spectra of a sample of stars with already confirmed BD companions both by radial velocity and astrometry.
Methods. We employ standard and automatic tools to perform an equivalent width (EW) based analysis and to derive chemical abundances from the CORALIE spectra of stars with BD companions.
Results. We compare these abundances with those of stars without detected planets and with low-mass and giant-mass planets. We find that stars with BDs do not have metallicities and chemical abundances similar to those of giant-planet hosts, but they resemble the composition of stars with low-mass planets. The distribution of mean abundances of α-elements and iron peak elements of stars with BDs exhibit a peak at about solar abundance, whereas for stars with low-mass and high-mass planets the [Xα/H] and [XFe/H] peak abundances remain at ~−0.1 dex and ~+0.15 dex, respectively. We display these element abundances for stars with low-mass and high-mass planets, and BDs versus the minimum mass, mC sin i, of the most massive substellar companion in each system, and we find a maximum in α-element as well as Fe-peak abundances at mC sin i ~ 1.35 ± 0.20 Jupiter masses.
Conclusions. We discuss the implications of these results in the context of the formation scenario of BDs compared with that of giant planets.
Key words: brown dwarfs / stars: abundances / planets and satellites: formation / planetary systems / stars: atmospheres
Appendices are available in electronic form at http://www.aanda.org
© ESO, 2014
1. Introduction
The most widespread convention places the mass range of brown dwarfs (BDs) at 13 − 80 MJ (being MJ the mass of Jupiter), having enough mass to burn deuterium but not for hydrogen fusion (Burrows et al. 1997), i.e., in between the heaviest giant planets and the lightest stars. Brown dwarfs above ~65 MJ are thought to fuse lithium and therefore the detection of the Li iλ 6708 Å can be used to identify BDs; this is the so-called lithium test (Rebolo et al. 1992).
Brown dwarfs were predicted by Kumar (1962) and Hayashi & Nakano (1963), but they were not empirically confirmed until 1995, when the first field BD was detected (Teide 1, Rebolo et al. 1995). This occurred the same year as the discovery of the first extra-solar planet (Mayor & Queloz 1995). The first BD companion to a M-dwarf star was also discovered that year (GJ 229B, Nakajima et al. 1995). During the following two decades high-precision radial velocity (RV) surveys have shown that close BDs around solar-type stars are rare (Grether & Lineweaver 2006, and references therein). Thus, at orbital separations of less than 10 AU, the frequency of BD companions remains below 1% (Marcy & Butler 2000), whereas it is ~7% for giant planets (Udry & Santos 2007; Mayor et al. 2011) and ~13% for stellar binaries (Halbwachs et al. 2003). The so-called Brown dwarf desert may be interpreted as the gap between the largest-mass objects that can be formed in protoplanetary disks, and the smallest-mass clumps that can collapse and/or fragment in the vicinity of a protostar (Ma & Ge 2014). The mass function, d, of close planetary and stellar companions drops away (α ~ − 1) towards the BD mass range (Grether & Lineweaver 2006). On the other hand, the mass function of isolated substellar objects is roughly flat or even with linear increase (α ~ 0) down to ~20 MJ (Chabrier 2002; Kirkpatrick et al. 2012). This may point to a different formation scenario for close BD companions and BDs in the field and clusters.
Sahlmann et al. (2011) presented the discovery of nine BD companions from a sample of 33 solar-type stars that exhibit RV variations caused by a companion in the mass range mC sin i ~ 13 − 80 MJ. They used Hipparcos astrometric data (Perryman et al. 1997) to confidently discard some of the BD candidates. Including literature data, these authors quoted 23 remaining potential BD candidates. From the CORALIE planet-search sample, they obtain an upper limit of 0.6% for the frequency of BD companions around Sun-like stars. Recently, Ma & Ge (2014) have collected all the BD candidates available in the literature including those in Sahlmann et al. (2011), some from the SDSS-III MARVELS survey (Ge et al. 2008) and some other RV surveys (e.g., Marcy & Butler 2000).
The metallicity of stars with BD companions have been briefly discussed in Sahlmann et al. (2011). They note that the sample is still too small to claim any possible metallicity distribution of stars hosting BDs. Ma & Ge (2014) extended the sample to roughly 65 stars with BD candidates, including dwarfs and giants, and stated that the mean metallicity of their sample is ⟨[Fe/H]⟩= − 0.04 (σ = 0.28), in other words remarkably lower than that of stars with giant planets (⟨[Fe/H]⟩= + 0.08, Sousa et al. 2008, 2011). On the other hand, stars with only detected small planets do not seem to require high metal content to form planets within planetary disks (Sousa et al. 2008, 2011; Adibekyan et al. 2012c). Sousa et al. (2011) study a sample of 107 stars with planets (97 giant and 10 small planets) and found an average metallicity of stars with small planets at about ⟨[Fe/H]⟩= − 0.11, very similar to that of stars without detected planets (Sousa et al. 2008). Hereafter “small” refers to low-mass planets, including super-Earths and Neptune-like planets, with mC sin i< 30 M⊕, whereas “giant” planet refers to high-mass planets, including Saturn-like and Jupiter-like planets, with 30 M⊕<mC sin i< 13 MJ (see Sect. 3).
Currently, there are two well-established theories for giant planet formation: the core-accretion scenario (Pollack et al. 1996) and disk gravitational instability (Boss 1997). The core-accretion model is more sensitive to the fraction of solids in a disk than is the disk-instability model. The formation of BDs has also been extensively studied. Two main mechanism have been proposed: molecular cloud fragmentation (Padoan & Nordlund 2004), and disk fragmentation (Stamatellos & Whitworth 2009). The latter mechanism, which requires a small fraction of Sun-like stars that host a massive extended disk, is able to explain most of the known BDs which may either remain bound to the primary star, or may be ejected into the field (Stamatellos & Whitworth 2009).
In this paper, we present a uniform spectroscopic analysis for a sample of stars with BD companions from Sahlmann et al. (2011) and we compare the results with those of a sample of stars with known giant and small planets from previous works (Adibekyan et al. 2012c). The aim of this work is to provide some information that could be useful to distinguish among the different and possible formation mechanisms of BD companions.
Stellar parameters of the CORALIE sample.
2. Observations
We analyze data for two different samples obtained with two different telescopes and instruments: stars with BDs with spectroscopic data at resolving power R ~ 500 000 taken at the 1.2 m Euler Swiss Telescope equipped with the CORALIE spectrograph (Udry et al. 2000), and stars with planetary companions observed with the HARPS spectrograph (Mayor et al. 2003) with R ~ 115 000 installed at the 3.6 m ESO telescope, both of them at La Silla Observatory (ESO) in Chile.
The individual spectra of each star were reduced in a standard manner, and later normalized within the package IRAF1, using low-order polynomial fits to the observed continuum.
3. Sample description and stellar parameters
3.1. Stars with BD companion candidates
Our stellar sample has been extracted mostly from F-, G-, and K-type main-sequence stars of the CORALIE RV survey (Udry et al. 2000). This sample consists of 15 stars with BD companion candidates reported in Sahlmann et al. (2011), for which the minimum mass, mC sin i, of the most massive companion is in the brown dwarf mass range (13 − 80 MJ). One of these 15 stars, HIP 103019, has been extracted from the HARPS RV survey (Mayor et al. 2003). In Table 1 we provide the minimum mass of these 15 BD candidates. Sahlmann et al. (2011) were also able to derive the orbital inclination, i, by using astrometric measurements from Hipparcos (Perryman et al. 1997; van Leeuwen 2007). This allowed them to confidently exclude as BD candidates eight stars from the initial sample of 15 stars because the current mass determinations, mC, place them in the M-dwarf stellar regime. Stellar parameters of the sample of 14 stars were collected from Sahlmann et al. (2011) and one star from Santos et al. (2005). Four additional stars from the CORALIE sample without detected BD companions were analyzed as a comparison/control sample (see Table 1). In Fig. 1 we show the histograms of Teff, log g, and [Fe/H] of our stellar sample. We only display those stars with available mC sin i values in Table 1.
We note that Ma & Ge (2014) also collected a sample 65 stars with BD candidates: 43 stars, 27 dwarfs, and 15 giants, with mC sin i values in the range 13−90 MJ. Two stars included in Ma & Ge (2014) as BD candidates, HD 30501 and HD 43848, were discarded by Sahlmann et al. (2011), probably because they have mC values above but close to the 80 MJ boundary. Figure 3 depicts the minimum mass of the BD companion mc sin i and orbital period against metallicity [Fe/H], including stars of our sample as well as those stars in the sample of Ma & Ge (2014) for comparison. We separate giant (log g< 4) and dwarfs (log g> 4) stars. The stars with BD candidates are spread over a wide range of orbital periods, a minimum masses of the BDs, and stellar metallicities.
3.2. Stars with planetary-mass companions
The HARPS subsample (HARPS-1 sample in Adibekyan et al. 2012c) used in this work contains 451 stars (Sousa et al. 2008; Neves et al. 2009), both with and without planetary companions. We collect the minimum mass of the most massive planet in each planetary system from the encyclopaedia of extra-solar planets2. The planetary-mass sample is separated into two groups: (i) small planets (SP; super-Earth-like and Neptune-like planets) with masses of mC sin i< ~ 0.094 MJ (~30 M⊕); and (ii) giant planets (GP; Saturn-like and Jupiter-like planets) with masses in the range 0.094 <mC sin i [MJ] < 13. Two of the stars with giant planets, HD 162020 and HD 202206, within the HARPS sample have companion masses above 13 MJ and will be considered as BDs hereafter.
Therefore, our final sample of confirmed BDs contains nine dwarf stars with companions in the mass range mC sin i ~ 13 − 80 MJ. In the following, we refer to “BD-host stars” as stars with confirmed BDs, i.e., with mC ~ 13 − 80 MJ, and “stars with discarded BDs” as those with mC sin i ~ 13 − 80 MJ but mC> 80 MJ, following Sahlmann et al. (2011). The sample of planet-host stars contains 25 stars with small planets and 78 stars with giant planets. In Tables A.4 and A.5 we provide the minimum mass of the most massive planet in each planetary system of the stars in the HARPS sample.
![]() |
Fig. 1 Histograms of the stellar parameters Teff, log g, and [Fe/H] of our CORALIE sample. |
4. Automatic codes for EW measurements: ARES versus TAME
![]() |
Fig. 2 ARES versus TAME graphics for the stars HD 89707 (S/N ~ 110) and HD 206505 (S/N ~ 70). Left panels: EW measured with TAME versus ARES. The 1:1 correspondence is shown as a dashed line. EW differences, EWTAME − EWARES, versus EWARES (middle panels) and spectral line wavelength (right panels). Dash-dotted lines define the mean value of the data points, and dotted lines define the mean plus the standard deviation. |
![]() |
Fig. 3 Minimum mass of the most massive companion, mC sin i, and orbital period, PC, versus the metallicity of stars with BDs. The CORALIE sample is depicted as filled symbols, and the sample in Ma & Ge (2014) is displayed as empty symbols, separated between giant and dwarf stars. |
![]() |
Fig. 4 High-resolution normalized spectra of late G-, K-type stars: CORALIE spectra of the coolest stars in the sample (black), HARPS spectra (red) of cool stars from the HARPS database (e.g., Sousa et al. 2008) and the CORALIE spectrum (blue) of HD 53680, and the HARPS spectrum (blue) of HIP 103019. The spectra are depicted following a sequence with decreasing Tefffrom top to bottom. |
We measure the equivalent widths (EWs) of spectral lines using the linelists in Sousa et al. (2008), Neves et al. (2009), and Adibekyan et al. (2012c), using automatic tools. We explore two different automatic codes for EW spectra analysis: the automatic code ARES3 based on C++ (Sousa et al. 2007) and a new code named TAME4 based on IDL (Kang & Lee 2012). In order to compare these two automatic codes, we measure the EWs of the CORALIE sample with the same input parameters to these two codes.
In Fig. 2, we compare the EWs measured using TAME, EWTAME, against those estimated using ARES, EWARES. The mean value of the EW differences EWTAME − EWARES is found at ~−1.2 mÅ and ~−1.5 mÅ for the stars HD 89707 star (S/N ~ 110) and HD 206505 (S/N ~ 70). The TAME code very slightly underestimates the EW compared to the ARES measurements, and the scatter of these comparisons is lower than ~1.5 mÅ. These EW differences do not exhibit any remarkable dependence on wavelength. We also tested whether the signal-to-noise ratio (S/N) from our stellar spectra is the source of this observed tendency and no trend has been found. The mean value of the EW differences fluctuates in the −2 mÅ to −1 mÅ range. The standard deviation of the EW differences improves slightly as the S/N increases, but it oscillates between 0.5 mÅ and 1.5 mÅ.
This analysis lead us to conclude that both programs show good agreement and their differences are not significant and do not have any relevant impact on the chemical abundance analysis, within the typical error bars of the EW-based chemical abundance analysis.
5. Chemical abundances
We compute the EWs using ARES for consistency with chemical abundance analysis in Adibekyan et al. (2012c). We also use version 2010 of the MOOG 5 code (Sneden 1973) together with the Kurucz ATLAS9 stellar model atmospheres (Kurucz 1993) for chemical abundance determination.
We first check the Fe i and Fe ii abundances, using the line list from Sousa et al. (2008). The high dispersion, σ, of the Fe abundances (see Table 2) of the stars HD 53680 and HIP 103019 suggests that these stars may have lower effective temperatures than those given in Table 1. In Fig. 4 we depict the normalized spectra of the coolest stars in the sample together with some spectra of late-G, K-type dwarfs from the HARPS database (e.g., Sousa et al. 2008). The Hα profiles of these two stars do not follow the sequence of temperatures, but they appear to be the coolest objects in Fig. 4. In addition, these stars were discarded as BD candidates by Sahlmann et al. (2011). We note the unusually high surface gravity estimated for the star HD 53680 which is likely spurious, consistent with the large scatter in the Fe abundances and the difference between Fe i and Fe ii abundances. This star is catalogued in the SIMBAD database as a K6V star in a visual binary system. The narrow Hα line and the large dispersion in Fe i and Fe ii abundances may indicate an even later spectral type. For these reasons, they may deserve further analysis and from this point on, we will not consider them.
Fe i and Fe ii abundances and standard deviations.
We use the linelist in Neves et al. (2009) on 17 stars of the CORALIE remaining sample (seven stars with confirmed BDs, six with discarded BD candidates, and four without detected BD companions) to derive element abundances of Na i, Mg i, Al i, Si i, Ca i, Sc i, Sc ii, Ti i, Ti ii, V i, Cr i, Cr ii, Mn i, Co i, and Ni i. Some lines of the original list are not considered: lines non-detected by ARES, lines discarded in the Adibekyan et al. (2012c), and lines whose abundances were out of the 3σ range. The abundance results are shown in Tables A.1−A.3. The chemical abundances of the HARPS sample were obtained from Adibekyan et al. (2012c) where they use exactly the same tools and model atmospheres.
6. Discussion
Gonzalez (1997); Gonzalez & Laws (2000) already noticed that giant planet hosts tend to be more metal-rich than stars without detected planets. Santos et al. (2001) provided supporting evidence of a metal-rich origin of giant-planet host stars and following studies confirmed this result (e.g., Santos et al. 2004, 2005; Valenti & Fischer 2005). Recent studies show that Neptune and super-Earth class planet hosts have a different metallicity distribution, more similar to stars without planets (e.g., Udry et al. 2006; Sousa et al. 2008, 2011; Ghezzi et al. 2010; Mayor et al. 2011; Buchhave et al. 2012).
In this section, we inspect the abundance ratios of different elements, [X/ Fe], as a function of the metallicity for our BD-companion stellar sample, and compare them with the planetary-companion stars analized by Adibekyan et al. (2012c). We also study the distributions of different element abundances in different samples and we compare them with other stars with and without planets as a function of the minimum mass of the most massive companion of each host star, mC sin i.
6.1. Galactic abundance trends
The abundances of the refractory elements in the CORALIE subsample exhibit similar behavior to other stars with and without planets analyzed in previous works (Neves et al. 2009; Adibekyan et al. 2012c) (see Fig. B.1).
Stars with BDs follow the galactic abundance trend except for some particular elements (Co, Si, Sc) where the abundances are slightly lower than expected. These exceptions may be due to the small number of lines needed to achieve a reliable mean abundance in certain elements (e.g., ScI for HD 167665). Stars with BD companions appear to be located at an intermediate range of metallicity, between stars with and without planets.
In Fig. 5 we display the mean abundance ratio of the α-elements Mg, Si, Ca, and Ti (with [Xα/H] computed as the sum of individual element abundances [X/H] divided by 4, and [Xα/Fe] = [Xα/H] − [Fe/H]) against the metallicity of stars with small and giant planets from (Adibekyan et al. 2012b) together with stars with BDs. The range of metallicities of stars with confirmed BDs seems to be narrower than that of the planet hosts, although this may not be statistically significant due to the small number of stars in the BD sample. Adibekyan et al. (2012b,a) remarked that stars with both small and giant planets at low metallicities [Fe/H] < − 0.3 dex tend to be α-enhanced and therefore to belong chemically to thick-disk populations. There is only one star with BD at these relatively low metallicities, at [Fe/H] ~ − 0.35 dex, showing a relatively low [Xα/Fe] ratio, and therefore the question whether this behaviour still holds for BD hosts remains open. The α-element abundance ratios [Xα/ Fe] of BD hosts seem to be consistent with the Galactic trend at higher metallicities, although there are some stars at metallicities below solar with relatively low [Xα/Fe] ratios, even below the trend described by the stars without detected planets of the HARPS sample. Stars with discarded BDs and without detected BD candidates also follow the general abundance trend.
![]() |
Fig. 5 α-element abundance ratios [Xα/ Fe] versus [Fe/H] for the samples of stars without planets (green empty squares), with small planets (black empty triangles), with giant planets (red empty triangles), BD-host stars (blue filled circles), stars without known BD companions from the CORALIE sample (blue empty circles), and stars with discarded BD candidates (violet filled circles). |
6.2. Element abundance distributions
The histograms displayed in Fig. B.2 allow us to study the abundance distribution [X/H] of our sample for each element. Stars without planets have a maximum abundance at ~−0.1 for most species. Stars with giant planets exhibit a more metal-rich maximum at ~0.2 − 0.3 dex, while the stars with small planets (whose maximum is at ~−0.1 dex) resemble the stars without planets. This general behaviour has already been noticed in previous papers (Neves et al. 2009; Adibekyan et al. 2012b) and agrees well with the so-called metallicity effect, i.e., the strong correlation between stellar metallicity and the likelihood of finding giant planets (Santos et al. 2001).
The abundance distribution of the sample of confirmed BD-host stars appears to be located in between the stars with small planets and stars with giant planets. Some element distributions (Na, Si, Mg, Mn) are more similar to those of giant planets, whereas for others elements (Ti, Cr, Co, Ni) the behaviour is closer to the stars without planets (see Fig. B.2).
![]() |
Fig. 6 Histogram of α-element abundances, [Xα/ H], (top panel) of the samples of stars without planets (green continuous line), stars with small planets (black dashed-triple-dotted line), stars with giant planets (red dashed-dotted line), stars with confirmed BD companions (blue dashed line), and stars with discarded BD candidates, i.e., binaries with low-mass M-dwarf companions (violet dotted line). The left y-axis in the top panel is labeled with the number of stars with and without small and giant planets, whereas the right y-axis shows the number of stars with BD-companion candidates. The lower panel shows cumulative histogram. |
In Figs. 6 and 7 we display the distribution and cumulative histograms of the α-element abundances (including the species Mg i, Si i, Ca i, Ti i), [Xα/ H], and the Fe-peak element abundances (including the species Cr i, Mn i, Co i, Ni i), [XFe/ H]. Here it is clear that confirmed BDs behave differently from giant planets. Although the sample is small, it may tentatively point to a bimodal distribution with two peaks, one at the position of the small planet distribution and one at the position of the giant planets. However, the mean values of the [Xα/ H] and [XFe/ H] abundances are roughly solar. The [Xα/ H] and [XFe/ H] cumulative histograms support the previous statement. Stars with small planets and without detected planets go together, whereas the stars with BDs exhibit a slightly different behaviour with a later growth of the cumulative histogram which resembles that of stars with small planets only at [X/H] > 0.1 dex. The cumulative histogram of stars with giant planets clearly manifest a later increase towards high metallicities reaching the saturation at [X/H] ~ 0.3 dex. We perform a K-S test to statistically evaluate the significance of this apparent different behaviour (see Table 3). This test provides a clear difference between the BD-host sample and the GP sample, but the SP and NP sample seems to be statistically very similar to the BD sample. The number of stars with confirmed BDs must be increased in order to be able to distinguish these populations. On the other hand, in Table 3 we also show the same K-S test for the stars with discarded BDs which are binaries hosting low-mass M dwarfs. Although again there are only six stars in this sample, it appears to be statistically different from all the samples of GP, SP, and NP, especially for the Fe-peak element abundances. However, the significance is lower than the comparison with confirmed BD-host stars.
Significance of the K-S test of the BD sample.
6.3. Abundance ratios [X/H] against companion mass
The mC sin i values of the BD candidates in our sample are higher than 13 MJ (the star with the lightest BD companion analyzed in this work is orbiting HD 211847, mC sin i ~ 19.2 MJ), but two of the most massive giant planets of HARPS sample exceed this value. The boundary between brown dwarfs and giant planets has been extensively investigated in the literature, and some works agree with the definition that brown dwarfs can burn the deuterium that is present when they form, and giant planets cannot (e.g., Burrows et al. 1997). In Bodenheimer et al. (2013), the borderline between giant planets and brown dwarfs is found to depend only slightly on different parameters, such as core mass, stellar mass, formation location, solid surface density in the protoplanetary disk, disk viscosity, and dust opacity. More than 50% of the initial deuterium is burned for masses above 11.6 − 13.6 MJ, in agreement with previous determinations that do not take the formation process into account. Thus, we use the mass ~13MJ to distinguish between giant planets and brown dwarfs.
The metal content of stars at birth surely affects the formation of the planetary-mass companions, but is this also true for stars with BD companions? To answer this question we show in Fig. B.3 the chemical element abundances, [X/H], as a function of the minimum mass of the most massive substellar companion, which could be a small planet, a giant planet, or a brown dwarf according to their mC sin i values. Qualitatively, these element abundances seem to progressively increase with the companion mass from small planets until reaching a maximum at about 1 MJ and then slightly decrease when entering the BD regime. The scatter in [X/H] may be due to the different intrinsic metallicities of the stars at every bin in mC sin i.
In Fig. B.4 we display the mean values of these element abundances, [X/H], in bins of mC sin i. The standard deviations from the mean values are different for every element, but they are around 0.15 − 0.2 dex. These mean element abundances keep roughly constant for small planets with masses lower than 0.04 MJ. From this point on, the abundances grow with the companion mass even in the giant-planet companion range, from low-mass up to Jupiter-mass planets, reaching a maximum in ~0.8 MJ (see e.g., Si i, Ti i, Cr i, and Ni i). For more massive giant-planet companions, these abundances start to decrease slowly with the companion mass towards high-mass BD companions. The stars with BD companions just follow the decreasing trend of the stars with giant planets.
![]() |
Fig. 8 Abundance ratios [Xα/ H] and [XFe/H] versus minimum mass of the most massive companion mC sin i of α-elements (top-left panel) and iron-peak elements (top-right panel), for stars with small planets (black empty triangles), with giant planets (red empty triangles), with confirmed BD companions (blue filled circles). The average values are shown as violet three-dotted-dashed lines. Lower panels show the mean abundances of stars with small and giant planets, and confirmed BD companions, in equal-sized bins appropriated for the logarithmic scale of companion masses as in Fig. B.4. Error bars represents the standard deviation of the mean divided by the square root of the number of stars in each bin. The green solid line depicts the parabolic fit of the data, whereas the violet three-dotted-dashed lines show the zero-order fits. |
We also decided to use the mean values of the α-elements, [Xα/ H], and Fe-peak elements, [XFe/ H], of each star in these samples. Thus, in Fig. 8 we display these abundances as a function of the minimum mass of the most massive companion, mC sin i. The scatter of these abundances is high because of the different global metallicity of different stars. The average values of these abundances are shown in Table 4 for the SP, GP, and BD samples. One can see in Fig. 8 the different levels of these three samples, although the BD sample shows an average value consistent with the SP sample within the error bars (see Table 4).
Mean abundance of SP, GP, and BD samples.
In the lower panels of Fig. 8, we depict the weighted average of these abundances at each mass bin and there the trend is clearer. We fit a parabolic function, using the IDL routine curvefit, to these mean values of the α-element and Fe-peak element abundances. The peaks of these trends have a maximum of abundance at ~0.15 ± 0.01 dex and companion mass of mC sin i ~ 1.42 ± 0.17 and 1.32 ± 0.18 MJ, respectively. The parabolic fit, with values of 3.2 and 2.3, respectively, provides a better representation of the average abundances than the linear fit (with
and 8.6). We perform an F-test to confirm that the parabolic curve fits this data set better than the linear fit. We use the IDL routine mpftest from the Markwardt library6. The result reveals a significance level of 1 × 10-5 for α-elements (F = 45.8) and 3 × 10-6 for iron-peak elements (F = 37.9), implying that the parabolic fit is significantly better than the linear.
We also perform a fit of a three zero-order function of three levels describing the SP, GP and BD samples and the values are given in Table 4, with values of 3.0 and 2.6 for α-element and Fe-peak element abundances, respectively. We compare this three-step model with the parabolic fit using an F-test, resulting in F values of −0.5 and 0.6 which gives a significance level of 0 and 0.5 for α-elements and Fe-peak elements, respectively. Therefore, the 3-step model provides a similar description of the data than the parabolic fit. We also check that a two-step model provides a worse fit than the three-step model (with
and 3.1 for α-elements and Fe-peak elements, respectively).
Finally, Sahlmann et al. (2011) noticed that there is a lack of BD companions with masses in the range mC sin i ~ 35−55 MJ. More recently, Ma & Ge (2014) have collected the known BD companions from different studies and confirm this gap for stars with periods shorter than 100 days. Although the statistics are still poor, these authors suggest that BD companions below this gap, i.e., with mC sin i< 42 MJ, may have formed in protoplanetary disks as giant planets, probably through the disk instability-fragmentation mechanism (Boss 1997; Stamatellos & Whitworth 2009), whereas BD companions with mC sin i> 42 MJ may have formed by molecular cloud fragmentation as stars (Padoan & Nordlund 2004; Hennebelle & Chabrier 2008).
In the bottom panels of Fig. 8 we can see that on average stars with BDs at masses mC sin i below 42 MJ have higher α-element and Fe-peak element abundances than stars above this mass. Stars with massive BD companions have more similar abundances to those of stars without planets. This might tentatively support the above statement with the low-mass BD companions being formed in protoplanetary disks as giant planets, and the high-mass BDs being formed by cloud fragmentation as stars.
7. Conclusions
We have analyzed a subsample of stars with candidate BD companions from the CORALIE radial velocity survey. We derive chemical abundances of several elements including α-elements and Fe-peak elements. A comparison with the chemical abundances of stars with giant planets shows that BD-host stars seem to behave differently. In particular, we compute the abundance histograms [Xα/H] and [XFe/H], revealing a mean abundance at about solar for the BD-host sample, whereas for stars without planets (NP) and with small planets (SP) it remains at −0.1 dex, and for stars with giant planets (GP) at roughly + 0.10 dex. The cumulative histograms of [Xα/H] and [XFe/H] abundances exhibit the same situation, with the stars without planets and with small planets going together, similarly to the stars with BDs. However, the stars with giant planets reach a later saturation at [X/H] ~ 0.3 dex. A Kosmogorov-Smirnov (K-S) test does not show a statisticallly significant difference between the cumulative distribution of SP, NP, and BD samples, but clearly separates the GP and BD samples.
Finally, we depict the [Xα/H] and [XFe/H] abundances versus the minimum mass of the most massive substellar companion, mC sin i, and we find a peak of these element abundances for a companion mass mC sin i ~ 1.3 − 1.4 MJ, with the abundances growing with the companion mass from small planets to Jupiter-like planets and after decresing towards massive BD companions. A three-step model also provides a similar description of the data with no statistically significant difference with the parabolic model. Recently, Sahlmann et al. (2011) and Ma & Ge (2014) have suggested that the formation mechanism may be different for BD companion below and above 42 MJ. We find that BDs below this mass tend to have higher abundances than those above this mass, which may support this conclusion and BDs with mC sin i< 42 MJ may form by disk instability-fragmentation, whereas high-mass BDs may form as stars by cloud fragmentation.
Online material
Appendix A: Tables of chemical abundances
In this Appendix we provide all the tables containing the element abundances of the 17 stars within the CORALIE sample, the first 13 stars with BD-companion candidates (7 confirmed BDs and 6 discarded BDs; for details see Sect. 3) and 4 stars as a comparison sample. We also provide tables with the minimum mass of the most massive substellar companion of stars in the HARPS sample.
Abundances of the CORALIE stellar sample.
Abundances of the CORALIE stellar sample.
Abundances of the CORALIE stellar sample.
Minimum mass of the most massive substellar companion of stars in the HARPS sample.
Minimum mass of the most massive substellar companion of stars in the HARPS sample.
![]() |
Fig. B.1 Element abundance ratios [X/Fe] versus [Fe/H] for the samples of stars without planets (green empty squares), stars with planets (red empty triangles), BD-host stars (blue filled circles), stars without known BD companions from the CORALIE sample (blue empty circles), and stars with discarded BD candidates (violet filled circles). At the top-left corner of each panel, the mean standard deviation of the abundance measurements are shown for the stars without planetary companion (σNP), stars with planets (σWP), and stars with confirmed BD companions (σBD). |
![]() |
Fig. B.2 Histogram of element abundances, [X/H], of the samples of stars without planets (green continuous line), stars with low-mass planets (black dashed-triple-dotted line), stars with giant planets (red dashed-dotted line), stars with confirmed BD companions (blue dashed line), and stars with discarded BD candidates (violet dotted line). The left y-axis of the top panel is labeled with the number of stars with and without small and giant planets, whereas the right y-axis shows the number of stars with BD companion candidates. |
![]() |
Fig. B.3 Element abundances, [X/H], versus the minimum mass of the most massive companion, mC sin i. Black empty triangles refer to stars with small planets, red empty triangles to stars with giant planets, blue filled circles to stars with confirmed BDs, and violet filled circles to stars with discarded BD candidates. At the bottom of each panel, standard deviations of the mean element abundances are shown for stars without planets (σNP), stars with planets (σWP), and stars with confirmed BD companions (σBD). |
![]() |
Fig. B.4 Mean element abundances versus the minimum mass of the most massive companion, mC sin i, computed in equal-size bins appropriated for the logarithmic scale of companion masses. Error bars represent the standard deviation of the mean. |
The ARES code can be downloaded at: http://www.astro.up.pt/
The TAME code can be downloaded at: http://astro.snu.ac.kr/~wskang/tame/
The MOOG code can be downloaded at: http://www.as.utexas.edu/~chris/moog.html
The Markwardt IDL library can be downloaded at: http://cow.physics.wisc.edu/~craigm/idl/idl.html
Acknowledgments
D.M.S. is grateful to the Spanish Ministry of Education, Culture and Sport for the financial support from a collaboration grant, and to the Ph.D. contract funded by Fundación La Caixa. J.I.G.H. and G.I. acknowledge financial support from the Spanish Ministry project MINECO AYA2011-29060, and J.I.G.H. also from the Spanish Ministry of Economy and Competitiveness (MINECO) under the 2011 Severo Ochoa Program MINECO SEV-2011-0187. N.C.S. acknowledges the support by the European Research Council/European Community under the FP7 through Starting Grant agreement number 239953, and the support in the form of a Investigador FCT contract funded by Fundação para a Ciência e a Tecnologia (FCT) /MCTES (Portugal) and POPH/FSE (EC). This research has made use of the IRAF facilities, and the SIMBAD database, operated at the CDS, Strasbourg, France.
References
- Adibekyan, V. Z., Delgado Mena, E., Sousa, S. G., et al. 2012a, A&A, 547, A36 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Adibekyan, V. Z., Santos, N. C., Sousa, S. G., et al. 2012b, A&A, 543, A89 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Adibekyan, V. Z., Sousa, S. G., Santos, N. C., et al. 2012c, A&A, 545, A32 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bodenheimer, P., D’Angelo, G., Lissauer, J. J., Fortney, J. J., & Saumon, D. 2013, ApJ, 770, 120 [NASA ADS] [CrossRef] [Google Scholar]
- Boss, A. P. 1997, Science, 276, 1836 [NASA ADS] [CrossRef] [Google Scholar]
- Buchhave, L. A., Latham, D. W., Johansen, A., et al. 2012, Nature, 486, 375 [NASA ADS] [Google Scholar]
- Burrows, A., Marley, M., Hubbard, W. B., et al. 1997, ApJ, 491, 856 [NASA ADS] [CrossRef] [Google Scholar]
- Chabrier, G. 2002, ApJ, 567, 304 [NASA ADS] [CrossRef] [Google Scholar]
- Ge, J., Mahadevan, S., Lee, B., et al. 2008, in Extreme Solar Systems, eds. D. Fischer, F. A. Rasio, S. E. Thorsett, & A. Wolszczan, ASP Conf. Ser., 398, 449 [Google Scholar]
- Ghezzi, L., Cunha, K., Smith, V. V., et al. 2010, ApJ, 720, 1290 [NASA ADS] [CrossRef] [Google Scholar]
- Gonzalez, G. 1997, MNRAS, 285, 403 [NASA ADS] [CrossRef] [Google Scholar]
- Gonzalez, G., & Laws, C. 2000, AJ, 119, 390 [NASA ADS] [CrossRef] [Google Scholar]
- Grether, D., & Lineweaver, C. H. 2006, ApJ, 640, 1051 [NASA ADS] [CrossRef] [Google Scholar]
- Halbwachs, J. L., Mayor, M., Udry, S., & Arenou, F. 2003, A&A, 397, 159 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Hayashi, C., & Nakano, T. 1963, Prog. Theor. Phys., 30, 460 [NASA ADS] [CrossRef] [Google Scholar]
- Hennebelle, P., & Chabrier, G. 2008, ApJ, 684, 395 [NASA ADS] [CrossRef] [Google Scholar]
- Kang, W., & Lee, S.-G. 2012, MNRAS, 425, 3162 [NASA ADS] [CrossRef] [Google Scholar]
- Kirkpatrick, J. D., Gelino, C. R., Cushing, M. C., et al. 2012, ApJ, 753, 156 [NASA ADS] [CrossRef] [Google Scholar]
- Kumar, S. S. 1962, AJ, 67, 579 [NASA ADS] [CrossRef] [Google Scholar]
- Kurucz, R. 1993, ATLAS9 Stellar Atmosphere Programs and 2 km s-1 grid. Kurucz CD-ROM No. 13 (Cambridge, Mass.: Smithsonian Astrophysical Observatory), 13 [Google Scholar]
- Ma, B., & Ge, J. 2014, MNRAS, 439, 2781 [NASA ADS] [CrossRef] [Google Scholar]
- Mayor, M., & Queloz, D. 1995, Nature, 378, 355 [NASA ADS] [CrossRef] [Google Scholar]
- Mayor, M., Pepe, F., Queloz, D., et al. 2003, The Messenger, 114, 20 [NASA ADS] [Google Scholar]
- Mayor, M., Marmier, M., Lovis, C., et al. 2011 [arXiv:1109.2497] [Google Scholar]
- Marcy, G. W., & Butler, R. P. 2000, PASP, 112, 137 [NASA ADS] [CrossRef] [Google Scholar]
- Nakajima, T., Oppenheimer, B. R., Kulkarni, S. R., et al. 1995, Nature, 378, 463 [NASA ADS] [CrossRef] [Google Scholar]
- Neves, V., Santos, N. C., Sousa, S. G., Correia, A. C. M., & Israelian, G. 2009, A&A, 497, 563 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Padoan, P., & Nordlund, Å. 2004, ApJ, 617, 559 [NASA ADS] [CrossRef] [Google Scholar]
- Perryman, M. A. C., Lindegren, L., Kovalevsky, J., et al. 1997, A&A, 323, L49 [NASA ADS] [Google Scholar]
- Pollack, J. B., Hubickyj, O., Bodenheimer, P., et al. 1996, Icarus, 124, 62 [NASA ADS] [CrossRef] [Google Scholar]
- Rebolo, R., Martin, E. L., & Magazzu, A. 1992, ApJ, 389, L83 [NASA ADS] [CrossRef] [Google Scholar]
- Rebolo, R., Zapatero Osorio, M. R., & Martín, E. L. 1995, Nature, 377, 129 [NASA ADS] [CrossRef] [Google Scholar]
- Sahlmann, J., Ségransan, D., Queloz, D., et al. 2011, A&A, 525, A95 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Santos, N. C., Israelian, G., & Mayor, M. 2001, A&A, 373, 1019 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Santos, N. C., Israelian, G., & Mayor, M. 2004, A&A, 415, 1153 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Santos, N. C., Israelian, G., Mayor, M., et al. 2005, A&A, 437, 1127 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Sneden, C. 1973, ApJ, 184, 839 [NASA ADS] [CrossRef] [Google Scholar]
- Sousa, S. G., Santos, N. C., Israelian, G., Mayor, M., & Monteiro, M. J. P. F. G. 2007, A&A, 469, 783 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Sousa, S. G., Santos, N. C., Mayor, M., et al. 2008, A&A, 487, 373 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [Google Scholar]
- Sousa, S. G., Santos, N. C., Israelian, G., Mayor, M., & Udry, S. 2011, A&A, 533, A141 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Stamatellos, D., & Whitworth, A. P. 2009, MNRAS, 392, 413 [Google Scholar]
- Udry, S., & Santos, N. C. 2007, ARA&A, 45, 397 [NASA ADS] [CrossRef] [Google Scholar]
- Udry, S., Mayor, M., Queloz, D., Naef, D., & Santos, N. 2000, in From Extrasolar Planets to Cosmology: The VLT Opening Symposium, eds. J. Bergeron, & A. Renzini (Springer Verlag), 571 [Google Scholar]
- Udry, S., Mayor, M., Benz, W., et al. 2006, A&A, 447, 361 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Valenti, J. A., & Fischer, D. A. 2005, in Protostars and Planets V, Proc. Conf., in Hilton Waikola, 8592 [Google Scholar]
- van Leeuwen, F. 2007, A&A, 474, 653 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
All Tables
Minimum mass of the most massive substellar companion of stars in the HARPS sample.
Minimum mass of the most massive substellar companion of stars in the HARPS sample.
All Figures
![]() |
Fig. 1 Histograms of the stellar parameters Teff, log g, and [Fe/H] of our CORALIE sample. |
In the text |
![]() |
Fig. 2 ARES versus TAME graphics for the stars HD 89707 (S/N ~ 110) and HD 206505 (S/N ~ 70). Left panels: EW measured with TAME versus ARES. The 1:1 correspondence is shown as a dashed line. EW differences, EWTAME − EWARES, versus EWARES (middle panels) and spectral line wavelength (right panels). Dash-dotted lines define the mean value of the data points, and dotted lines define the mean plus the standard deviation. |
In the text |
![]() |
Fig. 3 Minimum mass of the most massive companion, mC sin i, and orbital period, PC, versus the metallicity of stars with BDs. The CORALIE sample is depicted as filled symbols, and the sample in Ma & Ge (2014) is displayed as empty symbols, separated between giant and dwarf stars. |
In the text |
![]() |
Fig. 4 High-resolution normalized spectra of late G-, K-type stars: CORALIE spectra of the coolest stars in the sample (black), HARPS spectra (red) of cool stars from the HARPS database (e.g., Sousa et al. 2008) and the CORALIE spectrum (blue) of HD 53680, and the HARPS spectrum (blue) of HIP 103019. The spectra are depicted following a sequence with decreasing Tefffrom top to bottom. |
In the text |
![]() |
Fig. 5 α-element abundance ratios [Xα/ Fe] versus [Fe/H] for the samples of stars without planets (green empty squares), with small planets (black empty triangles), with giant planets (red empty triangles), BD-host stars (blue filled circles), stars without known BD companions from the CORALIE sample (blue empty circles), and stars with discarded BD candidates (violet filled circles). |
In the text |
![]() |
Fig. 6 Histogram of α-element abundances, [Xα/ H], (top panel) of the samples of stars without planets (green continuous line), stars with small planets (black dashed-triple-dotted line), stars with giant planets (red dashed-dotted line), stars with confirmed BD companions (blue dashed line), and stars with discarded BD candidates, i.e., binaries with low-mass M-dwarf companions (violet dotted line). The left y-axis in the top panel is labeled with the number of stars with and without small and giant planets, whereas the right y-axis shows the number of stars with BD-companion candidates. The lower panel shows cumulative histogram. |
In the text |
![]() |
Fig. 7 Same as Fig. 6, but for Fe-peak element abundances. |
In the text |
![]() |
Fig. 8 Abundance ratios [Xα/ H] and [XFe/H] versus minimum mass of the most massive companion mC sin i of α-elements (top-left panel) and iron-peak elements (top-right panel), for stars with small planets (black empty triangles), with giant planets (red empty triangles), with confirmed BD companions (blue filled circles). The average values are shown as violet three-dotted-dashed lines. Lower panels show the mean abundances of stars with small and giant planets, and confirmed BD companions, in equal-sized bins appropriated for the logarithmic scale of companion masses as in Fig. B.4. Error bars represents the standard deviation of the mean divided by the square root of the number of stars in each bin. The green solid line depicts the parabolic fit of the data, whereas the violet three-dotted-dashed lines show the zero-order fits. |
In the text |
![]() |
Fig. B.1 Element abundance ratios [X/Fe] versus [Fe/H] for the samples of stars without planets (green empty squares), stars with planets (red empty triangles), BD-host stars (blue filled circles), stars without known BD companions from the CORALIE sample (blue empty circles), and stars with discarded BD candidates (violet filled circles). At the top-left corner of each panel, the mean standard deviation of the abundance measurements are shown for the stars without planetary companion (σNP), stars with planets (σWP), and stars with confirmed BD companions (σBD). |
In the text |
![]() |
Fig. B.2 Histogram of element abundances, [X/H], of the samples of stars without planets (green continuous line), stars with low-mass planets (black dashed-triple-dotted line), stars with giant planets (red dashed-dotted line), stars with confirmed BD companions (blue dashed line), and stars with discarded BD candidates (violet dotted line). The left y-axis of the top panel is labeled with the number of stars with and without small and giant planets, whereas the right y-axis shows the number of stars with BD companion candidates. |
In the text |
![]() |
Fig. B.3 Element abundances, [X/H], versus the minimum mass of the most massive companion, mC sin i. Black empty triangles refer to stars with small planets, red empty triangles to stars with giant planets, blue filled circles to stars with confirmed BDs, and violet filled circles to stars with discarded BD candidates. At the bottom of each panel, standard deviations of the mean element abundances are shown for stars without planets (σNP), stars with planets (σWP), and stars with confirmed BD companions (σBD). |
In the text |
![]() |
Fig. B.4 Mean element abundances versus the minimum mass of the most massive companion, mC sin i, computed in equal-size bins appropriated for the logarithmic scale of companion masses. Error bars represent the standard deviation of the mean. |
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.