Evidence for candidate X-ray pulsations from the ultraluminous X-ray source NGC 7456 ULX-1

Resumo Técnico: Evidência de Pulsações de Raios X Candidatas em NGC 7456 ULX-1

Problema e Contexto
Fontes ultraluminosas de raios X (ULXs) são fontes extragalácticas com luminosidades que excedem o limite de Eddington ($>10^{39}$ erg s$^{-1}$). Embora um subconjunto destas, conhecidas como pulsares de ULX (ULXPs), tenham motores centrais de estrelas de nêutrons (NS) confirmados via pulsações coerentes, a vasta maioria dos mais de 1800 candidatos a ULX permanece não classificada. A detecção de pulsações em ULXs distantes é dificultada pelo baixo fluxo, modulação Doppler decorrente do movimento binário e evolução intrínseca rápida do spin devido à acreção super-Eddington. NGC 7456 ULX-1, uma ULX luminosa ($\sim 10^{40}$ erg s$^{-1}$) na galáxia espiral barrada NGC 7456, não apresentou detecção de pulsações anteriormente, deixando sua natureza (IMBH vs. NS acrecionando) incerta. Este estudo visa buscar pulsações candidatas em ULX-1 utilizando técnicas avançadas de temporização capazes de lidar com o desvio de sinal e a modulação orbital.

Metodologia
Os autores analisaram observações de arquivo do XMM-Newton de NGC 7456 de 2005, 2018 e 2023, com foco principal no conjunto de dados de 2023 (exposição de 123,9 ks) devido à sua estatística de fótons superior ($N_\gamma \approx 1,5 \times 10^4$ em 0,2–12,0 keV). A redução de dados seguiu os procedimentos padrão do SAS, utilizando o detector EPIC-pn para análise de temporização devido à sua resolução temporal superior.

A análise empregou uma abordagem de múltiplos estágios:

1. Buscas Cegas e de Aceleração: Espectros de potência iniciais foram gerados usando POWSPC. Para considerar o desvio de frequência, os autores utilizaram o algoritmo ACCELSEARCH (PRESTO) para realizar uma busca de aceleração no domínio de Fourier sobre a faixa de frequência 0,01–6,8 Hz, permitindo um desvio de frequência máximo de $|z| \le 50$ bins.
2. Análise $Z^2_n$ Direcionada: Uma busca refinada utilizando a tarefa HENZSEARCH (HENDRICS) foi conduzida sobre uma banda de frequência estreita (0,221–0,223 Hz) para confirmar candidatos identificados na busca cega.
3. Busca Coerente Desmodulada Orbitalmente: Para abordar a decoerência de fase causada pelo movimento binário, uma busca totalmente coerente foi realizada utilizando um algoritmo de otimização por enxame de partículas (PSO). Esta buscou um espaço de parâmetros de cinco dimensões ($a_1, P_{orb}, T_{asc}, \nu, \dot{\nu}$) para maximizar o estatístico $H$-test (incorporando até três harmônicos) após a aplicação de um modelo de atraso de órbita circular.
4. Validação Estatística: A probabilidade de falso-alarme (FAP) global foi estimada via 100.000 simulações de Monte Carlo de listas de eventos de apenas ruído, preservando os Intervalos de Bom Tempo (GTIs) e as contagens totais de eventos. A Teoria de Valores Extremos (EVT) foi aplicada para modelar a cauda da distribuição além do alcance simulado. Testes de controle foram realizados em uma ULX próxima (ULX-2) e regiões de fundo.

Resultos Principais

• Detecção de Candidato: Um sinal de pulsação candidato foi identificado em uma frequência de $\nu \approx 0,22$ Hz na observação de 2023 do XMM-Newton. O sinal foi detectado com alta significância na busca de aceleração ($p_{search} \sim 4 \times 10^{-3}$) e confirmado via análise $Z^2_n$ direcionada ($Z^2_3 \approx 47,9$).
• Solução Coerente: A busca coerente desmodulada orbitalmente rendeu um estatístico $H$-test máximo de $H_{max} = 64,6$ em $\nu \approx 0,2266$ Hz. Os parâmetros orbitais de melhor ajuste sugerem um semieixo maior projetado de $a_1 \approx 84$ segundos-luz e um período orbital de $P_{orb} \approx 1,40$ dias.
• Evolução do Spin: O candidato exibe uma derivada de frequência de $\dot{\nu} \approx 1,5 \times 10^{-9}$ Hz s$^{-1}$ (equivalente a $\dot{P} \approx -8,2 \times 10^{-8}$ s s$^{-1}$), indicando um rápido aumento de spin (spin-up) impulsionado pelo torque de acreção.
• Dependência de Energia: A fração pulsada aumenta monotonicamente com a energia, subindo de $\sim 0,15$ abaixo de 1 keV para $>0,25$ acima de alguns keV. O perfil de pulso é amplamente sinusoidal, tornando-se mais pontiagudo em energias mais altas, consistente com a emissão de uma coluna de acreção.
• Significância Estatística: A FAP empírica das simulações de Monte Carlo é $< 10^{-5}$. A extrapolação usando o modelo de Valor Extremo Generalizado (GEV) fornece uma FAP global estimada de $\approx 2,5 \times 10^{-6}$. Nenhum sinal significativo foi encontrado nos dados de 2018 (devido à menor estatística) ou em regiões de controle.

Interpretação Física e Significância
Os autores interpretam estes resultados como evidência de que a NGC 7456 ULX-1 abriga uma estrela de nêutrons acrecionando. Com base na taxa de spin-up observada e na suposição de acreção super-Eddington (onde o disco interno reside dentro do raio de esferização), os autores estimam que a força do campo magnético de dipolo superficial seja de $B \sim 10^{12}–10^{14}$ G. Isto coloca o objeto no regime de uma estrela de nêutrons altamente magnetizada passando por uma rápida acreção.

Alegações e Limitações
O artigo enquadra explicitamente estas descobertas como "evidência para um candidato" em vez de uma confirmação definitiva. Os autores observam que a linha de base de observação (124 ks) é comparável a, no máximo, um único ciclo orbital, o que significa que os parâmetros orbitais derivados representam uma solução de correção de fase efetiva em vez de uma efeméride física única. Consequentemente, os autores afirmam que a confirmação independente com conjuntos de dados ou instrumentos adicionais é necessária para verificar o sinal e restringir os parâmetros orbitais definitivamente. O estudo contribui para o crescente corpo de evidências sugerindo que uma fração substancial das ULXs pode abrigar estrelas de nêutrons em vez de buracos negros de massa intermediária.