$^3$H and $^3$He nuclei production in a combined thermal and coalescence framework for heavy-ion collisions in the few-GeV energy regime

Cet article présente un modèle combiné thermique et de coalescence pour les collisions d'ions lourds dans le régime de quelques GeV qui reproduit avec succès les rendements en protons, pions et deutérons, mais sous-estime la production de noyaux $^3$H et $^3$He d'un facteur deux par rapport aux données expérimentales.

L'essentiel

Imaginez une collision de particules à haute énergie comme une fête dansante massive et chaotique où des milliers de petites particules (protons et neutrons) sont créées dans un éclair. Pendant une fraction de seconde, elles forment une soupe d'énergie chaude et tourbillonnante. À mesure que la fête se calme, ces particules tentent de trouver des partenaires pour former des groupes stables, comme des couples ou de petites équipes de danse.

Ce document traite de la tentative de prédire la fréquence à laquelle ces particules forment des groupes plus larges et spécifiques : le Tritium (un noyau avec un proton et deux neutrons, écrit ³H) et l'Hélium-3 (deux protons et un neutron, écrit ³He).

Voici la décomposition de ce que les scientifiques ont fait, en utilisant des analogies simples :

1. La recette en deux étapes

Les auteurs ont combiné deux façons différentes de concevoir la formation de ces groupes :

• Étape 1 : Le Modèle Thermique (la phase de la « soupe chaude ») :
  D'abord, ils ont utilisé un « modèle statistique ». Imaginez que la collision crée un immense bol de soupe chaude. Les particules dans cette soupe se déplacent de manière aléatoire. Les scientifiques ont calculé combien de protons et de pions (un autre type de particule) flottent autour en fonction de la température et de la pression de cette soupe. Ils savaient déjà que cette méthode fonctionnait bien pour prédire le nombre de particules seules et de paires (comme le Deutérium, qui est juste un proton et un neutron se tenant la main) produites.

• Étape 2 : Le Modèle de Coalescence (la phase du « regroupement ») :
  Ensuite, ils se sont demandé : « Si ces particules sont assez proches, vont-elles s'assembler pour former un trio ? » C'est ce qu'on appelle la coalescence. Pensez à un jeu de chaises musicales. Si trois joueurs (nucléons) se trouvent par hasard très proches les uns des autres lorsque la musique s'arrête (lorsque le système se fige), ils se donnent la main pour former une équipe (un noyau). Le document utilise les mathématiques pour calculer les probabilités que trois joueurs spécifiques soient assez proches pour former une équipe.

2. La configuration : Une balle légèrement aplatie

Les scientifiques n'ont pas simplement supposé que la « soupe » était une sphère parfaite. Ils ont réalisé que l'explosion issue de la collision ressemble plutôt à une balle légèrement aplatie (un sphéroïde) qui se développe vers l'extérieur. Ils ont utilisé une forme plus réaliste pour cette expansion, ce qui les a aidés à obtenir de meilleurs chiffres pour les particules individuelles (protons et pions) avant de tenter de prédire les trios.

3. Prédiction vs Réalité

L'équipe a lancé ses calculs pour prédire combien de noyaux de Tritium et d'Hélium-3 devraient être créés lors de collisions or-or à un niveau d'énergie spécifique (2,4 GeV).

• Le résultat : Leurs calculs prédisaient qu'il y aurait environ la moitié de ces noyaux par rapport à ce que l'expérience HADES (un détecteur réel) a effectivement observé.

  • Pour le Tritium (³H), ils ont prédit environ 3,16, mais l'expérience a trouvé 8,65.
  • Pour l'Hélium-3 (³He), ils ont prédit environ 2,26, mais l'expérience a trouvé 4,55.

• La bonne nouvelle : Même s'ils se sont trompés d'un facteur deux, ils ont visé le bon ordre de grandeur. Dans le monde de la physique des particules, prédire que vous obtiendrez « quelques » unités plutôt que « zéro » ou « un million » est un succès significatif. Cela prouve que leur idée combinée « Soupe chaude + Regroupement » est sur la bonne voie.

4. Pourquoi cet écart ?

Les auteurs suggèrent que le facteur manquant de deux provient de la manière dont ils ont calculé le « taux de formation ».

