A&A 462, 345-353 (2007)
DOI: 10.1051/0004-6361:20066319
S. Desidera1 - M. Barbieri1,2
1 - INAF - Osservatorio Astronomico di Padova,
Vicolo dell' Osservatorio 5, 35122 Padova, Italy
2 -
Dipartimento di Fisica, Università di Padova, Italy
Received 29 August 2006 / Accepted 1 October 2006
Abstract
Aims. The statistical properties of planets in binaries were investigated. Any difference to planets orbiting single stars can shed light on the formation and evolution of planetary systems. As planets were found around components of binaries with very different separation and mass ratio, it is particularly important to study the characteristics of planets as a function of the effective gravitational influence of the companion.
Methods. A compilation of planets in binary systems was made; a search for companions orbiting stars recently shown to host planets was performed, resulting in the addition of two further binary planet hosts (HD 20782 and HD 109749). The probable original properties of the three binary planet hosts with white dwarfs companions were also investigated. Using this updated sample of planets in binaries we performed a statistical analysis of the distributions of planet mass, period, and eccentricity, fraction of multiplanet systems, and stellar metallicity for planets orbiting components of tight and wide binaries and single stars.
Results. The only highly significant difference revealed by our analysis concerns the mass distribution of short-period planets. Massive planets in short period orbits are found in most cases around the components of rather tight binaries. The properties of exoplanets orbiting the components of wide binaries are compatible with those of planets orbiting single stars, except for a possible greater abundance of high-eccentricity planets. The previously suggested lack of massive planets with P>100 days in binaries is not confirmed.
Conclusions. We conclude that the presence of a stellar companion with separation smaller than 100-300 AU is able to modify the formation and/or migration and/or the dynamical evolution history of giant planets while wide companions play a more limited role.
Key words: stars: planetary systems - stars: binaries: visual - stars: general - stars: individual: HD 20782, HD 109749, GL 86
The study of the frequency of planets in binary systems and of the
properties of these planets and of the binary systems themselves
is very important to improve our knowledge of planet formation and evolution.
On one hand, the frequency of planets
in binary systems has a strong effect on the global frequency of planets,
more than half of solar type stars being in binary or multiple systems
(Duquennoy & Mayor 1991).
On the other hand, the properties of planets in binaries, and any
difference to those of the planets orbiting single stars, would
shed light on the effects caused by the presence of the companions.
The search for planets in binaries can follow two complementary approaches. The first one is to perform dedicated surveys looking for planets in binary systems. Several programs are currently in progress focusing on different types of binaries. The SARG planet search has been studying about 50 wide pairs with similar components (Desidera et al. 2006) since late 2000. A similar study was recently started by the Geneva group using HARPS. About 100 spectroscopic binaries with orbital periods between 1.5 and 100 yr were monitored by Eggenberger et al. (2006a), while Konacki (2005a) is monitoring about 450 stars in binaries.
The second approach is to study the binarity of the hosts of planets discovered in general surveys, which include many binary stars in their lists in spite of some selection biases against them. Dedicated adaptive optics observations are required to reveal faint companions close to the star (e.g. Patience et al. 2002; Chauvin et al. 2006; Mugrauer et al. 2006a; see Table 1 for further references). However, the search for brighter companions at large separation can be performed by checking astronomical catalogs for nearby, comoving objects.
The first analysis of the properties of planets in binaries revealed differences with respect to those of planets orbiting single stars. Zucker & Mazeh (2002) showed some difference in the period-mass relation. A more complete statistical analysis by Eggenberger et al. (2004), hereafter E04, confirms that massive planets in close orbits are mostly found in binaries. They also propose other possible peculiar characteristics of planets in binaries: low eccentricities for orbital periods shorter than 40 days and a lack of massive planets with periods longer than 100 days in binaries. The estimated significance of these latter features is less than 97%, therefore requiring larger samples for confirmation.
The continuing observational efforts of planet searches and
dedicated follow-up to study the characteristics of planet hosts
lead to a significant increase of the total number of planets
in binaries.
We now know more than 40 planets in binaries or multiple systems,
an increase by more than a factor of two with respect to the
first compilation of the planets in binary systems assembled by
E04 (19 planets in 15 systems).
The recent publication of the Catalog of Nearby Exoplanets by Butler et al. (2006), hereafter B06, provides an updated and homogeneous assembly of planet and stellar properties (excluding binarity). This work is a useful starting point for a statistical comparison of the properties of planets orbiting the components of multiple systems and single stars. In this paper we present such a discussion.
The outline of the paper is the following. In Sect. 2 the samples of planets in binaries and in single stars are defined; in Sect. 3 the statistical analysis of the properties of planets in binaries and single stars is presented. In Sect. 4 the results and their implications are discussed. Appendices A-C address specific issues relevant for individual cases: in Appendix A we present the results of the search for companions of stars recently shown to host planets, resulting in two further binaries (HD 20782 and HD 109749), whose properties are discussed is detail; in Appendix B we discuss the probable original properties (mass ratio, separation) of the three systems including white dwarf (hereafter WD) companions; in Appendix C we provide details about individual objects.
We assembled the properties of the planets in multiple systems
and of the multiple system themselves (Table 1,
only available in electronic form).
For the planet hosts not included in the latest
study of multiplicity of planet hosts
(Raghavan et al. 2006, hereafter R06),
we performed a search in existing catalogs of binaries.
Details are presented in Appendix A.
The planetary properties, as well as mass and metallicity of the planet hosts,
were taken from B06.
The binary parameters were taken from the references listed in the
Table caption.
As we are interested in evaluating the actual dynamical effects of the
companions,
we listed individual masses, derived using the mass-luminosity relations by
Delfosse et al. (2000) or taken from the works cited in the
caption.
In a few cases, the complete
binary orbit is available, making feasible a more in depth analysis
of the possible formation and evolution of the planet
(see e.g. Thebault et al. 2004, for Cep).
The last column lists the critical semiaxis for dynamical stability
of planetary orbit derived
using Eq. (1) of Holman & Weigert (1999).
When the full binary orbit is not available,
this was calculated deriving the semimajor axis from the projected
separation and adopting an eccentricity of 0.35
(Fischer & Marcy 1992; Duquennoy & Mayor 1991).
Remarks on individual
objects are reported in Appendix C, while the properties
of the systems containing a white dwarf are discussed in
Appendix B.
The comparison sample of planets in single stars is that of B06, excluding the
objects in Table 1.
The mass limit adopted by B06
(
)
is different to those
adopted in R06.
One "super-planet'' in the range
was found in a binary,
the companion orbiting HD 41004B.
Figure 1 shows the mass ratio vs. semimajor axis for stars with planets in multiple systems. For hierarchical triple systems in which the planet orbits the isolated companion, the masses of the binary companions to the planet host are summed. It appears that planets might exist in binaries with very different properties. In some cases (e.g. very low mass companions at a projected separation larger than 1000 AU) the dynamical effects of the companion on the formation and evolution of the planetary system might be very limited, while in the cases of very tight binaries the presence of the planet represents a challenge for the current models of planet formation (Hatzes & Wuchterl 2005).
A simple look at Fig. 1 suggests a possible depopulated region for the binary semimajor axis between about 20 to 100 AU. A proper evaluation of its reality requires a detailed study of the selection effects of the original samples of radial velocity surveys and the discovery of stellar companions of planet hosts, which is postponed for a future study.
The dynamical effects of the companions on the circumstellar region
of the planet hosts are certainly very different for the binaries
in our sample.
To consider the effects of dynamical perturbation by the
stellar companion(s) we used the critical semiaxis for dynamical
stability of the planet
(Holman & Wiegert 1999).
This allows us to take into account the effects
of both the separation and the mass ratio.
The critical semiaxis for dynamical stability represents the limit
for stable planetary orbits around the planet hosts (calculated
for coplanar circular orbits).
It is also similar to the radius of tidal truncation of the circumstellar disk
(Pichardo et al. 2005; Pfahl & Muterspaugh 2006).
However, the area in which relative
impact velocities between planetesimals is expected to allow planet formation
might be significantly smaller than
(Thebault et al. 2006).
![]() |
Figure 1: Mass ratio vs. semimajor axis of the binary orbit for stars with planets in binary systems. Open circles represent the pairs for which binary orbit is available, open squares the pairs for which only the binary separation is available. |
Open with DEXTER |
The critical semimajor axis
was used to divide the sample according
to the relevance of the dynamical effects.
We define as "tight'' binaries those with
AU
and "wide'' binaries those with
AU. The limit corresponds to
a projected separation of about 200-300 AU depending on the mass ratio.
The choice was fixed as a compromise, to obtain
a sample of binaries wide enough to guarantee limited
dynamical perturbations for the portion of the disk
on which the formation of giant planets
is expected to occur
and to allow the inclusion of an adequate number of objects in the
"tight'' binary sample.
When relevant for the statistical analysis, the effect of changing
this limit is considered in the following discussion.
The statistical comparison to test the hypothesis that the parameters of planets (mass, period, eccentricity) in tight and wide binaries and in single stars can be drawn from the same parent distribution was estimated using the Kolmogorov-Smirnov (hereafter KS) test and the Mann-Whitney U (hereafter MWU) test (Babu & Feigelson 1996). We also used the hypergeometric distribution (see E04), that yields the probability of obtaining by chance a number of cases from a subsample of given size. The statistical analysis was performed including only planets with radial velocity (hereafter RV) semi-amplitude K>15 m/s, to exclude planets with uncertain orbital parameters and to ensure a more homogeneous detectability. This choice has little effect on the results.
Zucker & Mazeh (2002) showed that planets in binaries follow a different mass period relation than those orbiting single stars. In particular, high-mass, short-period planets appear to be present only in binary systems. Such a correlation was confirmed by E04 but further inclusion of new systems makes it weaker (Mugrauer et al. 2005). The addition of further low-mass hot-Jupiters in binaries from our work consolidates this latter trend.
The possibility that the Zucker & Mazeh correlation was simply due to the small number of planet in binaries known at that time can be considered. However, an alternative interpretation appears more convincing when one takes into account the binary properties of the massive hot-Jupiter hosts: these are mostly relatively tight binaries (Fig. 2). Indeed, Zucker & Mazeh (2002) considered in their statistical analysis only binaries whose projected separation is smaller than 1000 AU, then excluding very wide binaries included in increasingly larger number in later studies.
![]() |
Figure 2: Projected mass vs. orbital period of extrasolar planets. Open circles: single stars; filled circles: binary stars. The size of the symbol is proportional to the critical semimajor axis for dynamical stability (larger symbols refer to the tighter binaries). |
Open with DEXTER |
To test this hypothesis statistically, we selected from
Table 1 and B06 the planets with
P<40 days (about in the middle of the so-called "period
valley'' on the period distribution of planets, see
Udry et al. 2003 and Fig. 5).
We performed a KS test comparing
the mass distribution of planets with P<40 days orbiting
single stars, binaries with
75 AU and
binaries with
75 AU (34, 5, and 6 planets respectively)
(Fig. 3).
The hypothesis that the mass of close-in planets in tight
binaries can be drawn from the same parent distribution of close-in planets
orbiting members of wide binaries or single stars can be rejected
with a probability of 98.2% and 99.5%.
The MWU test yields an even higher significance (>99.9%).
The mass distribution of close-in planets orbiting single stars and
wide binaries is not statistically different.
![]() |
Figure 3:
Distribution of the projected mass ![]() ![]() ![]() |
Open with DEXTER |
We checked the effects of the assumptions made on our analysis
on the resulting significance.
The significance of the difference of the mass distribution between
planets in tight binaries and single stars is larger than 99.9% for
AU, it is about 97.5% for
AU and
becomes smaller moving the limit beyond 100 AU.
The exclusion of super-planets more massive than
has some effect but the significance remains high (>99%) for both
tests. Changing the period limit of close-in planets
in the range 20-200 days keeps the significance always larger than 99%.
The RV semi-amplitude selection limit plays little role.
For some of the close-in planets with
orbiting stars classified as single there are indications of the presence
of significant dynamical perturbations.
HD 118203 b, HD 68988 b and HIP 14810 b have an
unusually high eccentricity
(e=0.31, 0.14, and 0.15 respectively) for their orbital
period (
6 days).
A second planet is orbiting HIP 14810 (B06), while in
the other two cases
linear trends in the residuals from the short period orbit indicate
the presence of a further companion in a long period orbit.
Their mass and period will be revealed by the continuation of the
observations, allowing a proper classification of their hosts as binary systems
or multi-planet hosts.
A candidate companion close to HD 162020 was reported by
Chauvin et al. (2006), whose physical association remains to
be established.
For HD 73256 b there is currently no indication
of additional companions.
Figure 4 shows the distribution of of planets with P> 40 days.
The lack of massive planets with P> 100 days
in binaries noted by E04 is no longer valid: three planets with
are now in the sample. All of them are orbiting the components of very
wide binaries (projected separation larger than 800 AU).
The probability of having 3 planets with P>100 days and
derived using the hypergeometric distribution is 17% when considering only
wide binaries and 7.4% when considering all the binaries together.
However, this probability becomes larger than 10% when considering
a different mass limit (
)
and KS and MWU tests do not
reveal statistically significant differences.
The mass distribution of long-period planets around tight binaries (7 planets) is not statistically significantly different with respect to that of planets in single stars and wide binaries (83 and 19 planets respectively).
![]() |
Figure 4:
Distribution of the projected mass ![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 5:
Distribution of log (period) for the planets
orbiting single stars (continous line), planets orbiting binaries with
![]() ![]() |
Open with DEXTER |
The period distribution of short period planets shows a pile-up for periods
around 3 days (see e.g. B06).
Only one of the five close-in planets in tight binaries has a period
close to 3 days (
Boo b).
Two have a very short period (HD 41004Bb and HD 189733b
have the shortest and the third shortest period
in the B06 catalog), while GL86b and HD 195019b have
longer periods (15 and 18 days respectively).
The differences are not statistically significant using KS and MWU tests.
When considering the period distribution of planets with P>40 days, there are hints that the period distribution of planets in tight binaries is different, with a lack of planets with periods longer than 1000 days. However, the KS and MWU tests do not show statistically significant differences.
When considering the whole period range (Fig. 5), the difference between the period distribution of planets around single stars and tight binaries (117 and 12 planets respectively) is more significant (94.3%) for the MWU test while it remains not significant according to the KS (50%). Similar values of significance are derived for the comparison between the period distribution of planets in tight vs. planets in wide binaries (25 planets). The discrepant level of significance between the two tests can be understood considering that the period distribution of planets in tight binaries appears to be shifted with respect to those orbiting single stars for the whole period range, but without very large deviations. The maximum difference of the cumulative distributions, the estimator used by the KS test, is then not significant, while the global shift is marginally significant by the MWU-test, that is sensitive instead to differences of the median of the distributions.
As the significance is not large and the different surveys have their own period sensitivity and biases concerning binarity (see Sect. 3.7), we conclude that the evidence of the different period distribution for planets orbiting single stars and components of tight binaries needs confirmation.
The period distribution of planets orbiting single stars and components of wide binaries are instead remarkably similar considering both separate period ranges and for the whole period range.
The possible link between high planet eccentricity and binarity has been considered since the discovery of the planet orbiting 16 Cyg B (Holman et al. 1997; Mazeh et al. 1997). It is thus relevant to consider the eccentricity distribution for planets orbiting tight and wide binaries and single stars.
Figure 6 shows the period-eccentricity diagram. The only two planets with e>0.80 (HD 80606b and HD 20782b) are in binary systems. However, planets with e up to 0.75 appear to be present also around stars currently known as single.
![]() |
Figure 6: Eccentricity vs. orbital period for planets in binaries (filled circles) and orbiting single stars (empty circles). Different sizes of filled circles refer to different periastron of the binary orbit (larger sizes: closer orbits). |
Open with DEXTER |
The eccentricity distribution of planets in tight binaries with
periods longer than 40 days (Fig. 7) is not significantly
different to those
orbiting single stars.
The comparison for planets in wide binaries indicates no significant differences according to the KS test (74%), while the MWU test gives a significance of 91.2% that the median of the two samples are different.
E04 noted that all the five planets in binaries in their sample with periods shorter than 40 days have an unusually low eccentricity (e<0.05). They estimated a chance probability of 3.8% using the hypergeometric distribution. While a tendency to a low eccentricity of short-period planets in binaries seems to be still present (4 planets of 5 and 4 of 6 with e<0.05, for the tight and wide binaries respectively, compared with 17 of 34 for planets in single stars) a similar analysis in our sample (taking together wide and tight binaries) yields a probability of 12%, therefore not confirming the E04 result at a highly significant level. KS and MWU tests also do not suggest statistically significant differences. Selecting only planets with periods between 10 and 40 days (to exclude objects whose orbits are circularizated by tidal effects) does not increase the significance of a lower eccentricity for short-period planets in binaries.
![]() |
Figure 7:
Cumulative distribution of the planet eccentricity for the
planets with period longer than 40 days orbiting
single stars (continous line), planets
orbiting binaries with
![]() ![]() |
Open with DEXTER |
The occurrence of systems with more than one planet
around the components of wide binaries
(4 of 23,
%) is similar with respect to that of planets
orbiting single stars
(15 of 112,
%).
On the other hand, no multiple planets have been yet discovered around the
components of tight binaries (13 objects).
The closest binary with a known multi-planet system is
And (projected separation 750 AU).
However, the small number of tight binaries in the sample
makes the lack of multi-planet systems
not highly significant (probability of 15% of occurring by chance using
the hypergeometric distribution).
The similar fraction of multiplanet systems around single stars and
wide companions is a further indication that the presence of a wide companion
does not alter too much the process of planet formation.
The study of the metallicity of planet host in binaries can shed light on possible different planet formation mechanisms (see Sect. 4). The close-binary planet-hosts appears slightly more metal poor than single stars, while the wide-binary planet-hosts are slightly more metal rich. The comparison of cumulative distributions (Fig. 8) using KS and MWU tests indicates probabilities of 12.1% and 4.6% respectively that the metallicity of close and wide binary planet hosts can be drawn from the same parent population. The differences with respect to single stars are less significant (probability of 23.3% and 17.8% for the close binary vs. single stars and and wide binary vs. single stars respectively using the MWU test, a smaller significance using KS). Considering together all the binaries makes any difference not significant.
![]() |
Figure 8:
Distribution of the stellar metallicity [Fe/H] for the planets
orbiting single stars (continous line), planets orbiting
binaries with
![]() ![]() |
Open with DEXTER |
While intriguing, the possible differences of the metallicity of binary planet hosts should be taken with caution. Beside the rather large probability of a chance result, it should be considered that the small number of planets in close binaries makes the role of individual objects critical (e.g. HD 114762, whose companion might be a brown dwarf or even a low mass star, Cochran et al. 1991). Furthermore, indirect selection effects might cause spurious correlations (see Sect. 3.7).
The B06 catalog includes planets discovered by several different surveys. Different RV surveys typically have their own biases concerning binarity and other properties (e.g. stellar mass and evolutionary phase, metallicity, etc.). Furthermore, they have different precision and time spans, and thus different planet mass and period sensitivities.
All these effects cannot controlled in our study, and we cannot exclude that they play a role in producing some of the marginal differences we showed (e.g. in metallicity and the period distributions). However, the most significant result, the different mass distribution of close-in planets in tight binaries, should not be affected by these selection effects. In fact, these planets are the easiest to discover for any survey because of their short periods and large radial velocity amplitude.
An additional caveat is represented by the uncertainty in the
upper limit of planetary mass, which is assumed in this study
to be
following B06, and by the other choices performed
in the analysis (e.g. period limits, etc.)
The sample of binaries among planet hosts is certainly not complete. A large fraction of the planet hosts were scrutinized only for comoving companions at rather large separations (a few arcsec). The analysis was performed considering the actual binary configuration that in some cases can be significantly different to that at the time of the formation of the planets. In the case of binary systems including a white dwarf, the probable original configuration can be derived estimating the original mass of the (now) white dwarf (see Appendix B), but other cases are speculative (see Appendix A.1).
Exploiting the continuously expanding samples of extrasolar planets, we investigated the occurrence of a difference between planets in binaries and planets orbiting single stars. The size of the sample allowed us to divide the binaries in two sub-groups, "tight'' and "wide'' binaries. The separation was fixed at a critical semiaxis for dynamical stability equal to 75 AU. This corresponds to separations of about 200-300 AU depending on the mass ratio between the components.
The most significant feature we were able to show is that the mass distribution of close-in planets around close and wide binaries is different. Only the close-in planets in tight binaries have a mass distribution significantly different to those orbiting single stars. This result implies that the formation and/or migration and/or dynamical evolution processes acting in the presence of a sufficiently close external perturber are modified with respect to single stars.
The fact that the hosts of massive planets in close orbit are typically binaries somewhat resembles the results by Tokovinin et al. (2006) i.e. short-period spectroscopic binaries have in most cases a further companion. That a third member could act as an angular momentum sink resulting in a closer inner binary was first suggested by Mazeh & Shaham (1979). Tokovinin et al. (2006) favor generation of high orbital eccentricity by Kozai resonances, and subsequent tidal circularization, as the mechanism to explain the close spectroscopic binaries and their link with the presence of additional companions.
While it cannot be excluded that some of the massive close-in planets in binaries were moved into their current location by Kozai migration, some other massive planets (GL 86b and HD 195019b) have periastrons of 0.108 and 0.137 AU, probably excluding the occurrence of tidal circularization. This calls for a more general process to explain the presence of massive close-in planets.
A possible effect on migration and mass accumulation
timescales in a disk in the presence
of a companion
was predicted by the models by Kley (2000).
He considers the evolution of a Jupiter mass planet embedded in
a circumstellar disk in binary
systems with a separation of 50-100 AU, with eccentricity
,
and stellar masses of 1.0 and 0.5
.
The presence of the companion significantly enhances the growth rate
and make the migration timescale of the planet shorter.
As discussed by E04, this can qualitatively explain the
presence of massive close-in planets in binary systems.
A step forward here is that we have shown
that most of the binary systems that host the "anomalous'' massive,
close-in planets are indeed not too different
to those modeled by Kley (2000), while E04 were
forced to consider also binary systems with very
wide separation, for which the described model is likely not to work.
The models by Kley (2000)
have as an initial condition an already formed planet and thus
do not include the formation phase of giant planets.
As discussed by Thebault et al. (2006),
a possible way to achieve the formation of planets in a swarm
of planetesimals perturbed by a stellar companion might be
type II runaway accretion identified by
Kortenkamp et al. (2001).
More detailed modeling is required to achieve a full understanding
of the processes that cause the differences of planet
properties in tight binaries.
Additional possible differences of the properties of planets in tight binaries (to be confirmed as the statistical significance is not very high and selection effects might be at work) are the different period distributions (there are no planets with P>1000, the close-in planets probably do not pile-up at period close to three days) and the lack of multi-planet systems.
In a few cases,
the lack of long period planets can be understood
as the result of dynamical stability constrains: for
Cep and HD 41004A, the critical semimajor axis
for dynamical stability is just 1.9 and
3 times the planet semimajor
axis.
For the other cases, a more general mechanism should be at work
(enhanced migration; Kley 2000 or lack of suitable conditions for
planet formation in the
external parts of the disks, under the gravitational
influence of the companion).
The fact that properties of planets in tight binaries are different to those orbiting single stars can be exploited to distinguish between different scenarios for their origin. As discussed by e.g. Hatzes & Wurchterl (2005) and Thebault et al. (2004), the presence of planets in binaries with a separation of 20 AU or less represents a challenge for our current view of planet formation. Such problems can be overcome if the binary was initially wider, and its orbit was modified through dynamical interactions in star clusters after planet formation (Pfahl & Muterspaugh 2006). The different mass distribution of planets in close binaries implies that, if this is the case, the dynamical interactions should not only modify the binary orbit but also force a massive planet in external orbit to migrate inward, as massive hot Jupiters are rare around single stars. This point should be addressed by future dynamical simulations. Kozai oscillations and tidal circularization after the dynamical interaction cannot be excluded in some cases (large relative inclinations between the planet and binary orbits might be the result of the proposed dynamical encounters). However, less massive planets originally in external orbits should follow a similar fate and thus we expect that low mass planets should also be present in close orbits in tight binaries, at odds with current observational data. Furthermore, as discussed above, Kozai migration cannot explain planets with periods longer than a few days and low eccentricities.
Further clues to the origin of planets in tight binaries
can be derived by a determination of their frequency.
Pfahl & Muterspaugh (2006) predict that the fraction
of close binaries that dynamically acquire giant planets
is about 0.1%, with an uncertainty of about one order
of magnitude.
They note that the rough number of these planets seems
too high to be compatible with these predictions.
However, the number of tight binaries
included in the radial velocity samples is not well known.
A more detailed analysis is postponed to a forthcoming study.
The mentioned difficulties of the binary-interaction scenario to account for the occurrence of planets in close binaries leave open the alternative hypothesis, i.e. giant planets formed in binaries with small separation at the time of planet formation, possibly in a different way than planets around single stars.
The metallicity distribution of tight binaries with planets can give clues to the formation mechanism for these planets. The large effect of the stellar metallicity on the frequency of planetary companions is well established (Fischer & Valenti 2005; Santos et al. 2004) and can be understood in the framework of the core-accretion model for the formation of giant planets (Ida & Lin 2004; Robinson et al. 2006). If these planets formed in a different way than those orbiting single stars and the components of wide binaries, e.g. by disk instability triggered by dynamical perturbations (Boss 2006), we should expect a lower metallicity for planet hosts in close binaries, as the disk instability mechanism is expected to be rather insensitive to metallicity (Boss 2002). There is a marginal indication of a larger fraction of metal poor stars among tight binaries with planets. However, the distribution of metallicity remains on average larger than that if the solar vicinity, with a mean and median larger than the solar value. This indicates that high metallicity is a factor favoring planet formation also in tight binary systems.
The mass and period distributions of planets in wide binaries are not statistically significant different to those of planets orbiting single stars. The lack of long-period massive planets in binaries shown by E04 is not confirmed at a highly significant level in our study. The fraction of multi-planet systems is similar for planets orbiting single stars and components of multiple systems. These results indicate that a distant companion (separation >300-500 AU) does not significantly affect the process of planet formation. The only marginally significant difference between planets orbiting single stars and components of wide binaries concerns the planet eccentricity.
The eccentricity distribution of planets in wide binaries shows
a marginal excess of high-eccentricity ()
planets.
However, current results indicate that the high planet eccentricity is not
confined to planets in binaries,
and that the possible differences in eccentricity are limited
to the range
.
This indicates that there are mechanism(s) generating planet
eccentricity up to 0.4-0.5 that are independent of the binarity
of the planet host, and are characteristic of formation and
evolution of a planetary system
(e.g. disk-planet interactions: Tremaine & Zakamska 2004;
planet-planet scattering:
Marzari & Weidenschilling 2002;
Ford et al. 2003).
These probably act
during or shortly after planet formation.
Further eccentricity enhancements, possibly linked to the presence
of a companion, might take place
at later epochs.
In fact, Takeda et al. (2006) noted that most
high-eccentricity planets orbit old stars (ages >5 Gyr).
Mechanisms that require long time scales to modify
planetary orbits then seem favored.
A popular mechanism that might produce high eccentricity planets in binaries
is represented by Kozai oscillation.
Kozai oscillation is expected to be at work also for
low mass companions at large separation (the case of
most of the companions of planet hosts discovered up to now),
provided that a sufficiently long time is available
and that the planet has a sufficiently long period
to prevent the suppression of Kozai oscillation
by general relativistic effects
(Takeda & Rasio 2005; Wu & Murray 2003).
As an example, a brown dwarf companion such as HD 3651B
(mass about
,
projected separation about 500 AU)
might be able to induce the eccentricity oscillation for
a planet with an initial semimajor axis larger than about 3 AU.
Kozai oscillations are expected to produce large planet eccentricities. Wu & Murray (2003) considered in detail the possibility that the high eccentricity of HD 80606b (e=0.93) is due to Kozai oscillation. This was found to be possible only for a relative inclination of the planetary and binary orbit larger than 85 deg. In the case of the only other planet with e>0.8 (HD 20782b), the period of the Kozai eccentricity oscillation derived from Eq. (2) of Takeda & Rasio (2005) is much longer than the Hubble time. Therefore, the extremely high eccentricities of HD 20782b is unlikely to be due to the Kozai mechanism (see Appendix A.1), unless further companions exist in the system at smaller separation.
An alternative to Kozai oscillations is represented by chaotic evolution of planetary orbits induced by dynamical perturbations (Benest & Gonczi 2003; Marzari et al. 2005). This kind of evolution might also arise for the dynamical interactions within the planetary system, without requiring a stellar companion in a wide orbit.
To further investigate the origin of the high eccentricities, we consider the eccentricity distributions of planets residing in multiplanet systems and isolated planets (Fig. 9). There is a lack of high eccentricity planets in multi-planet systems (as noted also by Takeda et al. 2006). When considering planets with P>40 days, the comparison of the two distributions reveals a difference that is significant to about a 89.1% confidence level (according to the MWU test), similar to the results of the planets in wide binaries vs those in single stars. This can be understood considering that planetary systems with highly eccentric planets are in general more subject to close encounters, and then less stable, than a system populated with planets in low eccentricity orbits. If the chaotic evolution of a planetary system has as a final result a planet on a highly eccentricity orbit, it is probable that that planet remains alone. On the other hand, Kozai oscillations are more likely to take place in single-planet systems, as they are suppressed in the presence of significant mutual planet-planet interactions. Therefore, eccentricity enhancements due to Kozai oscillation are expected to occur for single-planets systems with a companion in a suitable binary configuration. A proper evaluation of the link between high eccentricities, binarity of planet hosts, and planet multiplicity requires larger samples to allow more statistically significant inferences.
![]() |
Figure 9: Distribution of the planet eccentricity for the single planets (continous line), and planets in planetary systems (dotted line). Upper panel: histogram; lower panel: cumulative distribution. |
Open with DEXTER |
The only highly significant difference shown in our analysis concerns the mass distribution of short-period planets in tight binaries. Massive planets in short period orbit are found in most cases around the components of rather tight binaries. Other possible peculiar features of planets in tight binaries with respect to planets orbiting single stars, such as a lack of long-period planets and multiple planets, are of marginal significance and need confirmation.
The properties of exoplanets orbiting the components of wide binaries are compatible with those of planets orbiting single stars, except for a possible greater abundance of planets on highly eccentric orbits. The previously suggested lack of massive planets with P>100 days in binaries is not confirmed.
This result indicates that the binary separation plays a role in affecting the properties of planets and that the simple classification of planets in binaries vs planets orbiting single stars is not adequate to describe the dynamical effects. This opens the perspective for more detailed modeling of the role of stellar companions on the formation and evolution of planetary systems.
Acknowledgements
This research has made use of the SIMBAD database, operated at CDS, Strasbourg, France, and of data products from the Two Micron All Sky Survey. We thank J.R. Hurley for providing the BSE program. We thank G. Takeda for useful information about the stellar ages of planet hosts. We thank R. Gratton for his comments and suggestions. We thank the anonymous referee for a prompt and detailed report. This work was funded by COFIN 2004 "From stars to planets: accretion, disk evolution and planet formation'' by the Ministero Univ. e Ricerca Scientifica Italy.
Table 1:
Properties of planets in multiple systems, of their host
stars and their companions: planet ,
period, semimajor axis,
eccentricity, radial velocity semi-amplitude,
metallicity and mass of the planet host (from B06); projected separation,
semimajor axis
and eccentricity (when available) and mass of the companions,
and critical semiaxis for
dynamical stability of planets. The asterisks in the last column mark systems
discussed individually in Appendix C.
As planet announcements arise quite frequently,
existing compilations quickly become incomplete.
Twenty-one stars are not considered in the latest
study of multiplicity of planet hosts
(R06). Most of them are
new planet hosts while some are stars with companions
with mass between 13 to 24 ,
not considered by R06
but included in B06 catalog.
We searched for companions of these stars in existing astronomical catalogs such as Hipparcos (ESA 1997), CCDM (Dommaget & Nys 2002), WDS (Mason et al. 2001), and ADS (Aitken 1932). For the following 16 stars there are no evidence of binarity in the literature: HD 4308, HIP 14810, HD 33283, HD 66428, HD 81040, HD 86081, HD 99109, HD 102195, HD 107148, HD 109143, HD 118203, GJ 581, HD 137510, HD 187085, HD 212301, HD 224693.
HD 33564 and HD 164922 are included in CCDM and WDS (the latter also in ADS) as triple systems, but the companions display different proper motion to the planet hosts and thus they are not physical companions. Two additional objects (HD 20782 and HD 109749) are binaries. The physical association of the companions and their physical properties are discussed in the next subsections. Finally, the binarity of HD 189733 is discussed in a dedicated paper by Bakos et al. (2006a).
HD 20782 is listed as a binary in the CCDM catalog (CCDM 03201-2850). The companion is HD 20781. The very large separation, 252 arcsec, that corresponds to 9080 AU at the distance of HD 20782 (36 pc) makes a detailed check of physical association mandatory. Table A.1 lists the stellar parameters of the two stars. Hipparcos parallaxes and proper motions, and Nordstrom et al. (2004) radial velocities, are fully compatible within errors for the two components. This strongly suggests physical association. Therefore, HD 20782 and HD 20781 form a very wide common proper motion pair. Considering the nominal radial velocity and proper motion differences, the pair could be bound. The observed separation is very large but not extreme for binaries in the solar neighborhood (Poveda & Allen 2004). Clearly, the orbit of the wide binary remains unconstrained.
Table A.1: Stellar properties of HD 20782 and HD 20781.
We tested the dynamical stability of the system via numerical integration and
we have performed 10 simulation of this system following the evolution
for 200 Myr (577 binary orbits) using a RADAU integrator
(Everhart et al. 1985).
The binary orbit was assumed with typical values of
AU and
,
the orbital parameter of the planet was chosen equal to the observed values,
and inclination was chosen slightly inclined (
)
with respect
to the binary orbital plane.
The planet results stable over all
the simulation,
thus we argue that the large binary separation ensures dynamical
stability of the planet
in spite of its very high eccentricity.
Only two extrasolar planets have orbits with eccentricities larger than 0.8: HD 20782b and HD 80606b. Both their host stars are member of wide common proper motion pairs.
While the link between the high eccentricity and binarity still needs statistical confirmation, nevertheless it is interesting to investigate the possible ways in which a distant companion might have induced extreme planet eccentricities.
As discussed in Sect. 4, the extremely high eccentricity of HD 20782b is unlikely to be due to the Kozai mechanism, as the period of the eccentricity oscillation is much larger than the Hubble time.
An interesting possibility to explain both the very large separation
of the binary and the very high eccentricity of the planet orbiting
the primary is a dynamical encounter of the binary
(or of an originally higher multiplicity system)
within a star cluster or with a passing star,
that might have perturbed both the binary and the planet orbit.
According to the simulations by Weinberg et al. (1987),
the probability of survival for a binary with an initial semimajor axis
of about 10 000 AU after 4 Gyr is less than 50%.
Furthermore, the actual semimajor axis
is probably larger than the observed projected separation by about 30% (Duquennoy & Mayor 1991; Fischer & Marcy 1992),
further decreasing the probability of survival, if this is the case
of HD 20782.
The galactic orbit of HD 20782 (Table A.1)
makes the star on average
closer to the Galactic center than the Sun,
further increasing the chance of stellar encounters, because
of the higher stellar density.
Therefore, at the age of the system (
7.1+1.9-4.1 Gyr;
Valenti & Fischer 2005;
;
Takeda et al. 2006),
the probability that the currently observed orbit is not the original one
but it was significantly modified by dynamical encounters in the Galactic
disk is rather high.
The process of destruction of wide binaries is typically the result of a large number of weak encounters (see e.g. Portegies Zwart et al. 1997). The changes in the binary orbit caused by one or more stellar encounters might alter the orbit of the planet(s) orbiting one of the components. A multi-planet system can be destabilized, leading to chaotic evolution of the planetary orbits. Dedicated dynamical modelling should be performed to check if the current configuration of the system might represent a possible outcome of a stellar encounter. This is postponed to future work. The hypothesis of stellar encounters appears less appealing to explain the extremely high eccentricity of HD 80606b, as in this case the projected separation of the binary is about 1200 AU, 7-8 times smaller than for HD 20782. However, the real orbit remains unknown.
The possibility that the high planet eccentricity is unrelated to the presence of the stellar companions and is caused by internal evolution of the planetary system cannot be excluded. Conclusive inferences on the origin of the very high eccentricity of HD 20782b and HD 86606b are not possible from the available data. Further clues will come from follow-up observations, such as the search for further planetary and stellar companions, and from the discovery of further planets with very highly eccentric orbit and the study of the binarity of their hosts.
HD 109749 is listed in Hipparcos and CCDM catalogs as a binary (CCDM J12373-4049). The Hipparcos magnitude of the secondary transformed to the standard system following Bessell (2000) yields V=10.76, about 2.5 mag fainter than the primary (V=8.09, Fischer et al. 2006; see Table A.2). Table A.3 shows the relative astrometry for the two components. The constancy of the projected separation and position angle, coupled with the rather large proper motion of HD 109749, makes the physical association between the two objects very probable. The projected separation of 8.35 arcsec corresponds to about 500 AU.
Table A.2: Stellar properties of HD 109749 and HD 109749B.
Table A.3: Relative astrometry of HD 109749.
![]() |
Figure A.1:
Absolute magnitudes and colors for the components of HD 109749
(filled circles). Left panel: MK vs. V-K;
right panel: ![]() |
To further check the physical association between the two stars, we identified
them in 2MASS catalog (Cutri et al. 2003).
The photometry is shown in Table A.2.
When placed in a color-magnitude diagram, both components lies along the same isochrone,
indicating a common distance and confirming the physical
association (Fig. A.1).
The masses derived from the isochrone shown in Fig. A.1
are 1.11 and 0.78 .
Three planets are orbiting stars with white dwarfs
companions.
For binaries wide enough to escape common envelope
evolution, the mass loss from the originally
more massive star during the RGB and AGB phases determines
a widening of the orbit.
The mass loss phase lasts several Myr, longer than the typical
binary orbital
periods.
In these conditions, the widening of the binary orbit
can be approximated by
,
and no eccentricity
changes are expected (Lagrange et al. 2006
and references therein).
In this section, we consider the three planet-host binaries
with a WD, placing constraints on the original
configuration of the system in terms of binary separation
and mass ratio.
The companion of Gl86 was discovered by
Els et al. (2001) and it was
classified as a brown dwarf on the basis of the near-infrared colors.
Mugrauer & Neuhauser (2005) showed that the companion is instead
a white dwarf.
Lagrange et al. (2006) further constrained the properties
of the system, using archive radial velocity and new photometry and
position measurements.
However, they also noted some inconsistency between the white dwarf
cooling age (1.2-1.8 Gyr depending on white dwarf mass and atmosphere
composition), the age of the star as derived from
chromospheric activity (about 2-3 Gyr), and the presence of the planet.
In fact, such a young age implies an original mass
of the companion that is too large to be compatible with
a dynamically stable orbit for the planet once the widening
of the binary orbit
due to mass loss from the originally more massive star
is taken into account.
A picture that might explain the properties
of the system is the following.
Kinematic parameters (U=-102 km s-1, V=-75 km s-1, W=-29 km s-1,
galactic orbit
kpc,
kpc, e=0.39,
kpc;
Nordstrom et al. 2004) indicate a very
old age
.
The kinematic criterium of Bensby et al. (2005) suggests
a probable membership to the thick disk.
Chemical abundances can also be used to constrain the galactic population.
At the metallicity of GL86 (
),
differences between thin and thick stars are fairly small
(the typical [
/Fe] of thin disk and thick stars is
about +0.05 and +0.15 respectively, Bensby et al. 2005).
High-resolution abundances analysis of several elements were published by
Allende Prieto et al. (2004), Valenti & Fischer (2005),
Santos et al. (2004), and
Gilli et al. (2006).
The uncertainty in the chemical abundances of GL 86
and the discrepancies between different studies do not allow
a firm classification, but they appear compatible
with thick disk membership or with thick disk-thin disk transition
objects
.
Thick disk stars with this metallicity have typical ages of
8 Gyr according to Bensby et al. (2005).
This is compatible with the lower limit (68% confidence level) of 8.5 Gyr
estimated by Takeda et al. (2006) from isochrone fitting.
If GL 86 is instead a thin disk star,
its kinematic parameters strongly
supports an old age.
Matching the cooling age of the white dwarf given by
Lagrange et al. (2006)
with the assumed age of about 8 Gyr for the system
requires an original mass of about
.
Older ages yield lower initial masses.
A
star ends its life as a
white dwarf. The amount of mass loss implies a widening of the binary
orbit of about 40-50% of its original separation.
To match the semimajor axis of 18.4 AU derived by
Lagrange et al. (2006),
the semimajor axis of the original orbit had to be about 13 AU.
A fraction of the mass loss by GL86 B during the
RGB and AGB evolution should have been
captured by the planet host.
A rough estimate based on the BSE binary evolution code by
Hurley et al. (2002), using the default
wind and accretion parameters, is
.
The
accreted mass corresponds to roughly half of
the mass of the convective envelope of GL86
(
;
Takeda et al. 2006).
The lack of peculiar abundance of neutron capture
elements ([Ba/Fe] = -0.12; Allende Prieto et al. 2004)
suggests an original mass for the AGB star
smaller than
(Busso et al. 1999),
in agreement with our estimate of
.
Along with mass, angular momentum should also have been transferred to the planet host. This was invoked to explain the high rotation rates of companions of hot white dwarfs such as the barium dwarf 2RE J0357+283 (Jeffries & Smalley 1996). The accretion of angular momentum by GL86 would explain its activity level typical of a much younger star.
This scenario would be able to explain the available observational constraints, including the stability of the planetary orbit. However, the formation of a planet with a companion at about 13 AU (and with an eccentricity of 0.4) represents a challenge for current models (Hatzes & Wuchterl 2005). The signature of accretion of angular momentum by GL86 is an indication that the binary system was in its current status at the time of the mass loss by GL86B, about 1.5 Gyr ago. Earlier modifications of the binary orbit due to e.g. dynamical interactions within the native star cluster (Pfahl & Muterspaugh 2006) cannot be excluded. The possible effects of the mass and angular momentum accretion on the mass and orbit of the planet remain to be investigated.
HD 27442 =
Ret has a companion at a
projected separation of about 13 arcsec = 240 AU.
According to Chauvin et al. (2006),
visual and near-infrared photometry are not compatible
with any main sequence star. They then argue that the
companion is a white dwarf.
The mass of the planet host yields a lower limit to
the original mass of the companion.
Valenti & Fischer (2005) determine a mass of
(the value also adopted by B06), with
an age of 3.5 Gyr, while Takeda et al. (2006)
list
with an age of
2.84+0.60-0.36 Gyr.
A lower limit to the WD mass
estimated using the BSE code (taking into account
the high metallicity of
Ret;
[Fe/H] = +0.42)
is
.
Theoretical models should be taken with caution at these
extremely high metallicities, as they fail to explain the
WD luminosity function of the old super-metal-rich open
cluster NGC 6791 (Bedin et al. 2005).
The relatively wide separation
from the primary (13 arcsec) and the moderately
bright magnitude (
,
WDS) open the perspective for
a more detailed characterization of
Ret B.
We can guess that the system
had lost more than one third of its original mass,
implying a substantial widening of the binary orbit.
The original separation might have been about 150 AU, making the
influence of the companion not negligible.
HD 147413 was classified as a member of the Ursa Major
moving group by Soderblom & Mayor (1983) while King
et al. (2003) classified the membership as unprobable,
on the basis of the discrepant kinematic parameters. The
activity level and the position on the color-magnitude diagram
are instead fully compatible with membership.
Its companion (WD 1620-39 = HIP 80300) is a widely studied
WD, for which Silvestri et al. (2001)
derived a mass of
and a cooling age of 20 Myr.
The planet host was claimed to be a barium star ([Ba/Fe] = +0.37)
by Porto De Mello & Da Silva (1997) but
Castro et al. (1999) showed that the high barium
abundance is shared by other stars in the UMa moving group
(and possibly by other young stars).
This points
toward a primordial abundance and not to accretion
by the WD companion, unexpected on the basis of the large
separation (more than 4000 AU).
We estimated the possible original mass of the (now) WD companion
using the BSE code (Hurley et al. 2002).
Adopting the most recent age estimate of the UMa moving group
(500 Myr; King et al. 2003)
and the cooling age by Silvestri et al. (2001) leads to
an original mass of
and a WD mass of
,
some
larger than that derived by Silvestri et al. (2001).
The WD mass and cooling age by Silvestri et al. (2001)
can be matched simultaneously by adopting an initial mass of
and a system age of 1 Gyr. This age is still marginally
compatible with the
photometry and activity level, if the star is not
a member of UMa moving group.
We conclude that the original separation of the binary was about a half
of the present one, as the system has lost about half of its original mass.
Considering the current very wide separation, this does not change the
classification of this planet host as a very wide binary.
We further discuss individual objects.
The criteria for the inclusion of planets in the
"Catalog of Nearby Exoplanets'' are different to
other compilations such as the Extrasolar Planet Encyclopedia
maintained by J. Schneider.
Three stars classified as planet hosts by the Extrasolar Planet
Encyclopedia and not included in the Catalog of Nearby Exoplanets
are members of multiple systems:
A few additional planet hosts not included in Table 1 have companion candidates whose identification or physical association needs confirmation: These are: