3 Numerical results – I. Glance at the raw data
In this section we briefly review the long-term stability of planetary orbital motion through some snapshots of raw numerical data. The orbital motion of planets indicates long-term stability in all of our numerical integrations: no orbital crossings nor close encounters between any pair of planets took place.
3.1 General description of the stability of planetary orbits
First， we briefly look at the general character of the long-term stability of planetary orbits. Our interest here focuses particularly on the inner four terrestrial planets for which the orbital time-scales are much shorter than those of the outer five planets. As we can see clearly from the planar orbital configurations shown in Figs 2 and 3， orbital positions of the terrestrial planets differ little between the initial and final part of each numerical integration， which spans several Gyr. The solid lines denoting the present orbits of the planets lie almost within the swarm of dots even in the final part of integrations (b) and (d). This indicates that throughout the entire integration period the almost regular variations of planetary orbital motion remain nearly the same as they are at present.
Vertical view of the four inner planetary orbits (from the z -axis direction) at the initial and final parts of the integrationsN±1. The axes units are au. The xy -plane is set to the invariant plane of Solar system total angular momentum.(a) The initial part ofN+1 ( t = 0 to 0.0547 × 10 9 yr).(b) The final part ofN+1 ( t = 4.9339 × 10 8 to 4.9886 × 10 9 yr).(c) The initial part of N?1 (t= 0 to ?0.0547 × 109 yr).(d) The final part ofN?1 ( t =?3.9180 × 10 9 to ?3.9727 × 10 9 yr). In each panel， a total of 23 684 points are plotted with an interval of about 2190 yr over 5.47 × 107 yr . Solid lines in each panel denote the present orbits of the four terrestrial planets (taken from DE245).
这章没有结束，请点击下一页继续阅读！