Formación de los planetas gigantes del Sistema Solar
Authors:   , Guilera Octavio Miguel 
Published: August 2019

El estudio de los sistemas planetarios es uno de los temas fundamentales de las ciencias astronómicas. El interés en nuestro Sistema Solar surgió en tiempos remotos. Hace casi dos décadas, Mayor y Queloz (1995) detectaron el primer planeta extrasolar en órbita alrededor de una estrella de tipo solar. Este notable descubrimiento ha hecho surgir desde entonces un enorme interés en el estudio de sistemas planetarios en general. Esta Tesis es una continuación natural de un proyecto de investigación que empezo hace más de diez años en el Grupo de Ciencias Planetarias del cual formo parte. En al año 2005, Benvenuto y Brunini desarrollaron un código para el cálculo de la formación de un planeta gigante basado en los códigos estandar de formación estelar. A partir de este nuevo código, la Dra. Andrea Fortier desarrolló su Tesis de Doctorado, profundizando sobre los fenómenos físicos que dan lugar a la formación de un planeta gigante y mejorando el código. Esta Tesis continúa esta línea de investigación enfocándose principalmente en la formación de los planetas gigantes del Sistema Solar. El objetivo de la misma es intentar contribuir a la comprensión de los procesos que dan lugar a la formación de sistemas planetarios como un todo. Se trabajó en la construcción de modelos de crecimiento planetario en el contexto de un disco protoplanetario realista que evoluciona en el tiempo. Específicamente, se estudió la interacción que surge entre dos (o más) embriones que crecen simultáneamente en un mismo disco protoplanetario. En relación a la formación de los planetas gigantes, actualmente se consideran dos modelos antagónicos para explicar la existencia de estos objetos, el modelo de inestabilidad gravitatoria y el modelo de inestabilidad nucleada o modelo de acreción del núcleo. El modelo de inestabilidad gravitatoria propone la formación de los planetas gigantes como consecuencia de inestabilidades hidrodinámicas en el disco protoplanetario. Estas inestabilidades producirían el colapso gravitatorio de una porción de la componente gaseosa del disco protoplanetario dando lugar a la formación de objetos con masas subestelares del orden de las de los planetas gigantes (la masa de Júpiter es un milésimo de la masa del Sol) en una escala de tiempo muy corta, de algunos miles de años. Esta corta escala de tiempo en la que se produce la formación planetaria es considerada la principal virtud de la teoría de inestabilidad gravitatoria. Por otro lado, el modelo de inestabilidad nucleada, actualmente el más aceptado por la comunidad científica, propone que el planeta comienza a formarse a través de la acreción de planetesimales (bloques fundamentales en el proceso de formación planetaria, con tamaños que van desde el metro a centenas de kilómetros). Inicialmente, el planeta tiene una masa pequeña, del orden de la masa de la Luna (aproximadamente un centésimo de la masa de la Tierra). El planeta aumenta su masa a expensas de los planetesimales, y al estar inmerso en un disco con una componente gaseosa, poco a poco comienza a ligar el gas circundante generando una envoltura gaseosa, la cual inicialmentre tiene una masa varios órdenes de magnitud menor que la del núcleo (llamaremos núcleo a la componente sólida, o de alta densidad, del planeta). Cuando el núcleo alcanza una masa del orden de diez veces la masa de la Tierra, se produce la inestabilidad nucleada. Las capas de la envoltura gaseosa ya no pueden ser sostenidas en equilbrio hidrostático, y se produce el colapso de la envoltura sobre el núcleo. De esta manera, el planeta acreta la mayor parte del gas que lo compone alcanzando su masa final en una escala de tiempo muy corta (de unos miles de años). La duración total del proceso es de algunos millones de años, la cual se suele citar como más larga que la escala de tiempo en que se ha observado que subsisten las nebulosas protoplanetarias (Mamajek, 2009). Esto constituye una evidente paradoja que tradicionalmente se ha planteado como la principal dificultad del modelo de inestabilidad nucleada. Sin embargo, una de las principales razones por las que se prefiere a este modelo frente al de inestabilidad gravitatoria está relacionada con la masa de los núcleos que predice. Los estudios, tanto teóricos como observacionales (estudio del pasaje de satélites artificiales en las proximidades de los planetas de nuestro Sistema Solar, estudio de la dinámica de los satélites naturales de los planetas gigantes del Sistema Solar) predicen que Júpiter, Saturno, Urano y Neptuno tienen núcleos con masas del orden de una decena de veces la terrestre (Podolak et al., 2000; Saumon y Guillot, 2004; Guillot, 2005), tal como predice la teoría de inestabilidad nucleada. En referencia a la evolución del disco protoplanetario, éste sufre una evolución de gran complejidad, la cual es necesario modelar en forma detallada ya que afecta la capacidad de crecimiento de los planetas inmersos en el mismo. Además, el mismo disco es el sistema físico a través del cual se produce la interación planeta – planeta. Aquí no nos referimos simplemente a la interacción gravitatoria sino a la modificación de las poblaciones de planetesimales como consecuencia de la presencia de varias masas planetarias. Dichas masas fuerzan la migración de planetesimales modificando su densidad superficial. Ésta, a su vez, es la que alimenta a los planetas restantes. Por lo tanto un planeta afecta la disponibilidad de materia de la que podrían alimentarse los planetas restantes de un sistema en formación. Cabe destacar que los trabajos de Fortier et al. (2007, 2009) son los primeros en los que se considera el régimen de acreción de planetesimales conocido como de crecimiento oligárquico (Ida y Makino, 1993; Kokubo e Ida, 1998, 2000, 2002). Las simulaciones numéricas más detalladas muestran que el régimen de acreción según el cual los planetas sufren la mayor parte de su crecimiento es el crecimiento oligárquico. Según este régimen, el planeta, luego de alcanzar una masa algo menor a la lunar, es capaz de perturbar su entorno de forma tal que su crecimiento se autolimita. De esta forma el proceso de formación planetaria se vuelve más lento que el predicho por el crecimiento rápido conocido como crecimiento en fuga tal como el considerado por Pollack et al. (1996) u otros autores (Hubicky et al., 2005; Alibert et al., 2005; Dodson-Robinson et al., 2009; Mordacini et al., 2009). Sin embargo, en todos estos trabajos se considera la formación aislada de cada planeta, en donde los posibles efectos que un embrión planetario en formación podría ejercer sobre otros, que crecen simultáneamente en el mismo disco protoplanetario, son despreciados. Esta configuración, la más sencilla posible, es poco realista y resulta insuficiente para comprender de manera más global la formación de los planetas gigantes del Sistema Solar. En esta Tesis desarrollamos un código, en base a los trabajos de Fortier et al. (2007, 2009), en donde se calcula por primera vez –para este tipo de modelos– la formación simultánea de un numéro arbitrario de planetas gigantes inmersos en un disco protoplanetario en evolución. Aplicamos el nuevo modelo para calcular cómo la formación aislada de Júpiter y Saturno se modifica cuando se considera que ambos planetas se forman simultáneamente (Guilera et al., 2010). En este trabajo mostramos que la formación aislada de un planeta gigante puede sufrir cambios significativos cuando la formación del mismo se produce simultáneamente en presencia de otros embriones. En lo que respecta al Sistema Solar, el Modelo de Niza (una triología de trabajos que lleva su nombre debido a que sus autores lo desarrollaron en el Observatorio de la ciudad de Niza, Francia: Tsiganis et al., 2005; Gomes et al., 2005; Morbidelli et al., 2005) cambió el paradigma acerca de su formación. La configuración inicial de este modelo representa la configuración inicial del Sistema Solar exterior cuando la nebulosa primordial fue disipada, y propone que los planetas gigantes estaban en una configuración orbital mucho más compacta que la actual. Particularmente, el modelo propone que los planetas gigantes del Sistema Solar, una vez disipada la nebulosa primordial, se encontraban en órbitas circulares y coplanares, entre ~5.5 UA y ~14 UA (Unidad Astronómica: representa la distancia media entre la Tierra y el Sol, ~150 millones de km). Otro aspecto importante en el Modelo de Niza es la existencia de un disco residual de planetesimales detras de las órbitas de los planetas gigantes. Este disco de planetesimales interactuaría gravitatoriamente con los planetas gigantes y causaría la migración de los mismos a sus posiciones actuales. El éxito del Modelo de Niza radica en que el mismo puede explicar cuantitatívamente muchos aspectos actuales del Sistema Solar: las órbitas, excentricidades e inclinaciones de los planetas gigantes del Sistema Solar (Tsiganis et al., 2005); la existencia de los Troyanos de Júpiter (Morbidelli et al., 2005); el origen del Gran Bombardeo Tardío del Sistema Solar (Gomes et al., 2005) y la formación de la Región Transneptuniana (Levison et al., 2008). Sin embargo, todos estos estudios mencionados consideran que los planetas gigantes del Sistema Solar ya estaban formados, y no plantean discusión alguna acerca de la formación de los mismos. El primero en investigar este problema fue Desch (2007). Utilizando la configuración inicial propuesta por el Modelo de Niza, Desch recalculó la nebulosa solar mínima y considerando una población de planetesimales de 100 m de radio y el régimen de crecimiento oligárquico para los embriones, estimó de manera simple el tiempo de formación de los núcleos de los planetas gigantes. Desch encontró que los mismos podrían formarse en una escala de tiempo compatible con la vida media observada de los discos protoplanetarios. Continuando esta idea, Benvenuto et al. (2009) calcularon de manera detallada la formación aislada de los cuatro planetas gigantes del Sistema Solar. Adoptando la configuración inicial propuesta por el Modelo de Niza, el nuevo modelo de nebulosa solar calculado por Desch y considerando una distribución de tamaños para los planetesimales, encontraron que los planetas gigantes del Sistema Solar pueden formarse en escalas de tiempo compatibles con la vida media observada de los discos protoplanetarios. Además, encontraron que las masas de los núcleos de los cuatros planetas están en perfecto acuerdo con las estimaciones teóricas y observacionales predichas para dichos cuerpos. No obstante, tanto Desch como Benvenuto et al. no tuvieron en cuenta un fenómeno muy importante: la migración de los planetesimales debido a la fricción gaseosa generada por el gas nebular. Thommes et al. (2003), Chambers (2006) y Brunini y Benvenuto (2008) mostraron que este fenómeno tiene una fuerte influencia en las escalas de tiempo de acreción, especialmente para los planetesimales pequeños (menores a 1 km de radio). Este fenómeno también es introducido por primera vez –para este tipo de modelos– en nuestro código, con el cual, continuando el trabajo de Benvenuto et al. (2009), calculamos por primera vez la formación simultánea de los planetas gigantes del Sistema Solar (Guilera et al., 2011). En este trabajo encontramos que el modelo de nebulosa solar propuesto por Desch no favorece la formación simultánea de los planetas gigantes del Sistemas Solar. Sin embargo, modelos de discos con perfiles de densidades superficiales más suaves (como los propuestos por los modelos de discos de acreción) si lo hacen. Aun más, para estos discos, encontramos que si la mayor parte de la masa del sistema está distribuida en planetesimales con radios menores o iguales a 1 km, la formación de los cuatro planetas gigantes puede darse en escalas de tiempo similares (siempre compatibles con la escala de tiempo de vida media observada para los discos protoplanetarios, condición necesaria que debe satisfacer todo modelo de formación de planetas gigantes). En este trabajo, mostramos que la evolución de la población de planetesimales juega un papel importante en el proceso de acreción y formación planetaria. En este estudio, dicha población evoluciona solo por la acreción de los embriones inmersos en el disco y por la migración debida al gas nebular. Sin embargo, la evolución de la población de planetesimales es un fenómeno complejo, y otros procesos pueden tener implicancias significativas, como por ejemplo la evolución colisional, la dispersión o la apertura de brechas en el disco de planetesimales (procesos complejos de incorporar en los modelos como en los que esta Tesis se basa). En lo que respecta a la evolución colisional de la población de planetesimales, a medida que los embriones crecen, debido a las excitaciones gravitatorias que producen, incrementan las velocidades relativas de los planetesimales. Este aumento en las velocidades relativas de los planetesimales causa la fragmentación de los mismos debido a las colisiones mutuas. Después de sucesivas colisiones destructivas los planetesimales van reduciendo sus tamaños. Inaba et al. (2003) y Kobayashi et al. (2010, 2011) encontraron que grandes cantidades de masa pueden perderse por la migración, debido a la fricción gaseosa del gas nebular, de los fragmentos pequeños productos de las colisiones entre planetesimales. Por otro lado, a medida que los embriones crecen, éstos comienzan a ligar el gas circundante. Inicialmente, estas envolturas son poco masivas pero relativamente bastante extendidas, y producen una fricción sobre los planetesimales, aumentando notablemente la sección eficaz de captura de los planetas. Los planetesimales más chicos de la distribución son quienes más sufren ambos efectos. Resulta entonces importante estudiar en forma detallada si la fragmentación de planetesimales y la generación de fragmentos pequeños favorece o inhibe la formación de un planeta gigante. La última etapa de esta Tesis estuvo basada en la generación de un modelo de fragmentación para incorporar en nuestro modelo global de formación planetaria. De esta manera, la población de planetesimales del disco protoplanetario evoluciona ahora por acreción de los embriones, migración orbital y fragmentación. Nuestros estudios (Guilera et al., 2014) muestran que el proceso de fragmentación de planetesimales inhibe fuertemente la formación de los planetas gigantes en un amplio rango de masas para los discos y para un amplio rango de tamaños para los planetesimales. Sin embargo, si la mayor parte de la masa que se pierde en las colisiones entre planetesimales se distribuye en los fragmentos más grandes producto de las mismas, el proceso de formación planetaria se ve favorecido siempre y cuando se considere una distrinbución inicial de planetesimales grandes (con radios del orden de 100 km). Finalmente, encontramos que en este caso, para planetesimales con tamaños menores o iguales a 10 km de radio, es imprescindible considerar un modelo más general en donde además de la fragmentación se tenga en cuenta la coagulación entre planetesimales.