arxiv: v3 [astro-ph.ep] 3 Oct 2013
|
|
- Belinda Dawson
- 5 years ago
- Views:
Transcription
1 Analysis of the motion of an extra-solar planet in a binary system arxiv:.3843v3 [astro-ph.ep] 3 Oct 03 Eva Plávalová Astronomical Institute Slovak Academy of Science Bratislava Slovak Republic plavala@slovanet.sk and Nina A. Solovaya Sternberg State Astronomical Institute Lomonosov Moscow State University Moscow Russia solov@sai.msu.ru ABSTRACT More than 0% of extra-solar planets EPs) orbit in a binary or multiple stellar system. We investigated the motion of planets revolving in binary systems in the frame of the particular case of the three body problem. We carried out an analysis of the motion an EP revolving in a binary system by following conditions; a) a planet in a binary system revolves around one of the components parent star) b) the distance between the star s components is greater than between the parent star and the orbiting planet ratio of the semi-major axes is a small parameter) c) the mass of the planet is smaller than the mass of the stars but is not negligible. The Hamiltonian of the system without short periodic terms was used. We expanded the Hamiltonian in terms of Legendre polynomial and truncated after the second order term depending on only one angular variable. In this case the solution of this system was obtained and the qualitative analysis of motion was produced. We have applied this theory to real EPs and compared to the numerical integration. Analyses of the possible regions of motion are presented. It is shown that the case of the stable and unstable motion of the EPs are possible. We applied our calculations to two binary systems hosting an EP and calculated the possible values for their unknown orbital elements. Subject headings: celestial mechanics planets and satellites: dynamical evolution and stability planets and satellites: individual 6 Cyb HD9994). Introduction More than half of all main sequence stars reside in binary or multiple systems Duquennoy & Mayor 99). At least 0% of the discovered extra-solar planets EPs) have been observed to be orbiting binary or multiple stellar systems Schneider et al. 0). Orbits in multiple stellar systems can have a miscellaneous in many cases unbelievable and improbable architecture. We targeted the binary stellar systems which are hosting EPs. For these systems Dvorak 986) defined three different possibilities for stable planetary orbit: S-orbit S-type) the Satellite-type orbit where an EP orbits one of the stars; P-orbit P-type) the Planet-type orbit where an EP surrounds both stars; L-orbit L-type) the Libration type orbit where an EP is librating around one of the triangular Lagrangian points. We have focused on an S-type orbit in this paper. The EP s motion in such system was considered in the frame of the general three body problem.. Setting up the problem The motion of an EP is investigated in the frame of the general three-body problem where the planet in the binary system revolves around one of the components a parent star) and the mass of the planet is much smaller than the mass of the stars but is not negligi-
2 ble. The distance between the binary components is greater than between the parent star and the orbiting planet the ratio the semi-major axes is a small parameter). The motion is considered in the Jacobian coordinate system and the invariable plane is taken as the reference plane. We used the Delaunay canonical elements L i G i H i l i g i and h i i= for the planet s orbit i= for the distant star s orbit). They can be expressed through the Keplerian elements as L i =β i ai G i = L i e i H i= G i cos I i l i = M i g i =ω i h i =Ω i. Where β = k m 0 m m0 + m = kµ ) β = k m 0+ m ) m m0 + m + m = kµ. ) In the previous expression the notation has the usual meaning: m 0 m the masses of the stars m the mass of the planet k the Gaussian constant a i the semi-major axis e i the eccentricity M i the mean anomaly I i ω i Ω i are the angular variables to the observation s plane andg i the argument of the pericenter in the invariable plane this plane is perpendicular to the angular momentum of the system). The eccentricities of the star s and planet s orbits can have any value from 0<e i <. In the general case the motion is defined by the masses of components and by the six initial values of the Keplerian elements of the planet and the distant star. The solution of a task using the Hamiltonian without short-periodic terms was obtained in the hyper elliptic integrals by the Hamilton-Jacobi method Orlov & Solovaya 988). The short-periodic terms are small and have no significant influence on the dynamic evolution of the system the values of which are less than ±0 3 which is less than the precision capabilities of the observations. The secular and essential long periodic terms are included into the Hamiltonian. Hamiltonian expanded in the terms of the Legendre polynomials and truncated after the second order terms carries the following form: F= γ + γ L L 6 γ 3 L 4 L 3 G3 [ 3q ) 5 3η ) 5 q ) η ) cosg ) ] 3) where the coefficientsγ γ andγ 3 depend on mass as follows: γ = β4 γ = β4 γ 3 = k β 6 µ µ 4) µ µ β 4 the parameter η= e 5) and the cosine of the angle between the plane of the planet s orbit and the plane of the distant star s orbit q= c G G G G. 6) Where c is the constant of the angular momentum of the system and c= C x + C y + C z 7) where C x = C y = C z = ) β i a i e i sin Ii sinω i i= ) β i a i e i sin Ii cosω i i= ) β i a i e i cos Ii. 8) i= The Hamilton s equation truncated after the second order terms depends on the one angular variableg only. It is necessary to integrate two equations from the canonical system of the differential equations. dg dt = F g dg dt Then the other equations will be as follows: dl i = 0 dt dg = 0 dt dc dt = 0 = F G. 9) dl i = F dt L i dg dt = F G dh = F dt c 0) where h=h and h = h+80. If system 9) is integrated the angular variables l i g and h will be defined from 0) in quadratures. For the solution of
3 system 9) we used the Hamilton-Jacobi method. According to this method we must find the complete integral W= W l l g g h A A A 3 A 4 A 5 ) ) of the equation in the partial derivations W +Φ g W W W W W ) = 0 ) t l l g g h whereφ= F and A A A 3 A 4 A 5 are the arbitrary constants. The general solution of the system is presented as follows: L i = W l i B i = W A i G i = W g i c= W h. 3) So as the Hamiltonian does not depend on l i g and h the integral ) is as follows: W=ε t t 0 )+A l +A l +A 4 g +A 5 h+w g ) 4) whereεis a new constant and W g ) is a new function which we must define. When we substituted the function W g ) defined by 4) to ) we found the equation to satisfy the function W : ε+φ g A A W g ) A 4 A 5 ) = 0 5) where W g ) is the derivation. Using equation 3) we rewrote previous equation to ε= γ A + γ A 6 γ 3 A 3 A 3 6) A3 4 where the constant A 3 has the following form: A 3 = 3 A 5 A 4 ) W 5 3 W 4 A 4 W A 5 A 5 A 4 ) W 4 A 4 W W cos g A ). 7) The equation 7) is an ordinary differential equation of the first order. If the solution is found then the general solution of the canonical system is: A 4 L = A L = A B = ε A t t 0 )+l + W A B = ε A t t 0 )+l + W A G = W g B 3 = ε A 3 t t 0 )+ W A 3 G = A 4 c=a 5 B 4 = ε A 4 t t 0 )+g + W A 4 B 5 = ε A 5 t t 0 )+h+ W A 5. 8) In the third line of the system 8) one can see the dependence between ξ and time t where ξ= W A = G. 9) A After the differentiation and the algebraic operations we obtained the following equation connecting ξ and t: G ξ ξ dξ = B 3 + A 6 γm e ) 3 n t t 0 ). 0) Where m = n /n and n n the mean motions of the planet and the distant starγ=m /m 0 + m + m ) while A and B 3 are the constants of integration. In this approximation we obtained an exact solution. We used this solution as an intermediary orbit for the planet s motion in which the second order perturbations are included. The expressions for a and a do not contain secular terms and hence are restricted in their time evolution. The change of the eccentricity of a planet is described by the expression e = ξ. In the case when the eccentricity of a planet s orbit can change to almost as much as the planet s pericenter is near to or beyond the Roche limit and large perturbations and tidal forces lead to the destruction of the planet. Such a system will be dynamically unstable. On the left side of the equation 0) under the integral in the denominator is the square root from the polynomial of the fifth order. This polynomial can be presented as the product of the two polynomials the second and the third orders = f ξ) f 3 ξ) where ξ= e and e is the eccentricity of the planet s orbit. The determination of the regions of the motion are possible when the roots of the equations f ξ)=0 f 3 ξ)=0 are found and the signs of the functions defined. 3
4 3. Equations of the second and third order We investigated the roots of the equations f ξ)=0 and f 3 ξ)=0with following forms: ) + f ξ)=ξ c + 3G ξ+ c G) A 3) G 3.. The equation of the second order We rewrote equation of the second order ) in the following form: ) f ξ)=ξ c + 3G ξ+ [ ) + + c + 3G + ] 3 G h 6) and f 3 ξ)=ξ 3 c + G + 5 ) ξ + 4 [ 5 ) + c + G + c G ) ] 6 G 0+ A 3) ξ 5 4 c G ). ) where h is difference between values A 3 and A 3crit Then h=a 3 A 3crit. 7) [ G h=3 η ) 8G 4G qη η + +0G ) ] q sin g. 8) Where and We note that c= c G > L ) a e G = G L = β β A 3 = 6η 0 q 0 6 η 0) ) a 3) [ 5 q 0) sin g 0 ]. 4) η 0 = e 0 5) where e 0 is the initial value of the eccentricity of the planet q 0 is the initial value of the cosine of the mutual inclination between the orbits of the planet and the distant star andg 0 is the initial value of the argument of the pericenter of the planet s orbit in the invariable plane. So the orbit is assumed to be elliptic 0<e < then 0 < ξ <. Value ofξ = corresponds to a circular motion. Such a case was not considered in this paper. For the determination of the boundaries of the regions of the possible motion it is necessary to find the roots of the equations ) and ) and to define the signs of the function in the interval between the roots. We named the meaning of the constant A 3 as A 3crit for which the square equation has the root equal to one. We found A 3crit when f ξ)=0 andξ=. We note the roots of the equation 6) asǫ andǫ andǫ <ǫ. The coefficient by ξ is always less than zero. This equation has no negative roots and the free term is also greater then 0. The roots are ) ǫ = c + 3G ) ± c + 3G + 3 G h 9) andǫ >> always. Derivation dǫ )/dh) is always positive. Soǫ is a growing function of h. When h=0 thenǫ = we conclude for h<0 thenǫ < and for h>0 thenǫ >. For the definition of the branch of the parabola we found the extreme of the function in the point which is validξ=c + 3G. The value f ξ)<0 and the branches of this parabola are directed up. 3.. Equation of the third order We rewrote equation ) using h in the form: ) 5+8c + 4G ξ + f 3 ξ)=ξ 3 [ c + G ] 6 G h ξ 5 4 c G ) c G ). 30) 4
5 Then multiplied this equation by 4 and identified in its left part the termξ f ξ). It is possible to rewrite equation 30) to take the following form: [ ) ξ f 3 ξ)=5 ξ) c+ G ][ξ ) ] c G ξ f ξ). 3) For the determination of the qualitative characteristics of the motion it is necessary to define how the roots of the equation of the third order f 3 ξ)=0 are located in relation to the points of the axes: ξ=0 ξ=ǫ ξ= ξ=ǫ. 3). If h<0 then we have: a)ξ=0 f 3 0)<0. b)ξ=ǫ f 3 ǫ )>0. c)ξ= f 3 )>0. d)ξ=ǫ f 3 ǫ )>0. From Sturm theorem e.g. Dörrie 965) it follows that between 0 andǫ lay at least one root of the third order equation.. If h>0 then we have: a)ξ=0 f 3 0) 0. b)ξ= f 3 )<0. c)ξ=ǫ f 3 ǫ )<0. d)ξ=ǫ f 3 ǫ )>0. The behaviours of the roots of the second and third order equations depend on the sign of h and valid:. for h<0 is valid 0 ǫ 3 <ǫ < <ǫ 4 <ǫ 5 <ǫ.. for h>0 is valid 0 ǫ 3 <ǫ 4 < <ǫ <ǫ 5 <ǫ. We have two roots which are less than therefore we identified the roots in ascending order as ξ ξ ξ 3 ξ 4 ξ 5. This means that e can change in limits from e min = ξ to e max = ξ. We showed the behaviour of these functions for the hypothetical value c G and h on the Figure for h<0 and on Figure for h>0. 4. The investigation of the value h Consider the value h which has the following form: h= 3 η ) G [ 8G 4G qη η + +0G ) ] q sin g. 33) When the valueg = 0 then value h reaches minimum h min = 3 η ) [ 8G ] 4G qη η. 34) G The value of the h min is always negative. In this case we have the two roots valued less than the one root of the equation of the third order and one root of the square equation. Ifg =π/ then h has a maximum value h max = 3 η )[ η 4G qη+ G +G 0G q ]. 35) The value of h max can either be negative or positive. For the caseg =π/ we rewrote equations ) and f Ξ for h Ξ f Ξ for h Ξ Fig.. Behaviour of the functions f ξ) dashed line) and f 3 ξ) solid line) for h>0. Fig.. Behaviour of the functions f ξ) dashed line) and f 3 ξ) solid line) for h<0. 5
6 ) in the following form: f ξ)= ) ξ η 0 ξ η 0 4G η 0 q 0 + 8G + 4G ) + 5q 0 ) η 0 ) + 6G 36) and f 3 ξ)= ) [ ξ η 0 ξ 5 ) ξ η 4 0) ] 3G ξ+ G η 0 q 0 5 4ξ)+5G q 0. 37) When one of the roots of the equation 37) carries the value ofξ = η 0 the other roots have to satisfy the equation ξ 5 ) ξ η 4 0) 3 G ξ+ G η 0 q 0 5 4ξ)+5 G q 0 = 0. 38) For h max the initial valueη =η 0 is one of the boundary limits for change in the value ofξ. We should establish thatξ=η 0 is the least root and the second root has a value less than. We therefore substitute the valueξ= η 0 in 38). If the obtained expression is negative then η 0 should become the second root in equation 37) if positive then it should become the least root. The left part of the equation 38) forξ=η 0 is [ 5q 0 3η 0 + ) ] η 0 q 0 5 4η 0. 39) G Equating this expression to zero and solving the obtained equation concerning q 0 we find the value of q 0 of whichη 0 is the root of expression39). So [ η 0 4η 0 60G 5 + ) ] 5 4η 0 q=. 40) 0G We denote the meaning of roots 40) as q 0 for the minus sign before the root term and the q 0 for plus sign. If the value of q 0 lays inside the region of q 0 < q 0 < q 0 expression 39) is negative and the value of the root lays between the other two roots. If either of the conditions q 0 < q 0 or q 0 > q 0 is valid thenξ=η 0 is the least root of equation 38). In this case e max = η 0 and the maximum value of the eccentricity of the planet s orbit can not exceed the initial value of eccentricity. The orbit of the EP may be dynamical stable. When the starting value of the cosine of the mutual inclination is q 0 = η 0 /G then from Eg.6) c G 0 and e. 5. Application of the theory on real extra-solar planets The knowledge of the six pairs of Keplerian elements of the orbits of the system allow us to investigate the character of the evolution of the planet s orbit and the possible conditions of stability. They may be presented by the orbital parameters which we obtain from the analytical theory. They are the angle of the mutual inclination between orbits of the planet and star the angular moment of the star and the maximum value of the eccentricity of the planet s orbit. The growth of the eccentricity of the orbit could lead to the destruction of the orbit in pericenter from tidal forces. In the catalogue of EPs the data for the longitude of the ascending node and the value of inclination is generally absent. The existing observational techniques do not allow estimation of these two elements unambiguously. Our theory allows us to define a range of possible values for these unknown elements by which the planet s eccentricity does not increase. In our application of the theory we have selected two binary stellar systems with a planet revolving around one of the components. The first is system 6 Cyg with an interesting planet 6 Cyg Bb which has the argument of pericenter close to 90. The second system is HD Cyg 6 Cyg A HD 86408) and 6 Cyg B HD 8647) are both members of a well-known wide binary system with stars spectral types G.5V and G3V. For the calculations we used orbital elements used by Hauser & Marcy 999) which were published in the Sixth Catalogue of Orbits of Visual Binary Stars Mason & Hartkopf 00). The semi-major axis a was calculated to be AU using the values.4 pc Fuhrmann et al. 998) for the distance and 35.4 for the angular separation. We made a revision calculation for this semi-major axis using the third Keplerian law and the value of 35.7 year for a period P Mason & Hartkopf 00). We reached a value of a = AU. There is an obvious difficulty in determining orbital parameters for periods greater than 0000 years. Only a small fraction of the orbit has transpired since astrometric measurements were first 6
7 Table : Initial orbital elements of the system 6 Cyg. star 6 Cyg A planet 6 Cyg Bb Mass [M S un ].53 Mass [M Jup sin I ].68±0.07 Semi-major axis [AU] ±0.07 Eccentricity ±0.0 Inclination Ascending node Argument of pericenter ±. Period 35.7 years days Note. For the calculations we copied the orbital elements for 6 Cyg A used by Hauser & Marcy 999) and for the planet we used elements from Schneider et al. 0). We used the planet s semi-major axis value.693 AU for the calculations. The mass of the parent star 6 Cyg B is.0±0.04m S un. Table : Our proposal for possible orbital elements for planet 6 Cyg Bb. Prograde orbit Retrograde orbit Mass [M Jup ].38± ±0.04 Semi-major axis [AU] Eccentricity 0.689± ±0.0 Inclination 45 ± 35 ± Argument of pericenter 83.4±. 83.4±. Period [day] obtained in the early 800s. So we used the value AU derived from the distance for the following calculations because we considered the determination period less precise than the determination of the values for angular separation and distance. The secondary star in the 6 Cyg binary system 6 Cyg B parent star) with a mass of.0±0.04 M S un Fuhrmann et al. 998) is known to host a giant gas planet 6 Cyg Bb with taxonomy class J0.W7 Plávalová 0). The EP has a minimum mass M Jup sin I =.68± 0.07M Jup in an orbit with high eccentricity e = Wittenmyer et al. 007). The initial conditions for investigating this planet we took from The Extrasolar Planets Encyclopaedia Schneider et al. 0). The values of the inclination and ascending note for the planet are unknown. In Table are shown the orbital elements for the planet and distant star. We made a revision calculation for the semi-major axis of the planet using the third Keplerian law and values for the orbital period of the planet which is days. It is known the minimum mass of the planet m =.68 M Jup. We varied the elements I from 0 to 80 at a time andω from 0 to 360 also for each degree. We recognized the value of the planet s inclination I for which the maximum value of the planet s eccentricity e is close to the initial value According to our calculation we received the value.693 AU for the semi-major axis.38m Jup for the mass and 45 or 35 for the inclination of the planet. These values are presented in Table and were used for our next calculations. The value of the ascending node of the planet s orbitω is unknown. But the orientation of the orbit is defined by the node therefore the different ways of the dynamical evolution are possible. We varied the ascending node of the planet from 0 to 360 by step. We recognized the region of the values ofω for which the planet s orbit would be stable or unstable. For the stability criteria we used the Roche limit. This is the minimum distance which a planet can approach its parent star without being torn apart by tidal forces. For the calculation of the Roche limit d R in our paper we used the equation published by Eggleton 983). 7
8 e max 0.9 Roche limit e max 0.9 Roche limit 0.7 I =45 I = q q I =5 I = q 0 I = q I = I =5 0.8 q 0 I = Ω Fig. 3. The planet 6 Cyg Bb. The evolution of the maximum value of the planet s eccentricity e and the cosine of the mutual inclination between the planet s orbit and the distant star s orbit q versus the ascending node of the planetω for the prograde I = 45 and retrograde I = 35 planet orbits. The Roche limit is plotted by a dashed line. 0.8 q Ω Fig. 4. The planet HD9994 b. The evolution of the maximum value of the planet s eccentricity e and the cosine of the mutual inclination between the planet s orbit and the distant star s orbit q versus the ascending node of the planetω for the prograde I = 65 and retrograde I = 5 planet orbits. The Roche limit is plotted by a dashed line. 0.49µ 3 d R = 0.6µ 3+ ln ) 4) +µ 3 whereµ=m /m 0. The results of our calculations are presented in Figure 3. and 5. If the ascending node isω ) for prograde motion I = 45 andω ) for retrograde motion I = 35 the planet reaches the Roche limit in its pericentre d R = AU) with an eccentricity of e In such a pericenter large perturbations and tidal forces drastically effect the planet and lead to its destruction. For the values I = 45 andω = 35 the maximum and minimum values of eccentricity e are equal to the initial value e = In this case q 0 = < q 0 = With these elements the planet s orbit would be stable. When the eccentricity of the planetary orbit grows close to the Roche limit large perturbations in the planet s pericenter are affected. In such cases measurements based on point mass bodies can be inaccurate rather a dynamic theory for real dimensional bodies must be used. 5.. HD 9994 The binary system HD Ceti ADS 406 AB) contains an A component: a yellow-white dwarf with a mass of.34 M S un and a B component: a red dwarf with a mass of 0.37 M S un. For the distant star we used Keplerian elements published by Hale 994). To define its semi-major axis a we used two methods. Firstly we derived this value using stellar parallax where the distance was.38 pc Schneider et al. 0) and 6.77 for the angular separation Hale 994) and resulted in a semi-major axis of 5.5 AU. Secondly with the application of Keplerian law using the values from Table 3 we calculated the value to be 5.37 AU. The difference between these two values is negligible we applied the value a = 5.5 AU in 8
9 Table 3: Initial orbital elements of the system HD9994. star HD9994 planet HD9994 b Mass [M S un ] 0.37 Mass [M Jup sin I ].68 Semi-major axis [AU] Eccentricity ±0.04 Inclination 4. Ascending node 84.3 Argument of pericenter ±8 Period 40 years 535.7±3. days Note. For the distant star we used Keplerian elements published by Hale 994). For our calculations we used the orbital elements published by Mayor et al. 004). We used the planet s semi-major axis value.473 AU for the calculations. The mass of the parent star HD9994 is.34m S un. Table 4: Our proposal for possible orbital elements for planet HD9994 b. Prograde orbit Retrograde orbit Mass [M Jup ].86± ±0.045 Semi-major axis [AU] Eccentricity 0.300± ±0.04 Inclination 65 ±3 5 ±3 Argument of pericenter 4 ±8 4 ±8 Period [day] 535.7± ±3. the following calculations. The planet HD 9994 b was discovered in 000 Queloz et al. 00) and is orbiting the A component. The taxonomy class is J0.G3 Plávalová 0) and the minimum mass is.68 M Jup. This planet orbiting with a semi-major axis of.4 AU with quite a high eccentricity 0.3 ± For our calculations we used the orbital elements published by Mayor et al. 004). In Table 3 the initial Kepler orbital elements for the planet and its distant star are shown. As with the first system we revised the calculation for the semi-major axis of the planet using Keplerian law. With the values listed in Table 3 we valued the planet s semi-major axis at a =.47 AU. If we varied the mass of the planet to 5 M Jup we would get a value of.48 AU. We decided to use the value.47 AU for the semi-major axis and.68 M Jup for the mass of the planet in our calculations. As in the case of 6 Cyg b we varied the elements I from 0 to 80 at a time andω from 0 to 360 also for each degree. We recognized the value of the planet s inclination I for which the maximum value of the planet s eccentricity e is close to the initial value According to our calculation we received the value a =.47 AU for the semi-major axis.86m Jup for the mass and 65 or 5 for the inclination of the planet. These values are presented in Table 4 and were used for our next calculations. The results of our calculations are presented in Figure 4. and 6. If the ascending node isω ) orω ) for prograde motion I = 65 andω 7 0 ) orω ) for retrograde motion I = 5 the planet reaches the Roche limit in its pericentre d R = 0.05 AU) with an eccentricity of e In such a pericenter large perturbations and tidal forces drastically effect the planet and lead to its destruction. For the values I = 65 andω = 63 the maximum and minimum values of the eccentricity e are equal to the initial value e = In this case q 0 = 0.999<q 0 = With these elements the planet s orbit would be stable. 9
10 e Ω =85 Ω =0 Roche limit Ω =34 e Ω =65 Ω =30 Roche limit Ω =63 r p [AU] Ω =0 Ω =34 0. Ω =85 Roche limit t[myr] Fig. 5. The planet 6 Cyg Bb. The evolution of the planet s eccentricity e and the pericenter distance r p within the interval of years. The curves are the result of the numerical integration. For the case if I = 45 andω = 0 or I = 45 andω = 34 the planet s orbit does not reach the Roche limit. For the case if I = 45 andω = 85 the planet stays within the Roche limit. The boundaries of the gray zones for all three cases were computed from the our theory. For the case ofω = 34 the gray zone is comparable with the line. The Roche limit is plotted by a dashed line. r p [AU] Ω =30 Ω =63 0. Ω =65 Roche limit t[myr] Fig. 6. The planet HD9994 b. The evolution of the planet s eccentricity e and the pericenter distance r p within the interval of years. The curves are the result of the numerical integration. For the case if I = 65 andω = 30 or I = 65 andω = 63 the planet s orbit does not reach the Roche limit. For the case I = 65 andω = 65 the planet stays within the Roche limit. The boundaries of the gray zones for all three cases were computed from the our theory. For the case ofω = 63 the gray zone is comparable with the line. The Roche limit is plotted by a dashed line. 6. Comparison of our theoretical results with numerical integration For confirmation of the obtained analytic results we compared them with the results of numerical integration. The equations of the motion of the systems were numerically integrated from the initial date 000 January forward within the interval of years using the Everhart s integrator Everhart 985). We have integrated the three variants ofω for the prograde and retrograde orbits for both systems. For 6 Cyg Bb we used the initial value of the planet s inclination i = 45 and the planet s ascending nodeω = 0 Ω = 34 andω = 85. The results obtained by the theory and by the numerical integration agree quite well. For an illustration the evolution of the eccentricity e and the pericenter distance r p = a e ) of the planet 6 Cyg Bb are presented in Figure 5. For HD 9994 b see Figure 6) we used as the initial values i = 65 andω = 30 Ω = 65 andω = 63. As in the case of 6 Cyg Bb the results of the numerical integration are the same as the results we obtained by the analytical theory. We obtained identical results from the numerical integration and analytical theory for the retrograde orbits of the both systems too. 7. Comparison with previous works In our investigation we have used the methods of the classical celestial mechanics developed by Hamilton. We come from the Principle of Determinacy PD)Lidov 00) that if the initial conditions of each object in a mechanical system are defined at any instant of time then their further behavior is defined expressly. Note that Hamilton s equations are first order in the time derivative which makes them more convenient for computation. The Hamiltonian 3) of the system 9) and 0) without short-periodic terms permits the solution in which the secular and the long-periodic 0
11 terms are taken into account in the intermediate orbit and allows the close approach of the EP to the star to be established. The Long-periodic stability of planets in a binary system was investigated by Holman & Wiegert 999). The planets are taken to be test particles moving in the field of an eccentric binary system. This study investigated the orbital stability numerically with the elliptic restricted three body problem. We made a comparison with their results concerning planet 6 Cyg Bb. We used the elements of this planet for the application of our theory. Holman & Wiegert 999) integrated differential equations of the motion usingω = 0 and obtained instability between 0 7 and 0 9 years. We varied this element from 0 till 360 and obtained instability aroundω = 0 too. The secular dynamics of massless particles orbiting a central star and perturbed by a secondary star component with high eccentricity have been investigated Giuppone et al. 0). They used the Lie series perturbation scheme restricted to second order in the small parameter. Their results have shown that the second-order secular dynamic reproduces the behaviour of a planet with good precision. We eliminated the short-periodic terms by von Zeipel s method in the general three-body problem. This method allows us to estimate values for the short periodic terms; these values are less than±0 3. The short-periodic terms are small and do not influence the evolution and stability. Our condition of stability is valid at any time interval. Our theory may also be used for EPs with a large mass. In the Sun asteroid Jupiter problem the Hamiltonian describing the motion of the massless asteroid in the heliocentric coordinates was used by Innanen et al. 997). They used the so-called Kozai mechanism which shows the behaviour of the eccentricity as a function of the initial inclination. In the restricted three-body problem with the averaging over the mean anomalies the perturbation function there are the integrals: e sin i=const. and ) e 5 sin i sin ω = const. 4) These were obtained at nearly the same time by Lidov 96) and Kozai 96). The integrals permit us to execute qualitative analysis of the family of the phase trajectories. In our case the connection between the eccentricity the inclination and the argument of pericenter is more complex depends on all Keplerian elements and defined by the equation 33). It allows us to arrive at the some conclusion similar to the results of Innanen et al. 997). The increase of the eccentricity occurs even if the third body is very distant and the perturbation is small. The time scale needed to observe this phenomenon is quite lengthy for example the sudden increase in the eccentricity of Neptune is observed at 00 Myr. Innanen et al. 997). But what is clear is that instability only occurs when the mutual inclination of a planet s orbit and the distant star s orbit is high. So we suppose for the stability of a planet s orbit the conditions are demanded: the initial value of the cosine of the mutual inclination q 0 < q 0 or q 0 > q 0 where q 0 and q 0 are roots of expression 40) and in this case c G > Conclusion We have shown that an EP revolving in a binary system around one of the components may move on a stable or unstable orbits. The conditions for a stable motion depend on the orbital parameters which can be calculated from the formulas of the previous sections. There are the angle of the mutual inclination between orbits the angular momentum of the distant star and the maximum value of the eccentricity of the planet s orbit. When the value of the planet s eccentricity is close to one in the pericenter the planet reaches the Roche limit. The tidal forces of the star destroy the planet. We have suggested the possible values for the unknown elements with which the orbit of the planet would remain stable. The results of our calculations are presented in Figures 6. The results of the numerical integration support the results obtained by the analytical theory. For 6 Cyg Bb there are three possible regions of stability. Firstly for the prograde orbit I = 45 and the ascending nodeω [4 49 ]. Secondly for the retrograde orbit I = 35 and the ascending node Ω [0 74 ] orω [ ]. We proposed the planet s mass to be.38±0.04 M Jup and the value of its semi-major axis equal to.693 AU for which the third Keplerian law is valid. For the second system HD9994 the values of
12 inclination and the ascending node of the planet for which the motion is stable are: for the prograde motion I = 65 and the ascending node Ω [3 46 ] orω [ ]; for the retrograde motion I = 5 and the ascending node Ω [0 35 ] or Ω [358 0 ] or Ω [0 7 ]. For the planet s mass we proposed the value of.86±0.045 M Jup and for the planet s semimajor axis a value of.47 AU. These values are in accordance with the third Keplerian law. This work was supported by the Slovak Research and Development Agency under the contract No. APVV REFERENCES Dörrie H Great Problems of Elementary Mathematics: Their History and Solution Eggleton P. P. 983 ApJ Mason B.D. & Hartkopf W.I. 00 The U.S. Naval Observatory Double Star CD 00.0 American Astronomical Society Meeting Abstracts 33 3 Mayor M. Udry S. Naef D. et al. 004 A&A Orlov A. A. & Solovaya N.A. 988 Astrophysics and Space Science Library Vol: 40 IAU Colloq. 96: The Few Body Problem ed. Valtonen M.J. 43 Plávalová E. 0 Astrobiology 36 Queloz D. Mayor M. Udry S. et al. 00 The Messenger 05 Schneider J. Dedieu C. Le Sidaner P. Savalle R. & Zolotukhin I. 0 A&A 53 A79 Wittenmyer R. A. Endl M. & Cochran W. D. 007 ApJ Everhart E. 985 Dynamics of Comets: Their Origin and Evolution Proceedings of IAU Colloq. 83 held in Rome Italy June Edited by Andrea Carusi and Giovanni B. Valsecchi. Dordrecht: Reidel Astrophysics and Space Science Library. Vol Duquennoy A. and Mayor M. 99 A&A Dvorak R. 986 A&A Fuhrmann K. Pfeiffer M. J. & Bernkopf J. 998 A&A Giuppone C.A. Leiva A. M. Correa-Otto J. & Beaugé C. 0 A&A 530 A03 Hale A. 994 AJ Hauser H.M. & Marcy G.W. 999 PASP 3 Holman M.J. & Wiegert P.A. 999 AJ 7 6 Innanen K.A. Zheng J. Q. Mikkola S. & Valtonen M. J. 997 AJ 3 95 Kozai Y. 96 AJ Lidov M.L. 96 Planet. Space Sci Lidov M.L. 00 Course of lectures on theoretical mechanics Kurs lektsiy po teoreticheskoy mekhanike) In Russian Fizmatlib Moscow 33 ISBN This -column preprint was prepared with the AAS L A TEX macros v5..
Kozai-Lidov oscillations
Kozai-Lidov oscillations Kozai (1962 - asteroids); Lidov (1962 - artificial satellites) arise most simply in restricted three-body problem (two massive bodies on a Kepler orbit + a test particle) e.g.,
More informationarxiv: v1 [astro-ph.ep] 3 Apr 2018
Astronomy& Astrophysics manuscript no. zanardi_ c ESO 28 September 2, 28 The role of the general relativity on icy body reservoirs under the effects of an inner eccentric Jupiter M. Zanardi, 2, G. C. de
More informationJoseph Castro Mentor: Nader Haghighipour
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
More informationDynamical properties of the Solar System. Second Kepler s Law. Dynamics of planetary orbits. ν: true anomaly
First Kepler s Law The secondary body moves in an elliptical orbit, with the primary body at the focus Valid for bound orbits with E < 0 The conservation of the total energy E yields a constant semi-major
More informationResearch Paper. Trojans in Habitable Zones ABSTRACT
ASTROBIOLOGY Volume 5, Number 5, 2005 Mary Ann Liebert, Inc. Research Paper Trojans in Habitable Zones RICHARD SCHWARZ, 1 ELKE PILAT-LOHINGER, 1 RUDOLF DVORAK, 1 BALINT ÉRDI, 2 and ZSOLT SÁNDOR 2 ABSTRACT
More informationInduced Eccentricities of Extra Solar Planets by Distant Stellar Companions
Induced Eccentricities of Extra Solar Planets by Distant Stellar Companions T. Mazeh Wise Observatory, Tel Aviv University, Tel Aviv, Israel Abstract. This paper explores the assumption that the high eccentricities
More informationUNIVERSITY of CALIFORNIA SANTA CRUZ
UNIVERSITY of CALIFORNIA SANTA CRUZ NUMERICAL AND ANALYTIC RESEARCH INTO THE FORMATION OF THE HD 80606B PLANETARY SYSTEM A thesis submitted in partial satisfaction of the requirements for the degree of
More informationISIMA lectures on celestial mechanics. 3
ISIMA lectures on celestial mechanics. 3 Scott Tremaine, Institute for Advanced Study July 2014 1. The stability of planetary systems To understand the formation and evolution of exoplanet systems, we
More informationAnalysis of Radial Velocity Measurements
Analysis of Radial Velocity Measurements Sistemas Planetários - 2nd project - 2011/2012 João Faria 1. Introduction The existence of solar systems other than our own has been speculated for a long time,
More informationPervasive Orbital Eccentricities Dictate the Habitability of Extrasolar Earths
Pervasive Orbital Eccentricities Dictate the Habitability of Extrasolar Earths Ryosuke Kita, 1 Frederic Rasio, 1,2 and Genya Takeda 1 arxiv:1003.0633v2 [astro-ph.ep] 2 Mar 2010 1 Department of Physics
More informationarxiv:astro-ph/ v3 2 Mar 2006
Mon. Not. R. Astron. Soc. 000, 1?? (2006) Printed 5 February 2008 (MN LATEX style file v2.2) Pluto s moon system: survey of the phase space I. I. Nagy, Á. Süli and B. Érdi Department of Astronomy, Eötvös
More informationOrbital and Celestial Mechanics
Orbital and Celestial Mechanics John P. Vinti Edited by Gim J. Der TRW Los Angeles, California Nino L. Bonavito NASA Goddard Space Flight Center Greenbelt, Maryland Volume 177 PROGRESS IN ASTRONAUTICS
More informationImportance of the study of extrasolar planets. Exoplanets Introduction. Importance of the study of extrasolar planets
Importance of the study of extrasolar planets Exoplanets Introduction Planets and Astrobiology (2017-2018) G. Vladilo Technological and scientific spin-offs Exoplanet observations are driving huge technological
More informationSecular dynamics of extrasolar-systems
Secular dynamics of extrasolar-systems Marco Sansottera [1] [1] Università degli Studi di Milano [2] Namur Center for Complex Systems (naxys) Based on a research project in collaboration with A.-S. Libert
More informationarxiv: v1 [astro-ph] 24 Dec 2008
Mon. Not. R. Astron. Soc. 000, 1 7 (0000) Printed 6 June 2018 (MN LATEX style file v2.2) High Inclination Planets and Asteroids in Multistellar Systems P. E. Verrier 1 and N. W. Evans 1 1 Institute of
More informationEVOLUTIONS OF SMALL BODIES IN OUR SOLAR SYSTEM
EVOLUTIONS OF SMALL BODIES IN OUR SOLAR SYSTEM Dynamics and collisional processes Dr. Patrick MICHEL Côte d'azur Observatory Cassiopée Laboratory, CNRS Nice, France 1 Plan Chapter I: A few concepts on
More informationEccentricity pumping of a planet on an inclined orbit by a disc
Mon. Not. R. Astron. Soc. 44, 49 414 21) doi:1.1111/j.1365-2966.21.16295.x Eccentricity pumping of a planet on an inclined orbit by a disc Caroline Terquem 1,2 and Aikel Ajmia 1 1 Institut d Astrophysique
More informationOrbital Evolution in Extra-solar systems
Orbital Evolution in Extra-solar systems George Voyatzis Section of Astronomy, Astrophysics and Mechanics, Department of Physics, Aristotle University of Thessaloniki, Greece. Abstract Nowadays, extra-solar
More informationDynamic Exoplanets. Alexander James Mustill
Dynamic Exoplanets Alexander James Mustill Exoplanets: not (all) like the Solar System Exoplanets: not (all) like the Solar System Solar System Lissauer et al 14 Key questions to bear in mind What is role
More informationTERRESTRIAL TROJAN PLANETS IN EXTRASOLAR SYSTEMS
TERRESTRIAL TROJAN PLANETS IN EXTRASOLAR SYSTEMS Richard Schwarz Institute of Astronomy University of Vienna Türkenschanzstrasse 17 A-1180 Vienna, Austria schwarz@astro.univie.ac.at Abstract Keywords:
More informationSTABILITY OF HYPOTHETICAL TROJAN PLANETS IN EXOPLANETARY SYSTEMS
STABILITY OF HYPOTHETICAL TROJAN PLANETS IN EXOPLANETARY SYSTEMS Bálint Érdi, Georgina Fröhlich, Imre Nagy and Zsolt Sándor Department of Astronomy Loránd Eötvös University Pázmány Péter sétány 1/A H-1117
More informationTHE exploration of planetary satellites is currently an active area
JOURNAL OF GUIDANCE, CONTROL, AND DYNAMICS Vol. 9, No. 5, September October 6 Design of Science Orbits About Planetary Satellites: Application to Europa Marci E. Paskowitz and Daniel J. Scheeres University
More informationDynamics of possible Trojan planets in binary systems
Mon. Not. R. Astron. Soc. 398, 20852090 (2009) doi:1111/j.1365-2966.2009.15248.x Dynamics of possible Trojan planets in binary systems R. Schwarz, 1 Á. Süli 1 and R. Dvorak 2 1 Department of Astronomy,
More informationarxiv:astro-ph/ v1 16 Sep 1996
The stability of planets in the Alpha Centauri system Paul A. Wiegert 1, Matt Holman 2 1 Department of Astronomy, University of Toronto, Toronto, Canada 2 Canadian Institute for Theoretical Astrophysics,
More informationPlanetary system dynamics Part III Mathematics / Part III Astrophysics
Planetary system dynamics Part III Mathematics / Part III Astrophysics Lecturer: Prof. Mark Wyatt (Dr. Amy Bonsor on 9,11 Oct) Schedule: Michaelmas 2017 Mon, Wed, Fri at 10am MR11, 24 lectures, start Fri
More informationThe stability of planets in the Alpha Centauri system
The stability of planets in the Alpha Centauri system Paul A. Wiegert 1, Matt Holman 2 1 Department of Astronomy, University of Toronto, Toronto, Canada 2 Canadian Institute for Theoretical Astrophysics,
More informationarxiv: v1 [astro-ph.ep] 20 Nov 2018
Astronomy & Astrophysics manuscript no. A_A_33097_Busetti c ESO 2018 November 21, 2018 Stability of planets in triple star systems F. Busetti 1, H. Beust 2, and C. Harley 1 1 Univ. of the Witwatersrand,
More informationA study upon Eris. I. Describing and characterizing the orbit of Eris around the Sun. I. Breda 1
Astronomy & Astrophysics manuscript no. Eris c ESO 2013 March 27, 2013 A study upon Eris I. Describing and characterizing the orbit of Eris around the Sun I. Breda 1 Faculty of Sciences (FCUP), University
More informationCHAPTER 8 PLANETARY MOTIONS
1 CHAPTER 8 PLANETARY MOTIONS 8.1 Introduction The word planet means wanderer (πλάνητες αστέρες wandering stars); in contrast to the fixed stars, the planets wander around on the celestial sphere, sometimes
More informationSatellite Communications
Satellite Communications Lecture (3) Chapter 2.1 1 Gravitational Force Newton s 2nd Law: r r F = m a Newton s Law Of Universal Gravitation (assuming point masses or spheres): Putting these together: r
More informationSupplementary Materials for
www.sciencemag.org/content/345/622/1317/suppl/dc1 Supplementary Materials for Chaotic dynamics of stellar spin in binaries and the production of misaligned hot Jupiters Natalia I. Storch, Kassandra R.
More informationA SYMPLECTIC MAPPING MODEL FOR THE STUDY OF 2:3 RESONANT TRANS-NEPTUNIAN MOTION
3 A SYMPLECTIC MAPPING MODEL FOR THE STUDY OF 2:3 RESONANT TRANS-NEPTUNIAN MOTION K.G. HADJIFOTINOU and JOHN D. HADJIDEMETRIOU Department of Physics, University of Thessaloniki, 540 06 Thessaloniki, Greece
More informationWhite Paper. Terrestrial and Habitable Planet Formation in Binary and Multi-star Systems
Authors White Paper (Submitted to ExoPlanet Task Force) Terrestrial and Habitable Planet Formation in Binary and Multi-star Systems Nader Haghighipour Institute for Astronomy and NASA Astrobiology Institute,
More informationANNEX 1. DEFINITION OF ORBITAL PARAMETERS AND IMPORTANT CONCEPTS OF CELESTIAL MECHANICS
ANNEX 1. DEFINITION OF ORBITAL PARAMETERS AND IMPORTANT CONCEPTS OF CELESTIAL MECHANICS A1.1. Kepler s laws Johannes Kepler (1571-1630) discovered the laws of orbital motion, now called Kepler's laws.
More informationarxiv: v1 [astro-ph] 4 Jul 2007
Astronomy & Astrophysics manuscript no. pap5 c ESO 28 February, 28 Dispersal of planets hosted in binaries, transitional members of multiple star systems F. Marzari and M. Barbieri 2 arxiv:77.59v [astro-ph]
More informationASTB01 Exoplanets Lab
ASTB01 Exoplanets Lab Author: Anders Johansen Revision date: $Date: 2015/08/28 14:55:59 $ Planets orbiting stars other than the Sun are called exoplanets. Stellar light reflected off
More informationarxiv: v1 [astro-ph.ep] 3 Jul 2018
Astronomy & Astrophysics manuscript no. _main c ESO 2018 July 5, 2018 Centaurs potentially in retrograde co-orbit resonance with Saturn Miao Li 1, Yukun Huang 1, and Shengping Gong 1, School of Aerospace
More informationAST111, Lecture 1b. Measurements of bodies in the solar system (overview continued) Orbital elements
AST111, Lecture 1b Measurements of bodies in the solar system (overview continued) Orbital elements Planetary properties (continued): Measuring Mass The orbital period of a moon about a planet depends
More informationPlanetenbewegung in Sternsystemen. The Effect of Resonances. Part 2
Planetenbewegung in Sternsystemen The Effect of Resonances Part 2 Topics overview 1. Definition and examples of resonances 2. Disturbing function 3. Mean-motion resonance (MMR) 4. Secular resonance (SR)
More informationDynamical Stability of Terrestrial and Giant Planets in the HD Planetary System
Dynamical Stability of Terrestrial and Giant Planets in the HD 155358 Planetary System James Haynes Advisor: Nader Haghighipour ABSTRACT The results of a study of the dynamical evolution and the habitability
More informationTowards stability results for planetary problems with more than three bodies
Towards stability results for planetary problems with more than three bodies Ugo Locatelli [a] and Marco Sansottera [b] [a] Math. Dep. of Università degli Studi di Roma Tor Vergata [b] Math. Dep. of Università
More informationLecture 13. Gravity in the Solar System
Lecture 13 Gravity in the Solar System Guiding Questions 1. How was the heliocentric model established? What are monumental steps in the history of the heliocentric model? 2. How do Kepler s three laws
More information2 Ford, Rasio, & Yu. 2. Two Planets, Unequal Masses
2 Ford, Rasio, & Yu unlikely to have developed such a large eccentricity, since dissipation in the disk tends to circularize orbits. Dynamical instabilities leading to the ejection of one planet while
More informationAdvisor : Prof. Hagai Perets. Israel Institute of Technology, Haifa
Secular Dynamics in Hierarchical Three-Body Systems with Mass Loss and MassErez Michaely Advisor : Prof. Hagai Perets Israel Institute of Technology, Haifa Context Systematically exploring triple systems
More informationOn the Secular Evolution of Extrasolar Planetary Systems
On the Secular Evolution of Extrasolar Planetary Systems Marco Sansottera [a] [a] Namur Center for Complex Systems (naxys) Based on a research work in collaboration with Anne-Sophie Libert [a] JISD 2012
More informationCALCULATION OF POSITION AND VELOCITY OF GLONASS SATELLITE BASED ON ANALYTICAL THEORY OF MOTION
ARTIFICIAL SATELLITES, Vol. 50, No. 3 2015 DOI: 10.1515/arsa-2015-0008 CALCULATION OF POSITION AND VELOCITY OF GLONASS SATELLITE BASED ON ANALYTICAL THEORY OF MOTION W. Góral, B. Skorupa AGH University
More informationChaos in some young asteroid families
Chaos in some young asteroid families Rosaev A., Plavalova E Research and Educational Center of Nonlinear Dynamics, Yaroslavl State University, Yaroslavl, Russia Astronomical Institute Slovak Academy of
More informationObservations of Extrasolar Planets
Observations of Extrasolar Planets Hamilton 2005 Shay Zucker Observations of Extrasolar Planets Spectroscopic detection of exoplanets Emerging properties of the sample Transiting planets Future prospects
More informationSECULAR EVOLUTION OF HIERARCHICAL PLANETARY SYSTEMS Man Hoi Lee and S. J. Peale
The Astrophysical Journal, 59:0 6, 003 August # 003. The American Astronomical Society. All rights reserved. Printed in U.S.A. SECULAR EVOLUTION OF HIERARCHICAL PLANETARY SYSTEMS Man Hoi Lee and S. J.
More informationTheory of mean motion resonances.
Theory of mean motion resonances. Mean motion resonances are ubiquitous in space. They can be found between planets and asteroids, planets and rings in gaseous disks or satellites and planetary rings.
More informationarxiv:astro-ph/ v1 14 Sep 2005
For publication in Bayesian Inference and Maximum Entropy Methods, San Jose 25, K. H. Knuth, A. E. Abbas, R. D. Morris, J. P. Castle (eds.), AIP Conference Proceeding A Bayesian Analysis of Extrasolar
More informationA REGION VOID OF IRREGULAR SATELLITES AROUND JUPITER
The Astronomical Journal, 136:909 918, 2008 September c 2008. The American Astronomical Society. All rights reserved. Printed in the U.S.A. doi:10.1088/0004-6256/136/3/909 A REGION VOID OF IRREGULAR SATELLITES
More informationCelestial Mechanics III. Time and reference frames Orbital elements Calculation of ephemerides Orbit determination
Celestial Mechanics III Time and reference frames Orbital elements Calculation of ephemerides Orbit determination Orbital position versus time: The choice of units Gravitational constant: SI units ([m],[kg],[s])
More informationStability Dynamics Habitability of Planets in Binary Systems
Stability Dynamics Habitability of Planets in Binary Systems Facts and Questions: Binary and Multi-Star Systems in the Solar neighbourhood: >60% Do the gravitational perturbations allow planetary formation?
More informationarxiv: v2 [astro-ph.ep] 30 Nov 2013
Extreme orbital evolution from hierarchical secular coupling of two giant planets Jean Teyssandier,, Smadar Naoz 2,3, Ian Lizarraga 4, Frederic A. Rasio 3,5 arxiv:3.548v2 [astro-ph.ep] 3 Nov 23 ABSTRACT
More informationResonancias orbitales
www.fisica.edu.uy/ gallardo Facultad de Ciencias Universidad de la República Uruguay Curso Dinámica Orbital 2018 Orbital motion Oscillating planetary orbits Orbital Resonances Commensurability between
More informationPlanetary system dynamics. Planetary migration Kozai resonance Apsidal resonance and secular theories Mean motion resonances Gravitational scattering
Planetary system dynamics Planetary migration Kozai resonance Apsidal resonance and secular theories Mean motion resonances Gravitational scattering How should the planets of a typical planetary system
More informationRestricted three body problems in the Solar System: simulations
Author:. Facultat de Física, Universitat de Barcelona, Diagonal 645, 0808 Barcelona, Spain. Advisor: Antoni Benseny i Ardiaca. Facultat de Matemàtiques, Universitat de Barcelona, Gran Via de les Corts
More informationLong period comet encounters with the planets: an analytical approach. G. B. Valsecchi IAPS-INAF, Roma, Italy IFAC-CNR, Sesto Fiorentino, Italy
Long period comet encounters with the planets: an analytical approach G. B. Valsecchi IAPS-INAF, Roma, Italy IFAC-CNR, Sesto Fiorentino, Italy The engine of cometary orbital evolution The orbital evolution
More informationThree objects; 2+1 problem
Three objects; 2+1 problem Having conquered the two-body problem, we now set our sights on more objects. In principle, we can treat the gravitational interactions of any number of objects by simply adding
More informationA dwarf planet is a celestial body that (a) is in orbit around the Sun, (b) has sufficient mass for its selfgravity
Celestial Mechanics - 2013-14 1 Carlo Nipoti, Dipartimento di Fisica e Astronomia, Università di Bologna 20/5/2014 6. Planetary systems 6.1 The Solar System 6.1.1 Structure of the Solar System [R05 1.2;
More informationExoplanets: a dynamic field
Exoplanets: a dynamic field Alexander James Mustill Amy Bonsor, Melvyn B. Davies, Boris Gänsicke, Anders Johansen, Dimitri Veras, Eva Villaver The (transiting) exoplanet population Solar System Hot Jupiters:
More informationPlanetary Perturbations on the 2 : 3 Mean Motion Resonance with Neptune
PASJ: Publ. Astron. Soc. Japan 54, 493 499, 2002 June 25 c 2002. Astronomical Society of Japan. Planetary Perturbations on the 2 : 3 Mean Motion Resonance with Neptune Tetsuharu FUSE Subaru Telescope,
More informationLecture 20: Planet formation II. Clues from Exoplanets
Lecture 20: Planet formation II. Clues from Exoplanets 1 Outline Definition of a planet Properties of exoplanets Formation models for exoplanets gravitational instability model core accretion scenario
More informationarxiv: v1 [astro-ph.ep] 15 Apr 2016
An Earth-Like Planet in GJ 832 System S. Satyal 1, J. Griffith 1, Z. E. Musielak 1,2 1 The Department of Physics, University of Texas at Arlington, Arlington TX 76019 2 Kiepenheuer-Institut für Sonnenphysik,
More informationSecular Evolution of Extrasolar Planetary Systems:
Secular Evolution of Extrasolar Planetary Systems: an Extension of the Laplace-Lagrange Secular Theory Marco Sansottera [a] [a] Namur Center for Complex Systems (naxys) Based on a research work in collaboration
More informationarxiv:astro-ph/ v1 22 Sep 2005
Dynamical Stability and Habitability of γ Cephei Binary-Planetary System Nader Haghighipour arxiv:astro-ph/0509659 v1 22 Sep 2005 Institute for Astronomy and NASA Astrobiology Institute, University of
More informationORBITAL RESONANCES IN PLANETARY SYSTEMS
ORBITAL RESONANCES IN PLANETARY SYSTEMS Renu Malhotra Lunar & Planetary Laboratory, The University of Arizona, Tucson, AZ, USA Keywords: planets, orbits, solar system, Hamiltonian, action-angle variables,
More informationAST111 PROBLEM SET 2 SOLUTIONS. RA=02h23m35.65s, DEC=+25d18m42.3s (Epoch J2000).
AST111 PROBLEM SET 2 SOLUTIONS Home work problems 1. Angles on the sky and asteroid motion An asteroid is observed at two different times. The asteroid is located at RA=02h23m35.65s, DEC=+25d18m42.3s (Epoch
More informationas the orbits of distant planetoids are expected to be randomized over billions of year by the gravity of the four giant planets.
Dynamics of Extreme Outer Solar System Objects Principle investigators/supervisors: Denis Gaidashev (Matematiska Institution); Erik Zackrisson (Institutionen fo r fysik och astronomi) 1. Purpose and aims
More informationRESONANCE OVERLAP IS RESPONSIBLE FOR EJECTING PLANETS IN BINARY SYSTEMS
ACCEPTED FOR PUBLICATION IN APJ Preprint typeset using LATEX style emulateapj v. 11/26/04 RESONANCE OVERLAP IS RESPONSIBLE FOR EJECTING PLANETS IN BINARY SYSTEMS LAWRENCE R. MUDRYK AND YANQIN WU Department
More informationarxiv: v1 [astro-ph] 30 May 2007
Submitted to ApJ Preprint typeset using L A TEX style emulateapj v. 10/09/06 SHRINKING BINARY AND PLANETARY ORBITS BY KOZAI CYCLES WITH TIDAL FRICTION Daniel Fabrycky 1 and Scott Tremaine 1,2 Submitted
More informationarxiv: v1 [astro-ph.ep] 11 Aug 2016
Noname manuscript No. (will be inserted by the editor) Secular and tidal evolution of circumbinary systems. Alexandre C. M. Correia Gwenaël Boué Jacques Laskar arxiv:168.3484v1 [astro-ph.ep] 11 Aug 216
More informationCelestial Mechanics I. Introduction Kepler s Laws
Celestial Mechanics I Introduction Kepler s Laws Goals of the Course The student will be able to provide a detailed account of fundamental celestial mechanics The student will learn to perform detailed
More informationObservational Astronomy - Lecture 4 Orbits, Motions, Kepler s and Newton s Laws
Observational Astronomy - Lecture 4 Orbits, Motions, Kepler s and Newton s Laws Craig Lage New York University - Department of Physics craig.lage@nyu.edu February 24, 2014 1 / 21 Tycho Brahe s Equatorial
More informationThe two-body Kepler problem
The two-body Kepler problem set center of mass at the origin (X = 0) ignore all multipole moments (spherical bodies or point masses) define r := r 1 r 2,r:= r,m:= m 1 + m 2,µ:= m 1 m 2 /m reduces to effective
More informationEXTENDING NACOZY S APPROACH TO CORRECT ALL ORBITAL ELEMENTS FOR EACH OF MULTIPLE BODIES
The Astrophysical Journal, 687:1294Y1302, 2008 November 10 # 2008. The American Astronomical Society. All rights reserved. Printed in U.S.A. EXTENDING NACOZY S APPROACH TO CORRECT ALL ORBITAL ELEMENTS
More informationHABITABLE EXTRASOLAR PLANETARY SYSTEMS, THE CASE OF 55 CNC
HABITABLE EXTRASOLAR PLANETARY SYSTEMS, THE CASE OF 55 CNC Desiree Cotto-Figueroa University of Puerto Rico at Humacao Institute for Astronomy, University of Hawaii Mentor : Nader Haghighipour ABSTRACT
More informationSymplectic Correctors for Canonical Heliocentric N-Body Maps
Symplectic Correctors for Canonical Heliocentric N-Body Maps J. Wisdom Massachusetts Institute of Technology, Cambridge, MA 02139 wisdom@poincare.mit.edu Received ; accepted 2 ABSTRACT Symplectic correctors
More informationarxiv: v2 [astro-ph] 18 Dec 2008
Planet influence on the shape of the hosting star - ellipsoidal variations of tau Bootis W. Dimitrov Astronomical Observatory of Adam Mickiewicz University ul. S loneczna 36, 60-286 Poznań, Poland dimitrov@amu.edu.pl
More informationHD10647 and the Distribution of Exoplanet Properties with Semi-major Axis
Extrasolar planets : Today and Tomorrow ASP Conference Series, Vol. 321, 2004 J.-P. Beaulieu, A. Lecavelier des Etangs, and C. Terquem HD10647 and the Distribution of Exoplanet Properties with Semi-major
More informationPADEU. Pulsating zero velocity surfaces and capture in the elliptic restricted three-body problem. 1 Introduction
PADEU PADEU 15, 221 (2005) ISBN 963 463 557 c Published by the Astron. Dept. of the Eötvös Univ. Pulsating zero velocity surfaces and capture in the elliptic restricted three-body problem F. Szenkovits
More informationThe planet search program at the Coudé Echelle. spectrometer
A&A 374, 675 681 (2001) DOI: 10.1051/0004-6361:20010723 c ESO 2001 Astronomy & Astrophysics The planet search program at the ESO Coudé Echelle spectrometer II. The α Centauri system: Limits for planetary
More informationSpectroscopic orbits of 10 nearby solar-type dwarfs
doi:10.1093/mnras/stu743 Spectroscopic orbits of 10 nearby solar-type dwarfs N. A. Gorynya 1,2 and A. Tokovinin 3 1 Institute of Astronomy of Russian Academy of Science, 48 Pyatnitskaya Str, 109017 Moscow,
More informationA Kozai-resonating Earth quasi-satellite
MNRASL 437, L85 L89 (2014) Advance Access publication 2013 November 13 doi:10.1093/mnrasl/slt147 A Kozai-resonating Earth quasi-satellite M. Connors Athabasca University Observatories, Athabasca AB T9S
More informationISIMA lectures on celestial mechanics. 1
ISIMA lectures on celestial mechanics. 1 Scott Tremaine, Institute for Advanced Study July 2014 The roots of solar system dynamics can be traced to two fundamental discoveries by Isaac Newton: first, that
More informationarxiv: v1 [astro-ph.sr] 11 Sep 2014
18th Cambridge Workshop on Cool Stars, Stellar Systems, and the Sun Proceedings of Lowell Observatory (9-13 June 2014) Edited by G. van Belle & H. Harris BinHab: A Numerical Tool for the Calculation of
More informationThe Heliocentric Model of Copernicus
Celestial Mechanics The Heliocentric Model of Copernicus Sun at the center and planets (including Earth) orbiting along circles. inferior planets - planets closer to Sun than Earth - Mercury, Venus superior
More informationThe Hill stability of the possible moons of extrasolar planets
Mon. Not. R. Astron. Soc. 406, 98 94 00) doi:0./j.65-966.00.6796.x The Hill stability of the possible moons of extrasolar planets J. R. Donnison Astronomy Unit, School of Mathematical Sciences, Queen Mary,
More informationFundamentals of Astrodynamics and Applications
Fundamentals of Astrodynamics and Applications Third Edition David A. Vallado with technical contributions by Wayne D. McClain Space Technology Library Published Jointly by Microcosm Press Hawthorne, CA
More informationSatellite meteorology
GPHS 422 Satellite meteorology GPHS 422 Satellite meteorology Lecture 1 6 July 2012 Course outline 2012 2 Course outline 2012 - continued 10:00 to 12:00 3 Course outline 2012 - continued 4 Some reading
More informationORBITAL RESONANCES IN PLANETARY SYSTEMS 1
ORBITAL RESONANCES IN PLANETARY SYSTEMS 1 Renu Malhotra Lunar & Planetary Laboratory, The University of Arizona, Tucson, AZ, USA Keywords. planets, orbits, solar system, Hamiltonian, action-angle variables,
More informationarxiv: v1 [astro-ph] 28 May 2008
Submitted to the The Astrophysical Journal Letters Possibility of Detecting Moons of Pulsar Planets Through Time-of-Arrival Analysis arxiv:0805.4263v1 [astro-ph 28 May 2008 Karen M. Lewis School of Mathematical
More informationOrbit Characteristics
Orbit Characteristics We have shown that the in the two body problem, the orbit of the satellite about the primary (or vice-versa) is a conic section, with the primary located at the focus of the conic
More informationLecture 11: Binary stars
Lecture 11: Binary stars Senior Astrophysics 2018-04-11 Senior Astrophysics Lecture 11: Binary stars 2018-04-11 1 / 31 Outline 1 Observed characteristics 2 Interacting binaries 3 Gravity in a rotating
More informationAy 1 Lecture 2. Starting the Exploration
Ay 1 Lecture 2 Starting the Exploration 2.1 Distances and Scales Some Commonly Used Units Distance: Astronomical unit: the distance from the Earth to the Sun, 1 au = 1.496 10 13 cm ~ 1.5 10 13 cm Light
More informationarxiv: v2 [astro-ph.sr] 3 Jan 2018
Mon. Not. R. Astron. Soc. 000, 1 7 (2018) Printed 8 January 2018 (MN LATEX style file v2.2) Spectroscopic orbits of nearby solar-type dwarfs. II. arxiv:1712.06389v2 [astro-ph.sr] 3 Jan 2018 N. A. Gorynya
More informationObservational Cosmology Journal Club
Observational Cosmology Journal Club 07/09/2018 Shijie Wang 1. Heller, R. (2018). Formation of hot Jupiters through disk migration and evolving stellar tides. Retrieved from arxiv.1806.06601 2. Rey, J.,
More informationGravitation and the Waltz of the Planets
Gravitation and the Waltz of the Planets Chapter Four Guiding Questions 1. How did ancient astronomers explain the motions of the planets? 2. Why did Copernicus think that the Earth and the other planets
More informationGravitation and the Waltz of the Planets. Chapter Four
Gravitation and the Waltz of the Planets Chapter Four Guiding Questions 1. How did ancient astronomers explain the motions of the planets? 2. Why did Copernicus think that the Earth and the other planets
More information