The dispersal of compact protoplanetary discs

Resumen Técnico: La Dispersión de Discos Protoplanetarios Compactos

Planteamiento del Problema
Los discos protoplanetarios son los sitios de formación planetaria, pero sus mecanismos de dispersión siguen siendo objeto de investigación activa. La evidencia observacional, particularmente de estudios de evolución color-color de poblaciones estelares jóvenes, indica que la dispersión de los discos procede preferentemente de adentro hacia afuera, creando una cavidad central antes de que el disco exterior sea erosionado. Sin embargo, los modelos teóricos estándar suelen asumir discos con grandes radios exteriores (cientos de unidades astronómicas), muy más allá del radio gravitacional ($R_g$) donde la fotoevaporación es más eficiente. Recientes sondeos (por ejemplo, ALMA) revelan que una fracción sustancial de discos son "compactos", con radios exteriores comparables o menores que $R_g$. Los modelos existentes, que típicamente asumen perfiles de pérdida de masa no restringidos que se extienden a grandes radios, no logran reproducir con precisión las trayectorias evolutivas y las escalas de tiempo de dispersión de estos sistemas compactos. Específicamente, las prescripciones estándar a menudo predicen una erosión de afuera hacia adentro para discos compactos, contradiciendo los diagnósticos observacionales.

Metodología
Los autores emplearon un enfoque de múltiples etapas que combina simulaciones hidrodinámicas de radiación de alta resolución con modelos de síntesis de poblaciones viscosos 1D:

1. Simulaciones Hidrodinámicas de Radiación: Utilizando el código PLUTO acoplado al marco de transferencia radiativa MOCASSIN, los autores simularon la fotoevaporación impulsada por rayos X para una estrella de $1 M_\odot$. Exploraron cuatro radios de corte iniciales del disco ($R_{cut}$) distintos: 10, 50, 100 y 200 ua. Las simulaciones rastrearon la estructura térmica, la densidad y los campos de velocidad hasta alcanzar un estado cuasi-estacionario. El radio gravitacional se adoptó como $R_g = 5$ ua.
2. Derivación de Relaciones de Escala: A partir de las simulaciones hidrodinámicas, los autores analizaron las distribuciones acumuladas de la tasa de pérdida de masa. Determinaron cómo depende la tasa total de pérdida de masa del viento del radio de corte exterior del disco.
3. Síntesis de Poblaciones: Estos hallazgos se incorporaron en un código de evolución viscosa 1D (SPOCK). Los autores realizaron una síntesis de poblaciones de 4.096 discos sintéticos utilizando una estrategia de muestreo Cuasi-Monte Carlo (secuencias de Sobol) para explorar un espacio de parámetros de 9 dimensiones (masa estelar, luminosidad en rayos X, masa inicial del disco, radio de corte, viscosidad, campo FUV externo, etc.).
4. Prescripciones Probadas: Los modelos compararon dos escenarios:
   • Prescripción Estándar: Tasas de pérdida de masa por fotoevaporación interna aplicadas sin tener en cuenta el radio exterior finito del disco (asumiendo efectivamente que el disco se extiende hasta el infinito o más allá del pico de pérdida de masa).
   • Prescripción Dependiente del Corte: La tasa de pérdida de masa se trunca y normaliza en función del radio de corte específico del disco, derivada de las relaciones de escala hidrodinámicas.
   • Fotoevaporación Externa: Algunos modelos incluyeron fotoevaporación externa FUV leve (utilizando la cuadrícula FRIED) para simular estrellas OB cercanas.

Contribuciones y Resultados Clave

• Escala Hidrodinámica: Las simulaciones revelaron que los perfiles de pérdida de masa superficial son casi independientes del tamaño del disco; la morfología del viento y las tasas locales de pérdida de masa en las regiones internas permanecen consistentes independientemente del corte exterior. Sin embargo, la tasa total de pérdida de masa integrada se reduce significativamente para discos compactos porque el viento se trunca en el radio de corte. La tasa de pérdida de masa acumulada escala con el radio, creciendo como $\dot{M}_w \propto R^3$ cerca del radio gravitacional y aplanándose a $\dot{M}_w \propto R^{0.5}$ a radios mayores.
• Trayectorias Evolutivas:
  • Prescripción Estándar: Cuando se aplica a discos compactos ($R_{cut} < 20-30$ ua), el modelo estándar predice que la tasa de pérdida de masa del viento excede la tasa de expansión viscosa fuera del radio de corte. Esto conduce a una dispersión de afuera hacia adentro, donde los bordes exteriores se erosionan primero, contradiciendo la evidencia observacional.
  • Prescripción Dependiente del Corte: Cuando la tasa de pérdida de masa se escala al radio de corte, la tasa total del viento se reduce lo suficiente como para permitir que la expansión viscosa domine inicialmente. Esto resulta en una dispersión de adentro hacia afuera, donde se abre una brecha cerca del radio gravitacional y se expande hacia el exterior, consistente con los diagnósticos observacionales.
• Resultados de la Síntesis de Poblaciones:
  • La prescripción dependiente del corte, combinada con una tasa de fotoevaporación interna reducida (escalada por un factor de 10 para alinearse con modelos químicos recientes), reproduce las fracciones de discos observadas y las tasas de acreción en función de la edad del cúmulo sin ajuste de parámetros.
  • La vida media del disco en el modelo mejorado es de $2.70 \pm 0.13$ Myr, lo que se alinea bien con las estimaciones observacionales ($2.24 \pm 0.35$ Myr) y es significativamente más larga que las predicciones de modelos sin la reducción de la tasa ($1.13 \pm 0.05$ Myr).
  • La inclusión de fotoevaporación externa leve no impide la expansión del disco cuando se utiliza la prescripción del radio de corte, pero resulta en un aumento más suave de los radios del disco con el tiempo y en una mayor fracción de discos compactos que persisten hasta edades más avanzadas.

Significado
El artículo argumenta que la dispersión de discos protoplanetarios compactos no puede capturarse con precisión mediante modelos de fotoevaporación estándar no escalados. Los autores proponen una relación de escala simple y motivada físicamente que tiene en cuenta el radio exterior finito del disco. Al incorporar esta dependencia, los modelos reconcilian con éxito las predicciones teóricas con la limpieza observada de adentro hacia afuera de las poblaciones estelares jóvenes y la prevalencia de discos compactos. El estudio concluye que las prescripciones de dispersión deben incluir explícitamente la dependencia del radio de corte del disco para predecir correctamente la evolución y la vida útil de los discos compactos, proporcionando una corrección necesaria a los modelos actuales de síntesis de poblaciones.