ON THE POSSIBILITY OF ADDITIONAL PLANETS IN THE γ CEPHEI BINARY-PLANETARY SYSTEM Joseph Castro Mentor: Nader Haghighipour

ABSTRACT Results of the simulations of the dynamical stability of additional hypothetical planets in the γ Cephei binary-planetary system are presented. Using the reported orbital parameters of the existing bodies (Hatzes et al. 2003), simulations were carried out for different values of the , eccentricity, and inclination of an additional Jupiter-like plantet to test the stability of this object at a distance of 3 to 5 AU from the primary . Simulations were also carried out for fictitious Saturn- and Neptune-like planets. Results indicate that the system appears to be dynamically full and no stable of any additional planets were found.

1. INTRODUCTION In 1992, the first discovery of an extra-solar planet was made around the pulsar PSR 1257 (Wolszczan and Frail 1992). It was not until three later on October 6, 1995, that an extra-solar planet was discovered orbiting a -like star (Mayor and Queloz 1995). Since then, astronomers have been finding extra-solar planets at a maddening rate, and the tally is now up to over 190. Among these planets, approximately 25% are members of binary or multistar systems. These systems, for the most part, have very large stellar separations, generally between 250 and 6500 AU; distances at which the gravitational effect of the farther companion is insignificant. When the separation is less than 100 AU, however, the farther companion can have profound effects on the dynamics of the system. Among the three known systems of this type, this paper focuses on the binary-planetary system of γ Cephei. Gamma Cephei is a spectroscopic binary located approximately 50 light-years from Earth. This system has a 1.59 solar-mass K1 IV as its primary star, and a 0.34-0.92 solar-mass as its secondary. The semimajor axis and eccentricity of this binary are 18.5 ± 1.1 AU and 0.361 ± 0.023, respectively (Hatzes et al. 2003). A Jupiter- like planet with at least 1.7 Jupiter-, orbits the system’s primary star in an with a semimajor axis of 2.13 ± 0.05 AU and an eccentricity of 0.12 ± 0.05 (Hatzes et al. 2003). Table 1 shows the six orbital parameters of the γ Cephei system as reported by Hatzes et al. (2003). Gamma Cephei is an interesting case in that, it is not only the first of three moderately close planet-hosting binary/multistar systems, it also could have been the very first extra-solar planetary system ever discovered. In 1988, a Canadian team led by Bruce Campbell and Gordon Walker, suspected the existence of a third body in the binary system of γ Cephei while studying the variations of this system. In 1992, however, the possibility that the RV variations observed by the authors were due to the long-period RV variability of most K giants (under which γ Cephei is classified) was brought about, and their claim of the presence of an extra-solar planet was retracted. Finally, ten years later, using new radial velocity measurements made by the McDonald Obseravtory, the existence of the planet was confirmed Studies have already been carried out on the stability of the planet of γ Cephei system. In 2003, a study by Dvorak et al. (2003) confirmed the dynamical stability of this planet for a period of approximately 108 years. The stability of this object was also confirmed separately by Haghighipour (2006). Using a range of values for the eccentricity and semimajor axis of the binary as well as various inclinations of the Jupiter-like planet, Haghighipour (2006) found that the orbit of this planet is stable when the value of the binary eccentricity is smaller than 0.5, and the inclination of the planet is between 0◦ and 40◦. Studies have also been carried out to determine whether additional planets and asteroids can exist within the γ Cephei system. Verrier and Evans (2006) performed numerical simulations to study the stability of test particles in this system and were able to identify possible locations where additional planets and asteroids could be stable. These authors also studied the stability of objects around each star (interior and exterior to the planet) as well as the circumbinary case. Around the primary star, these authors identified three zones of stability for at least 1-Myr, from 0.5 to 0.75 AU, 1.05 to 1.15 AU, and also 1.2 to 1.3 AU. Around the secondary, they identified a region of stability at approximately 0.5 to 1.5 AU from this star. Finally, in a region around both stars, they determined a zone of stability

