Temporal evolution of the circumstellar disk orientation in the transient X-ray pulsar GRO J1008-57

Resumen Técnico: Evolución temporal de la orientación del disco circumestelar en el púlsar de rayos X transitorio GRO J1008−57

Planteamiento del Problema
Las binarias Be/rayos X (BeXRB) exhiben dos tipos distintos de erupciones de rayos X: erupciones periódicas de Tipo I (normales), que típicamente ocurren cerca del periastrón, y erupciones esporádicas de Tipo II (gigantes). Mientras que las erupciones de Tipo I se atribuyen generalmente al acercamiento más cercano de la estrella de neutrones al disco de acreción circumestelar de la estrella Be, las fases orbitales de estas erupciones no siempre permanecen estables. En algunos sistemas, como EXO 2030+375, la fase orbital de las erupciones de Tipo I exhibe variaciones complejas a largo plazo, incluidas alternancias periódicas entre retrasos y adelantos de fase, a menudo desencadenadas por erupciones de Tipo II. Se cree que estas variaciones reflejan la evolución estructural y física del disco Be. Sin embargo, los mecanismos que impulsan estos desplazamientos de fase —ya sea precesión del disco, oscilaciones de excentricidad (Kozai-Lidov) o perturbaciones de densidad— siguen siendo objeto de debate. GRO J1008−57, una BeXRB con un largo período orbital (249.48 días) y una historia de múltiples erupciones de Tipo II, sirve como un caso de estudio crítico para investigar si los mecanismos de evolución de fase observados en sistemas de período más corto se aplican a sistemas con diferentes escalas de tiempo dinámicas.

Metodología
Los autores analizaron datos de monitoreo de rayos X a largo plazo de GRO J1008−57 obtenidos de los instrumentos Swift/BAT (15–150 keV) y MAXI/GSC (2–30 keV), cubriendo el período de 2005 a 2026.

1. Análisis de Curvas de Luz: Las curvas de luz diarias se utilizaron para identificar 31 erupciones de Tipo I (Swift) y 24 (MAXI), así como cuatro erupciones principales de Tipo II ocurridas en 2012, 2015, 2017 y 2020.
2. Modelado de Temporización: Los perfiles temporales de las erupciones de Tipo I se modelaron utilizando una función simplificada de ascenso rápido y decaimiento exponencial (FRED) para determinar los tiempos precisos de inicio ($t_{start}$) y pico ($t_{peak}$).
3. Cálculo de Fase: Las fases orbitales se calcularon con respecto a los tiempos de paso por el periastrón derivados de la efeméride de Kühnel et al. (2013).
4. Modelado Evolutivo: La evolución a largo plazo de las fases orbitales de los inicios de las erupciones se ajustó contra tres modelos: deriva lineal, deriva cuadrática y un modelo tipo escalón. Se empleó el Criterio de Información de Akaike (AIC) para seleccionar estadísticamente el modelo de mejor ajuste.
5. Energetica: La liberación total de energía (fluencia) de cada erupción de Tipo I se estimó integrando las curvas del modelo ajustado sobre una ventana de $\pm0.15$ períodos orbitales alrededor del pico.

Resultados Clave

1. Evolución de Fase Tipo Escalón: El análisis revela que las fases orbitales de las erupciones de Tipo I en GRO J1008−57 no siguen una deriva lineal o cuadrática continua. En su lugar, exhiben una evolución tipo escalón. Las fases permanecen estables durante muchos períodos orbitales pero experimentan saltos abruptos de pequeña amplitud (que oscilan entre $\sim -0.003$ y $-0.016$ en fase) que coinciden temporalmente con cada una de las cuatro erupciones de Tipo II observadas.
2. Selección de Modelo: El modelo tipo escalón proporcionó un ajuste estadísticamente superior en comparación con los modelos lineal y cuadrático para los datos de Swift y MAXI, con valores de $\Delta$AIC que indican una significancia superior a $5\sigma$.
3. Ciclos Energéticos: La liberación total de energía de las erupciones de Tipo I muestra un patrón sistemático: un aumento gradual previo a una erupción de Tipo II, seguido de un declive rápido y una recuperación gradual posterior. Esto sugiere ciclos de acumulación de masa del disco, agotamiento y reconstrucción impulsados por eventos de Tipo II.
4. Estabilidad entre Saltos: A diferencia del sistema EXO 2030+375, que muestra desplazamientos de fase de gran amplitud y períodos de recuperación a largo plazo (cientos de días), GRO J1008−57 no muestra evidencia de recuperación de fase prolongada ni de deriva continua entre eventos de Tipo II. El sistema regresa a un nuevo estado de fase estable casi inmediatamente después de que concluye la erupción de Tipo II.

Significado y Mecanismo Propuesto
El artículo argumenta que la evolución tipo escalón observada en GRO J1008−57 desafía los modelos existentes propuestos para otras BeXRB, como el mecanismo de oscilación Kozai-Lidov (KL), las oscilaciones globales de disco de un solo brazo o la precesión a largo plazo del disco.

• Rechazo de las Oscilaciones KL: La escala de tiempo estimada para las oscilaciones KL en GRO J1008−57 es de aproximadamente 79 años, muy por encima de la base de observación. Además, los mecanismos KL predicen típicamente grandes desplazamientos de fase y crecimiento de la excentricidad, lo cual no se observa aquí.
• Rechazo de Precesión/Oscilaciones: La ausencia de períodos de recuperación de fase de larga duración (observados en EXO 2030+375) descarta la precesión a largo plazo del disco o las perturbaciones de densidad globales como el impulsor principal.

Explicación Propuesta:
Los autores proponen que el comportamiento único de GRO J1008−57 está impulsado por su largo período orbital (249.48 días), lo que resulta en interacciones infrecuentes entre la estrella de neutrones y el disco circumestelar.

• Ciclo de Reconstrucción del Disco: Las erupciones de Tipo II se interpretan como eventos que agotan el disco Be. Debido al largo período orbital, el sistema tiene tiempo suficiente entre interacciones para que el disco reconstruya completamente su densidad y restaure una estructura geométrica similar a su estado previo a la erupción.
• Saltos Pequeños: En consecuencia, después de cada evento de Tipo II, el disco se asienta en una configuración ligeramente diferente pero estable, causando un desplazamiento pequeño y abrupto en la fase orbital de las erupciones de Tipo I subsiguientes. El sistema permanece luego en este nuevo estado estable hasta el próximo evento de agotamiento.
• Contraste con EXO 2030+375: Por el contrario, EXO 2030+375 tiene un período orbital corto ($\sim46$ días), lo que lleva a interacciones frecuentes que impiden que el disco se relaje completamente, permitiendo así que las perturbaciones se acumulen y impulsen grandes desplazamientos de fase continuos.

Conclusión
Este estudio amplía la comprensión de los mecanismos de erupción de las BeXRB al demostrar que la evolución temporal de las fases de las erupciones de Tipo I no es universal, sino que depende de la dinámica orbital del sistema. La evolución tipo escalón en GRO J1008−57 proporciona evidencia de un ciclo de agotamiento y reconstrucción del disco impulsado por erupciones de Tipo II, donde el largo período orbital permite que el disco reestablezca una geometría estable entre eventos, resultando en saltos de fase discretos en lugar de una deriva continua.