Semi-analytic studies of accretion disk and magnetic field geometry in M87*

Resumo Técnico: Estudos semianalíticos do disco de acreção e da geometria do campo magnético em M 87*

Problema e Contexto
O comportamento dos campos magnéticos próximos a buracos negros (BHs) e seu papel nos fluxos de acreção e na formação de jatos permanecem questões centrais na astrofísica de buracos negros. Embora as observações do Event Horizon Telescope (EHT) de M 87* e Sagitário A* tenham fornecido insights sem precedentes sobre estruturas na escala do horizonte, interpretar essas observações requer um arcabouço teórico que ligue a emissão polarizada às propriedades do sistema. Simulações de Magnetohidrodinâmica Relativística Geral (GRMHD) são a ferramenta padrão para isso, mas são computacionalmente dispendiosas e dependem de condições iniciais incertas e suposições de MHD ideal que podem não capturar totalmente a física do plasma colisional. Por outro outro lado, modelos simplificados carecem da consistência física para fazer previsões quantitativas. Este artigo aborda essa lacuna empregando um modelo de Fluxo de Acreção de Baixa Eficiência Radiativa (RIAF) estacionário e semianalítico para investigar sistematicamente como a geometria do campo magnético, o spin do buraco negro, a dinâmica do plasma e a espessura do disco afetam a emissão de sincrotron observável em M 87*.

Metodologia
Os autores construíram um modelo RIAF estacionário e axissimétrico em um espaço-tempo de Kerr usando o código de traçado de raios de código aberto ipole. O modelo assume que a emissão se origina próximo ao plano equatorial, negligenciando a contribuição da base do jato para a morfologia primária do anel.

• Fluxo de Acreção: A densidade numérica de elétrons ($n_e$) e a temperatura ($T_e$) seguem leis de potência do raio esférico. A quatro-velocidade do fluxo é interpolada entre o movimento orbital Kepleriano e a queda livre geodésica, parametrizada por $\kappa_K$ (Kepleriano) e $\kappa_{ff}$ (inflow radial).
• Geometria do Campo Magnético: Sete configurações distintas de campo magnético foram testadas: Radial (monopolo dividido), Vertical, Dipolo, Quadrupolo, Parabólico, Combinado (solução de Wald + radial) e Toroidal. Estas foram definidas via potenciais vetoriais para garantir um campo livre de divergência.
• Parâmetros: O estudo variou o spin adimensional do BH ($a^*$) de -0,94 a +0,94, a espessura geométrica do disco ($H$) e a dinâmica do fluxo. O modelo fiduciário utilizou $H=0,5$ e um perfil de velocidade misto ($\kappa_{ff}=0,5, \kappa_K=0,5$).
• Observáveis Sintéticos: Os autores geraram imagens sintéticas a 230 GHz e extraíram métricas restringidas pelo EHT: diâmetro angular do anel ($\alpha$), assimetria de brilho ($I_r$), polarização fracionária líquida ($|m|_{net}$), polarização fracionária resolvida ($\langle|m|\rangle$) e o segundo modo da decomposição de Fourier da imagem polarizada ($\beta_2$).
• Comparação: Os resultados semianalíticos foram comparados com restrições observacionais do EHT e uma biblioteca de simulações GRMHD de tempo médio (estados MAD e SANE).

Resultados Principais

1. Geometria do Campo Magnético: O estudo confirma que a geometria do campo magnético dita fortemente a estrutura da imagem polarizada. Campos toroidais produzem padrões de Ângulo de Posição do Vetor de Polarização Elétrica (EVPA) radiais, enquanto campos poloidais produzem padrões azimutais. O argumento do parâmetro $\beta_2$ ($\arg \beta_2$) serve como um discriminador robusto: ele está próximo de $0^\circ$ para dominância azimutal (toroidal) e próximo de $180^\circ$ para dominância poloidal. Os dados do EHT, com $\arg \beta_2$ na faixa intermediária ($-163^\circ$ a $-127^\circ$), sugerem uma geometria mista ou dominada por campos poloidais, em vez de uma puramente toroidal.
2. Dinâmica do Plasma e Assimetria: A dinâmica do fluxo impacta significativamente a assimetria de brilho ($I_r$). Modelos Keplerianos (movimento puramente azimutal) produzem alta assimetria impulsionada pelo efeito Doppler, frequentemente excedendo as restrições do EHT. Em contraste, modelos com inflow radial significativo (momento angular zero) são excessivamente simétricos. A assimetria observada em M 87* é melhor reproduzida por modelos com velocidades sub-Keplerianas e spins positivos de moderados a altos.
3. Spin do Buraco Negro: Embora o spin tenha um efeito secundário no diâmetro do anel, ele é um motor primário da assimetria de brilho. Spins positivos ($a^* > 0$) alinham o efeito Doppler e os efeitos relativísticos para aumentar a assimetria, enquanto spins negativos frequentemente resultam em imagens muito simétricas para corresponder às observações. Os dados favorecem valores de spin positivos de moderados a altos.
4. Espessura do Disco: A espessura geométrica do disco ($H$) tem um efeito limitado na maioria dos observáveis, validando o uso de suposições de emissão equatorial para modelos globais dominados pelo disco, exceto em casos específicos como campos dipolares, onde regiões fora do equador contribuem significativamente para as assinaturas polarimétricas.
5. Discrepância de Polarização: Um achado notável é que os modelos RIAF semianalíticos permanecem "Faraday-thin", resultando em alta polarização resolvida ($\langle|m|\rangle$) e altas amplitudes de $\beta_2$ que são inconsistentes com as observações do EHT. Em contraste, simulações GRMHD (particularmente com altas razões de temperatura íon-elétron, $R_h$) exibem despolarização de Faraday significativa. Os autores sugerem que a despolarização observada em M 87* provavelmente surge de estruturas turbulentas não axissimétricas ou do plasma da base do jato fora do equador, efeitos não capturados pela estrutura estacionária semianalítica.

Significância e Conclusões
O artigo demonstra que os observáveis na escala do horizonte dependem sistemática e fisicamente da geometria do campo magnético, da velocidade do plasma e do spin do BH. Os autores concluem que o sistema M 87* é mais consistente com um fluxo dominado por campos magnéticos poloidais com inflow parcialmente radial. Embora o spin do BH permaneça elusivo, os dados desfavorecem spins zero ou levemente negativos, preferindo valores positivos de moderados a altos.

O estudo destaca as limitações dos modelos puramente semianalíticos em reproduzir a complexidade total dos dados do EHT, especificamente no que diz respeito à despolarização de Faraday, o que sugere que características turbulentas, não axissimétricas ou contribuições do jato são críticas. O trabalho serve como uma ponte entre modelos de brinquedo intuitivos e o exaustivo computacionalmente GRMHD, oferecendo um arcabouço paramétrico para interpretar restrições polarimétricas. Os autores sugerem que trabalhos futuros explorem configurações de campo mistas (toroidal-poloidal) e espaços-tempos não-Kerr dentro do arcabouço RIAF para refinar ainda mais essas restrições.