Southern massive stars at high angular resolution. Physical separations and mass ratios

Resumen Técnico: Estrellas Masivas del Sur a Alta Resolución Angular: Separaciones Físicas y Razones de Masa

Problema y Contexto
La evolución de las estrellas masivas está regulada fundamentalmente por su multiplicidad, ya que las interacciones como la transferencia de masa y las fusiones alteran significamente sus trayectorias evolutivas y los productos de su final de vida. Mientras que los sondeos espectroscópicos han caracterizado binarias estrechas ($<1$ UA) y las técnicas de imagen han identificado sistemas amplios ($>10^3$ UA), el régimen de separación intermedia (aproximadamente de 1 milisegundo de arco a 8 segundos de arco) ha permanecido en gran medida inexplorado. El sondeo Southern Massive Stars at High Angular Resolution (smash+), utilizando interferometría (PIONIER) y enmascaramiento de apertura dispersa/óptica adaptativa (NACO), detectó previamente una alta fracción de compañeras en este dominio. Sin embargo, los resultados iniciales del sondeo se presentaron en unidades observacionales (separación angular y contraste de magnitud). Para comprender la naturaleza física de estos sistemas y sus mecanismos de formación, estos parámetros observacionales deben convertirse en unidades físicas: separaciones físicas proyectadas (en UA) y razones de masa ($q = M_2/M_1$).

Metodología
Este artículo detalla la conversión de los datos observacionales de smash+ en parámetros físicos mediante un proceso de calibración y modelado de múltiples pasos:

1. Determinación de la Distancia: Reconociendo que las paralaxas de Gaia pueden ser poco fiables para estrellas brillantes, calientes o binarias debido a desplazamientos del fotocentro y sistemáticas, los autores desarrollaron una calibración de distancia espectrofotométrica. Este método utiliza calibraciones de magnitud absoluta basadas en el tipo espectral (SpT) y la clase de luminosidad (LC) derivadas de Martins & Plez (2006).

   • Las magnitudes para clases de luminosidad intermedias (IV, II) se derivaron mediante interpolación lineal.
   • Las magnitudes aparentes se corrigieron por la contribución de flujo de las compañeras resueltas dentro de la función de dispersión de punto de 2MASS ($4''$) y las contribuciones estimadas de binarias internas no resueltas (donde se conocían los tipos espectrales).
   • Se aplicaron correcciones de extinción utilizando curvas de Fitzpatrick (1999).
   • Las distancias resultantes se validaron contra valores de literatura de cúmulos y el catálogo basado en Gaia Alma Luminous Star III (ALS III), revelando una sobreestimación sistemática de $\sim5\%$ en las distancias espectrofotométricas, la cual se consideró insignificante para el análisis logarítmico de las separaciones.

2. Determinación de la Masa:

   • Primarias: Las masas se derivaron utilizando relaciones masa-luminosidad calibradas a partir de modelos de Martins et al. (2005) y Brott et al. (2011), ajustadas para clases de luminosidad específicas.
   • Compañeras: Dado que las clases de luminosidad de las compañeras son generalmente desconocidas, los autores asumieron una clase de luminosidad de enana (LC V) para la mayoría de las compañeras, justificado por la expectativa de que las compañeras de menor masa evolucionan más lentamente. Para las compañeras más brillantes donde una suposición de enana implicaría un $q > 1$ (físicamente imposible si la primaria está más evolucionada), los autores adoptaron la clase de luminosidad menos evolucionada que resultara en un $q < 1$ para establecer un límite superior de masa.

3. Sensibilidad y Límites de Detección: Los autores emplearon simulaciones de Monte Carlo para mapear las probabilidades de detección en el espacio de la razón de masa ($q$) frente al semieje mayor ($a$). Estas simulaciones tuvieron en cuenta los límites instrumentales específicos del sondeo (ángulos de trabajo interno/externo y límites de contraste) y las orientaciones orbitales aleatorias (excentricidad, inclinación, argumento del periastro).

Contribuciones Clave y Resultados

• Separaciones Físicas: Las separaciones físicas proyectadas de las compañeras siguen una distribución uniforme en la separación logarítmica (una ley de Öpik), consistente con una distribución plana en $\log(a)$. Los autores señalan un potencial exceso de compañeras alrededor de 80–100 AU, aunque esto puede ser un artefacto de la transición entre modos instrumentales.
• Distribuciones de la Razón de Masa: La distribución de la razón de masa sigue una ley de potencia, $f_q \propto q^\kappa$. El exponente $\kappa$ varía con la separación:
  • Para la muestra completa, $\kappa \approx -0.7$.
  • Para separaciones $<100$ AU, $\kappa_{<100} = -0.6^{+0.9}_{-0.7}$. Este valor es más plano que el $\kappa = -1.4 \pm 0.4$ reportado por Moe & Di Stefano (2017) para una submuestra más pequeña, pero es compatible con las distribuciones planas ($\kappa \approx 0$) encontradas en binarias espectroscópicas.
  • Para separaciones $>1000$ AU, $\kappa_{>1000} = -1.6^{+1.0}_{-1.2}$, indicando una distribución más pronunciada que favorece razones de masa más bajas.
• Correlación Separación vs. Razón de Masa: Se observa una tendencia distinta donde la razón de masa máxima detectable declina a medida que aumenta la separación. Más allá de $\approx 100$ AU, hay una clara falta de compañeras de masa casi igual ($q > 0.8$). En separaciones $>1000$ AU, no se detectan compañeras con $q > 0.7$, y por encima de $10^4$ AU, el límite cae a $q \approx 0.5$. Los autores derivan un envolvente superior empírico para esta "zona de exclusión" (Ec. 6).
• Completitud: El sondeo se considera casi completo para binarias con $1 < \log(a/\text{AU}) < 4$ y $q > 0.2$. Dentro de estos límites, la probabilidad de detección es uniforme, lo que permite que las distribuciones observadas sean tratadas como representativas de la población real sin correcciones de sesgo complejas.

Significancia y Reivindicaciones
El artículo afirma que el sondeo smash+ proporciona un censo casi completo de la multiplicidad de estrellas masivas en el régimen de separación intermedia. Al convertir los datos angulares a unidades físicas, el estudio confirma que las estrellas masivas se forman casi exclusivamente en sistemas múltiples, siendo la mayoría triples o sistemas de orden superior.

La principal contribución científica radica en caracterizar la distribución de la razón de masa como una función de la separación. Los resultados sugieren una transición en la formación o historia dinámica:

• En separaciones cercanas ($<100$ AU), la distribución de la razón de masa es relativamente plana, lo que potencialmente indica una formación correlacionada o interacciones dinámicas tempranas.
• En separaciones anchas ($>1000$ AU), la distribución se vuelve pronunciada, y la ausencia de pares de alta razón de masa sugiere que los componentes a estas distancias se formaron de manera casi independiente, consistente con la captura o fragmentación de núcleos separados.

Los autores concluyen que, si bien sus valores de $\kappa$ derivados son compatibles tanto con las binarias espectroscópicas como con las distribuciones más pronunciadas de Moe & Di Stefano (2017), se requiere un conjunto de datos más grande para resolver definitivamente la distribución de la razón de masa para binarias de periodo intermedio. Esta distribución es crítica para los modelos de síntesis de población, particularmente respecto a los resultados de la transferencia de masa de Caso B. Además, la calibración espectrofotométrica desarrollada en este trabajo sirve como una herramienta para identificar estrellas con historias evolutivas anómalas (por ejemplo, productos de fusión o remanentes de transferencia de masa) al compararlas contra estimaciones independientes de distancia o masa.