Issue |
A&A
Volume 535, November 2011
|
|
---|---|---|
Article Number | A50 | |
Number of page(s) | 6 | |
Section | Planets and planetary systems | |
DOI | https://doi.org/10.1051/0004-6361/201116907 | |
Published online | 03 November 2011 |
A binary merger origin for inflated hot Jupiter planets
1
Centro de Astrobiologia CSIC – INTA, Carretera Torrejón-Ajalvir km
4,
28850
Madrid,
Spain
2
Max-Planck-Institut für Astrophysik, 85748
Garching,
Germany
e-mail: henk@mpa-garching.mpg.de
3
Instituto de Astrofisica de Canarias,
calle Via Lactea, 38200 La
Laguna, Tenerife,
Spain
Received: 16 March 2011
Accepted: 27 August 2011
We hypothesize that hot Jupiters with inflated sizes represent a separate planet formation channel, the merging of two low-mass stars. We show that the abundance and properties of W UMa stars and low mass detached binaries are consistent with their being possible progenitors. The degree of inflation of the transiting hot Jupiters correlates with their expected spiral-in life time by tidal dissipation, and this could indicate youth if the stellar dissipation parameter Q∗' is sufficiently low. Several Jupiter-mass planets can form in the massive compact disk formed in a merger event. Gravitational scattering between them can explain the high incidence of excentric, inclined, and retrograde orbits. If the population of inflated planets is indeed formed by a merger process, their frequency should be much higher around blue stragglers than around T Tauri stars.
Key words: planets and satellites: formation / binaries: eclipsing / planets and satellites: dynamical evolution and stability / blue stragglers / methods: statistical
© ESO, 2011
1. Introduction
Ever since the first discoveries of Jupiter-mass planets around solar-type stars (Mayor & Queloz 1995; Butler & Marcy 1996), the origins of hot Jupiters have inspired a variety of ideas. Formation scenarios invoked in the literature include orbital migration, radiative stripping, planet scattering, and secular chaos (Ford & Rasio 2008; Ida & Lin 2010; Naoz et al. 2011). Each of these mechanisms can explain some of observational properties of hot Jupiters, but not all of them. The main properties are summarized as follows: a) a pile up of orbital periods around three days, b) a wide range of obliquity, c) a lack of close companion stars, d) a mass function different than for other types of exoplanets, e) a limited mass range of host star masses, f) a steep rise in frequency for super-solar metallicities, g) a frequency insensitive to stellar metallicity from solar down to very metal poor stars, and h) a wide range of planetary sizes for a given planetary mass. Together, all these observational patterns constitute a substantial set of constrains that no single theory has been able to explain, though planet population synthesis (Mordasini et al. 2009) has met with success in accounting for several of the observed features. It is possible that different pathways exist to the formation of hot Jupiters. Hopefully they would leave different signatures on the planet properties that could be used to identify their origin. In this paper we hypothesize a new scenario, where they form in excretion disks produced by low-mass binary mergers, which, we argue, could be the cause of the large sizes of some transiting hot Jupiters.
The radius of brown dwarfs and giant planets is determined by the interior temperatures which, for single objects, depend on the total mass, age, and the opacity. At a given planetary mass and composition, all planets are expected to converge to a constant radius (R0) with time, which is determined by electron degeneracy pressure, Coulomb pressure and exchange pressure, in the interior (Zapolsky & Salpeter 1969; Stevenson 1991). Since the stellar hosts of exoplanets are usually on the main sequence, it is difficult to estimate their ages, but typically they are thought to be several Gyr old, and if their planets have the same age, they should have had enough time to evolve to their degeneracy dominated configuration where R → R0. The advent of wide area transit surveys from the ground (such as WASP) and in space (CoRoT and Kepler) has brought about the discovery of over a hundred transiting exoplanets. The combination of detailed photometric transit observations and high-precision radial velocity measurements have provided mass and radius determinations for those planets. With all this wealth of data, it has become clear that our understanding of the mass-radius relationship in exoplanets is far from satisfactory. In particular two properties of transiting exoplanets are especially troubling: (1) the wide spread of planetary radii for masses around 1 Jupiter mass (from 10% to as much a factor 2), and (2) the large size of some exoplanets, which has been termed the “radius anomaly” of hot Jupiters (Bodenheimer et al. 2001; Guillot et al. 2006).
Several mechanisms have been invoked to account for the radius anomaly, such as diffusion of stellar irradiation into the planet interior, enhanced opacities, inefficient convection, tidal heating (Leconte et al. 2009), Ohmic heating (Batygin & Stevenson 2010), burial of heat by turbulence (Youdin & Mitchell 2010), and strong winds (Chabrier et al. 2010). The wide spread in radius excess (more than a factor of 2 for planets with masses around that of Jupiter) makes it unlikely that the cause of the radius anomaly lies in systematic errors in the physics of interior structure (EOS, opacities) or other factors that are just functions of intrinsic properties of the planet.
Recent Spitzer observations of the transiting planet WASP-12b have indicated that the atmospheric composition may be highly non-solar (Madhusudhan et al. 2011). Then it is possible that the scatter in the mass-radius relationship is at least partly due to variations in the atmospheric composition of hot Jupiters and in their cloud thickness (Burrows et al. 2011). However, to produce hot Jupiters as inflated as observed with enhanced atmospheric opacities a factor 10 above the solar mixture is needed for cooling and contraction to be slowed down significantly (Guillot 2008). The conventional planet formation models (core accretion) have indeed predicted that Jupiter-mass planets could have different chemical compositions from those of the nebulosity from which their parent star formed (Pollack & Bodenheimer 1989). However, in these models, the higher mass planets that form in the protoplanetary disk have abundances closer to that of the central star. Thus, it is unclear how planets larger than Jupiter and even brown dwarfs may form with chemical abundances radically different from those of the host stars. Furthermore, the degree of inflation of hot Jupiters seems to be quite insensitive to the metallicity of the host star (Laughlin et al. 2011).
The transiting brown dwarf CoRoT-15b has shown that the radius anomaly extends into the brown dwarf mass domain, making it more unlikely that the correct explanation is stellar irradiation (Bouchy et al. 2011). The age of the stellar host has been estimated to lie in the range 1.14−3.35 Ga using evolutionary models, while the cooling age of the brown dwarf is around 0.5 Ga, although high-metallicity models could recover the coevality of the system (Burrows et al. 2011). This example illustrates that the radius anomaly could be interpreted either as an age discrepancy between the stars and their close substellar-mass companion or as a difference in chemical composition.
2. Hot Jupiter formation in binary mergers
A close binary with a system mass in the range of planet hosts consists of stars with convective envelopes. The magnetic activity of these stars causes them to lose angular momentum by “magnetic braking”. Tidal coupling between the stars transmits (part of) this loss to the orbit, the binary orbit becomes narrower, and the angular momentum loss speeds up. Eventually, the stars go through a merger process. If in this process a substantial amount of mass remains in orbit around the primary, it would form a disk in which planets could form. Such a scenario for “second generation” planets around single main sequence stars has already been discussed by Tutukov et al. (1991, 2004), but not in the context of hot Jupiters. The low angular momentum of the post-merger disk compared with a normal protoplanetary disk means that the orbits of these planets would naturally be close. Tidal interaction with the star would eventually cause the innermost planets to be accreted. Other examples of “late planet formation” scenarios in circumbinary excretion disks are those proposed to explain the existence of planets around pulsars and highly evolved stars (Banit et al. 1993; Wolszczan & Kuchner 2010; Perets 2011).
Direct observational evidence that contact binaries merge to become single stars has been elusive, but one case is now known. OGLE photometric monitoring has shown that the progenitor of Nova Sco 2008 (V1309 Sco) was not the expected cataclysmic variable but a contact binary with an orbital period of 1.4 d (Tylenda et al. 2010). Over the years leading up to the outburst, the period was observed to decrease ever more rapidly with time, as expected in a merger event. The nova-like behavior of the event shows that energetic mass ejection can accompany a W UMa merger (at least, at the high end of the W UMa mass distribution represented by V1309 Sco). An example of circumstellar dust formation around a stellar merger event could be the mid-infrared flux variability observed in V838 Mon (Wisniewski et al. 2008), although the nature of this object prior to its 2002 outburst is still a matter of debate.
If a hot Jupiter planet is formed in the excretion disk produced by a binary merger, the cooling age of the planet would reflect the time elapsed since the merging event. Calculations of the life times of a contact binaries before their final merger indicate that they can be rather long (6−10 Gyr, Stępień 2011). Thus, hot Jupiters produced in mergers would look much younger than their host stars. Their host population is a few Gy old, as judged from its kinematics and their inferred evolutionary status in individual cases (such as Lee et al. 2011).
An obvious indicator of a young cooling age for any planet is that it has a radius larger than an older planet with the same mass. Recognizably inflated hot Jupiters (e.g., 20% above the nominal mass-radius relation) comprise about 1/one-third of the known sample, and would need to be about 108 yrs or younger to explain their sizes with standard evolutionary models (Baraffe et al. 2003). To examine the likelihood of this scenario we proceed by checking whether there is a suitable population of progenitor binaries, and by looking for a connection between planet inflation and tidal evolution.
2.1. Progenitor population
If the progenitors are recognizable as such for a fraction f of the age of the host stars (3 Gy, say), the progenitors must thus be some 30f times more abundant than the inflated planet population and should be easily identifiable. The frequency of transiting planets found by CoRoT around G-type main sequence stars is around 1 per 2000 stars, so the recognizably inflated fraction of 1/3 has an abundance of about 1.6 × 10-4 per G star. Their progenitors must then have an abundance of ~5 × 10-3f. An obvious candidate population is contact binaries (W UMa stars).
The abundance of stars classified as W UMa on the basis of their (EW type) light curves is reasonably well established. From a complete sample of bright stars (Rucinski 2002) and ASAS light curves (Rucinski 2006), an abundance of 2 × 10-3 per main sequence star of the same spectral type has been derived. Since they are eclipsing binaries, the selection effects of orbital inclination are probably similar to those of the inflated transiting planets. In the presently most plausible interpretation of the evolution of W UMas (Stępień 2006, 2009), they are long-lived (a few Gy), as for the estimated age of the planet host population, i.e. f ≈ 1. On the basis of these numbers, contact binaries thus appear to be a possible progenitor population, although they may not be the only relevant possibility as discussed below.
The outcome of close binary evolution (e.g., Ritter 1996) is different depending on the mass ratio, total mass, and the evolutionary status of the components. The primary star can be brought (or kept) in synchronous rotation with the orbit if the mass ratio q = M2/M1 of the secondary is large enough, and the stars close enough for tidal interaction to be effective. Ongoing angular momentum loss from the system by magnetic braking narrows the orbit until one of the stars fills its Roche lobe and mass transfer starts. If, on the other hand, the orbital angular momentum of the secondary is too small to spin the primary up to corotation with the orbit, the orbit will shrink by tidal interaction (even if magnetic braking were ineffective) until one of the stars fills its Roche lobe. If the mass ratio is below Darwin’s stability limit (“tidal instability”, Darwin 1879; Hut 1980), this will happen even if the stars initially rotate in synchrony with the orbit. For low mass main sequence stars, this limit is on the order of q = 0.08 (Rasio 1995).
After the binary orbit has shrunk until one of the stars fills its Roche lobe, the further evolution depends on the nature of the lobe-filling star, since its internal structure determines how its size responds to loss of mass from its surface. For main sequence stars the star that first fills its lobe is the more massive one. The simplest case is when it is a fully convective, low-mass star. Its adiabatic mass-radius exponent is then negative; i.e., sudden mass loss causes it to expand. Mass loss then causes it to overfill its Roche lobe, and the mass loss rate increases exponentially, with the final stages happening on a dynamical (orbital) time scale (cf. Ritter 1988). Since the receiving star does not fill its Roche lobe yet and, being of lower mass, is also fully convective with a negative mass-radius exponent, it can receive a large amount of mass from the primary before also filling its Roche lobe. The SPH simulations by Rasio & Shapiro (1995) and grid based simulations (D’Souza et al. 2006) show that the final merger of such binaries happens on a time scale of some ten orbits.
The situation is more complicated when the mass-radius exponent of the primary is positive or becomes positive after an initial phase of dynamical mass transfer (Stępień 2006). A longer phase of mass transfer on the thermal time scale then takes place, until the mass ratio of the system has reversed, as in Algol type binaries. For primary masses on the order of 1 M⊙, the star is likely to have evolved a (small) helium core during the magnetic braking period that brought the binary into contact. Stępień (2006) shows that this allows binaries in the mass range of W UMa stars to settle into a stable contact configuration, which lasts for a few Gy until angular momentum loss by magnetic braking causes the stars to finally merge. This solves a longstanding puzzle in the theory of contact binaries and explains their high abundance.
The nature and abundance of W UMas makes them a plausible candidate population for making hot Jupiters, but they are probably not the only ones. The distribution of (total) mass of W UMas in the catalog of Gazeas & Stępień (2008) peaks at 1.5−2 M⊙, the distribution of known host masses of transit planets around 1.1 M⊙. In Fig. 1 we compare the distribution of total masses of W UMa stars with that of transiting hot Jupiter host stars. The overlap between the distributions is in fact only modest. The V1309 Sco event shows that some mass can be lost in the final merger, although loss as large as several tenths of a solar mass seems unlikely. It is thus worth exploring other channels. We propose here that the additional channel is in fact the direct merger of low-mass binaries (M1 + M2 ≲ 1.2 M⊙) on a dynamical time scale as discussed above. There is no obvious reason why such systems, at masses of 0.8−1.2 M⊙, would be formed at rates much less than the 1−2.5 M⊙ binaries that end up becoming W UMas. Their lower luminosity, and the absence of the extended contact phase that makes W UMas stand out, would make them a much less prominent population. Detached binary systems that could be such progenitors are known (e.g., Coughlin et al. 2011).
Fig. 1 Masses of the host stars of 103 transiting planets for which mass and radius of star and planet are given in the “Extrasolar Planets Encyclopaedia” compared with the masses of W UMa binaries from Gazeas & Stępień (2008). |
2.2. Spiral-in time scales
Inflated hot Jupiters orbit at distances close to their host star, typically within 0.05 AU (Fig. 2). At such distances, friction in the tides they raise on the star cause their orbits to circularize. The functional dependence of this process on system parameters is well established (Goldreich 1963; Jackson et al. 2009). The associated time scale is proportional to a stellar tidal dissipation parameter Q∗' whose value is less well known since it depends on details of the interaction of the tides with the flows in the convective envelope of the host star. By fitting the observed distribution of excentricities of planet orbits, Jackson et al. (2008) deduce a value of × 105 for this dimensionless parameter in the host stars and × 106 in the planets.
A much higher value of stellar dissipation parameter is predicted, however, by a quantitative theory for the dissipation of internal waves in the radiative interior of the solar type stars by Barker & Ogilvie (2010), although it does not include the dissipation due to interaction of the tides with flows in the convective envelope, which is a more complicated process. On the other hand, a lower value, Q′ = 3.6 × 104, has been inferred for Jupiter from astrometric analysis of the orbital motions of the Galilean moons by Lainey et al. (2009), and a value was adopted for exoplanets by Weidner & Horne (2010), who noted that the constraints of detections of moons around hot Jupiters could be used to infer a lower limit on the planetary dissipation parameter.
Fig. 2 Radius excess δ vs. orbital distance from the host star of 103 planets for which mass and radius of star and planet are given in the “Extrasolar Planets Encyclopaedia”. |
The tides raised on the star are a significant sink of orbital angular momentum. If the mass of the companion is low, its angular momentum can be too small to keep the host star corotating with the orbit (Darwin instability), and the planet will eventually spiral into the star. The time scale for this is also known (cf. Debes & Jackson 2010), subject again to the uncertainty of the parameter . If mass and radius of both the planet and the host star are known, the expected life time of the planet before spiral-in can be calculated. The required data are available for a sample of 103 transiting planets listed in the “The Extrasolar Planets Encyclopaedia” (Schneider 2011). The time to spiral in from a distance a0 is given by (from Eq. (4) in Debes & Jackson) (1)We define the excess δ of the planet’s radius Rp relative to its equilibrium radius R0(Mp) as (2)For R0(Mp) we take the mass-radius relation in Fortney et al. (2011). Figure 3 shows the correlation between spiral-in life times t thus calculated and the radius excess δ if a value is assumed. The correlation between these two quantities does not depend on the value of Q∗' if the dissipation parameter is a constant.
Fig. 3 Radius excess δ vs. spiral-in time scale for the sample of Fig. 2 for an assumed dissipation parameter . |
The spiral-in time could be an indicator of the planet’s age by the statistical argument that a population of observed objects that are about to disappear in the near future is probably not very old. In Fig. 4 we show the age distribution obtained with a constant value of . This distribution would suggest that about half of the hot Jupiters have formed within the last 1 Gyr of the age of the Milky Way. However, the usefulness of spiral-in time as an age indicator is limited by its sensitivity to the value of , which is highly uncertain. On the other hand, the correlation shown in Fig. 3 could come from other effects, such as the trend that planets with stronger irradiation tend to have larger radii as discussed by Laughlin et al. (2011), but this would not eliminate the conclusion that inflated planets would be young if has values of about 106 or lower.
Spiral-in times are compared with standard cooling curves of irradiated planets in Fig. 5. Since spiral-in time measures age only in an average sense, we have binned the data from Fig. 3, sorting the sample by the value of δ in groups of 4 and averaging their values of δ. The ages of the bins are calculated as geometric means. Curves show the predicted radius excess as a function of age from the irradiated models in Fig. 6 of Baraffe et al. (2003). The overall agreement between the time scales derived from the tidal evolution and from the cooling tracks is poor for the adopted value of . A good agreement would require a value of × 104. This value is much lower than the single-value time-averaged estimates from ab initio calculations by Ogilvie & Lin (2007). More details on the dependence of tidal evolution with Q∗' can be found in Carone & Patzold (2007).
Fig. 5 Data of Fig. 3 binned and compared with planet cooling tracks from Baraffe et al. (2003). Solid: track for 10 MJ, dotted 1 MJ. |
2.3. Host star lithium and rotation
If the most inflated hot Jupiters are indeed the youngest planets in the sample, there could be correlations between radius excess and other age indicators such as the lithium abundance or rotation period of the host stars. The possibility that hot Jupiters could be young has been discussed by Gandolfi et al. (2010). They find that the size of CoRoT-11b fits standard cooling models if its age is only 12 Myr. They reject this possibility because CoRoT-11 has depleted its lithium. The depletion of lithium is expected, however, if the star has been formed in a binary merger because simulations of stellar collisions indicate that there is a substantial amount of chemical mixing (Trac et al. 2007). The lithium of each of the stars in the contact binary would have been depleted during the pre-main sequence and main sequence evolution of the stars prior to the merging event (cf. Martin 1997). Observations of cataclysmic binaries have indicated that lithium depletion is efficient during mass transfer evolution (Martin et al. 1995). Thus, the binary merger scenario leads to the conclusion that the host stars of hot Jupiters should have enhanced lithium depletion, and this is consistent with observational data (Israelian et al. 2004; Takeda et al. 2007).
After the merger event, the resulting single star would be a fast, magnetically active rotator surrounded by a disk. If the most inflated hot Jupiters are young, in the age range from 10 to 100 Myr, it could be expected that their host stars would be fast rotators because they still did not have enough time to slow down due to magnetic braking. Some examples of fast rotation have been noted in the literature, such as CoRoT-11 (Gandolfi et al. 2010) and HD 15082 (Collier Cameron et al. 2010). To check whether there is a correlation between planet radius excess and host star rotation we have plotted the data in Fig. 6. No clear correlation is seen, but we note that all the planets with radius excess below about 20% have host stars with rotation periods longer than five days.
Though the hosts rotate more rapidly than normal stars of the same inferred age (Pont 2009; Hartman 2010), the difference is not as large as might be expected for post-merger ages of 10−100 My. We tentatively suggest that the difference is due to the compact nature of a post-merger disk. Its specific angular momentum is much lower than in the large protostellar disks of T Tauris, and with a mass on the order of a tenth of a solar mass, it might be more efficient at braking the host star rotation through interaction with its magnetosphere.
Fig. 6 Radius excess δ vs. rotational period (days) for the host stars of 103 planets for which mass and radius of star and planet are given in the “Extrasolar Planets Encyclopaedia”. Open diamonds: stars earlier than F5 with planets with null orbital eccentricity; filled diamonds: stars earlier than F5 with significantly eccentric planet orbits; open circles: stars later than F5 with planets with zero orbital eccentricity; filled circles stars later than F5 that host planets with significant eccentricities. |
2.4. Blue stragglers
The expected eventual merging of a W UMa star would produce a single star that would look younger than a star of the same age as the original binary system. Such mergers are the default origin of the blue stragglers (BSS) in globular clusters (main sequence stars above the turn-off point), as has been proposed already by McCrea (1964) and studied by Mateo et al. (1990) and Stępień (1995, 2009). Some BSS are found to be fast rotators, but the distribution of rotation velocities peaks at a modest Vrotsini = 7 km s-1 in the 47 Tuc globular cluster (Ferraro et al. 2006). In Fig. 7 we compare the distribution of vsini values measured for BSS in 47 Tuc with that of the transiting hot Jupiters for the overlapping 0 mass range between 1.0 and 1.3 solar masses. A multiplicative correction factor was applied to the BSS data to account for their unknown inclination angle following the prescription of Chandrasekar & Munch (1950). A two-sample Kolmogorov-Smirnov test indicates that the probability that these two populations are identical is 88.3%. The statistical similarity is striking, although we note that there seems to be trend for higher vsini values among the BSS in 47 Tuc. Thus, in line with the conclusions from the previous section, we identify the slow rotation of some hot Jupiter host stars as the main pitfall of the binary merger hypothesis.
Fig. 7 Comparison of Vsini distributions for host star of transiting hot Jupiters (dashed line) and BSS (solid line) in the 47 Tuc globular cluster for the stellar mass range between 1.0 and 1.3 solar mass. A correction factor has been applied to account for the unknown inclination angle of the BSS. |
In stellar mergers, accretion of chemically fractionated material may be enhanced over normal stars as shown by Desidera et al. (2007) in the case of HD 113984 AB, a wide binary system where the primary is a BSS with an iron content 0.25 dex lower than the secondary star. A subpopulation of C and O depleted BSS has been identified in the 47 Tuc cluster by Ferraro et al. (2006), indicating again that stellar mergers can sometimes produce chemical abundance anomalies.
3. Final remarks
We have presented the working hypothesis that hot Jupiter planets could form in binary mergers. The final merger of W UMa stars is a possibility, but the direct merger of lower mass binaries that have not gone through an extended W UMa-type contact phase is likely to contribute as well. We estimate that the frequency of these binaries in the Milky Way is consistent with what would be needed to account for the observed frequency of hot Jupiters.
To test our hypothesis we looked for a correlation between the radius anomaly and the life time of the planet before spiral-in by tidal interaction with the host star. Using a value for the stellar dissipation constant we found that about half of the hot Jupiters could have ages of about 1 Gyr or younger. However, the usefulness of this age indicator is limited by its strong dependence on the value of , which is highly uncertain.
A merger scenario agrees with the peculiar distribution of orbital distances of transiting planets. The concentration near 0.05 AU, with a rapid decline at larger distances, indicates an origin close to the host star. Merger events are likely to produce planets even closer in as well, but these would have disappeared rapidly from the population owing to the steep dependence of spiral-in time on distance.
To bring their radii in agreement with standard cooling curves (Baraffe et al. 2003), the most inflated planets would need to be younger than about 0.1 Gyr. A correlation between degree of inflation and the rotation rate of the host star would be expected from the merger scenario. It is found to be only weakly present in the observations, overshadowed by large scatter. In agreement with the scenario, however, the average rotation rate of the hosts of inflated planets is greater than in normal stars of the same (inferred) age, as noted previously (Pont 2009; Hartman 2010).
The inclined and retrograde orbits identified through the Rossiter-McLaughlin effect (cf. Triaud et al. 2010) require explanation in the merger interpretation, since the angular momentum of a disk formed in a merger would have the same direction as that of the star. We suggest here that planet formation in the massive, compact disk resulting from merger on a dynamical time scale, would have produced a number of planets in closely packed orbits. Gravitational interaction in such compact systems of planets causes their orbits to be secularly unstable. This would scatter them into stochastic orbits (e.g., Chatterjee et al. 2008), including retrograde ones. (The total angular momentum vector and a corresponding preponderance of prograde orbits would be conserved.) The process is also consistent with the high frequency of excentric orbits.
Main sequence GK stars with super-solar metallicity have a markedly increased frequency of hot Jupiters (Gonzalez et al. 2001; Santos et al. 2004; Fischer & Valenti 2005) and a similar trend is seen in M dwarfs (Johnson & Apps 2009). Another intriguing observation is the presence of a massive planet around a very metal-poor star (Setiawan et al. 2011). These strange observational results could be explained in the context of stellar mergers, which may lead to unexpected abundance patterns in the resulting host stars and will form planets even around the most metal-poor stars in the Universe.
The frequency of stellar companions around contact binaries has been found to be 31% ± 6% using adaptive optics assisted imaging (Ruckinski et al. 2007). in the separation range of 3 to 100 AU. The survival of these companions will depend on their mass, their orbital separation, and the amount of mass loss from the system during the merger event. The evolution of these triple systems needs to be modeled to be able to compare with the binary properties of hot Jupiter host stars so that a useful constrain can be derived.
Our main conclusion is that the hypothesis that hot Jupiters may result from binary mergers has some advantages over other competing explanations of the radius anomaly, and it may also account for other odd observational trends. However, it has a clear issue with the slow rotation of host stars of inflated planets. We note that it would be worthwhile to study the detailed mechanisms of angular momentum loss and disk formation in binary mergers.
As observational tests of the binary merger scenario for the formation of hot Jupiters we propose the following: a) the frequency of these planets among BSS should be significantly higher than among T Tauri stars and b) the age of some exoplanet hosts inferred from nucleocosmochronology could be older than the age estimated from standard evolutionary models for single stars.
Acknowledgments
We thank the referee for detailed and constructive criticism. This work has been supported by the Spanish Plan Nacional de Astronomia y Astrofisica under grants AYA2007-2357 and AYA 2010-21308-C03-02, by RoPACS, a Marie Curie Initial Training Network funded by the European Commission’s Seventh Framework Program, and by the CONSOLIDER-INGENIO GTC project. We thank Michele Montgomery for detailed comments on an earlier version of the text. Part of this work was carried out while EM was a visiting research professor at the Geological Sciences Department of the University of Florida in Gainesville.
References
- Banit, M., Ruderman, N. A., Shaham, J., & Applegate, J. H. 1993, ApJ, 415, 779 [NASA ADS] [CrossRef] [Google Scholar]
- Baraffe, I., Chabrier, G., Barman, T. S., Allard, F., & Hauschildt, P. H. 2003, A&A, 402, 701 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Barker, A. J., & Ogilvie, G. I. 2010, MNRAS, 404, 1849 [NASA ADS] [Google Scholar]
- Batygin, K., & Stevenson, D. J. 2010, ApJ, 714, L238 [NASA ADS] [CrossRef] [Google Scholar]
- Boenheimer, P., Lin, D. N. C., & Mardling, R. A. 2001, ApJ, 548, 466 [NASA ADS] [CrossRef] [Google Scholar]
- Bouchy, F., Deleuil, M., Guillot, T., et al. 2011, A&A, 525, A68 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Burrows, A., Heng, K., & Nampaisarn, T. 2011, ApJ, 736, 47 [NASA ADS] [CrossRef] [Google Scholar]
- Butler, R. P., & Marcy, G. 1996, ApJ, 464, L156 [Google Scholar]
- Carone, L., & Patzold, M. 2007, Planet. Space Sci., 55, 643 [NASA ADS] [CrossRef] [Google Scholar]
- Chabrier, G., Leconte, J., & Baraffe, I. 2010, in The Astrophysics of Planetary Systems, ed. A. Sozzetti, et al., IAU Symp., 276, in press [arXiv:1012.0765] [Google Scholar]
- Chandrasekhar, S., & Munch, G. 1950, ApJ, 111, 142 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
- Chatterjee, S., Ford, E. B., Matsumura, S., & Rasio, F. A. 2008, ApJ, 686, 580 [NASA ADS] [CrossRef] [Google Scholar]
- Cameron, A. C., Guenther, E., Smalley, B., et al. 2010, MNRAS, 407, 507 [NASA ADS] [CrossRef] [Google Scholar]
- Coughlin, J. L., López-Morales, M., Harrison, T. E., Ule, N., & Hoffman, D. I. 2011, AJ, 141, 78 [NASA ADS] [CrossRef] [Google Scholar]
- Darwin, G. H. 1879, Roy. Soc. London Proc. Ser. I, 29, 168 [Google Scholar]
- Debes, J. H., & Jackson, B. 2010, ApJ, 723, 1703 [NASA ADS] [CrossRef] [Google Scholar]
- Desidera, S., Gratton, R. G., Lucatello, S., Endl, M., & Udry, S. 2007, A&A, 462, 1039 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- D’Souza, M. C. R., Motl, P. M., Tohline, J. E., & Frank, J. 2006, ApJ, 643, 381 [NASA ADS] [CrossRef] [Google Scholar]
- Ferraro, F. R., Sabbi, E., Gratton, R., et al. 2006, ApJ, 647, L53 [NASA ADS] [CrossRef] [Google Scholar]
- Fischer, D. A., & Valenti, J. 2005, ApJ, 622, 1102 [NASA ADS] [CrossRef] [Google Scholar]
- Ford, E. B., & Rasio, F. A. 2008, ApJ, 686, 621 [NASA ADS] [CrossRef] [Google Scholar]
- Fortney, J. J., Ikoma, M., Nettelman, N., Guillot, T., & Marley, M. S. 2011, ApJ, 729, 32 [NASA ADS] [CrossRef] [Google Scholar]
- Fusi Pecci, F., Ferraro, F. R., Corsi, C. E., Cacciari, C., & Buonnano, R. 1992, AJ, 104, 1831 [NASA ADS] [CrossRef] [Google Scholar]
- Gandolfi, D., Hébrard, G., Alonso, R., et al. 2010, A&A, 524, A55 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Gazeas, K., & Stępień, K. 2008, MNRAS, 390, 1577 [NASA ADS] [Google Scholar]
- Goldreich, P. 1963, MNRAS, 126, 257 [NASA ADS] [Google Scholar]
- Gonzalez, G., Laws, C., Tyagi, S., & Reddy, B. E. 2001, AJ, 121, 432 [NASA ADS] [CrossRef] [Google Scholar]
- Guillot, T. 2008, Phys. Scripta, 130, 014023 [Google Scholar]
- Guillot, T., Santos, N. C., Pont, F., et al. 2006, A&A, 453, L21 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Hartman, J. D. 2010, ApJ, 717, L138 [NASA ADS] [CrossRef] [Google Scholar]
- Hut, P. 1980, A&A, 92, 167 [NASA ADS] [Google Scholar]
- Ida, S., & Lin, D. 2010, ApJ, 719, 810 [NASA ADS] [CrossRef] [Google Scholar]
- Israelian, G., Santos, N. C., Mayor, M., & Rebolo, R. 2004, A&A, 414, 601 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Jackson, B., Greenberg, R., & Barnes, R. 2008, ApJ, 678, 1396 [NASA ADS] [CrossRef] [Google Scholar]
- Jackson, B., Barnes, R., & Greenberg, R. 2009, ApJ, 698, 1357 [NASA ADS] [CrossRef] [Google Scholar]
- Johnson, J. A., & Apps, K. 2009, PASP, 699, 933 [Google Scholar]
- Lainey, V., Arlot, J.-E., Karatekin, O., & van Hoolst, T. 2009, Nature, 459, 957 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
- Laughlin, G., Crismani, M., & Adams, F. C. 2011, ApJ, 729, L7 [NASA ADS] [CrossRef] [Google Scholar]
- Leconte, J., Baraffe, I., Chabrier, G., Barman, T., & Levrard, B. 2009, A&A, 506, 385 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Lee, B. L., Ge, J., Fleming, S. W., et al. 2011, ApJ, 728, 32 [NASA ADS] [CrossRef] [Google Scholar]
- Madhusudhan, N., Harrington, J., Stevenson, K. B., et al. 2011, Nature, 469, 64 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
- Martin, E. L. 1997, Mem. Soc. Ast. It., 68, 905 [Google Scholar]
- Martin, E. L., Casares, J., Charles, P. A., & Rebolo, R. 1995, A&A, 303, 785 [NASA ADS] [Google Scholar]
- Mateo, M., Harris, H. C., Nemec, J., & Olszewski, E. W. 1990, AJ, 100, 469 [NASA ADS] [CrossRef] [Google Scholar]
- Mayor, M., & Queloz, D. 1995, Nature, 378, 355 [NASA ADS] [CrossRef] [Google Scholar]
- McCrea, W. H. 1964, MNRAS, 128, 147 [NASA ADS] [CrossRef] [Google Scholar]
- Mordasini, C., Alibert, Y., Benz, W., & Naef, D. 2009, A&A, 501, 1161 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Naoz, S., Farr, W. M., Lithwick, Y., Rasio, F. A., & Teyssandier, J. 2011, Nature, 473, 187 [NASA ADS] [CrossRef] [Google Scholar]
- Ogilvie, G. I., & Lin, D. N. C. 2007, ApJ, 661, 1180 [NASA ADS] [CrossRef] [Google Scholar]
- Perets, H. B. 2011, in Planetary Systems Beyond the Main Sequence, AIP Conf. Proc., 1331, 56 [NASA ADS] [CrossRef] [Google Scholar]
- Pollack, J. B., & Bodenheimer, P. 1989, in Origin and Evolution of Planetary and Satellite Atmospheres, ed. S. K. Atreya, J. B. Pollack, & M. S. Matthews (Tucson: Univ. Arizona Press), 564 [Google Scholar]
- Pont, F. 2009, MNRAS, 396, 1789 [NASA ADS] [CrossRef] [Google Scholar]
- Rasio, F. A. 1995, ApJ, 444, L41 [NASA ADS] [CrossRef] [Google Scholar]
- Rasio, F. A., & Shapiro, S. L. 1995, ApJ, 438, 887 [NASA ADS] [CrossRef] [Google Scholar]
- Ritter, H. 1988, A&A, 202, 93 [NASA ADS] [Google Scholar]
- Ritter, H. 1996, Evolutionary Processes in Binary Stars, NATO ASIC Proc., 477, 223 [Google Scholar]
- Rucinski, S. M. 2002, PASP, 114, 1124 [NASA ADS] [CrossRef] [Google Scholar]
- Rucinski, S. M. 2006, MNRAS, 368, 1319 [NASA ADS] [CrossRef] [Google Scholar]
- Rucinski, S. M., Pribulla, Th., & van Kerkwijk, M. H. 2006, AJ, 134, 2353 [Google Scholar]
- Stępień, K. 1995, MNRAS, 274, 1019 [NASA ADS] [Google Scholar]
- Stępień, K. 2006, Acta Astron., 56, 347 [NASA ADS] [Google Scholar]
- Stępień, K. 2009, MNRAS, 397, 857 [NASA ADS] [CrossRef] [Google Scholar]
- Stępień, K. 2011, Acta Astron., 61, 139 [NASA ADS] [Google Scholar]
- Stevenson, D. J. 1991, ARA&A, 29, 163 [NASA ADS] [CrossRef] [Google Scholar]
- Takeda, Y., Kawanomoto, S., Honda, S., Ando, H., & Sakurai, T. 2007, A&A, 468, 663 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Trac, H., Sills, A., & Penn, U.-E. 2007, MNRAS, 377, 977 [NASA ADS] [CrossRef] [Google Scholar]
- Triaud, A. H. M. J., Collier Cameron, A., Queloz, D., et al. 2010, A&A, 524, A25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Tutukov, A. V. 1991, AZh, 68, 837 [NASA ADS] [Google Scholar]
- Tutukov, A. V., Dremova, G. N., & Svechnikov, M. A. 2004, Astron. Rep., 48, 219 [Google Scholar]
- Tylenda, R., Hajduk, M., Kamiński, T., et al. 2010, A&A, 528, A114 [Google Scholar]
- Vidotto, A. A., Jardine, M., & Helling, C. 2011, MNRAS, 411, L46 [NASA ADS] [CrossRef] [Google Scholar]
- Weidner, C., & Horne, K. 2010, A&A, 521, A76 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Wisniewski, J. O., Clampin, M., Bjorkman, K., & Barry, R. 2008, ApJ, 683, L171 [NASA ADS] [CrossRef] [Google Scholar]
- Wolszczan, A., & Kuchner, M. 2010, Exoplanets, ed. S. Seager (Tucson: Univ. of Arizona Press), 175 [Google Scholar]
- Youdin, A. N., & Mitchell, J. L. 2010, ApJ, 721, 1113 [NASA ADS] [CrossRef] [Google Scholar]
- Zanni, C., & Ferreira, J. 2011, ApJ, 727, L22 [NASA ADS] [CrossRef] [Google Scholar]
- Zapolski, H. S., & Salpeter, E. E. 1969, ApJ, 158, 809 [NASA ADS] [CrossRef] [Google Scholar]
All Figures
Fig. 1 Masses of the host stars of 103 transiting planets for which mass and radius of star and planet are given in the “Extrasolar Planets Encyclopaedia” compared with the masses of W UMa binaries from Gazeas & Stępień (2008). |
|
In the text |
Fig. 2 Radius excess δ vs. orbital distance from the host star of 103 planets for which mass and radius of star and planet are given in the “Extrasolar Planets Encyclopaedia”. |
|
In the text |
Fig. 3 Radius excess δ vs. spiral-in time scale for the sample of Fig. 2 for an assumed dissipation parameter . |
|
In the text |
Fig. 4 Distribution of spiral-in times obtained from the results presented in Fig. 3. |
|
In the text |
Fig. 5 Data of Fig. 3 binned and compared with planet cooling tracks from Baraffe et al. (2003). Solid: track for 10 MJ, dotted 1 MJ. |
|
In the text |
Fig. 6 Radius excess δ vs. rotational period (days) for the host stars of 103 planets for which mass and radius of star and planet are given in the “Extrasolar Planets Encyclopaedia”. Open diamonds: stars earlier than F5 with planets with null orbital eccentricity; filled diamonds: stars earlier than F5 with significantly eccentric planet orbits; open circles: stars later than F5 with planets with zero orbital eccentricity; filled circles stars later than F5 that host planets with significant eccentricities. |
|
In the text |
Fig. 7 Comparison of Vsini distributions for host star of transiting hot Jupiters (dashed line) and BSS (solid line) in the 47 Tuc globular cluster for the stellar mass range between 1.0 and 1.3 solar mass. A correction factor has been applied to account for the unknown inclination angle of the BSS. |
|
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.