Bridging Scales in Black Hole Accretion and Feedback: Subgrid Prescription from First Principles

Resumo Técnico: Prescrição de Subgrade para Acumulação de Buracos Negros a partir de Primeiros Princípios

1. O Problema

A compreensão de como os buracos negros supermassivos (BNS) acoplam-se às suas galáxias hospedeiras através de uma vasta gama de escalas espaciais e temporais permanece um dos principais desafios na evolução galáctica.

• Desafio Computacional: Simulações de Magnetohidrodinâmica Relativística (GRMHD) em escala de horizonte resolvem a física de jatos e acreção, mas geralmente alcançam estados estacionários apenas até algumas centenas de raios gravitacionais ($r_g$). Por outro lado, simulações cosmológicas e de formação de galáxias operam em escalas de kiloparsecs, mas dependem de "prescrições de subgrade" (modelos aproximados) para a acreção e feedback não resolvidos.
• Incerteza Atual: Existe uma lacuna crítica sobre a confiabilidade de extrapolar resultados de simulações de GRMHD de pequena escala (domínio limitado) para escalas galácticas. Modelos atuais frequentemente assumem regimes de equilíbrio que podem não ser alcançados em simulações de curta duração ou ignoram a variabilidade temporal intrínseca.
• Contexto Físico: O foco deste trabalho são os fluxos de acreção "quentes" e radiativamente ineficientes (ADAFs), típicos de buracos negros com baixas razões de Eddington ($f_{Edd} \lesssim 10^{-3}$), onde o feedback é impulsionado principalmente pelo mecanismo de Blandford-Znajek (jatos relativísticos).

2. Metodologia

Os autores utilizam o método "Multizone" (multizona), desenvolvido por seu grupo, para realizar simulações de GRMHD de longa duração que conectam o horizonte de eventos ao raio de Bondi ($R_B$), cobrindo uma faixa dinâmica de ~8 ordens de magnitude.

• Código: Utilização do código KHARMA (Kokkos-based High-Accuracy Relativistic Magnetohydrodynamics with Adaptive mesh refinement).
• Parâmetros Variados:
  • Rotação do Buraco Negro ($a_*$): Variação de $0$a$0.9$ (incluindo um caso de $0.97$).
  • Raio de Bondi ($R_B$): Variação de $4 \times 10^2$ a $2 \times 10^6$ $r_g$.
  • Condições Iniciais: Fluxo de Bondi esférico (sem rotação inicial do gás), mas com testes adicionais para rotação co-rotante, contra-rotante e campos magnéticos fracos.
• Técnica Multizona: O método itera sequencialmente entre passos de tempo pequenos (perto do horizonte) e grandes (longe do horizonte), permitindo que o sistema atinja um equilíbrio dinâmico global em escalas de tempo cosmológicas, algo impossível com métodos convencionais de GRMHD.
• Duração: Simulações de até $700 t_B$ (escalas de tempo de Bondi), o que corresponde a dezenas de milhões de anos físicos para sistemas como M87.

3. Contribuições Principais

O trabalho fornece um conjunto de prescrições de subgrade derivadas de primeiros princípios, prontas para implementação em simulações cosmológicas e modelos semi-analíticos (SAMs). As contribuições dividem-se em duas abordagens:

1. Prescrição Média Temporal (Determinística): Fórmulas analíticas compactas para a taxa de acreção média e potência de feedback, normalizadas à taxa de Bondi clássica.
2. Prescrição Estocástica (Distribuição): Caracterização das distribuições de probabilidade da taxa de acreção e eficiência de feedback, reconhecendo a alta variabilidade temporal.
3. Evolução do Spin: Uma nova relação para a evolução do spin do buraco negro em regimes de acreção quente, demonstrando que o spin é efetivamente "congelado" em baixas taxas de acreção.

4. Resultados Chave

