Bolometric correction factor and radiative efficiency for the super-Eddington accretion flow in tidal disruption events

Resumo Técnico: Fator de Correção Bolométrica e Eficiência Radiativa para Fluxo de Acreção Super-Eddington em Eventos de Disrupção de Maré

Enunciado do Problema
A estimativa da luminosidade bolométrica e da eficiência radiativa é crítica para a compreensão dos fluxos de acreção ao redor de buracos negros (BHs) supermassivos. Em Eventos de Disrupção de Maré (TDEs), existe uma discrepância significativa entre a massa accrecionada ($\Delta M$) derivada das observações e a expectativa teórica de $\sim 0,5 M_\odot$ para a disrupção de uma estrela do tipo solar. Este problema da "energia ausente" surge porque os cálculos padrão frequentemente assumem um fator de correção bolométrica de raios-X constante ($k_{bol} \approx 1$) e uma alta eficiência radiativa ($\eta = 0,1$). Estudos anteriores sugerem que essas suposições podem ser falhas, particularmente durante a fase inicial super-Eddington, onde a emissão é dominada pela banda extrema ultravioleta (EUV) não observável, em vez de raios-X, e onde a eficiência radiativa é teoricamente esperada como sendo menor devido à advecção.

Metodologia
Para abordar essas incertezas, os autores realizaram simulações de hidrodinâmica de radiação axissimétricas bidimensionais utilizando o código ATHENA++. O estudo modelou o fluxo de acreção resultante de uma estrela do tipo solar ($M_* = M_\odot, R_* = R_\odot$) em uma órbita parabólica disrompida por BHs de massas de $10^6, 10^{6,5}$ e $10^7 M_\odot$.

• Configuração da Simulação: Os autores utilizaram um potencial pseudo-newtoniano para mimetizar efeitos de BH de Schwarzschild. A massa foi injetada continuamente no raio de circularização ($R_c = 2R_t$) com um perfil de taxa de queda (fallback rate) $\dot{M}_{inject} \propto t^{-5/3}$. O domínio computacional abrangeu $2R_S \le r \le 10^5 R_S$ com uma resolução de $768 \times 256$.
• Pós-processamento: Em $t = 16$ dias pós-injeção (um tempo escolhido quando a acreção super-Eddington estável é estabelecida), os autores realizaram cálculos de transferência de radiação de Monte Carlo usando PYTHON para gerar espectros emergentes.
• Parâmetros: Os espectros foram calculados para ângulos de visão ($\theta$) variando de $0^\circ$ (face a face) a $90^\circ$ (de borda). A banda de raios-X foi definida como 0,3–2 keV.
• Cálculos: A luminosidade de raios-X equivalente isotrópica ($L_X$) e a luminosidade bolométrica ($L_{bol}$) foram integradas a partir dos espectros para derivar o fator de correção bolométrica ($k_{bol} = L_{bol}/L_X$) e a eficiência radiativa aparente ($\eta = L_{bol}/\dot{M}_{in}c^2$).

Resultados Principais

1. Dependência de Massa e Ângulo de Visão: Tanto $k_{bol}$ quanto $\eta$ são fortemente dependentes da massa do BH ($M_{BH}$) e do ângulo de visão ($\theta$).
   • $k_{bol}$: Os valores variam de algumas dezenas a vários milhares. Para $M_{BH} = 10^6 M_\odot$, $k_{bol}$ aumenta de $\sim 25$ em $0^\circ$ para $\sim 1000$ em $90^\circ$. Para massas mais altas ($10^7 M_\odot$), $k_{bol}$ permanece relativamente constante em torno de $\sim 1000$ em todos os ângulos. Esses valores são substancialmente maiores do que o valor comumente adotado de unidade.
   • $\eta$: A eficiência radiativa derivada varia de $\sim 10^{-3}$ a $10^{-1}$. Especificamente, para $10^6 M_\odot$, $\eta$ é $\sim 0,001-0,02$; para $10^7 M_\odot$, atinge até $\sim 0,13$, mas geralmente permanece abaixo de 0,1 para a maioria dos ângulos. $\eta$ diminui com o aumento do ângulo de visão e é geralmente menor para massas de BH mais baixas (devido a taxas de acreção escalonadas por Eddington mais altas).
2. Características Espectrais: Os espectros emergentes são dominados pela banda EUV. À medida que $\theta$ aumenta, a emissão de raios-X é cada vez mais obscurecida e reprocessada pelo fluxo (outflow) opticamente espesso para as bandas óptica e UV. Consequentemente, $L_X$ cai por ordens de magnitude (até $\sim 1000\times$ para $10^6 M_\odot$) conforme $\theta$ vai de $0^\circ$ para $90^\circ$, enquanto $L_{bol}$ diminui de forma mais branda ($\sim 15-20\times$).
3. Razão do Fator de Correção: A razão $(k_{bol}/1,1) / (\eta/0,1)$, que escala a massa accrecionada inferida em relação às estimativas anteriores, varia de $\sim 90$ a $8 \times 10^4$ através dos parâmetros simulados. Isso indica que as estimativas anteriores de $\Delta M$ foram provavelmente subestimadas por 2 a 5 ordens de magnitude.

Aplicação a TDEs Observados
Os autores aplicaram seus valores derivados de $k_{bol}$ e $\eta$ a vários TDEs de raios-X específicos (ex: NGC 5905, NGC 3599, SDSS J132341.97+482701.3, ASASSN-14li). Ao selecionar a $M_{BH}$ e $\theta$ apropriados (tipicamente inclinações pequenas consistentes com o modelo unificado de TDEs), a massa accrecionada recalculada $\Delta M$ para esses eventos cai dentro da faixa de $0,25 - 0,5 M_\odot$. Este resultado alinha-se com a expectativa teórica para a disrupção de uma estrela de sequência principal do tipo solar, oferecendo assim uma solução para o problema da "energia ausente" para esses casos específicos.

Significância e Limitações
O artigo afirma que os valores de $k_{bol}$ e $\eta$ derivados fornecem uma correção motivada fisicamente que pode aliviar significativamente o problema da energia ausente em TDEs, desde que o ângulo de visão do observador seja pequeno. Os autores enfatizam que sua solução é consistente com o modelo unificado de TDEs, onde eventos selecionados por raios-X são observados em baixas inclinações.

No entanto, os autores mantêm uma postura modesta quanto à universalidade de suas descobertas:

• Incerteza do Ângulo de Visão: A solução é fortemente dependente de $\theta$, que não pode ser determinado com precisão apenas através de observações.
• Evolução Temporal: Os cálculos baseiam-se em $t=16$ dias. Embora os autores argumentem que a fase super-Eddington é estável por centenas de dias, a evolução de longo prazo de $k_{bol}$ e $\eta$ conforme a taxa de acreção cai abaixo dos limites de Eddington requer mais estudos.
• Simplificações: O estudo assume uma estrela do tipo solar, um BH de Schwarzschild (usando potencial pseudo-newtoniano) e ignora o spin do BH e campos magnéticos. Fatores como o fator de penetração, a estrutura estelar e efeitos relativísticos não foram testados e podem influenciar os resultados.

Em conclusão, este trabalho fornece um conjunto específico de parâmetros derivados de simulações que, quando aplicados a TDEs de raios-X com pequenos ângulos de visão, reconcilia a energia radiada observada com os limites teóricos de disrupção de massa, embora reconheça que uma resolução completa para o problema da energia ausente requer investigação adicional de todo o espaço de parâmetros da física de TDEs.