exoALMA XX: Tomographic Detection of Embedded Planets in Protoplanetary Disks

1. Planteamiento del Problema

Las observaciones de alta resolución angular con ALMA han revelado una gran variedad de subestructuras en discos protoplanetarios (huecos, anillos, espirales). Si bien muchas de estas estructuras sugieren la formación de planetas, otras pueden ser causadas por inestabilidades del disco (como la inestabilidad de cizalladura vertical - VSI, la inestabilidad gravitacional - GI, o la inestabilidad magnetorrotacional - MRI).

El desafío principal reside en distinguir inequívocamente entre las señales cinemáticas provocadas por planetas incrustados y aquellas generadas por inestabilidades fluidas. Los métodos tradicionales basados únicamente en la cinemática del centroide de la línea (desviaciones del movimiento kepleriano) a menudo son ambiguos o insuficientes para localizar planetas con precisión, especialmente en discos con baja inclinación o cuando las perturbaciones son sutiles.

2. Metodología

El estudio emplea un enfoque tomográfico que va más allá del análisis de centroides, utilizando simulaciones hidrodinámicas y observaciones sintéticas para desarrollar y validar nuevos diagnósticos:

• Simulaciones y Observaciones Sintéticas:
  • Se utilizaron simulaciones de interacción planeta-disco (código fargo3d) con planetas de masas $0.1\%$ y $0.2\%$ de la masa estelar.
  • Se incluyeron simulaciones de inestabilidades del disco (VSI, MRI, GI) usando el código PLUTO.
  • Se generaron observaciones sintéticas de líneas de CO ($^{12}$CO y $^{13}$CO $J=3\to2$) replicando las condiciones del programa exoALMA (resoluciones de $0.15''$ y $0.30''$, espaciado de canal de $0.1$ km/s).
• Análisis Tomográfico y "Folded Maps":
  • Se aplicó un procedimiento de plegado (folding) a los mapas de velocidad y anchura de línea. Esto implica sumar o restar los lados corrimiento al rojo y al azul del disco para suprimir flujos simétricos (como la rotación kepleriana global) y aislar perturbaciones no simétricas.
  • Anchura de Línea (Line Width): Se utilizó el plegado sustractivo para buscar aumentos localizados en la dispersión de velocidad (ensanchamiento de línea), que son menos sensibles a la posición azimutal del planeta que las señales de velocidad.
  • Asimetría de Línea: Se analizó la forma del perfil de la línea (sesgo o skewness) para distinguir la naturaleza de los flujos de velocidad.
  • Algoritmos de Detección: Se emplearon algoritmos de agrupamiento (K-means) y estimadores estadísticos (desviación absoluta mediana - MAD) para identificar picos significativos ($>3\sigma$) en los mapas residuales.
• Aplicación a Datos Reales:
  • La metodología se aplicó a los datos de exoALMA de los discos HD 135344B y MWC 758, comparando los residuos con modelos keplerianos suaves (discminer).

3. Contribuciones Clave

1. Nueva Firma de Detección: Demostración de que el ensanchamiento de línea localizado es un indicador más robusto y directo de la presencia de planetas que las perturbaciones de velocidad del centroide, especialmente para planetas de masa intermedia ($>0.1\% M_*$).
2. Diagnóstico de Asimetría: Identificación de que las líneas ópticamente gruesas cerca de un planeta muestran una asimetría pronunciada (cola hacia velocidades bajas y corrimiento al azul), causada por flujos meridionales convergentes. En contraste, las inestabilidades del disco producen un ensanchamiento más isotrópico y simétrico.
3. Robustez Observacional: Se establece que la detección de planetas es más dependiente de la relación señal-ruido (SNR) que de la resolución angular o espectral. Degradar la resolución angular para aumentar el SNR puede mejorar la detección de planetas de baja masa.
4. Distinción Mecanismo-Origen: Desarrollo de un criterio claro para diferenciar entre señales planetarias y de inestabilidad del disco basándose en la morfología del perfil de la línea y la coherencia de los flujos de velocidad.

