Planetary Perturbations on the 2 : 3 Mean Motion Resonance with Neptune

Size: px
Start display at page:

Download "Planetary Perturbations on the 2 : 3 Mean Motion Resonance with Neptune"

Transcription

1 PASJ: Publ. Astron. Soc. Japan 54, , 2002 June 25 c Astronomical Society of Japan. Planetary Perturbations on the 2 : 3 Mean Motion Resonance with Neptune Tetsuharu FUSE Subaru Telescope, National Astronomical Observatory of Japan, 650 North A ohoku Place, Hilo, HI 96720, U.S.A. tetsu@subaru.naoj.org (Received 2001 December 31; accepted 2002 April 2) Abstract A mean motion resonance with Neptune is one of the major dynamical characteristics of the Edgeworth Kuiper Belt Objects that spread over the outside of Neptune s orbit. We investigated the influence of Jupiter, Saturn, and Uranus on the dynamical structure of the 2 : 3 mean motion resonance with Neptune. Since the distance between the planets and objects has great influence on the dynamical structure, numerical integrations were carried out for test particles for 5 Myr under five fictitious solar-system models having different configurations. In discussion about the results, we concentrate on the size of the EKBO eccentricities, e, because they have an especially close relation to the orbital evolution; for example, larger e objects have more opportunities to encounter other objects or planets. By comparing the results of each model, it was found that the resonance structure is not largely affected by Jupiter and Saturn, and that the largest e of the resonant objects depends on the distance between Uranus and the resonance. The distance also influences the structure around the resonance because of variation in the secular resonance positions. We further found that the largest e values of the non-resonant objects at each semi-major axis do not depend on the distance from Uranus; however, they depend on the distance from Neptune. Key words: celestial mechanics Kuiper Belt solar system: formation 1. Introduction As Edgeworth (1949) and Kuiper (1951) predicted, approximately 500 objects, called Edgeworth Kuiper Belt Objects (EKBOs), have been discovered in a region outside of Neptune s orbit (called the Edgeworth Kuiper Belt; EKB) since August of 1992 (Jewitt, Luu 1993). The EKBOs are thought to be remnants of planetesimals in the early stage of the solar system formation as well as the origin of short-period comets (Everhart 1972; Fernández 1980). Some of the discovered EKBOs show two major dynamical characteristics; a mean motion resonance with Neptune is the first example. When the ratio of the orbital periods between Neptune and an EKBO can be expressed as (p + q):p, where p and q are integers, the EKBO will exhibit the (p + q) :p mean motion resonance with Neptune, and a resonant angle, called a critical argument, σ, librates around a fixed angle. The definition of σ is σ (p + q)λ N + pλ + qϖ =(p + q)[ l N +(ϖ ϖ N )] + pl, (1) where λ, ϖ,andl are the mean longitude, longitude of the perihelion, and mean anomaly, respectively, and N indicates Neptune. Among the discovered EKBOs, 39 multi-opposition EKBOs have been found in the 2 : 3 mean motion resonance with Neptune (Marsden 2002). They are called Plutinos because the motion is like Pluto and σ librates around 180.We henceforth designate a mean motion resonance with Neptune as a resonance. A secular resonance is also an important dynamical characteristic of EKBOs. When there is a commensurability between the frequencies of the precession rates of either ϖ (a longitude of perihelion) or Ω (a longitude of ascending node) among two objects, a secular resonance occurs. A ν 8 secular resonance takes place if the frequencies of the precession rates of an EKBO s and Neptune s ϖ are commensurable. When an EKBO exhibits the ν 8 secular resonance, the eccentricity, e, of the EKBO is excited. On the other hand, the ν 18 secular resonance occurs when the frequencies of the precession rates of an EKBO s and Neptune s Ω are resonant, and the EKBO s inclination, i, is excited. Knežević et al. (1991) analytically investigated the locations of secular resonances between 2 and 50 AU from the Sun, and found that a few secular resonances are present in the EKB. Gallardo and Ferraz-Mello (1997) discussed the σ libration of 2 : 3, 3 : 5 and 1 : 2 resonant objects based on a time frequency analysis. It is understood that the 2 : 3 resonance is complex because of the ν 8 and ν 18 secular resonances and the Kozai resonance (Morbidelli et al. 1995; Morbidelli 1997); the Kozai resonance couples the evolution of the EKBO s i and e, and keeps the value of a(1 e 2 )cosi constant, where a is the semi-major axis of the EKBO. Carrying out semianalytical and numerical integration methods, Gallardo and Ferraz-Mello (1998) studied stable regions in the 2 : 3 resonance. Yu and Tremaine (1999) also used three-body planar semi-analytical and numerical integration methods for the 2 : 3 resonance, and found that the stability of the resonant orbits is determined by the size of the object s e. NesvornỳandRoig (2000, 2001) investigated the 2 : 3, 1 : 2, and 3 : 4 resonances with numerical integration methods, and showed chaotic and regular motion of the resonant objects. The purpose of the present work was to extensively research into the 2 : 3 resonance and to investigate the influence of perturbing planets on the dynamical structure in and around the 2 : 3 resonance. We carried out numerical integrations of EKBO test objects under five fictitious solar system models

