Rates of tidal disruption events from constrained cosmological simulations of the local Universe: population properties and implications for transient surveys

Resumo Técnico: Taxas de Eventos de Disrupção de Maré a partir de Simulações Cosmológicas Restritas do Universo Local

Problema
Eventos de disrupção de maré (TDEs) servem como sondas críticas para o estudo de populações de buracos negros supermassivos (SMBHs) quiescentes. No entanto, existe uma discrepância significativa entre as previsões teóricas e as medições observacionais. Modelos idealizados de relaxação de dois corpos preveem taxas de $\gtrsim 10^{-4} \text{ yr}^{-1} \text{ gal}^{-1}$, enquanto restrições observacionais de levantamentos ópticos, de raios X e de infravermelho médio tipicamente resultam em taxas um ordem de magnitude menores ($\sim 10^{-5} \text{ yr}^{-1} \text{ gal}^{-1}$). Além disso, abordagens analíticas tradicionais frequentemente dependem de extrapolações 2D e suposições simplificadas sobre perfis de densidade galáctica (ex: esferas isotérmicas singulares), falhando em considerar a interação complexa da formação de estruturas cosmológicas, histórias de fusão e a diversidade morfológica específica (cusp de densidade ou núcleo/core) das galáxias hospedeiras. À medida que o Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) se aproxima, são necessárias estimativas de taxas de TDE espacialmente resolvidas e refinadas para interpretar os próximos dados de transientes com precisão.

Metodologia
Os autores utilizam o "Simulation of the LOcal Web" (SLOW), uma simulação cosmológica restrita de última geração do Universo local, focando especificamente em seis regiões de alta resolução ("zoom-in"): os superaglomerados Coma, Hercules, Shapley, Virgo e Perseus, e o aglomerado de galáxias Fornax. Estas regiões abrangem extensões radiais de 55–92 Mpc, cobrindo ambientes que variam de escalas de grupos a superaglomerados massivos ($M_{500} \approx 3.5 \times 10^{13} - 1.7 \times 10^{15} M_\odot$).

A metodologia procede através de várias etapas rigorosas de filtragem e cálculo:

1. Extração de Buracos Negros: Os SMBHs são extraídos dos volumes de simulação, com suas massas ($M_\bullet$) e parâmetros de spin adimensional ($\tilde{a}$) registrados.
2. Filtragem de Sequência Principal: Para garantir a validade dos cálculos de dinâmica de cone de perda (loss-cone), um esquema de filtragem estrito isola buracos negros de "sequência principal" dinamicamente estáveis. Isso exclui sistemas não ligados, dominados por rotação ou submetidos a disrupção de maré via quatro critérios:
   • Corte de Deslocamento (Offset Cut): Os BHs devem estar dentro de 1,5 kpc do centro de seu subhalo.
   • Corte de Velocidade: A velocidade relativa deve ser inferior à dispersão de velocidade estelar local.
   • Corte de Anisotropia Fracionária (FA Cut): O ambiente estelar local deve exibir $FA < 0,4$, garantindo uma esfera aproximadamente isotrópica e suportada por pressão, necessária para os formalismos padrão de cone de perda.
   • Corte de Razão de Matéria Escura: É necessária uma razão de massa de matéria escura para massa estelar $> 1$ para excluir núcleos despidos/esvaziados.
3. Classificação Morfológica: As galáxias hospedeiras são classificadas como "cuspy" (com pico de densidade, $g \ge 2$) ou "cored" (com núcleo, $g \le 1$) com base no declive da densidade estelar 3D medido dentro de 1 kpc. Este raio é escolhido para evitar artefatos de suavização numérica, servindo como um proxy para a estrutura nuclear não resolvida.
4. Correção de Eficiência Relativística: O estudo aplica um fator de eficiência relativística ($\eta$) derivado de Kesden (2012). Este fator quantifica a fração do cone de perda que resulta em um TDE observável em vez de uma captura direta, contabilizando o spin do buraco negro e o limite de massa de Hills.
5. Cálculo de Taxa: As taxas volumétricas de TDE ($\Gamma_{vol}$) são calculadas usando o arcabouço de Stone & Metzger (2016), que aplica taxas de linha de base distintas para ambientes de cusp e core, ponderadas pela eficiência relativística $\eta$.

