Bayesian Inference of Stellar r-Process Abundances

Résumé Technique : Inférence Bayésienne des Abondances Stellaires du Processus r

Énoncé du Problème
Le processus de capture neutronique rapide (processus r) est responsable de la synthèse des éléments lourds dans l'univers. Bien que les profils d'abondance des éléments lourds (A > 130) dans les étoiles pauvres en métaux enrichies par le processus r et les résidus du processus r solaire présentent une robustesse remarquable, il existe une variabilité étoile par étoile significative pour les éléments lourds plus légers (Sr, Y, Zr, jusqu'au second pic). Les tentatives précédentes pour expliquer ces profils ont suggéré une variété de mécanismes, incluant des processus primaires, la fission, des processus r faibles ou des contributions de multiples sites astrophysiques. Cependant, déterminer les conditions astrophysiques spécifiques (fraction électronique, entropie, échelle de temps d'expansion) requises pour reproduire les abondances stellaires observées reste un défi. Les études existantes reposent souvent sur l'ajustement d'un nombre fixe de composantes ou de trajectoires spécifiques sans un cadre statistique rigoureux pour quantifier les incertitudes ou déterminer le nombre nécessaire de conditions contributives.

Méthodologie
Les auteurs présentent un cadre d'inférence bayésienne utilisant l'échantillonnage par chaîne de Markov de Monte Carlo (MCMC) pour inférer les conditions astrophysiques du processus r. La méthodologie intègre les composants suivants :

