Constraints on cyclotron features and accretion regime in the high-mass X-ray binary 4U 1700-37 from NuSTAR

Resumen Técnico: Restricciones en características de ciclotrón y régimen de acreción en el binario de rayos X de alta masa 4U 1700–37 mediante NuSTAR

1. Problema de Investigación

El sistema 4U 1700–37 es un binario de rayos X de alta masa (HMXB) alimentado por vientos, compuesto por una estrella supergigante O6.5 Iaf+ (HD 153919) y un objeto compacto, presumiblemente una estrella de neutrones (EN). A pesar de ser una fuente brillante y persistente, existen incertidumbres fundamentales sobre sus propiedades físicas:

• Campo magnético: No se ha medido directamente debido a la falta de detección firme de pulsaciones coherentes y de características de resonancia de dispersión de ciclotrón (CRSFs).
• Inconsistencias previas: Estudios anteriores con NuSTAR reportaron candidatos a CRSFs marginales (alrededor de ~16 keV), pero la estabilidad, energía y existencia real de estas líneas son debatidas.
• Dependencia del modelo: Las mediciones previas se basaron en épocas únicas, suposiciones de un solo continuo y calibraciones antiguas, lo que dificultaba evaluar la robustez de las características reportadas frente a diferentes modelos espectrales.

El objetivo principal fue realizar un análisis homogéneo de todas las observaciones disponibles de NuSTAR para caracterizar las propiedades temporales y espectrales, probar la robustez de los candidatos a CRSFs y restringir el campo magnético y el régimen de acreción del objeto compacto.

2. Metodología

Los autores realizaron un análisis temporal y espectral homogéneo de dos observaciones de NuSTAR (2016 y 2022) y datos complementarios de Swift/BAT:

• Reducción de Datos: Se utilizaron NuSTARDAS y la base de datos de calibración más reciente (v.20240229). Se aplicaron correcciones bariocéntricas y filtros estrictos para eliminar eventos de la Anomalía del Atlántico Sur (SAA). Se modelaron explícitamente las características de absorción instrumental (bordes K de Zn y Cd) presentes en los detectores.
• Análisis Temporal: Se examinaron curvas de luz en bandas suave (3–8 keV) y dura (8–79 keV), se calcularon espectros de densidad de potencia (PDS) y periodogramas de Lomb-Scargle. Se buscó pulsaciones en el rango de 0.01–5000 Hz utilizando la herramienta HENaccelsearch y se establecieron límites superiores a la fracción pulsada (PF) mediante simulaciones.
• Análisis Espectral: Se ajustaron espectros promediados en el tiempo y espectros resueltos por intensidad (dividiendo la observación en bandas de baja, media y alta tasa de conteo) utilizando XSPEC.
  • Se probaron cinco modelos de continuo diferentes (incluyendo bbody+cutoffpl, nthcomp, thcomp, highecut, npex).
  • Se buscaron CRSFs utilizando el modelo gabs (y cyclabs) sobre los continuos óptimos.
  • Validación de Significancia: Se realizaron más de $10^5$ simulaciones espectrales para cuantificar la significancia de las líneas y su dependencia del modelo de continuo, evitando falsos positivos derivados de un ajuste global deficiente.

3. Contribuciones Clave

• Análisis Homogéneo Multi-época: Es el primer estudio que compara sistemáticamente dos épocas de NuSTAR con calibraciones actualizadas y múltiples modelos de continuo para 4U 1700–37.
• Caracterización de la Dependencia del Modelo: Se demuestra que la detección y la energía de las características de absorción candidatas son altamente dependientes del modelo de continuo subyacente.
• Restricciones Cuantitativas: Se establecen límites superiores rigurosos a la fracción pulsada y se derivan rangos para el campo magnético y el periodo de equilibrio de rotación, incluso sin una detección firme de CRSF.
• Reconciliación de Resultados Previos: Se explica por qué estudios anteriores reportaron líneas a diferentes energías (~16 keV vs ~38 keV), atribuyendo las discrepancias a la elección del modelo de continuo y a la variabilidad de la fuente.

4. Resultados Principales

• Pulsaciones: No se detectaron pulsaciones coherentes en el rango de 0.01–5000 s. Se establecieron límites superiores de fracción pulsada (PF) del 1.36% (FPMA) y 1.47% (FPMB) al 90% de confianza.
• Espectro Continuo: El espectro se describe mejor con un modelo de cuerpo negro absorbido más una ley de potencia con corte (bbody+cutoffpl).
  • Temperatura del cuerpo negro: ~1.64 keV.
  • Índice de fotones: ~1.05.
  • Energía de corte: ~17.1 keV.
  • Se detectó una línea de emisión de Fe K$\alpha$ en ~6.4 keV.
• Candidatos a CRSFs:
  • Se identificaron residuos de absorción alrededor de ~19.5 keV y ~52 keV.
  • Bajo el modelo de continuo óptimo (bbody+cutoffpl), la significancia de estas características es < 3$\sigma$ (probabilidades de ~95%), por lo que no constituyen detecciones firmes.
  • Sin embargo, con modelos de continuo alternativos (con peores ajustes globales), la significancia puede superar el 99.9%.
  • La energía aparente del centro de la línea muestra un desplazamiento con el flujo (espectroscopía resuelta por intensidad).
• Parámetros Físicos Derivados (Interpretación):
  • Si se interpretan las características como líneas de ciclotrón, el campo magnético de la estrella de neutrones se restringe a $B \sim (1.7 - 4.4) \times 10^{12}$ G.
  • Utilizando el modelo de acreción cuasi-esférica subsónica (Shakura et al. 2012) y el campo magnético derivado, se estima un periodo de equilibrio de rotación ($P_{eq}$) de ~1.9 ks (1900 segundos).
  • La luminosidad observada ($L_X \approx 1.7 \times 10^{36}$ erg s$^{-1}$) es consistente con un régimen de acreción subsónica estable.

5. Significancia e Implicaciones

• Naturaleza del Objeto: Los resultados apoyan la clasificación de 4U 1700–37 como un sistema con una estrella de neutrones de campo magnético "normal" (no magnetar, no campo muy bajo) en un régimen de acreción por viento subsónico.
• Posición en el Diagrama de Corbet: El periodo de equilibrio estimado (~1.9 ks) sitúa al sistema en la región de HMXBs con supergigantes del diagrama de Corbet, lo cual es consistente con la falta de detección de pulsaciones (el periodo de rotación real podría estar cerca del equilibrio, reduciendo el torque y la visibilidad de pulsaciones).
• Metodología Futura: El estudio subraya que en sistemas alimentados por vientos como 4U 1700–37, la afirmación de la presencia de líneas de ciclotrón requiere un escrutinio cuidadoso de los residuos espectrales más allá de las estadísticas globales ($\chi^2$). Se necesitan observaciones más profundas y de banda ancha para confirmar las características con un umbral de >3$\sigma$.
• Resolución de Discrepancias: El trabajo explica las mediciones previas dispares de energía de línea como un artefacto de la modelización del continuo y la variabilidad de la fuente, proporcionando un marco unificado para interpretar los datos de rayos X duros.

En conclusión, aunque no se logra una detección definitiva de CRSFs, el análisis proporciona las restricciones más sólidas hasta la fecha sobre el campo magnético y la física de acreción de 4U 1700–37, demostrando la importancia crítica de la modelización espectral adecuada en la astrofísica de rayos X.