Interaction between disk and extended corona in a general relativistic framework

Resumo Técnico: Interação entre o Disco e a Coroa Estendida em um Framework de Relatividade Geral (Parte I)

Definição do Problema
Modelos atuais de fluxos de acreção de duas fases em binárias de raios X (XRBs) frequentemente têm dificuldade em explicar as observações de estado duro, onde o disco opticamente espesso se estende próxima ao buraco negro (BH), iluminado por uma coroa de emissão de raios X duros. Estudos anteriores (ex: Haardt & Maraschi 1991; Poutanen et al. 2018) estabeleceram que uma coroa de placa estática acima de um disco possui um limite superior na potência Comptonizada que pode produzir em equilíbrio energético. Se a coroa se tornar excessivamente espessa ou quente, ocorre um efeito cascata descontrolado onde a iluminação aumentada do disco incrementa os fótons semente, impulsionando ainda mais a saída da coroa, impedindo um equilíbrio autoconsistente. Além disso, modelos existentes frequentemente assumem um equilíbrio de energia local, negligenciando variações radiais e efeitos de relatividade geral (GR), ou dependem de códigos que não incorporam o feedback da coroa para o disco. Este trabalho aborda a lacuna no alcance de uma solução de equilíbrio globalmente autoconsistente para um sistema disco-coroa dentro de um framework de GR completo, testando especificamente a viabilidade de uma geometria de placa estática para estados duros.

Metodologia
Os autores atualizaram o código de transferência radiativa de Monte Carlo MONK (Zhang et al. 2019) para incorporar a interação mútua entre o disco de acreção e uma coroa estendida. Os principais avanços metodológicos incluem:

1. Equilíbrio Energético Global Iterativo: O código agora computa iterativamente a Comptonização pela coroa e a iluminação do disco até que um equilíbrio global seja alcançado. A restrição fundamental é que a luminosidade total observada ($L_{tot}$) deve ser igual à potência de acreção total disponível ($L_{acc} = \eta \dot{M}c^2$).
2. Albedo e Mecanismo de Feedback: Um parâmetro de albedo é introduzido para definir a fração do fluxo coronal incidente refletido pelo disco (modelado como uma atmosfera de elétrons semi-infinita usando o formalismo de Chandrasekhar). A fração restante é absorvida, aumentando a temperatura efetiva e a dissipação intrínseca do disco.
3. Parâmetro de Dissipação Viscosa ($\alpha$): Diferente de trabalhos anteriores que assumiam que toda a energia de acreção era dissipada na coroa, este modelo ajusta a fração $\alpha$ da potência viscosa dissipada intrinsecamente no disco. O código itera $\alpha$ de modo que $L_{tot} = L_{acc}$. Se $L_{tot} > L_{acc}$, $\alpha$ é reduzido para diminuir o fluxo de fótons semente e vice-versa.
4. Framework de Relatividade Geral: A propagação de fótons é rastreada via geodésicas nulas em geometria de Kerr. O código contabiliza o desvio da luz (light bending), o redshift gravitacional e as quatro trajetórias possíveis dos fótons: atingir o infinito, atingir a coroa, auto-irradiar o disco ou cair no BH.
5. Análise Espectral: O índice de fóton ($\Gamma$) do componente Comptonizado é derivado ajustando o espectro emergente a uma lei de potência ($N(E) \propto E^{-\Gamma}$) dentro de uma faixa de energia definida pelo overlap de múltiplas ordens de espalhamento.

Resultos Principais
Utilizando parâmetros de referência (massa do BH $10 M_\odot$, spin $a=0$, temperatura da coroa $T_e = 120$ keV, albedo $= 0.5$), o estudo produz os seguintes achados:

• Soluções de Equilíbrio: Soluções de equilíbrio global existem tanto para estados duros quanto intermediários. Para uma corona de placa estática cobrindo totalmente o disco, o índice de fóton mais baixo alcançável é aproximadamente 1.7–1.8 sob condições típicas.
• Limite de Espessura Óptica: Existe um limite superior estrito na espessura óptica da coroa ($\tau_{max} \approx 0.16$ para o caso de referência) além do qual nenhuma solução de equilíbrio física existe (exigindo $\alpha < 0$). À medida que $\tau$ aumenta, a fração de dissipação intrínseca do disco $\alpha$ deve diminuir para manter o equilíbrio energético, atingindo valores próximos de zero onde a emissão do disco é inteiramente reprocessada pela iluminação da coroa.
• Dependências de Parâmetros:
  • Spin do BH: Um spin mais alto ($a=0.998$) leva a espectros mais duros ($\Delta \Gamma \sim -0.15$) devido ao aumento do desvio da luz e da iluminação do disco, o que força um $\alpha$ menor.
  • Temperatura da Coroa: Uma $T_e$ mais alta resulta em espectros mais duros.
  • Albedo: Um albedo mais alto endurece significativamente os espectros. Para um albedo de 1.0 (reflexão completa), o sistema pode sustentar espessuras ópticas muito maiores ($\tau \approx 0.6$), permitindo que o índice de fóton caia para $\Gamma \sim 1.0$. Isso ocorre porque o alto albedo suprime os fótons semente suaves reprocessados, desacoplando o loop de feedback.
• Forma Espectral: O espectro Comptonizado frequentemente desvia de uma simples lei de potência, exibindo uma "quebra de anisotropia" em altas energias, uma característica atribuída à injeção unidirecional de fótons semente e aos efeitos de GR.
• Fluxo de Energia Radial: Sob a suposição de uma corona de placa estendida e espacialmente uniforme, alcançar o equilíbrio local em cada anel radial requer que a energia flua dos anéis internos para os externos, além de satisfazer o equilíbrio global.

Significância e Conclusões
Os autores concluem que, com o código MONK atualizado, é possível alcançar o equilíbrio de energia global para sistemas disco-coroa, fornecendo um framework robusto para restringir os parâmetros do sistema. No entanto, o estudo constata que uma corona de placa estática não é uma geometria favorável para explicar os espectros mais duros observados em XRBs, mesmo quando o equilíbrio de energia global é imposto. A incapacidade da placa estática de produzir índices de fóton significativamente abaixo de 1.7–1.8 (a menos que o albedo seja unidade) sugere que geometrias alternativas são necessárias. Os autores afirmam que trabalhos futuros explorarão geometrias de disco truncado e coronas de fluxo de saída (outflowing) como potenciais soluções para essas limitações. A abordagem também possui potencial para restringir a geometria da coroa usando dados de polarização, além do fluxo.