2 494 T. Fuse [Vol. 54, with developed integrator (section 2) and showed variations in the orbital elements of the objects in and around the 2 : 3 resonance (section 3). Comparing the results of the integrations, we concentrate on the eccentricity of EKBOs and discuss the influence of the perturbing planets on the dynamical structure in and around the resonance in section Methods of Analysis Symplectic integrators (SI) are good algorithms because solutions are area-preserving, the discretization error in the energy integral does not have a secular term, and it can treat a highly eccentric orbit without changing the step-size (Forest, Ruth 1990). It is convenient to use mixed-variable symplectic integrators (MVSI) if we can write the Hamiltonian as F = F 0 (p,q) +εr(q), where F 0 is an integrable Hamiltonian and εr is a perturbation with a small parameter ε (Kinoshita et al. 1991; Wisdom, Holman 1991). This is because we can use a longer step-size in MVSI than that in SI. We generally have to choose a step-size which is adequate for the smallest orbital period when N-body numerical integrations are carried out. Since we know that the orbital period of EKBOs is the longest in the Sun planets EKBO system, we developed MVSI to MVSI with Files (MVSI F ; see Appendix). The positions of the perturbing planets were calculated in advance with a small step-size, and we saved the data into a file. We then integrated EKBO test objects using the file as the perturbing planet positions. MVSI F can include Jupiter, Saturn, Uranus, and Neptune as the perturbing planets (the masses of the four inner planets are added to that of the Sun). We did not consider galactic tides and a passing star s perturbation, because they finally begin to have dynamical effects on objects whose a is 1000 AU (Duncan et al. 1987). The step-size was 300 days to make the planet positions file and 1500 days for the EKBO test objects. Integrations using a step-size of 300 days for the test objects were also carried out, and it was confirmed that there is no global difference between 1500 and 300 days and CPU time can be saved several dozen percent. The calculation span of the numerical integrations was 5 Myr, which is almost equal to the long-period variations in the EKBO s e and ϖ. The error of the Neptune s longitude was estimated to be 5 for 5 Myr. The initial positions, velocities, and masses of the planets were obtained from JPL DE403 planetary ephemerides at the epoch of JD Since MVSI F could not treat a close approach between the test objects and planets, the calculation was stopped when the test objects entered any active spheres of the planets given in Danby (1988). 3. Results Not only the existence of perturbing planets, but also the distance between a planet and an EKBO is important to study the influence on the orbital stability of the objects. Since Uranus is especially the closest planet to the resonant EKBOs with Neptune, the distance between them should have an effect on the dynamical structure of the EKBOs. Therefore, we integrated EKBO test objects placed in and around the 2:3 resonance for 5 Myr assuming five fictitious solar system models: [1] the perturbing planets are Jupiter, Saturn, Uranus, and Neptune; the semi-major axis of Uranus a U is obtained from DE403 ( a U0 ) : Model J N; [2] the perturbing planets are Uranus and Neptune, and a U = a U0 : Model U + N; [3] the perturbing planets are Uranus and Neptune, and a U = a U0 3AU : Model U +N; [4] the perturbing planets are Uranus and Neptune, and a U = a U0 + 3AU : Model U + +N;and [5] the perturbing planet is Neptune only : Model N. The initial conditions of the test objects were as follows: AU a AU ( a =0.01 AU), 0.01 e ( e =0.0025), i =1, which refer to invariable planes consisting of the Sun four planets for [1], of the Sun Uranus Neptune for [2], [3], and [4], and of the Sun Neptune for [5], ϖ ϖ N =0, Ω = Ω N +90,andσ = 180. The two initial secular angles were chosen from the viewpoint of close pericenters and far nodes. The initial a and e sets are shown on an a e plane in figure 1. We also draw the 2 : 3 resonance borders using the semi-analytical methods developed based on Yoshikawa (1990, 1991). The two arcs crossing the resonance borders indicate Neptune-crossing and Uranus-crossing curves. [1] Model J N We show the initial a and e values of the test objects that survive for 5 Myr in figure 2. The colors of σ stand for the amplitude of σ,definedas σ σ 180, (2) because σ of the 2 : 3 resonant objects librates around 180, like Pluto. It is clear that the resonant σ distribution forms a layer structure at the center of a =39.3AU, and that the inner region has a smaller amplitude than the outer. This is consistent with Gallardo and Ferraz-Mello (1998) as well as Nesvornỳ and Roig (2000). The smallest value of the resonant σ is approximately 6, which means that the amplitude of σ is 174 σ 186. The total number of the surviving objects is 8102 compared with the initial objects. We can see difference between the Neptune crossing-curve and the largest non-resonant initial e, which is found to be approximately Gladman et al. (2001) also discuss the boundary between the stable and chaotic orbits on the initial a e plane. They empirically imply that it is q = (a 30), where q is the object s perihelion distance. Because of the relation q = a(1 e), the difference between the Neptunecrossing curve and the stable maximum e was calculated at around the 2 : 3 resonance. The value in figure 2 supports the empirical rule by Gladman et al. (2001). In order to investigate the dynamical structure of the resonance on the a e plane, we give the following time-average values of a and e: N / N / a = a i N, e = e i N, (3) i=1 i=1 where N stands for the number of output steps for the integration time. The colors indicate σ defined by equation (2).

3 No. 3] The 2 : 3 Mean Motion Resonance with Neptune Fig. 1. Initial a and e conditions on the a e plane in and around the 2 : 3 resonance for numerical integrations. a is 0.01 AU and e is The arcs show the Neptune-crossing and Uranus-crossing curves. Fig. 2. Initial a and e of test objects that survived for 5 Myr for [1] Model J N. The colors of σ denote σ = σ 180. Fig. 3. Time-averaged a and e of test objects that survived for 5 Myr for [1] Model J N. The colors of σ denote σ = σ 180. The test objects in the 2 : 3 resonance form a line. Fig. 4. Initial a and e of test objects which survived for 5 Myr for [2] Model U + N. The colors of σ denote σ = σ 180. Fig. 5. Time-averaged a and e of test objects which survived for 5 Myr for [2] Model U + N. The colors of σ denote σ = σ 180. We can see that the a and e values in figure 3 form a line at a AU and this is the 2 : 3 resonance center. Since figure 2 shows the osculating orbital elements for the initial conditions, there is a difference in the positions of the layer and line centers in figures 2 and 3. The maximum resonant e at the center is approximately 0.40, and the value is approximately 1.67-times larger than the Neptune-crossing curve. This is because the ν8 secular resonance is not around the resonance center (Morbidelli 1997), and the 2 : 3 resonant objects can be stable at a large-e region for a long time. Gallardo and Ferraz-Mello (1998) reported that the maximum e is 0.37 in their integrations. The slight difference might be due to the different initial conditions. According to a list of EKBO orbital elements at the website of the Minor Planet Center ( on 2002/2/24), the largest e of the 2 : 3 resonant EKBO is of 1996 TP66. It is consistent with the results in figure 3. Furthermore, several lines are, for example, seen at a = 38.20, 38.48, 39.98, 40.42, and AU. According to the Kepler s third law, they correspond to the 7 : 10 (38.20 AU), 9 : 13 (38.48 AU), 9 : 14 (40.42 AU), and 7 : 11 (40.70 AU) resonances with Neptune, and the 1 : 3 (39.98 AU) 495

