The jumping-Jupiter scenario specifies an evolution of giant-planet migration described by the Nice model, in which an ice giant (Uranus, Neptune, or an additional Neptune-mass planet) encounters first Saturn and then Jupiter, causing the step-wise separation of their orbits. The jumping-Jupiter scenario was proposed by Ramon Brasser, Alessandro Morbidelli, Rodney Gomes, Kleomenis Tsiganis, and Harold Levison after their studies revealed that the smooth divergent migration of Jupiter and Saturn resulted in an inner Solar System significantly different from the current Solar System. The sweeping of secular resonances through the inner Solar System during the migration excited the eccentricities of the terrestrial planets beyond current values and left an asteroid belt with an excessive ratio of high- to low-inclination objects. The step-wise separation of Jupiter and Saturn described in the jumping-Jupiter scenario allows these resonances to quickly cross the inner Solar System without altering orbits excessively. The jumping-Jupiter scenario also results in a number of other differences with the original Nice model. The fraction of lunar impactors from asteroid belt during the Late Heavy Bombardment is significantly reduced, most of the Jupiter trojans are captured via an alternative mechanism, and Jupiter acquires its population of irregular satellites via the same process as the other planets. The frequent ejection of an ice giant during simulations of the jumping-Jupiter scenario has led some to propose an additional giant planet in the early Solar System.
- 1 Background
- 2 Description
- 3 Implications for the early Solar System
- 4 See also
- 5 References
Original Nice model
The original Nice model begins with the giant planets in a compact configuration with nearly circular orbits. Initially, interactions with planetesimals originating in an outer disk drive the slow divergent migration of the giant planets. This planetesimal-driven migration of the giant planets continues until Jupiter and Saturn cross their mutual 2:1 resonance. The resonance crossing excites the eccentricity of Jupiter and Saturn. The increased eccentricity creates perturbations on Uranus and Neptune, increasing their eccentricities until the system becomes chaotic and orbits begin to intersect. Gravitational encounters between the planets scatter Uranus and Neptune outward into the planetesimal disk. The disk is disrupted, scattering many of the planetesimals onto planet-crossing orbits. A rapid phase of divergent migration of the giant planets is initiated and continues until the disk is depleted. Dynamic friction during this phase dampens the eccentricities of Uranus and Neptune stabilizing the system. In numerical simulations of the original Nice model the final orbits of the giant planets are similar to the current Solar System.
Updated initial conditions
Later versions of the Nice model begin with the giant planets in a series of resonances. This change reflects hydrodynamic models of the early Solar System. In these models, interactions between the giant planets and the gas disk results in the giant planets migrating toward the central star, in some cases becoming hot Jupiters. However, in a multiple-planet system, this inward migration may be halted or reversed if a more rapidly migrating smaller planet is captured in an outer orbital resonance. The Grand Tack hypothesis, with Jupiter's migration being reversed at 1.5 AU following the capture of Saturn in a 3:2 resonance, is an example of this type of orbital evolution. The resonance in which Saturn is captured, a 3:2 or a 2:1 resonance, and the extent of the outward migration (if any) are dependent on the physical properties of the gas disk. The capture of Uranus and Neptune into further resonances during this outward migration results in a quadruply resonant system, with several stable combinations having been identified. Following the dissipation of the gas disk, the quadruple resonance is eventually broken due to interactions with planetesimals from the outer disk. Numerical simulations from this point resembling the original Nice model with an instability beginning shortly after the quadruple resonance is broken or following a period of planetesimal-driven migration when planets cross a different resonance. However, they differ in that there is no slow approach to the 2:1 resonance as Jupiter and Saturn begin in this resonance or cross it rapidly during the instability. There are still some open questions about the damping of the orbital eccentricity (via interactions with planetesimals) of Jupiter and Saturn after the dispersal of disk's gas. In fact, capture of these two planets into the 2:1 mean-motion resonance is far more likely than capture into the 3:2 resonance and evolution in the 2:1 resonance can lead to substantial eccentricity growth during the disk-gas phase.
Alternative instability trigger
The addition of the gravitational interactions between planetesimals to the models revealed another mechanism for triggering a late instability in the outer Solar System. During numerical simulations which included these interactions a transfer of energy between the planetesimals and the giant planets was observed. This energy transfer led to the migration of the planets toward the Sun and occurred even when there were no encounters between planetesimals and the planets. Closer investigation indicated that the energy transfer was due to a coupling between the average eccentricity of the planetesimal disk and the semi-major axes of the outer planets. Because the planets are locked in resonance the inward migration also resulted in an increase in the eccentricity of the inner ice giant. The increased eccentricity changes the precession frequency of the inner ice giant, leading to the crossing of secular resonances during this migration. Eventually, the breaking of the quadruple resonance of the outer planets occurs during one of these secular-resonance crossings. Gravitational encounters begin shortly afterward due to the close proximity of the planets in the previously resonant configuration. Numerical simulations indicate that the timing of the instability caused by this mechanism is fairly independent of the distance between the outer planet and the planetesimal disk and typically occurs after several hundred million years. In combination with the updated initial conditions, this alternative mechanism for triggering a late instability has been called the Nice 2 model.
Solar System constraints
Ramon Brasser, Alessandro Morbidelli, Rodney Gomes, Kleomenis Tsiganis, and Harold Levison published a series of three papers analyzing the orbital evolution of the Solar System during giant planet migration. The impact this migration had on the eccentricities of Jupiter and Saturn, the orbits of the terrestrial planets, and the orbital distribution of the asteroid belt allowed them to identify several constraints on the evolution of the outer Solar System. A number of these were found to be incompatible with smooth planetesimal-driven migration of Jupiter after the 2:1 resonance crossing.
Jupiter and Saturn have modest eccentricities that oscillate out of phase, with Jupiter reaching maximum eccentricity when Saturn reaches its minimum and vice versa. A smooth migration of the planets without resonance crossing results in very small eccentricities. Resonance crossings excites their mean eccentricities, with the 2:1 resonance crossing reproducing Jupiter's current eccentricity, but these do not generate the oscillations in their eccentricities. Recreating both requires either a combination of resonance crossings and an encounter between Saturn and an ice giant, or the encounters of an ice giant with both gas giants.
During the smooth migration of the giant planets the ν5 secular resonance sweeps through the inner Solar System, exciting the eccentricities of the terrestrial planets. For the original Nice model, these eccentricities can reach levels that destabilize the inner Solar System, leading to collisions between planets or the ejection of Mars, during the slow approach of Jupiter and Saturn to their 2:1 resonance. In later versions of the Nice model, Jupiter's and Saturn's divergent migration across (or from) the 2:1 resonance is more rapid and the nearby ν5 resonance crossings are brief, hence not resulting in the excessive excitation of the orbits of Earth and Mars. The authors proposed that the last two resonance crossings were avoided, also preventing the excessive excitation of the eccentricities of Mercury and Venus. This would occur if gravitational encounters between an ice giant and both gas giants caused the Jupiter–Saturn period ratio to jump from below 2.1 to beyond 2.3. The authors named this alternative evolution the jumping-Jupiter scenario.
A smooth planetesimal-driven migration of the giant planets does not result in an orbital distribution that resembles that of the current asteroid belt. The ν6 secular resonance removes low-inclination asteroids and the ν16 secular resonance excites asteroid inclinations, resulting in a ratio of high- to low-inclination asteroids that is too large. The interaction of the ν6 secular resonance with the 3:1 mean-motion resonance also leaves a prominent clump in the semi-major-axis distribution. A giant-planet migration that includes a jumping-Jupiter scenario, in contrast, does not significantly alter the inclination distribution, yielding an asteroid belt with a final orbital distribution that is similar to its initial distribution.
The jumping-Jupiter scenario replaces the smooth separation of Jupiter and Saturn with a series of jumps, thereby avoiding the sweeping of secular resonances through the inner Solar System as their period ratio crosses from 2.1-2.3. In the jumping-Jupiter scenario an ice giant is scattered inward by Saturn onto a Jupiter-crossing orbit and then scattered outward by Jupiter. Saturn's semi-major axis is increased in the first gravitational encounter and Jupiter's reduced by the second with the net result being an increase in their period ratio. In numerical simulations the process can be much more complex: while the trend is for Jupiter's and Saturn's orbits to separate, depending on the geometry of the encounters, individual jumps of Jupiter's and Saturn's semi-major axes can be either up and down. In addition to numerous encounters with Jupiter and Saturn, the ice giant can encounter other ice giant(s) and in some cases cross significant parts of the asteroid belt. The gravitational encounters occur over a period of 10,000–100,000 years, and end when dynamical friction with the planetesimal disk dampens the ice giant's eccentricity, raising its perihelion beyond Saturn's orbit; or when the ice giant is ejected from the Solar System. A jumping-Jupiter scenario occurs in a subset of numerical simulations of the Nice model, including some done for the original Nice model paper. The chances of Saturn scattering an ice giant onto a Jupiter-crossing orbit increases when the initial Saturn–ice giant distance is less than 3 AU, and with the 35-Earth-mass planetesimal belt used in the original Nice model, typically results in the ejection of the ice giant.
Implications for the early Solar System
In addition to preserving low eccentricities of the terrestrial planets and maintaining the pre-migration orbital distribution of the asteroid belt the jumping-Jupiter scenario results in a number of other differences with the original Nice model. These include the source regions for lunar impactors during the Late Heavy Bombardment, constraints on the formation of the asteroid belt, the capture mechanisms for Jupiter's irregular satellites and trojans, and the possibility of additional giant planets in the early Solar System.
Late Heavy Bombardment
The jumping-Jupiter scenario yields a much smaller flux of impactors from the asteroid belt during the Late Heavy Bombardment. Numerical simulations of the asteroid belt during a jumping-Jupiter scenario revealed that roughly 50% of the asteroids were removed. For comparison, 90% of the asteroids were removed during planetary migration in the original Nice model. The mass of asteroids impacting the inner planets is reduced by roughly an order of magnitude in the jumping-Jupiter scenario, potentially leaving comets as the dominant source of the impactors. This conclusion, however, conflicts with evidence indicating that the impactors that created the lunar basins were not comets. A recent paper that offers asteroids from a primordial extension of the asteroid belt to Mars-crossing orbits as the primary source of the LHB of the Moon may resolve this issue. These E-belt asteroids were in stable orbits until the outer planets reached their current configuration at the end of the giant planet migration. Their orbits were then destabilized by the ν6 secular resonance driving them onto Earth-crossing orbits as part of the Late Heavy Bombardment.
The distribution of the orbital elements and the population of the asteroid belt is largely preserved during the planetary migration in the jumping-Jupiter scenario. This places a number of constraints on models of the early Solar System. The excitation and depletion of the asteroid belt and the partial mixing of its taxonomical classes must have occurred during the planetary formation era. If this was the result of gravitational encounters with planetary embryos, bodies larger than Ceres, few if any could have remained in the asteroid belt at the end of that era. Furthermore, because fossil Kirkwood gaps formed while Jupiter was on its initial orbit are not observed, Jupiter must have ended that era on a low-eccentricity orbit. Simulations conducted by Fernando Roig and David Nesvorny indicated that the inclination distribution at the end of the planetary formation era resembled that of the Grand Tack model, but that the eccentricity distribution was less excited. The difference may indicate that higher-eccentricity asteroids left by the Grand Tack were depleted over time by the terrestrial planets, which were not included in their simulation. The rapid migration in the jumping-Jupiter scenario also favors the survival of asteroid collisional families formed during the Late Heavy Bombardment. One example of these is the Hilda collisional family, a subset of the Hilda group, which some researches have suggested was formed at that time because of the current low collision rate. Asteroids in mean-motion resonance with Jupiter, such as those in the 2:1 resonance, can be jump-captured as the semi-major axis of the resonance jumps during Jupiter's encounters. A few asteroids from the asteroid belt may be captured as Hildas, but most originate from the outer Solar System.
In the jumping-Jupiter scenario, most of the Jupiter trojans are captured shortly after a gravitational encounter between Jupiter and an ice giant via a process called jump-capture. During these encounters Jupiter's semi-major axis can jump by as much as 0.2 AU, displacing the L4 and L5 points radially, and releasing any existing Jupiter trojans. New Jupiter trojans are captured from the population of planetesimals with semi-major axes matching Jupiter's new semi-major axis. Following its final encounter with Jupiter the ice giant may pass through one of Jupiter's trojan swarms, scattering the captured bodies and reducing its population, potentially explaining the asymmetry of the L4 and L5 groups. Later, some of the trojans are lost and others captured during weak-resonance crossings as the co-orbital region becomes temporarily chaotic. In numerical simulations the orbital distribution of Jupiter trojans captured was similar to that of the current Solar System and independent of the particular history of the jumping-Jupiter scenario. The capture efficiency was such that the current population could be reproduced starting with a planetesimal belt consistent with models reproducing other aspects of the outer Solar System.
The encounters between the ice giant and Jupiter dynamically perturb the orbits of the Galilean satellites. The perturbations excite the eccentricities and inclinations of the orbits and alter the semi-major axes, potentially breaking the Laplace resonance of Io, Europa and Ganymede. These possibilities were investigated Rogerio Deienno, David Nesvorny, David Vokrouhlicky, and Tadashi Yokoyama using numerical simulations. Although encounters beyond 0.05 AU had little effect, those inside 0.02 AU were found to disrupt the orbits of the satellites, potentially leading to collisions between or the ejections of the satellites. The Laplace resonance of the inner satellites was disrupted in some simulations but was often restored by tidal interactions. The orbit of Callisto was the most affected by the encounters. Its inclination, which is not damped by tidal interactions, places the strongest constraints on the encounters, limiting those between 0.02 AU and 0.03 AU to a few.
Jupiter and Saturn undergo many gravitational encounters with an ice giant in the jumping-Jupiter scenario, allowing Jupiter to capture a population of irregular satellites and increasing the relative size of Saturn's population. During these encounters, the hyperbolic orbits of planetesimals around a gas giant are perturbed by the presence of the ice giant. If the geometry and velocities are right, these three-body interactions leave the planetesimal in a bound orbit when the ice giant recedes. Although loosely bound satellites can also be ejected, tightly bound satellites remain, and over a series of encounters the number of irregular satellites tends to increase. The population and orbits of the irregular satellites captured by Jupiter in the simulations are consistent with its current population. The more-frequent encounters between Saturn and the ice giant increase the size of its population relative to Uranus and Neptune when compared to the original Nice model, producing a closer match with observations.
Fifth giant planet
The early Solar System may have begun with five giant planets. In numerical simulations of the jumping-Jupiter scenario the ice giant is often ejected following its gravitational encounters with Jupiter and Saturn, leaving planetary systems that begin with four giant planets with only three. Although beginning with a higher-mass planetesimal disk was found to stabilize four-planet systems, these simulations end with Jupiter and Saturn too far apart. This problem led David Nesvorný to investigate planetary systems beginning with five giant planets. After conducting thousands of simulations he reported that simulations beginning with five giant planets were 10 times as likely to reproduce the current Solar System.
The most difficult aspect of the current Solar System to reproduce in simulations has been Jupiter's eccentricity. A follow-up study by David Nesvorný and Alessandro Morbidelli reported that even for the best combination of initial conditions this constraint was met in only 7% of simulations. The simulations with the best results began with a significant migration of Neptune into the planetesimal disk. This disrupted the planetesimal disk and drove the divergent migration of Jupiter and Saturn until an instability was triggered. The inner ice giant then began its encounters with Jupiter and Saturn. With a smaller mass of planetesimals remaining less dampening of Jupiter's eccentricity and post-encounter migration of Jupiter and Saturn occurred. Although this evolution yields a good match with the current Solar System the authors noted that a wide variety of outcomes were produced by the jumping-Jupiter scenario and that this case should be considered neither the typical nor the expected result.
A separate study by Konstantin Batygin and Michael Brown also found a low probability of reproducing the current Solar System. However, their study yielded similar probabilities for planetary systems beginning with four and five giant planets. This is in part due to using different criteria to judge success, such as retaining a primordial cold classical Kuiper belt. Their results suggest that preserving a cold classical belt would require the ice giant to be ejected in 10,000 years.
The slow outward migration of Neptune into the planetesimal disk seen in the simulations that best reproduce Jupiter's eccentricity also reproduces the large range of inclinations of objects in the Kuiper belt. Numerical simulations conducted by David Nesvorný using a variety of starting points and timescales revealed that the inclination distribution of the plutinos and of the hot classical Kuiper belt objects were best produced when Neptune migrated smoothly from 24 AU to 30 AU over a timescale of 30 million years. This slow migration provides time for encounters with Neptune to excite the eccentricities and inclinations of the planetesimals. It is also needed for a fraction of these objects to be captured onto stable orbits with sizable inclinations in a three-step process first described by Rodney Gomes. The objects are first scattered from the planetesimal disk onto orbits with larger semi-major axes. After some of these objects are captured in mean-motion resonances with Neptune, their inclinations and eccentricities evolve through processes such as the Kozai mechanism, reducing their eccentricities and increasing their eccentricities. At those lower eccentricities some objects then escape from the mean-motion resonance onto stable orbits while Neptune is migrating. The simulations conducted by Nesvorný ended with an excess of objects in mean-motion resonances when compared to observation. Nesvorný speculated that including the encounters between Neptune and a fifth giant planet, which would cause objects previously captured into the resonances to be released, would resolve this issue.
An outward jump in the semi-major axis of Neptune caused by a gravitational encounter with another planet after Neptune migrated several AU may be responsible for the kernel of the cold classical Kuiper belt objects. The kernel is a concentration of Kuiper belt objects with small eccentricities and inclinations, and semi-major axes of 44–44.5 AU identified by the Canada-France Ecliptic Plane Survey. In numerical simulations conducted by David Nesvorny planetesimals were captured by Neptune's 2:1 mean-motion resonance as Neptune migration out to 28 AU. These object then escaped from this resonance when Neptune's semi-major axis was made to jump, as it would if Neptune had a gravitational encounter with another ice giant, leaving a clump of objects with a semi-major axes near 44 AU. If the jump in Neptune's semi-major axis was outward these object remained in these orbits as the resonance moved away.
- Brasser, R.; Morbidelli, A.; Gomes, R.; Tsiganis, K.; Levison, H.F. (2009). "Constructing the secular architecture of the Solar System II: The terrestrial planets". Astronomy and Astrophysics 507 (2): 1053–1065. arXiv:0909.1891. Bibcode:2009A&A...507.1053B. doi:10.1051/0004-6361/200912878.
- Morbidelli, Alessandro; Brasser, Ramon; Gomes, Rodney; Levison, Harold F.; Tsiganis, Kleomenis (2010). "Evidence from the asteroid belt for a violent past evolution of Jupiter’s orbit". The Astronomical Journal 140 (5): 1391–1401. arXiv:1009.1521. Bibcode:2010AJ....140.1391M. doi:10.1088/0004-6256/140/5/1391.
- Nesvorný, David; Vokrouhlický, David; Morbidelli, Alessandro (2013). "Capture of Trojans by Jumping Jupiter". The Astrophysical Journal 768 (1): 45. arXiv:1303.2900. Bibcode:2013ApJ...768...45N. doi:10.1088/0004-637X/768/1/45.
- Nesvorný, David (2011). "Young Solar System's Fifth Giant Planet?". The Astrophysical Journal Letters 742 (2): L22. arXiv:1109.2949. Bibcode:2011ApJ...742L..22N. doi:10.1088/2041-8205/742/2/L22.
- Batygin, Konstantin; Brown, Michael E.; Betts, Hayden (2012). "Instability-driven Dynamical Evolution Model of a Primordially Five-planet Outer Solar System". The Astrophysical Journal Letters 744 (1): L3. arXiv:1111.3682. Bibcode:2012ApJ...744L...3B. doi:10.1088/2041-8205/744/1/L3.
- Morbidelli, Alesandro (2010). "A coherent and comprehensive model of the evolution of the outer Solar System". Comptes Rendus Physique 11 (9–10): 651–659. arXiv:1010.6221. Bibcode:2010CRPhy..11..651M. doi:10.1016/j.crhy.2010.11.001.
- Lin, D. N. C.; Bodenheimer, P.; Richardson, D. C. (1996). "Orbital migration of the planetary companion of 51 Pegasi to its present location". Nature 380 (6575): 606–607. doi:10.1038/380606a0.
- Masset, F.; Snellgrove, M. (2001). "Reversing type II migration: resonance trapping of a lighter giant protoplanet". Monthly Notices of the Royal Astronomical Society 320 (4): L55–L59. arXiv:astro-ph/0003421. Bibcode:2001MNRAS.320L..55M. doi:10.1046/j.1365-8711.2001.04159.x.
- Walsh, Kevin J.; Morbidelli, Alessandro; Raymond, Sean N.; O'Brien, David P.; Mandell, Avi M. "A low mass for Mars from Jupiter's early gas-driven migration". Nature 475 (7335): 206–209. arXiv:1201.5177. doi:10.1038/nature10201.
- Pierens, A.; Nelson, R. P (2008). "Constraints on resonant–trapping for two planets embedded in a protoplanetary disc". Astronomy and Astrophysics 482 (1): 333–340. arXiv:0802.2033. Bibcode:2008A&A...482..333P. doi:10.1051/0004-6361:20079062.
- D'Angelo, G.; Marzari, F. (2012). "Outward Migration of Jupiter and Saturn in Evolved Gaseous Disks". The Astrophysical Journal 757 (1): 50. arXiv:1207.2737. Bibcode:2012ApJ...757...50D. doi:10.1088/0004-637X/757/1/50.
- Pierens, Arnaud; Raymond, Sean N; Nesvorny, David; Morbidelli, Alessandro. "Outward Migration of Jupiter and Saturn in 3:2 or 2:1 Resonance in Radiative Disks: Implications for the Grand Tack and Nice models". The Astrophysical Journal Letters 795 (1): L11. arXiv:1410.0543. doi:10.1088/2041-8205/795/1/L11.
- Morbidelli, Alessandro; Tsiganis, Kleomenis; Crida, Aurélien; Levison, Harold F.; Gomes, Rodney (2007). "Dynamics of the Giant Planets of the Solar System in the Gaseous Protoplanetary Disk and Their Relationship to the Current Orbital Architecture". The Astronomical Journal 134 (5): 1790–1798. arXiv:0706.1713. Bibcode:2007AJ....134.1790M. doi:10.1086/521705.
- Batygin, Konstantin; Brown, Michael E. (2010). "Early Dynamical Evolution of the Solar System: Pinning Down the Initial Conditions of the Nice Model". The Astrophysical Journal 716 (2): 1323–1331. arXiv:1004.5414. Bibcode:2010ApJ...716.1323B. doi:10.1088/0004-637X/716/2/1323.
- Levison, Harold F.; Morbidelli, Alessandro; Tsiganis, Kleomenis; Nesvorný, David; Gomes, Rodney (2011). "Late Orbital Instabilities in the Outer Planets Induced by Interaction with a Self-gravitating Planetesimal Disk" (PDF). The Astronomical Journal 142 (5): 152. Bibcode:2011AJ....142..152L. doi:10.1088/0004-6256/142/5/152.
- Morbidelli, A.; Brasser, R.; Tsiganis, K.; Gomes, R.; Levison, H. F (2009). "Constructing the secular architecture of the Solar System I. The giant planets". Astronomy and Astrophysics 507 (2): 1041–1052. Bibcode:2009A&A...507.1041M. doi:10.1051/0004-6361/200912876.
- Agnor, Craig B.; Lin, D. N. C. (2012). "On the Migration of Jupiter and Saturn: Constraints from Linear Models of Secular Resonant Coupling with the Terrestrial Planets". The Astrophysical Journal 745 (2): 143. arXiv:1110.5042. Bibcode:2012ApJ...745..143A. doi:10.1088/0004-637X/745/2/143.
- Roig, Fernando; Nesvorný, David. "The evolution of asteroids in the jumping-Jupiter migration model". eprint arXiv.
- Tsiganis, K.; Gomes, R.; Morbidelli, A.; Levison, H. F. (2005). "Origin of the orbital architecture of the giant planets of the Solar System". Nature 435 (7041): 459–461. Bibcode:2005Natur.435..459T. doi:10.1038/nature03539. PMID 15917800.
- Gomes, R.; Levison, H. F.; Tsiganis, K.; Morbidelli, A. (2005). "Origin of the cataclysmic Late Heavy Bombardment period of the terrestrial planets". Nature 435 (7041): 466–469. Bibcode:2005Natur.435..466G. doi:10.1038/nature03676. PMID 15917802.
- Kring, David A.; Cohen, Barbara A. (2002). "Cataclysmic bombardment throughout the inner Solar System 3.9–4.0 Ga" (PDF). Journal of Geophysical Research (Planets) 107 (E2): 4–1. doi:10.1029/2001JE001529.
- Strom, Robert G.; Malhotra, Renu; Ito, Takashi; Yoshida, Fumi; Kring, David A. (2005). "The Origin of Planetary Impactors in the Inner Solar System" (PDF). Science 309 (5742): 1847–1860. arXiv:astro-ph/0510200. Bibcode:2005Sci...309.1847S. doi:10.1126/science.1113544.
- Bottke, William F.; Vokrouhlický, David; Minton, David; Nesvorný, David; Morbidelli, Alessandro; Brasser, Ramon; Simonson, Bruce; Levison, Harold F. (2012). "An Archaean heavy bombardment from a destabilized extension of the asteroid belt". Nature 485 (7396): 78–81. Bibcode:2012Natur.485...78B. doi:10.1038/nature10967.
- Brož, M.; Morbidelli, A.; Bottke, W. F.; Rozehnal, J.; Vokrouhlický, D.; Nesvorný, D. (2013). "Constraining the cometary flux through the asteroid belt during the late heavy bombardment". Astronomy & Astrophysics 551: 16. arXiv:1301.6221. Bibcode:2013A&A...551A.117B. doi:10.1051/0004-6361/201219296. A177.
- Brož, M.; Vokrouhlický, D.; Morbidelli, A.; Nesvorný, D.; Bottke, W. F. (2011). "Did the Hilda collisional family form during the late heavy bombardment?". Monthly Notices of the Royal Astronomical Society 414 (3): 2716–2727. arXiv:1109.1114. Bibcode:2011MNRAS.414.2716B. doi:10.1111/j.1365-2966.2011.18587.x.
- Chrenko, O.; Brož, M.; Nesvorný, D.; Tsiganis, K.; Skoulidou, D. K. "The origin of long-lived asteroids in the 2:1 mean-motion resonance with Jupiter". Monthly Notices of the Royal Astronomical Society 451 (3): 2399–2416. doi:10.1093/mnras/stv1109.
- Morbidelli, A.; Levison, H. F.; Tsiganis, K.; Gomes, R. (2005). "Chaotic capture of Jupiter's Trojan asteroids in the early Solar System". Nature 435 (7041): 462–465. Bibcode:2005Natur.435..462M. doi:10.1038/nature03540. PMID 15917801.
- Deienno, Rogerio; Nesvorný, David; Vokrouhlický, David; Yokoyama, Tadashi (2014). "Orbital Perturbations of the Galilean Satellites during Planetary Encounters". The Astronomical Journal 148 (2): 25. arXiv:1405.1880. Bibcode:2014AJ....148...25D. doi:10.1088/0004-6256/148/2/25.
- Nesvorný, David; Vokrouhlický, David; Morbidelli, Alessandro (2007). "Capture of Irregular Satellites during Planetary Encounters". The Astronomical Journal 133 (5): 1962–1976. Bibcode:2007AJ....133.1962N. doi:10.1086/512850.
- Nesvorný, David; Vokrouhlický, David; Deienno, Rogerio. "Capture of Irregular Satellites at Jupiter". The Astrophysical Journal 784 (1): 22. arXiv:1401.0253. Bibcode:2014ApJ...784...22N. doi:10.1088/0004-637X/784/1/22.
- Jewitt, David; Haghighipour, Nader. "Irregular Satellites of the Planets: Products of Capture in the Early Solar System". Annual Review of Astronomy & Astrophysics 45 (1): 261–295. arXiv:astro-ph/0703059. doi:10.1146/annurev.astro.44.051905.092459.
- Stuart, Colin. "Was a giant planet ejected from our Solar System?". Physics World. Retrieved 16 January 2014.
- Nesvorný, David; Morbidelli, Alessandro (2012). "Statistical Study of the Early Solar System's Instability with Four, Five, and Six Giant Planets". The Astronomical Journal 144 (4): 117. arXiv:1208.2957. Bibcode:2012AJ....144..117N. doi:10.1088/0004-6256/144/4/117.
- Nesvorný, David. "Evidence for Slow Migration of Neptune from the Inclination Distribution of Kuiper Belt Objects" (PDF). The Astronomical Journal 150 (3): 73. doi:10.1088/0004-6256/150/3/73.
- Gomes, Rodney. "The origin of the Kuiper Belt high-inclination population". Icarus 161 (2): 404–418.
- Petit, J.-M.; Gladman, B.; Kavelaars, J. J.; Jones, R. L.; Parker, J. (2011). "Reality and origin of the Kernel of the classical Kuiper Belt" (PDF). EPSC-DPS Joint Meeting (2–7 October 2011).
- Nesvorný, David. "Jumping Neptune Can Explain the Kuiper Belt Kernel" (PDF). The Astronomical Journal 150 (3): 68. doi:10.1088/0004-6256/150/3/68.