Cyclotron lines in subcritical X-ray pulsars: Monte Carlo simulations reveal the origin of the observed variability

Resumo Técnico: Linhas Ciclotrônicas em Púlsares de Raios X Subcríticos

Declaração do Problema
Púlsares de raios X (XRPs) exibem características de espalhamento ressonante ciclotrônico (CRSFs) em seus espectros, que servem como diagnósticos diretos da intensidade do campo magnético da estrela de nêutrons (NS). Observações revelam variabilidade significativa na energia central da CRSF ($E_{CRSF}$) e na largura da linha ($\sigma_{CRSF}$) em função da luminosidade de raios X ($L$) e da fase do pulso. No regime de luminosidade subcrítica ($L \lesssim 10^{37}$ erg s$^{-1}$), uma correlação positiva entre $E_{CRSF}$ e $L$ é frequentemente observada, contrastando com a correlação negativa vista em fontes supercríticas. Embora modelos qualitativos sugiram que isso surja do desvio Doppler de fótons espalhados por plasma em queda no funil de acreção, um tratamento quantitativo e rigoroso de transferência radiativa dessa variabilidade tem estado ausente. Especificamente, a interação entre espalhamento ressonante, efeitos da pressão de radiação na velocidade do fluxo e geometria de observação requer simulação detalhada para explicar as tendências observadas.

Metodologia
Os autores desenvolveram um código de Monte Carlo (MC) relativístico para realizar cálculos detalhados de transferência radiativa dentro do funil de acreção acima do ponto quente da NS. O modelo adota o quadro físico proposto por Mushtukov et al. (2015b) para fontes subcríticas:

• Geometria: Um canal de acreção cilíndrico semi-infinito alinhado com o eixo magnético ($B = B\hat{z}$).
• Perfis de Plasma: Expressões analíticas para densidade eletrônica ($n_e$) e velocidade de fluxo ($\beta$) são derivadas da conservação de massa e momento, incorporando a desaceleração do plasma em queda devido à pressão de radiação. O perfil de velocidade depende da luminosidade relativa a uma luminosidade crítica $L^*$.
• Transferência Radiativa: Fótons semente são injetados isotropicamente a partir do ponto quente com um espectro de corpo negro. Pacotes de fótons propagam-se através do plasma, sofrendo espalhamento Compton.
• Física de Espalhamento: O código emprega seções de choque de espalhamento magnético ressonante aproximadas e médias em polarização (Loudas et al. 2021) válidas para $b \ll 1$ e $E \ll m_ec^2$. O espalhamento é tratado no referencial de repouso do elétron, com transformações de Lorentz aplicadas para contabilizar o movimento de fluxo do plasma.
• Condições de Contorno: Dois esquemas foram testados: "Reflexão" (reflexão geométrica na base, conservando energia) e "Termalização" (absorção e reemissão). Os resultados principais utilizam o esquema de Reflexão.
• Parâmetros: Simulações foram executadas para uma NS fiducial ($M=1.4 M_\odot$, $R=10^6$ cm, $B \approx 10^{12}$ G) e aplicadas ao caso específico do XRP GX 304−1. A luminosidade de acreção foi variada no intervalo $L/L^* \in (0, 1)$.

Contribuições e Resultados Principais
O estudo fornece a primeira verificação quantitativa via MC do modelo de Mushtukov et al. (2015b), produzindo os seguintes resultados:

1. Morfologia Espectral: Os espectros emergentes exibem uma CRSF proeminente e assimétrica acompanhada por uma asa azul ampla. A CRSF é sistematicamente desviada para o vermelho em relação à energia ciclotrônica clássica local ($E_c$). Este desvio para o vermelho é causado pelo efeito Doppler entre o referencial de laboratório e o referencial do elétron em queda.
2. Dependência da Luminosidade: Tanto $E_{CRSF}$ quanto $\sigma_{CRSF}$ mostram uma correlação positiva com a luminosidade para todos os ângulos de observação. À medida que $L$ aumenta, a pressão de radiação desacelera o fluxo em queda mais efetivamente, reduzindo a velocidade de fluxo ($\beta$). Essa redução na velocidade diminui o desvio Doppler para o vermelho, deslocando $E_{CRSF}$ mais próximo de $E_c$. Simultaneamente, o aumento da profundidade óptica em luminosidades mais altas aprofunda a linha e alarga a asa azul.
3. Dependência do Ângulo de Observação: As propriedades da CRSF variam significativamente com o ângulo de observação $\theta$ (relativo ao eixo magnético).
   • $E_{CRSF}$ diminui à medida que $\cos \theta$ aumenta (ou seja, à medida que a visão se torna mais "frontal"). Isso é devido à dependência do desvio Doppler em relação ao ângulo entre o momento do fóton e a velocidade do fluxo.
   • $\sigma_{CRSF}$ aumenta com $\cos \theta$. Visões frontais amostram uma faixa mais ampla de gradientes de velocidade ao longo da linha de visada, levando a um maior alargamento Doppler.
   • Consequentemente, o modelo prevê uma anticorrelação entre $E_{CRSF}$ e $\sigma_{CRSF}$ durante o ciclo do pulso (variabilidade resolvida em fase), whereas espectros médios em fase mostram uma correlação positiva com a luminosidade.
4. Aplicação ao GX 304−1: Quando aplicado a dados observacionais do GX 304−1, o modelo reproduz com sucesso as correlações positivas observadas de $E_{CRSF}$ e $\sigma_{CRSF}$ com a luminosidade ao longo de quase uma ordem de magnitude. O melhor ajuste é alcançado para ângulos de observação onde o funil de acreção é visto predominantemente de perfil ($\cos \theta \approx 0.0-0.1$), consistente com configurações geométricas específicas (por exemplo, um rotador quase alinhado em alta inclinação ou um rotador quase ortogonal em baixa inclinação).

Significado e Afirmações
O artigo afirma que o espalhamento ressonante da radiação propagando-se através do plasma acrecente no funil magnético fornece uma explicação natural e quantitativa para a variabilidade observada de CRSF em XRPs subcríticos. Os autores enfatizam que seus resultados:

• Validam a hipótese de que o canal de acreção é a região formadora da linha em fontes subcríticas.
• Demonstram que o movimento de fluxo e o desvio Doppler resultante são os principais impulsionadores da variabilidade dependente da luminosidade, sem exigir a formação de um choque sem colisões (CS) como o local primário de formação da linha.
• Destacam o papel crítico da geometria de observação na interpretação de dados resolvidos em fase, prevendo anticorrelações específicas entre energia e largura da linha que diferem das tendências médias em fase.
• Distinguem seu cenário do modelo de CS (onde a formação da linha ocorre entre o choque e a superfície), observando que o cenário CS não preveria a forte variabilidade dependente da fase observada em suas simulações.

Os autores notam limitações, incluindo o uso de geometria cilíndrica simplificada, perfis analíticos 1D e a negligência de efeitos de relatividade geral e polarização. No entanto, argumentam que essas simplificações não alteram as conclusões fundamentais sobre a origem da variabilidade. O trabalho serve como uma prova de conceito para o uso de espalhamento ressonante no funil de acreção para modelar espectros de XRP subcríticos.