Deep optical spectroscopic monitoring of the pulsating ULX NGC 1313 X-2 with longslit Gemini observations

Resumo Técnico: Monitoramento espectroscópico óptico profundo da ULX pulsante NGC 1313 X-2 com observações de fenda longa do Gemini

Problema e Contexto
Fontes ultraluminosas de raios X (ULXs) são hipotetizadas como sistemas binários de alta massa (HMXBs) extremos que excedem o limite de Eddington. Embora as pulsações em várias ULXs tenham confirmado a presença de estrelas de nêutrons, a natureza de suas estrelas doadoras e os parâmetros orbitais da binária permanecem mally restritos para muitos sistemas. NGC 1313 X-2 é uma ULX pulsante conhecida, localizada nos arredores da galáxia espiral barrada NGC 1313. Estudos anteriores debateram sua história de formação (ex: interações de maré vs. formação estelar in-situ) e o tipo espectral de sua companheira, com sugestões variando de tipos O/B a supergigantes vermelhas, frequentemente baseadas em ajustes de distribuição de energia espectral (SED) de banda larga ou dados ópticos limitados. O principal desafio abordado neste trabalho é a identificação do tipo espectral da estrela companheira e a subsequente derivação de restrições orbitais usando espectroscopia óptica de alta qualidade e profundidade.

Metodologia
Os autores utilizaram dados espectroscópicos de fenda longa arquivados do observatório Gemini-South (GMOS-S), coletados ao longo de dez noites entre outubro e dezembro de 2009.

• Redução de Dados: Dez quadros de ciência por noite foram reduzidos usando rotinas pyraf. A redução incluiu subtração de viés (bias subtraction), correção de campo plano (flat-fielding), calibração de comprimento de onda (usando lâmpadas de arco CuAr) e subtração de fundo de céu. Para mitigar as lacunas do detector, foi empregado um padrão de dithering de três posições.
• Empilhamento: Espectros calibrados em fluxo de todas as dez noites foram co-adicionados usando empilhamento de mediana com clipping de $3\sigma$ para excluir raios cósmicos, resultando em um espectro 2D e 1D de alta razão sinal-ruído.
• Análise Espectral:
  • Linhas de Emissão: Linhas de emissão características ([O III] $\lambda$4959, [O III] $\lambda$5007 e [He II] $\lambda$4686) foram modeladas usando uma combinação de funções Lorentziana e Gaussiana para determinar comprimentos de onda centrais, largura total à meia altura (FWHM) e dispersões de velocidade.
  • Ajuste de Continum: O continuum foi ajustado utilizando modelos estelares do atlas Castelli-Kurucz. O processo de ajuste permitiu que o raio estelar, gravidade superficial ($\log g$), temperatura efetiva ($T_{\text{eff}}$) e metalicidade ([M/H]) variassem. A extinção foi corrigida usando $E(B-V) = 0,13 \pm 0,01$ derivado de espectroscopia MUSE.
  • Análise de Variabilidade: Os autores realizaram a correlação cruzada dos dados ópticos com dados de raios X do Swift-XRT quase simultâneos e fluxos de [He II] para investigar correlações e buscar atrasos de tempo (time lags) usando espectros de atraso-frequência.
  • Modelagem Orbital: Assumindo que a companheira preenche seu lobo de Roche, os autores restringiram a separação orbital, a massa da companheira e a excentricidade usando o raio estelar derivado e a fórmula de Eggleton (1983).
• Identificação da Companheira: O espectro empilhado revela uma queda distinta de fluxo abaixo de $\sim$4000 Å. Ao ajustar o espectro na faixa de 3300–6500 Å (incluindo esta queda), os melhores modelos estelares apontam consistentemente para uma supergigante do tipo A ($T_{\text{eff}} \approx 7000\text{--}8250$ K, $\log g \approx +0,5$, $R_* \approx 45\text{--}75 R_\odot$). Os autores observam que, se os comprimentos de onda abaixo de 4000 Å forem excluídos, companheiras do tipo B tardio também são estatisticamente possíveis, mas a presença da quebra de Balmer favorece fortemente a classificação do tipo A.
• Restrições Orbitais: Assumindo que a supergigante do tipo A preenche seu lobo de Roche, o estudo deriva restrições para o sistema binário. As separações orbitais são restritas ao intervalo $50 R_\odot < a < 200 R_\odot$. A análise exclui soluções onde o raio estelar é significativamente menor que 0,9 $R_L$ ou maior que 1,1 $R_L$.
• Cinemática de Linhas de Emissão: As linhas de [O III] mostram velocidades e dispersões constantes, consistentes com a nebulosa de bolha circundante. Em contraste, a linha [He II] $\lambda$4686 exibe variabilidade significativa tanto em fluxo quanto em velocidade, apoiando sua origem nas proximidades imediatas da ULX (provavelmente o fluxo de acreção). Nenhuma periodicidade foi detectada na linha [He II] ou no continuum óptico, impedando uma medição direta do período orbital via curvas de velocidade radial.
• Disco de Acreção e Variabilidade: O espectro de atraso-frequência entre os dados de raios X e ópticos mostra alguma estrutura sugestiva de atrasos de propagação em escalas de tempo de dezenas de dias. No entanto, devido à amostragem esparsa e aos requisitos de interpolação, atrasos de tempo precisos não puderam ser robustamente restringidos. Os autores observam que um argumento simples de tempo de viagem da luz para o atraso observado implicaria um tamanho de disco não fisicamente possível ($\sim 3 \times 10^5 R_\odot$), sugerindo que o atraso pode ser impulsionado por processos viscosos ou que a estimativa é pouco confiável.
• Excesso Azul: Um excesso de fluxo para o azul de 3500 Å permanece inexplicado apenas pelo modelo da supergigante do tipo A. Os autores sugerem que isso poderia resultar de um disco de acreção irradiado, reprocessamento em um vento opticamente espesso ou um disco circumbinário, mas não identificam a fonte definitivamente.

Significância e Alegações
O artigo afirma fornecer a primeira identificação tentativa do tipo estelar para a companheira de NGC 1313 X-2 como uma supergigante do tipo A, baseada na detecção de uma quebra de Balmer em espectroscopia profunda do Gemini. Esta descoberta atualiza as restrições anteriores sobre a natureza binária e os parâmetros orbitais do sistema.

Os autores argumentam que a presença de uma supergigante do tipo A (idade $\lesssim$ 1 Myr) é consistente com simulações de síntese de população para jovens ULXs contendo estrelas de nêutrons. Além disso, a localização da ULX nos arredores da galáxia, combinada com a idade jovem da nebulosa de bolha circundante, sugere que o sistema provavelmente se formou via formação estelar in-situ desencadeada pela acreção de gás do meio intergaláctico, em vez de através de recuo gravitacional do centro galáctico. O estudo reforça a visão de que NGC 1313 X-2 é uma HMXB extrema padrão, consistente com a ausência de um corte na função de luminosidade de raios X das ULXs.

O trabalho também corrobora achados anteriores sobre as taxas de expansão da bolha de gás e o poder cinético do vento/jato em relação ao poder de acreção, utilizando perfis de linhas de [O III]. Embora os autores reconheçam as limitações em restringir o tamanho do disco de acreção e a origem exata do excesso UV/azul, as restrições espectroscópicas sobre a estrela doadora representam um passo significativo na caracterização da história de formação e do estado evolutivo desta ULX pulsante.