Radiation-hydrodynamics of star-disc collisions: From system parameters to outflows and lightcurves

Resumo Técnico: Hidrodinâmica de radiação de colisões estrela–disco: Dos parâmetros do sistema aos fluxos de saída e curvas de luz

Problema e Contexto
As erupções quase periódicas (QPEs) são transientes nucleares associados a buracos negros supermassivos (SMBHs), caracterizados por flares de raios-X suaves recorrentes. Embora existam vários modelos, o cenário de colisão estrela–disco — onde uma estrela em uma órbita apertada, excêntrica e inclinada intersecta repetidamente um disco de acreção — é um candidato principal. Este modelo explica naturalmente os tempos de recorrência, os intervalos alternados "longo–curto" e os padrões de flare "forte–fraco" observados em fontes como GSN 069. No entanto, conectar os parâmetros físicos locais dessas colisões às propriedades observáveis dos flares (luminosidade, duração, assimetria) requer a compreensão da evolução acoplada de choques, geometria de fluxo de saída e transporte de radiação em um meio em expansão e opticamente espesso. Estudos anteriores exploraram partes deste problema usando abordagens 1D, 2D ou não radiativas, mas um estudo sistemático de hidrodinâmica de radiação em 3D dedicado, que ligue os parâmetros do sistema às curvas de luz, era inexistente.

Metodologia
Os autores realizaram um conjunto de simulações de hidrodinâmica de radiação tridimensionais locais utilizando o código de hidrodinâmica de partículas suavizadas (smoothed-particle hydrodynamics - SPH) Phantom. As simulações foram conduzidas em um referencial cartesiano local centrado no plano médio do disco, negligenciando a dinâmica orbital global, a relatividade geral e a autogravidade do gás, uma vez que os tempos de escala das colisões são muito menores do que os tempos de escala orbitais.

O estudo utilizou um modelo fiduciário baseado em uma estrela de sequência principal ($R_\star \approx R_\odot$) colidindo com um disco-$\alpha$ opticamente espesso e dominado pela pressão de radiação em torno de um SMBH de $10^6 M_\odot$. A estrela foi modelada como uma esfera rígida movendo-se a $v_\star = 0,1c$. Os autores variaram sistematicamente cinco parâmetros principais em relação ao modelo fiduciário, variando um por vez:

1. Velocidade estelar ($v_\star$): Variando de $0,5$a$1,5$ vezes o valor fiduciário.
2. Densidade superficial do disco ($\Sigma_d$): Variando de $1$a$100$ vezes o valor fiduciário, abrangendo desde discos finos padrão até discos densos pós-evento de destruição de maré (TDE).
3. Raio estelar ($R_\star$): Variando de $0,5$a$4$ vezes o valor fiduciário, variando efetivamente a razão entre a escala de altura do disco e o raio estelar ($H/R_\star$).
4. Perfil de densidade vertical do disco: Modelado como uma distribuição Gaussiana $\rho(z) \propto e^{-z^2/2\sigma^2}$, variando a concentração de massa em direção ao plano médio ($\sigma/R_\odot$).
5. Ângulo de colisão local ($i$): O ângulo entre o vetor de velocidade relativa e o plano médio do disco, variando de $90^\circ$ (perpendicular) a $30^\circ$ (oblíquo).

O transporte de radiação foi tratado usando a aproximação de difusão de limite de fluxo com opacidade de espalhamento de elétrons. A luminosidade bolométrica emergente foi calculada integrando os fluxos de difusão e advecção através de uma fotosfera determinada dinamicamente.

Principais Contribuições e Resultados
O estudo fornece escalonamentos empíricos para a luminosidade de pico ($L_{\text{peak}}$) e a duração do flare ($\Delta t$) como funções dos parâmetros do sistema. As principais descobertas são:

