From Common Envelope Evolution to Luminous Red Novae I: A One-dimensional Radiation Hydrodynamic Model

1. Il Problema Scientifico

L'evoluzione della comune envelope (CEE) è una delle fasi meno comprese nell'evoluzione delle stelle binarie. Durante la fase di "plunge-in" (immersione), i nuclei binari rilasciano energia orbitale che viene convertita in energia cinetica, termica, magnetica e radiativa. Se questa energia viene depositata rapidamente in un ambiente otticamente spesso, la pressione di radiazione può superare quella del gas, promuovendo l'espulsione dell'involucro.
Tuttavia, esiste un divario significativo tra la teoria e le osservazioni (come le Novae Rosse Luminose o LRNe). I modelli tridimensionali (3D) attuali spesso trascurano il trasporto radiativo per limiti computazionali, assumendo un'evoluzione adiabatica. Questa semplificazione fallisce quando l'involucro si espande e diventa luminoso, ignorando il ruolo cruciale della pressione di radiazione e del trasporto di energia radiativa nel guidare l'espulsione e determinare le firme osservative.

2. Metodologia

L'autore ha sviluppato e utilizzato un modello idrodinamico radiativo (RHD) unidimensionale (1D) implementato nel codice Guangqi.

• Equazioni Governative: Il modello risolve le equazioni RHD nel sistema di riferimento di laboratorio utilizzando l'approssimazione della diffusione limitata dal flusso (FLD - Flux-Limited Diffusion). Le equazioni includono la conservazione di massa, quantità di moto ed energia, accoppiando l'energia del gas e quella della radiazione.
• Fisica Inclusa:
  • Equazione di Stato (EoS): Include l'energia interna del gas (idrogeno ed elio in varie fasi di ionizzazione) e l'energia di ricombinazione ($H \to H^+$, $H \to H_2$).
  • Trasporto Radiativo: Utilizza tabelle di opacità combinate (polvere e gas) e calcola l'accelerazione dovuta alla pressione di radiazione ($a_{rad}$) in modo auto-consistente.
  • Geometria: Simmetria sferica con un dominio computazionale che va da $r_{in} = 450 R_\odot$ a $r_{out} = 1.5 \times 10^4 R_\odot$.
• Condizioni al Contorno:
  • Esterno: Condizioni di flusso libero (otticamente sottile).
  • Interno: Durante la fase di espansione, le condizioni corrispondono alle proprietà fisiche dell'involucro iniettato. Successivamente, si passa a una condizione di solo afflusso (inflow-only) per permettere il ricadimento (fallback) del materiale.
• Indagine Parametrica: È stata eseguita una vasta survey di parametri per quantificare l'impatto su:
  • Rapporto tra energia interna radiativa e del gas ($E/e_g \in [0.2, 3.2]$).
  • Rapporto tra l'energia totale dell'involucro e l'energia di legame gravitazionale ($\zeta \in [0.54, 2.87]$).
  • Tasso di iniezione di massa ($\dot{M} \in [2.5, 10] M_\odot/\text{yr}$).
  • Regimi di espansione subsonici e supersonici ($v_{ej}/v_{esc} \in [0.3, 0.6]$).

3. Contributi Chiave

Il lavoro introduce un approccio 1D auto-consistente che integra la pressione di radiazione in modo dinamico, superando le limitazioni dei modelli puramente adiabatici o puramente idrodinamici.

• Modellazione Auto-consistente: A differenza di studi precedenti (es. CI24), questo modello calcola l'accelerazione radiativa e il trasporto di energia in modo accoppiato, permettendo di studiare come la radiazione influenzi la dinamica del gas in tempo reale.
• Analisi di Convergence e Errori: Il paper include un'analisi rigorosa della conservazione dell'energia e dei test di convergenza numerica, dimostrando che il modello è robusto durante la fase dinamica, sebbene vi siano piccole sottostime di energia in regimi dominati dalla radiazione o subsonici.

4. Risultati Principali

Dall'analisi delle simulazioni emergono tre conclusioni fondamentali:

1. Accelerazione Dominata dalla Pressione di Radiazione:
   Nella regione ad alta opacità e alta luminosità immediatamente sotto il fronte di ricombinazione, la pressione di radiazione diventa il motore principale dell'espulsione di massa. L'accelerazione radiativa ($a_{rad}$) supera significativamente la gravità in queste zone, specialmente quando il rapporto $E/e_g$ è alto. Questo meccanismo è responsabile di una frazione di massa non legata ($\eta$) maggiore rispetto ai modelli puramente idrodinamici.

2. Relazione Non Lineare tra $\eta$ e $\zeta$:
   La frazione di massa espulsa ($\eta$) mostra una dipendenza non lineare dal parametro energetico $\zeta$ (rapporto energia totale/energia di legame).

   • Per $\zeta < 0.7$, l'espulsione è inefficiente.
   • Per $\zeta \in [0.7, 1.05]$, l'accelerazione diventa efficace.
   • Un tasso di iniezione di massa ($\dot{M}$) più alto porta a un $\eta$ maggiore, poiché l'alta opacità di un involucro denso intrappola la radiazione, trasferendo più momento al gas. Al contrario, in involucri tenui, la radiazione si disaccoppia e sfugge più facilmente.

3. Plateau Secondario nelle Curve di Luce:
   La ricombinazione dell'idrogeno atomico in molecolare ($H \to H_2$) rilascia calore latente che sostiene un secondo plateau nelle curve di luce bolometriche a tempi tardivi (tra 140 e 280 giorni). Questo fenomeno è correlato all'aumento della frazione di massa di $H_2$ e potrebbe spiegare le emissioni infrarosse osservate in eventi transitori come gli "SPRITEs".

5. Significato e Implicazioni

• Colmare il Divario Teoria-Osservazione: Il modello fornisce un collegamento quantitativo tra i parametri fisici della CEE e le firme osservative delle LRNe, aiutando a interpretare la diversità delle curve di luce osservate.
• Raffinamento del Formalismo $\alpha_{CE}$: La relazione non lineare trovata tra l'energia depositata e la massa espulsa suggerisce che il parametro di efficienza $\alpha$ nel formalismo standard di Webbink potrebbe non essere costante, ma dipendere da fattori come il tasso di iniezione di massa e il trasporto radiativo.
• Limiti e Futuro: Sebbene il modello 1D sia efficace per studiare la termodinamica e il trasporto radiativo, non può catturare il trasporto di momento angolare o le asimmetrie geometriche (che richiedono modelli 2D/3D). Tuttavia, questo lavoro stabilisce una base fisica solida per future indagini multidimensionali, sottolineando che la pressione di radiazione è un fattore critico spesso sottostimato nei modelli di CEE.

In sintesi, questo studio dimostra che l'inclusione auto-consistente del trasporto radiativo è essenziale per comprendere correttamente la dinamica di espulsione dell'involucro e le conseguenti firme osservative delle Novae Rosse Luminose.