Evolution of wide O star binaries through their LBV stage. Population synthesis with mass-ejection-driven orbital evolution

1. Il Problema Scientifico

Il lavoro affronta una discrepanza significativa tra le previsioni della sintesi stellare e le osservazioni astronomiche riguardanti le binarie di stelle Wolf-Rayet (WR) a lungo periodo.

• Il Paradosso: I modelli evolutivi standard prevedono che le binarie di stelle O massicce e inizialmente ampie, dopo aver subito trasferimento di massa, dovrebbero evolvere in sistemi WR+O con periodi orbitali lunghi. Tuttavia, tali sistemi sono estremamente rari nelle osservazioni (specialmente nelle Nubi di Magellano).
• Conseguenze: Questa scarsità osservativa pone vincoli severi sui canali evolutivi che portano alla formazione di buchi neri binari (BBH) che si fondono, in particolare il canale dell'evoluzione dell'involucro comune (CEE), che richiede progenitori ampi per sopravvivere alla fase di CEE.
• Il Caso WR 140: Il sistema galattico WR 140 (periodo di 2895 giorni, eccentricità $e \approx 0.9$) presenta caratteristiche (ampia separazione e alta eccentricità) difficili da spiegare con i modelli di trasferimento di massa classici, che tendono a circularizzare le orbite.

2. Metodologia e Ipotesi Fondamentali

Gli autori propongono un nuovo meccanismo di evoluzione orbitale guidato dall'eiezione di massa durante la fase di Variabile Blu Luminosa (LBV) della stella primaria, invece del classico trasferimento di massa tramite overflow del lobo di Roche (RLOF).

• Ipotesi Principale: Nelle binarie ampie ed eccentriche, la stella primaria massiccia può espandersi fino a raggiungere il limite di Eddington prima di iniziare il trasferimento di massa. Invece di un flusso stabile, l'involucro gonfiato e debolmente legato viene espulso istantaneamente al passaggio pericentro.
• Modellazione Analitica:
  • Gli autori derivano equazioni analitiche per l'evoluzione orbitale post-eiezione, assumendo che l'eiezione avvenga istantaneamente e isotropicamente al pericentro.
  • Dimostrano che l'eiezione di massa al pericentro conserva la direzione del vettore di Laplace-Runge-Lenz, mantenendo l'orientamento dell'orbita, ma modifica drasticamente i parametri orbitali.
  • Le equazioni mostrano che l'eiezione di massa aumenta il periodo orbitale ($P_{orb}$) e l'eccentricità ($e$), e conferisce una "spinta" (kick) al centro di massa del sistema, generando velocità spaziali misurabili.
• Sintesi di Popolazione:
  • Vengono eseguiti calcoli di sintesi di popolazione per le stelle WR nella Piccola Nube di Magellano (SMC).
  • Il modello utilizza una griglia densa di modelli evolutivi binari dettagliati (da lavori precedenti degli stessi autori) con metallicità SMC.
  • Vengono introdotti criteri critici: la massa primaria deve superare $M_{crit} = 30 M_\odot$ e il periodo iniziale deve essere $P_{crit} > 40$ giorni per innescare la fase LBV e l'eiezione di massa.
  • Vengono considerate diverse distribuzioni di eccentricità iniziale (basate su Sana et al. 2012) e parametri critici variabili.

3. Contributi Chiave e Risultati

A. Dinamica dell'Evoluzione

• Aumento di Periodo ed Eccentricità: A differenza del trasferimento di massa (che spesso riduce il periodo e circularizza l'orbita), l'eiezione di massa al pericentro allarga l'orbita e ne aumenta l'eccentricità.
• Velocità Spaziali: Il sistema riceve un impulso di quantità di moto, generando velocità di fuga misurabili per le stelle risultanti (sia le binarie legate che le stelle singole "runaway").
• Distruzione del Sistema: Se la massa espulsa supera una certa soglia (dipendente dall'eccentricità iniziale), la binaria si disgrega, producendo stelle WR singole ad alta velocità.

B. Applicazione a WR 140

• Il modello riesce a riprodurre coerentemente le proprietà di WR 140 (periodo $\sim$2900 giorni, $e \approx 0.9$).
• I risultati indicano che il progenitore di WR 140 aveva un periodo orbitale iniziale inferiore a $\sim 100$ giorni e una massa di circa 30-45 $M_\odot$.
• Il modello predice una velocità spaziale attuale per WR 140 di 50-75 km/s, compatibile con le osservazioni del moto proprio ($\sim 34$ km/s), suggerendo che il sistema sia un risultato tipico di questa evoluzione.

C. Risultati della Sintesi di Popolazione (SMC)

Il modello "Simulation1" (con $M_{crit}=30 M_\odot$, $P_{crit}=40$ d, e distribuzione di eccentricità iniziale) predice:

• 5.3 binarie WR+O strette (periodo < 20 giorni), formate tramite trasferimento di massa classico.
• 3.8 binarie WR+O ampie (periodo > 20 giorni, picco a $\sim 3000$ giorni), formate tramite il canale di eiezione di massa. Queste hanno alte eccentricità ($0.4 - 1.0$).
• 1.9 stelle WR "runaway" (singole), risultanti dalla distruzione della binaria.
• Confronto con le osservazioni: Questi numeri si allineano bene con il campione osservato nella SMC (4 binarie WR+O strette e 7 stelle WR apparentemente singole, di cui 2 runaway).
• Implicazione Osservativa: Le binarie WR+O a lungo periodo predette sono difficili da rilevare con le ricerche di velocità radiale tradizionali a causa delle loro alte eccentricità e lunghi periodi, il che spiega la loro apparente assenza.

D. Applicazione ai Sistemi Gaia BH

Il modello viene applicato ai sistemi di buchi neri scoperti da Gaia (Gaia BH1 e BH2), che mostrano orbite ampie ed eccentriche.

• Il meccanismo di eiezione di massa LBV offre un canale evolutivo plausibile per spiegare l'origine di questi sistemi, permettendo di "pompare" l'eccentricità e allargare l'orbita prima del collasso in buco nero.

4. Significato e Implicazioni

1. Risoluzione del "Missing Link": Il modello spiega perché le binarie WR a lungo periodo siano rare nelle osservazioni: non sono inesistenti, ma sono nascoste a causa delle loro orbite altamente eccentriche e dei lunghi periodi, oppure sono state distrutte diventando stelle runaway.
2. Impatto sui Buchi Neri Binari (BBH):
   • Se l'evoluzione guidata dall'eiezione di massa è dominante nelle binarie ampie, il numero di progenitori per la fase di CEE (che richiede binarie WR+O strette e circolari) potrebbe essere inferiore alle attese.
   • Le binarie che sopravvivono alla fase LBV tendono ad essere molto eccentriche. La formazione di BBH che si fondono attraverso il canale CEE richiederebbe che le maree riescano a circularizzare queste orbite eccentriche prima del collasso in buco nero, una condizione non banale.
3. Nuovi Canali Evolutivi: Il lavoro suggerisce che l'interazione mareale durante il passaggio pericentro in sistemi LBV è un meccanismo cruciale, spesso trascurato, che può dominare l'evoluzione di binarie massicce ampie, portando a sistemi con proprietà orbitali uniche (come WR 140) e a stelle singole ad alta velocità.

In conclusione, questo studio offre un quadro teorico coerente che collega l'evoluzione delle binarie O massicce, la scarsità osservata di binarie WR a lungo periodo, la presenza di stelle WR runaway e la formazione di sistemi di buchi neri eccentrici come quelli rilevati da Gaia.