• Grille de Nucléosynthèse : L'analyse repose sur un large sondage indépendant du site (KAR25) contenant 120 000 calculs de nucléosynthèse paramétriques. Ces calculs utilisent le réseau de réactions WinNet et le modèle de masse FRDM2012, paramétrés par la fraction électronique initiale ($Y_e$), l'entropie initiale ($s$) et l'échelle de temps d'expansion ($\tau$).
• Modèle Statistique : Le profil d'abondance stellaire observé est modélisé comme une superposition linéaire de $N_{comp}$ composantes de nucléosynthèse. Chaque composante est définie par quatre paramètres : $Y_{e,i}$, $s_i$, $\tau_i$, et un poids $w_i$. L'abondance totale est donnée par $Y_{sample}(Z) = \sum w_i \cdot Y_{model,i}$.
• Configuration de l'Inférence : Les auteurs utilisent le package Python emcee pour échantillonner la distribution de probabilité a posteriori des paramètres. La fonction de vraisemblance est définie via une métrique de chi-deux réduit comparant le modèle aux abondances observées (normalisées par rapport au Strontium).
• Priors (Lois A Priori) : Deux priors distincts sont testés :
  1. Prior Uniforme Agnostique : Assigne une probabilité égale à tous les paramètres dans les limites physiques ($Y_e \in [0,005, 0,5]$, plages spécifiques pour $s$ et $\tau$).
  2. Prior Informé par la Simulation : Construit à l'aide d'une estimation de la densité par noyau gaussien (KDE) basée sur plus de 29 000 trajectoires hydrodynamiques provenant de 22 simulations différentes (couvrant les fusions d'étoiles à neutrons, les éjectas de disques et les supernovas magnétorotationnelles). Ce prior pondère les conditions selon leur fréquence dans les simulations actuelles.
• Identifiabilité : Pour éviter le double comptage de composantes identiques, l'échantillonneur MCMC impose une contrainte d'ordonnancement où les composantes doivent avoir des fractions électroniques croissantes ($Y_{e,1} \le \dots \le Y_{e,N_{comp}}$).

Résultats Clés
Le cadre a été appliqué à l'étoile template quasi complète du processus r HD222925, aux étoiles limitées-r (HD128279, HD122563) et aux résidus du processus r solaire.

1. Dominance de deux composantes : Pour HD222925, deux composantes suffisent pour reproduire le profil d'abondance observé avec une grande fidélité.

   • Composante H : Une condition de $Y_e$ plus faible ($Y_e \approx 0,16-0,18$) produisant les éléments du second au troisième pic (processus r lourd).
   • Composante L : Une condition de $Y_e$ plus élevée ($Y_e \approx 0,32-0,33$) produisant les éléments du premier au second pic (processus r léger).
   • Augmenter le nombre de composantes ($N_{comp} = 3, 4, 5$) n'améliore pas significativement l'accord médian avec les observations, car les composantes additionnelles échantillonnent largement les mêmes régions de l'espace des paramètres que la solution à deux composantes.

2. Universalité des Conditions : Les deux mêmes ensembles de conditions astrophysiques (composantes H et L) dominent les profils du processus r dans HD222925, HD128279, HD122563 et les résidus du processus r solaire. Les différences entre ces étoiles sont principalement attribuées au poids relatif ($w_2/w_1$) de ces deux composantes plutôt qu'à l'existence de sites nucléosynthétiques fondamentalement différents.

   • HD222925 : $w_2/w_1 \approx 1,7$
   • HD128279 : $w_2/w_1 \approx 2,7$
   • HD122563 : $w_2/w_1 \approx 7,4$
   • R solaire : $w_2/w_1 \approx 4,9$

3. Écarts du Processus r Solaire : Bien que le modèle à deux composantes s'ajuste bien à HD222925, il échoue à reproduire les éléments de capture neutronique plus légers ($Z \approx 30-35$) dans les résidus du processus r solaire, qui sont significativement enrichis par rapport à HD222925. Une troisième composante avec un $Y_e$ plus élevé ($\approx 0,42$) et une échelle de temps d'expansion plus longue est nécessaire pour ajuster le premier pic solaire.

4. Écarts Résiduels : Les désaccords restants entre les modèles et les observations sont concentrés dans des éléments spécifiques :

   • Systématiques d'Observation : Des éléments comme Ag et Cd ($Z=47, 48$) présentent des écarts probablement dus à des effets 3D et de non-équilibre thermodynamique local (NLTE).
   • Physique Nucléaire : Les écarts autour du troisième pic (Re, Os, Ir) et pour Th suggèrent une sensibilité aux modèles de masse nucléaire et aux barrières de fission.
   • Ga, Ge, As : Les abondances solaires de ces éléments sont plus élevées que celles de HD222925, suggérant des canaux de production additionnels dans le Soleil ou des problèmes dans l'extraction du processus r solaire.

5. Comparaison des Priors : Les résultats dérivés du prior uniforme agnostique et du prior informé par la simulation sont presque identiques. Cela indique que les conditions requises pour reproduire HD222525 sont déjà communes dans les simulations hydrodynamiques actuelles et ne nécessitent pas de conditions exotiques en dehors de l'espace des paramètres simulés.

Signification et Revendications
L'article affirme démontrer la puissance de la combinaison des observations d'abondances stellaires avec de larges calculs de nucléosynthèse via l'inférence bayésienne pour contraindre les conditions du processus r. La principale signification réside dans le fait de :

• Établir que la superposition de seulement deux conditions astrophysiques distinctes (une composante légère et une lourde) peut expliquer les profils robustes du processus r observés dans les étoiles pauvres en métaux et le système solaire, soutenant ainsi l'« universalité » du processus r.
• Fournir un cadre statistique qui quantifie les incertitudes et détermine le nombre nécessaire de composantes sans hypothèses arbitraires.
• Souligner que les écarts restants sont probablement dus à des systématiques d'observation (effets NLTE, 3D) et aux entrées de physique nucléaire (masses, fission) plutôt qu'à un manque de compréhension des sites astrophysiques fondamentaux.
• Motiver la nécessité d'échantillons plus larges et analysés de manière homogène d'étoiles enrichies par le processus r, ainsi que l'amélioration des données nucléaires/atomiques pour affiner davantage ces contraintes.

Les auteurs ne prétendent pas avoir identifié le site astrophysique spécifique (par exemple, fusion d'étoiles à neutrons vs supernova), mais plutôt les conditions thermodynamiques requises, notant que ces conditions se retrouvent dans divers scénarios de simulation actuels.