RABBITS - III. Modelling relativistic accretion discs around spinning black holes in galaxy formation simulations

Título: RABBITS – III. Modelagem de Discos de Acreção Relativísticos ao Redor de Buracos Negros em Rotação em Simulações de Formação de Galáxias

1. O Problema

As simulações hidrodinâmicas de formação de galáxias, embora essenciais para entender a evolução cósmica, enfrentam limitações significativas na modelagem da acreção de gás em Buracos Negros Supermassivos (BNS).

• Subgrid e Estocasticidade: Devido à resolução espacial e temporal limitada, os processos de acreção em pequena escala são tratados por modelos "subgrid" (sub-rede), frequentemente baseados na fórmula de Bondi-Hoyle-Lyttleton. Esses métodos são inerentemente estocásticos e não resolvem a física do disco de acreção.
• Física Newtoniana vs. Relativística: A maioria dos modelos existentes utiliza a prescrição de Shakura & Sunyaev (1973), que é newtoniana e negligencia efeitos relativísticos cruciais, como a dependência do spin do buraco negro na eficiência radiativa e no transporte de momento angular.
• Falta de Conexão com Observações: Sem a estrutura física explícita do disco, é difícil extrair perfis de temperatura e fluxo radiais realistas para gerar Distribuições de Energia Espectral (SEDs) comparáveis diretamente com observações de quasares e AGNs.

2. Metodologia

Os autores desenvolveram e implementaram um modelo de disco de acreção geometricamente fino e relativístico que evolui autoconsistentemente a massa e o vetor de spin do buraco negro. O modelo integra três componentes principais:

1. Modelo de Disco Relativístico: Baseado nas soluções locais de Compère & Oliveri (2017) e Mummery & Balbus (2023). O disco é construído combinando "peças" de soluções analíticas locais que são válidas em diferentes regimes físicos, definidos pela dominância de pressão (gás vs. radiação) e opacidade (espalhamento de elétrons vs. absorção livre-livre).
   • Regimes: O modelo identifica e combina três regimes principais: Gas-ES (Pressão de Gás + Espalhamento de Elétrons), Rad-ES (Pressão de Radiação + Espalhamento de Elétrons) e Gas-FF (Pressão de Gás + Absorção Livre-Livre).
   • Limites: O disco estende-se do raio de fóton ($R_{ph}$) até o raio de estabilidade gravitacional ($R_{outer}$), calculado via o parâmetro de Toomre relativístico. A região interna ao ISCO (Órbita Circular Estável Interna) é tratada com soluções analíticas suaves que conectam a região de queda livre.
2. Integrador Ketju: Utilizado para resolver a dinâmica de pequena escala ao redor dos buracos negros, incluindo termos pós-newtonianos para capturar a interação precisa entre o buraco negro e as partículas estelares/gás.
3. Simulação SPH (Sphgal): O modelo é implementado no código de formação de galáxias Gadget-3 (com o integrador Mstar e hidrodinâmica Sphgal), permitindo simular a evolução de galáxias isoladas com feedback de AGN.

Evolução do Sistema:

• Alinhamento: O modelo calcula explicitamente o momento angular total do disco para determinar se ele está alinhado ou contra-alinhado com o spin do buraco negro (critério de King et al., 2005), influenciando a taxa de acreção e a evolução do spin.
• Feedback: A eficiência radiativa ($\epsilon_r$) é calculada dinamicamente baseada no spin e na orientação do disco, incluindo correções de torque na ISCO. O feedback térmico é injetado no gás circundante com base nessa eficiência.

3. Contribuições Principais

• Primeira Implementação em Simulações de Galáxias: É a primeira vez que um modelo de disco de acreção relativístico completo, que evolui o spin e a estrutura do disco, é acoplado a simulações de formação de galáxias em larga escala.
• Código Público: Os autores lançaram publicamente o modelo em duas versões:
  • Versão C: Para uso como módulo em simulações hidrodinâmicas.
  • Versão Python: Para pós-processamento de simulações existentes e construção de discos de acreção fictícios ("mock discs").
  • Atualização do Relagn: Uma versão atualizada do modelo Relagn capaz de gerar SEDs a partir de parâmetros físicos do disco.
• Redução de Estocasticidade: Ao derivar a taxa de acreção da estrutura física do disco (primeiros princípios) em vez de uma fórmula estocástica de Bondi, o modelo reduz o ruído numérico e fornece taxas de acreção mais físicas e estáveis.
• Previsão de SEDs e Raios: A capacidade de extrair perfis de temperatura e fluxo diretamente da simulação permite a geração de SEDs realistas e a previsão de raios de meia-luz, conectando a física do disco com observações.

4. Resultados

• Estrutura do Disco e Parâmetros: A exploração do espaço de parâmetros (viscosidade, spin, massa e taxa de acreção) mostrou que:
  • O raio externo do disco é sensível à massa do buraco negro e à taxa de acreção, mas menos sensível ao spin.
  • Discos com maior viscosidade suportam mais massa nas bordas antes de se tornarem gravitacionalmente instáveis.
  • A contra-rotação (disco girando no sentido oposto ao spin do BH) altera significativamente o raio da ISCO e a estrutura global do disco, tornando-o mais frio e menos denso.
• Simulação de Galáxia Isolada:
  • O modelo capturou episódios de acreção cíclicos, onde o disco é alimentado por partículas de gás capturadas e depois drena para o buraco negro.
  • Alinhamento e Spin: Simulações com acreção "caótica" (orientação aleatória do disco) resultaram em buracos negros com spins mais altos do que em cenários de acreção "coerente". Isso ocorre porque discos grandes, mesmo com orientação aleatória inicial, tendem a ter momento angular total suficiente para satisfazer o critério de alinhamento, evitando o "spin-down" eficiente que ocorre em discos pequenos contra-alinhados.
  • Feedback: Buracos negros com spins mais altos (no caso caótico) tiveram eficiências radiativas maiores, suprindo mais fortemente o gás ao redor e resultando em um crescimento final de massa ligeiramente menor em comparação com o caso coerente.
• Comparação Observacional:
  • Raios de Meia-Luz: Os raios de meia-luz previstos pelo modelo são sistematicamente maiores (fator de alguns) do que as observações de quasares microlenteados, mas seguem a tendência de escala com a massa do buraco negro ($R \propto M^{0.8}$).
  • SEDs: As SEDs geradas pelo modelo relativístico mostram desvios significativos em relação ao modelo padrão de Novikov & Thorne (1973), especialmente na cauda de baixa frequência (UV/Óptico) e no pico de alta frequência, devido às diferenças na estrutura térmica e na extensão radial do disco.

5. Significância

Este trabalho representa um avanço crucial na astrofísica computacional ao fechar a lacuna entre a física de acreção em pequena escala (relativística) e a evolução de galáxias em grande escala.

• Precisão Física: Permite que simulações de formação de galáxias incorporem a física real dos discos de acreção, melhorando a previsão de taxas de crescimento de buracos negros, evolução de spin e feedback energético.
• Teste de Teorias: Oferece uma ferramenta para testar teorias de acreção e gravidade em ambientes cosmológicos realistas.
• Preparação para Ondas Gravitacionais: Ao modelar a evolução do spin e a acreção de gás versus fusões, o modelo contribui para previsões mais precisas de sinais de ondas gravitacionais de fusões de buracos negros supermassivos.
• Reprodutibilidade: A disponibilização do código permite que a comunidade científica utilize e expanda este modelo, facilitando a comparação direta entre simulações e dados observacionais de quasares e AGNs.