Optical Super-orbital Modulation of SMC X-1: Disk Precession and a Revised Pulsar Mass

Resumen Técnico: Modulación Óptica Superorbital de SMC X-1: Precesión del Disco y una Masa Revisada del Púlsar

Planteamiento del Problema
La determinación observacional del límite inferior de masa para las estrellas de neutrones es crítica para restringir la física de las supernovas de colapso del núcleo (CCSNe). El sistema binario de rayos X de masa alta SMC X-1 es un objeto pivotal en este contexto, ya que las estimaciones de masa anteriores situaban la masa de su púlsar cerca o potencialmente por debajo del límite teórico inferior de $\approx 1.1 M_\odot$. Estas estimaciones previas se basaron en determinar la relación de masas binaria y la inclinación orbital, principalmente mediante mediciones espectroscópicas de velocidad radial (RV) de la estrella donante y modelado de la modulación elipsoidal óptica. Sin embargo, SMC X-1 exhibe modulaciones ópticas superorbitales complejas (variaciones sincronizadas con el ciclo superorbital de rayos X de $\sim 40\text{--}60$ días) que no se tuvieron en cuenta completamente en los modelos estáticos anteriores. Las observaciones desde tierra carecían de la precisión fotométrica y el muestreo temporal necesarios para resolver estas variaciones dentro de los ciclos orbitales individuales, lo que llevó a posibles errores sistemáticos en los parámetros binarios derivados.

Metodología
Los autores analizaron curvas de luz ópticas de alta cadencia de SMC X-1 obtenidas por el Satélite de Sondeo de Exoplanetas en Tránsito (TESS) a través de seis sectores, combinadas con datos de monitoreo de rayos X a largo plazo del Monitor de Imagen de Rayos X de Todo el Cielo (MAXI).

1. Análisis Observacional: El estudio identificó variaciones sistemáticas en las curvas de luz orbitales ópticas sincronizadas con la modulación superorbital de rayos X. Específicamente, la profundidad del mínimo en la conjunción inferior y la asimetría de la curva de luz de doble pico variaban dependiendo del estado de rayos X (ascendente, decreciente, activo o inactivo).
2. Desarrollo del Modelo: Se desarrolló un código de modulación elipsoidal modificado en Python (utilizando NumPy, SciPy, Astropy y CuPy para la paralelización). A diferencia de los modelos estáticos, este código incorpora un disco de acreción precesante que altera:
   • La geometría de la irradiación de rayos X en la superficie de la estrella donante.
   • La geometría del reprocesamiento óptico/UV (OUV) en la superficie del disco.
   • El ocultamiento mutuo entre la donante, el disco y el púlsar.
     El modelo calcula la temperatura y la radiación de cuerpo negro de los elementos de superficie (SE) tanto en la donante como en el disco, teniendo en cuenta el oscurecimiento gravitacional, el calentamiento por rayos X y la irradiación OUV.
3. Estimación de Parámetros: Se realizó un muestreo de Cadena de Markov Monte Carlo (MCMC) para ajustar el modelo a las curvas de luz ópticas de TESS y de rayos X de MAXI. El análisis probó dos escenarios para la relación de masas ($q = M_{NS}/M_s$): uno basado en las RV observadas ($\log q = -1.17$) y un segundo donde $q$ se aumentó un 20% ($\log q = -1.09$) para probar la hipótesis de que la RV de la donante observada está subestimada.

Resultados Clave

• Correlaciones: El estudio confirmó fuertes correlaciones entre el flujo superorbital de rayos X y las características ópticas: el mínimo de conjunción inferior es superficial durante los estados "inactivos" de rayos X y profundo durante los estados "activos"; la asimetría del doble pico se correlaciona con la derivada del flujo de rayos X (por ejemplo, el pico izquierdo es más alto durante la fase ascendente de rayos X).
• Validación del Modelo: El modelo elipsoidal modificado reprodujo con éxito las curvas de luz ópticas y de rayos X observadas, proporcionando evidencia observacional de que la modulación superorbital es impulsada por un disco de acreción precesante. El análisis sugiere que la precesión del disco es retrógrada.
• Sesgo en la Velocidad Radial: La intensa irradiación de rayos X de la estrella donante desplaza el centro de emisión óptica lejos del centro gravitacional. El modelo predice que este efecto hace que la amplitud de la velocidad radial observada esté subestimada en aproximadamente un 20% (la curva de RV aparente parece "recortada" en comparación con una onda sinusoidal verdadera).
• Estimación de Masa Revisada:
  • Utilizando la relación de masas estándar ($\log q = -1.17$), la masa del púlsar derivada es $M_{NS} = 1.11 \pm 0.07 M_\odot$, consistente con estudios anteriores.
  • Corrigiendo por la subestimación del 20% en la RV (aumentando $q$ a $\log q = -1.09$), la masa revisada del púlsar es $M_{NS} = 1.35 \pm 0.10 M_\odot$.
• Discrepancias: El modelo requirió una luminosidad de rayos X ($L_X \approx 4 \times 10^{39}$ erg s$^{-1}$) significativamente más alta que los valores observados típicos para reproducir la intensidad de la modulación óptica. Los autores sugieren que esto podría implicar emisión anisotrópica de rayos X (haz dirigido hacia la donante) o desalineación espín-órbita, aunque esta última es difícil de conciliar con la órbita circular del sistema a menos que la acreción de masa durante una fase de envoltura común fuera ineficiente.

Significado
El artículo afirma que las variaciones sistemáticas en la curva de luz óptica, impulsadas por el disco precesante, no se habían modelado previamente de una manera que afectara la determinación de la masa. Al desarrollar un modelo que tenga en cuenta la geometría cambiante de la irradiación, los autores demuestran que la velocidad radial de la donante probablemente ha sido subestimada debido a desplazamientos inducidos por la irradiación en el centro de emisión. Corregir este efecto arroja una masa revisada del púlsar de $\approx 1.35 M_\odot$. Este valor sitúa al púlsar SMC X-1 bien dentro del rango típico para las estrellas de neutrones, eliminando la tensión con los límites inferiores teóricos para las supernovas de colapso del núcleo y sugiriendo que las estimaciones de masa anteriores estaban sesgadas por efectos de irradiación no corregidos.