Planet formation at the inner edge of the dead zone -- I. the interplay between accretion outbursts and dust growth

Título: Formación de planetas en el borde interior de la zona muerta: I. La interacción entre erupciones de acreción y el crecimiento de polvo.

Contexto y Problema:
Las observaciones de objetos estelares jóvenes (YSO) revelan variaciones de luminosidad episódicas, conocidas como erupciones de tipo FU Orionis. Un mecanismo teórico prominente para explicar esto es la inestabilidad térmica (TI) en el borde interior de la "zona muerta" (DZIE) de los discos protoplanetarios. En esta región, la transición entre una zona interna ionizada y activa (donde opera la inestabilidad magnetorrotacional, MRI) y una zona externa fría y menos turbulenta (zona muerta) provoca la acumulación de masa. Cuando el calentamiento viscoso supera el enfriamiento radiativo, se desencadena una erupción.

El problema central abordado en este trabajo es que los modelos anteriores han tratado el polvo de manera simplificada (distribución de tamaños fija y acoplada al gas). Sin embargo, el polvo es el principal portador de opacidad en la zona muerta y su dinámica (crecimiento, fragmentación, arrastre) tiene efectos térmicos y dinámicos críticos que no se han explorado completamente. El objetivo es investigar cómo la evolución dinámica del polvo influye en las erupciones de acreción y en la formación de subestructuras ricas en polvo que podrían dar lugar a planetesimales.

Metodología:

Los autores realizaron simulaciones de hidrodinámica radiativa utilizando el código PLUTO en un marco axisimétrico 1D (integración vertical).

1. Física del Gas y Polvo:

   • Se resolvieron las ecuaciones de Navier-Stokes para el gas y un componente de polvo sin presión.
   • Se incluyó calentamiento viscoso, irradiación estelar, enfriamiento radiativo (superficial y difusión en el plano) y dinámica de polvo-gas.
   • Modelo de Evolución de Polvo: Se implementó el método TriPoD (Pfeil et al. 2024), que rastrea la evolución de poblaciones de granos "pequeños" y "grandes", así como el tamaño máximo ($a_{max}$), permitiendo coagulación y fragmentación.
   • Opacidad: Se utilizó una opacidad dependiente del tamaño del grano, la distribución de tamaños ($q_{dust}$) y la temperatura, calculada con el módulo growpacity. Se modeló la sublimación del polvo a altas temperaturas ($T_{subl} \approx 2000$ K).

2. Configuración Numérica:

   • Se modeló un disco alrededor de una estrella de $1 M_{\odot}$.
   • Se varió el coeficiente de viscosidad en la zona muerta ($\alpha_{DZ}$) entre $10^{-3}$ (modelo fiducial) y $10^{-4}$ para estudiar el papel de la difusión viscosa.
   • Se compararon tres escenarios: modelo fiducial (con evolución de polvo), modelo sin evolución de polvo (distribución fija) y modelo con difusión radiativa en el plano.

Contribuciones Clave:

1. Modelado Autoconsistente del Polvo: Es el primer estudio que integra la evolución dinámica completa del polvo (coagulación/fragmentación) y su acoplamiento térmico-dinámico con el gas en simulaciones de erupciones de acreción en la zona muerta.
2. Método Semi-Analítico: Se desarrolló un método para predecir los estados pre- y post-erupción basándose en el equilibrio térmico y la viscosidad, validado contra las simulaciones numéricas completas.
3. Análisis de la Formación de Planetesimales: Se evaluó si las estructuras de polvo resultantes de las erupciones cumplen los criterios de inestabilidad de flujo (streaming instability, SI) para formar planetesimales.

Resultados Principales:

• Impacto de la Evolución del Polvo en las Erupciones:

  • La fragmentación del polvo durante la fase de calentamiento aumenta drásticamente la opacidad (al haber más granos pequeños), lo que intensifica la erupción y permite que el frente de calentamiento penetre más profundamente en la zona muerta (hasta $\sim 0.9$ au en el modelo fiducial, frente a $\sim 0.6$ au sin evolución de polvo).
  • La inclusión de la evolución dinámica del polvo aumenta la masa acrecida por un factor de $\sim 1.5$ en comparación con modelos estáticos.
  • La difusión radiativa en el plano tiene un efecto secundario menor en comparación con la evolución del polvo.

• Formación de Estructuras Anulares:

  • Las erupciones generan múltiples anillos de gas y polvo que se extienden hasta $\sim 1$ au dentro de la zona muerta.
  • Estos anillos contienen masas de polvo de hasta $\sim 1.6 M_{\oplus}$ (masas terrestres).
  • Durante la fase de erupción, el polvo se sublima y fragmenta; tras el paso del frente de calor, el polvo se recondensa y crece rápidamente por coagulación, alcanzando tamaños de milímetros en escalas de tiempo comparables a la duración de la erupción.

• Evolución a Largo Plazo y Formación de Planetesimales:

  • Tras la erupción, el disco entra en una fase cuiescente donde los anillos de presión formados actúan como trampas para el polvo.
  • En modelos con menor viscosidad ($\alpha_{DZ} = 10^{-4}$), la difusión viscosa es más lenta, permitiendo que los granos crezcan más (hasta $\sim 3$ mm) y se acumulen con mayor eficiencia en los máximos de presión.
  • Simulaciones de Formación: Utilizando criterios de inestabilidad de flujo (SI) recientes (Lim et al. 2026), se encontró que solo en el caso más favorable ($\alpha_{DZ} = 10^{-4}$) se supera el umbral crítico de concentración de polvo. En este escenario, tras 50 kyr de evolución, se forman aproximadamente 1.6 $M_{\oplus}$ de planetesimales, distribuidos entre los picos de acumulación post-erupción y el borde interior.

• Parámetros Temporales:

  • Las erupciones duran $\sim 85$ años, seguidas de fases cuiescentes de $\sim 3200$ años (para $\alpha_{DZ} = 10^{-3}$) hasta $\sim 100$ kyr (para $\alpha_{DZ} = 10^{-4}$).

Significado e Implicaciones:

1. Importancia Térmica y Dinámica del Polvo: El estudio demuestra que ignorar la evolución del polvo lleva a subestimar la intensidad y extensión de las erupciones de acreción. El polvo no es un espectador pasivo; su fragmentación y crecimiento modulan la opacidad y, por tanto, la termodinámica del disco.
2. Vía de Formación de Planetas: Se propone un mecanismo viable para la formación de planetas rocosos (super-Tierras y mini-Neptunos) cerca de la estrella ($\lesssim 1$ au). Las erupciones periódicas crean trampas de presión eficientes donde el polvo se acumula y crece, facilitando la inestabilidad de flujo y la formación de planetesimales in situ.
3. Limitaciones y Futuro: Aunque los resultados son prometedores, el modelo 1D no puede capturar inestabilidades no axisimétricas como la inestabilidad de ondas de Rossby (RWI), que podrían erosionar los frentes de erupción o crear vórtices que afecten la acumulación de polvo. Además, los criterios de formación de planetesimales aún tienen incertidumbres significativas.

En conclusión, este trabajo establece que el tratamiento autoconsistente de la dinámica del polvo es crucial para entender la evolución de los discos protoplanetarios y sugiere que las erupciones de acreción en el borde de la zona muerta podrían ser un motor fundamental para la formación temprana de planetas en el sistema solar interior.