Can photoevaporation open gaps in protoplanetary discs?

Planteamiento del problema
La dispersión de los discos protoplanetarios es un factor crítico para determinar el cronograma y las condiciones de la formación planetaria. Si bien la fotoevaporación —la eliminación de gas mediante vientos impulsados térmicamente lanzados por la radiación estelar de alta energía— es ampliamente reconocida como un mecanismo de dispersión clave, la base teórica de cómo abre y mantiene brechas sigue siendo incierta. Los modelos estándar suelen depender de cálculos de evolución viscosa unidimensionales (1D) donde la tasa de pérdida de masa por fotoevaporación, $\dot{\Sigma}(r)$, se trata como un término de sumidero fijo derivado de simulaciones bidimensionales (2D) estáticas. Este enfoque asume implícitamente que el perfil de pérdida de masa permanece constante incluso cuando la estructura del disco evoluciona. Los autores argumentan que esto ignora un bucle de retroalimentación crucial: a medida que comienza a formarse una brecha, la estructura local del disco cambia, alterando potencialmente las condiciones de lanzamiento del viento, reduciendo las tasas locales de pérdida de masa y permitiendo que el flujo viscoso rellene la región deprimida. Si la fotoevaporación es autolimitada de esta manera, puede ser menos eficiente en la limpieza del gas en los radios donde se forman los planetas gigantes de lo que se suponía anteriormente, con implicaciones significativas para los modelos de acreción de núcleo y migración.

Metodología
Para investigar si la fotoevaporación por sí sola puede abrir y mantener brechas cuando el flujo responde dinámicamente a la estructura del disco, los autores emplearon un enfoque computacional de dos fases:

1. Fase I (Evolución 1D): Utilizaron el código de evolución viscosa spock para simular un disco sometido a fotoevaporación por rayos X hasta el inicio de la formación de la brecha. La simulación evolucionó un perfil de densidad superficial de gas ($\Sigma$) durante 2.1 Myr, incorporando un perfil de pérdida de masa de Picogna et al. (2019). Esta fase estableció las condiciones iniciales para las simulaciones 2D, específicamente un estado donde se había desarrollado una depresión parcial en la densidad superficial.
2. Fase II (Radiación-Hidrodinámica 2D): El perfil de densidad superficial de la simulación 1D a los 2.1 Myr se mapeó en un sistema de coordenadas esféricas axisimétrico 2D ($r-\theta$) utilizando el código pluto. Estas simulaciones modelaron el flujo de salida fotoevaporativo y su interacción con el disco.
   • Modelo p-g: Utilizó la densidad superficial evolucionada de la Fase I, que incluía una depresión de densidad preexistente.
   • Modelo p-ng: Utilizó un perfil de densidad superficial inicial modificado donde la depresión fue suprimida (reemplazada por una ley de potencia) para evaluar el impacto de la morfología inicial en la apertura de la brecha.
   • Las simulaciones emplearon una aproximación de calentamiento en la línea de visión y asumieron equilibrio radiativo, desacoplando la transferencia radiativa de la hidrodinámica para reducir el costo computacional.

Resultos clave

• Formación de brechas autolimitadas: En las simulaciones 2D, una vez que se formó una depresión de densidad, la tasa local de pérdida de masa disminuyó drásticamente dentro de la región de la brecha. Esto ocurrió porque la región deprimida permaneció en la sombra del disco interno, reduciendo la exposición a la irradiación estelar. En consecuencia, la brecha no se profundizó hasta convertirse en una cavidad totalmente evacuada.
• Relleno parcial: El flujo viscoso y el transporte de masa radial a lo largo de la superficie del disco (y a través del viento) actuaron para rellenar parcialmente la región deprimida. El sistema evolucionó hacia una zona de depresión parcial persistente en lugar de un hueco interno limpio.
• Insensibilidad a las condiciones iniciales: La comparación de los modelos con (p-g) y sin (p-ng) una depresión inicial mostró que el resultado evolutivo a largo plazo fue ampliamente independiente de la morfología inicial del disco. En ambos casos, el sistema se autorreguló hacia un estado de depresión parcial.
• Limitaciones de las prescripciones 1D: Los autores probaron una prescripción de primer orden para aproximar el comportamiento 2D en modelos 1D escalando el perfil de pérdida de masa primordial basado en la depresión de la densidad superficial local. Aunque esto mejoró la concordancia respecto a las prescripciones estáticas, los modelos 1D todavía evolucionaron demasiado rápido y produjeron brechas más profundas que las simulaciones 2D. Esta discrepancia resalta la incapacidad de los modelos 1D para capturar el transporte de masa radial a lo largo de la superficie del disco y la compleja interacción de sombreado.
• Firmas observacionales: Aunque la brecha de gas permanece parcialmente llena, el máximo de presión en el borde exterior de la depresión aún puede actuar como una trampa eficiente para granos de polvo. Esto podría producir firmas observacionales que asemejen a los discos de transición (cavidades internas en la emisión de continuo) incluso sin una limpieza completa del gas.

Significancia y Afirmaciones
El artículo desafía el paradigma estándar de que la fotoevaporación talla eficientemente cavidades internas limpias y produce directamente discos de transición. En su lugar, los autores proponen que la fotoevaporación puede conducir a un estado cuasi-estacionario caracterizado por una zona de depresión parcial persistente. Esto sugiere que los discos que exhiben una depresión moderada sin una limpieza interna completa pueden representar sistemas influenciados por la fotoevaporación pero no totalmente dispersados.

Los autores proporcionan una prescripción modesta de primer orden para aproximar este comportamiento autorregulado en modelos de evolución de discos 1D, la cual es adecuada para su uso en estudios de formación planetaria y síntesis de poblaciones. Sin embargo, declaran explícamente que esta prescripción sigue siendo aproximada y no puede capturar plenamente la naturaleza multidimensional del flujo, particularmente la redistribución de masa radial. El artículo concluye que son necesarias más simulaciones multidimensionales y la exploración del espacio de parámetros para derivar recetas robustas para modelos globales de discos y poblaciones de planetas. El trabajo subraya la necesidad de ir más allá de los tratamientos de pérdida de masa estáticos para comprender la interacción entre los vientos internos y las vías de formación planetaria.