On the Origin of the Wide-Orbit Circumbinary Giant Planet HD 106906 a Dynamical Scenario and Its Impact on the Disk
Total Page:16
File Type:pdf, Size:1020Kb
Astronomy & Astrophysics manuscript no. Article c ESO 2018 September 11, 2018 On the origin of the wide-orbit circumbinary giant planet HD 106906 A dynamical scenario and its impact on the disk L. Rodet1, H. Beust1, M. Bonnefoy1, A.-M. Lagrange1, P. A. B. Galli1; 2, C. Ducourant3, and R. Teixeira2 1 Univ. Grenoble Alpes, CNRS, IPAG, F-38000 Grenoble, France e-mail: [email protected] 2 Instituto de Astronomia, Geofísica e Ciências Atmosféricas, Universidade de São Paulo, Rua do Matão, 1226, Cidade Universitária, 05508-900, São Paulo - SP, Brazil 3 Laboratoire d’Astrophysique de Bordeaux, Univ. Bordeaux, CNRS, B18N, allée Geoffroy Saint-Hilaire, 33615 Pessac, France Received 16 December 2016 / Accepted 4 March 2017 ABSTRACT Context. A giant planet has been recently resolved at a projected distance of 730 au from the tight pair of young (∼ 13 Myr) intermediate-mass stars HD 106906AB in the Lower Centaurus Crux (LCC) group. The stars are surrounded by a debris disk which displays a ring-like morphology and strong asymmetries at multiple scales. Aims. We aim to study the likelihood of a scenario where the planet formed closer to the stars in the disk, underwent inward disk- induced migration, and got scattered away by the binary star before being stabilized by a close encounter (fly-by). Methods. We performed semi-analytical calculations and numerical simulations (Swift_HJS package) to model the interactions be- tween the planet and the two stars. We accounted for the migration as a simple force. We studied the LCC kinematics to set constraints on the local density of stars, and therefore on the fly-by likelihood. We performed N-body simulations to determine the effects of the planet trajectories (ejection and secular effects) onto the disk morphology. Results. The combination of the migration and mean-motion resonances with the binary star (often 1:6) can eject the planet. Nonethe- less, we estimate that the fly-by hypothesis decreases the scenario probability to less than 10−7 for a derived local density of stars of 0.11 stars/pc3. We show that the concomitant effect of the planet and stars trajectories induce spiral-features in the disk which may correspond to the observed asymmetries. Moreover, the present disk shape suggests that the planet is on an eccentric orbit. Conclusions. The scenario we explored is a natural hypothesis if the planet formed within a disk. Conversely, its low probability of occurrence and the fact that HD 106906 b shares some characteristics with other systems in Sco-Cen (e.g., HIP 78530, in terms of mass ratio and separation) may indicate an alternative formation pathway for those objects. Key words. celestial mechanics – methods: numerical – planet and satellites: dynamical evolution and stability – stars: planetary systems – planet-disk interactions 1. Introduction spread of about 10 Myr (Pecaut et al. 2012). In recent years, high contrast imaging has revealed the circumstellar environment of More than 3500 exoplanets have been found in the last three 1 HD 106906AB: an 11 ± 2 MJ planet located at 732 ± 30 au in decades , but few among them have been detected to be hun- projected separation (Bailey et al. 2013) and an asymmetric de- dreds of astronomical units (au) from their star. As the devel- bris disk nearly viewed edge-on, imaged by SPHERE (Lagrange opment of direct imaging reveals more of those wide planetary- et al. 2016a), GPI and HST (Kalas et al. 2015) and MagAO (Wu mass companions, classical theories of planet formation fail at et al. 2016). More recently, the binary nature of HD 106906 was explaining their origin. In the two scenarios, core accretion (Pol- inferred thanks to observations with the instruments HARPS and lack et al. 1996) and gravitational instability (Boss 1997), the PIONIER (Lagrange et al. 2016b). It turns out to be a 13±2 Myr planets form within the primordial gas disk. However, the limited old SB2 binary consisting of two F5 V-type stars with very sim- extent of the disk (see eg Fig. 5 in Lieman-Sifry et al. 2016) does ilar masses. Table 1 summarizes the key characteristics of the not enable the formation of a giant planet far away from its star. system components. No further information is known about the arXiv:1703.01857v2 [astro-ph.EP] 23 Mar 2017 Thus, when the star around which orbits the very wide and mas- orbit of the planet, which must have an orbital period of at least sive HD 106906AB b turned out to be a binary star (Lagrange 3000 years. The binary orbit is also not much constrained yet, but et al. 2016b), it has been suggested that dynamical interactions given its short orbital period (< 100 days), it will presumably be could account for the current position of the planet (Lagrange et better known in the near future. al. 2016b, Wu et al. 2016). The planet HD 106906 (or also HIP 59960) is located at a The edge-on debris disk has an unusual shape: its luminous distance of 103 ± 4 pc (Van Leeuwen 2007) and belongs to the intensity has a very asymmetric profile. The longest peak, point- Lower Centaurus Crux (LCC) group, which is a subgroup of the ing westward, extends up to 550 au, while the east edge reaches Scorpius–Centaurus (Sco-Cen) OB association (de Zeeuw et al. 120 au only (see Figs. 1 and 3 of Kalas et al. 2015). Conversely, 1999). The LCC group has a mean age of 17 Myr, with an age- below 120 au, the disk is more luminous on its east side than on its west side. This reversed asymmetry might suggest the pres- 1 http://exoplanet.eu ence of a spiral density wave extending over the whole disk, and Article number, page 1 of 16 A&A proofs: manuscript no. Article Table 1: Key characteristics of the HD 106906 system. a very stable orbit. On the other hand, if it migrates too close to the binary, it undergoes a close encounter with the stars and can be ejected. System component Mass Projected separation a a The binary is surrounded by a chaotic zone where no sta- HD 106906 AB ∼ 1:34 and ∼ 1:37 M 0.36-0.58 au b b ble circumbinary orbit is possible. Dvorak (1986) uses a semi- HD 106906 b 11 ± 2 MJ 732 ± 30 au Disk 0:067 M c from 65 ± 3 to ∼ 550 au d e analytical approach to compute the upper critical orbit (lower Moon radius of the stable zone) and lower critical orbit (upper radius of the chaotic zone) around two stars of same masses for differ- ent binary eccentricities, and found that this gap size typically References: a – Lagrange et al. 2016b, b – Bailey et al. 2013, c – Chen et al. 2011, d – Kalas et al. 2015, e – Lagrange et al. 2016a. ranges between two and three times the semi-major axis of the binary orbit. Numerical results for this mass ratio are missing, so that we computed the limits of the gap with our Swift_ HJS viewed edge-on from the Earth. Finally, a large cavity splits the package and compared them to the semi-analytical approach in disk into two debris belts. Chen et al. (2014) modeled the stars’ Fig 1. In each simulation, the evolutions of 10,000 test parti- excess emission and suggested 13.7 and 46 au for the radii of cles have been studied during 105 orbital periods of the binary. the belts. The latter likely corresponds to the one imaged by La- The particles have been randomly chosen with semi-major axes grange et al. (2016a) and Kalas et al. (2015) at ∼ 50 au. between 1.5 and four times the binary semi-major axis aB, ec- Despite the richness of the observations, the geometry and centricities between 0 and 0.1, and inclinations with respect to kinematics of the whole system are strongly underconstrained. the binary orbital plan between 0 and 3 ◦. The time step has been If the actual planet-binary distance is less than 1,000 au, then chosen to be 1/20 of the binary orbital period. the orbit inclination with respect to the plane of the disk must be significant (20 degrees). However, a coplanar configuration can- not be excluded, but the separation should then be around 3,000 au. In any case, the large separation between the planet and the central binary, as well as the possible misalignment between the planet orbit and the debris disk, challenges classical mechanisms of planet formation. According to current theories, planet formation takes place in the primordial gaseous disk. However, as we mentioned above, forming a giant planet via core accretion or gravitational instabil- ity at 700 au or more from any central star appears very unlikely, first due to the lack of circumstellar gas at that distance, and sec- ond because the corresponding formation timescale would ex- ceed the lifetime of the gaseous disk. The disk asymmetries (in particular the suspected spiral structure) indicate strong ongoing dynamical interaction with the dust. This may suggest that the planet did not form where it resides today, but may have formed Fig. 1: Chaotic zone (in dark gray) as a function of the binary inside and be scattered afterwards. The recently discovered bi- eccentricity, for binary components of same masses. The lighter nary nature of HD 106906AB is indeed a source of potentially part designates a critical zone, where some test particles can sur- strong dynamical perturbations that could trigger planet scatter- vive. The red lines represent the lower and upper critical orbit ing. parabolic fits found by Dvorak (1986) in its study of circumbi- The purpose of this paper is to investigate both analytically nary planet stability.