• L'analogie : Imaginez essayer de prédire combien de personnes vont former un groupe. Si vous supposez que tout le monde se tient en un cercle parfait, vous risquez de faire erreur dans vos calculs. Les scientifiques ont utilisé une forme simplifiée pour déterminer où se trouvent les particules (une sphère rigide). Ils soupçonnent que s'ils avaient utilisé une « fonction d'onde » plus complexe et plus réaliste (une meilleure carte de l'emplacement exact des particules), leur prédiction se rapprocherait des chiffres réels.

5. La forme des données

Bien que le nombre total de noyaux ait été sous-estimé, les scientifiques ont vérifié la forme des données (comment les particules sont réparties selon différentes vitesses et directions).

• Ils ont constaté que la forme de leur modèle était légèrement trop « abrupte » par rapport aux données expérimentales.
• Cependant, s'ils multipliaient simplement leur prédiction par un facteur d'échelle (comme en montant le volume), la forme de leur courbe correspondait très bien aux données expérimentales. Cela suggère que la physique de leur formation est correcte, même si le décompte exact nécessite un léger ajustement.

Résumé

Ce document est une tentative réussie de mélanger deux théories (soupe thermique et regroupements par coalescence) pour expliquer comment les noyaux lourds se forment lors de collisions de particules.

• Ce qui a fonctionné : Le modèle a correctement prédit l'ampleur générale de l'effet ainsi que la forme de la distribution des particules.
• Ce qui doit être amélioré : Le modèle prédit environ la moitié du nombre réel de noyaux trouvés dans les expériences. Les auteurs pensent que cela est dû au fait que leur « carte » mathématique de l'emplacement des particules est un peu trop simple, et qu'une carte plus détaillée corrigerait le décompte.

Ils concluent que leur cadre constitue une base solide pour comprendre ces petites équipes nucléaires, même si le décompte final nécessite un petit ajustement.

Résumé technique

Résumé technique : Production de noyaux $^3$H et $^3$He dans un cadre combiné thermique et de coalescence

Énoncé du problème
Ce travail traite de la production de noyaux légers, spécifiquement le tritium ($^3$H) et l'hélium-3 ($^3$He), dans les collisions d'ions lourds dans le régime d'énergie de quelques GeV. Alors que les modèles d'hadronisation statistique décrivent avec succès les spectres des protons et des pions, et que le modèle de coalescence est souvent utilisé pour les noyaux légers, ces approches sont fréquemment traitées comme des alternatives exclusives. Les auteurs visent à construire un cadre cohérent qui combine une description thermique de la production initiale de particules avec un mécanisme de coalescence ultérieur pour prédire les rendements et les spectres de $^3$H et $^3$He dans les collisions Au-Au à $\sqrt{s_{NN}} = 2,4$ GeV.

Méthodologie
Les auteurs emploient une approche hybride intégrant un modèle d'hadronisation statistique avec un cadre de coalescence de nucléons :

1. Modèle de dégel thermique (Thermal Freeze-out) :

   • L'étude utilise un modèle d'hadronisation statistique avec une expansion hydrodynamique sphéroïdale, ce qui est plus réaliste que la symétrie sphérique supposée dans les modèles de type "blast-wave" antérieurs.
   • Les conditions de dégel sont contraintes par les données expérimentales de la collaboration HADES (centralité de 10 %). L'analyse sélectionne une température de dégel plus élevée ($T = 70,3$ MeV) plutôt qu'une alternative plus basse ($T = 49,6$ MeV), car la première reproduit mieux le rendement expérimental du deutéron.
   • Le système est caractérisé par un petit rayon ($R \approx 6$ fm), un potentiel chimique baryonique ($\mu_B = 876$ MeV) et un potentiel chimique d'isospin ($\mu_{I3} = -21,5$ MeV).
   • Le modèle suppose que tous les protons détectés (tant libres que liés dans des noyaux légers) proviennent d'un système thermique unique. Par conséquent, les spectres de protons utilisés comme entrée pour la coalescence sont mis à l'échelle pour tenir compte de la fraction de protons se trouvant dans des états liés.