Contribuições Principais e Resultados

• Taxas Volumétricas e por BH: O estudo deriva uma taxa volumétrica média de TDE de $\approx 600 \text{ Gpc}^{-3} \text{ yr}^{-1}$ através de todos os seis ambientes. A taxa média de TDE por buraco negro é encontrada como estável em $\approx 4,5 \times 10^{-5} \text{ yr}^{-1}$, um valor consistente com as restrições observacionais atuais, mas significativamente menor que as previsões teóricas idealizadas.
• Distribuição Espacial vs. Massa Total: O artigo demonstra que os rendimentos de TDE são impulsionados pela demografia dos buracos negros e concentração espacial, em vez da massa total do aglomerado. Enquanto o orçamento absoluto de TDE é esmagadoramente dominado por galáxias satélites cuspy em halos de aglomerados estendidos, a taxa volumétrica ($\Gamma_{vol}$) é fortemente influenciada pela densidade numérica de buracos negros ($n_\bullet$).
• Depleção de Núcleo e Eficiência: As taxas centrais de núcleo são fortemente reduzidas devido à depleção dinâmica (fusões destruindo subhalos de baixa massa) e restrições de captura direta (limite de Hills). Consequentemente, o orçamento total de TDE não é dominado pelas galáxias centrais massivas, mas pela população de satélites de menor massa e tipo cuspy nas periferias dos aglomerados.
• Dependência Ambiental:
  • Hercules: Apesar de ser um superaglomerado massivo e em assembleia ativa, com uma alta densidade numérica de buracos negros, apresenta a menor eficiência por buraco negro ($\tilde{\Gamma}$) devido ao crescimento da massa do buraco negro via fusões, o que empurra muitos BHs em direção ao limite de massa de Hills, onde $\eta \to 0$.
  • Fornax: Este ambiente de baixa massa é o mais eficiente por buraco negro ($\tilde{\Gamma} \approx 1,07$) porque abriga uma fração desproporcional de buracos negros de baixa massa que não sofreram fusão, evitando o corte de Hills e mantendo altas eficiências de cone de perda. Crucialmente, Fornax exibe a menor densidade numérica de buracos negros ($n_\bullet$) e o maior raio de meio número ($R_{1/2,\bullet} = 31,26$ Mpc) entre os seis ambientes, indicando a distribuição espacialmente mais dispersa de buracos negros.
• Discrepância com Extrapolações 2D: O estudo destaca que as extrapolações analíticas 2D tradicionais superestimam as densidades de superfície centrais. As densidades de taxa de superfície 2D simuladas são aproximadamente uma ordem de magnitude menores do que os valores da literatura derivados de projeções 2D, visto que o verdadeiro orçamento de TDE está distribuído pelo halo estendido, em vez de concentrado no núcleo.

Significado
O artigo afirma preencher a lacuna entre as previsões teóricas e as medições observacionais ao utilizar um arcabouço cosmológico 3D completo que resolve os ambientes estelares locais. Ao ir além das extrapolações 2D simplificadas e perfis de densidade idealizados, os autores fornecem um "orçamento de TDE" mais realista que contabiliza a história de montagem complexa do Universo local.

Os resultados sugerem que a taxa observada de TDE de $\sim 3 \times 10^{-5} \text{ yr}^{-1} \text{ gal}^{-1}$ é robustamente sustentada por simulações cosmológicas quando se considera a depleção do núcleo dinâmico e os limites de captura relativística. A taxa derivada de $4,5 \times 10^{-5} \text{ yr}^{-1}$ é apresentada como um limite teórico superior robusto para os próximos levantamentos de transientes, reconhecendo que a física de subgrade da simulação pode produzir levemente em excesso galáxias anãs de baixa massa, que dominam o rendimento absoluto. O estudo enfatiza que as taxas de TDE não são constantes universais, mas são fundamentalmente ditadas pelas propriedades macroscópicas da teia cósmica, especificamente a concentração espacial de buracos negros e a história de fusão de seus ambientes hospedeiros.