DiskMINT-GARDEN: Self-consistent Models to Estimate Disk Masses

Resumen Técnico: DiskMINT-GARDEN: Modelos autoconsistentes para estimar las masas de los discos

Planteamiento del problema
Determinar la masa de gas de los discos protoplanetarios es fundamental para comprender la formación planetaria, pero sigue siendo un desafío. El monóxido de carbono (CO) es el trazador principal del gas en los discos, pero los modelos termoquímicos tempranos encontraban frecuentemente que la emisión de los isotopólogos del CO observados era significativamente más débil de lo predicho. Esta discrepancia se interpretaba a menudo como evidencia de una sustracción sustancial de CO en fase gaseosa (hasta dos órdenes de magnitud) debido al secuestro en planetesimales o al procesamiento químico. Sin embargo, tales factores de sustracción introducen una incertidumbre significativa en las estimaciones de masa ($M_{\text{gas}}$). Además, los modelos existentes suelen depender de supuestos simplificados, como estructuras de densidad vertical isotérmicas o factores de sustracción fijos, que pueden no contabilizar de manera autoconsistente el acoplamiento entre la estructura, la temperatura y la química del disco. Aunque existen restricciones independientes vía mediciones dinámicas (limitadas a discos masivos) o observaciones de deuteruro de hidrógeno (HD) (limitadas por la sensibilidad y la disponibilidad de instalaciones), se necesita un método rápido, reproducible y autoconsistente para estimar las masas de los discos a partir de datos de ALMA comúnmente disponibles.

Metodología
Los autores presentan DiskMINT-GARDEN (Grid of Astrochemical Radiative Disk EmissioN), un marco que combina una cuadrícula exhaustiva de modelos autoconsistentes con una herramienta de inferencia de aprendizaje automático.

• Construcción de la cuadrícula del modelo: La cuadrícula se construye sobre el marco DiskMINT, que acopla la estructura hidrostática del disco, la transferencia radiativa de continuo/líneas y una red química reducida. Los componentes físicos clave incluyen:

  • Estructura vertical autoconsistente: El modelo resuelve el equilibrio hidrostático vertical, vinculando la masa del plano medio con las capas emisoras, en lugar de asumir un perfil de densidad isotérmico fijo.
  • Química: La red incluye la congelación de CO, la conversión en la superficie de los granos (específicamente de CO a CO$_2$) y la fotodisociación selectiva de isótopos. Notablemente, no asume factores de sustracción elemental ad hoc.
  • Espacio de parámetros: La cuadrícula abarca masas estelares ($0.1 - 2.0 M_\odot$), masas de gas del disco ($10^{-5} - 10^{-1} M_\star$), relaciones polvo-gas ($0.003 - 0.1$) y radios característicos ($10 - 300$ au).
  • Observables sintéticos: Para cada modelo, el marco computa la luminosidad del continuo de milimétrico ($L_{\text{mm}}$ a $\sim 1.3$ mm), las luminosidades de las líneas de C$^{18}$O (2–1) y (3–2), y métricas de tamaño de polvo ($R_{\text{dust,90}}$).

• Herramienta de inferencia: Se entrena un modelo de regresión de aprendizaje automático supervisado basado en árboles de decisión de gradiente aumentado (XGBoost) sobre la cuadrícula del modelo. Esta herramienta mapea el vector de observables $\mathcal{O} = (L_{\text{C}^{18}\text{O}}, L_{\text{mm}}, R_{90,\text{dust}})$ al vector de parámetros físicos $\Theta = (M_{\text{gas}}, \epsilon, R_c)$ para una masa estelar dada. La herramienta está diseñada para predicciones rápidas y deterministas dentro del dominio de la cuadrícula.

• Validación e incertidumbre: Las incertidumbres estadísticas se estiman propagando los errores observacionales y la dispersión de la interpolación (típicamente $\lesssim 25\%$). Las incertidumbres sistemáticas derivadas de los supuestos de los modelos subyacentes (estructura, química) se estiman conservadoramente como un factor de $\sim 2$. El marco incluye controles de dominio para señalar objetivos que se encuentran fuera de los límites de la cuadrícula de entrenamiento.

Resultados clave
Los autores aplicaron DiskMINT-GARDEN a una muestra de archivo de 34 discos protoplanetarios con detecciones de continuo y C$^{18}$O de ALMA.

1. Acuerdo con restricciones independientes: Las masas de gas inferidas muestran un buen acuerdo (dentro de un factor de $\lesssim 2$) con las estimaciones de masa dinámica independientes para discos masivos y con las estimaciones de masa basadas en HD para los pocos discos donde se detecta HD. Esto valida la capacidad del marco para reproducir masas en regímenes donde las restricciones dinámicas son robustas.
2. Comparación con modelos DALI: Al compararse con las masas derivadas utilizando el marco termoquímico DALI (que típicamente requiere grandes factores de sustracción de CO para coincidir con las observaciones de N$_2$H$^+$ y CO), DiskMINT-GARDEN produce masas de gas totales similares sin invocar ningún factor de sustracción.
   • En DALI, la falta de química en la superficie de los granos y una estructura de densidad isotérmica fija requieren un factor de sustracción (frecuentemente $\sim 10\times$) para coincidir con las observaciones.
   • En DiskMINT, la inclusión de la química en la superficie de los granos (convirtiendo el hielo de CO en CO$_2$) y una estructura vertical autoconsistente (situando la línea de nieve de CO a temperaturas más altas, $\sim 35$ K) reproduce naturalmente la emisión de C$^{18}$O observada y las restricciones de N$_2$H$^+$ sin necesidad de sustracciones ad hoc.
3. Implicaciones para la evolución de los discos: El hecho de que DiskMINT-GARDEN recupere masas consistentes con DALI sin sustracción sugiere que los datos actuales no requieren un procesamiento químico significativo (como el secuestro de CO en planetesimales) para explicar la emisión observada. El CO "faltante" en interpretaciones previas puede explicarse mediante la estructura física y la química de superficie de granos implementada en DiskMINT.

Significación y afirmaciones
El artículo afirma que DiskMINT-GARDEN proporciona un método rápido, de código abierto y homogéneo para estimar las masas de gas de los discos a partir de datos de ALMA ampliamente disponibles (continuo y C$^{18}$O), evitando la necesidad de observaciones de tiempo intensivo de trazadores adicionales como N$_2$H$^+$.

• Eficiencia: Permite estudios a nivel de población de la evolución de los discos, extendiendo las estimaciones de masa a discos de menor masa donde las restricciones dinámicas son inviables.
• Consistencia física: Al demostrar que modelos autoconsistentes con química de superficie de granos pueden reproducir las observaciones sin grandes factores de sustracción, este trabajo desafía la necesidad de asumir una sustracción significativa de CO en los modelos evolutivos de los discos.
• Escalabilidad: El marco sirve como una herramienta de inferencia de primer orden que puede guiar el modelado detallado y espacialmente resuelto de objetivos individuales. Los autores señalan que, aunque la cuadrícula actual asume propiedades de polvo fijas, el marco puede extenderse para incorporar trazadores adicionales (por ejemplo, HD de futuras instalaciones) y datos de continuo de múltiples bandas para restringir aún más los parámetros del disco.

Los autores concluyen que los datos existentes son consistentes con poco procesamiento químico debido a procesos evolutivos de los discos, y que las discrepancias en las estimaciones de masa anteriores a menudo surgen de las diferencias en cómo se acoplan la estructura y la química del disco, en lugar de una sustracción elemental real.