The Nice (//) model is a scenario for the dynamical evolution of the Solar System. It is named for the location of the Observatoire de la Côte d'Azur, where it was initially developed, in Nice, France. It proposes the migration of the giant planets from an initial compact configuration into their present positions, long after the dissipation of the initial protoplanetary gas disk. In this way, it differs from earlier models of the Solar System's formation. This planetary migration is used in dynamical simulations of the Solar System to explain historical events including the Late Heavy Bombardment of the inner Solar System, the formation of the Oort cloud, and the existence of populations of small Solar System bodies including the Kuiper belt, the Neptune and Jupiter trojans, and the numerous resonant trans-Neptunian objects dominated by Neptune. Its success at reproducing many of the observed features of the Solar System means that it is widely accepted as the current most realistic model of the Solar System's early evolution, although it is not universally favoured among planetary scientists. One of its limitations that it does not reproduce the outer-system satellites and the Kuiper belt (see below).
The original core of the Nice model is a triplet of papers published in the general science journal Nature in 2005 by an international collaboration of scientists: Rodney Gomes, Hal Levison, Alessandro Morbidelli and Kleomenis Tsiganis. In these publications, the four authors proposed that after the dissipation of the gas and dust of the primordial Solar System disk, the four giant planets (Jupiter, Saturn, Uranus and Neptune) were originally found on near-circular orbits between ~5.5 and ~17 astronomical units (AU), much more closely spaced and compact than in the present. A large, dense disk of small, rock and ice planetesimals, their total about 35 Earth masses, extended from the orbit of the outermost giant planet to some 35 AU.
Scientists understand so little about the formation of Uranus and Neptune that Levison states, "...the possibilities concerning the formation of Uranus and Neptune are almost endless." However, it is suggested that this planetary system evolved in the following manner. Planetesimals at the disk's inner edge occasionally pass through gravitational encounters with the outermost giant planet, which change the planetesimals' orbits. The planets scatter the majority of the small icy bodies that they encounter inward, exchanging angular momentum with the scattered objects so that the planets move outwards in response, preserving the angular momentum of the system. These planetesimals then similarly scatter off the next planet they encounter, successively moving the orbits of Uranus, Neptune, and Saturn outwards. Despite the minute movement each exchange of momentum can produce, cumulatively these planetesimal encounters shift (migrate) the orbits of the planets by significant amounts. This process continues until the planetesimals interact with the innermost and most massive giant planet, Jupiter, whose immense gravity sends them into highly elliptical orbits or even ejects them outright from the Solar System. This, in contrast, causes Jupiter to move slightly inward.
The low rate of orbital encounters governs the rate at which planetesimals are lost from the disk, and the corresponding rate of migration. After several hundreds of millions of years of slow, gradual migration, Jupiter and Saturn, the two inmost giant planets, cross their mutual 1:2 mean-motion resonance. This resonance increases their orbital eccentricities, destabilizing the entire planetary system. The arrangement of the giant planets alters quickly and dramatically. Jupiter shifts Saturn out towards its present position, and this relocation causes mutual gravitational encounters between Saturn and the two ice giants, which propel Neptune and Uranus onto much more eccentric orbits. These ice giants then plough into the planetesimal disk, scattering tens of thousands of planetesimals from their formerly stable orbits in the outer Solar System. This disruption almost entirely scatters the primordial disk, removing 99% of its mass, a scenario which explains the modern-day absence of a dense trans-Neptunian population. Some of the planetesimals are thrown into the inner Solar System, producing a sudden influx of impacts on the terrestrial planets: the Late Heavy Bombardment.
Eventually, the giant planets reach their current orbital semi-major axes, and dynamical friction with the remaining planetesimal disc damps their eccentricities and makes the orbits of Uranus and Neptune circular again.
In some 50% of the initial models of Tsiganis and colleagues, Neptune and Uranus also exchange places. An exchange of Uranus and Neptune would be consistent with models of their formation in a disk that had a surface density that declined with distance from the Sun, which predicts that the masses of the planets should also decline with distance from the Sun.
Solar System features
Running dynamical models of the Solar System with different initial conditions for the simulated length of the history of the Solar System will produce the various populations of objects within the Solar System. As the initial conditions of the model are allowed to vary, each population will be more or less numerous, and will have particular orbital properties. Proving a model of the evolution of the early Solar System is difficult, since the evolution cannot be directly observed. However, the success of any dynamical model can be judged by comparing the population predictions from the simulations to astronomical observations of these populations. At the present time, computer models of the Solar System that are begun with the initial conditions of the Nice scenario best match many aspects of the observed Solar System.
The Late Heavy Bombardment
The crater record on the Moon and on the terrestrial planets is part of the main evidence for the Late Heavy Bombardment (LHB): an intensification in the number of impactors, at about 600 million years after the Solar System's formation. In the initial Nice model icy planetesimals are scattered onto planet-crossing orbits when the outer disc is disrupted by Uranus and Neptune causing a sharp spike of impacts by icy objects. The migration of outer planets also causes mean-motion and secular resonances to sweep through the inner Solar System. In the asteroid belt these excite the eccentricities of the asteroids driving them onto orbits that intersect those of the terrestrial planets causing a more extended period of impacts by stony objects and removing roughly 90% of its mass.
In this initial model the number of planetesimals that would reach the Moon is consistent with the crater record from the LHB. Impacts onto Jupiter's moons are sufficient to trigger Ganymede's differentiation but not Callisto's. However, significant differences with observations have also been identified. The orbital distribution of the surviving asteroid belt has an excess of high inclination objects due to the secular resonances exciting inclinations and also removing low inclination objects. During Jupiter and Saturn's slow approach to the 2:1 resonance leading up to the LHB Mars's eccentricity can reach a level that results in the destabilization of the inner Solar System. The eccentricities of the other terrestrial planets can also be excited beyond current levels by sweeping secular resonances. In the outer Solar System the impacts of icy planetesimals onto Saturn's inner moons are sufficient to vaporize their ice. These issue provided motivations for modifications to the initial Nice model.
Trojans and the asteroid belt
During the period of orbital disruption following Jupiter and Saturn reaching the 2:1 resonance, the combined gravitational influence of the migrating giant planets would have quickly destabilized any existing Trojan groups in the L4 and L5 Lagrange points of Jupiter and Neptune. During this time, the Trojan co-orbital region is termed "dynamically open". Under the Nice model, the planetesimals leaving the disrupted disk cross this region in large numbers, temporarily inhabiting it. After the period of orbital instability ends, the Trojan region is "dynamically closed", capturing planetesimals present at the time. The present Trojan populations are then these acquired scattered planetesimals of the primordial asteroid belt. This simulated population matches the libration angle, eccentricity and the large inclinations of the orbits of the Jupiter Trojans. Their inclinations had not previously been understood. This mechanism of the Nice model similarly generates the Neptune trojans.
A large number of planetesimals would have also been captured in Jupiter's mean motion resonances as Jupiter migrated inward. Those that remained in a 3:2 resonance with Jupiter form the Hilda family. The eccentricity of other objects declined while they were in a resonance and escaped onto stable orbits in the outer asteroid belt, at distances greater than 2.6 AU as the resonances moved inward. These captured objects would then have undergone collisional erosion, grinding the population away into smaller fragments that can then be acted on by the Yarkovsky effect, causing small objects to drift into unstable resonances, and Poynting–Robertson drag causing smaller grains to drift toward the sun. These processes remove more than 90% of the origin mass implanted into the asteroid belt according to Bottke and colleagues. The size frequency distribution of this simulated population following this erosion are in excellent agreement with observations. This suggests that the Jupiter Trojans, Hildas and some of the outer asteroid belt, all spectral D-type asteroids, are the remnant planetesimals from this capture and erosion process. It has also been suggested that the dwarf planet Ceres was captured via this process.
Any original populations of irregular satellites captured by traditional mechanisms, such as drag or impacts from the accretion disks, would be lost during the interactions of the planets at the time of global system instability. In the Nice model, large numbers of planetesimals interact with the outer planets at this time, and some are captured during three-way interactions with those planets. The probability for any planetesimal to be captured by an ice giant is relatively high, a few 10−7. These new satellites could be captured at almost any angle, so unlike the regular satellites of Saturn, Uranus and Neptune, they do not necessarily orbit in the planets' equatorial planes. Some irregulars may have even been exchanged between planets. The resulting irregular orbits match well with the observed populations' semimajor axes, inclinations and eccentricities. Subsequent collisions between these captured satellites may have created the suspected collisional families seen today. These collisions are also required to erode the population to the present size distribution.
Triton, the largest moon of Neptune, can be explained if it was captured in a three-body interaction involving the disruption of a binary planetoid. Such binary disruption would be more likely if Triton was the smaller member of the binary. However, Triton's capture would be more likely in the early Solar System when the gas disk would damp relative velocities, and binary exchange reactions would not in general have supplied the large number of small irregulars.
There were not have been enough interactions with Jupiter to explain Jupiter's retinue of irregulars in the initial Nice model simulations that reproduced other aspects of the outer Solar System. This suggests either that a second mechanism was at work for that planet, or that the early simulations did not reproduce the evolution of the giant planets orbits.
Formation of the Kuiper belt
The migration of the outer planets is also necessary to account for the existence and properties of the Solar System's outermost regions. Originally, the Kuiper belt was much denser and closer to the Sun, with an outer edge at approximately 30 AU. Its inner edge would have been just beyond the orbits of Uranus and Neptune, which were in turn far closer to the Sun when they formed (most likely in the range of 15–20 AU), and in opposite locations, with Uranus farther from the Sun than Neptune.
Gravitational encounters between the planets scatter Neptune outward into the planetesimal disk with a semi-major axis of ~28 AU and an eccentricity as high as 0.4. Neptune's high eccentricity causes it mean-motion resonances to overlap and orbits in the region between Neptune and its 2:1 mean motion resonances to become chaotic. The orbits of objects between Neptune and the edge of the planetesimal disk at this time can evolve outward onto stable low eccentricity orbits within this region. When Neptune's eccentricity is damped by dynamical friction they become trapped on these orbits. These objects form a dynamically cold belt since their inclinations remain small during the short time they interact with Neptune. Later, as Neptune migrates outward on a low eccentricity orbit, objects that have been scattered outward are captured into its resonances can have their eccentricities decline and their inclinations increase due to the Kozai mechanism, allowing them to escape onto stable higher inclination orbits. Other objects remain captured in resonance, forming the plutinos and other resonant populations. These two population are dynamically hot, with higher inclinations and eccentricities, due to their being scattered outward and the longer period these objects interact with Neptune.
This evolution of Neptune's orbit produces both resonant and non-resonant populations, an outer edge at Neptune's 2:1 resonance, and a small mass relative to the original planetesimal disk. The excess of low inclination plutinos in other models is avoided due to Neptune being scattered outward, leaving its 3:2 resonance beyond the original edge of the planetesimal disk. The differing initial locations, with the cold classical objects originating primarily from the outer disk, and capture processes offer explanations for the bi-modal inclination distribution and its correlation with compositions. However, this evolution of Neptune's orbit fails to account for some of the characteristics of the orbital distribution. It predicts a greater average eccentricity in classical Kuiper belt object orbits than is observed (0.10–0.13 versus 0.07) and it does not produce enough higher inclination objects. It also cannot explain the apparent complete absence of gray objects in the cold population, although it has been suggested that color differences arise in part from surface evolution processes rather than entirely from differences in primordial composition.
The shortage of the lowest eccenticty objects predicted in the Nice model may indicate that the cold population formed at its location. The hot and cold populations not only possess different orbits, but different colors; the cold population is markedly redder than the hot, suggesting it has a different composition and formed in a different region. The cold population also includes a large number of binary objects with loosely bound orbits that would be unlikely to survive close encounter with Neptune. If the cold population formed at it current location preserving it would require that Neptune's eccentricity remained small, or that it perihelion precessed rapidly due to a strong interaction between it and Uranus.
Scattered disc and Oort cloud
Objects scattered outward by Neptune onto orbits with semi-major axis greater than 50 AU can be captured in resonances forming the resonant population of the scattered disc, or if their eccentricities are reduced while in resonance they can escape from the resonance onto stable orbits in the scattered disc while Neptune is migrating. When Neptune's eccentricity is large its aphelion can reach well beyond it current orbit. Objects that attain perihelia close to or larger than Neptune's at this time can become detached from Neptune when its eccentricity is damped reducing its aphelion, leaving them on stable orbits in the scattered disc.
Objects scattered outward by Uranus and Neptune onto larger orbits (roughly 5,000 AU) can have their perihelion raised by the galactic tide detaching them from the influence of the planets forming the inner Oort cloud with moderate inclinations. Others that reach even larger orbits can be perturbed by nearby stars forming the outer Oort cloud with isotropic inclinations. Objects scattered by Jupiter and Saturn are typically ejected from the Solar System Several percent of the initial planetesimal disc can be deposited in these reservoirs.
The Nice model has undergone significant modification since its initial publication. The initial conditions of the model have been changed as a result of investigations of the behavior of planets orbiting in a gas disk to a quadruple resonant configuration with Jupiter and Saturn in their 3:2 resonance. The gravitational stirring of the outer planetesimal disk by Pluto-sized objects has been shown to result in breaking of the quadruple resonance via a mechanism that is not sensitive to the distance between the outer planet and the planetesimal disk. This mechanism for triggering the late instability of resonant planets similar to that in the original Nice model has been referred to as the Nice 2 model.
The smooth divergent migration of Jupiter and Saturn following the instability that resulted in the excitation of the eccentricities of the terrestrial planets beyond their current values and an asteroid belt with an excessive ratio of high- to low-inclination objects has been replaced by a step-wise separation of their orbits called the jumping-Jupiter scenario. This is driven by an ice giant encountering Saturn, causing it orbit to expand, then encountering Jupiter, causing its orbit to shrink. The encounters between ice giant and Jupiter in this model also allow Jupiter to acquire its own irregular satellites.
Five-planet Nice model
The frequent ejection of the ice giant encountering Jupiter has led David Nesvorný to hypothesize an early Solar System with five giant planets, one of which was ejected during the instability. In recent work this version of the Nice model begins with a 20 Earth-mass planetesimals disk orbiting beyond the five giant planets that are in a resonant chain. Following the breaking of the resonant chain Neptune first migrates outward into the planetesimal disk reaching 28 AU before encounters between planets begin. This initial migration reduces the mass of the outer disk enabling Jupiter's eccentricity to be preserved and produces a Kuiper belt with an inclination distribution that matches observations. The lower mass planetesimal belt in combination with the excitation of inclinations and eccentricities by the Pluto-massed objects also significantly reduce the loss of ice by Saturn's inner moons. A study published in 2015 based on a large ensemble of n-body simulations suggests that the ejection of a hypothetical fifth giant planet is statistically unlikely to produce the observed orbits of the terrestrial planets. The study concluded that, if there was a giant-planet instability leading to the ejection of one or more additional ice giants, it had to occur prior to the formation of the terrestrial planets and that such an instability could not be the source of the Late Heavy Bombardment.
- "Solving solar system quandaries is simple: Just flip-flop the position of Uranus and Neptune". Press release. Arizona State University. 11 Dec 2007. Retrieved 2009-03-22.
- Desch, S. (2007). "Mass Distribution and Planet Formation in the Solar Nebula". The Astrophysical Journal. 671 (1): 878–893. Bibcode:2007ApJ...671..878D. doi:10.1086/522825.
- Crida, A. (2009). "Solar System formation". Reviews in Modern Astronomy. 21: 3008. arXiv: . Bibcode:2009RvMA...21..215C. doi:10.1002/9783527629190.ch12.
- R. Gomes; H. F. Levison; K. Tsiganis; A. Morbidelli (2005). "Origin of the cataclysmic Late Heavy Bombardment period of the terrestrial planets". Nature. 435 (7041): 466–9. Bibcode:2005Natur.435..466G. doi:10.1038/nature03676. PMID 15917802.
- Tsiganis, K.; Gomes, R.; Morbidelli, A.; F. Levison, H. (2005). "Origin of the orbital architecture of the giant planets of the Solar System" (PDF). Nature. 435 (7041): 459–461. Bibcode:2005Natur.435..459T. doi:10.1038/nature03539. PMID 15917800.
- Morbidelli, A.; Levison, H.F.; Tsiganis, K.; Gomes, R. (2005). "Chaotic capture of Jupiter's Trojan asteroids in the early Solar System" (PDF). Nature. 435 (7041): 462–465. Bibcode:2005Natur.435..462M. doi:10.1038/nature03540. OCLC 112222497. PMID 15917801.
- G. Jeffrey Taylor (21 August 2001). "Uranus, Neptune, and the Mountains of the Moon". Planetary Science Research Discoveries. Hawaii Institute of Geophysics & Planetology. Retrieved 2008-02-01.
- Hansen, Kathryn (June 7, 2005). "Orbital shuffle for early solar system". Geotimes. Retrieved 2007-08-26.
- Levison HF, Morbidelli A, Van Laerhoven C, Gomes RS, Tsiganis K (2007). "Origin of the Structure of the Kuiper Belt during a Dynamical Instability in the Orbits of Uranus and Neptune". Icarus. 196 (1): 258–273. arXiv: . Bibcode:2008Icar..196..258L. doi:10.1016/j.icarus.2007.11.035.
- T. V. Johnson; J. C. Castillo-Rogez; D. L. Matson; A. Morbidelli; J. I. Lunine. "Constraints on outer Solar System early chronology" (PDF). Early Solar System Impact Bombardment conference (2008). Retrieved 2008-10-18.
- Baldwin, Emily. "Comet impacts explain Ganymede-Callisto dichotomy". Astronomy Now. Retrieved 23 December 2016.
- 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–1501. arXiv: . Bibcode:2010AJ....140.1391M. doi:10.1088/0004-6256/140/5/1391.
- 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: . Bibcode:2009A&A...507.1053B. doi:10.1051/0004-6361/200912878.
- Nimmo, F.; Korycansky, D. G. (2012). "Impact-driven ice loss in outer Solar System satellites: Consequences for the Late Heavy Bombardment". Icarus. 219 (1): 508–510. Bibcode:2012Icar..219..508N. doi:10.1016/j.icarus.2012.01.016.
- Levison, Harold F.; Shoemaker, Eugene M.; Shoemaker, Carolyn S. (1997). "Dynamical evolution of Jupiter's Trojan asteroids". Nature. 385 (6611): 42–44. Bibcode:1997Natur.385...42L. doi:10.1038/385042a0.
- Levison, Harold F.; Bottke, William F.; Gounelle, Matthieu; Morbidelli, Alessandro; Nesvorny, David; Tsiganis, Kleomeis (2009). "Contamination of the asteroid belt by primordial trans-Neptunian objects". Nature. 460 (7253): 364–366. doi:10.1038/nature08094.
- Bottke, W. F.; Levison, H. F.; Morbidelli, A.; Tsiganis, K. (2008). "The Collisional Evolution of Objects Captured in the Outer Asteroid Belt During the Late Heavy Bombardment". 39th Lunar and Planetary Science Conference. 39th Lunar and Planetary Science Conference. 39 (LPI Contribution No. 1391): 1447. Bibcode:2008LPI....39.1447B.
- William B. McKinnon (2008). "On The Possibility Of Large KBOs Being Injected Into The Outer Asteroid Belt". Bulletin of the American Astronomical Society. 40: 464. Bibcode:2008DPS....40.3803M.
- Turrini & Marzari, 2008, Phoebe and Saturn's irregular satellites: implications for the collisional capture scenario
- Nesvorný, D.; Vokrouhlický, D.; Morbidelli, A. (2007). "Capture of Irregular Satellits during Planetary Encounters". The Astronomical Journal. 133 (5): 1962–1976. Bibcode:2007AJ....133.1962N. doi:10.1086/512850.
- Nesvorný, David; Beaugé, Cristian; Dones, Luke (2004). "Collisional Origin of Families of Irregular Satellites". The Astronomical Journal. 127 (3): 1768–1783. Bibcode:2004AJ....127.1768N. doi:10.1086/382099.
- Bottke, William F.; Nesvorný, David; Vokrouhlick, David; Morbidelli, Alessandro (2010). "The Irregular Satellites: The Most Collisionally Evolved Populations in the Solar System". The Astronomical Journal. 139 (3): 994–1014. Bibcode:2010AJ....139..994B. doi:10.1088/0004-6256/139/3/994.
- Agnor, Craig B.; Hamilton, Douglas B. (2006). "Neptune's capture of its moon Triton in a binary-planet gravitational encounter". Nature. 441 (7090): 192–194. doi:10.1038/nature04792. PMID 16688170.
- Vokrouhlický, David; Nesvorný, David; Levison, Harold F. (2008). "Irregular Satellite Capture by Exchange Reactions". The Astronomical Journal. 136 (4): 1463–1476. Bibcode:2008AJ....136.1463V. doi:10.1088/0004-6256/136/4/1463.
- Levison, Harold F.; Morbidelli, Alessandro; VanLaerhoven, Christa; Gomes, Rodney S. (2008-04-03). "Origin of the structure of the Kuiper belt during a dynamical instability in the orbits of Uranus and Neptune". Icarus. 196 (1): 258–273. arXiv: . Bibcode:2008Icar..196..258L. doi:10.1016/j.icarus.2007.11.035. Retrieved 2012-05-26.
- Morbidelli, Alessandro (2006). "Origin and dynamical evolution of comets and their reservoirs". arXiv: [astro-ph].
- Lovett, Rick (2010). "Kuiper Belt may be born of collisions". Nature. doi:10.1038/news.2010.522.
- Wolff, Schuyler; Dawson, Rebekah I.; Murray-Clay, Ruth A. (2012). "Neptune on Tiptoes: Dynamical Histories that Preserve the Cold Classical Kuiper Belt". The Astrophysical Journal. 746 (2): 171. arXiv: . Bibcode:2012ApJ...746..171W. doi:10.1088/0004-637X/746/2/171.
- Batygin, Konstantin; Brown, Michael E.; Fraser, Wesley (2011). "Retention of a Primordial Cold Classical Kuiper Belt in an Instability-Driven Model of Solar System Formation". The Astrophysical Journal. 738 (1): 13. arXiv: . Bibcode:2011ApJ...738...13B. doi:10.1088/0004-637X/738/1/13.
- Dones, L.; Weissman, P. R.; Levison, H. F.; Duncan, M. J. (2004). "Oort cloud formation and dynamics". Comets II: 153–174. Bibcode:2004ASPC..323..371D.
- Brasser, R.; Morbidelli, A. (2013). "Oort cloud and Scattered Disc formation during a late dynamical instability in the Solar System". Icarus. 225 (1): 40.49. arXiv: . Bibcode:2013Icar..225...40B. doi:10.1016/j.icarus.2013.03.012.
- 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: . Bibcode:2007AJ....134.1790M. doi:10.1086/521705.
- 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.
- Nesvorný, David; Vokrouhlický, David; Deienno, Rogerio. "Capture of Irregular Satellites at Jupiter". The Astrophysical Journal. 784 (1): 22. arXiv: . Bibcode:2014ApJ...784...22N. doi:10.1088/0004-637X/784/1/22.
- Nesvorný, David. "Young Solar System's Fifth Giant Planet?". The Astrophysical Journal Letters. 742 (2): L22. arXiv: . Bibcode:2011ApJ...742L..22N. doi:10.1088/2041-8205/742/2/L22.
- Siegel, Ethan. "Jupiter May Have Ejected A Planet From Our Solar System". Starts With a Bang. forbes.com. Retrieved 20 December 2015.
- Nesvorný, David (2015). "Jumping Neptune Can Explain the Kuiper Belt Kernel". The Astronomical Journal. 150 (3): 68. arXiv: . Bibcode:2015AJ....150...68N. doi:10.1088/0004-6256/150/3/68.
- 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: . Bibcode:2012AJ....144..117N. doi:10.1088/0004-6256/144/4/117.
- Nesvorný, David (2015). "Evidence for Slow Migration of Neptune from the Inclination Distribution of Kuiper Belt Objects". The Astronomical Journal. 150 (3): 73. arXiv: . Bibcode:2015AJ....150...73N. doi:10.1088/0004-6256/150/3/73.
- Dones, L.; Levison, H. L. "The Impact Rate on Giant Planet Satellites During the Late Heavy Bombardment" (PDF). 44th Lunar and Planetary Science Conference (2013).
- Kaib, Nathan A.; Chambers, John E. (2016). "The fragility of the terrestrial planets during a giant-planet instability". Monthly Notices of the Royal Astronomical Society. 455 (4): 3561–3569. arXiv: . Bibcode:2016MNRAS.455.3561K. doi:10.1093/mnras/stv2554.
|Wikimedia Commons has media related to Nice Model.|