A HINSA view of cosmic-ray ionization in IC 348 and NGC 1333: evidence for a strong low-energy cosmic-ray disparity

Résumé technique : Une perspective HINSA sur l'ionisation par les rayons cosmiques dans IC 348 et NGC 1333

Énoncé du problème
Le taux d'ionisation par les rayons cosmiques (CRIR, $\zeta_{ion}$) est un paramètre fondamental qui régit l'évolution chimique et dynamique des nuages moléculaires, notamment en régulant le couplage entre les champs magnétiques et le gaz. Bien que des observations récentes aient révélé des valeurs de CRIR élevées dans les régions de formation d'étoiles massives et à proximité des protostars, les origines physiques et les mécanismes d'accélération des rayons cosmiques de basse énergie (LECR, $E < 1$ GeV) restent incertaines. De plus, une disparité existe dans le complexe de nuages moléculaires de Persée tout proche : alors que NGC 1333 présente un déficit de résidus de rayons $\gamma$ diffus sub-GeV, sa voisine IC 348 présente un résidu positif et un spectre d'énergie de CR plus dur, malgré le fait que NGC 1333 possède une activité de formation stellaire plus élevée. La mesure dans laquelle les sources locales de CR supplantent le fond galactique dans ces régions, ainsi que la nature de la disparité entre elles, nécessite d'être clarifiée.

Méthodologie
Cette étude présente de nouvelles estimations du CRIR vers IC 348 et NGC 1333 en utilisant la technique HINSA (H I narrow self-absorption). La méthodologie implique :

1. Observations : Utilisation de données H I provenant du télescope FAST (NGC 1333) et du relevé GALFA-H I (IC 348), combinées à des données archivistiques de $^{13}$CO ($J=1-0$ et $3-2$) issues du relevé COMPLETE et du JCMT, ainsi qu'à des données de continuum de poussière du Herschel Gould Belt Survey.
2. Analyse HINSA : Dérivation de la colonne de densité d'hydrogène atomique $N(\text{HINSA})$ via l'ajustement par dérivée de second ordre des caractéristiques d'absorption H I par rapport au champ de vitesse de $^{13}$CO. La fraction HINSA ($f_{\text{HINSA}} = N(\text{HINSA})/N(\text{H}_2)$) est calculée pour les régions où $N(\text{H}_2) > 3 \times 10^{21}$ cm$^{-2}$ afin de minimiser les effets de photodissociation.
3. Contraintes de densité : Estimation de la densité volumique du gaz $n(\text{H}_2)$ par modélisation de transfert radiatif non-LTE (RADEX) contrainte par les intensités de raies de $^{13}$CO et un échantillonnage par chaîne de Markov de Monte Carlo (MCMC).
4. Calcul du CRIR : Dérivation de $\zeta_{ion}$ à partir de l'équilibre en régime permanent entre la formation de H$_2$ sur les grains de poussière et sa destruction par les rayons cosmiques (incluant l'ionisation, la dissociation et les processus UV induits par les CR).
5. Modélisation spectrale : Modélisation du spectre de LECR dans un cadre d'atténuation de lame finie. Le modèle incorpore des contraintes issues des spectres de CR de haute énergie déduits des observations de rayons $\gamma$ de Fermi pour ajuster les profils de $\zeta_{ion}$ observés, en faisant varier la pente spectrale de basse énergie ($\alpha$) et la rupture de quantité de mouvement ($p_0$).

Résultats clés

• Disparité de CRIR : L'étude trouve une différence prononcée de CRIR entre les deux nuages. Dans IC 348, la valeur moyenne pondérée de $\zeta_{ion}$ varie de $(5,35 \pm 1,49) \times 10^{-17}$ s$^{-1}$ à $(2,33 \pm 0,55) \times 10^{-16}$ s$^{-1}$. En revanche, NGC 1333 présente des valeurs environ un ordre de grandeur plus faibles, allant de $(3,94 \pm 2,46) \times 10^{-18}$ s$^{-1}$ à $(1,68 \pm 0,25) \times 10^{-17}$ s$^{-1}$.
• Tendance de la colonne de densité : Dans les deux nuages, $\zeta_{ion}$ diminue avec l'augmentation de la colonne de densité de H$_2$, ce qui est cohérent avec les modèles d'atténuation dépourvus de sources d'accélération internes.
• Modélisation spectrale : Les meilleurs spectres de CR de basse énergie reproduisent les profils observés. Le modèle indique que le facteur de $\sim$10 de différence dans les taux d'ionisation est piloté par une différence correspondante dans le flux de protons de basse énergie ($\sim$10 MeV). Le spectre d'IC 348 est nettement plus dur et plus intense que celui de NGC 1333, qui présente un déficit par rapport au fond galactique aux énergies du GeV.
• Analyse de la densité et de la structure : Les auteurs examinent les biais potentiels provenant de l'estimation de la densité du gaz et de la structure des nuages (via des simulations 3D-PDR). Bien que les incertitudes sur $n(\text{H}_2)$ et les structures fractales des nuages puissent biaiser les valeurs absolues de $\zeta_{ion}$ (potentiellement en les sous-estimant à faible $N(\text{H}_2)$ et en les surestimant à haute $N(\text{H}_2)$), l'analyse conclut que ces facteurs ne peuvent expliquer la disparité d'un ordre de grandeur entre les deux régions.

Signification et affirmations
L'article affirme que la disparité observée de CRIR est intrinsèque aux propriétés physiques des nuages plutôt qu'un artefact des échelles de temps évolutives ou des biais de mesure.

• Sources d'accélération externes : Le CRIR élevé dans IC 348, couplé à son spectre de rayons $\gamma$ plus dur, suggère la présence de sources d'accélération de LECR locales externes aux cœurs denses des nuages, probablement distinctes de l'activité protostellaire qui domine NGC 1333. Les auteurs proposent que les chocs de terminaison des vents stellaires dans l'amas d'IC 348, plus évolué et moins enfoui (contenant plus d'étoiles de type B), sont une source plausible.
• Contribution protostellaire : Les résultats suggèrent que si les jets protostellaires peuvent augmenter le CRIR à petite échelle (comme observé dans la cartographie de NGC 1333), ils ne pilotent pas le CRIR à grande échelle dans ces environnements, compte tenu notamment du fait que NGC 1333 possède une activité de formation stellaire plus élevée tout en ayant un CRIR global plus faible.
• Implications pour la formation stellaire : Le CRIR globalement plus élevé dans IC 348 implique une fraction d'ionisation plus élevée, conduisant à un couplage plus fort entre le gaz et le champ magnétique. Les auteurs notent que cela peut accentuer le freinage magnétique lors de l'effondrement des cœurs, ce qui pourrait expliquer les observations selon lesquelles les disques protostellaires dans IC 348 sont en moyenne plus petits que ceux de NGC 1333.
• Validation méthodologique : L'étude valide la technique HINSA comme une sonde pour le CRIR dans les régions de formation d'étoiles de faible masse, soulignant la nécessité de futures modélisations 3D pour affiner les valeurs absolues tout en confirmant la robustesse des disparités relatives.