Dust Formation in Common Envelope Binary Interactions -- III. Lightcurves

Résumé technique : Formation de poussière dans les interactions binaires d'enveloppe commune — III. Courbes de lumière

Énoncé du problème
Les Novae Rouges Lumineuses (LRN) sont des événements transitoires résultant d'interactions d'enveloppes communes (CE), pourtant la prédiction de leurs courbes de lumière observables à partir de simulations hydrodynamiques reste un défi. Les tentatives précédentes pour modéliser ces événements traitaient souvent la formation de poussière comme un processus instantané ou reposaient sur des tables d'opacité qui ne résolvaient pas explicitement l'échelle de temps de la nucléation. De plus, combler efficacement le fossé entre les simulations hydrodynamiques 3D et les données observationnelles nécessite des méthodes de transfert radiatif post-traitement robustes, capables de rendre compte de la nature complexe, asymétrique et évolutive des éjectas de la CE, en particulier la formation d'une coquille de poussière optiquement épaisse.

Méthodologie
Les auteurs effectuent des calculs de courbes de lumière par post-traitement sur deux simulations 3D de type SPH (Smoothed Particle Hydrodynamics) d'interactions de CE, décrites initialement dans Bermúdez-Bustamante et al. (2024). Les simulations modélisent l'interaction entre un compagnon compact ($0,6 M_\odot$) et des étoiles de la branche asymptotique des géantes (AGB) de deux masses distinctes : $1,7 M_\odot$ et $3,7 M_\odot$. Crucialement, ces simulations incluent une nucléation explicite de poussière à base de carbone.

La méthodologie de dérivation des courbes de lumière comprend :

1. Transfert radiatif : Application de l'équation du transfert radiatif le long des rayons projetés depuis le domaine de simulation vers un observateur. L'intensité spécifique est calculée en utilisant une approximation d'opacité grise supposant un équilibre thermodynamique local (LTE).
2. Traitement de l'opacité : L'opacité totale est la somme de l'opacité du gaz (dérivée des tables MESA pour $T > 1450$ K) et de l'opacité de la poussière (calculée via la théorie de Mie pour la poussière de carbone). Pour les températures inférieures à 1450 K, là où la poussière se forme, l'opacité du gaz est fixée à une valeur constante faible, permettant à l'opacité de la poussière de dominer.
3. Calcul de la luminosité : La luminosité bolométrique est dérivée en intégrant l'intensité spécifique sur une grille de $256 \times 256$ pixels. La fonction source est traitée comme un corps noir gris.
4. Analyse de la résolution et des incertitudes : Les auteurs quantifient explicitement les incertitudes découlant de la faible résolution de surface (où peu de particules SPH définissent la photosphère) et de l'hypothèse adiabatique (absence de refroidissement radiatif). Ils introduisent une métrique de résolution ($N_{res}$) basée sur le nombre de particules contribuant à l'intégration de la luminosité par rayon pour définir les barres d'erreur.

Résultats clés
Les simulations, exécutées pendant 44 ans sous conditions adiabatiques, produisent les caractéristiques de courbe de lumière suivantes :

