El Modelo de Niza es un escenario teórico analizado a partir de un modelo computacional que permite describir de forma dinámica los procesos que pudieron llevar a la formación y evolución del sistema solar en su configuración actual. Se llama así por la ubicación del Observatorio de la Costa Azul, donde se desarrolló inicialmente a partir del año 2005 en Niza, Francia.[1][2][3]
Investigaciones posteriores revelaron una serie de diferencias entre las predicciones del modelo de Niza original y las observaciones del sistema solar actual, por ejemplo, en las órbitas de los planetas terrestres y los asteroides, lo que llevó a su modificación.
Simulación mostrando los planetas exteriores y el cinturón de planetesimales: a) Configuración inicial, después de que Júpiter y Saturno alcanzaran la resonancia 2:1 ; b) Dispersión de los planetesimales en el sistema solar interior después del cambio orbital de Neptuno (azul oscuro) y Urano (azul claro); c) Tras la eyección de los planetesimales por los planetas.[4]
Descripción
El núcleo original del modelo de Niza es un triplete de artículos publicados en la revista de ciencia general Nature en 2005 gracias a la colaboración de un grupo internacional de científicos: Rodney Gomes, Hal Levison, Alessandro Morbidelli y Kleomenis Tsiganis.[4][5][6] En estas publicaciones, los cuatro autores propusieron que después de la disipación del gas y el polvo del disco del sistema solar primordial, los cuatro planetas gigantes (Júpiter, Saturno, Urano y Neptuno) se encontrarían originalmente en órbitas casi circulares de entre ~5.5 y ~17 unidades astronómicas (UA), espaciados de una forma mucho más compacta que en el presente. Un disco grande y denso de pequeños planetesimales de hielo y roca, con un total de aproximadamente 35 masas terrestres, se extendería desde la órbita del planeta gigante más externo hasta aproximadamente 35 UA.
Los científicos comprenden tan poco acerca de la formación de Urano y Neptuno, que Levison afirma que "... las posibilidades con respecto a la formación de Urano y Neptuno son casi infinitas".[7]
Sin embargo, se sugiere que este sistema planetario evolucionó a partir de la interacción entre el cinturón externo de planetesimales y los planetas gigantes. Así, los planetesimales situados en el borde interno del disco experimentarían ocasionalmente el tirón gravitatorio del planeta gigante más exterior, lo que cambiaría sus órbitas. Los planetas dispersarían a la mayoría de los pequeños cuerpos helados que encontrasen hacia el interior, intercambiando momento angular con los objetos dispersos para que los planetas se movieran hacia afuera en respuesta, preservando el momento angular del sistema. Estos planetesimales se dispersarían de forma similar al encontrar el siguiente planeta, desplazándose hacia el exterior de forma sucesiva las órbitas de Urano, Neptuno y Saturno.[7] El intercambio de impulso se traduciría en un desplazamiento (migración) de los planetesimales. Este proceso continuaría hasta que interactuasen con el planeta gigante más interno y masivo, Júpiter, cuya inmensa gravedad los enviaría a órbitas altamente elípticas o incluso los expulsaría directamente del sistema solar. Esto, en contraste, haría que Júpiter se moviera ligeramente hacia adentro.[8]
La baja tasa de encuentros orbitales gobierna la tasa con la que los planetesimales salen del disco y la tasa de migración correspondiente. Después de varios cientos de millones de años de migración lenta y gradual, Júpiter y Saturno, los dos planetas más grandes, se acoplan mutuamente en resonancia orbital 1:2. Esta resonancia aumenta su excentricidad orbital, desestabilizando todo el sistema planetario. La disposición de los planetas gigantes se alteraría rápida y drásticamente.[9] Júpiter desplazaría a Saturno hacia su posición actual, y esta reubicación provocaría encuentros gravitacionales mutuos entre Saturno y los dos gigantes helados, impulsando a Neptuno y Urano a órbitas mucho más excéntricas. Estos gigantes de hielo posteriormente se introducirían en el disco planetesimal, dispersando decenas de miles de planetesimales de sus órbitas anteriormente estables hacia el sistema solar exterior. Esta interrupción dispersa casi por completo el disco primordial, eliminando el 99% de su masa, un escenario que explica la ausencia actual de una densa población de Objetos transneptunianos.[5] Algunos de los planetesimales se impulsarían hacia el sistema solar interno, produciendo una repentina afluencia de impactos sobre los planetas terrestres: el bombardeo intenso tardío.[4]
Finalmente, los planetas gigantes alcanzarían sus semiejes mayores orbitales actuales, y la fricción dinámica con el disco planetesimal restante amortiguaría sus excentricidades, obligando a que las órbitas de Urano y Neptuno fuesen circulares nuevamente.[10]
En aproximadamente el 50% de los modelos iniciales de Tsiganis y sus colegas, Neptuno y Urano también intercambian lugares.[5] Un intercambio de Urano y Neptuno sería consistente con los modelos de su formación en un disco que tuviese una densidad de superficie que se redujera con la distancia al Sol, lo que conllevaría que las masas de los planetas también deberían disminuir con la distancia desde el Sol.[1]
Ejemplo de simulación del modelo de Niza de la migración de los cuatro planetas gigantes, según su distancia al Sol.
Características del sistema solar
La ejecución de modelos dinámicos con diferentes condiciones iniciales para una duración simulada de la historia del sistema solar producirá diversas poblaciones de objetos en su interior. Como las condiciones iniciales del modelo pueden variar, cada población será más o menos numerosa y tendrá propiedades orbitales particulares. Demostrar rigurosamente la validez de un modelo de la evolución del sistema solar primigenio es difícil, ya que este proceso no se puede observar directamente.[9] Sin embargo, el éxito de cualquier modelo dinámico se puede evaluar comparando las predicciones de población de las simulaciones con las observaciones astronómicas de estas poblaciones.[9] En la actualidad, los modelos informáticos del sistema solar se inician con las condiciones de partida del escenario de Niza, que en la práctica se ha visto que produce resultados coincidentes con muchos de los aspectos observados del sistema solar actual.[11]
El registro de cráteres en la Luna y en los planetas terrestres es parte de la evidencia principal del bombardeo intenso tardío: un incremento en el número de cuerpos de impacto, unos 600 millones de años después de la formación del sistema solar. En el modelo de Niza, los planetesimales helados se dispersan en las órbitas de los cruces de planetas cuando Urano y Neptuno rompen el disco exterior, lo que provoca un aumento agudo de los impactos de los objetos helados. La migración de planetas externos también hace que los movimientos medios asociados a la resonancia secular se propaguen a través del sistema solar interno. En el cinturón de asteroides, aumentan las excentricidades, lo que los conduce a órbitas que se intersecan con las de los planetas terrestres, causando un período prolongado de impactos de objetos rocosos, eliminando aproximadamente el 90 % de su masa.[4] El número de planetesimales que alcanzaría la luna es consistente con el registro de cráteres del bombardeo intenso tardío.[4] Sin embargo, la distribución orbital de los asteroides restantes no coincide con las observaciones.[12] En el sistema solar exterior, los impactos en las lunas de Júpiter son suficientes para desencadenar la diferenciación de Ganimedes pero no la de Calisto.[13] Sin embargo, la cantidad de impactos de planetesimales helados en las lunas internas de Saturno es excesiva, lo que se traduciría en la vaporización de su hielo.[14]
Troyanos y el cinturón de asteroides
Después de que Júpiter y Saturno entrasen en resonancia 2:1, su influencia gravitacional combinada desestabilizaría la región coorbital de los troyanos existentes en los puntos de Lagrange L4 y L5 de Júpiter y de Neptuno, propiciando que se escapasen y a la vez que se capturasen nuevos objetos del disco planetesimal externo.[15]
Los objetos en la región coorbital troyana experimentan una libración que se desplaza cíclicamente en relación con los puntos L4 y L5. Cuando Júpiter y Saturno están cerca pero no en resonancia, la ubicación donde Júpiter pasa Saturno en relación con su perihelio circula lentamente. Si el período de esta circulación cae en resonancia con el período en el que los troyanos pueden aumentar el alcance de sus libraciones, finalmente llegarían a escaparse.[6] Cuando esto ocurriese, la región coorbitaria del troyano quedaría "dinámicamente abierta" y los objetos podrían tanto escapar como entrar en esta zona. Los troyanos primordiales se escaparían, y una fracción de los numerosos objetos del disco planetesimal interrumpido lo habitarían temporalmente.[3] Más adelante, cuando Júpiter y Saturno estuvieran más alejados, la región de los troyanos se "cerraría dinámicamente", y se capturarían numerosos planetesimales en esta región, como se observa en la actualidad.[6] Los troyanos capturados tienen una amplia gama de inclinaciones, que no se habían entendido previamente, debido a sus repetidos encuentros con los planetas gigantes.[3] El ángulo de libración y la excentricidad de la población simulada también coincide con las observaciones de las órbitas de los asteroides troyanos de Júpiter.[6] Este mecanismo del modelo de Niza genera de manera similar los asteroides troyanos de Neptunos.[3]
Un gran número de planetesimales también se habrían capturado en las resonancias de movimiento medias de Júpiter a medida que migraba hacia adentro. Los que permanecieron en una resonancia 3:2 con Júpiter forman el grupo de Hilda. La excentricidad de otros objetos disminuyó mientras se encontraban en resonancia, escapando a órbitas estables en el cinturón de asteroides exterior, a distancias superiores a 2.6 UA cuando las resonancias se extendieron hacia adentro.[16] Estos objetos capturados habrían sufrido erosión por colisión, produciendo fragmentos más pequeños capaces de ser afectados por el efecto Yarkovsky, causando que los objetos pequeños se desvíen hacia resonancias inestables, y por el efecto Poynting-Robertson, provocando que los granos más pequeños se desplacen hacia el Sol. Estos procesos eliminan más del 90 % de la masa de origen implantada en el cinturón de asteroides según Bottke y sus colegas.[17] La distribución de frecuencia de tamaños de esta población simulada después de esta erosión coincide plenamente con las observaciones.[17] Esto sugiere que los troyanos de Júpiter, los asteroides Hildas y algunos de los cinturones de asteroides externos, y todos los asteroides del tipo D espectrales, son los planetesimales remanentes de este proceso de captura y erosión.[17] También se ha sugerido que el planeta enano Ceres se capturó mediante un proceso similar.[18] Han sido descubiertos algunos asteroides del tipo D con semiejes mayores de menos de 2.5 UA, más cercanos que los que se capturarían en el modelo original de Niza.[19]
Satélites del sistema exterior
Cualquier población original de satélites irregulares capturada por los mecanismos tradicionales, como el arrastre o los impactos de los discos de acreción,[20] se perdería durante los encuentros entre los planetas en el momento de la inestabilidad global del sistema.[5] En el modelo de Niza, los planetas exteriores encuentran un gran número de planetesimales después de que Urano y Neptuno entren en el disco planetesimal, interrumpiéndolo. Una fracción de estos planetesimales son capturados por estos planetas a través de la interacción de tres cuerpos durante los encuentros entre planetas. La probabilidad de que un gigante helado capture cualquier planetesimal es relativamente alta, del orden de 10−7.[21] Estos nuevos satélites podrían capturarse en casi cualquier ángulo, por lo que, a diferencia de los satélites regulares de Saturno, Urano y Neptuno, no necesariamente orbitan en los planos ecuatoriales de los planetas. Algunos satélites irregulares pueden incluso haber sido intercambiados entre planetas.
Las órbitas irregulares resultantes coinciden bien con los ejes, inclinaciones y excentricidades de los semiejes mayores de las poblaciones observadas.[21] Se sospecha que las colisiones posteriores entre estos satélites capturados pueden haber creado las familias de colisión observadas en la actualidad.[22] Estas colisiones también son necesarias para erosionar la población hasta la distribución de tamaño actual.[23]
La configuración de Tritón, la luna más grande de Neptuno, puede explicarse si se supone que se capturó en una interacción de tres cuerpos que involucra la interrupción de un planetoide binario.[24] Tal interrupción binaria sería más probable si Tritón fuera el miembro más pequeño del sistema binario.[25] Sin embargo, la captura de Tritón sería más probable en el sistema solar primoigenio, cuando el disco de gas amortiguaría las velocidades relativas, y las reacciones de intercambio binario en general no habrían suministrado la gran cantidad de pequeños satélites irregulares.[25]
No hubo suficientes interacciones entre Júpiter y los otros planetas para explicar el séquito de irregulares de Júpiter en las simulaciones iniciales del modelo de Niza que reproducían otros aspectos del Sistema Solar exterior. Esto sugiere que un segundo mecanismo estaba funcionando para ese planeta o que las simulaciones iniciales no reproducían la evolución de las órbitas de los planetas gigantes.[21]
Formación del cinturón de Kuiper
La migración de los planetas exteriores también es necesaria para explicar la existencia y las propiedades de las regiones más externas del sistema solar.[10] Originalmente, el cinturón de Kuiper era mucho más denso y estaba más cercano al Sol, con un borde exterior de aproximadamente 30 UA. Su borde interior habría estado más allá de las órbitas de Urano y Neptuno, que a su vez estaban mucho más cerca del Sol cuando se formaron (muy probablemente en el rango de 15 a 20 UA), y en ubicaciones opuestas, con Urano más alejado del Sol que Neptuno.[4][10]
Los encuentros gravitacionales entre los planetas desplazarían a Neptuno hacia el disco planetesimal, con un eje semi mayor de ~ 28 UA y una excentricidad tan alta como 0.4. La gran excentricidad de Neptuno hace que sus resonancias de movimiento medio se superpongan y las órbitas en la región entre Neptuno y sus resonancias de movimiento medio de 2:1 se vuelven caóticas. Las órbitas de los objetos entre Neptuno y el borde del disco planetesimal en este momento pueden evolucionar hacia afuera en órbitas estables de baja excentricidad dentro de esta región. Cuando la excentricidad de Neptuno se ve amortiguada por la fricción dinámica, quedan atrapadas en estas órbitas. Estos objetos forman un cinturón dinámico poco activo, ya que sus inclinaciones permanecen pequeñas durante el corto tiempo que interactúan con Neptuno. Más tarde, a medida que Neptuno migra hacia el exterior en una órbita de baja excentricidad, los objetos que se han dispersado hacia fuera se capturan en sus resonancias y pueden disminuir sus excentricidades y aumentar sus inclinaciones debido al mecanismo de Kozai, lo que les permite escapar a órbitas estables de mayor inclinación. Otros objetos quedan capturados en resonancia, formando los plutinos y otras poblaciones resonantes. Estas dos poblaciones son dinámicamente calientes, con mayores inclinaciones y excentricidades; debido a que están dispersos hacia afuera y a un período más largo, estos objetos interactúan con Neptuno.[10]
Esta evolución de la órbita de Neptuno produce poblaciones tanto resonantes como no resonantes, un borde exterior en la resonancia 2:1 de Neptuno y una pequeña masa en relación con el disco planetesimal original. El exceso de plutinos de baja inclinación en otros modelos se evita debido a que Neptuno se dispersa hacia afuera, dejando su resonancia 3:2 más allá del borde original del disco planetesimal. Las diferentes ubicaciones iniciales, con los objetos clásicos fríos originarios principalmente del disco externo, y los procesos de captura, ofrecen explicaciones para la distribución de inclinación bimodal y su correlación con las composiciones.[10] Sin embargo, esta evolución de la órbita de Neptuno no tiene en cuenta algunas de las características de la distribución orbital. Predice una excentricidad promedio mayor en las órbitas de objetos del cinturón de Kuiper clásico de lo que se observa (0.10–0.13 contrs 0.07) y no produce suficientes objetos con mayor inclinación. Tampoco puede explicar la aparente ausencia completa de objetos grises en la población fría, aunque se ha sugerido que las diferencias de color surgen en parte de los procesos de evolución de la superficie y no completamente de las diferencias en la composición primordial.[26]
La escasez de los objetos de menor excentricidad predichos en el modelo de Niza puede indicar que la población fría se formó in situ. Además de sus diferentes órbitas, las poblaciones frías y calientes tienen diferentes colores. La población fría es marcadamente más roja que la caliente, lo que sugiere que tiene una composición diferente y está formada en una región diferente.[26][27] La población fría también incluye un gran número de objetos binarios con órbitas poco ligadas que probablemente no sobrevivirían en un encuentro cercano con Neptuno.[28] Si la población fría se formara en su ubicación actual, preservarla requeriría que la excentricidad de Neptuno se mantuviera en valores reducidos,[29] o que su perihelio se eliminara rápidamente debido a una fuerte interacción entre esta y Urano.[30]
Discos dispersos y la nube de Oort
Los objetos dispersados hacia afuera por Neptuno en órbitas con un semi eje mayor que 50 UA pueden ser capturados en resonancias que forman la población resonante del disco disperso, o si sus excentricidades se reducen mientras están en resonancia, entonces pueden escapar de la resonancia a órbitas estables en el disco disperso mientras Neptuno está migrando. Cuando la excentricidad de Neptuno es grande, su afelio puede llegar mucho más allá de su órbita actual. Los objetos que alcanzan el perihelio cerca o por encima que el de Neptuno en este momento pueden desprenderse del planeta cuando se modera su excentricidad, reduciendo su afelio, dejándolos en órbitas estables en el disco disperso.[10]
Los objetos dispersados hacia afuera por Urano y Neptuno en órbitas más grandes (aproximadamente 5000 UA) pueden elevar su perihelio por la marea galáctica separándolos de la influencia de los planetas que forman la nube de Oort interna con inclinaciones moderadas. Otros que alcanzan órbitas aún más grandes pueden ser perturbados por estrellas cercanas que forman la nube externa de Oort con inclinaciones isotrópicas. Los objetos dispersados por Júpiter y Saturno normalmente se expulsan fuera del sistema solar.[31] Varios porcentajes del disco planetesimal inicial pueden depositarse en estos reservorios.[32]
El modelo de Niza ha experimentado varias modificaciones desde su publicación inicial, ya que la comprensión de la formación del sistema solar ha avanzado y se han identificado diferencias significativas entre sus predicciones y observaciones. Los modelos hidrodinámicos del sistema solar temprano indican que las órbitas de los planetas gigantes convergen, dando lugar a su captura en resonancias.[33] Durante el acercamiento lento de Júpiter y Saturno a la resonancia 2:1, Marte pudo capturarse en una resonancia secular, excitando su excentricidad a un nivel que desestabilizara el sistema solar interior. Las excentricidades de los otros planetas terrestres también pueden excitarse más allá de los niveles actuales, barriéndose las resonancias seculares después de la inestabilidad.[34] La distribución orbital del cinturón de asteroides también deja un exceso de objetos de alta inclinación debido a las resonancias seculares que excitan las inclinaciones y eliminan los objetos de baja inclinación.[12] Otras diferencias entre las predicciones y las observaciones incluyeron la captura de algunos satélites irregulares por Júpiter, la vaporización del hielo de las lunas internas de Saturno, la escasez de objetos de alta inclinación capturados en el cinturón de Kuiper y el reciente descubrimiento de asteroides del tipo D en el cinturón interior de asteroides.
Las primeras modificaciones al modelo de Niza fueron las posiciones iniciales de los planetas gigantes. Las investigaciones sobre el comportamiento de los planetas que orbitan en un disco de gas utilizando modelos hidrodinámicos revelan que los planetas gigantes migrarían hacia el Sol. Si la migración continuara, habría dado como resultado que Júpiter orbitaría cerca del Sol como los exoplanetas conocidos como Júpiteres calientes. Sin embargo, la captura de Saturno en una resonancia con Júpiter evitaría esta situación, y la captura posterior de los otros planetas da como resultado una configuración resonante cuádruple con Júpiter y Saturno en una resonancia 3:2.[33] Una inestabilidad tardía a partir de esta configuración es posible si el disco exterior contiene objetos con masa similar a la de Plutón. La agitación gravitatoria del disco planetesimal externo por estos objetos de la masa de Plutón aumentaría sus excentricidades y también contribuiría a la migración hacia el interior de los planetas gigantes. La resonancia cuádruple de los planetas gigantes se rompe cuando las resonancias seculares se cruzan durante la migración hacia el interior. A continuación sigue una inestabilidad tardía similar a la del modelo original de Niza. Sin embargo, contradiciendo al modelo inicial, el momento de esta inestabilidad no es sensible a la distancia entre el planeta exterior y el disco planetesimal. La combinación de las órbitas planetarias resonantes y la inestabilidad tardía provocada por estas interacciones de larga distancia se han denominado Modelo de Niza 2.[35]
La segunda modificación fue el requisito de que uno de los gigantes de hielo se encuentre con Júpiter, lo que provocaría que su semieje mayor se modificara bruscamente. En este escenario de Júpiter saltador, un gigante de hielo se encontraría con Saturno y se dispersaría hacia el interior en una órbita que cruzase la de Júpiter, lo que a su vez haría que la órbita de Saturno se expandiera; luego se encontraría con Júpiter y se dispersaría hacia afuera, lo que provocaría que la órbita de Júpiter se redujera. Esto da como resultado una separación gradual de las órbitas de Júpiter y Saturno en lugar de una suave migración divergente.[34] La separación gradual de las órbitas de Júpiter y Saturno evita el barrido lento de las resonancias seculares a través del Sistema Solar interior que pudo producir el incremento de las excentricidades de las órbitas de los planetas terrestres[34] y un cinturón de asteroides con una proporción excesiva de objetos de inclinación de alta a baja.[12] Los encuentros entre el gigante de hielo y Júpiter en este modelo permiten a Júpiter adquirir sus propios satélites irregulares.[36]
Los troyanos de Júpiter también se capturan después de estos encuentros, cuando su semieje mayor salta y, si el gigante de hielo pasa a través de uno de los puntos de libración que dispersan los troyanos, una población se agota en relación con la otra.[37] El desplazamiento más rápido de las resonancias seculares a través del cinturón de asteroides limita la pérdida de estos de su núcleo. La mayoría de los impactadores rocosos del bombardeo intenso tardío en su lugar se originan en una extensión interna que se interrumpe cuando los planetas gigantes alcanzan sus posiciones actuales, quedando un remanente como asteroides de Hungaria.[38] El modelo contempla que algunos asteroides del tipo D están incrustados en el cinturón interior, dentro de 2.5 UA, cuando es cruzado por el gigante de hielo durante sus encuentros.[39]
La frecuente expulsión del gigante de hielo al encontrarse con Júpiter, ha llevado a David Nesvorný y otros a la hipótesis de un sistema solar temprano con cinco planetas gigantes, uno de los cuales fue expulsado durante un episodio de inestabilidad.[40][41] Este modelo de Niza de cinco planetas comienza con los planetas gigantes en una cadena de resonancias 3:2, 3:2, 2:1 y 3:2, con un disco planetesimal en órbita más allá de ellos.[42] Tras la ruptura de la cadena de resonancias, Neptuno migraría primero hacia el interior del disco planetesimal, alcanzando 28 UA antes de que comiencen los encuentros entre planetas.[43] Esta migración inicial reduce la masa del disco externo, lo que permite preservar la excentricidad de Júpiter[44] y produce un cinturón de Kuiper con una distribución de inclinación que coincide con las observaciones si 20 masas terrestres permanecieran en el disco planetesimal cuando comenzó la migración.[45]
La excentricidad de Neptuno puede permanecer pequeña durante la inestabilidad, ya que solo se encuentra con el gigante de hielo expulsado, lo que permite que se conserve un cinturón frío clásico in situ.[43] El cinturón planetesimal de masa inferior en combinación con el incremento de inclinaciones y excentricidades de los objetos del tamaño de Plutón también reduce significativamente la pérdida de hielo por las lunas internas de Saturno.[46] La combinación de una ruptura tardía de la cadena de resonancia y una migración de Neptuno a 28 UA antes de la inestabilidad es improbable con el modelo de Niza 2. Esta brecha se puede salvar mediante una migración lenta impulsada por el polvo durante varios millones de años después de un escape temprano de la resonancia.[47]
Un estudio reciente encontró que el modelo de Niza de cinco planetas tiene una probabilidad estadísticamente pequeña de reproducir las órbitas de los planetas terrestres. Aunque esto implica que la inestabilidad ocurrió antes de la formación de los planetas terrestres y no pudo ser la fuente del bombardeo intenso tardío,[48][49] la ventaja de una inestabilidad temprana se reduce por los saltos considerables en el semieje mayor de Júpiter y de Saturno necesarios para preservar el cinturón de asteroides.[50][51]
↑ abG. Jeffrey Taylor (21 de agosto de 2001). «Uranus, Neptune, and the Mountains of the Moon». Planetary Science Research Discoveries. Hawaii Institute of Geophysics & Planetology. Consultado el 1 de febrero de 2008.
↑Joseph M. Hahn; Malhotra, Renu (13 de julio de 2005). «Neptune's Migration into a Stirred–Up Kuiper Belt: A Detailed Comparison of Simulations to Observations». Astronomical Journal130 (5): 2392-2414. Bibcode:2005AJ....130.2392H. arXiv:astro-ph/0507319. doi:10.1086/452638.
↑ abcdefLevison, Harold F.; Morbidelli, Alessandro; Van Laerhoven, Christa; Gomes, Rodney S.; Tsiganis, Kleomenis (2007). «Origin of the Structure of the Kuiper Belt during a Dynamical Instability in the Orbits of Uranus and Neptune». Icarus196 (1): 258-273. Bibcode:2008Icar..196..258L. arXiv:0712.0553. 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». Early Solar System Impact Bombardment conference (2008). Consultado el 18 de octubre de 2008.
↑ abcMorbidelli, 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 Journal140 (5): 1391-1501. Bibcode:2010AJ....140.1391M. arXiv:1009.1521. doi:10.1088/0004-6256/140/5/1391.
↑Levison, Harold F.; Shoemaker, Eugene M.; Shoemaker, Carolyn S. (1997). «Dynamical evolution of Jupiter's Trojan asteroids». Nature385 (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». Nature460 (7253): 364-366. Bibcode:2009Natur.460..364L. PMID19606143. doi:10.1038/nature08094.
↑ abcBottke, 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 Society40: 464. Bibcode:2008DPS....40.3803M.
↑ abMorbidelli, 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 Journal134 (5): 1790-1798. Bibcode:2007AJ....134.1790M. arXiv:0706.1713. doi:10.1086/521705.
↑ abcBrasser, 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 Astrophysics507 (2): 1053-1065. Bibcode:2009A&A...507.1053B. arXiv:0909.1891. doi:10.1051/0004-6361/200912878.
↑Batygin, Konstantin; Brown, Michael E.; Betts, Hayden (2012). «Instability-driven Dynamical Evolution Model of a Primordially Five-planet Outer Solar System». The Astrophysical Journal Letters744 (1): L3. Bibcode:2012ApJ...744L...3B. arXiv:1111.3682. doi:10.1088/2041-8205/744/1/L3.
↑Deienno, Rogerio; Morbidelli, Alessandro; Gomes, Rodney S.; Nesvorny, David (2017). «Constraining the giant planets' initial configuration from their evolution: implications for the timing of the planetary instability». The Astronomical Journal153: 153. Bibcode:2017AJ....153..153D. arXiv:1702.02094. doi:10.3847/1538-3881/aa5eaa.