Simulation based parameter space for shock in transonic, sub-Keplerian accretion flow onto non-rotating black holes

Resumo Técnico: Espaço de Parâmetros Baseado em Simulação para Choque em Fluxo de Acreção Transônico e Sub-Kepleriano

Definição do Problema
Fluxos de acreção transônicos, sub-keplerianos e não dissipativos sobre buracos negros são governados por dois parâmetros conservados: energia específica ($\epsilon$) e momento angular específico ($l$). Embora modelos analíticos baseados em análise de Rankine-Hugoniot (RH) e equilíbrio vertical tenham identificado um espaço de parâmetros limitado onde choques estacionários podem se formar, simulações numéricas multidimensionais anteriores sugeriram que os choques podem ocorrer fora desses limites teóricos. Essas simulações também indicaram que os choques podem ser dinâmicos (oscilatórios) em vez de estacionários. No entanto, um mapeamento abrangente do espaço de parâmetros de choque baseado em simulações carecia de sistematização. Este trabalho aborda a lacuna entre as previsões analíticas e a realidade numérica ao escanear sistematicamente o espaço de parâmetros $(l, \epsilon)$ para identificar onde ocorre a formação de choque, caracterizar a natureza desses choques (estacionários vs. oscilatórios) e analisar a dinâmica resultante do fluxo.

Metodologia
Os autores realizaram mais de 280 simulações hidrodinâmicas ideais bidimensionais em um domínio cilíndrico $R-Z$, assumindo a axisymmetricidade. As simulações utilizaram um potencial pseudo-newtoniano ($\Phi(r) = -1/2(r-1)$) para modelar um buraco negro não rotativo.

• Esquema Numérico: O código empregou um método de volume finito com precisão de segunda ordem, utilizando um limitador van Leer para reconstrução espacial, um solver de Riemann HLL para fluxos interfaciais e um esquema Runge-Kutta de dois estágios com estabilidade forte preservada (SSP) para integração temporal.
• Condições de Contorno: A entrada foi imposta na fronteira radial externa ($R_{out}$) com variáveis primitivas constantes derivadas de soluções analíticas para pares específicos de $(l, \epsilon)$, garantindo injeção supersônica para evitar o feedback acústico. Condições de saída foram aplicadas na fronteira vertical superior ($Z_{out}$). Condições de reflexão foram usadas no eixo de rotação ($R=0$) e no plano equatorial ($Z=0$). A matéria dentro de um raio $r_{BH} = 1,8$ foi absorvida.
• Varredura de Parâmetros: O estudo escaneou um domínio retangular de $l \in [1,4; 1,9]$ e $\epsilon \in [0,0001; 0,0061]$ usando 286 pontos. Os domínios de simulação foram ajustados para cada ponto com base na localização do ponto sônico do tipo X mais externo. As simulações foram executadas até $t=50000$ para garantir comportamento estacionário ou periódico.

Principais Contribuições e Resultados
A principal contribuição é a identificação de um espaço de parâmetros de formação de choque que é significativamente maior do que a região prevista analiticamente. Os autores categorizaram o espaço de parâmetros em cinco regiões distintas baseadas nas características das soluções:

1. Região A (Choques Estacionários): Esta região contém parâmetros onde se formam choques estacionários. O choque cria uma camada limite suportada por barreira centrífuga (CENBOL). Os autores observam que o limite identificado numericamente para esta região estende-se a valores de $l$ mais baixos do que as previsções analíticas. Essa extensão é atribuída aos movimentos verticais na simulação que acumulam matéria próximo ao equador, aumentando a pressão e facilitando a formação do choque mesmo quando as suposições de equilíbrio vertical em modelos analíticos o proibiam.
2. Região E (Sem Choque): Localizada em valores de $l$ mais baixos, estas soluções exibem uma barreira centrífuga que causa um acúmulo de densidade e um salto supersônico-para-supersônico, mas nenhum choque verdadeiro (transição supersônico-para-subsônico) se forma. O fluxo permanece supersônico durante todo o processo e acreta sobre o buraco negro.
3. Regiões B e C (Choques Oscilatórios): Estas regiões contêm parâmetros onde os choques se formam, mas oscilam.
   • Região B: Caracterizada por oscilações de alta frequência e baixa amplitude, confinadas principalmente à parte inferior da superfície do choque (próximo ao equador). A parte superior do CENBOL permanece estática.
   • Região C: Caracterizada por oscilações de grande amplitude e baixa frequência que envolvem toda a superfície do choque. Estas soluções exibem vórtices turbulentos significativos na região pós-choque.
4. Região D (Dominada por Outflow): Em valores de $l$ elevados, a barreira centrífuga é tão forte que impede a acreção. O fluxo desenvolve um choque de propagação externa que atinge o limite do domínio, deixando para trás um fluxo subsônico. Mais de 86% da matéria injetada é ejetada como outflow.

Dependência de Parâmetros

• Localização do Choque: Para um $\epsilon$ fixo, a localização do choque geralmente se move para mais longe do buraco negro conforme $l$ aumenta. A dependência de $\epsilon$ para um $l$ fixo é não monotônica; em alguns casos, a localização do choque aumenta com $\epsilon$, enquanto em outros (particularmente em $l$ mais baixos), permanece constante ou diminui ligeiramente.
• Outflow: O outflow é gerado em quase todas as simulações devido aos efeitos combinados da barreira centrífuga e dos gradientes de pressão térmica. A razão entre outflow e inflow aumenta significativamente conforme $l$ aumenta, transitando de uma absorção dominante (Regiões A e E) para uma dominância de outflow (Região D).

Significância e Alegações
O artigo afirma que simulações numéricas multidimensionais revelam um espaço de parâmetros para a formação de choque mais amplo do que o compreendido anteriormente através de métodos analíticos. Especificamente, a inclusão da dinâmica vertical permite que choques se formem em momentos angulares mais baixos do que o previsto por modelos de equilíbrio vertical.

Os autores destacam a relevância astrofísica destas descobertas:

• Formação de Coroa: O fluxo pós-choque comprimido e geometricamente espesso (CENBOL) atua como uma região quente e opticamente tênue capaz de realizar o espalhamento Compton inverso de fótons de baixa energia vindos de um disco Kepleriano (ou auto-Comptonização se o disco estiver ausente). Isso fornece um mecanismo natural para a formação da coroa observada em binárias de raios-X e núcleos galácticos ativos em estados espectrais duros.
• Oscilações Quasi-Periódicas (QPOs): A identificação de regiões de choque oscilatório (B e C) apoia a ligação teórica entre oscilações de choque e as QPOs observadas em radiação de alta energia.
• Outflows: As simulações produzem outflows de forma consistente, os quais são frequentemente observados em estados espectrais duros, sugerendo que o mecanismo de choque é um motor unificado tanto para a coroa quanto para o outflow.

Os autores mantêm a modéstia quanto ao seu escopo, observando que os resultados são específicos para buracos negros não rotativos usando um potencial pseudo-newtoniano e hidrodinâmica ideal. Eles reconhecem que incluir viscosidade, perdas radiativas, campos magnéticos ou efeitos relativísticos gerais para buracos negros rotativos poderia alterar os limites do espaço de parâmetros. No entanto, afirmam que o presente trabalho estabelece um mapa fundamental, baseado em simulações, da formação de choque que se estende além dos limites analíticos atuais.