Black Hole Spin-down in Collapsars in 3D Neutrino Transport GRMHD Simulations

Título: Evolução do Spin de Buracos Negros em Colapsares em Simulações 3D de GRMHD com Transporte de Neutrinos

1. O Problema

Os collapsars (estrelas massivas cujos núcleos colapsam diretamente em buracos negros) são os principais candidatos a motores centrais de Gamma-Ray Bursts de longa duração (LGRBs). A potência desses jatos relativísticos depende criticamente da força do campo magnético e do spin (rotação) do buraco negro (BH).

• Desafio Anterior: Trabalhos anteriores com simulações de discos de acreção adiabáticos (não radiativos) mostraram que os jatos magnéticos (estado MAD - Magnetically Arrested Disk) tendem a frear rapidamente o spin do BH para valores de equilíbrio muito baixos ($a_{eq} \approx 0.035 - 0.07$).
• Conflito: Esses spins baixos resultam em jatos de potência extremamente reduzida, que podem não conseguir escapar da estrela progenitora. Além disso, observa-se uma discrepância entre esses valores teóricos e as distribuições de spin inferidas por ondas gravitacionais (LIGO/Virgo/KAGRA), que sugerem spins ligeiramente mais altos, embora ainda baixos.
• Fator Faltante: Em collapsars, os discos de acreção são densos e quentes, emitindo neutrinos que resfriam o disco. O resfriamento por neutrinos reduz a espessura do disco e altera o transporte de momento angular, mas este efeito não havia sido incluído em simulações de evolução de spin de collapsars devido ao alto custo computacional do transporte de neutrinos.

2. Metodologia

Os autores realizaram uma série de simulações 3D de Magnetohidrodinâmica Relativística Geral (GRMHD) com transporte de neutrinos (esquema M1 de duas momentos, "cinza").

• Código: Utilizaram o código νH-AMR, acelerado por GPU, que implementa o transporte de neutrinos para três espécies ($\nu_e, \bar{\nu}_e, \nu_x$).
• Configuração Inicial: Simularam dois tipos de perfis de densidade de estrelas progenitoras:
  1. Modelos p#: Perfis de lei de potência ($\rho \propto r^{-1.5}$), representando collapsars típicos com taxas de acreção moderadas ($\dot{m} \sim 0.1-1 M_\odot/s$).
  2. Modelos c#: Núcleos de densidade constante seguidos por perfis mais íngremes, resultando em taxas de acreção muito altas no início ($\dot{m} \gtrsim 10 M_\odot/s$).
• Variação de Parâmetros: Foram testados diferentes spins iniciais do BH ($a = 0.0$ a $0.8$) e comparados casos com e sem transporte de neutrinos.
• Análise: Mediram fluxos de massa, momento angular e energia no horizonte de eventos para calcular o parâmetro de "spin-up" ($s$) e determinar o spin de equilíbrio ($a_{eq}$).

3. Contribuições Principais

• Primeira Simulação 3D: Esta é a primeira vez que a evolução do spin de BHs em collapsars é simulada com transporte de neutrinos completo em 3D.
• Modelo de Evolução de Spin: Desenvolveram um modelo semi-analítico que integra os fluxos hidrodinâmicos, eletromagnéticos e de neutrinos para prever a evolução temporal do spin do BH em escalas de tempo de GRB (dezenas de segundos).
• Quantificação do Efeito de Resfriamento: Demonstraram quantitativamente como o resfriamento por neutrinos altera a termodinâmica do disco e, consequentemente, a eficiência do freio magnético no BH.

4. Resultados Chave

• Spin de Equilíbrio Mais Alto: Ao contrário dos discos adiabáticos, os discos resfriados por neutrinos atingem um spin de equilíbrio significativamente maior: $a_{eq} \approx 0.13$.
  • Isso é 2 a 4 vezes maior que os valores encontrados em discos não radiativos ($0.035 - 0.07$).
  • É consistente com as distribuições de spin inferidas por observações de ondas gravitacionais.
• Potência dos Jatos: O aumento do spin de equilíbrio resulta em jatos 4 a 16 vezes mais potentes do que os previstos por modelos não radiativos. Isso torna os jatos muito mais capazes de perfurar e escapar da estrela progenitora.
• Independência do Perfil Inicial: O valor de $a_{eq} \approx 0.13$ é robusto e mantém-se consistente em uma ampla faixa de estruturas de progenitores e taxas de acreção ($\dot{m} \sim 0.1 - 10 M_\odot/s$).
• Papel dos Neutrinos:
  • A contribuição direta dos neutrinos para o torque de spin é desprezível.
  • No entanto, o resfriamento por neutrinos é crucial: ele reduz a espessura do disco ($h/r$), o que modifica a eficiência do acoplamento entre o disco e o BH, permitindo que o spin se mantenha mais alto antes de atingir o equilíbrio.
• Dinâmica de Baixo Spin: Em BHs com spin inicial muito baixo ($a \approx 0.1$), os jatos podem se tornar instáveis (instabilidade de kink), causando desalinhamento do disco e variabilidade na potência do jato, especialmente em progenitores com perfis de densidade mais rasos.
• Evolução Temporal: Para durações típicas de GRB ($t \gtrsim 30$ s), o BH tem tempo suficiente para acrecionar massa suficiente e atingir o estado de equilíbrio $a_{eq} \approx 0.13$, independentemente do spin inicial (desde que $0 < a_0 < 1$).

5. Significado e Implicações

• Resolução de Tensões Observacionais: O valor de $a_{eq} \approx 0.13$ reconcilia melhor as simulações teóricas com as observações de ondas gravitacionais de fusões de buracos negros binários, sugerindo que collapsars podem ser progenitores de uma fração significativa desses eventos.
• Limites de Luminosidade: O freio rápido do spin (dentro de ~5 segundos) limita a potência máxima dos jatos. Isso pode explicar a existência de um "corte superior" na luminosidade observada de GRBs, impedindo que jatos excessivamente energéticos (que ocorreriam com spins iniciais máximos) sejam produzidos.
• Viabilidade de Jatos: Os jatos mais potentes resultantes de spins ligeiramente mais altos ($a \approx 0.13$) têm maior probabilidade de escapar da estrela, resolvendo o problema de jatos "fracos" que falham em sair da estrela em modelos anteriores.
• Futuro: O trabalho estabelece uma nova base para modelagem de collapsars, destacando a necessidade de incluir transporte de neutrinos para previsões precisas de energia de GRBs e evolução de spin de buracos negros.

Em resumo, o estudo demonstra que o resfriamento por neutrinos é um fator físico crítico que eleva o spin de equilíbrio dos buracos negros em collapsars, tornando-os motores mais eficientes e consistentes com as observações astrofísicas atuais.