4 496 T. Fuse [Vol. 54, resonance with Uranus. The minimum e values around 40 AU are larger than those in other regions, because e is excited by the ν 8 secular resonance pointed out by Knežević et al. (1991) and Morbidelli (1997). Most of the largest e values of nonresonant objects (blue marks) are approximately 0.05 smaller than the Neptune-crossing curve for any a-region. This means surviving objects that have a larger e than the curve must be in the 2 : 3 resonance. On the other hand, the value is inconsistent with figure 2 and Gladman et al. (2001), because of the difference between the initial conditions and time average values. When we discuss a more realistic configuration of the outer solar system, we should use the value of 0.05 led by figure 3. The i values of the test objects were also calculated, and the maximum and minimum i were found to be 7. 5and 0. 5, respectively. In the left region of the resonance borders, most i values are smaller than approximately 3,andthe i values are between 1 and 5 in the right region of the resonance borders (around 40 AU). The phenomena are consistent with the presence of the ν 18 secular resonance (Knežević et al. 1991; Morbidelli 1997). Because of the existence of the ν 8 and ν 18 secular resonances outside the 2 : 3 resonance borders, the small-e and -i EKBOs should not be there in a real EKB. [2] Model U +N Figure 4 describes the initial a and e values of objects that survive during integrations under the Sun Uranus Neptune system. We can see that the shape of the resonant σ structure in figure 4 is fairly similar to that in figure 2. The total surviving objects are 8357, which is also close to the results of the model [1]. The a and e values under the model [2] are shown in figure 5. The maximum resonant e is approximately 0.39 at a AU, where the a and e values form a line. The difference between this value and that of the model [1] might be quite small, and is approximately We can say from these similarities that the influence of Jupiter and Saturn on the dynamical structure of the 2 : 3 resonance on the a e plane is small. It is seen that the e values outside the resonance borders are fairly larger than those of the model [1]. The difference should be due to the change in the positions of the ν 8 secular resonance because of the absence of Jupiter and Saturn, or the direct effects of the two planets on the objects. We need a further investigation to conclude that the difference is due to the direct or indirect influence of Jupiter and Saturn. On the other hand, most maximum e values of the non-resonant objects are approximately 0.05 smaller than the Neptune-crossing curve at any a values. Since this model does not include Jupiter and Saturn, it seems that the influence of the two planets on these values is also small. [3] Model U +N The next model contains two perturbing planets, the same as the model [2], but a U is by 3 AU smaller than a U0. It is expected that the effect of Uranus on the objects becomes weak compared with the model [2], and the resonance structure on the a e plane may be different from those of the two models, [1] and [2]. In figure 6, we indicate the initial a and e of the test objects that survive for 5 Myr. The Uranus-crossing curve disappears from the figure because of the change in a U. According to figure 6, the resonant σ distribution is fairly different from those in figures 2 and 4; i.e. more resonant objects that initially have large-e remain stably in the resonance. It was found that objects survive, and the number is approximately 1.3-times larger than those of models [1] and [2]. The results reflect the weakness of the influence of Uranus. We can see the a and e values in figure 7. The maximum resonant e is approximately 0.42 at a AU, which is obviously larger than those of models [1] and [2]. Because of the distance from Uranus, the maximum resonant e should be large. On the other hand, most objects outside the resonance borders stay near the initial positions, expect for some resonant objects. These show that the distance between Uranus and the resonance must be important for determining the largest-e value of the resonant objects and the dynamical structure around the 2 : 3 resonance. The largest e values of the non-resonant objects are approximately 0.05 smaller than the Neptune-crossing curve, which is the same as those in models [1] and [2]. From the results, we can say that the values are not determined by the distance from Uranus, but from Neptune. In figure 6, objects forming islands are observed around the initial a =38.20 AU and e = ; lines can also be seen near the same positions in figure 7. Such objects have far larger e values than the Neptune-crossing curve; however, they are stable because of being in the 7 : 10 resonance with Neptune. [4] Model U + +N Although Uranus and Neptune are included in this model in the same way as in models [2] and [3], a U is 3 AU larger than a U0. It is possible that we will see a large effect on the structure on the a e plane in and around the 2 : 3 resonance Figure 8 contains the initial a and e values of the objects that survive during the integrations. Since a U is larger, the Uranuscrossing curve comes down in figure 8. The smaller e region of the resonant σ structure around a 39.4AU is similar to those of models [1], [2], and [3] in figures 2, 4, and 6, respectively. It is, however, seen that the maximum initial e value of the surviving resonant objects is significantly smaller than those of other models. The total number of the surviving objects is 4856, which is the minimum value, and the smallest distance between Uranus and the resonance must make the differences. We show a and e on the a e plane in figure 9. The maximum resonant e has the smallest value among the previous four models, and is approximately 0.29 at a 39.45AU. The structure outside the resonance borders is entirely different from the others, and there are especially few surviving objects at the small-e region at 38AU a 39AU. Comparing these figures, we can say that the existence of Uranus determines the structure in and around the 2 : 3 resonance on the a e plane. Most of the values of the maximum non-resonant e are, in contrast, the same as those in models [1] [3] and it is approximately It becomes clear that it does not depend on the distance from Uranus.

5 No. 3] The 2 : 3 Mean Motion Resonance with Neptune 497 Fig. 7. Time-averaged a and e of test objects that survived for 5 Myr for [3] Model U + N. The colors of σ denote σ = σ 180. Fig. 8. Initial a and e of test objects that survived for 5 Myr for [4] Model U + + N. The colors of σ denote σ = σ 180. Fig. 9. Time-averaged a and e of test objects that survived for 5 Myr for [4] Model U + + N. The colors of σ denote σ = σ 180. Fig. 10. Initial a and e of test objects that survived for 5 Myr for [5] Model N. The colors of σ denote σ = σ 180. Fig. 11. Time-averaged a and e of the test objects that survived for 5 Myr for [5] Model N. The colors of σ denote σ = σ 180. Fig. 6. Initial a and e of test objects that survived for 5 Myr for [3] Model U + N. The colors of σ denote σ = σ 180.

