Circumbinary disk formation through AGB star winds

Sintesi Tecnica: Formazione di Dischi Circumbinari attraverso i Venti di Stelle AGB

Problema
I dischi circubinari sono frequentemente osservati attorno a sistemi binari post-AGB (post-asymptotic giant branch), tuttavia i meccanismi che ne governano la formazione e l'evoluzione a lungo termine rimangono scarsamente compresi. Sebbene studi precedenti abbiano esplorato il trasferimento di massa tramite fasi di involucro comune (common-envelope) o il superamento del limite di Roche (RLOF) in sistemi semi-distaccati, il ruolo specifico delle interazioni vento-binaria nel generare questi dischi su lunghi periodi temporali non è stato pienamente caratterizzato. Le simulazioni esistenti spesso operano per meno di $\sim100$ periodi orbitali, limitando la capacità di catturare l'evoluzione secolare della morfologia del disco. Inoltre, la transizione tra diversi regimi di trasferimento di massa — specificamente l'accrescimento di Bondi-Hoyle (BH) rispetto al Wind Roche-Lobe Overflow (WRLOF) — e il modo in cui i parametri binari (eccentricità e rapporto di massa) influenzano la struttura del disco risultante richiedono un'indagine sistematica.

Metodologia
Gli autori utilizzano simulazioni di idrodinamica a particelle smoothed (SPH) utilizzando il codice phantom per modellare l'interazione tra il vento di una stella AGB e una compagna binaria. Il sistema consiste in una stella AGB ($M_1$) e una compagna ($M_2$), trattate come particelle "sink" con raggi di accrescimento definiti ($R_{s,1} = 1.0$ AU, $R_{s,2} = 0.4$ AU). Le simulazioni si basano sui parametri del sistema AC Her, ma sono scalate in modo da essere indipendenti dalla scala (scale-free) in base all'energia specifica delle particelle iniettate.

Le caratteristiche metodologiche chiave includono:

• Iniezione del Vento: Particelle di gas vengono iniettate a $R_1 = 1.1$ AU dalla stella AGB. Lo studio investiga due regimi di iniezione definiti dal rapporto di energia specifica $\eta$:
  • $\eta > 1$: venti veloci dove tutte le particelle sono non legate (regime BH).
  • $\eta < 0$: venti lenti dove le particelle sono legate, portando al WRLOF.
  • $0 \leq \eta \leq 1$: regimi intermedi.
• Geometria di Iniezione: Vengono testate distribuzioni di iniezione sia uniformi che gaussiane (con bias equatoriale), con il modello gaussiano adottato per gli studi dettagliati sul disco per migliorare la risoluzione vicino al piano medio.
• Termodinamica: Viene utilizzata principalmente un'equazione di stato isotermica locale ($H/R = 0.1$), con una simulazione comparativa che utilizza un'equazione di stato isotermica binaria per testare le assunzioni termiche.
• Viscosità: Viene imposto un parametro di viscosità efficace di tipo Shakura-Sunyaev $\alpha \sim 0.1$, sebbene lo studio analizzi esplicitamente se le strutture fisiche (shock/spirali) dominino il trasporto del momento angolare rispetto a questa viscosità imposta.
• Durata: Le simulazioni durano per oltre 600 periodi orbitali binari ($\sim4500$ anni), superando significativamente la durata dei lavori precedenti.

Risultati Chiave
Le simulazioni rivelano che la morfologia dell'outflow è determinata dall'interazione tra la velocità del vento e la massa della compagna, categorizzata dal rapporto di energia specifica $\eta$:

1. Regime di Bondi-Hoyle (Venti Veloci/Compagne di Bassa Massa): Quando $\eta \gg 1$, il vento rimane ampiamente non legato. La compagna focalizza gravitazionalmente il vento, creando un persistente pattern di densità a spirale. Il flusso di gas è dominato dalla velocità radiale e il materiale tende a sfuggire direttamente dal sistema binario senza formare un disco circubinare stabile.
2. Regime WRLOF (Venti Lenti/Compagne di Alta Massa): Quando $\eta < 0$ (o intermedio), il vento lento e denso riempie il lobo di Roche del donatore. Un disco circumsingolare si forma attorno alla compagna, aumentando significativamente la sezione d'urto effettiva per la cattura della massa. Questo disco facilita un efficiente trasferimento di momento angolare, permettendo al materiale di sfuggire attraverso il punto di Lagrange esterno ($L_2$) e accumularsi in un disco circubinare coerente.
3. Impatto dell'Eccentricità Binaria ($e_b$):
   • Circolare/Moderatamente Eccentrica ($e_b = 0, 0.2$): I dischi si formano a piccoli raggi ($\sim10$ AU) ed espandono verso l'esterno. Il disco mantiene un'eccentricità moderata ($e_d \approx 0.2$) che si circolarizza gradualmente nel tempo a causa della competizione tra la dissipazione viscosa e i torques binari.
   • Altamente Eccentrica ($e_b = 0.4$): La compagna si avvicina alla stella AGB al perielio, permettendo a una grande frazione di materiale di essere iniettata direttamente in orbite circubinari altamente eccentriche. Ciò risulta in un disco concentrato a semiasse maggiore più elevato ($\sim30$ AU) con alta eccentricità ($e_d \approx 0.6$).
4. Impatto del Rapporto di Massa ($q$): Rapporti di massa inferiori (compagne più leggere) riducono l'efficienza della cattura del vento e della formazione del disco. Sebbene i dischi circubinari risultanti siano meno densi e crescano più lentamente, possono comunque formare strutture estese. Notevolmente, anche quando il disco circumsecondario non è risolto (nel caso con $q$ più basso), si sviluppa comunque un disco circubinare esteso e quasi circolare, suggerendo che un'alta densità iniziale non è strettamente necessaria per la formazione di dischi ampi.
5. Trasporto del Momento Angolare: L'analisi dello stress effettivo ($\alpha_{eff}$) indica che il trasporto del momento angolare all'interno del disco circubinare è dominato dalle strutture a spirale e dagli shock generati dall'interazione binaria-vento. Questi stress fisici sono comparabili o superiori alla viscosità artificiale imposta, implicando che la viscosità imposta non domina l'evoluzione del disco.

Significato e Rivendicazioni
L'articolo sostiene che le interazioni vento-binaria, in particolare nel regime WRLOF, possono naturalmente generare le diverse morfologie dei dischi circubinari osservate nei sistemi AGB e post-AGB senza richiedere una fase di involucro comune. Lo studio dimostra che:

• Sono necessarie simulazioni a lungo termine per catturare la transizione dalle strutture a spirale ai dischi stabili e per osservare l'evoluzione secolare dell'eccentricità e dell'estensione del disco.
• I parametri binari (eccentricità e rapporto di massa) controllano sensibilmente l'estensione radiale e l'eccentricità del disco, offrendo una potenziale spiegazione per la varietà di strutture di disco post-AGB osservate.
• Il momento angolare del risultante disco circubinare è fornito interamente dalla compagna, portando naturalmente alla formazione di dischi coplanari.

Gli autori notano esplicitamente che, sebbene il loro modello spieghi i dischi coplanari, esso non tiene conto dei dischi circubinari allineati ai poli osservati in sistemi come AC Her, lasciando aperta la questione dell'origine di tale momento angolare polare. Inoltre, lo studio riconosce che le semplificazioni nella termodinamica (assunzione isotermica) e l'omissione dell'accoppiamento esplicito polvere-radiazione potrebbero influenzare l'evoluzione dettagliata, sebbene il meccanismo globale di formazione del disco si ritenga robusto.