• Phase précoce (Montée vers le pic) : Les deux modèles présentent un pic chaud et brillant ($L \sim 2-4 \times 10^6 L_\odot$, $T_{eff} \sim 6000-7000$ K) survenant 1 à 5 ans après le début du débordement du lobe de Roche. Cette montée est principalement pilotée par l'expansion rapide de la photosphère due à l'entraînement (inspiral) et aux instabilités de surface, plutôt que par des changements de température significatifs. Le timing et la forme du pic varient selon l'angle de vue et la masse de la progénitrice ; le modèle plus massif ($3,7 M_\odot$) présente un pic plus large et plus irrégulier en raison d'un rapport de masse plus faible ($q=0,16$) et d'une formation de poussière retardée.
• Formation de poussière et déclin : La nucléation de la poussière commence environ 1 à 3 ans après le début de l'interaction. Peu après, la luminosité bolométrique subit un déclin brutal à mesure que la coquille de poussière devient optiquement épaisse, obscurcissant la photosphère de gaz chaud. Ceci est suivi d'une récupération partielle et d'une phase de plateau avec une température effective d'environ 400 K.
• Évolution de la photosphère : La photosphère « géométrique » (définie par $\tau = 2/3$) s'étend linéairement, atteignant $\sim 250$ au à la fin de la simulation. Cependant, la photosphère « observée » (dérivée de l'ajustement SED) révèle une évolution complexe où une composante plus froide, dominée par la poussière, finit par surpasser la composante de gaz plus chaude.
• Dispersion à long terme : L'extrapolation de l'expansion libre de la coquille de poussière suggère que la poussière deviendra optiquement mince aux longueurs d'onde visibles entre 100 et 200 ans après l'interaction, révélant à nouveau la photosphère interne plus chaude.
• Comparaison avec les observations : Les auteurs comparent leurs résultats à OGLE-2002-BLG-360 et AT 2025abao. Le modèle de $3,7 M_\odot$ offre une correspondance plus étroite avec AT 2025abao en termes de masse de la progénitrice et de luminosité de pic, bien que le temps de montée observé pour AT 2025abao soit nettement plus rapide que ce que prédit la simulation. Le modèle de $1,7 M_\odot$ partage certaines caractéristiques avec OGLE-2002, notamment concernant la masse de poussière et le rayon de la photosphère à des temps tardifs.

Limites et incertitudes
Le papier identifie deux sources principales d'incertitude qui limitent la précision quantitative des courbes de lumière :

1. Résolution de surface : Les 1 à 2 premières années des courbes de lumière sont fortement affectées par une résolution SPH insuffisante à la surface de l'étoile. La résolution actuelle (1,3 million de particules) résulte en seulement 1 à 2 particules définissant la photosphère à $t=0$, ce qui entraîne de grandes incertitudes sur la luminosité initiale et le temps de montée. Les auteurs notent que l'obtention d'une résolution convergée nécessiterait des centaines de millions, voire des milliards de particules.
2. Hypothèse adiabatique : Les simulations sont adiabatiques, manquant de refroidissement radiatif. Par conséquent, la luminosité du plateau à un stade tardif ($\sim 10^5 L_\odot$) est probablement surestimée d'un ordre de grandeur par rapport aux attentes physiques, car le système ne peut pas rayonner l'excès d'énergie interne.

Signification
Ce travail représente une première étape pour combler le fossé entre les simulations hydrodynamiques de CE en 3D avec nucléation explicite de poussière et les courbes de lumière observables. Bien que les auteurs reconnaissent que leurs modèles ne sont pas encore totalement convergés ou physiquement complets (en raison de l'hypothèse adiabatique), ils démontrent avec succès que :

• La formation de poussière dans les interactions de CE est un processus non instantané qui modifie fondamentalement la morphologie de la courbe de lumière, créant un déclin et une phase de plateau distincts.
• L'asymétrie de l'interaction entraîne des dépendances significatives vis-à-vis de l'angle de vue observé.
• Le timing de la formation de la poussière et la transition subséquente d'une photosphère dominée par le gaz vers une photosphère dominée par la poussière peuvent être suivis, fournissant un cadre pour interpréter les futures observations de LRN.

Le papier conclut que, bien que les limites actuelles empêchent une correspondance quantitative précise avec des transitoires spécifiques, les caractéristiques qualitatives des courbes de lumière (pic, déclin, plateau et échelles de temps de dispersion de la poussière) offrent un modèle précieux pour comprendre l'évolution des événements de CE et la formation des nébuleuses planétaires. Les améliorations futures nécessiteront un raffinement adaptatif des particules pour résoudre la surface et l'implémentation du refroidissement radiatif pour corriger l'énergie à des stades tardifs.