Dynamics and detectability of long-lived non-accretion phases for massive black hole binaries in cold, thermally regulating disks

Resumo Técnico: Dinâmica e Detectabilidade de Fases de Não-Acreção de Longa Duração para Binários de Buracos Negros Massivos

Enunciado do Problema
Simulações hidrodinâmicas anteriores de binários de buracos negros supermassivos (SMBHBs) interagindo com discos circumbinários (CBDs) basearam-se amplamente em suposições localmente isotérmicas ou discos com razões de aspecto relativamente altas ($h/r \sim 0,1$). Esses estudos sugeriram que, embora o potencial binário tramite o disco interno, o gás ainda acreta via fluxos de maré. No entanto, SMBHBs realistas devem possuir discos significativamente mais finos ($h/r \sim 10^{-2} - 10^{-3}$), onde a redução do suporte de pressão pode funcionalmente cessar a acreção. Uma questão aberta crítica é se as fases de "não-acreção em fuga" observadas em simulações de discos isotérmicos e finos sobrevivem quando a termodinâmica do disco é evoluída de forma autoconsistente. Especificamente, não está claro se o aquecimento viscoso e hidrodinâmico no disco interno pode aquecer o gás o suficiente para restaurar o suporte de pressão e permitir que a acreção retome, ou se o sistema permanece em um estado de acreção suprimida de longa duração.

Metodologia
Os autores apresentam simulações hidrodinâmicas de alta resolução baseadas em grade, utilizando o código de volume finito Sailfish. O estudo resolve as equações de Navier-Stokes compressíveis 2D, verticalmente integradas, para um disco de gás em torno de um binário de massas iguais. As principais características metodológicas incluem:

• Equação de Energia: Diferente dos estudos isotérmicos anteriores, os autores evoluem uma equação de energia que rastreia de forma autoconsistente o aquecimento viscoso e hidrodinâmico, acoplado ao resfriamento de corpo negro radiativo. Assume-se que o disco é opticamente espesso com opacidade dominada por espalhamento de elétrons.
• Configuração do Sistema: As simulações modelam um binário com massa total $M = 10^6 M_\odot$ e separação $a = 0,001$ pc (período orbital $\sim 3$ anos) dentro de um disco $\alpha$ dominado pela pressão do gás ($\alpha=0,1$).
• Espaço de Parâmetros: Os autores exploram três números de Mach ($M_a = \{12, 50, 80\}$), correspondendo a frações de Eddington efetivas variando de regimes super-Eddington ($\sim 10^5$) a quase-Eddington ($\sim 10$).
• Prescrição de Sink: Os autores implementam uma prescrição de "sink sem torque adiabático" modificada. Esta atualização corrige perdas espúrias de energia interna associadas à remoção de massa sem torque em simulações anteriores de evolução de energia, garantindo soluções de disco de buraco negro único estáveis fora do raio de sink.
• Teste de Convergência: As simulações foram executadas em múltiplas resoluções ($\Delta x = 0,015a, 0,0075a, 0,005a$) para verificar a convergência do parâmetro de torque e das taxas de acreção.

Principais Resultados

1. Persistência de Fases de Não-Acreção: O estudo confirma que fases de não-acreção de longa duração persistem mesmo com uma equação de energia autoconsistente. Embora o gás se acumule e aqueça na borda distante (apocentro) da cavidade circumbinária, o material que lança fluxos de acreção próximo ao pericentro da cavidade permanece comparativamente frio. Consequentemente, esses fluxos carecem de suporte de pressão interna suficiente para resistir à compactação de maré pelo potencial binário, levando a uma supressão sustentada da taxa de acreção do binário ($\dot{M}_b \ll \dot{M}_\infty$).
2. Moderação da Não-Acreção em Fuga: A inclusão da equação de energia enfraquece, mas não elimina, o problema da "não-acreção em fuga" encontrado em soluções isotérmicas. O parâmetro de torque ($\ell_0$), que governa a supressão da acreção, ainda cresce exponencialmente com o número de Mach, embora a uma taxa mais lenta ($\ln \ell_0 \propto 0,085 M_a$) em comparação com soluções isotérmicas ($\ln \ell_0 \propto 0,21 M_a$).
3. Estrutura Térmica e Formação de Fluxos: A cavidade exibe uma assimetria térmica distinta. O apocentro é aquecido por acúmulos de densidade e geração de fluxos, enquanto o pericentro resfria eficientemente via expansão adiabática conforme elementos de fluido retornam em órbitas excêntricas. Este mecanismo de resfriamento é dominante sobre as perdas radiativas na região de lançamento de fluxos.
4. Propriedades Eletromagnéticas: Apesar das baixas taxas de acreção, esses sistemas não são eletromagneticamente escuros. A borda interna aquecida do CBD permanece luminosa nas frequências óptica e do infravermelho próximo (NIR). No entanto, os sistemas são intrinsecamente fracos em raios-X devido à truncagem do disco interno e à falta de um contínuo de alta energia e fotoionizante.
5. Tempos de Fusão vs. Ativação: Para números de Mach elevados ($M_a \gtrsim 90$), o tempo necessário para o binário "ligar" (atingir um estado de acreção estacionário) pode exceder o tempo de vida de fusão do binário impulsionado por ondas gravitacionais. Isso sugere que binários massivos em discos dominados pela pressão do gás podem fundir-se enquanto ainda estão em uma fase de não-acreção.

Significância e Alegações
O artigo afirma que o problema da "não-acreção em fuga" é uma característica robusta de discos finos dominados pela pressão do gás, sobrevivendo mesmo quando o feedback térmico é incluído. Isso tem implicações profundas para a identificação eletromagnética de SMBHBs:

• Detectabilidade: Binários em fase de não-acreção podem ser detectáveis em futuros levantamentos de campo amplo (ex: LSST, Telescópio Roman) como fontes periodicamente variáveis nas bandas óptica e NIR, apesar de sua baixa luminosidade bolométrica.
• Assinaturas Espectrais: Esses sistemas podem aparecer como Núcleos Galácticos Ativos (AGNs) intrinsecamente fracos em raios-X, com linhas de emissão fracas ou ausentes, potencialmente ligando-os à população de quasares de linhas fracas.
• Fenômenos de Mudança de Aspecto (Changing-Look): O alimentação intermitente de minidiscos, impulsionada pela complexa dinâmica da cavidade truncada, oferece um mecanismo plausível para o comportamento de AGNs do tipo "changing-look", onde os fluxos ionizantes variam significamente ao longo do tempo.

Os autores concluem que, embora seus resultados sejam baseados em modelos de média vertical, não magnetizados e dominados pela pressão do gás com resfriamento simplificado, eles fornecem uma linha de base crítica para entender como a autoconsistência termodinâmica afeta a dinâmica de acreção binária. Eles sugerem que trabalhos futuros devem incorporar campos magnéticos, pressão de radiação e geometria 3D para caracterizar totalmente a transição para o estado de acreção estacionária e a fase final de fusão.