A Parameterized YSO Accretion Disk Model with Increasing Accretion Rate: Predicted Outburst Lightcurves

Sintesi Tecnica: Un Modello di Disco di Accrescimento per Oggetti Stellari Giovani (YSO) Parametrizzato con Tasso di Accrescimento Crescente

Enunciato del Problema
Gli Oggetti Stellari Giovani (YSO), in particolare la sottoclasse FU Orionis (FU Ori), subiscono drammatici scoppi di accrescimento caratterizzati da aumenti di luminosità di diversi ordini di grandezza nell'arco di mesi o anni, seguiti da decadimenti che durano decenni. Sebbene la transizione dall'accrescimento magnetosferico nello stato di bassa luminosità all'accrescimento durante lo scoppietto nello stato di alta luminosità sia un fenomeno riconosciuto, comprendere l'evoluzione specifica dei componenti del sistema – fotosfera stellare, shock magnetosferici, dischi di gas viscosi e dischi di polvere passivi – durante questi eventi rimane una sfida. I modelli esistenti si basano spesso su simulazioni numeriche complesse e computazionalmente costose (ad es. MHD) o su assunzioni di stato stazionario che potrebbero non catturare i cambiamenti morfologici dipendenti dal tempo osservati nelle curve di luce multibanda. È necessario un framework parametrizzato che accoppi questi componenti fisici per simulare l'evoluzione temporale delle distribuzioni spettrali di energia (SED) e delle curve di luce in risposta a tassi di accrescimento variabili.

Metodologia
Gli autori presentano una pipeline di modellazione parametrizzata costruita sul codice YSOpy, progettata per simulare fotometrie in serie temporali e spettri per stelle giovani in accrescimento. A differenza delle simulazioni numeriche di trasferimento radiativo completo o autoconsistenti, questo approccio utilizza un'infrastruttura modulare per calcolare gli spettri per quattro distinti componenti del sistema separatamente, per poi sommarli:

1. Fotosfera Stellare: Modellata come un corpo nero o una lastra atmosferica.
2. Shock di Accrescimento Magnetosferico: Modellati utilizzando un'approssimazione di lastra di idrogeno, dove la temperatura ($T_{hotspot}$), la densità elettronica ($n_e$) e la profondità ottica ($\tau$) sono derivate in modo autoconsistente dal tasso di accrescimento ($\dot{M}$) e dalla velocità di caduta libera ($v_{ff}$).
3. Disco di Gas Riscaldato Viscosamente: Modellato utilizzando un profilo di temperatura modificato di Shakura & Sunyaev (1973) per il disco interno, impiegando griglie atmosferiche BT-Settl per temperature $>1400$ K e corpi neri per le regioni più fredde.
4. Disco di Polvere Riscaldato Passivamente: Modellato utilizzando un profilo di temperatura di Chiang & Goldreich (1997) per la polvere irradiata, con un confine interno determinato dal raggio di sublimazione della polvere ($R_{sub}$).

La pipeline itera attraverso un profilo di tasso di accrescimento dipendente dal tempo ($\dot{M}(t)$). Ad ogni passo temporale, il modello calcola dinamicamente i raggi critici, inclusi il raggio di troncamento magnetosferico ($R_{trunc}$), il raggio di co-rotazione ($R_{co-rot}$) e il raggio interno del disco di polvere ($R_{dust,inner}$). Il raggio interno del disco di gas ($R_{in}$) è determinato dall'equilibrio tra le pressioni magnetiche e di accrescimento, transitando da $R_{co-rot}$ a bassi $\dot{M}$ a $R_{trunc}$ (o $R_*$) ad alti $\dot{M}$. Il modello tiene conto dello "spegnimento" dell'accrescimento magnetosferico poloidale quando il disco raggiunge la superficie stellare ($R_{in} \approx R_*$). Infine, gli spettri ad alta risoluzione risolti nel tempo sono convoluti con le funzioni di trasmissione dei filtri standard (ad es. Gaia, 2MASS, WISE) per generare curve di luce e curve di colore multibanda.

Contributi Chiave

• Framework Parametrizzato Integrato: Il documento introduce una pipeline snella che accoppia i componenti stellari, degli shock, del gas e della polvere in un unico modello evolutivo nel tempo, consentendo la simulazione di interi cicli di scoppietto senza il costo computazionale delle simulazioni MHD 3D.
• Fisica degli Shock Autoconsistente: Il modello deriva le temperature delle macchie calde e le densità elettroniche direttamente dal tasso di accrescimento e dalla velocità di caduta libera, garantendo la conservazione dell'energia ($L_{slab} = L_{blackbody}$) tramite un fattore di scala $\chi$.
• Confini del Disco Dinamici: Il modello traccia esplicitamente l'evoluzione del confine interno del disco di gas e del fronte di sublimazione della polvere in funzione di $\dot{M}$, catturando la transizione da un disco limitato dalla co-rotazione a un disco troncato magneticamente e, infine, a un disco che tocca la superficie stellare.
• Strumenti Diagnostici Multibanda: L'infrastruttura genera curve di luce e curve di colore sintetiche nelle bande ottiche e infrarosse, fornendo uno strumento diagnostico per interpretare le morfologie degli scoppietti osservati.