2. Cadre de coalescence :

   • La production de $^3$H et $^3$He est calculée en utilisant le formalisme de coalescence standard, où la distribution de quantité de mouvement du noyau est dérivée du produit des distributions de quantité de mouvement des protons et des neutrons à une fraction de la quantité de mouvement du noyau.
   • Distribution spatiale : Contrairement aux travaux précédents qui adoptent souvent un profil spatial gaussien pour la source, ce papier adopte une distribution spatiale uniforme des nucléons à l'intérieur d'une sphère de rayon $R$ à un temps de laboratoire fixe. Ce choix est cohérent avec les hypothèses du modèle thermique.
   • Fonctions d'onde : Les auteurs utilisent des approximations gaussiennes pour les fonctions d'onde de $^3$H et $^3$He en coordonnées de Jacobi.
   • Coefficient de taux de formation : Le coefficient de taux de formation du noyau ($A_{FR}^X$) est calculé en intégrant le produit de la fonction de distribution spatiale et du carré de la fonction d'onde sur l'espace des phases. Ce calcul produit des valeurs numériques élevées pour les coefficients ($A_{FR}^{^3H} \approx 5,51 \times 10^{11}$ MeV$^6$).

3. Implémentation numérique :

   • La formule de Cooper-Frye est utilisée pour générer les spectres de quantité de mouvement invariants pour les protons et les neutrons basés sur les paramètres d'expansion sphéroïdale.
   • Ces spectres de nucléons servent d'entrées dans l'équation de coalescence (Éq. 6 et Éq. 7) pour prédire les spectres de rapidité et de masse transversale des noyaux légers.

Résultats clés

• Rendements : Le modèle prédit des rendements totaux pour $^3$H et $^3$He qui sont environ 2 à 2,7 fois inférieurs aux valeurs expérimentales rapportées par la collaboration HADES. Spécifiquement, le modèle prédit $N \approx 3,16$ pour $^3$H et $N \approx 2,26$ pour $^3$He, contre des valeurs expérimentales de 8,65 et 4,55, respectivement.
• Formes de spectres : Malgré la divergence de normalisation, le modèle reproduit avec succès l'ordre de grandeur des multiplicités. Les formes prédites des spectres de rapidité et de masse transversale sont comparées aux données préliminaires de HADES.
• Facteurs d'échelle : Pour mieux correspondre aux formes des spectres de rapidité expérimentaux, les coefficients de taux de formation ont été mis à l'échelle par des facteurs ($\alpha$). Pour $^3$H, un facteur $\alpha \approx 2,4$ a fourni le meilleur ajustement à la forme spectrale, tandis que pour $^3$He, $\alpha \approx 1,9$ était optimal. Ces facteurs sont légèrement différents de ceux requis pour correspondre aux rendements totaux ($\alpha \approx 2,7$ pour $^3$H et $\alpha \approx 2,0$ pour $^3$He).
• Incertitudes : Les auteurs identifient la source principale d'incertitude comme étant le calcul des facteurs de taux de formation, qui repose sur des fonctions d'onde et des distributions spatiales simplifiées.

Signification et affirmations
Le papier affirme que sa contribution principale est de démontrer que les mécanismes de production thermique et de coalescence peuvent coexister et se compléter mutuellement dans un scénario théorique cohérent pour les collisions d'ions lourds à basse énergie.

• Les auteurs affirment que leur approche reproduit correctement l'ordre de grandeur des rendements des noyaux légers, ce qui est considéré comme un succès non trivial compte tenu de la complexité du calcul impliquant de grands et de petits nombres.
• Le travail étend la description réussie de la production de deutéron (établie précédemment dans la Réf. [4]) aux noyaux suivants les plus légers ($A=3$).
• Les auteurs concluent modestement que, bien que les rendements absols soient sous-estimés, le cadre fournit une base viable pour de futures comparaisons avec les données expérimentales. Ils suggèrent que la divergence des rendements pourrait être résolue en employant des fonctions d'onde plus réalistes pour $^3$H et $^3$He dans les calculs futurs du taux de formation.
• Le papier ne propose pas de nouveaux dispositifs expérimentaux mais souligne que ses prédictions de spectres peuvent être directement comparées avec de futures données expérimentales pour valider davantage le tableau combiné thermique-coalescence.