6 498 T. Fuse [Vol. 54, Table 1. Surviving object ratios and maximum resonant e values. Model Survived (%) Max e [1] Model J N [2] Model U+N [3] Model U + N [4] Model U + + N [5] Model N Note. Given are the ratios of the surviving object numbers to the initial object numbers, and the maximum resonant e values for the numerical integrations under the five fictitious solar system models. [5] Model N This model consists only of Neptune as the perturbing planet. Because the test object is massless, a system comprising the Sun and Neptune is a two-body problem and the orbital elements of Neptune do not change, except for the mean anomaly, l N. Figure 10 indicates the initial a and e of the objects that survive for 5 Myr; we can see that the resonant σ layer structure is the same as in the other models. The smallest resonant σ is approximately Since Uranus is not included in the system, the motion of the objects is not affected by Uranus, and the maximum initial e of the surviving objects is the largest initial value of ; however, a further investigation is required to confirm the real stable maximum e in this model. The total number of surviving objects is 13045, which is the maximum value among the five fictitious solar system models. The largest resonant e is approximately 0.47 at a AU in figure 11, which is also the maximum value among the five models. It can be seen that the resonance in this model is quite stable, and the resonant objects will be able to have large e values. On the other hand, the values between the maximum non-resonant e and the Neptune-crossing curve is approximately 0.05; we should notice that this value is the same as that for the other models. Therefore, we conclude that Jupiter, Saturn, and Uranus do not have any influence on the largest e values of non-resonant objects; however, the distance from Neptune determines them. As is the case with figures 6 and 7, we can see some islands and lines above the Neptune-crossing curve. They correspond to the 7:10 (38.20 AU), 9:13 (38.48 AU), and 7:11 (40.70 AU) resonances with Neptune; also, the resonant objects are stable even if they have large e values. 4. Discussion We summarize the results of our numerical integrations under the five solar system models for the 2 : 3 resonance in table 1. Comparing the maximum resonant e in table 1, we can say that the values are determined by the distance between Uranus and the resonance; for example, a further distance from Uranus makes the resonant maximum e larger. We also tabulate the ratio of the surviving object number to the initial number (= 50960) in table 1. It can be seen that the distance Fig. 12. Distribution of a and e of the discovered EKBOs between 37.8 AU and 41 AU. The orbital elements are based on the website of the Minor Planet Center on 2002/2/24. between Uranus and the resonance has a close relation to the ratio; if it is larger, the ratio becomes smaller. There is no large difference in both the resonance structure and the surviving number ratio between [1] Model J N and [2] Model U + N. This indicates that the influence of Jupiter and Saturn on the resonance structure is small, and we need a further investigation to confirm that it has direct or indirect effects of the planets. We cannot see the dependence of the distance from Uranus on the largest e values of the non-resonant objects; however, it is found that the whole structure outside the resonance would also be affected by Uranus according to the results of [2] Model U + N, [3] Model U + N, and [4] Model U + +N. This is because the motion of Neptune should be different if Uranus s orbit differs, and the positions of the ν 8 secular resonance will change. Therefore, it is concluded that the structure around the resonance is determined by the distance between Uranus and the resonance, while the maximum e values of the nonresonant objects are restricted by the distance from Neptune. It is thought that planets may have migrated outward during the solar system formation (Malhotra 1995). From the above discussion, if Uranus migrated during the solar system formation, the structure of the resonance was affected by the variation in the distance between Uranus and the resonance. This means that the evolution of Uranus will be determined if the resonance structure is revealed. According to the orbital distribution of the observed EKBOs between AU in figure 12, we can presume that the 2 : 3 resonant EKBOs having larger e values should exist, and that non-resonant small-e EKBOs ought to be around both 38.7 AU and 39.8 AU if Uranus did not migrate. Further survey observations for those undiscovered EKBOs are expected to clarify the real dynamical structure in and around the resonance and to understand the possibility of the migration of Uranus. 5. Conclusion We carried out numerical integrations for EKBO test objects to study the dynamical structure in and around the 2 : 3 resonances under five fictitious solar system models. We reached the following conclusions based on our investigations:

7 No. 3] The 2 : 3 Mean Motion Resonance with Neptune 499 Comparing the results of [1] Model J N and [2] Model U + N, we can say that the effect of Jupiter and Saturn on the dynamical structure of the 2 : 3 resonance is small. The maximum e fairly depends on the distance between Uranus and the 2 : 3 resonance. The distance between Uranus and the 2 : 3 resonance also has an effect on the dynamical structure around the 2 : 3 resonance, because of the variation in the positions of the ν 8 secular resonance. The values between the largest non-resonant objects e and the Neptune-crossing curve are approximately 0.05 under the five solar system models. They are determined by the distance from Neptune. If Uranus did not migrate, undiscovered larger-e 2:3 resonant and small-e non-resonant EKBOs should exist. The orbital evolution of Uranus will be confirmed when we reveal the dynamical structure in and around the 2 : 3 resonance by further observations. I thank Hiroshi Kinoshita and Hiroshi Nakai for fruitful discussions. I am also grateful Makoto Yoshikawa for useful comments. Workstations at the Astronomical Data Analysis Center of the National Astronomical Observatory of Japan were used for the data analysis. Symplectic Integrator for Hamiltonian includ- Appendix. ing Time We assume a Hamiltonian system, H, that includes time, t, explicitly: H (q,p,t)=k(p) +U(q,t), (A1) where q and p are a canonical set. After introducing T,which is a momentum conjugated to t, and the new Hamiltonian F Danby, J. M. A. 1988, Fundamentals of Celestial Mechanics, (Virginia: Willmann-Bell), 352 Duncan, M., Quinn, T., & Tremaine, S. 1987, AJ, 94, 1330 Edgeworth, K. E. 1949, MNRAS, 109, 600 Everhart, E. 1972, Astrophys. Lett., 10, 131 Fernández, J. A. 1980, MNRAS, 192, 481 Forest, E., & Ruth, R. D. 1990, Physica, 43D, 105 Gallardo, T., & Ferraz-Mello, S. 1997, AJ, 113, 863 Gallardo, T., & Ferraz-Mello, S. 1998, Planet. Space Sci., 46, 945 Gladman, B., Holman, M., Grav, T., Kavelaars, J., Nicholson, P., Aksnes, K., & Petit, J.-M. 2001, astro-ph/ Jewitt, D., & Luu, J. 1993, Nature, 362, 730 Kinoshita, H., Yoshida, H., & Nakai, H. 1991, Celest. Mech. Dyn. Astron., 50, 59 can be written as F = T + H = K(p) +T + U(q,t) = K (p,t )+U(q,t). The second order SI is q 1 = q 0 + c 1 τ K p (p 0,T 0 )=q 0 + c 1 τ K p (p 0), References t 1 = t 0 + c 1 τ K T (p 0,T 0 )=t 0 + c 1 τ, p 1 = p 0 d 1 τ U q (q 1,t 1 ), T 1 = T 0 d 1 τ U t (q 1,t 1 ), (A2) (A3) q 2 = q 1 + c 2 τ K p (p 1)=t 0 + c 1 τ, t 2 = t 1 + c 2 τ = t 0 +(c 1 + c 2 )τ, p 2 = p 1 d 2 τ U (A4) q (q 2,t 2 ), T 2 = T 1 d 2 τ U t (q 2,t 2 ), where τ is the step size, and (q 0,p 0 )and(q 2,p 2 ) are the initial values and numerical solutions after τ, respectively. We note that we do not have to evaluate T 1 and T 2 because they do not actually appear in the system. If we do not calculate the positions of perturbing planets simultaneously, the Hamiltonian includes the time explicitly, as above. The positions of the perturbing planets are, in practice, calculated with a small time step; they are saved into a data file beforehand for integrating for test objects. Knežević, Z., Milani, A., Farinella, P., Froeschlé, Ch., & Froeschlé, Cl. 1991, Icarus, 93, 316 Kuiper, G. P. 1951, Astrophysics, ed. J. A. Hynek (New York: McGraw-Hill), 357 Malhotra, R. 1995, AJ, 110, 420 Marsden, B. G. 2002, Minor Planet Electron. Circ., 2002-C08 Morbidelli, A. 1997, Icarus, 127, 1 Morbidelli, A., Thomas, F., & Moons, M. 1995, Icarus, 118, 322 Nesvornỳ, D., & Roig, F. 2000, Icarus, 148, 282 Nesvornỳ, D., & Roig, F. 2001, Icarus, 150, 104 Wisdom, J., & Holman, M. 1991, AJ, 102, 1528 Yoshikawa, M. 1990, Icarus, 87, 78 Yoshikawa, M. 1991, Icarus, 92, 94 Yu, Q., & Tremaine, S. 1999, AJ, 118, 1873

The orbit evolution of 32 plutinos over 100 million year

The orbit evolution of 32 plutinos over 100 million year A&A 368, 700 705 (2001) DOI: 10.1051/0004-6361:20010056 c ESO 2001 Astronomy & Astrophysics The orbit evolution of 32 plutinos over 100 million year X.-S. Wan 1 and T.-Y. Huang 1,2 1 Astronomy Department,

More information

A SYMPLECTIC MAPPING MODEL FOR THE STUDY OF 2:3 RESONANT TRANS-NEPTUNIAN MOTION

A 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 information

FROM THE SCATTERED DISK TO THE OORT CLOUD The Extended Scattered Disk

FROM THE SCATTERED DISK TO THE OORT CLOUD The Extended Scattered Disk FROM THE SCATTERED DISK TO THE OORT CLOUD The Extended Scattered Disk Julio A. Fernández Departamento de Astronomía, Facultad de Ciencias, Montevideo, URUGUAY Adrián Brunini, Tabaré Gallardo, Rodney Gomes

More information

ON THE ORIGIN OF THE HIGH-PERIHELION SCATTERED DISK: THE ROLE OF THE KOZAI MECHANISM AND MEAN MOTION RESONANCES

ON THE ORIGIN OF THE HIGH-PERIHELION SCATTERED DISK: THE ROLE OF THE KOZAI MECHANISM AND MEAN MOTION RESONANCES ON THE ORIGIN OF THE HIGH-PERIHELION SCATTERED DISK: THE ROLE OF THE KOZAI MECHANISM AND MEAN MOTION RESONANCES RODNEY S. GOMES 1, TABARE GALLARDO 2, JULIO A. FERNA NDEZ 2 AND ADRIA N BRUNINI 3 1 GEA/OV/UFRJ

More information

THE PLANE OF THE KUIPER BELT

THE PLANE OF THE KUIPER BELT The Astronomical Journal, 127:2418 2423, 2004 April # 2004. The American Astronomical Society. All rights reserved. Printed in U.S.A. THE PLANE OF THE KUIPER BELT Michael E. Brown Division of Geological

More information

Dynamical behaviour of the primitive asteroid belt

Dynamical behaviour of the primitive asteroid belt Mon. Not. R. Astron. Soc. 293, 405 410 (1998) Dynamical behaviour of the primitive asteroid belt Adrián Brunini Observatorio Astronómico de La Plata, Profoeg, Paseo del Bosque, (1900) La Plata, Argentina

More information

Forming the Kuiper Belt by the Outward Transport of Objects During Neptune s Migration

Forming the Kuiper Belt by the Outward Transport of Objects During Neptune s Migration Forming the Kuiper Belt by the Outward Transport of Objects During Neptune s Migration Harold F. Levison and Alessandro Morbidelli Southwest Research Institute, Boulder, CO 80302, USA Observatoire de la

More information

arxiv: v1 [astro-ph.ep] 3 Jul 2018

arxiv: 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 information

The Mars 1:2 Resonant Population. Tabaré Gallardo.

