Cyclotron lines in subcritical X-ray pulsars: Monte Carlo simulations reveal the origin of the observed variability

Resumen Técnico: Líneas Ciclotrón en Púlsares de Rayos X Subcríticos

Enunciado del Problema
Los púlsares de rayos X (XRPs) exhiben características de dispersión resonante ciclotrón (CRSF) en sus espectros, las cuales sirven como diagnósticos directos de la intensidad del campo magnético de la estrella de neutrones (NS). Las observaciones revelan una variabilidad significativa en la energía del centroide de la CRSF ($E_{CRSF}$) y en el ancho de la línea ($\sigma_{CRSF}$) en función de la luminosidad en rayos X ($L$) y de la fase del pulso. En el régimen de luminosidad subcrítica ($L \lesssim 10^{37}$ erg s$^{-1}$), a menudo se observa una correlación positiva entre $E_{CRSF}$ y $L$, en contraste con la correlación negativa observada en fuentes supercríticas. Si bien los modelos cualitativos sugieren que esto surge del desplazamiento Doppler de fotones que se dispersan en el plasma en caída dentro del embudo de acreción, ha faltado un tratamiento riguroso y cuantitativo de transferencia radiativa de esta variabilidad. Específicamente, la interacción entre la dispersión resonante, los efectos de la presión de radiación sobre la velocidad del flujo y la geometría de observación requiere una simulación detallada para explicar las tendencias observadas.

Metodología
Los autores desarrollaron un código de Monte Carlo (MC) relativista para realizar cálculos detallados de transferencia radiativa dentro del embudo de acreción por encima del punto caliente de la NS. El modelo adopta el marco físico propuesto por Mushtukov et al. (2015b) para fuentes subcríticas:

• Geometría: Un canal de acreción cilíndrico semi-infinito alineado con el eje magnético ($B = B\hat{z}$).
• Perfiles de Plasma: Expresiones analíticas para la densidad de electrones ($n_e$) y la velocidad de flujo ($\beta$) se derivan de la conservación de masa y momento, incorporando la desaceleración del plasma en caída debido a la presión de radiación. El perfil de velocidad depende de la luminosidad relativa a una luminosidad crítica $L^*$.
• Transferencia Radiativa: Se inyectan fotones semilla isotrópicamente desde el punto caliente con un espectro de cuerpo negro. Los paquetes de fotones se propagan a través del plasma, experimentando dispersión Compton.
• Física de Dispersión: El código emplea secciones eficaces de dispersión magnética resonante aproximadas y promediadas en polarización (Loudas et al. 2021) válidas para $b \ll 1$ y $E \ll m_ec^2$. La dispersión se trata en el marco de reposo del electrón, aplicando transformaciones de Lorentz para tener en cuenta el movimiento de flujo del plasma.
• Condiciones de Contorno: Se probaron dos esquemas: "Reflexión" (reflexión geométrica en la base, conservando energía) y "Termalización" (absorción y reemisión). Los resultados principales utilizan el esquema de Reflexión.
• Parámetros: Las simulaciones se ejecutaron para una NS fiducial ($M=1.4 M_\odot$, $R=10^6$ cm, $B \approx 10^{12}$ G) y se aplicaron al caso específico del XRP GX 304−1. La luminosidad de acreción se varió en el rango $L/L^* \in (0, 1)$.

Contribuciones y Resultados Clave
El estudio proporciona la primera verificación cuantitativa mediante MC del modelo de Mushtukov et al. (2015b), obteniendo los siguientes resultados:

1. Morfología Espectral: Los espectros emergentes exhiben una CRSF prominente y asimétrica acompañada por un ala azul amplia. La CRSF está sistemáticamente desplazada al rojo con respecto a la energía ciclotrón clásica local ($E_c$). Este desplazamiento al rojo es causado por el efecto Doppler entre el marco de laboratorio y el marco del electrón en caída.
2. Dependencia de la Luminosidad: Tanto $E_{CRSF}$ como $\sigma_{CRSF}$ muestran una correlación positiva con la luminosidad para todos los ángulos de visión. A medida que aumenta $L$, la presión de radiación desacelera el flujo en caída de manera más efectiva, reduciendo la velocidad de flujo ($\beta$). Esta reducción de velocidad disminuye el desplazamiento al rojo Doppler, desplazando $E_{CRSF}$ más cerca de $E_c$. Simultáneamente, el aumento de la profundidad óptica a mayores luminosidades profundiza la línea y amplía el ala azul.
3. Dependencia del Ángulo de Visión: Las propiedades de la CRSF varían significativamente con el ángulo de visión $\theta$ (relativo al eje magnético).
   • $E_{CRSF}$ disminuye a medida que aumenta $\cos \theta$ (es decir, a medida que la vista se vuelve más "de frente"). Esto se debe a la dependencia del desplazamiento Doppler del ángulo entre el momento del fotón y la velocidad del flujo.
   • $\sigma_{CRSF}$ aumenta con $\cos \theta$. Las vistas de frente muestrean un rango más amplio de gradientes de velocidad a lo largo de la línea de visión, lo que conduce a un mayor ensanchamiento Doppler.
   • En consecuencia, el modelo predice una anticorrelación entre $E_{CRSF}$ y $\sigma_{CRSF}$ durante el ciclo del pulso (variabilidad resuelta en fase), mientras que los espectros promediados en fase muestran una correlación positiva con la luminosidad.
4. Aplicación a GX 304−1: Cuando se aplica a datos observacionales de GX 304−1, el modelo reproduce exitosamente las correlaciones positivas observadas de $E_{CRSF}$ y $\sigma_{CRSF}$ con la luminosidad en casi un orden de magnitud. El mejor ajuste se logra para ángulos de visión donde el embudo de acreción se ve predominantemente de perfil ($\cos \theta \approx 0.0-0.1$), consistente con configuraciones geométricas específicas (por ejemplo, un rotador casi alineado a alta inclinación o un rotador casi ortogonal a baja inclinación).

Significado y Afirmaciones
El artículo afirma que la dispersión resonante de la radiación que se propaga a través del plasma acreciendo en el embudo magnético proporciona una explicación natural y cuantitativa para la variabilidad observada de las CRSF en XRPs subcríticos. Los autores enfatizan que sus resultados:

• Validan la hipótesis de que el canal de acreción es la región formadora de líneas en fuentes subcríticas.
• Demuestran que el movimiento de flujo y el desplazamiento Doppler resultante son los impulsores principales de la variabilidad dependiente de la luminosidad, sin requerir la formación de una onda de choque sin colisiones (CS) como el sitio principal de formación de líneas.
• Destacan el papel crítico de la geometría de observación en la interpretación de datos resueltos en fase, prediciendo anticorrelaciones específicas entre la energía y el ancho de la línea que difieren de las tendencias promediadas en fase.
• Distinguen su escenario del modelo de CS (donde la formación de líneas ocurre entre la onda de choque y la superficie), señalando que el escenario de CS no predeciría la fuerte variabilidad dependiente de la fase observada en sus simulaciones.

Los autores señalan limitaciones, incluido el uso de una geometría cilíndrica simplificada, perfiles analíticos 1D, y el descuido de efectos relativistas generales y polarización. Sin embargo, argumentan que estas simplificaciones no alteran las conclusiones fundamentales sobre el origen de la variabilidad. El trabajo sirve como una prueba de concepto para el uso de la dispersión resonante en el embudo de acreción para modelar espectros de XRP subcríticos.