TABLE 1 Orbital Parameters of the γ Cephei System (Hatzes et al. 2003)

Parameter Planet Binary a (AU) 2.13 ± 0.05 18.5 ± 1.1 e 0.12 ± 0.05 0.361 ± 0.023 i (deg) 0 0 Ω (deg) 0 0 ω (deg) 49.6 ± 25.6 158.76 ± 1.2 M (deg) 104 353 2

TABLE 2

Parameter Known Planet Extra Planet Binary

Mass 1.7 MJ 0.25-2.5 MJ 0.4 M⊙ a (AU) 2.13 3-5 18.5 e 0.12 0.0-0.2 0.361 i (deg) 0 0-10 0 Ω (deg) 0 0 0 ω (deg) 49.6 0 158.76 M (deg) 104 0 353 extending out from near 65 AU. Instead of using test particles, Haghighipour (2006) simulated the dynamics of an actual Earth-like planet. Con- sidering a range of 0.3 to 4.0 AU for the semimajor axis of this object, and several different values of its inclination, this author was able to show that a region of stability only exists at distances of 0.3 to 0.8 AU from the system’s primary star. These results slightly overlap with those presented by Verrier and Evans (2006). Haghighipour (2006) also showed that the habitable zone of the primary at 3.05 to 3.7 AU is unstable. The focus of this paper is on the existence of possible additional planets in the γ Cephei system. The objects of interest are those with masses equal to those of Jupiter, Saturn and Neptune, on orbits beyond the location of the actual planet of the system.

2. STABILITY OF AN ADDITIONAL JUPITER-LIKE PLANET In order to find out whether the γ Cephei system can harbor an additional large planet, one has to study the dynamical stability of such an object at different distances from the primary star. However, it is important to choose the initial value of the semimajor axis of this planet to be outside the Hill sphere of other objects. The Hill sphere is the sphere of gravitational influence of an astronomical body in the face of the perturbations of the body it orbits around. One can calculate the Hill radius (radius of the Hill sphere) of a body from equation

m 1/3 RH = a( ) . (1) 3M∗ In this equation, a is the semimajor axis of the object, m is its mass, and M∗ is the mass of the central body. In the case of γ Cephei, the Hill radius of the known Jupiter-like planet is approximately at 0.15 AU, and the Hill radius of the binary companion is at 8.1 AU (Verrier and Evans 2006). Taking into account the uncertainty on the semimajor axis of the known planet, and the furthest distance it can orbit due to its eccentricity, any additional planets must be placed at least 2.65 AU from the primary star. As for the binary, one can calculate its closest approach via the formula r = ab(1 − eb). (2)

Given that ab=18.5 AU and eb=0.361, the value of r is approximately 12 AU. This means that at certain points in its orbit, the binary’s effect can stretch to less than 4 AU from the primary star. Any additional planets should technically not be placed any further out than 4 AU, but since the binary can influence at this distance only every 57 years (one period of its orbit), there is some leeway. Numerical integrations were carried out using the hybrid symplectic and also the Bulirsch-Stoer integrators of the MERCURY 6.2 Integrator Package (Chambers 1999). The timestep of each simulation was ten days. The orbital parameters of the known planet of the system and those of the binary were taken from Table 1. The initial longitude of ascending node (Ω), argument of pericentre (ω), and mean anomaly (M) of the additional planet were set to zero and the semimajor axis, eccentricity, mass and inclination of this planet were varied. Table 2 shows an overview of these values. Figure 1 shows the time-variations of the semimajor axis, eccentricity, and inclination of the three bodies of the system for one of our simulations. As shown here, the system becomes unstable after approximately 2 million years. It is important to note the jump in semimajor axis, eccentricity, and inclination of the binary and Jupiter-like planet after the additional planet becomes unstable and leaves the system. To conserve angular momentum, when the additional object leaves the system, the other two bodies move to different orbits. A more pronounced representation of this jump can be seen in Figure 2. The instability of the additional planet was observed in all our simulations. Figure 3 shows the results of numerical integrations performed for a Jupiter-mass planet with 0 eccentricity and no inclination. As shown here, the longest lifetime of the Jupiter-mass planet is no more than 150,000 years, and corresponds to a distance of 3.5 AU. This increased ’stability’ at the distance of 3.5 AU from the primary star was also reported by Verrier and Evans (2006). We also carried out a few simulations at 2.75 AU, 5.25 AU and 5.5 AU. As expected, results indicate that the additional planet became unstable within a few tens of thousands of years. It is expected that decreasing or increasing the semimajor axes will futher cut short the lifetime of the additional planet due to its venturing into the Hill sphere of the other planet or the binary. In addition to investigating the dynamical stability of the additional planet at various semimajor axes, the stability of the system for various initial eccentricities of the additional planet was also examined. The value of the planet’s 3