The Mars 1:2 Resonant Population. Tabaré Gallardo. 1 The Mars 1:2 Resonant Population Tabaré Gallardo Departamento de Astronomía, Facultad de Ciencias, Iguá 4225, 11400 Montevideo, Uruguay gallardo@fisica.edu.uy Icarus, in press May 2007 2 Abstract An

More information

EVOLUTIONS OF SMALL BODIES IN OUR SOLAR SYSTEM

EVOLUTIONS 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 information

as the orbits of distant planetoids are expected to be randomized over billions of year by the gravity of the four giant planets.

as 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 information

Dynamic Exoplanets. Alexander James Mustill

Dynamic 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 information

arxiv:astro-ph/ v1 27 Sep 1995

arxiv:astro-ph/ v1 27 Sep 1995 THE PHASE SPACE STRUCTURE NEAR NEPTUNE RESONANCES IN THE KUIPER BELT arxiv:astro-ph/9509141v1 27 Sep 1995 RENU MALHOTRA Lunar and Planetary Institute, 3600 Bay Area Blvd., Houston, TX 77058. Electronic

More information

A REGION VOID OF IRREGULAR SATELLITES AROUND JUPITER

A 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 information

Chaotic diffusion of small bodies in the Solar System

Chaotic diffusion of small bodies in the Solar System Chaotic diffusion of small bodies in the Solar System Kleomenis Tsiganis (tsiganis@obs-nicefr) and Alessandro Morbidelli Observatoire de la Côte d Azur - Nice, CNRS, BP 4229, 06304 Nice Cedex 4, France

More information

LONG-TERM DYNAMICS AND THE ORBITAL INCLINATIONS OF THE CLASSICAL KUIPER BELT OBJECTS Marc J. Kuchner 1. Michael E. Brown. and Matthew Holman

LONG-TERM DYNAMICS AND THE ORBITAL INCLINATIONS OF THE CLASSICAL KUIPER BELT OBJECTS Marc J. Kuchner 1. Michael E. Brown. and Matthew Holman The Astronomical Journal, 124:1221 1230, 2002 August # 2002. The American Astronomical Society. All rights reserved. Printed in U.S.A. LONG-TERM DYNAMICS AND THE ORBITAL INCLINATIONS OF THE CLASSICAL KUIPER

More information

Investigation of evolution of orbits similar to that of (4179) Toutatis

Investigation of evolution of orbits similar to that of (4179) Toutatis Investigation of evolution of orbits similar to that of (4179) Toutatis Wm. Robert Johnston May 1994 Abstract This project compared the orbital evolution of asteroid (4179) Toutatis to that of several

More information

Resonance In the Solar System

Resonance In the Solar System Resonance In the Solar System Steve Bache UNC Wilmington Dept. of Physics and Physical Oceanography Advisor : Dr. Russ Herman Spring 2012 Goal numerically investigate the dynamics of the asteroid belt

More information

Long simulations of the Solar System: Brouwer s Law and chaos

