Turbulent infall onto class 0 disks as cause of CAI brief condensation episode in the solar system

Resumen Técnico: Aporte Turbulento a Discos de Clase 0 como Causa del Breve Episodio de Condensación de CAI en el Sistema Solar

Enunciado del Problema
Las inclusiones ricas en calcio y aluminio (CAI) encontradas en meteoritos condritas carbonáceas representan los sólidos más antiguos del sistema solar, caracterizados por un episodio de condensación notablemente breve que duró aproximadamente $10^5$ años. Un desafío central en la teoría de la formación del sistema solar es conciliar esta corta escala temporal de formación con la contaminación de larga duración (millones de años) de las regiones de formación estelar por radionúclidos de vida corta (RLC) provenientes de fuentes externas como supernovas o estrellas de la rama asintótica gigante. Aunque modelos anteriores, como el modelo de viento en X, intentaron explicar la formación de CAI mediante vientos magnetocentrífugos, la datación isotópica precisa ha restringido la formación de CAI a una ventana estrecha que desafía los modelos que dependen de una irradiación solar prolongada. Además, aunque se han propuesto diversos mecanismos de calentamiento (por ejemplo, erupciones de FU Orionis, supernovas cercanas), a menudo luchan por explicar simultáneamente la extensión espacial (varias UA), la naturaleza transitoria y la intensidad del calentamiento requeridos para sublimar y re-condensar granos refractarios. Los autores investigan si el aporte dinámico de gas a discos protostelares de Clase 0 puede generar las condiciones térmicas necesarias para este breve episodio de alta temperatura.

Metodología
El estudio emplea simulaciones hidrodinámicas de alta resolución utilizando el código GIZMO en modo de Masa Finita sin Malla (MFM). Este esquema lagrangiano se selecciona por su superior captura de choques, disipación numérica reducida y excelente conservación del momento angular, lo cual es crítico para modelar discos torcidos y aporte turbulento.

• Condiciones Iniciales: Las simulaciones utilizan condiciones iniciales derivadas de un estudio de síntesis de poblaciones de formación de cúmulos estelares (Bate 2012, 2019). Los autores rastrean la evolución de estrellas de masa solar hasta la etapa de Clase 0, capturando historias de aporte realistas, dependientes del tiempo y anisotrópicas desde corrientes de acreción filamentarias.
• Configuración de la Simulación: Se realizan dos simulaciones principales:
  • Ejecución 1: Presenta una fase de aporte rápida y caótica con corrientes de momento angular significativamente desalineadas impactando el disco.
  • Ejecución 2: Un caso de control con una configuración inicial similar pero sin aporte significativo de gas vertical desde corrientes, sirviendo como línea base para la evolución de un disco aislado.
• Modelos Físicos:
  • Transferencia Radiativa: En lugar de una ecuación de estado barotrópica simple, los autores emplean un método de enfriamiento radiativo aproximado (Young et al. 2024) que hibrida estimadores de densidad de columna (Stamatellos et al. 2007; Lombardi Jr et al. 2015) para manejar tanto geometrías esféricas como de disco. Esto permite una evolución térmica realista donde el calentamiento viscoso domina sobre la irradiación en el plano medio ópticamente grueso.
  • Resolución: Un estudio de resolución (Apéndice A) confirma que las características físicas clave, específicamente las ondas de densidad espiral y el calentamiento por choque asociado, convergen en conteos de partículas de $\sim 10$ millones. Las simulaciones finales utilizan 10 millones de partículas.

Resultados Clave
Las simulaciones revelan una evolución distintiva de tres fases en la Ejecución 1 (el escenario de aporte), contrastando marcadamente con la evolución estable de la Ejecución 2.

1. Evolución Global y Reorientación del Disco:

   • Ejecución 1: El sistema experimenta una "fase principal de calentamiento" ($t \approx 6000-9000$ años) impulsada por un aporte vigoroso y desalineado. La interacción entre las corrientes de aporte y el disco rotatorio causa una reorientación global del disco de casi $30^\circ$. Esto conduce a la formación de torceduras y ondas de densidad espiral globales.
   • Ejecución 2: Sin aporte externo significativo, el disco permanece geométricamente delgado y estable, con solo leves inclinaciones ($10^\circ-20^\circ$) y sin reorientación a gran escala.

2. Evolución Térmica:

   • Ejecución 1: La fase principal de calentamiento se caracteriza por un rápido aumento de temperatura. La temperatura máxima del disco alcanza $\sim 2000$ K, suficiente para sublimar granos de CAI. Este calentamiento se mantiene durante varios miles de años. El mecanismo de calentamiento se identifica como disipación de choque en la interfaz disco-envoltura y dentro del plano medio del disco debido a choques espirales. El disco denso y verticalmente engrosado durante el pico de aporte reduce la eficiencia del enfriamiento radiativo, atrapando el calor.
   • Ejecución 2: La temperatura del disco permanece significativamente más baja, alcanzando un máximo solo alrededor de $400-600$ K. La ausencia de aporte externo impide la generación del intenso calentamiento por choque observado en la Ejecución 1.

3. Mecanismos de Calentamiento:

   • Choques Superficiales: Cortes verticales muestran flujos de aporte convergentes colisionando con la superficie del disco, creando fuertes choques y zonas de alta entropía.
   • Choques Espirales: El análisis del plano medio revela que las ondas de densidad espiral (impulsadas por inestabilidad gravitacional con Toomre $Q \approx 1$) generan saltos de entropía localizados y calentamiento. En el disco interior ($r \lesssim 3$ UA), las temperaturas alcanzan $\sim 2000$ K, mientras que los choques espirales en el rango de 2–3 UA crean puntos calientes que superan los 1500 K.
   • Momento Angular: La desalineación de los vectores de momento angular de aporte impulsa la reorientación y la turbulencia, lo que a su vez aumenta la disipación viscosa y la formación de choques.

Significado y Afirmaciones
El artículo afirma que el aporte turbulento a discos de Clase 0 proporciona un mecanismo viable para el breve episodio de condensación de alta temperatura requerido para formar CAI. Específicamente:

• Calentamiento Transitorio: El calentamiento impulsado por el aporte es transitorio, durando $\sim 10^3-10^4$ años, lo cual se alinea con las restricciones de edad radiogénica de las CAI.
• Extensión Espacial: El calentamiento se extiende hasta $\sim 2-3$ UA, cubriendo la región donde se cree que se forman y preservan las CAI.
• Re-condensación: El modelo sugiere un ciclo donde los granos cargados con CAI "pre-solares" en nubes de aporte son sublimados por el calentamiento por choque y posteriormente re-condensan a medida que el sistema se enfría durante la fase de subsistencia del aporte.
• Dicotomía de Sistemas: La comparación entre la Ejecución 1 y la Ejecución 2 implica que la formación de un reservorio de CAI de población única no es un resultado inevitable de toda formación de discos. Depende de condiciones de contorno específicas, como la presencia de corrientes de aporte caóticas y desalineadas. Esto puede explicar por qué las CAI son abundantes en condritas carbonáceas (CC) pero escasas en condritas no carbonáceas (NC), sugiriendo que solo una fracción de los sistemas planetarios experimenta la dinámica de aporte específica necesaria para reiniciar globalmente el reloj de las CAI.

Los autores concluyen que, aunque sus simulaciones cubren la etapa muy temprana de la formación del disco (Clase 0, masa central $\sim 0.2 M_\odot$), el mecanismo de calentamiento impulsado por el aporte ofrece una explicación robusta para las condiciones térmicas necesarias para la formación de CAI, vinculando directamente la historia dinámica de la envoltura protostelar con el registro meteorítico.