20 15 Binary 10 Known Planet Additional Planet 5 0

Semimajor Axis (AU) 1 0.8 0.6 0.4 0.2 Eccentricity 0

10

5

0 Inclination (Deg.) 0 2 4 6 8 10 Millions of Years Fig. 1.— Plots of semimajor axis, eccentricity, and inclination versus time for a typical simulation. The initial semimajor axis of the additional planet is 3.5AU, its initial eccentricity is 0, and its initial inclination is 10◦.

0.37 0.365 0.36 0.355 0.35 0.345 0.34 18.1 18.2 18.3 18.4 18.5 18.6 1 0.8 0.6 0.4 0.2 Additional Planet Eccentricity 0 2 3 4 5 6 7 0.3 0.25 0.2

0.15 Known Planet 0.1 2 2.05 2.1 2.15 2.2 2.25 2.3 Semimajor Axis (AU) Fig. 2.— Eccentricity versus semimajor axis for the system of Figure 1. 4

1e+08

1e+07

1e+06

1e+05

1e+04

1e+03 Lifetime (Years) 1e+02

1e+01

1e+00 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 Initial semimajor axis (AU) Fig. 3.— Lifetime of additional Jupiter-mass planet with 0 eccentricity and 0◦ inclination. Black indicates the additoinal planet, red indicates the known Jupiter-sized object, and green indicates the single case where the known planet left the system first.

1e+07

1e+06 3.5 AU 3.75 AU 4.0 AU

1e+05

1e+04

1e+03

Lifetime (Years) 1e+02

1e+01

1e+00 0 0.05 0.1 0.15 0.2 0 0.05 0.1 0.15 0.2 0 0.05 0.1 0.15 0.2 Initial eccentricity Fig. 4.— Stability of the system as the initial eccentricity of the additional planet vaires. Other than 0 eccentricity, no other eccentricity appears to be preferred over another for the entire range of simulations. The initial semimajor axis of the planet was 3.5 AU in all cases shown here. eccentricity were chosen to be 0, 0.05, 0.1, 0.15, 0.2. Simulations were carried out for each value of the semimajor axis with the aforementioned eccentricities. Figure 4 shows the lifetime of the additional planet in terms of its initial eccentricity at three semimajor axes of 3.5 AU, 3.75 AU, and 4.0 AU. As shown in this figure, can play an important role in the stablitiy a planet. For most of our cases, an eccentricity of zero seemed to be preferred, but sometimes an eccentricity of 0.05 could increase the additional planet’s lifetime by a considerable amount. Some of our simulations, however, resulted in an increased lifetime by increasing the eccentricity beyond 0.05. In general, it appeared that the additional planet was least stable at 0.2 eccentrcity, leaving the system within a few thousand years. It is expected that increasing the eccentricity further will cause the system to become unstable. We also simulated the dynamics of the system with and additional planet with a mass of 0.25 to 2.5 MJ . The semimajor axis and eccentricity of this object were varied for each mass. Initially, it appeared that the lifetime of the additional planet increased as the mass decreased. Figure 5, however, shows that there is no apparent trend to the lifetime of this object. This figure considers the stability of various masses, in cirular orbits, and with an initial semimajor axes of 3.5 AU. In addition to varying the semimajor axis, eccentricity, and mass of the additional Jupiter-like planet, the planet’s inclination was also varied from 0◦ to 10◦, with a 2◦ increment. In these simulations, no stable orbits were found beyond 2 million years. Figure 6 shows the results of some simulations for a few masses in circular orbits, initially at 3.5 AU. Here, it seems, no inclination is preferred over another. The lifetime of the additional planet does not increase or decrease as a function of the initial inclination. Therefore, it is expected that inclinations beyond 10◦ will not produce better results. Since a 0.25 Jupiter-mass planet appeared to have the greatest probability of being stable, 5

