TDEs on FIRE: Illuminating the Cosmic Evolution of Tidal Disruption Rates

Resumo Técnico: TDEs em FIRE

Definição do Problema
Embora os eventos de ruptura de maré (TDEs) tenham sido extensamente caracterizados no universo local, sua prevalência e evolução cósmica em altos redshifts ($z > 1$) permanecem amplamente não restringidas, tanto observacional quanto teoricamente. As amostras observacionais atuais são limitadas pelo fluxo e predominantemente locais, com apenas alguns candidatos jetados detectados em $z \sim 1$. As previsões teóricas enfrentam incertezas significativas devido à falta de perfis estelares resolvidos em galáxias de alto redshift e à baixa restrição da função de massa de buracos negros (BHMF) para buracos negros de baixa massa, que dominam a população de TDEs. Consequentemente, não há uma previsão robusta para a evolução do redshift das taxas de TDE (TDR) ou sua dependência das propriedades das galáxias hospedeiras ao longo do tempo cósmico.

Metodologia
Este estudo utiliza as simulações de zoom cosmológico FIRE-2 (Feedback In Realistic Environments), especificamente o Massive Halo suite (simulações A1, A2, A4 e A8), para computar taxas intrínsecas de TDE na faixa de redshift de $z = 1–10$. A metodologia integra três componentes principais:

1. Estrutura de Simulação: As simulações FIRE-2 empregam o código GIZMO com alta resolução espacial (alcançando escalas de parsec) e modelagem abrangente de formação estelar, evolução e feedback estelar. O Massive Halo suite inclui um modelo sofisticado para a formação e crescimento de buracos negros massivos (MBHs), semeando-os com $M_{seed} = 1,4 \times 10^4 M_\odot$ quando os halos de matéria escura excedem um limiar de massa específico. Os buracos negros são tratados como partículas de sumidouro (sink particles) colisionais com massa dinâmica artificial para garantir que permaneçam no centro da galáxia.
2. Análise de Perfil Estelar: Os autores extraem perfis de densidade estelar radial dentro de 1 kpc dos centros das galáxias (em vez da posição do BH, para considerar potenciais deslocamentos por fricção dinâmica). Esses perfis são ajustados com uma função de lei de potência dupla. Sistemas com inclinações internas não físicas (mais íngremes que as externas) ou ruído numérico são filtrados. A amostra é restrita a galáxias com $M_{gal} > 10^7 M_\odot$ e $M_{BH} < 10^8 M_\odot$ (o limiar acima do qual as estrelas são engolidas inteiras em vez de sofrerem ruptura de maré).
3. Cálculo de Taxa: As taxas de TDE são calculadas usando a teoria do cone de perda (loss cone theory) e o código de acesso público phaseflow. O código resolve a função de distribuição estelar $f(E)$ via inversão de Eddington a partir dos perfis de densidade ajustados. A taxa total de consumo é derivada pela integração do espalhamento de estrelas para dentro do cone de perda devido ao relaxamento de dois corpos, assumindo uma população estelar monocromática do tipo Sol ($m_\star = 1 M_\odot, r_\star = 1 R_\odot$).

Principais Contribuições e Resultados

• Evolução do Redshift da TDR: O estudo fornece a primeira previsão da evolução da TDR através de $z = 1–10$. A TDR média aumenta no universo primitivo, atinge o pico em aproximadamente $4 \times 10^{-4} \text{ yr}^{-1}$ perto de $z \sim 2,5$, e declina para $\sim 10^{-5} \text{ yr}^{-1}$ em $z = 1$. Uma adequação de lei de potência para $z \ge 2,5$ é fornecida: $\log_{10}(\text{TDR}) \approx z^{-0,2} - 4$.
• Correlação com a Formação Estelar: Uma forte correlação é encontrada entre a TDR e a taxa de formação estelar (SFR) da galáxia hospedeira. A TDR acompanha a história da SFR cósmica, atingindo o pico em $z \sim 2,5$. A relação de escala é derivada como $\text{TDR} \propto \text{SFR}^{0,69}$.
• Dependência de Massas ($M_{BH}$ e $M_{gal}$):
  • Massa do Buraco Negro: A TDR aumenta com $M_{BH}$ no regime de buraco negro de massa intermediária (IMBH), atinge o pico em $M_{BH} \sim 1,8 \times 10^5 M_\odot$ e depois declina para buracos negros supermassivos (SMBHs). Essa tendência de "virada" (turnover) coincide qualitativamente com as descobertas no universo local, apesar das diferenças em normalização e inclinações.
  • Massa da Galáxia: Da mesma forma, a TDR aumenta com $M_{gal}$ para galáxias pequenas, atinge o pico em $M_{gal} \sim 5 \times 10^9 M_\odot$ e declina para sistemas mais massivos.
• Direcionadores Estruturais Estelares: A TDR é fortemente impulsionada pelas propriedades estelares centrais. Densidades estelares internas ($\rho_0$) mais altas e inclinações de densidade interna ($\gamma$) mais íngremes correlacionam-se com taxas mais altas. O estudo destaca que as maiores TDRs ocorrem em sistemas com $\rho_0$ alto, $\gamma$ íngreme e pequenos raios de influência de buraco negro ($r_{inf}$).
• Galáxias Satélites: Galáxias satélites (hospedeiras não primárias) exibem TDRs comparáveis às galáxias primárias. Sua contribuição fracionária para a TDR total aumenta significativamente em altos redshifts, sugerindo que elas são cruciais para sondar IMBHs e a montagem precoce de galáxias. A maioria das TDEs de satélites na amostra recai dentro dos limites de resolução angular de instalações de próxima geração, como o Telescópio Espacial Nancy Grace Roman, o Observatório Vera Rubin e o ULTRASAT.

Significância e Alegações
O artigo afirma que simulações cosmológicas, especificamente o FIRE-2, oferecem um caminho promissor para restringir a evolução cósmica das taxas de TDE. O trabalho demonstra que as tendências qualitativas da TDR com relação à massa de buraco negro e de galáxia, observadas no universo local, persistem do alto redshift até o presente dia, sugerindo relações de escala consistentes entre buraco negro e galáxia através do tempo cósmico.

Os autores reconhecem modestamente as limitações, incluindo a ausência de feedback de AGN no Massive Halo suite (que pode aumentar artificialmente as SFRs e densidades centrais em $z \lesssim 2$), a incapacidade de resolver totalmente os aglomerados estelares nucleares devido à resolução finita, e a restrição da simulação a $z=1$ (impedindo restrições em $z=0$). Apesar dessas ressalvas, o estudo estabelece um arcabouço para comparar previsões teóricas com observações futuras de levantamentos de ampla área e alta cadência, visando validar modelos de crescimento de buracos negros, dinâmica estelar e coevolução de galáxias.