Proton Acceleration by Collisionless Shocks in Supermassive Black Hole Coronae: Implications for High-Energy Neutrinos

Resumo Técnico: Aceleração de Prótons por Choques Sem Colisões em Coronas de Buracos Negros Supermassivos

Problema
Observações recentes do Observatório de Neutrinos IceCube detectaram um excesso significativo de neutrinos de alta energia provenientes de núcleos galácticos ativos (AGNs) radio-silenciosos próximos, notadamente a galáxia Seyfert NGC 1068. Estas observações apresentam uma discrepância: embora neutrinos de alta energia sejam detectados, o fluxo correspondente de raios gama de alta energia é significativamente menor do que o previsto pelos modelos hadrônicos padrão. Isso sugere que os neutrinos são produzidos em uma região opaca aos raios gama, provavelmente na coroa quente que circunda o buraco negro supermassivo (SMBH) central, onde os raios gama são atenuados via produção de pares $\gamma\gamma$ com fótons de raios X.

Para explicar o espectro de neutrinos observado, prótons parentais devem ser acelerados a energias de $\sim 100$ TeV dentro desta coroa. No entanto, o mecanismo de aceleração dominante na coroa de SMBHs permanece não estabelecido. Embora a Aceleração por Choque Difusiva (DSA) seja um candidato principal, sua viabilidade é incerta sob condições de coroa. Especificamente, choques de coroa são esperados ter números de Mach sônicos baixos ($M_s \sim 2-10$) e velocidades trans-relativísticas, regimes tradicionalmente considerados ineficientes para a aceleração de partículas. Além disso, o comportamento de plasmas de duas temperaturas ($T_i \gg T_e$) e a eficiência específica da aceleração de prótons versus elétrons nesses ambientes carecem de restrições baseadas em primeiros princípios.

Metodologia
Para abordar essas incertezas, os autores empregam simulações de partícula-em-célula (PIC) unidimensionais de primeiros princípios usando o código SMILEI. O estudo foca em choques quase-paralelos ($\theta = 30^\circ$) ajustados para os parâmetros físicos de coroas de SMBH.

O levantamento de simulações explora sistematicamente um amplo espaço de parâmetros motivado por restrições observacionais para galáxias Seyfert:

• Número de Mach Sônico ($M_s$): Variando de trans-sônico ($M_s \approx 2$) a valores mais altos ($M_s \approx 30$).
• Número de Mach Alfvénico ($M_A$): Variando de $4$a$46$.
• Velocidade do Choque ($v_{sh}$): Velocidades trans-relativísticas ($v_{sh}/c \approx 0,17 - 0,66$).
• Razão de Temperatura Íon-Elétron ($T_i/T_e$): Variando de $0,01$a$50$, refletindo os estados térmicos desacoplados esperados em coroas colisionais.

As simulações utilizam uma razão de massa íon-elétron reduzida ($m_i/m_e = 100$) para garantir a viabilidade computacional enquanto preservam parâmetros adimensionais fundamentais. As execuções foram evoluídas por tempo suficiente ($t \gtrsim 400 \Omega_{ci}^{-1}$) para permitir a saturação das instabilidades eletromagnéticas e o desenvolvimento de caudas de potência não térmicas.

Resultados Principais

1. Aceleração de Prótons Robusta: O estudo descobre que a DSA é um mecanismo robusto e eficiente para a aceleração de prótons em coroas de SMBH. Apesar dos baixos números de Mach sônicos, os choques consistentemente canalizam aproximadamente 10% da energia cinética do choque para íons não térmicos. Esta eficiência permanece estável em todo o levantamento de parâmetros, incluindo casos com $M_s \approx 2$, desafiando a noção de que choques de baixo $M_s$ são aceleradores ineficazes.
2. Dinâmica de Aceleração: A energia máxima dos prótons acelerados aumenta linearmente com o tempo após a estrutura do choque atingir um estado de quase-equilíbrio. A taxa de aceleração segue a escala prevista pela teoria DSA, $t_{acc} \propto \gamma \beta_{sh}^{-2} \Omega_{ci}^{-1}$. Os autores derivam um coeficiente de aceleração normalizado, $C^{-1} \equiv d(\gamma \beta_{sh}^{-2})/d(t\Omega_{ci})$, que exibe um valor universal de $C^{-1} \approx 0,1$, variando por menos de um fator de dois em todos os parâmetros testados.
3. Ineficiência da Aceleração de Elétrons: Em contraste com os prótons, a aceleração de elétrons é altamente ineficiente, com frações de energia de elétrons não térmicos permanecendo abaixo de 1% (tipicamente $\sim 10^{-3}$ para razões $T_i/T_e$ do tipo coronal). A eficiência da aceleração de elétrons diminui conforme a razão $T_i/T_e$ aumenta. Em casos trans-relativísticos específicos com amplificação de campo magnético forte, a cauda de potência de elétrons falha em se estender para energias mais altas ao longo do tempo, sugerindo que a DSA sustentada é inibida quando campos amplificados tornam o choque efetivamente superluminal para os elétrons.
4. Mecanismos de Instabilidade: A turbulência magnética upstream necessária para o espalhamento é impulsionada principalmente pela instabilidade não ressonante (NRI, ou instabilidade de Bell) na maior parte do espaço de parâmetros. A instabilidade ressonante (RI) torna-se dominante apenas em números de Mach alfvênicos muito baixos ($M_A \lesssim 8$).

Implicações para NGC 1068 e Produção de Neutrinos
Aplicando estes resultados a NGC 1068, os autores demonstram que:

• Escalas de Tempo: Prótons podem ser acelerados a $\sim 100$ TeV em escalas de tempo menores que as escalas de resfriamento relevantes (Bethe-Heitler, fotomessônica, $pp$) e dinâmicas (queda livre).
• Orçamento de Energia: A potência não térmica de prótons necessária para explicar a luminosidade de neutrinos observada é estimada em $\sim 3\%$ da potência total do choque. Isto está bem dentro dos $\sim 10\%$ de eficiência de aceleração medidos nas simulações.
• Déficit de Raios Gama: A baixa eficiência da aceleração de elétrons e a alta opacidade da coroa aos raios gama explicam naturalmente o déficit observado de raios gama de alta energia, permitindo simultaneamente que os neutrinos escapem.

Significância e Alegações
O artigo afirma fornecer as primeiras restrições de primeiros princípios e autossistentes sobre a DSA no contexto específico de coroas de SMBH. A principal significância reside em estabelecer que choques sem colisões são um mecanismo viável e eficiente para produzir os prótons de alta energia necessários para explicar o sinal de neutrinos do IceCube de NGC 1068, mesmo sob as condições desafiadoras de baixos números de Mach sônicos e altas razões de temperatura íon-elétron.

Os autores declaram explicitamente que, embora seus resultados apoiem fortemente o modelo hadrônico para a produção de neutrinos em coroas de AGN, a aplicação de resultados de simulações cinéticas locais a um objeto astrofísico complexo como uma coroa de AGN requer consideração cuidadosa. A análise do orçamento de energia é apresentada como uma avaliação de ordem de magnitude baseada na microfísica estabelecida no artigo, fazendo a ponte entre a física de plasma microscópica e os fenômenos macroscópicos de alta energia. O estudo não pretende resolver todas as incertezas relativas à energia máxima do próton ou a transição detalhada das instabilidades em 3D, observando que estas requerem investigação adicional.