Long simulations of the Solar System: Brouwer s Law and chaos ANZIAM J. 46 (E) pp.c1086 C1103, 2005 C1086 Long simulations of the Solar System: Brouwer s Law and chaos K. R. Grazier W. I. Newman James M. Hyman Philip. W. Sharp (received 28 October 2004, revised 18

More information

Resonance Capture. Alice Quillen. University of Rochester. Isaac Newton Institute Dec 2009

Resonance Capture. Alice Quillen. University of Rochester. Isaac Newton Institute Dec 2009 Resonance Capture Alice Quillen University of Rochester Isaac Newton Institute Dec 2009 This Talk Resonance Capture extension to nonadiabatic regime Astrophysical settings: Dust spiraling inward via radiation

More information

Dynamical evolution of asteroid fragments originating near the ν 6 resonance

Dynamical evolution of asteroid fragments originating near the ν 6 resonance Dynamical evolution of asteroid fragments originating near the ν 6 resonance Takashi Ito and Renu Malhotra Lunar & Planetary Laboratory, The University of Arizona, Tucson, AZ 85721 0092, USA Correspondence

More information

Long-term integrations and stability of planetary orbits in our Solar system

Long-term integrations and stability of planetary orbits in our Solar system Mon. Not. R. Astron. Soc. 336, 483 500 (2002) Long-term integrations and stability of planetary orbits in our Solar system Takashi Ito and Kiyotaka Tanikawa Data Analysis Center and Theoretical Astrophysics

More information

Planetary Embryos Never Formed in the Kuiper Belt

Planetary Embryos Never Formed in the Kuiper Belt Icarus 157, 241 248 (2002) doi:10.1006/icar.2002.6832, available online at http://www.idealibrary.com on Planetary Embryos Never Formed in the Kuiper Belt A. Morbidelli and C. Jacob Observatoire de la

More information

Inclination Mixing in the Classical Kuiper Belt

Inclination Mixing in the Classical Kuiper Belt Submitted to ApJ; preprint: Jan 11, 211. Inclination Mixing in the Classical Kuiper Belt Kathryn Volk and Renu Malhotra Lunar and Planetary Laboratory, University of Arizona, Tucson, AZ 85721, USA. kvolk@lpl.arizona.edu

More information

Asteroid 2014 YX 49 : a large transient Trojan of Uranus

Asteroid 2014 YX 49 : a large transient Trojan of Uranus Preprint 6 February 2017 Compiled using MNRAS LATEX style file v3.0 Asteroid 2014 YX 49 : a large transient Trojan of Uranus C. de la Fuente Marcos and R. de la Fuente Marcos Universidad Complutense de

More information

arxiv: v2 [astro-ph.ep] 1 Feb 2010

arxiv: v2 [astro-ph.ep] 1 Feb 2010 Noname manuscript No. (will be inserted by the editor) Where are the Uranus Trojans? R. Dvorak Á. Bazsó L.-Y. Zhou arxiv:911.2366v2 [astro-ph.ep] 1 Feb 21 Received: date / Accepted: date Abstract The area

More information

Dynamical properties of the Solar System. Second Kepler s Law. Dynamics of planetary orbits. ν: true anomaly

Dynamical 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 information

arxiv:astro-ph/ v2 14 Jan 1999

arxiv:astro-ph/ v2 14 Jan 1999 DYNAMICS OF THE KUIPER BELT RENU MALHOTRA Lunar and Planetary Institute arxiv:astro-ph/9901155v2 14 Jan 1999 MARTIN DUNCAN Queen s University and HAROLD LEVISON Southwest Research Institute ABSTRACT Our

More information

Observational Astronomy - Lecture 4 Orbits, Motions, Kepler s and Newton s Laws

Observational 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 information

DYNAMICS OF THE KUIPER BELT

DYNAMICS OF THE KUIPER BELT To appear in Protostars and Planets IV Preprint (Jan 1999): http://www.lpi.usra.edu/science/renu/home.html DYNAMICS OF THE KUIPER BELT Renu Malhotra (1), Martin Duncan (2), and Harold Levison (3) (1) Lunar

More information

Resonancias orbitales

Resonancias 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 information

Corralling a distant planet with extreme resonant Kuiper belt objects

Corralling a distant planet with extreme resonant Kuiper belt objects Corralling a distant planet with extreme resonant Kuiper belt objects Renu Malhotra 1 Kathryn Volk Xianyu Wang 2 arxiv:1603.02196v3 [astro-ph.ep] 2 Jun 2016 Lunar and Planetary Laboratory, The University

More information

Planetenbewegung in Sternsystemen. The Effect of Resonances. Part 2

Planetenbewegung 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 information

Theory of mean motion resonances.

Theory 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 information

A few points on the dynamical evolution of the young solar system. Renu Malhotra The University of Arizona

A few points on the dynamical evolution of the young solar system. Renu Malhotra The University of Arizona A few points on the dynamical evolution of the young solar system Renu Malhotra The University of Arizona A few points on the dynamical evolution of the young solar system Renu Malhotra The University

More information

The MATROS project: Stability of Uranus and Neptune Trojans. The case of 2001 QR322

The MATROS project: Stability of Uranus and Neptune Trojans. The case of 2001 QR322 A&A 41, 725 734 (23) DOI: 1.151/4-6361:231275 c ESO 23 Astronomy & Astrophysics The MATROS project: Stability of Uranus and Neptune Trojans. The case of 21 QR322 F. Marzari 1, P. Tricarico 1, and H. Scholl

More information

The primordial excitation and clearing of the asteroid belt Revisited

The primordial excitation and clearing of the asteroid belt Revisited Icarus 191 (2007) 434 452 www.elsevier.com/locate/icarus The primordial excitation and clearing of the asteroid belt Revisited David P. O Brien a,b,, Alessandro Morbidelli a, William F. Bottke c a Observatoire

More information

The dynamics of Plutinos

The dynamics of Plutinos The dynamics of Plutinos Qingjuan Yu and Scott Tremaine Princeton University Observatory, Peyton Hall, Princeton, NJ 08544-1001, USA ABSTRACT Plutinos are Kuiper-belt objects that share the 3:2 Neptune

More information

arxiv:astro-ph/ v3 2 Mar 2006

arxiv: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 information

Numerous Weak Resonances Drive Asteroids toward Terrestrial Planets Orbits

Numerous Weak Resonances Drive Asteroids toward Terrestrial Planets Orbits Icarus 139, 295 308 (1999) Article ID icar.1999.6097, available online at http://www.idealibrary.com on Numerous Weak Resonances Drive Asteroids toward Terrestrial Planets Orbits A. Morbidelli and D. Nesvorný

More information

The Scattered Disk: Origins, Dynamics and End States

The Scattered Disk: Origins, Dynamics and End States The Scattered Disk: Origins, Dynamics and End States Rodney S. Gomes Observatório Nacional, Rua General José Cristino 77, CEP 2921-, Rio de Janeiro, RJ, Brazil Julio A. Fernández, Tabaré Gallardo Departamento

More information

DIFFUSION OF ASTEROIDS IN MEAN MOTION RESONANCES

DIFFUSION OF ASTEROIDS IN MEAN MOTION RESONANCES DIFFUSION OF ASTEROIDS IN MEAN MOTION RESONANCES KLEOMENIS TSIGANIS and HARRY VARVOGLIS Section of Astrophysics, Astronomy and Mechanics, Department of Physics, University of Thessaloniki, 540 06 Thessaloniki,

More information

Evidence for an extended scattered disk.

Evidence for an extended scattered disk. Evidence for an extended scattered disk. arxiv:astro-ph/0103435v1 27 Mar 2001 B. Gladman 1, M. Holman 2, T. Grav 3, J. Kavelaars 4, P. Nicholson 5, K. Aksnes 3, J-M. Petit 1 1 Observatoire de la Côte d

More information

Frequency map analysis of the 3/1 resonance between planets b and c in the 55 Cancri system ABSTRACT

Frequency map analysis of the 3/1 resonance between planets b and c in the 55 Cancri system ABSTRACT A&A 442, 359 364 (2005) DOI: 151/0004-6361:20053164 c ESO 2005 Astronomy & Astrophysics Frequency map analysis of the 3/1 resonance between planets b and c in the 55 Cancri system F. Marzari 1,H.Scholl

More information

Orbital Structure and Dynamical Evolution of. TNOs. Patryk Sofia Lykawka ( )

Orbital Structure and Dynamical Evolution of. TNOs. Patryk Sofia Lykawka ( ) Orbital Structure and Dynamical Evolution of TNOs Patryk Sofia Lykawka ( ) patryksan@gmail.com Outline I: Introduction and motivation II: III: IV: Dynamical stability and planet migration Stable TNO populations

More information

The Determination of Asteroid Proper Elements

The Determination of Asteroid Proper Elements Knezevic et al.: Determination of Asteroid Proper Elements 63 The Determination of Asteroid Proper Elements Z. Knezevic Astronomical Observatory, Belgrade A. Lemaître Facultes Universitaires Notre-Dame

More information

arxiv:astro-ph/ v1 20 Dec 1999

arxiv:astro-ph/ v1 20 Dec 1999 POPULATION OF THE SCATTERED KUIPER BELT 1 Chadwick A. Trujillo Institute for Astronomy, 2680 Woodlawn Drive, Honolulu, HI 96822 arxiv:astro-ph/9912428v1 20 Dec 1999 chad@ifa.hawaii.edu David C. Jewitt

More information

Inner solar system dynamical analogs of plutinos

Inner solar system dynamical analogs of plutinos Inner solar system dynamical analogs of plutinos Martin Connors a,*, R. Greg Stacey a,b, Ramon Brasser c, and Paul Wiegert d a Centre for Science, Athabasca University, 1 University Drive, Athabasca AB,

More information

arxiv: v2 [astro-ph.ep] 16 Jun 2014

arxiv: v2 [astro-ph.ep] 16 Jun 2014 Mon. Not. R. Astron. Soc. 000, 000 000 (2014) Printed 17 June 2014 (MN LATEX style file v2.2) Extreme trans-neptunian objects and the Kozai mechanism: signalling the presence of trans-plutonian planets

More information

arxiv: v1 [astro-ph.ep] 30 Nov 2011

arxiv: v1 [astro-ph.ep] 30 Nov 2011 Mon. Not. R. Astron. Soc. 000, 000 000 (0000) Printed 1 December 2011 (MN LATEX style file v2.2) An Oort cloud origin for the high-inclination, high-perihelion Centaurs arxiv:1111.7037v1 [astro-ph.ep]

More information

AST111, 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 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 information

Prospects for unseen planets beyond Neptune

Prospects for unseen planets beyond Neptune Prospects for unseen planets beyond Neptune Renu Malhotra The University of Arizona, Tucson, Arizona, USA; renu@lpl.arizona.edu Abstract. Several anomalous features in the orbital distribution of distant

More information

Publ. Astron. Obs. Belgrade No. 90 (2010), DYNAMICAL CHARACTERISTICS OF HUNGARIA ASTEROIDS 1. INTRODUCTION

Publ. Astron. Obs. Belgrade No. 90 (2010), DYNAMICAL CHARACTERISTICS OF HUNGARIA ASTEROIDS 1. INTRODUCTION Publ. Astron. Obs. Belgrade No. 9 (21), 11-18 Invited Lecture DYNAMICAL CHARACTERISTICS OF HUNGARIA ASTEROIDS Z. KNEŽEVIĆ1, B. NOVAKOVIĆ2, A. MILANI 3 1 Astronomical Observatory, Volgina 7, 116 Belgrade,

More information

On the observability of resonant structures in planetesimal disks due to planetary migration. R. Reche, H. Beust, J.-C. Augereau, and O.

On the observability of resonant structures in planetesimal disks due to planetary migration. R. Reche, H. Beust, J.-C. Augereau, and O. A&A 480, 551 561 (2008) DOI: 10.1051/0004-6361:20077934 c ESO 2008 Astronomy & Astrophysics On the observability of resonant structures in planetesimal disks due to planetary migration R. Reche, H. Beust,

More information

THE STRUCTURE OF THE DISTANT KUIPER BELT IN A NICE MODEL SCENARIO

THE STRUCTURE OF THE DISTANT KUIPER BELT IN A NICE MODEL SCENARIO Accepted to the Astronomical Journal: January 21, 2017 Preprint typeset using L A TEX style AASTeX6 v. 1.0 THE STRUCTURE OF THE DISTANT KUIPER BELT IN A NICE MODEL SCENARIO R. E. Pike 1,2, S. Lawler 3,

More information

Survival of Trojan-Type Companions of Neptune During Primordial Planet Migration

Survival of Trojan-Type Companions of Neptune During Primordial Planet Migration Survival of Trojan-Type Companions of Neptune During Primordial Planet Migration Stephen J. Kortenkamp,, Renu Malhotra, and Tatiana Michtchenko ABSTRACT We investigate the survivability of Trojan-type

More information

The Dynamics of Known Centaurs 1

The Dynamics of Known Centaurs 1 The Dynamics of Known Centaurs 1 Matthew S. Tiscareno, Renu Malhotra Lunar and Planetary Laboratory, University of Arizona, Tucson, AZ 85721; matthewt@lpl.arizona.edu, renu@lpl.arizona.edu ABSTRACT We

More information

A Kozai-resonating Earth quasi-satellite

A 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 information

Kozai-Lidov oscillations

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 information

A survey of near-mean-motion resonances between Venus and Earth

A survey of near-mean-motion resonances between Venus and Earth Celest Mech Dyn Astr (2) 7:63 76 DOI 7/s569--9266-6 ORIGINAL ARTICLE A survey of near-mean-motion resonances between Venus and Earth Á. Bazsó R. Dvorak E. Pilat-Lohinger V. Eybl Ch. Lhotka Received: 3

More information

Dynamical evolution of asteroid 2009 HE60, a quasi-satellite of Earth.

Dynamical evolution of asteroid 2009 HE60, a quasi-satellite of Earth. Dynamical evolution of asteroid 2009 HE60, a quasi-satellite of Earth. P. Caspari ª ª Swinburne Astronomy Online, Student 6593658 Not yet received Not yet available online Abstract This paper announces

More information

PLANETARY MIGRATION AND PLUTINO ORBITAL INCLINATIONS R. S. GOMES

PLANETARY MIGRATION AND PLUTINO ORBITAL INCLINATIONS R. S. GOMES THE ASTRONOMICAL JOURNAL, 120:2695È2707, 2000 November ( 2000. The American Astronomical Society. All rights reserved. Printed in U.S.A. PLANETARY MIGRATION AND PLUTINO ORBITAL INCLINATIONS R. S. GOMES

More information

arxiv: v1 [astro-ph.ep] 25 Jun 2009

arxiv: v1 [astro-ph.ep] 25 Jun 2009 Two Dynamical Classes of Centaurs Brenae L. Bailey a, Renu Malhotra b arxiv:0906.4795v1 [astro-ph.ep] 25 Jun 2009 a Program in Applied Mathematics, 617 N. Santa Rita, The University of Arizona, Tucson,

More information

ISIMA lectures on celestial mechanics. 3

ISIMA 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 information

CHAOS AND THE EFFECTS OF PLANETARY MIGRATION ON THE ORBIT OF S/2000 S5 KIVIUQ

CHAOS AND THE EFFECTS OF PLANETARY MIGRATION ON THE ORBIT OF S/2000 S5 KIVIUQ The Astronomical Journal, 128:1899 1915, 2004 October # 2004. The American Astronomical Society. All rights reserved. Printed in U.S.A. CHAOS AND THE EFFECTS OF PLANETARY MIGRATION ON THE ORBIT OF S/2000

More information

Paweł Kankiewicz, Ireneusz Włodarczyk

Paweł Kankiewicz, Ireneusz Włodarczyk Meeting on Asteroids and Comets in Europe Vienna, Austria, 12-14.05.2006 The stability of orbits of effective Mars Crossers Paweł Kankiewicz, Ireneusz Włodarczyk Astrophysics Division, Institute of Physics,

More information

Kuiper Belt Dynamics and Interactions

Kuiper Belt Dynamics and Interactions Kuiper Belt Dynamics and Interactions Minor Planet Center Ruth Murray-Clay Harvard-Smithsonian Center for Astrophysics Kuiper belt µm ejected by radiation pressure larger grains migrate in via PR drag

More information

Two Dynamical Classes of Centaurs

Two Dynamical Classes of Centaurs Two Dynamical Classes of Centaurs Brenae L. Bailey a, Renu Malhotra b a Program in Applied Mathematics, 617 N. Santa Rita, The University of Arizona, Tucson, AZ 85721 b Lunar and Planetary Laboratory,

More information

1. INTRODUCTION. THE ASTRONOMICAL JOURNAL, 121:1730È1735, 2001 March ( The American Astronomical Society. All rights reserved. Printed in U.S.A.

1. INTRODUCTION. THE ASTRONOMICAL JOURNAL, 121:1730È1735, 2001 March ( The American Astronomical Society. All rights reserved. Printed in U.S.A. THE ASTRONOMICAL JOURNAL, 121:1730È1735, 2001 March ( 2001. The American Astronomical Society. All rights reserved. Printed in U.S.A. ON THE SIZE DEPENDENCE OF THE INCLINATION DISTRIBUTION OF THE MAIN

More information

arxiv: v1 [astro-ph.ep] 22 Feb 2016

arxiv: v1 [astro-ph.ep] 22 Feb 2016 Neptune s Orbital Migration Was Grainy, Not Smooth David Nesvorný 1, David Vokrouhlický 2 arxiv:1602.06988v1 [astro-ph.ep] 22 Feb 2016 (1) Department of Space Studies, Southwest Research Institute, 1050

More information

Towards stability results for planetary problems with more than three bodies

Towards 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 information

Meteorite transport Revisited

Meteorite transport Revisited Meteoritics & Planetary Science 52, Nr 8, 1660 1668 (2017) doi: 10.1111/maps.12876 Meteorite transport Revisited Jack WISDOM * Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA

More information

ORBITAL RESONANCES IN PLANETARY SYSTEMS

ORBITAL 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 information

Dynamical Evolution of Ecliptic Comets

Dynamical Evolution of Ecliptic Comets Duncan et al.: Dynamical Evolution of Ecliptic Comets 193 Dynamical Evolution of Ecliptic Comets Martin Duncan Queen s University Harold Levison and Luke Dones Southwest Research Institute Ecliptic comets

More information

A STUDY OF CLOSE ENCOUNTERS BETWEEN MARS AND ASTEROIDS FROM THE 3:1 RESONANCE. Érica C. Nogueira, Othon C. Winter

A STUDY OF CLOSE ENCOUNTERS BETWEEN MARS AND ASTEROIDS FROM THE 3:1 RESONANCE. Érica C. Nogueira, Othon C. Winter A STUDY OF CLOSE ENCOUNTERS BETWEEN MARS AND ASTEROIDS FROM THE 3: RESONANCE Érica C. Nogueira, Othon C. Winter Grupo de Dinâmica Orbital e Planetologia UNESP -- Guaratinguetá -- Brazil Antonio F.B. de

More information

Linking NEAs to their main-belt source regions

Linking NEAs to their main-belt source regions Near-Earth (NEAs) Department of Astronomy, University of Belgrade Stardust ITN: Opening Training School 21st November 2013, Glasgow, Scotland Near-Earth (NEAs) Table of contents 1 Main phases during the

More information

The asteroids. Example for the usage of the Virtual Observatory

The asteroids. Example for the usage of the Virtual Observatory Example for the usage of the Virtual Observatory The asteroids Florian Freistetter, ZAH, Heidelberg florian@ari.uni-heidelberg.de Asteroids in the solar system There are not only planets in our solar system.

More information

Pluto Data: Numbers. 14b. Pluto, Kuiper Belt & Oort Cloud. Pluto Data (Table 14-5)

Pluto Data: Numbers. 14b. Pluto, Kuiper Belt & Oort Cloud. Pluto Data (Table 14-5) 14b. Pluto, Kuiper Belt & Oort Cloud Pluto Pluto s moons The Kuiper Belt Resonant Kuiper Belt objects Classical Kuiper Belt objects Pluto Data: Numbers Diameter: 2,290.km 0.18. Earth Mass: 1.0. 10 22 kg

More information

Secular resonance sweeping of the main asteroid belt during planet migration

Secular resonance sweeping of the main asteroid belt during planet migration 27 pages, 7 figures; submitted to ApJ on May 27, 2010 Secular resonance sweeping of the main asteroid belt during planet migration David A. Minton Southwest Research Institute and NASA Lunar Science Institute

More information

Asteroids in the 2 : 1 resonance with Jupiter: dynamics and size distribution

Asteroids in the 2 : 1 resonance with Jupiter: dynamics and size distribution Mon. Not. R. Astron. Soc. 335, 417 431 (2002) Asteroids in the 2 : 1 resonance with Jupiter: dynamics and size distribution F. Roig, 1 D. Nesvorný 2 and S. Ferraz-Mello 1 1 Instituto Astronômico e Geofísico,

More information

Planetary migration and the Late Heavy Bombardment (better late than never)

Planetary migration and the Late Heavy Bombardment (better late than never) Planetary migration and the Late Heavy Bombardment (better late than never) S. Goldman April 2007 The early history of the Solar System, though still shrouded in mystery, is slowly being pieced together

More information

The behavior of the AZ and QZ solutions in Figure 3 shows only dierences of

The behavior of the AZ and QZ solutions in Figure 3 shows only dierences of years; the latter is magnied by about a factor of 50 for better visualization. The reduced solution has a high and a low eccentricity state with regular and fast transitions between them. Its phase is

More information

Dynamics of two planets in the 3/2 mean-motion resonance: application to the planetary system of the pulsar PSR B

Dynamics of two planets in the 3/2 mean-motion resonance: application to the planetary system of the pulsar PSR B Celestial Mech Dyn Astr (2006) 94:381 397 DOI 10.1007/s10569-006-9002-4 ORIGINAL ARTICLE Dynamics of two planets in the 3/2 mean-motion resonance: application to the planetary system of the pulsar PSR

More information

SOLAR SYSTEM 2019 SAMPLE EXAM

SOLAR SYSTEM 2019 SAMPLE EXAM SOLAR SYSTEM 2019 SAMPLE EXAM Team Name: Team #: No calculators are allowed. All questions are of equal weight unless otherwise noted. Turn in all materials when you have completed the test! Make sure

More information

Origin of the Structure of the Kuiper Belt during a Dynamical Instability in the Orbits of Uranus and Neptune

Origin of the Structure of the Kuiper Belt during a Dynamical Instability in the Orbits of Uranus and Neptune Origin of the Structure of the Kuiper Belt during a Dynamical Instability in the Orbits of Uranus and Neptune Harold F. Levison Southwest Research Institute, USA Alessandro Morbidelli Observatoire de la

More information

arxiv: v1 [astro-ph] 24 Dec 2008

arxiv: 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 information

Yes, inner planets tend to be and outer planets tend to be.

Yes, inner planets tend to be and outer planets tend to be. 1. Planet Density Make some general comments about inner and outer planets density Inner Planets Density Outer Planets Density Is there a pattern or a trend in planet density? Yes, inner planets tend to

More information

arxiv: v1 [astro-ph.ep] 3 Apr 2018

arxiv: 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 information

Habitability in the Upsilon Andromedae System

Habitability in the Upsilon Andromedae System Habitability in the Upsilon Andromedae System Adrienne Dove University of Missouri Columbia Institute for Astronomy University of Hawaii Mentor: Nader Haghighipour ABSTRACT We investigate the habitability

More information

arxiv:astro-ph/ v1 16 Sep 1996

arxiv: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 information

Discovery of a candidate inner Oort cloud planetoid

Discovery of a candidate inner Oort cloud planetoid submitted to ApJ Letters, March 16, 2004 Discovery of a candidate inner Oort cloud planetoid Michael E. Brown Division of Geological and Planetary Sciences, California Institute of Technology, Pasadena,

More information

Missing: a source of short-period comets

Missing: a source of short-period comets Missing: a source of short-period comets Bill Worraker Astronomers have now discovered more than 700 objects accompanying Pluto in the region of space known as the Kuiper Belt, which lies beyond the orbit

More information

Exercise 4.0 PLANETARY ORBITS AND CONFIGURATIONS

Exercise 4.0 PLANETARY ORBITS AND CONFIGURATIONS Exercise 4.0 PLANETARY ORBITS AND CONFIGURATIONS I. Introduction The planets revolve around the Sun in orbits that lie nearly in the same plane. Therefore, the planets, with the exception of Pluto, are

More information

THE CURIOUSLY WARPED MEAN PLANE OF THE KUIPER BELT

THE CURIOUSLY WARPED MEAN PLANE OF THE KUIPER BELT Draft version September 17, 217 Typeset using LATEX default style in AASTeX61 THE CURIOUSLY WARPED MEAN PLANE OF THE KUIPER BELT Kathryn Volk 1 and Renu Malhotra 1 1 Lunar and Planetary Laboratory, The

More information

arxiv: v2 [astro-ph.ep] 4 Mar 2011

arxiv: v2 [astro-ph.ep] 4 Mar 2011 32 pages, 7 figures. Accepted for publication in ApJ on March 1, 2011 arxiv:1102.3131v2 [astro-ph.ep] 4 Mar 2011 Secular resonance sweeping of the main asteroid belt during planet migration David A. Minton

More information

Exploring the orbital evolution of planetary systems

Exploring the orbital evolution of planetary systems Exploring the orbital evolution of planetary systems Tabaré Gallardo Instituto de Física, Facultad de Ciencias, Udelar. (Dated: January 22, 217) The aim of this paper is to encourage the use of orbital

More information

The Kuiper Belt, filled with icy worlds

The Kuiper Belt, filled with icy worlds The Kuiper Belt, filled with icy worlds By NASA.gov, adapted by Newsela staff on 11.16.16 Word Count 956 Level 1140L TOP: Artist's concept showing the exploration of the Kuiper Belt so far. New Horizons

More information

A Survey of the Trans-Neptunian Region

A Survey of the Trans-Neptunian Region Next Generation Space Telescope Ad-Hoc Science Working Group Design Reference Mission Proposal Program contacts: Philip D. Nicholson, Brett J. Gladman Scientific category: SOLAR SYSTEM Instruments: OPT/CAM,

More information

ORBITAL RESONANCES IN PLANETARY SYSTEMS 1

ORBITAL 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 information

The Planet Pluto. & Kuiper Belt. The Search for PLANET X Pluto Discovered. Note how Pluto Moved in 6 days. Pluto (Hades): King of the Underworld

The Planet Pluto. & Kuiper Belt. The Search for PLANET X Pluto Discovered. Note how Pluto Moved in 6 days. Pluto (Hades): King of the Underworld X The Planet Pluto & Kuiper Belt Updated May 9, 2016 The Search for PLANET X Recall Neptune was predicted from observed changes in orbit of Uranus Lowell & Pickering suggest small changes in Neptune s

More information