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

Resumo Técnico: Modulação Óptica Superorbital de SMC X-1: Precessão do Disco e uma Massa de Pulsar Revisada

Enunciação do Problema
A determinação observacional do limite inferior de massa para estrelas de nêutrons é crítica para restringir a física das supernovas de colapso do núcleo (CCSNe). O binário de raios-X de alta massa SMC X-1 é um objeto pivotal neste contexto, porque estimativas anteriores de massa colocavam a massa do seu pulsar próxima ou potencialmente abaixo do limite inferior teórico de $\approx 1.1 M_\odot$. Essas estimativas anteriores baseavam-se na determinação da razão de massas binária e da inclinação orbital, principalmente através de medições espectroscópicas de velocidade radial (RV) da estrela doadora e modelagem da modulação elipsoidal óptica. No entanto, SMC X-1 exibe modulações ópticas superorbitais complexas (variações sincronizadas com o ciclo superorbital de raios-X de $\sim 40\text{--}60$ dias) que não foram totalmente contabilizadas em modelos estáticos anteriores. Observações terrestres careciam de precisão fotométrica e amostragem temporal para resolver essas variações dentro de ciclos orbitais individuais, levando a possíveis erros sistemáticos nos parâmetros binários derivados.

Metodologia
Os autores analisaram curvas de luz ópticas de alta cadência de SMC X-1 obtidas pelo Satélite de Levantamento de Exoplanetas em Trânsito (TESS) em seis setores, combinadas com dados de monitoramento de raios-X de longo prazo do Monitor de Imagem de Raios-X de Todo o Céu (MAXI).

1. Análise Observacional: O estudo identificou variações sistemáticas nas curvas de luz orbitais ópticas sincronizadas com a modulação superorbital de raios-X. Especificamente, a profundidade do mínimo na conjunção inferior e a assimetria da curva de luz de duplo pico variavam dependendo do estado de raios-X (subindo, decaindo, ligado ou desligado).
2. Desenvolvimento do Modelo: Um código de modulação elipsoidal modificado foi desenvolvido em Python (utilizando NumPy, SciPy, Astropy e CuPy para paralelização). Ao contrário de modelos estáticos, este código incorpora um disco de acreção precessante que altera:
   • A geometria da irradiação de raios-X na superfície da estrela doadora.
   • A geometria do reprocessamento óptico/UV (OUV) na superfície do disco.
   • O ocultamento mútuo entre a doadora, o disco e o pulsar.
     O modelo calcula a temperatura e a radiação de corpo negro de elementos de superfície (ES) tanto na doadora quanto no disco, contabilizando o escurecimento gravitacional, o aquecimento por raios-X e a irradiação OUV.
3. Estimação de Parâmetros: A amostragem de Cadeia de Markov Monte Carlo (MCMC) foi realizada para ajustar o modelo às curvas de luz ópticas do TESS e de raios-X do MAXI. A análise testou dois cenários para a razão de massas ($q = M_{NS}/M_s$): um baseado em RVs observadas ($\log q = -1.17$) e um segundo onde $q$ foi aumentado em 20% ($\log q = -1.09$) para testar a hipótese de que a RV da doadora observada é subestimada.

Resultados Principais

• Correlações: O estudo confirmou fortes correlações entre o fluxo superorbital de raios-X e características ópticas: o mínimo na conjunção inferior é raso durante estados "desligados" de raios-X e profundo durante estados "ligados"; a assimetria do duplo pico correlaciona-se com a derivada do fluxo de raios-X (por exemplo, o pico esquerdo é mais alto durante a fase de subida de raios-X).
• Validação do Modelo: O modelo elipsoidal modificado reproduziu com sucesso as curvas de luz ópticas e de raios-X observadas, fornecendo evidência observacional de que a modulação superorbital é impulsionada por um disco de acreção precessante. A análise sugere que a precessão do disco é retrógrada.
• Viés de Velocidade Radial: A intensa irradiação de raios-X da estrela doadora desloca o centro de emissão óptica para longe do centro gravitacional. O modelo prevê que esse efeito faz com que a amplitude da velocidade radial observada seja subestimada em aproximadamente 20% (a curva de RV aparente parece "recortada" em comparação com uma onda senoidal verdadeira).
• Estimativa de Massa Revisada:
  • Usando a razão de massas padrão ($\log q = -1.17$), a massa do pulsar derivada é $M_{NS} = 1.11 \pm 0.07 M_\odot$, consistente com estudos anteriores.
  • Corrigindo a subestimação de 20% da RV (aumentando $q$ para $\log q = -1.09$), a massa do pulsar revisada é $M_{NS} = 1.35 \pm 0.10 M_\odot$.
• Discrepâncias: O modelo exigiu uma luminosidade de raios-X ($L_X \approx 4 \times 10^{39}$ erg s$^{-1}$) significativamente maior que os valores observados típicos para reproduzir a intensidade da modulação óptica. Os autores sugerem que isso pode implicar emissão anisotrópica de raios-X (feixes direcionados à doadora) ou desalinhamento spin-órbita, embora este último seja difícil de conciliar com a órbita circular do sistema, a menos que a acreção de massa durante uma fase de envelope comum tenha sido ineficiente.

Significância
O artigo afirma que as variações sistemáticas na curva de luz óptica, impulsionadas pelo disco precessante, não foram previamente modeladas de uma maneira que afetasse a determinação da massa. Ao desenvolver um modelo que contabiliza a geometria de irradiação em mudança, os autores demonstram que a velocidade radial da doadora provavelmente foi subestimada devido a deslocamentos induzidos por irradiação no centro de emissão. A correção desse efeito resulta em uma massa de pulsar revisada de $\approx 1.35 M_\odot$. Este valor coloca o pulsar SMC X-1 bem dentro da faixa típica para estrelas de nêutrons, removendo a tensão com os limites inferiores teóricos para supernovas de colapso do núcleo e sugerindo que estimativas anteriores de massa foram enviesadas por efeitos de irradiação não corrigidos.