1e+08

1e+07

1e+06

1e+05

1e+04

1e+03 Lifetime (Years) 1e+02

1e+01

1e+00 0 0.25 0.5 0.75 1 1.25 1.5 1.75 2 2.25 2.5 2.75 Mass of additional planet (Jupiter-masses) Fig. 5.— Additional planet lifetimes when the mass was varied. The planet was given an initial semimajor axis of 3.5 AU and a circular orbit. The longest lifetime occured for a 0.25 MJ planet, and the shorest lifetime occurued for a 1.75 MJ planet.

1e+07 0.5 Jup. Mass 1.5 Jup. Mass 2.5 Jup. Mass

1e+06

1e+05

1e+04

1e+03

Lifetime (Years) 1e+02

1e+01

1e+00 0 2 4 6 8 10 0 2 4 6 8 10 0 2 4 6 8 10 Inclination (degrees) Fig. 6.— Additional planet lifetimes when the inclination was varied for three different masses.

TABLE 3

Parameter Known Planet Extra Planet Binary

Mass 1.7 MJ 0.5-1.0 MS 0.4 M⊙ a (AU) 2.13 3-5 18.5 e 0.12 0.0-0.2 0.361 i (deg) 0 0-10 0 Ω (deg) 0 0 0 ω (deg) 49.6 0 158.76 M (deg) 104 0 353 more inclinations were tried for that case. The results indicated that the system was unstable.

3. STABILITY OF AN ADDITIONAL SATURN-LIKE PLANET Simulations were also carried out for an additional Saturn-like planet in the γ Cephei system. Table 3 shows the orbital parameters used in these simulations. Numerical simulations were performed for a Saturn-mass object as well as a 0.5 Saturn-mass planet. The initial semimajor axis of each planet was varied between 3 and 5 AU, with an eccentricity ranging from 0.0 to 0.2. The inclination was not varied in these simulations, and was kept constant at 0◦. The results mirror those of the Jupiter-like planet. That is, no stable orbits were found and the longest lifetimes 6

1e+08

1e+07

1e+06

1e+05

1e+04

1e+03 Lifetime (Years) 1e+02

1e+01

1e+00 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 Initial semimajor axis (AU) Fig. 7.— Lifetime of an additional Saturn-mass planet with 0 eccentricity and 0◦ inclination. Black indicates the additoinal planet, and red indicates the known Jupiter-sized object. The longest Saturn lifetime was found at 3.5 AU.

TABLE 4

Parameter Known Planet Extra Planet Binary

Mass 1.7 MJ 1.0-1.5 MN 0.4 M⊙ a (AU) 2.13 3-5 18.5 e 0.12 0.0-0.2 0.361 i (deg) 0 0-10 0 Ω (deg) 0 0 0 ω (deg) 49.6 0 158.76 M (deg) 104 0 353 were at 3.5 AU from the primary star. Figure 7 shows how the planet’s stability varies with semimajor axis.

