Formation, Dynamical Evolution, and Habitability of Planets in Binary Star Systems
A survey of currently known planet-hosting stars indicates that approximately 25% of extrasolar planetary systems are within dual-star environments. Several of these systems contain stellar companions on moderately close orbits, implying that studies of the formation and dynamical evolution of giant and terrestrial planets, in and around binary star systems have now found realistic grounds. With the recent launch of the space telescope CoRoT, and the launch of NASA’s Kepler satellite in 2009, the number of such dynamically complex systems will soon increase and many more of their diverse and interesting dynamical characteristics will soon be discovered. It is therefore, both timely and necessary, to obtain a deep understanding of the history and current status of research on planets in binary star systems. This chapter will serve this purpose by reviewing the models of the formation of giant and terrestrial planets in dual-star environments, and by presenting results of the studies of their dynamical evolution and habitability, as well as the mechanisms of delivery of water and other volatiles to their terrestrial-class objects. In this chapter, the reader is presented with a comprehensive, yet relatively less technical approach to the study of planets in and around binary stars, and with discussions on the differences between dynamical characteristics of these systems and planetary systems around single stars.
The concept of a “world with two suns” has been of interest to
astronomers for many years. Many scientists tried to understand
whether planets could form in binary star systems, and whether
the notion of habitability, as we know it, could be extended
to such environments. Although as a result of their respective works,
many dynamical features of binary-planetary systems
The discovery of extrasolar planets during the past decade has, however,
changed this trend. Although the candidate planet-hosting stars have
been routinely chosen to be single, or within wide (100 AU)
The first detection of a planet in a binary system was reported by Campbell, Walker & Yang in 1988. In an attempt to identify planetary objects outside our solar system, these authors measured the variations in the radial velocities of a number of stars, and reported the possibility of the presence of a Jovian-type body around the star Cephei (figure 9.1, Campbell, Walker & Yang, 1988). This star, that is a K1 IV sub-giant with a mass of 1.59 solar-masses (Fuhrmann, 2004), is the primary of a binary system with a semimajor axis of 18.5 AU and an eccentricity of 0.36 (Griffin, Carquillat, Ginestet, 2002; Hatzes et al., 2003). The secondary of this system is an M dwarf with a mass of 0.44 solar-masses (Neuhäuser et al., 2007; Torres, 2007). Initial radial velocity measurements of Cephei implied that this star may be host to a giant planet with a probable mass of 1.7 Jupiter-masses, in an orbit with a semimajor axis of 1.94 AU (Campbell, Walker & Yang, 1988).
Unfortunately, the discovery of the first binary-planetary system, which could have also marked the detection of the first planet outside our solar system, did not withstand skepticism. In an article in 1992, Walker and his colleagues attributed their measured variations of the radial velocity of Cephei to the chromospheric activities of this star, and announced that the possibility of the existence of a giant planet around Cephei may be none (Walker et al., 1992). It took observers an additional 12 years to monitor Cephei and measure its radial velocity to arrive at the conclusion that the previously observed variations were not due to stellar activities and were in fact representative of a planetary companion (figure 9.2). It was the discovery of a giant planet in the binary system of GL 86 (Queloz et al., 2000), and the (re-)announcement of the detection of a giant planet in Cephei system (Hatzes et al., 2003) that opened a new chapter in the theoretical and observational studies of extrasolar planetary systems.
The fact that giant planets exist in moderately close ( AU) binary systems has confronted dynamicists with many new challenges. Questions such as, how are these planets formed, can binary-planetary systems host terrestrial and/or habitable planets, how are habitable planets formed in such dynamically complex environments, and how do such planets acquire the ingredients necessary for life, are among major topics of research in this area. This chapter is devoted to review these issues and present the current status of research on the formation of planet in dual-star systems and habitability of terrestrial bodies in and around binary stars.
The chapter begins with a review of the dynamics of a planet in a binary star system. In general, prior to constructing a theory for the formation of planets, it proves useful to study whether the orbit of a planet around its host star would be stable. In a binary system, such studies are of quite importance since in these systems the perturbation of the stellar companion may dictate the possibility of the formation of planetary bodies by affecting the stability and dynamics of smaller objects.
The formation of planets in binary star systems is reviewed in the third section. Although the study of the dynamics of planets in and around binary stars dates back to approximately forty year ago, the formation of planets in these systems is an issue that is still unresolved. In spite the observational evidence that indicate majority of main and pre-main sequence stars are formed in binaries or clusters (Abt, 1979; Duquennoy & Mayor, 1991; Mathieu, 1994; Mathieu et al., 2000; White Ghez, 2001), and in spite the detection of potentially planet-forming environments in and around binary stars (figure 9.3, also see Mathieu, 1994; Akeson, Koerner Jensen, 1998; Rodriguez et al., 1998; White et al., 1999; Silbert et al., 2000; Mathieu et al., 2000), planet formation theories are still unclear in explaining how planets may form in multi-star environments. The focus of section 9.3 is on discussing the formation of giant and terrestrial planets in moderately close binary-planetary systems, and reviewing the current status of planet formation theories in this area.
The habitability of a binary system is presented in section 9.4. The notion of habitability is defined based on the habitability of Earth and life, as we know it. Such a definition requires a habitable planet to have the capability of retaining liquid water in its atmosphere and on its surface. The latter is determined by the luminosity of the central star, the size of the planet, and also the distribution of water in the protoplanetary disk from which terrestrial-class objects are formed. In section 9.4, a review of the current status of the models of habitable planet formation in and around binary systems are presented, and their connections to
models of terrestrial planet formation and water-delivery around single stars are discussed. Finally the chapter ends by discussing the future prospects of research in the field of planets in binaries.
0.2 Dynamical Evolution and Stability
In general, one can consider orbital stability synonymous with the capability of an object in maintaining its orbital parameters (i.e., semimajor axis, eccentricity, and inclination) at all times. In other words, an object is stable if small variations in its orbital parameters do not progress exponentially, but instead vary sinusoidally. Instability occurs when the perturbative forces create drastic changes in the time variations of these parameters and result in either the ejection of the object from the system (i.e., leaving the system’s gravitational field), or its collision with other bodies.
The concept of stability, as explained above, although simple, has been defined differently by different authors. A review by Szebehely in 1984 lists 50 different definitions for the stability of an object in a multi-body system (Szebehely, 1984). For instance, Harrington (1977) considered the orbit of an object stable if, while numerically integrating the object’s orbit, its semimajor axis and orbital eccentricity would not undergo secular changes. Szebehely (1980), and Szebehely McKenzie (1981), on the other hand, considered integrals of motion and curves of zero velocity to determine the stability of a planet in and around binary stars. In this chapter, the values of the orbital eccentricity of an object and its semimajor axis are used to evaluate its stability. An object is considered stable if, for the duration of the integration of its orbit, the value of its orbital eccentricity stays below unity, it does not collide with other bodies, and does not leave the gravitational field of its host star.
The study of the stability of a planetary orbit in dual stars requires a detailed analysis of the dynamical evolution of a three-body system. Such an analysis itself is dependent upon the type of the planetary orbit. In general, a planetary-class object may have three types of orbit in and around a binary star system. Szebehely (1980) and Dvorak (1982) have divided these orbits into three different categories. As indicated by Szebehely (1980), a planet may be in an inner orbit, where it revolves around the primary star, or it may be in a satellite orbit, where it revolves around the secondary star. A planet may also be revolving the entire binary system in which case its orbit is called an outer orbit. As classified by Dvorak (1982), on the other hand, a study of the stability of resonant periodic orbits in a restricted, circular, three-body system indicates that a planet may have an S-type orbit, where it revolves around only one of the stars of the binary, or may be in a P-type orbit, where it revolves the entire binary system (figure 9.4). A planet may also be in an L-type orbit where it librates in a stable orbit around the or Lagrangian points.
The rest of this section is devoted to a review of the studies of the stability of planets in S-type and P-type orbits. Given that in a binary star system, a planet is subject to the gravitational attraction of two massive bodies (i.e., the stars), it would be important to understand how the process of the formation of planetary objects, and the orbital dynamics of small bodies would be affected by the orbital characteristics of the binary’s stellar components. In general, except for some special cases for which analytical solutions may exist, such studies require numerical integrations of the orbits of all the bodies in the system. In the past, prior to the invention of symplectic integrators (Wisdom & Holman, 1991), which enabled dynamicists to extend the studies of the stability of planetary systems to several hundred million years, and before the development of fast computers, majority of such studies were either limited to those special cases, or were carried out numerically for only a small number of binary’s orbital period. Examples of such studies can be found in articles by Graziani Black (1981), Black (1982), and Pendleton Black (1983), in which the authors studied the orbital stability of a planet in and around a binary star. By numerically integrating the equations of the motion of the planet, this authors showed that, when the stars of the binary have equal masses, the orbital stability of the planet is independent of its orbital inclination (also see Harrington, 1977). Their integrations also indicate that, when the mass of one of the stellar components is comparable to the mass of Jupiter, planetary orbits with inclinations higher than 50 tend to become unstable.
The invention of symplectic integrators, in particular routines that have
been designed specifically for the purpose of integrating orbits
of small bodies in dual-star systems
0.2.1 Stability of S-type orbits
As mentioned above, instability occurs when the perturbative effects cause the semimajor axis and orbital eccentricity of a planet change in such a way that either the object leaves the gravitational field of the system, or it collides with another body. For a planet in an S-type orbit, the gravitational force of the secondary star is the source of these perturbations. That implies, a planet at a large distance from the secondary, i.e., in an orbit closer to its host star, may receive less perturbation from the binary companion and may be able to sustain its dynamical state for a longer time (Harrington, 1977). Since the perturbative effect of the stellar companion varies with its mass, and the eccentricity and semimajor axis of the binary (which together determine the closest approach of the secondary to the planet), it is pos-
sible to estimate an upper limit for the planet’s distance to the star beyond which the orbit of the planet would be unstable. As shown by Rabl Dvorak (1988) and Holman Wiegert (1999), the maximum value that the semimajor axis of a planet in an S-type orbit can attain and still maintain its orbital stability is a function of the mass-ratio and orbital elements of the binary, and is given by (Rabl Dvorak, 1988; Holman Wiegert, 1999)
In this equation, is critical semimajor axis , , and are the semimajor axis and eccentricity of the binary, and and are the masses of the primary and secondary stars, respectively. The signs in equation (9.1) define a lower and an upper value for the critical semimajor axis , and set a transitional region that consists of a mix of stable and unstable systems. Such a dynamically gray area, in which the state of a system changes from stability to instability, is known to exist in multi-body environments, and is a characteristic of any dynamical system.
Equation (9.1) is an empirical formula that has been obtained by numerically integrating the orbit of a test particle (i.e., a massless object) at different distances from the primary of a binary star (Rabl Dvorak, 1988; Holman Wiegert, 1999). Figure 9.5 shows this in more detail. Similar studies have been done by Moriwaki & Nakagawa (2004), and Fatuzzo et al. (2006)
who obtained critical semimajor axes slightly larger than given by equation (9.1).
Since the mass of a Jovian-type planet is approximately three
orders of magnitude smaller than the mass of a star, such a
test particle approximation yields
results that not only are
applicable to the stability of giant planets, but
can also be used in identifying regions where
smaller bodies, such as terrestrial-class objects
(Quintana et al., 2002; Quintana & Lissauer, 2006; Quintana et al., 2007)
and dust particles (Trilling et al., 2007),
can have long term stable orbits
The stability of S-type systems has been studied by many authors (Benest, 1988, 1989, 1993, 1996; Wiegert Holman, 1997; Pilat-Lohinger Dvorak, 2002; Dvorak et al., 2003, 2004; Pilat-Lohinger et al., 2004; Musielak et al., 2005). In a recent article, Haghighipour (2006) extended such studied to the dynamical stability of the Jupiter-like planet of the Cephei planetary system. By numerically integrating the orbit of this object for different values of and (the orbital inclination of the giant planet relative to the plane of the binary), Haghighipour (2006) has shown that the orbit of this planet is stable for the values of the binary eccentricity within the range . Figure 9.6 shows the results of such integrations for a coplanar system with and for different values of the binary eccentricity. The initial value of the semimajor axis of the binary was chosen to be
21.5 AU. Integrations also indicated that the binary-planetary system of Cephei becomes unstable in less than a few thousand years when the initial value of the binary eccentricity exceeds 0.5.
Interesting results were obtained when the Cephei system was integrated for different values of . The results indicated that for the above-mentioned range of orbital eccentricity, the planet maintains its orbit for all values of inclination less than 40. Figure 9.7 shows the semimajor axes and orbital eccentricities of the system for and for =5, 10, and 20. For orbital inclinations larger than 40, the system becomes unstable in a few thousand years.
An interesting dynamical phenomenon that may occur in an S-type binary, and has also been observed in the numerical simulations of a few of these systems, is the Kozai resonance (Haghighipour, 2003, 2005; Verrier & Evans, 2006; Takeda & Rasio, 2006; Malmberg, Davies & Chambers, 2007). As demonstrated by Kozai (1962), in a three-body system with two massive objects and a small body (e.g., a binary-planetary system), the exchange of angular momentum between the planet and the secondary star, can cause the orbital eccentricity of the planet to reach high values at large inclinations. Averaging the equations of motion of the system over mean anomalies, one can show that in this case, the averaged system is integrable when the ratio of distances are large (the Hill’s approximation, Kozai, 1962). The Lagrange equations of motion in this case, indicate that, to the first order of planet’s eccentricity, the longitude of the periastron of this object, , librates around a fix value. Figure 9.8 shows this for the giant planet of Cephei. As shown here, librates around 90 (Haghighipour, 2003, 2005).
In a Kozai resonance, the longitude of periastron and the orbital eccentricity of the small body are related to its orbital inclination as (Innanen, 1997)
From equation (9.3), one can show that the Kozai resonance may occur if the orbital inclination of the small body is larger than 39.23. For instance, as shown by Haghighipour (2003, 2005), in the system of Cephei, Kozai resonance occurs at . For the minimum value of , the maximum value of the planet’s orbital eccentricity, as given by equation (9.4), is equal to 0.764. Figure 9.8 also shows that stays below this limiting value at all times.
As shown by Kozai (1962) and Innanen (1997), in a Kozai resonance, the disturbing function of the system, averaged over the mean anomalies, is independent of the longitudes of ascending nodes of the small object (the planet) and the perturbing body (the stellar companion). As a result, the quantity (shown as the “Reduced Delaunay Momentum” in figure 9.8) becomes a constant of motion. Since the eccentricity and inclination of the planet vary with time, the
fact that the quantity above is a constant of motion implies that the time-variations of these two quantities have similar periods and, at the same time, they vary in such a way that when reaches its maximum, reaches its minimum and vice versa. Figure 9.9 shows this clearly.
0.2.2 Stability of P-type Orbits
Numerical simulations have also been carried out for the stability of P-type orbits in binary-planetary systems (Ziglin, 1975; Szebehely McKenzie, 1981; Dvorak, 1984, 1986; Dvorak, Froeschlé, Froeschlé, 1989; Kubala, Black & Szebehely; Holman Wiegert, 1999; Broucke, 2001; Pilat-Lohinger, Funk & Dvorak, 2003; Musielak et al., 2005). Similar to S-type orbits, in order for a P-type planet to be stable, it has to be at a safe distance from the two stars so that it would be immune from their perturbative effects. That is, planets at large distances from the center of mass of a binary will have a better chance of begin stable. This distance, however, cannot be too large because at very large distances, other astronomical effects, such as galactic perturbation, and perturbations due to passing stars, can render the orbit of a planet unstable.
To determine the critical value of the semimajor axis of a P-type planet in a stable orbit, preliminary attempts were made by Dvorak (1984), who numerically integrated the orbit of a circumbinary planet in a circular orbit around an eccentric binary system and showed that planets at distances 2-3 times the separation of the binary have stable orbits. Subsequent studies by Dvorak (1986), Dvorak, Froeschlé, Froeschlé (1989), and Holman Wiegert (1999) complemented Dvorak’s results of 1984 and showed that the orbit of a P-type planet will be stable as long as the semimajor axis of the planet stays larger than the critical value given by (figure 9.10)
Similar to equation (9.1), equation (9.4) represents a transitional region with a lower boundary below which the orbit of a P-type planet will be certainly unstable, and an upper boundary beyond which the orbit of the planet will be stable. The mixed zone between these two boundaries represents a region where a planet, depending on its orbital parameters, and the orbital parameters and the mass-ratio of the binary, may or may not be stable. Recently, by applying the stability criteria of equation (9.4) to their observational results, Trilling et al. (2007) have confirmed the presence of stable dust band, possibly resulted from collision of planetesimal, around close binary star systems.
A dynamically interesting feature of the stable region around the stars of a binary is the appearance of islands of instability. As shown by Holman Wiegert (1999) islands of instability may develop beyond the inner boundary of the mixed zone, which correspond to the locations of mean-motion resonances. The appearance of these unstable regions have been reported by several authors under various circumstances (Hénon Guyot, 1970; Dvorak, 1984; Rabl Dvorak, 1988; Dvorak, Froeschlé, Froeschlé, 1989). Extensive numerical simulations would be necessary to determine whether the overlapping of these resonances would result in stable P-type binary-planetary orbits.
0.3 Planet Formation in Binaries
Despite a wealth of articles on planets in binary star systems, the process of the formation of these objects is still poorly understood. The current theories of planet formation focus only on the formation of planets in a circumstellar disk around a single star, and their extensions to binary environments are limited to either the Sun-Jupiter system, where the focus is on the effect of Jupiter on the formation of inner planets of our solar system (Heppenheimer, 1974, 1978; Drobyshevski, 1978; Diakov & Reznikov, 1980; Whitmire et al., 1998; Kortenkamp, Wetherill & Inaba), or binaries resembling some of extrasolar planets, in which the secondary star has a mass in the brown dwarf regime (Whitmire et al., 1998). Although attempts have been made to extend such studies to binaries with comparable-mass stellar components (Marzari Scholl, 2000; Nelson, 2000; Barbieri, Marzari Scholl, 2002; Quintana et al., 2002; Lissauer et al., 2004), the extent of the applicability of the results of these studies has been only to hypothetical cases since, until recently, there had been no observational evidence on the existence of such binary-planetary systems.
In general, it is believed that planet formation proceeds through the following four stages (figure 9.11):
1) coagulation of dust particles and their growth to centimeter-sized objects,
2) collisional growth of centimeter-sized particles to kilometer-sized bodies (planetesimals),
3) formation of Moon- to Mars-sized protoplanets (also known as planetary embryos) through the collision and coalescence of planetesimals, and
4) collisional growth of planetary embryos to terrestrial-sized objects.
The latter is a slow process that may take a few hundred million years. During the first few million years of this process, at larger distances from the star, planetesimals and planetary embryos may form planetary cores several times more massive than Earth, and may proceed to form giant planets.
In a binary star system with a moderate to small separation, the secondary star will have significant effects on the efficiency of each of these processes. As shown by Boss (2006), a binary companion can alter the structure of a planet-forming nebula, and create regions where the densities of the gas and dust are locally enhanced (figure 9.12). Also, as shown by Artymowicz & Lubow (1994), and
Pichardo, Sparke & Aguilar (2005), a stellar component on an eccentric orbit can truncate the circumprimary disk of embryos to smaller radii and remove material that may be used in the formation of terrestrial planets (figure 9.13). As a result, it used to be believed that circumstellar disks around the stars of a binary may not be massive enough to form planets. However, observations by Mathieu (1994), Akeson, Koerner Jensen (1998), Rodriguez et al. (1998), and Mathieu et al. (2000) have indicated that potentially planet-forming circumstellar disks can indeed exist around the stars of a binary system, implying that planet formation in binaries may be as common as around single stars (figure 9.14). Among these circumstaller disks, the two well-separated disks of the system L1551 retain equivalent of approximately 0.03 to 0.06 solar-masses of their original circumstellar materials in a region with an outer radius of 10 AU (figure 9.14, Rodriguez et al., 1998). The masses of these disks are comparable to the minimum solar-mass model of the primordial nebula of our solar system (Weidenschilling, 1977; Hayashi, 1981), implying that, planet formation in dual-star systems can begin and continue in the same fashion as around our Sun.
Despite the observational evidence in support of the existence of planet-forming environments in moderately close binary star systems, the perturbative effect of the
binary companion may not always favor planet formation. For instance, as shown by Nelson (2000), giant planet formation cannot proceed through the disk instability mechanism (Boss, 2000) around the primary of a binary star system with separation of AU. Also, when forming planetary embryos, as shown by Heppenheimer (1978), Whitmire et al. (1998), and Thébault et al. (2004), the perturbation of the secondary star may increase the relative velocities of planetesimals and cause their collisions to result in breakage and fragmentation (figure 9.15). Results of the studies by these authors suggest that planetesimal accretion will be efficient only in binaries with large separation [50 AU as indicated by Heppenheimer (1978), 26 AU as shown by Whitmire et al. (1998), and 100 AU as reported by Mayer et al. (2005)]. Finally, in a binary star system, the stellar companion may create unstable regions where the building blocks of planets will not maintain their orbits and, as a result, planet formation will be inhibited (Whitmire et al., 1998).
Interestingly, despite all these difficulties, numerical simulations have shown that it may indeed be possible to form giant and/or terrestrial planets in and around a dual-star system. Recent simulations by Boss (2006), and Mayer, Boss & Nelson (2007) indicate that Jupiter-like planets can form around the primary of a binary star system via gravitational instability in a marginally unstable circumprimary disk (figure 9.16). On the other hand, as shown by Thébault et al. (2004), the core accretion mechanism may also be able to form giant planets around the primary of a binary star. However, as the results of their simulations for planet formation in
the Cephei system indicate, the semimajor axis of the final gas-giant planet may be smaller than its observed value.
In regard to the formation of terrestrial planets in binary systems, in a series of articles, Quintana and her colleagues integrated the orbits of a few hundred Moon- to Mars-sized objects and showed that terrestrial-class objects can form in and around binaries (Quintana et al., 2002; Lissauer et al., 2004; Quintana & Lissauer, 2006; Quintana et al., 2007). Figure 9.17 shows the results of some of their simulations. As shown here, depending on the mass-ratio of the binary and the initial values of its orbital parameters, in a few hundred million years, terrestrial planets can form around a close (0.01 to 0.1 AU) binary star system.
Quintana and colleagues also studied terrestrial planet formation in binaries with larger separations (Quintana et al., 2007). Figure 9.18 shows the results of their simulations for a binary with a separation of 20 AU. Similar to figure 9.17, terrestrial-type objects are formed around the primary of the binary in a few hundred million years. Statistical analysis of their results, as shown in figures 9.19,
9.20, and 9.21 indicate that, as expected in binaries with larger perihelia, where disk truncation has been smaller and more planet-forming material is available, terrestrial planet formation is efficient and the number of final terrestrial planets is large. The results of simulations by Quintana et al. (2007) also indicate that in a binary with a periastron distance larger than 10 AU, terrestrial planet formation can proceed efficiently in a region within 2 AU of the primary star. In binaries with periastra smaller than 5 AU, this region may be limited to inside 1 AU (figure 9.20).
Despite the destructive role of the binary companion in increasing the relative velocities of planetesimals, which causes their collisions to result in erosion, this efficiency of terrestrial planet formation in binary systems may be attribute to the fact that the effect of the binary companion on increasing the relative velocities of planetesimals can be counterbalanced by dissipative forces such as gas-drag and dynamical friction (Marzari et al., 1997; Marzari Scholl, 2000). The combination of the drag force of the gas and the gravitational force of the secondary star may result in the alignment of the periastra of planetesimals and increases the efficiency of their accretion by reducing their relative velocities. This is a process that is more effective when the sizes of the two colliding planetesimals are comparable and small. For colliding bodies with different sizes, depending on the size distribution of small objects, and the radius of each individual planetesimal, the process of the alignment of periastra may instead increase the relative velocities of the two objects, and cause their collisions to become eroding (figure 9.22, also see Thébault et al., 2006).
Simulations of terrestrial planet formation have also been extended to binaries with larger separations (20-40 AU) that also host a giant planet (Haghighipour & Raymond, 2007). As discussed in the next section, by numerically integrating the orbits of the binary, its giant planet, and a few hundred planetary embryos, these authors have shown that it is possible to form Earth-like objects, with considerable amount of water, in the habitable zone of the primary of a moderately close binary-planetary system.
It is widely accepted that a planet capable of developing life (as we know it), has to be able to continuously maintain liquid water on its surface and in its atmosphere. The capability of a planet in retaining water depends on its size and the processes involving its interior dynamics and atmospheric circulation. It also depends on its orbital parameters (i.e., its semimajor axis and orbital eccentricity) and the brightness of the central star at the location of the planet. While a dynamic interior and atmospheric circulation are necessary for a habitable planet to develop CO cycle and generate greenhouse effect (which helps the planet to maintain a uniform temperature), a long-term stable orbit, at a right distance from the star, is essential to ensure that the planet will receive the amount of radiation that enables it to maintain liquid water on its surface. These seemingly unrelated characteristics of a potential habitable planet, have strong intrinsic correlations, and combined with the luminosity of the star, determine the system’s habitable zone.
The inner and outer boundaries of a habitable zone vary with the star’s luminosity and the planet’s atmospheric circulation models [see Menou & Tabachnik (2003), Jones, Underwood, Sleep (2005), and Jones, Sleep, & Underwood (2006) for a table of distances of the inner and outer boundaries of the habitable zones of exoplanetary systems]. A conservative estimate of the habitable zone of a star can be made by assuming that its inner edge is located at a distance closer than which
water on the surface of the planet evaporates due to runaway greenhouse effect, and its outer edge is at a distance where, in the absence of CO clouds, runaway glaciation will freeze the water and creates permanent ice on the surface of the planet. As shown by Kasting, Whitmire, Reynolds (1993), such a definition of a habitable zone results in a habitable region between 0.95 AU and 1.15 AU for the Sun (figure 9.23). This is a somewhat conservative estimate of the Sun’s habitable zone and as noted by Jones, Underwood, Sleep (2005), the outer edge of this region may, in fact, be farther away (Forget Pierrehumbert, 1997; Williams Kasting, 1997; Mischna et al., 2000).
Since the notion of habitability is based on life on Earth, it is possible to calculate the location of the boundaries of the habitable zone of a star by comparing its luminosity with that of the Sun. For a star with the surface temperature and radius , the luminosity is given by
where is Boltzmann constant. Using equation (9.5) and the fact that Earth is in the habitable zone of the Sun, the radial distances of the inner and outer edges of the habitable zone of a star can be obtained from (Haghighipour 2006)
The quantities and in equation (9.6) are the surface temperature and radius of the Sun, respectively, and represent the radial distance of Earth from the Sun. Equation (9.6) implies that a habitable zone can be defined as a region around a star where an Earth-like planet can receive the same amount of radiation as Earth receives from the Sun, so that it can develop and maintain similar habitable conditions as those on Earth.
As mentioned above, the orbit of a potential habitable planet in the habitable zone of a star has to be stable over long durations of time. As shown in section 9.2, the stability of the orbit of a planet in a binary system is strongly affected by the orbital motion of the binary companion. In binary systems where the primary hosts other planetary bodies (e.g., giant planets), the dynamics of a habitable planet will also be affected by the gravitational perturbations of these objects. It is therefore important to determine under what conditions a terrestrial-class object will have a long-term stable orbit in the habitable zone of a binary system, prior to construction a theory for the formation of Earth-like planets in such environments.
Since a terrestrial-class planet is approximately two orders of magnitude less massive than a Jovian-type object, it will not have significant effect on the motion of the stars and the giant planets of a binary system. Therefore, as explained in sections 9.2.1 and 9.2.2, if a binary system does not contain a Jupiter-like planet, any dynamical criterion that is obtained for the stability or instability of a general planetary body,
can also be applied to the dynamics of a terrestrial planet. Equation (9.1) and the stability conditions presented by figure 9.5 can be used to determine the long-term stability of an Earth-like object in a binary system.
If a binary contains giant planets, however, the situation is different. The gravitational perturbations of the latter objects will have significant effects on the motion and dynamics of terrestrial planets in the system. As shown by Haghighipour (2006), an Earth-size object, in a region between the giant planet and the primary of Cephei binary system, can maintain its stability
only in orbits close to the primary star and outside
the influence zone
As mentioned above, the instability of an Earth-like planet in the habitable zone of Cephei can be attributed to the interaction between this object and the giant planet of the system. When the Earth-like planet is outside the giant planet’s influence zone (e.g., at closer distances to the primray star) it can maintain its orbit for several hundred million years. Figure 9.24 suggest that, in order for a binary-planetary system to be habitable, its habitable zone has to be outside the influence region of its giant planet. In an S-type binary-planetary system, this implies a primary with a close-in habitable region. In a recent article, Haghighipour & Raymond (2007) have studies the habitability of such a system. By considering a binary with a Sun-like primary star and a Jupiter-sized planet in a circular orbit at 5 AU, and by adopting the model of Morbidelli et al. (2000), which is based on the assumption that water-carrying objects, in the Sun’s asteroid belt, were the primary source of the delivery of water to Earth, these authors integrated the orbits of a few hundred protoplanetary (Moon- to Mars-sized) objects, and showed that it is indeed possible to form Earth-sized planets, with substantial amount of water, in the habitable zone of the primary star (figure 9.25). As shown by these authors, the mass and orbital parameters of the secondary star play important roles in the radial mixing of protoplanetary objects and the delivery of water to the habitable zone of the primary star. The giant planet of the system also plays the important role of transferring angular momentum from the secondary star to the disk of protoplanets, and enhancing the radial mixing of these objects. As shown in figure 9.26, water delivery is less efficient in binaries with smaller perihelia since in such systems, the close approach of the binary companion to the giant planet increases its eccentricity, which in turn results in stronger interaction between this object and the disk of protoplanets, causing them to become unstable in very short time. The results of the simulations by Haghighipour & Raymond (2007) indicate that binary-planetary systems with giant planets at 5-10 AU, and binary perihelion distances of approximately 20 AU to 25 AU, will be more efficient in forming and hosting habitable planets (figure 9.27).
0.5 Future Prospects
The discovery of planets in binary star systems is one the interesting surprises of modern astronomy. Despite the long history of the study of planets in such environments, the recent detection of planets in moderately close binaries has confronted astronomers with many new challenges. Many aspects of the formation process of these planets are still unresolved, and questions regarding their frequency and detection techniques demand more detailed investigation.
The habitability of binary systems is also an open question. Although recent simulations of the late stage of terrestrial planet formation in binary systems have indicated that water-carrying planets can form in the habitable zone of a binary-planetary system, more studies are necessary to understand how protoplanetary objects can develop and evolve in such an environment. Such studies have implications for investigating the habitability of extrasolar planets, and tie directly to several of near-future NASA missions, in particular, the space mission Kepler.
- A binary-planetary system is a dual-star system that also hosts planetary bodies.
- As shown by Norwood Haghighipour (2002), the perturbative effect of the stellar companion on the dynamics of a planetary system around a star becomes important when the binary separation is smaller than 100 AU. At the present, approximately 25% of extrasolar planetary systems detected by radial velocity technique are in binary systems with separations ranging from 250 to 6000 AU.
- Symplectic integrators, as they were originally developed by Wisdom & Holman (1991), are not suitable for numerically integrating the orbits of small bodies in the gravitational fields of two massive objects. These integrators have been designed to integrate the orbits of planetary or smaller objects when they revolve around only one massive central body. Recently Chambers et al. (2002) have developed a version of a symplectic integrator that is capable of integrating the motion of a small object in the gravitational fields of two stellar bodies.
- In applying equation (9.1) to the stability of dust particles, one has to note that this equation does not take into account the effects of non-gravitational forces such as gas-drag or radiation pressure. The motion of a dust particles can be strongly altered by the effects of these forces.
- The influence zone of a planetary object with a mass around a star with a mass is defined as the region between and , where is the semimajor axis of the planet, and is its Hill radius
- Abt, H. A. 1979, The Frequencies of Binaries on the Main Sequence, Astron. J., 84, 1591
- Akeson, R. L., Koerner, D. W. & Jensen, E. L. N. 1998, A Circumstellar Dust Disk around T Tauri N: Subarcsecond Imaging at lambda =3 Millimeters, Astrophys. J., 505, 358
- Artymowicz, P. & Lubow, S. H. 1994, Dynamics of Binary-Disk Interaction. 1: Resonances and Disk Gap Sizes, Astrophys. J., 421, 651
- Barbieri, M., Marzari, F. & Scholl, H. 2002, Formation of Terrestrial Planets in Close Binary Systems: The Case of Alpha Centauri A, Astron. Astrophys., 396, 219
- Benest, D. 1988, Planetary Orbits in the Elliptic Restricted Problem. I-The Alpha Centauri System, Astron. Astrophys., 206, 143
- Benest, D. 1989, Planetary Orbits in the Elliptic Restricted Problem. II-The Sirius System, Astron. Astrophys., 223, 361
- Benest, D. 1993, Stable Planetary Orbits Around One Component in Nearby Binary Stars. II, Celest. Mech., 56, 45
- Benest, D. 1996, Planetary Orbits in the Elliptic Restricted Problem. III. The Coronae Borealis System, Astron. Astrophys., 314, 983
- Black, D. C. 1982, A Simple Criterion for Determining the Dynamical Stability of Three-Body Systems, Astron. J., 87, 1333
- Bodenheimer, P., Hubickyj, O. & Lissauer, J. J. 2000, Models of the in Situ Formation of Detected Extrasolar Giant Planets, Icarus, 143, 2
- Boss, A. P. 2000, Possible Rapid Gas Giant Planet Formation in the Solar Nebula and Other Protoplanetary Disks, Astrophys. J., 536, L101
- Boss, A. P. 2006, Gas Giant Protoplanets Formed by Disk Instability in Binary Star Systems, Astrophys. J., 641, 1148
- Broucke, R. A. 2001, Stable Orbits of Planets of a Binary Star System in the Three-Dimensional Restricted Problem, Celest. Mech. Dynamic. Astron., 81, 321
- Campbell, B., Walker, G. A. H. & Yang, S. 1988, A Search for Substellar Companions to Solar-Type Stars, Astrophys. J., 331, 902
- Chambers, J. E., Quintana, E. V., Duncan, M. J. & Lissauer, J. J., 2002, Symplectic Integrator Algorithms for Modeling Planetary Accretion in Binary Star Systems, Astrophys. J., 123, 2884
- Diakov, B. B. & Reznikov, B. I. 1980, Computer Simulation of Planet Formation in a Binary Star System Terrestrial Planets, Moon. Planet., 23, 429
- Drobyshevski, E. M. 1978, The Origin of the Solar System - Implications for Transneptunian Planets and the Nature of the Long-Period Comets, Moon. Planet., 18, 145
- Duquennoy, A. & Mayor, M. 1991, Multiplicity Among Solar-Type Stars in the Solar Neighborhood. II - Distribution of the Orbital Elements in an Unbiased Sample, Astron. Astrophys., 248, 485
- Dvorak, R. 1982, Österr. Akad. d. Wiss., Math-Nat. Klasse n 191, 423
- Dvorak, R. 1984, Numerical Experiments on Planetary Orbits in Double Stars, Celest. Mech., 34, 369
- Dvorak, R. 1986, Critical Orbits in the Elliptic Restricted Three-Body Problem, Astron. Astrophys., 167, 379
- Dvorak, R., Froeschlé, Ch. & Froeschlé, Ci. 1989, Stability of Outer Planetary Orbits (P-Types) in Binaries, Astron. Astrophys., 226, 335
- Dvorak, R., Pilat-Lohinger, E., Funk, B. & Freistetter, F. 2003, Planets in Habitable Zones:. A Study of the Binary Gamma Cephei, Astron. Astrophys., 398, L1
- Dvorak, R., Pilat-Lohinger, E., Bois, E., Funk, B., Freistetter, F. & Kiseleva-Eggleton, L. 2004, Planets in Double Stars: the Cephei System, RevMexAA (Series de Conferencias), 21, 222
- Fatuzzo, M., Adams, F. C., Gauvin, R. & Proszkow, E. M. 2006, A Statistical Stability Analysis of Earth-like Planetary Orbits in Binary Systems, Pub. Astron. Soc. Pacific., 118, 1510
- Forget, F. & Pierrehumbert, R. T. 1997, Warming Early Mars with Carbon Dioxide Clouds That Scatter Infrared Radiation, Science, 278, 1273
- Fuhrmann, K. 2004, Astron. Nachr., 325, 3
- Graziani, F. & Black, D. C. 1981, Orbital Stability Constraints on the Nature of Planetary Systems, Astrophys. J., 251, 337
- Griffin, R. F., Carquillat, J. M. & Ginestet, N. 2002, Spectroscopic Binary Orbits from Photoelectric Radial Velocities, Paper 163: HD 213503/4 and HD 220636/7 with a Conjecture Concerning Cephei, The Observatory, 122, 90
- Haghighipour, N. 2004, On the Dynamical Stability of Cephei; An S-type Binary Planetary System, In: AIP Conference Proceedings 713: The Search for Other Worlds, ed by S. S. Holt, & D. Deming, Melville, New York, pp 269-272
- Haghighipour, N. 2005, Dynamical Stability and Habitability of Cephei Binary-Planetary System, astro-ph/0509659
- Haghighipour, N. 2006, Dynamical Stability and Habitability of the Cephei Binary-Planetary System, Astrophys. J., 644, 543
- Haghighipour, N. & Raymond, S. N. 2007, Habitable Planet Formation in Binary-Planetary Systems, Astrophys. J., 666, 436
- Hagle, J. & Dvorak, R. 1988, An Analytical Study of Stable Planetary Orbits in the Circular Restricted Problem, Celest. Mech., 42, 355
- Harrington, R. S. 1977, Planetary Orbits in Binary Stars, Astron. J., 82, 753
- Hatzes, A. P., Cochran, W. D., Endl, M., McArthur, B., Paulson, D. B., Walker, G. A. H., Campbell, B. & Yang, S. 2003, A Planetary Companion to Cephei A, Astrophys. J., 599, 1383
- Hayashi, C. 1981, Structure of the Solar Nebula, Growth and Decay of Magnetic Fields and Effects of Magnetic and Turbulent Viscosities on the Nebula, Prog. Theor. Phys. Suppl., 70, 35
- Hénon, M. & Guyot, M. 1970, Stability of Periodic Orbits in the Restricted Problem, In: Periodic Orbits, Stability and Resonances, ed by G. E. O. Giacaglia, D. Reidel Pub. Comp., Netherlands, pp 349
- Heppenheimer, T. A. 1974, Outline of a Theory of Planet Formation in Binary Systems, Icarus, 22, 436
- Heppenheimer, T. A. 1974, On the Formation of Planets in Binary Star Systems, Astron. Astrophys., 65, 421
- Holman, M. J. & Wiegert, P. A. 1999, Long-Term Stability of Planets in Binary Systems, Astron. J., 117, 621
- Innanen, K. A., Zheng, J. Q., Mikkola, S. & Valtonen, M. J. 1997, The Kozai Mechanism and the Stability of Planetary Orbits in Binary Star Systems, Astron. J., 113, 1915
- Jones, B. W., Underwood, D. R. & Sleep, P.N. 2005, Prospects for Habitable “Earths” in Known Exoplanetary Systems, Astrophys. J., 622, 1091
- Jones, B. W., Sleep, P.N. & Underwood, D. R. 2006, Habitability of Known Exoplanetary Systems Based on Measured Stellar Properties, Astrophys. J., 649, 1010
- Kasting, J. F., Whitmire, D. P. & Reynolds, R. T. 1993, Habitable Zones Around Main Sequence Stars, Icarus, 101, 108
- Kortenkamp, S. J., Wetherill, G. W. & Inaba, S. 2001, Runaway Growth of Planetary Embryos Facilitated by Massive Bodies in a Protoplanetary Disk, Science, 293, 1127
- Kozai, Y. 1962, Secular Perturbations of Asteroids with High Inclination and Eccentricity, Astron. J., 67, 591
- Krist, J. E., Stapelfeldt, K. R., Golimowski, D. A., Ardila, D. R., Clampin, M., Martel, A. R., Ford, H. C., Illingworth, G. D. & Hartig, G. F. 2005, Hubble Space Telescope ACS Images of the GG Tauri Circumbinary Disk, Astron. J., 130, 2778
- Kubala, A., Black, D. C. & Szebehely, V. 1993, Stability of Outer Planetary Orbits Around Binary Stars - A Comparison of Hill’s and Laplace’s Stability Criteria, Celest. Mech. Dynamic. Astron., 56, 51
- Lissauer, J. J., Quintana, E. V., Chambers, J. E., Duncan, M. J. & Adams, F. C. 2004, Terrestrial Planet Formation in Binary Star Systems, RevMexAA (Series de Conferencias), 22, 99
- Malmberg, D., Davies, M. B. & Chambers, J. E. 2007, The Instability of Planetary Systems in Binaries: How the Kozai Mechanism Leads to Strong Planet-Planet Interactions, Mon. Not. Roy. Ast. Soc., 377, L1
- Marzari, F., Scholl, H., Tomasella, L. & Vanzani, V. 1997, Gas Drag Effects on Planetesimals in the 2:1 Resonance With Proto-Jupiter, Planet. Space. Sci., 45, 337
- Marzari, F. & Scholl, H. 2000, Planetesimal Accretion in Binary Star Systems, Astrophys. J., 543, 328
- Mathieu, R. D. 1994, Pre-Main-Sequence Binary Stars, Ann. Rev. Astron. Astrophys. 32, 465
- Mathieu, R. D., Ghez, A. M., Jensen, E. L. & Simon, M. 2000, Young Binary Stars and Associated Disks. In: Protostars and Planets IV, ed by V. Mannings, A. P. Boss & S. S. Russell, Univ. Arizona Press, Tucson, pp 703
- Mayer, L., Boss, A. P. & Nelson, A. F. 2007, Gravitational Instability in Binary Protoplanetary Disks, arXiv0705.3182M
- Mayer, L., Wadsley, J., Quinn, T. & Stadel, J. 2005, Gravitational Instability in Binary Protoplanetary Discs: New Constraints on Giant Planet Formation, Mon. Not. Roy. Ast. Soc., 363, 641
- Menou, K. & Tabachnik, S. 2003, Dynamical Habitability of Known Extrasolar Planetary Systems, Astrophys. J., 583, 473
- Mischna, M. A., Kasting, J. F., Pavlov, A. & Freedman, R. 2000, Influence of Carbon Dioxide Clouds on Early Martian Climate, Icarus, 145, 546
- Morbidelli, A., Chambers, J., Lunine, J. I., Petit, J. M., Robert, F., Valsecchi, G. B. & Cyr, K. E. 2000, Source Regions and Time Scales for the Delivery of Water to Earth, Meteorit. Planet. Sci., 35, 1309
- Moriwaki, K. & Nakagawa, Y. 2004, A Planetesimal Accretion Zone in a Circumbinary Disk, Astrophys. J., 609, 1065
- Musielak, Z. E., Cuntz, M., Marshall, E. A. & Stuit, T. D. 2005, Stability of Planetary Orbits in Binary Systems, Astron. Astrophys., 434, 355
- Nelson, A. F. 2000, Planet Formation is Unlikely in Equal-Mass Binary Systems with a 50 AU, Astrophys. J., 537, L65
- Neuhäuser, R., Mugrauer, M., Fukagawa, M., Torres, G. & Schmidt, T. 2007, Direct Detection of Exoplanet Host Star Companion Cep B and Revised Masses for Both Stars and the Sub-stellar Object, Astron. Astrophys., 462, 777
- Norwood, J. W. & Haghighipour, N. 2002, On the Stability of And. Extrasolar Planetary System: An S-Type Binary-Planetary System with more Than One Planet, Bul. Am. Ast. Soc., 34, 892
- Pendleton, Y. J. & Black, D. C. 1983, Further Studies on Criteria for the Onset of Dynamical Instability in General Three-Body Systems, Astron. J., 88, 1415
- Pichardo, B., Sparke, L. S. & Aguilar, L. A. 2005, Circumstellar and Circumbinary Discs in Eccentric Stellar Binaries, Mon. Not. Roy. Ast. Soc., 359, 521
- Pilat-Lohinger, E. & Dvorak, R. 2002, Stability of S-type Orbits in Binaries, Celest. Mech. Dynamic. Astron., 82, 143
- Pilat-Lohinger, E., Funk, B. & Dvorak, R. 2003, Stability Limits in Double Stars. A Study of Inclined Planetary Orbits, Astron. Astrophys., 400, 1085
- Pilat-Lohinger, E., Dvorak, R., Bois, E. & Funk, B. 2004, Stable Planetary Motion in Double Stars In: APS Conference Series 321: Extrasolar Planets: Today and Tomorrow, ed by J. P. Beaulieu, A. Lecavelier des Etangs, & C. Terquem, pp 410
- Queloz, D., Mayor, M., Weber, L., Blécha, A., Burnet, M., Confino, B., Naef, D., Pepe, F., Santos, N., & Udry, S. 2000, The CORALIE Survey for Southern Extrasolar Planets. I. A Planet Orbiting the Star Gliese 86, Astron. Astrophys., 354, 99
- Quintana, E. V., Lissauer, J. J., Chambers, J. E. & Duncan, M. J. 2002, Terrestrial Planet Formation in the Alpha Centauri System, Astrophys. J., 576, 982
- Quintana, E. V. & Lissauer, J. J. 2006, Terrestrial Planet Formation Surrounding Close Binary Stars, Icarus, 185, 1
- Quintana, E. V., Adams, F. C., Lissauer, J. J. & Chambers, J. E. 2007, Terrestrial Planet Formation Around Individual Stars Within Binary Star Systems, Astrophys. J., 660, 807
- Rabl, G. & Dvorak, R. 1988, Satellite-Type Planetary Orbits in Double Stars - A Numerical Approach, Astron. Astrophys., 191, 385
- Rodriguez, L. F., D’Alessio, P., Wilner, D. J., Ho, P. T. P., Torrelles, J. M., Curiel, S., Gomez, Y., Lizano, S., Pedlar, A., Canto, J. & Raga, A. C. 1998, Compact Protoplanetary Disks Around the Stars of a Young Binary System, Nature, 395, 355
- Silbert, J., Gledhill, T., Duchéne, G. & Ménard, F. 2000, Near-Infrared Imaging Polarimetry of the GG Tauri Circumbinary Ring, Astrophys. J., 536, L89
- Szebehely, V. 1980, Stability of Planetary Orbits in Binary Systems, Celest. Mech., 22, 7
- Szebehely, V. & McKenzie, R. 1981, Stability of Outer Planetary Systems, Celest. Mech., 23, 3
- Szebehely, V. 1984, Review of Concepts of Stability, Celest. Mech., 34, 49
- Takeda, G. & Rasio, F. A. 2006, Eccentricities of Planets in Binary Systems, Astrophy. Sp. Sci., 304, 239
- Thébault, P., Marzari, F., Scholl, H., Turrini, D. & Barbieri, M. Planetary Formation in the Cephei System, Astron. Astrophys., 427, 1097
- Thébault, P., Marzari, F. & Scholl, H. 2006, Relative Velocities Among Accreting Planetesimals in Binary Systems: The Circumprimary Case, Icarus, 183, 193
- Trilling, D. E., Stansberry, J. A., Stapelfeldt, K. R., Rieke, G. H., Su, K. Y. L., Gray, R. O., Corbally, C. J., Bryden, G., Chen, C. H., Boden, A. & Beichman, C. A. 2007, Debris Disks in Main-Sequence Binary Systems, Astrophys. J., 658, 1264
- Torres, G. 2007, The Planet Host Star Cephei: Physical Properties, the Binary Orbit, and the Mass of the Substellar Companion, Astrophys. J., 654, 1095
- Verrier, P. E. & Evans, N. W. 2006, Planets and Asteroids in the Cephei System, Mon. Not. Roy. Ast. Soc., 368, 1599
- Walker, G. A. H., Bohlender, D. A., Walker, A. R., Irwin, A. W., Yang, S. L. S. & Larson, A. 1992, Gamma Cephei - Rotation or Planetary Companion? Astrophys. J., 396, L91
- Weidenschilling, S. J. 1977, The Distribution of Mass in the Planetary System and Solar Nebula, Astrophys. Sp. Sci., 51, 153
- White, R. J., Ghez, A. M., Reid, I. N. & Schultz, G. 1999, A Test of Pre-Main-Sequence Evolutionary Models across the Stellar/Substellar Boundary Based on Spectra of the Young Quadruple GG Tauri, Astrophys. J., 520, 811
- White, R. J. & Ghez, A. M. 2001, Observational Constraints on the Formation and Evolution of Binary Stars, Astrophys. J., 556, 265
- Whitmire, D. P., Matese, J. L., Criswell, L. & Mikkola, S. 1998, Habitable Planet Formation in Binary Star Systems, Icarus, 132, 196
- Wiegert, P. A. & Holman, M. J. 1997, The Stability of Planets in the Alpha Centauri System, Astrophys. J., 113, 1445
- Williams, D. M. & Kasting, J. F. 1997, Habitable Planets with High Obliquities, Icarus, 129, 254
- Wisdom, J. & Holman, M. J. 1991, Symplectic Maps for the N-body Problem, Astrophys. J., 102, 1528
- Ziglin, S. L. 1975, Secular Evolution of the Orbit of a Planet in a Binary-Star System, Sov. Astron. Let., 1, 194