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

Résumé technique : Un modèle de disque d'accrétion d'objet stellaire jeune paramétré avec un taux d'accrétion croissant

Énoncé du problème
Les objets stellaires jeunes (YSO), en particulier la sous-classe des FU Orionis (FU Ori), subissent des éruptions d'accrétion dramatiques caractérisées par des augmentations de luminosité de plusieurs ordres de grandeur sur des mois à des années, suivies de déclins sur des décennies. Bien que la transition de l'accrétion magnétosphérique en état faible à l'accrétion d'éruption en état élevé soit un phénomène reconnu, la compréhension de l'évolution spécifique des composants du système — photosphère stellaire, chocs magnétosphériques, disques de gaz visqueux et disques de poussière passifs — durant ces événements reste difficile. Les modèles existants reposent souvent sur des simulations numériques complexes et coûteuses en calcul (par exemple, MHD) ou sur des hypothèses d'état stationnaire qui peuvent ne pas capturer les changements morphologiques dépendants du temps observés dans les courbes de lumière multi-bandes. Il existe un besoin d'un cadre paramétré couplant ces composants physiques pour simuler l'évolution temporelle des distributions d'énergie spectrale (SED) et des courbes de lumière en réponse à des taux d'accrétion variables.

Méthodologie
Les auteurs présentent un pipeline de modélisation paramétré construit sur le code YSOpy, conçu pour simuler la photométrie et les spectres en série temporelle pour les étoiles jeunes accrétant. Contrairement aux simulations numériques complètes de transfert radiatif ou auto-cohérentes, cette approche utilise une infrastructure modulaire pour calculer les spectres de quatre composants distincts du système séparément avant de les sommer :

1. Photosphère stellaire : Modélisée comme un corps noir ou une couche atmosphérique.
2. Chocs d'accrétion magnétosphériques : Modélisés à l'aide d'une approximation de couche d'hydrogène où la température ($T_{hotspot}$), la densité électronique ($n_e$) et l'épaisseur optique ($\tau$) sont déduites de manière auto-cohérente du taux d'accrétion ($\dot{M}$) et de la vitesse de chute libre ($v_{ff}$).
3. Disque de gaz chauffé par viscosité : Modélisé à l'aide d'un profil de température modifié de Shakura & Sunyaev (1973) pour le disque interne, utilisant des grilles atmosphériques BT-Settl pour les températures $>1400$ K et des corps noirs pour les régions plus froides.
4. Disque de poussière chauffé passivement : Modélisé à l'aide d'un profil de température de Chiang & Goldreich (1997) pour la poussière irradiée, avec une limite intérieure déterminée par le rayon de sublimation de la poussière ($R_{sub}$).

Le pipeline itère sur un profil de taux d'accrétion dépendant du temps ($\dot{M}(t)$). À chaque pas de temps, le modèle calcule dynamiquement des rayons critiques, notamment le rayon de truncation magnétosphérique ($R_{trunc}$), le rayon de co-rotation ($R_{co-rot}$) et le rayon interne du disque de poussière ($R_{dust,inner}$). Le rayon interne du disque de gaz ($R_{in}$) est déterminé par l'équilibre entre les pressions magnétique et d'accrétion, passant de $R_{co-rot}$ à faible $\dot{M}$ à $R_{trunc}$ (ou $R_*$) à fort $\dot{M}$. Le modèle prend en compte la « coupure » de l'accrétion magnétosphérique polaire lorsque le disque atteint la surface stellaire ($R_{in} \approx R_*$). Enfin, les spectres haute résolution résolus dans le temps sont convolués avec des fonctions de transmission de filtres standards (par exemple, Gaia, 2MASS, WISE) pour générer des courbes de lumière et des courbes de couleur multi-bandes.

Contributions clés

• Cadre paramétré intégré : L'article introduit un pipeline rationalisé qui couple les composants stellaire, de choc, gazeux et de poussière en un seul modèle évoluant dans le temps, permettant la simulation de cycles complets d'éruption sans le coût computationnel des simulations MHD 3D.
• Physique de choc auto-cohérente : Le modèle dérive les températures des points chauds et les densités électroniques directement à partir du taux d'accrétion et de la vitesse de chute libre, assurant la conservation de l'énergie ($L_{slab} = L_{blackbody}$) via un facteur d'échelle $\chi$.
• Limites de disque dynamiques : Le modèle suit explicitement l'évolution de la limite interne du disque de gaz et du front de sublimation de la poussière en fonction de $\dot{M}$, capturant la transition d'un disque limité par la co-rotation à un disque tronqué magnétiquement, et enfin à un disque touchant la surface stellaire.
• Outils de diagnostic multi-bandes : L'infrastructure génère des courbes de lumière et des courbes de couleur synthétiques dans les bandes optiques et infrarouges, fournissant un outil de diagnostic pour interpréter les morphologies d'éruption observées.