4. STABILITY OF AN ADDITIONAL NEPTUNE-LIKE PLANET Neptune-like planets were also placed in the γ Cephei system. Numerical integrations were carried out for a 1.0 and 1.5 Neptune-mass object. The object was placed between 3 to 5 AU and was given an eccentricity between 0.0 and 0.2 (Table 4). Similar to the simulations of the Saturn-like planet, no inclinations other than 0◦ were used. Results from these integrations were very similar to those seen in the case of the Saturn-like planet. The only difference is that the Neptune-like planet actually appeared to be less stable than the Saturn-like object. Figures 8 and 9 show plots of the semimajor axes and eccentricities of γ Cephei system when a Neptune-mass planet is present. Note the periodicity in the graphs of the known planet of the system (Figure 8). This is due to the closest approach between the binary and the additional object.

5. CONCLUSION We presented the results of over 200 numerical integrations of the three body system of γ Cephei, assuming additional Jupiter to Neptune-sized planets. Simulations indicated that no such large bodies can exist in the system for a period larger than a few million years. Results also indicate that a region of increased lifetime exists at 3.5 AU from the primary star. The first part of this study focused on the possibility of the presence of additional Jupiter-like planets in the γ Cephei binary-planetary system. Previous studies report the possbility of terrestrial planets close to the primary star, as well as possible planets/asteroids orbiting the secondary star. However, this study, focused on the possibility of an extra object in the region between 3 to 5 AU from the primary star. Simulations were run for ficticious planets of masses of 0.25 to 2.5 Jupiter-masses, with eccentricities ranging from 0.0 to 2.0 and inclinations 0◦ to 10◦. Numerical integrations resulted in no long-term stable orbits. Also, due to the complex dynamics of the system, no stability trends were found through mass, eccentricity, or inclination variations. After all integrations were complete, the longest stable orbit was found to be approximately 2,400,000 years corresponding to a 0.25 MJ planet initially at 3.5 AU and with an initial eccentricity of zero and inclination of 10◦. Similar results followed for the study of possible Saturn- and Neptune-like planets in the system. No long-term stable orbits were found. Comparing the results from these two studies show that the Saturn-like object may have a lifetime greater than that of the Neptune-like body. Even so, it is not expected that smaller bodies, such as Uranus-like 7

18.6 Binary 18.5 18.4 18.3 18.2 0 1 2 3 4 5 6 7 8 9 10

4.2 4.0 Additional Planet 3.8 3.6 3.4 0 1 2 3 4 5 6 7 8 9 10 Semimajor Axis (AU)

2.13 Known Planet 2.11 2.09 2.07 0 1 2 3 4 5 6 7 8 9 10 1E5 Years Fig. 8.— Semimajor axis vs time for the binary, Jupiter-like planet, and a Neptune-like object. The initial semimajor axis of the Neptune-like planet was 3.5 AU, with an eccentricity of 0 and 0◦ inclination.

0.37 0.365 Binary 0.36 0.355 0.35 0.345 0.34 0 1 2 3 4 5 6 7 8 9 10

0.6 0.5 Additional Planet 0.4 0.3 0.2 0.1 Eccentricity 0 0 1 2 3 4 5 6 7 8 9 10

0.3 0.25 Known Planet 0.2 0.15 0.1 0 1 2 3 4 5 6 7 8 9 10 1E5 Years Fig. 9.— Eccentricity vs time for the binary, Jupiter-like planet, and Neptune-like planet. The initial semimajor axis of the Neptune-like planet was 3.5 AU, with an eccentricity of 0 and 0◦ inclination. 8 planets, would have even less stable orbits.

REFERENCES Chambers, J. E. 1999, MNRAS, 304, 793. Rivera, E., and Haghighipour, N. 2006, MNRAS submitted (astro- Dvorak, R. Pilar-Lohinger, E., Funk, B. and, Freistetter, F. 2003, ph/0406429). A&A 398, L1-L4. Verrier, P. E., and Evans, N. W. 2006, MNRAS 368, 1599. Haghighipour, N. 2006, APJ, 644, 543. Wolszczan, A., and Frail, D. 1992, Nature, 355, 145. Hatzes, A. Cochran, W., Endl, M., McArthur, B., Paulson, D., Walker, G., Campbell, B. and Yang, S. 2003, APJ, 599, 1383. Mayor, M., and Queloz, D. 1995, Nature 378, 355.