Risultati
Gli autori hanno applicato il modello a un tasso di accrescimento che aumenta linearmente e a quattro profili di accrescimento specifici che mimano oggetti FU Ori osservati: HBC 722 (lineare a tre pezzi), V890 Aur (tangente iperbolica), V960 Mon (esponenziale a due lati) e Gaia 17bpi (logistica $\times$ Gaussiana). I risultati chiave includono:

• Evoluzione della Dominanza dei Componenti:

  • Basso $\dot{M}$: Il sistema è dominato dalla fotosfera stellare e dal disco di polvere passivo; gli shock magnetosferici sono trascurabili o più freddi della fotosfera.
  • $\dot{M}$ in Ascesa: All'aumentare di $\dot{M}$, le macchie calde magnetosferiche si illuminano rapidamente, dominando prima le bande UV e ottiche blu (ad es. banda u).
  • Picco/Alto $\dot{M}$: Il riscaldamento viscoso del disco di gas interno diventa la fonte di flusso dominante, in particolare nell'infrarosso vicino (bande Gaia, 2MASS). Se $\dot{M}$ supera $\sim 10^{-5} - 10^{-4} M_\odot$ yr$^{-1}$, il raggio di troncamento del disco si sposta verso l'interno fino alla superficie stellare, spegnendo l'accrescimento poloidale e l'emissione delle macchie calde.
  • Comportamento nell'Infrarosso Medio: Le bande dell'infrarosso medio (WISE) rispondono principalmente al riscaldamento e all'espansione del fronte di sublimazione della polvere, spesso illuminandosi prima delle bande ottiche mentre il raggio interno del disco di polvere si sposta verso l'esterno.

• Morfologia delle Curve di Luce:

  • Le curve di luce ottiche rosse e nell'infrarosso vicino tracciano generalmente il profilo del tasso di accrescimento in ingresso.
  • Le curve di luce nell'infrarosso medio sono più sensibili alla geometria e al riscaldamento del disco di polvere interno, mostrando talvolta un "pre-cursore" di illuminazione o una tempistica distinta rispetto alle bande ottiche.
  • L'evoluzione del colore varia: i sistemi spesso diventano "più blu quando più luminosi" durante la fase di ascesa a causa del riscaldamento delle macchie calde e del disco interno, ma possono arrossarsi durante il decadimento mentre il disco viscoso si raffredda e riemerge la fotosfera.

• Sensibilità ai Parametri:

  • Intensità del Campo Magnetico: Campi più forti ritardano l'avvicinamento del disco alla superficie stellare, prolungando la fase di accrescimento magnetosferico e creando curve di luce multibanda nell'UV/ottico.
  • Periodo di Rotazione: Le variazioni del periodo di co-rotazione influenzano principalmente il raggio interno del disco nello stato di bassa luminosità, ma hanno un impatto trascurabile sulle curve di luce durante la fase di scoppietto ad alto accrescimento.

Significato e Affermazioni
Gli autori collocano questo lavoro come uno strumento illustrativo per interpretare i meccanismi fisici alla base degli scoppietti FU Ori, piuttosto che come un modello predittivo per singoli oggetti specifici. Affermano che il loro approccio parametrizzato dimostra con successo come diversi componenti del sistema dominino il flusso a varie lunghezze d'onda durante un ciclo di scoppietto.

Il documento sottolinea che la tempistica dell'illuminazione attraverso le lunghezze d'onda offre potere diagnostico: l'illuminazione UV/ottica blu indica l'inizio del riscaldamento degli shock magnetosferici, mentre l'illuminazione nell'infrarosso medio segnala l'espansione del fronte di sublimazione della polvere. Gli autori suggeriscono che future strutture come LSST e UVEX potrebbero utilizzare queste firme previste per identificare scoppietti in fase iniziale. Inoltre, il modello fornisce un framework per inferire i profili di accrescimento sottostanti e i tassi di accrescimento nello stato di bassa luminosità dalle curve di luce e di colore osservate, colmando il divario tra descrizioni fenomenologiche semplici e simulazioni fisiche complesse. Gli autori riconoscono che, sebbene il loro modello semplifichi la ristrutturazione del campo magnetico e gli strati limite di accrescimento equatoriale, cattura la fisica essenziale della transizione dagli stati di accrescimento simili ai T-Tauri a quelli simili alle FU Ori.