Résultats
Les auteurs ont appliqué le modèle à un taux d'accrétion augmentant linéairement et à quatre profils d'accrétion spécifiques imitant des objets FU Ori observés : HBC 722 (linéaire en trois parties), V890 Aur (tangente hyperbolique), V960 Mon (exponentielle bilatérale) et Gaia 17bpi (logistique × Gaussienne). Les principaux résultats incluent :

• Évolution de la dominance des composants :

  • Faible $\dot{M}$ : Le système est dominé par la photosphère stellaire et le disque de poussière passif ; les chocs magnétosphériques sont négligeables ou plus froids que la photosphère.
  • $\dot{M}$ croissant : À mesure que $\dot{M}$ augmente, les points chauds magnétosphériques s'illuminent rapidement, dominant d'abord les bandes UV et optique bleue (par exemple, la bande u).
  • Pic/Élevé $\dot{M}$ : Le chauffage visqueux du disque de gaz interne devient la source de flux dominante, en particulier dans l'infrarouge proche (bandes Gaia, 2MASS). Si $\dot{M}$ dépasse $\sim 10^{-5} - 10^{-4} M_\odot$ an$^{-1}$, le rayon de truncation du disque se déplace vers l'intérieur jusqu'à la surface stellaire, coupant l'accrétion polaire et l'émission des points chauds.
  • Comportement en infrarouge moyen : Les bandes infrarouges moyennes (WISE) répondent principalement au chauffage et à l'expansion du front de sublimation de la poussière, s'illuminant souvent avant les bandes optiques à mesure que le rayon interne du disque de poussière se déplace vers l'extérieur.

• Morphologie des courbes de lumière :

  • Les courbes de lumière optique rouge et infrarouge proche suivent généralement le profil du taux d'accrétion d'entrée.
  • Les courbes de lumière en infrarouge moyen sont plus sensibles à la géométrie et au chauffage du disque de poussière interne, présentant parfois un éclaircissement « précurseur » ou un timing distinct par rapport aux bandes optiques.
  • L'évolution des couleurs varie : les systèmes deviennent souvent « plus bleus quand ils sont plus brillants » durant la montée en raison du chauffage des points chauds et du disque interne, mais peuvent rougir durant le déclin à mesure que le disque visqueux se refroidit et que la photosphère réapparaît.

• Sensibilité aux paramètres :

  • Intensité du champ magnétique : Des champs plus forts retardent l'approche du disque vers la surface stellaire, prolongeant la phase d'accrétion magnétosphérique et créant des courbes de lumière à pics multiples dans l'UV/optique.
  • Période de rotation : Les variations de la période de co-rotation affectent principalement le rayon interne du disque en état faible, mais ont un impact négligeable sur les courbes de lumière durant la phase d'éruption à fort taux d'accrétion.

Signification et affirmations
Les auteurs positionnent ce travail comme un outil illustratif pour interpréter les mécanismes physiques derrière les éruptions FU Ori plutôt que comme un modèle prédictif pour des objets individuels spécifiques. Ils affirment que leur approche paramétrée démontre avec succès comment différents composants du système dominent le flux à différentes longueurs d'onde tout au long d'un cycle d'éruption.

L'article souligne que le timing de l'éclaircissement à travers les longueurs d'onde offre un pouvoir de diagnostic : l'éclaircissement UV/optique bleu indique le début du chauffage par choc magnétosphérique, tandis que l'éclaircissement en infrarouge moyen signale l'expansion du front de sublimation de la poussière. Les auteurs suggèrent que les futures installations comme LSST et UVEX pourraient utiliser ces signatures prédites pour identifier les éruptions à un stade précoce. De plus, le modèle fournit un cadre pour inférer les profils d'accrétion sous-jacents et les taux d'accrétion en état faible à partir des courbes de lumière et de couleur observées, comblant le fossé entre les descriptions phénoménologiques simples et les simulations physiques complexes. Les auteurs reconnaissent que, bien que leur modèle simplifie la restructuration du champ magnétique et les couches limites d'accrétion équatoriales, il capture la physique essentielle de la transition des états d'accrétion de type T-Tauri aux états de type FU Ori.