4. Resultados Principales

A. Resultados en Simulaciones

• Límites de Detección: Es posible detectar planetas con masas superiores al 0.1% de la masa estelar (aprox. 1 $M_{Jup}$ para una estrella solar) con solo unas horas de integración de ALMA.
• Precisión: El ensanchamiento de línea permite localizar planetas con una precisión de menos de media anchura de haz en radio y pocos grados en azimut.
• Influencia de la Inclinación: La detección es robusta para inclinaciones de $-15^\circ$ a $-30^\circ$. A mayor inclinación, las señales de velocidad se vuelven menos localizadas, pero el ensanchamiento de línea sigue siendo detectable.
• Distinción: Las inestabilidades (VSI, MRI, GI) producen perturbaciones localizadas, pero sus perfiles de línea carecen de la fuerte asimetría y los flujos meridionales coherentes observados en las interacciones planeta-disco.

B. Aplicación a HD 135344B

Se identificaron tres candidatos a planetas incrustados basados en perturbaciones localizadas de velocidad y anchura de línea:

1. Candidato C1: En $R = 41$ au, $\phi = 36^\circ$. Asociado a ensanchamiento de línea en capas superiores, posiblemente un planeta de baja masa ($<1.6 M_{Jup}$).
2. Candidato C2: En $R = 73$ au, $\phi = 15^\circ$. Presenta ensanchamiento de línea significativo en $^{12}$CO y $^{13}$CO, así como asimetrías en el perfil de línea consistentes con flujos meridionales.
3. Candidato DF (Doppler Flip): En $R = 95$ au, $\phi = -133^\circ$. Detectado mediante una inversión de signo de velocidad (Doppler flip) y un mínimo en la anchura de línea promediada azimutalmente. Este candidato explica grandes estructuras espirales y coincide con propuestas anteriores basadas en polvo.

C. Aplicación a MWC 758

• El disco está dominado por grandes estructuras espirales en velocidad y anchura de línea que se extienden desde el interior hasta $\sim200$ au.
• A diferencia de HD 135344B, no se encontraron perturbaciones localizadas (picos aislados) que sugieran un planeta masivo incrustado en las regiones externas.
• Las espirales de velocidad están dominadas por movimientos verticales (hacia arriba y hacia abajo).
• Interpretación: La morfología observada es más consistente con un disco excéntrico inducido por un compañero subestelar en las regiones internas, o con un disco con geometría torcida (warped), en lugar de un planeta masivo externo excitando ondas de densidad.

5. Significado e Impacto

Este trabajo representa un avance significativo en la astrofísica de discos protoplanetarios al proporcionar una metodología sistemática para "ver" planetas ocultos en el gas, más allá de las estructuras de polvo.

• Validación de la Formación Planetaria: Confirma que las interacciones planeta-disco dejan huellas cinemáticas específicas (ensanchamiento asimétrico) que pueden ser detectadas incluso cuando el planeta no es visible directamente.
• Optimización de Observaciones: Sugiere que para la búsqueda de exoplanetas en discos, es más eficiente priorizar la sensibilidad (SNR) sobre la resolución angular máxima, permitiendo detectar planetas de masa joviana con observaciones más rápidas y menos costosas.
• Resolución de Ambigüedades: Ofrece una herramienta crucial para distinguir entre la formación de planetas y las inestabilidades gravitacionales o magnéticas, un debate central en la comprensión de la evolución de los discos.
• Nuevos Objetivos: Proporciona candidatos concretos para futuras observaciones de alta resolución y seguimiento directo, especialmente en sistemas como HD 135344B donde la evidencia cinemática es fuerte.

En resumen, el estudio establece que la tomografía de líneas moleculares, combinando análisis de anchura, asimetría y flujos verticales, es una técnica poderosa y necesaria para mapear la población de planetas incrustados en los discos de formación estelar.