• Taxa de Acreção ($\dot{M}$):

  • A taxa de acreção é fortemente suprimida pelos campos magnéticos em comparação com a previsão clássica de Bondi.
  • A supressão segue uma lei de potência: $\dot{M}/\dot{M}_B \propto (R_B/r_g)^{-0.5}$.
  • Existe uma dependência moderada com o spin: buracos negros com spin alto ($a_* \approx 0.9$) têm taxas de acreção ligeiramente menores (fator de poucos) do que os não rotativos, devido a jatos mais fortes que regulam o influxo.
  • Fórmula Ajustada: $\frac{\dot{M}}{\dot{M}_B} \approx (-3|a_*| + 4.7) \left(\frac{R_B}{r_g}\right)^{-0.5}$.

• Potência de Feedback ($\dot{E}_{fb}$) e Eficiência ($\eta$):

  • A potência de feedback escala como $R_B^{-0.6}$.
  • Para spins altos ($a_* \gtrsim 0.7$), a eficiência de feedback média é constante e independente do raio de Bondi ($\eta \approx 0.3$ para $a_*=0.9$).
  • Para spins baixos, a eficiência depende levemente de $R_B$.
  • Fórmula Ajustada: $\frac{\dot{E}_{fb}}{\dot{M}_B c^2} \approx (0.98|a_*| + 0.13) \left(\frac{R_B}{r_g}\right)^{-0.6}$.
  • Nota-se uma dependência linear com o spin para a potência de feedback, diferindo de modelos teóricos de força livre que preveem dependência quadrática.

• Variabilidade e Estatística:

  • Tanto $\dot{M}$ quanto a eficiência $\eta$ seguem estatísticas lognormais.
  • A largura da distribuição (variabilidade) aumenta com o aumento do raio de Bondi ($R_B$).
  • Para spins altos e grandes $R_B$, a variabilidade pode exceder 2 ordens de magnitude, sugerindo que modelos puramente médios podem ser insuficientes para certas aplicações.

• Evolução do Spin:

  • Contrariando expectativas de simulações de curta duração, os buracos negros sofrem desaceleração (spindown) em todos os casos de acreção quente, independentemente da rotação inicial do gás.
  • O parâmetro de aceleração de spin é $s(a_*) \approx -3.7 a_*$.
  • O tempo de desaceleração é extremamente longo: $t_s \approx 12 (10^{-3}/f_{Edd})$ Gyr.
  • Conclusão: Em regimes de baixa acreção ($f_{Edd} \lesssim 10^{-3}$), o spin do buraco negro permanece efetivamente "congelado" em escalas de tempo cosmológicas. Mudanças significativas de spin ocorrem apenas durante fases de acreção rápida ou fusões.

• Independência das Condições Iniciais:

  • O estado de equilíbrio dinâmico final é robusto e independente da configuração inicial do gás (rotação co-rotante vs. contra-rotante) ou da força do campo magnético inicial. O sistema "esquece" a rotação inicial devido à frenagem magnética e instabilidades de intercâmbio.

5. Significado e Implicações

• Validação de Modelos de Subgrade: O estudo demonstra que extrapolar diretamente resultados de simulações de GRMHD de pequena escala (domínio limitado) para modelos de galáxias inteiras pode levar a erros sistemáticos. Simulações de curta duração não capturam o estado de variabilidade forte e o equilíbrio global alcançado em escalas de tempo galácticas.
• Aplicabilidade Cosmológica: As fórmulas fornecidas permitem que simulações cosmológicas e SAMs incorporem física de acreção e feedback baseada em primeiros princípios, sem necessidade de ajuste fino observacional, melhorando a previsão da evolução de galáxias "vermelhas e mortas".
• Mecanismo de Feedback: A descoberta de que o spin é congelado em baixas taxas de acreção implica que a evolução do spin e a energia dos jatos em galáxias elípticas massivas e aglomerados são governadas por episódios raros de alta acreção, e não por um processo contínuo de desaceleração.
• Limitações: As prescrições são válidas para fluxos de acreção quente e radiativamente ineficientes (ADAF). Em taxas de acreção mais altas (regime de disco de acreção radiativamente eficiente), modelos diferentes seriam necessários.

Em suma, este trabalho estabelece uma ponte física robusta entre a física do horizonte de eventos e a evolução galáctica, fornecendo ferramentas quantitativas essenciais para a próxima geração de simulações cosmológicas.