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

Resumo Técnico: Evolução temporal da orientação do disco circumestelar no pulsar X transitório GRO J1008−57

Enunciação do Problema
Binárias Be/X (BeXRBs) exibem dois tipos distintos de erupções em raios-X: erupções periódicas do Tipo I (normais), tipicamente ocorrendo próximas ao periastro, e erupções esporádicas do Tipo II (gigantes). Enquanto as erupções do Tipo I são geralmente atribuídas à aproximação mais próxima da estrela de nêutrons ao disco de ejeção circumestelar da estrela Be, as fases orbitais dessas erupções nem sempre permanecem estáveis. Em alguns sistemas, como EXO 2030+375, a fase orbital das erupções do Tipo I exibe variações complexas de longo prazo, incluindo alternâncias periódicas entre atrasos e adiantamentos de fase, frequentemente desencadeadas por erupções do Tipo II. Acredita-se que essas variações reflitam a evolução estrutural e física do disco Be. No entanto, os mecanismos que impulsionam essas mudanças de fase — seja precessão do disco, oscilações de excentricidade (Kozai-Lidov) ou perturbações de densidade — permanecem debatidos. GRO J1008−57, uma BeXRB com longo período orbital (249,48 dias) e histórico de múltiplas erupções do Tipo II, serve como um estudo de caso crítico para investigar se os mecanismos de evolução de fase observados em sistemas de período mais curto se aplicam a sistemas com diferentes escalas de tempo dinâmicas.

Metodologia
Os autores analisaram dados de monitoramento de raios-X de longo prazo de GRO J1008−57 obtidos pelos instrumentos Swift/BAT (15–150 keV) e MAXI/GSC (2–30 keV), cobrindo o período de 2005 a 2026.

1. Análise de Curvas de Luz: Curvas de luz diárias foram utilizadas para identificar 31 erupções do Tipo I (Swift) e 24 (MAXI), bem como quatro grandes erupções do Tipo II ocorridas em 2012, 2015, 2017 e 2020.
2. Modelagem de Temporização: Os perfis temporais das erupções do Tipo I foram modelados usando uma função simplificada de crescimento rápido e decaimento exponencial (FRED) para determinar tempos precisos de início ($t_{start}$) e pico ($t_{peak}$).
3. Cálculo de Fase: Fases orbitais foram calculadas relativas aos tempos de passagem pelo periastro derivados da efeméride de Kühnel et al. (2013).
4. Modelagem Evolutiva: A evolução de longo prazo das fases orbitais dos inícios das erupções foi ajustada contra três modelos: deriva linear, deriva quadrática e um modelo em degrau. O Critério de Informação de Akaike (AIC) foi empregado para selecionar estatisticamente o melhor modelo ajustado.
5. Energetica: A liberação total de energia (fluência) de cada erupção do Tipo I foi estimada integrando as curvas do modelo ajustado sobre uma janela de $\pm0,15$ períodos orbitais em torno do pico.

Principais Resultados

1. Evolução de Fase em Degraus: A análise revela que as fases orbitais das erupções do Tipo I em GRO J1008−57 não seguem uma deriva linear ou quadrática contínua. Em vez disso, elas exibem uma evolução em degraus. As fases permanecem estáveis ao longo de muitos períodos orbitais, mas sofrem saltos abruptos de pequena amplitude (variando de $\sim -0,003$ a $-0,016$ em fase) que coincidem temporalmente com cada uma das quatro erupções do Tipo II observadas.
2. Seleção de Modelo: O modelo em degrau forneceu um ajuste estatisticamente superior em comparação com os modelos linear e quadrático para ambos os dados Swift e MAXI, com valores de $\Delta$AIC indicando uma significância superior a $5\sigma$.
3. Ciclos Energéticos: A liberação total de energia das erupções do Tipo I exibe um padrão sistemático: um aumento gradual levando a uma erupção do Tipo II, seguido por um declínio rápido e uma recuperação gradual subsequente. Isso sugere ciclos de acúmulo de massa do disco, esgotamento e reconstrução impulsionados por eventos do Tipo II.
4. Estabilidade Entre Saltos: Ao contrário do sistema EXO 2030+375, que mostra grandes deslocamentos de fase de amplitude e períodos de recuperação de longo prazo (centenas de dias), GRO J1008−57 não mostra evidências de recuperação de fase prolongada ou deriva contínua entre eventos do Tipo II. O sistema retorna a um novo estado de fase estável quase imediatamente após o término da erupção do Tipo II.

Significado e Mecanismo Proposto
O artigo argumenta que a evolução em degraus observada em GRO J1008−57 desafia modelos existentes propostos para outras BeXRBs, como o mecanismo de oscilação Kozai-Lidov (KL), oscilações globais de disco de um braço ou precessão de disco de longo prazo.

• Rejeição de Oscilações KL: A escala de tempo estimada para oscilações KL em GRO J1008−57 é de aproximadamente 79 anos, muito superior à base de observação. Além disso, mecanismos KL tipicamente preveem grandes deslocamentos de fase e crescimento de excentricidade, o que não é observado aqui.
• Rejeição de Precessão/Oscilações: A ausência de períodos de recuperação de fase de longa duração (observados em EXO 2030+375) argumenta contra precessão de disco de longo prazo ou perturbações globais de densidade como o motor principal.

Explicação Proposta:
Os autores propõem que o comportamento único de GRO J1008−57 é impulsionado pelo seu longo período orbital (249,48 dias), o que resulta em interações infrequentes entre a estrela de nêutrons e o disco circumestelar.

• Ciclo de Reconstrução do Disco: Erupções do Tipo II são interpretadas como eventos que esgotam o disco Be. Devido ao longo período orbital, o sistema tem tempo suficiente entre interações para que o disco reconstrua totalmente sua densidade e restaure uma estrutura geométrica semelhante ao seu estado pré-erupção.
• Pequenos Saltos: Consequentemente, após cada evento do Tipo II, o disco se estabelece em uma configuração ligeiramente diferente, mas estável, causando um deslocamento pequeno e abrupto na fase orbital das erupções subsequentes do Tipo I. O sistema então permanece nesse novo estado estável até o próximo evento de esgotamento.
• Contraste com EXO 2030+375: Em contraste, EXO 2030+375 possui um período orbital curto ($\sim46$ dias), levando a interações frequentes que impedem o relaxamento completo do disco, permitindo assim que perturbações se acumulem e impulsionem grandes deslocamentos de fase contínuos.

Conclusão
Este estudo amplia a compreensão dos mecanismos de erupção em BeXRBs ao demonstrar que a evolução temporal das fases das erupções do Tipo I não é universal, mas depende da dinâmica orbital do sistema. A evolução em degraus em GRO J1008−57 fornece evidências de um ciclo de esgotamento e reconstrução do disco impulsionado por erupções do Tipo II, onde o longo período orbital permite que o disco reestabeleça uma geometria estável entre eventos, resultando em saltos de fase discretos em vez de deriva contínua.