Circumbinary disk formation through AGB star winds

Resumo Técnico: Formação de Discos Circumbinários através de Ventos de Estrelas AGB

Definição do Problema
Discos circumbinários são frequentemente observados em torno de sistemas binários pós-AGB (pós-ramo assintótico gigante), contudo, os mecanismos que governam sua formação e evolução a longo prazo permanecem mal compreendidos. Embora estudos anteriores tenham explorado a transferência de massa via fases de envelope comum ou transbordamento do lóbulo de Roche (RLOF) em sistemas semitotalmente separados, o papel específico das interações vento-binária na geração desses discos ao longo de longos períodos não foi totalmente caracterizado. Simulações existentes frequentemente rodam por menos de $\sim100$ períodos orbitais, limitando sua capacidade de capturar a evolução secular da morfologia do disco. Além disso, a transição entre diferentes regimes de transferência de massa — especificamente a acreção de Bondi-Hoyle (BH) versus o transbordamento do lóbulo de Roche por vento (WRLOF) — e como os parâmetros binários (excentricidade e razão de massa) influenciam a estrutura resultante do disco requerem investigação sistemática.

Metodologia
Os autores utilizam simulações de hidrodinâmica de partículas suavizadas (SPH) usando o código phantom para modelar a interação entre o vento de uma estrela AGB e uma companheira binária. O sistema consiste em uma estrela AGB ($M_1$) e uma companheira ($M_2$), tratadas como partículas de sumidouro (sink particles) com raios de acreção definidos ($R_{s,1} = 1,0$ UA, $R_{s,2} = 0,4$ UA). As simulações baseiam-se nos parâmetros do sistema AC Her, mas são escalonadas de forma livre à escala dependendo da energia específica das partículas injetadas.

Principais características metodológicas incluem:

• Injeção de Vento: Partículas de gás são injetadas a $R_1 = 1,1$ UA da estrela AGB. O estudo investiga dois regimes de injeção definidos pela razão de energia específica $\eta$:
  • $\eta > 1$: Ventos rápidos onde todas as partículas são não ligadas (regime BH).
  • $\eta < 0$: Ventos lentos onde as partículas são ligadas, levando ao WRLOF.
  • $0 \leq \eta \leq 1$: Regimes intermediários.
• Geometria de Injeção: Distribuições de injeção uniforme e gaussiana (com viés equatorial) são testadas, sendo o modelo gaussiano adotado para estudos detalhados de disco para aumentar a resolução próximo ao plano médio.
• Termodinâmica: Uma equação de estado isotérmica local é utilizada primordialmente ($H/R = 0,1$), com uma simulação comparativa usando uma equação de estado binária isotérmica para testar as suposições térmicas.
• Viscosidade: Um parâmetro de viscosidade efetiva de Shakura-Sunyaev de $\alpha \sim 0,1$ é imposto, embora o estudo analise explicitamente se estruturas físicas (choques/espirais) dominam o transporte de momento angular sobre esta viscosidade imposta.
• Duração: As simulações rodam por mais de 600 períodos orbitais binários ($\sim4500$ anos), excedendo significativamente a duração de trabalhos anteriores.

Resultos Principais
As simulações revelam que a morfologia do fluxo de saída é determinada pela interação entre a velocidade do vento e a massa da companheira, categorizada pela razão de energia específica $\eta$:

1. Regime Bondi-Hoyle (Ventos Rápidos/Companheiras de Baixa Massa): Quando $\eta \gg 1$, o vento permanece amplamente não ligado. A companheira foca gravitacionalmente o vento, criando um padrão de densidade espiral persistente. O fluxo de gás é dominado pela velocidade radial, e o material tende a escapar do sistema binário diretamente sem formar um disco circumbinário estável.
2. Regime WRLOF (Ventos Lentos/Companheiras de Alta Massa): Quando $\eta < 0$ (ou intermediário), o vento lento e denso preenche o lóbulo de Roche da doadora. Um disco circunestelar se forma ao redor da companheira, aumentando significativamente a seção de choque efetiva para captura de massa. Este disco facilita o transporte eficiente de momento angular, permitindo que o material escape através do ponto de Lagrange externo ($L_2$) e se acumule em um disco circumbinário coerente.
3. Impacto da Excentricidade Binária ($e_b$):
   • Circular/Levemente Excêntrica ($e_b = 0, 0,2$): Discos se formam em pequenos raios ($\sim10$ UA) e expandem-se para fora. O disco mantém uma excentricidade moderada ($e_d \approx 0,2$) que circulariza gradualmente ao longo do tempo devido à dissipação viscosa competindo com os torques binários.
   • Altamente Excêntrica ($e_b = 0,4$): A companheira aproxima-se da estrela AGB no periastro, permitindo que uma grande fração do material seja injetada diretamente em órbitas circumbinárias altamente excêntricas. Isso resulta em um disco concentrado em semieixos maiores ($\sim30$ UA) com alta excentricidade ($e_d \approx 0,6$).
4. Impacto da Razão de Massa ($q$): Razões de massa menores (companheiras mais leves) reduzem a eficiência de captura de vento e de formação de disco. Embora os discos circumbinários resultantes sejam menos densos e cresçam mais lentamente, eles ainda podem formar estruturas estendidas. Notavelmente, mesmo quando o disco circunsecundário não é resolvido (no caso de $q$ mais baixo), um disco circumbinário estendido e quase circular ainda se desenvolve, sugerindo que uma alta densidade inicial não é estritamente necessária para a formação de um disco amplo.
5. Transporte de Momento Angular: A análise do estresse efetivo ($\alpha_{eff}$) indica que o transporte de momento angular dentro do disco circumbinário é dominado por estruturas espirais e choques gerados pela interação vento-binária. Esses estresses físicos são comparáveis ou superiores à viscosidade artificial imposta, implicando que a viscosidade imposta não domina a evolução do disco.

Significância e Alegações
O artigo afirma que as interações vento-binária, particularmente no regime WRLOF, podem naturalmente gerar as diversas morfologias de discos circumbinários observadas em sistemas AGB e pós-AGB sem exigir uma fase de envelope comum. O estudo demonstra que:

• Simulações de longo prazo são necessárias para capturar a transição de estruturas espirais para discos estáveis e para observar a evolução secular da excentricidade e extensão do disco.
• Parâmetros binários (excentricidade e razão de massa) controlam sensivelmente a extensão radial e a excentricidade do disco, oferecendo uma explicação potencial para a variedade de estruturas de discos pós-AGB observadas.
• O momento angular do disco circumbinário resultante é fornecido inteiramente pela companheira, levando naturalmente à formação de discos coplanares.

Os autores observam explicitamente que, embora seu modelo explique discos coplanares, ele não contabiliza os discos circumbinários alinhados ao polo observados em sistemas como AC Her, deixando a origem de tal momento angular polar como uma questão em aberto. Adicionalmente, o estudo reconhece que simplificações na termodinâmica (suposição isotérmica) e a negligência do acoplamento poeira-radiação explícito podem influenciar a evolução detalhada, embora o mecanismo global de formação do disco deva permanecer robusto.