Coupled Time-Dependent Proton Acceleration and Leptonic-Hadronic Radiation in Turbulent Supermassive Black Hole Coronae

Resumo Técnico: Aceleração de Prótons Dependente do Tempo e Radiação Leptônica-Hadrônica Acopladas em Coroas de Buracos Negros Supermassivos Turbulentas

Definição do Problema
Coroas turbulentas que circundam buracos negros supermassivos (SMBHs) são locais potenciais para acelerar partículas não térmicas a altas energias, gerando radiação multi-mensageira observável (neutrinos e ondas eletromagnéticas). No entanto, modelar esses processos é desafiador devido às escalas de tempo comparáveis da aceleração de partículas, resfriamento radiativo e o desenvolvimento de cascata eletromagnéticas (EM). É necessário um tratamento auto-consistente para incorporar o feedback das cascata EM sobre o resfriamento e a aceleração de prótons, particularmente em regiões coronais compactas onde as densidades de fótons são elevadas. Embora existam modelos de estado estacionário, há uma carência de estruturas sistemáticas que acoplem a aceleração de prótons dependente do tempo com o feedback radiativo através de múltiplas escalas espaciais e temporais, especialmente para eventos transientes como eventos de ruptura de maré (TDEs).

Metodologia
Os autores apresentam um arcabouço numérico que acopla auto-consistentemente a aceleração de prótons dependente do tempo com a modelagem de radiação leptônica-hadrônica.

• Aceleração de Prótons: A evolução da função de distribuição de prótons $f(p, t)$ é governada pela equação de Fokker-Planck (FP) dependente do tempo. Esta equação incorpora termos de difusão de momento (aceleração), escape e resfriamento. Os autores empregam o método de Chang-Cooper para discretização espacial e o esquema de Crank-Nicolson para discretização temporal para garantir estabilidade e positividade, mesmo em regimes dominados por resfriamento.
• Modelagem de Radiação: O solver FP é acoplado ao AM3, um software de código aberto para modelagem de radiação lepto-hadrônica dependente do tempo. O AM3 calcula a evolução de partículas secundárias (fótons, neutrinos, $e^\pm$, $\mu^\pm$, $\pi^\pm$, $\pi^0$) resultantes de interações $p\gamma$, $pp$, $\gamma\gamma$, sincrotron de prótons e Bethe-Heitler.
• Mecanismo de Acoplamento: Para evitar a contagem dupla de resfriamento, a equação FP determina exclusivamente a evolução do espectro de prótons. Em cada passo de tempo, o AM3 calcula as taxas de resfriamento ($t_{cool}^{-1} = t_{p\gamma}^{-1} + t_{pp}^{-1} + t_{BH}^{-1} + t_{p-sy}^{-1}$) baseando-se na distribuição de prótons atual e nos campos de fótons. Essas taxas são retroalimentadas no solver FP para atualizar a evolução dos prótons.
• Validação: O código é validado contra soluções de estado estacionário (equações diferenciais ordinárias reduzidas) e testado com funções de injeção e perfis de aceleração/resfriamento específicos para garantir convergência e independência de detalhes de injeção no regime de alta energia.

Principais Aplicações e Resultados
O arcabouço é aplicado a dois cenários astrofísicos distintos:

1. Corona Estacionária em NGC 1068:

   • Configuração: O modelo simula uma corona magnetizada turbulenta ao redor do SMBH na galáxia Seyfert NGC 1068, uma conhecida fonte de neutrinos de alta energia.
   • Resultados: A simulação dependente do tempo reproduz o espectro de neutrinos do IceCube (com pico entre 1–100 TeV) enquanto satisfaz as restrições de raios gama impostas pelo Fermi-LAT. O modelo demonstra que o sistema converge para um estado estacionário dentro de $\sim 300$ escalas de tempo de queda livre ($t_{fs}$).
   • Descobertas: Neste ambiente de alta luminosidade de raios X, o resfriamento de prótons é dominado por interações com fótons de raios X da corona. O feedback das cascata EM sobre a distribuição de prótons é considerado negligenciável, e a solução dependente do tempo alinha-se estreitamente com resultados prévios de estado estacionário.

2. Corona Transiente em TDE AT 2019dsg:

   • Configuração: Um modelo fenomenológico é desenvolvido para uma corona transiente alimentada por um TDE, caracterizada por uma taxa de queda de massa $\propto t^{-5/3}$. O modelo vincula a taxa de acreção dependente do tempo às propriedades da corona (densidade, opacidade).
   • Resultados: A simulação prevê emissões multi-comprimento de onda e de neutrinos atrasadas. Especificamente, o feedback da cascata EM em estágios iniciais impacta significativamente o resfriamento de prótons em coronas mais fracas (onde a luminosidade de raios X é baixa).
   • Descobertas: O acúmulo de cascata EM, combinado com o acúmulo de prótons durante a fase super-Eddington, impulsiona as curvas de luz óptica/UV, raios X e neutrinos para atingir o pico com um atraso de aproximadamente $O(100)$ dias (especificamente $\sim 140$ dias no modelo). Isso oferece uma explicação potencial para a emissão de neutrinos com atraso temporal observada em TDEs.

Significância e Alegações
O artigo afirma fornecer a primeira descrição sistemática acoplando a aceleração de prótons dependente do tempo (via Fokker-Planck) com o feedback radiativo de maneira totalmente dependente do tempo.

• Versatilidade: O arcabouço é apresentado como uma ferramenta flexível capaz de modelar tanto fontes estáticas (ex: NGC 1068, outras galáxias Seyfert emissoras de neutrinos) quanto fenômenos transientes (ex: TDEs, surtos de raios gama, flares de AGN).
• Agnosticismo de Mecanismo: Embora a aplicação principal foque na aceleração turbulenta, os autores observam que o arcabouço pode acomodar outros mecanismos, como reconexão magnética, fluxos de cisalhamento e aceleração de choque, modificando os termos de difusão e deriva.
• Visão Multi-Mensageira: O trabalho une a dinâmica de plasma (microfísica de aceleração) e a modelagem de radiação, oferecendo uma ferramenta para interpretar fenômenos astrofísicos de alta energia no contexto de observatórios futuros como IceCube-Gen2 e KM3NeT.
• Limitações Modestas: Os autores reconhecem que o modelo atual de aceleração de prótons depende de parametrizações existentes (ex: coeficientes de difusão derivados de simulações PIC) e ainda não modela auto-consistentemente o feedback dos prótons acelerados sobre a cascata de turbulência em si, nem a origem detalhada dos raios X da corona. Eles sugerem estes pontos como direções para refinamentos futuros.