• Dinâmica de Colisão: A estrela impulsiona um choque de proa (bow shock), criando dois fluxos de saída (frente e trás) em lados opostos do disco. O fluxo de saída frontal é geralmente mais massivo e energético. A morfologia qualitativa do choque de proa e dos fluxos de saída permanece amplamente inalterada através das variações de $v_\star$ e $\Sigma_d$, embora a normalização termodinâmica mude.
• Velocidade Estelar ($v_\star$): Estrelas mais rápidas produzem flares mais brilhantes ($L_{\text{peak}} \propto v_\star^2$) devido a taxas de aquecimento de choque mais altas. Elas também produzem flares mais longos ($\Delta t \propto v_\star^{1,4-2,2}$) porque mais energia é injetada, embora a luminosidade decaia mais rapidamente devido a perdas adiabáticas mais fortes.
• Densidade Superficial do Disco ($\Sigma_d$): Discos mais densos aumentam a massa e o momento dos fluxos de saída, mas têm um efeito fraco na luminosidade de pico ($L_{\text{peak}} \propto \Sigma_d^{-0,1}$ para o fluxo frontal, $\propto \Sigma_d^{0,2}$ para o traseiro). Isso se deve a uma competição entre o aumento do aquecimento de choque e as perdas adiabáticas mais fortes em ejectas opticamente mais espessas. No entanto, discos mais densos prolongam significativamente a duração do flare ($\Delta t \propto \Sigma_d^{0,5-0,6}$) devido aos tempos de difusão de fótons mais longos.
• Raio Estelar ($R_\star$): Estrelas maiores interagem com mais gás, produzindo flares mais brilhantes ($L_{\text{peak}} \propto R_\star^{1,5}$ frontal, $\propto R_\star^{0,7}$ traseiro) e mais longos ($\Delta t \propto R_\star^{1,5}$). O aumento de $R_\star$ intensifica a assimetria de momento entre os fluxos de saída frontal e traseiro.
• Perfil de Densidade Vertical ($\sigma$): Discos mais centralizados (menor $\sigma$) resultam em flares mais brilhantes e curtos. As camadas centrais densas absorvem mais energia de choque, enquanto as camadas externas de baixa densidade permitem que a radiação escape mais eficientemente com perdas adiabáticas reduzidas.
• Ângulo de Colisão ($i$): Travessias mais oblíquas (menor $i$) reduzem a assimetria de momento e luminosidade entre os dois fluxos. O fluxo de saída traseiro torna-se significativamente mais brilhante ($L_{\text{peak}} \propto (\sin i)^{-0,8}$) e os flares tornam-se mais longos ($\Delta t \propto (\sin i)^{-0,9 \text{ a } -1}$) devido ao tempo de interação mais longo e à redistribuição mais simétrica do gás chocado.

Aplicação a GSN 069
Os autores aplicaram seus escalonamentos empíricos à fonte bem caracterizada GSN 069. Ao amostrar o espaço de parâmetros para corresponder às luminosidades de pico observadas ($\sim 3 \times 10^{42}$ erg s$^{-1}$), durações ($\sim 1$ h) e a razão de amplitude forte-fraco, eles identificaram duas regiões viáveis:

1. Solução Preferencial: Uma órbita retrógrada colidindo com um disco pós-TDE denso e verticalmente concentrado ($\Sigma_d \sim 45 \times$ fiduciante, $\sigma/R_\odot \sim 0,3$). Este modelo reproduz bem as propriedades observadas dos flares.
2. Solução Alternativa: Uma órbita retrógrada com densidade de disco fino padrão, mas com uma geometria de colisão muito rasante ($i \sim 6^\circ$). Esta solução é menos precisa e requer parâmetros extremos.

Os autores observam que a preferência por um disco denso e verticalmente concentrado é consistente com a evidência observacional de atividade do tipo TDE em GSN 069.

Significância e Limitações
O artigo afirma fornecer o primeiro estudo sistemático de hidrodinâmica de radiação em 3D ligando parâmetros locais de colisão às curvas de luz de QPEs. Os escalonamentos empíricos derivados servem como uma ponte entre simulações locais e modelos globais de QPE, permitindo a previsão de propriedades de flares baseadas nos parâmetros do sistema. Os autores enfatizam que suas descobertas sugerem que combinações específicas de parâmetros podem reproduzir as amplitudes, durações e ciclos de dever de trabalho características de QPEs.

O estudo reconhece modestamente limitações: os modelos assumem escalonamentos de um parâmetro separáveis, negligenciam a rotação e o cisalhamento do disco, tratam a estrela como uma esfera rígida e não consideram efeitos de ângulo de visão ou evolução espectral (ex: Comptonização). Consequentemente, os resultados são apresentados como regiões fisicamente motivadas do espaço de parâmetros, em vez de modelos únicos e definitivos de melhor ajuste. Os autores sugerem que trabalhos futuros incorporem estruturas estelares realistas, dinâmica de disco e luminosidades dependentes do ângulo de visão para quantificar totalmente a capacidade do modelo de reproduzir a diversidade das fontes observadas de QPE.