A Kinematic Study of Wolf-Rayet Stars at the Galactic Center I: Binary Candidates and Constraints on the Binary Fraction

Resumen Técnico: Un Estudio Cinemático de Estrellas Wolf-Rayet en el Centro Galáctico I: Candidatos Binarios y Restricciones sobre la Fracción Binaria

Planteamiento del Problema
El origen de la población estelar joven en el Centro Galáctico, particularmente las estrellas Wolf-Rayet (WR) ubicadas dentro de 0.5 pc de la agujero negro supermasivo Sgr A*, sigue siendo poco comprendido. Los modelos estándar sugieren que la cizalla de marea de Sgr A* debería impedir la formación estelar a estas distancias, lo que requiere escenarios de formación alternativos, como la formación in situ dentro de un disco de acreción. La fracción binaria de estas estrellas sirve como un diagnóstico crítico para probar estos modelos de formación; por ejemplo, la formación estelar en un disco de acreción predice una fracción binaria más alta debido a las altas tasas de fragmentación. Además, la fracción binaria influye en la evolución dinámica del cúmulo a través de fusiones y evaporación, y el movimiento binario puede sesgar las mediciones cinemáticas utilizadas para determinar la pertenencia estelar a estructuras específicas. Si bien la fracción binaria de la población estelar joven general en el Centro Galáctico ha sido estimada mediante fotometría (>71%), las restricciones espectroscópicas específicamente para la población WR han estado limitadas por tamaños de muestra pequeños (típicamente $\le$ 8 estrellas) y líneas de base temporales cortas.

Metodología
Este estudio presenta el estudio cinemático con la línea de base temporal más larga (1994–2024) de estrellas WR en el Centro Galáctico hasta la fecha.

• Ensamblaje de la Muestra: Los autores construyeron una muestra de 27 estrellas y candidatos WR dentro de 0.5 pc de Sgr A*, que comprende 18 estrellas WN y 9 estrellas WC. Esta muestra triplica el número de estrellas WR incluidos en búsquedas de binarias espectroscópicas previas en la región.
• Observaciones: El estudio utiliza 110 nuevas observaciones de unidad de campo integral tomadas entre 2006 y 2024, principalmente utilizando el espectrógrafo OSIRIS en el Observatorio W. M. Keck (101 observaciones) y el espectrógrafo NIFS en el telescopio Gemini North (9 observaciones). Estas se combinaron con 70 mediciones de velocidad radial ($v_z$) previamente publicadas en la literatura.
• Reducción y Análisis de Datos: Los espectros se redujeron utilizando pipelines estándar (OSIRIS DRP y Nifty4Gemini). Las velocidades radiales se derivaron mediante correlación cruzada contra plantillas empíricas construidas a partir de los propios espectros observados para minimizar la sensibilidad a parámetros atmosféricos específicos. Para la binaria de doble línea IRS 16SW, la $v_z$ se determinó mediante el ajuste gaussiano de la línea de absorción He I de 2.112 µm.
• Presupuesto de Error: Las incertidumbres estadísticas se derivaron mediante simulaciones de Monte Carlo. Las incertidumbres sistemáticas derivadas de la variabilidad espectral impulsada por el viento se añadieron en cuadratura: 15 km/s para las estrellas WN (incluyendo Ofpe/WN9 e híbridas) y 3 km/s para las estrellas WC, basadas en caracterizaciones previas de la variabilidad del viento.
• Identificación de Binarias: Los candidatos se identificaron como estrellas que exhiben una variación de $v_z$ entre dos épocas cualesquiera con una significancia mayor a $3\sigma$, tras corregir por la aceleración orbital de las estrellas alrededor de Sgr A*.
• Inferencia Estadística: Para inferir la fracción binaria intrínseca ($f_{bin}$), los autores realizaron 1000 simulaciones. Estas simulaciones asignaron aleatoriamente el estatus de binaria a las 27 estrellas basándose en una $f_{bin}$ de prueba, generaron propiedades físicas binarias (periodo, razón de masa, excentricidad) basadas en distribuciones observadas en estrellas WR de campo (Sana et al. 2012; Dsilva et al. 2022, 2023), y simularon mediciones de $v_z$ que coincidan con la cadencia observacional y las incertidumbres de los datos reales. El número de binarias simuladas detectadas con una significancia $>3\sigma$ se comparó con el número real de detecciones para derivar la distribución posterior de $f_{bin}$.

Contribuciones Clave y Resultados

• Nuevos Candidatos Binarios: El estudio identifica dos nuevos candidatos binarios entre las estrellas WR: S8-181 (WN8) e IRS 13E4 (WC9).
  • S8-181: Muestra una variación de $v_z$ de $139 \pm 42$ km/s ($3.4\sigma$) entre los espectros de 2022 y 2024. La variación excede los límites de variabilidad impulsada por el viento y la aceleración gravitacional de Sgr A*.
  • IRS 13E4: Muestra una variación de $v_z$ de $71 \pm 15$ km/s ($4.8\sigma$) sobre una línea de base de 14 años. La variación excede la variabilidad de viento esperada para estrellas WC ($\sim 10$ km/s) y la aceleración de Sgr A*.
• Binarias Confirmadas: El estudio confirma tres binarias previamente conocidas: IRS 16SW, IRS 16NE y S4-258.
• Fracción Binaria: Basándose en la detección de 5 binarias (3 conocidas, 2 nuevas) de las 27 estrellas de la muestra, y contabilizando la sensibilidad experimental mediante simulación, los autores infieren una fracción binaria espectroscópica para las estrellas WR del Centro Galáctico de $f_{bin} = 0.56 \pm 0.18$.
• Restricciones Combinadas: Al combinarse con los resultados de la fracción binaria fotométrica de Gautam et al. (2024) para la población estelar joven más amplia, el estudio deriva una fracción binaria combinada de $f_{bin} = 0.69 \pm 0.17$.

Significancia y Reivindicaciones
El artículo afirma que la fracción binaria derivada de $0.56 \pm 0.18$ para las estrellas WR del Centro Galáctico es consistente con:

1. Estudios fotométricos previos de estrellas jóvenes en el Centro Galáctico.
2. Las fracciones binarias de las estrellas WR de campo (específicamente estrellas WN) y estrellas O en cúmulos jóvenes de la Vía Láctea.
3. La sugerencia de que la fracción binaria aumenta con la distancia radial de Sgr A*, dado que la fracción medida es significativamente mayor que el límite superior de $f_{bin} \le 0.17$ encontrado para las estrellas tipo B dentro de 0.04 pc del agujero negro (Chu et al. 2023).

Los autores enfatizan que, si bien su resultado respalda la alta fracción binaria predicha por los modelos de formación de disco de acreción, la inferencia depende de supuestos respecto a la distribución de las propiedades físicas de las binarias (periodo, razón de masa) que no están independientemente restringidas para la población del Centro Galáctico. El estudio también señala que las estrellas WC productoras de polvo, a menudo asociadas con binarias de vientos colisionantes, no mostraron universalmente variaciones de velocidad en esta muestra, lo que llevó a los autores a excluir conservadoramente la producción de polvo como un indicador único de binaridad en su análisis. El trabajo representa la muestra espectroscópica más grande de estrellas WR del Centro Galáctico hasta la fecha, proporcionando una restricción cinemática robusta sobre la población binaria en este entorno extremo.