Probing in-medium effect via giant dipole resonance in the extended quantum molecular dynamics model

Cet article utilise une approche stochastique dans le modèle EQMD étendu pour démontrer que la largeur de la résonance dipolaire géante dans le ${}^{208}$Pb est sensible à l'énergie de symétrie et nécessite une réduction significative des sections efficaces nucléon-nucléon in-medium pour correspondre aux données expérimentales.

L'essentiel

Imaginez que le noyau d'un atome (comme celui du plomb) est un peu comme une grosse boule de pâte à modeler composée de deux types de billes : des billes rouges (les protons) et des billes bleues (les neutrons).

Normalement, ces billes sont bien tassées ensemble. Mais si vous donnez un petit coup sec à cette boule, les billes rouges et les billes bleues vont se mettre à osciller l'une contre l'autre, un peu comme les deux côtés d'un ressort qui vibre. En physique, on appelle ce phénomène la Résonance Dipolaire Géante (GDR). C'est une sorte de "battement de cœur" du noyau.

L'objectif de cet article est de comprendre comment ce battement de cœur s'arrête (comment il s'amortit) et de voir ce que cela nous apprend sur la "colle" qui maintient les billes ensemble.

Voici les points clés expliqués simplement :

1. Le problème de l'ancien modèle : La "boussole géométrique"

Pour simuler ce battement de cœur sur un ordinateur, les scientifiques utilisent des modèles mathématiques. L'ancien modèle (appelé EQMD) fonctionnait un peu comme un jeu de billard très strict :

• Il calculait si deux billes allaient se percuter en regardant si leurs trajectoires se croisaient exactement dans un petit espace de temps.
• Le problème : Cette méthode était trop rigide. Elle ne parvenait pas à reproduire la vitesse à laquelle le battement de cœur du noyau s'arrêtait réellement. C'était comme si le modèle prédisait que la boule de pâte à modeler vibrait pendant des heures, alors qu'en réalité, elle s'arrête très vite.

2. La nouvelle solution : Le "hasard intelligent"

Les auteurs de l'article ont décidé de changer de méthode. Au lieu de regarder les trajectoires avec une règle, ils ont utilisé une approche stochastique (basée sur le hasard et les probabilités).

• L'analogie : Imaginez que vous êtes dans une foule dense. Au lieu de calculer exactement à quelle seconde vous allez heurter quelqu'un, vous dites : "Plus je suis proche de quelqu'un et plus je me déplace vite, plus j'ai de chances de le bousculer".
• Cette nouvelle méthode permet de mieux simuler les collisions entre les billes à l'intérieur du noyau.

3. La découverte majeure : L'effet "miroir" du milieu

En utilisant cette nouvelle méthode, les chercheurs ont découvert deux choses fascinantes :

• La "colle" change de force : La force qui maintient les protons et les neutrons ensemble (appelée énergie de symétrie) détermine la hauteur du battement de cœur. En ajustant cette force dans leur modèle, ils ont pu faire correspondre la fréquence du battement à celle observée dans la réalité.
• Le frein invisible (l'effet de milieu) : C'est le point le plus important. Ils ont réalisé que pour que le battement de cœur s'arrête aussi vite que dans la réalité, les billes doivent se heurter beaucoup plus souvent qu'on ne le pensait.
  • L'analogie : Imaginez courir dans un couloir vide (l'espace libre). Vous courez vite. Maintenant, imaginez courir dans une piscine remplie de gens (le milieu nucléaire). Vous allez beaucoup moins vite car vous heurtez tout le monde.
  • Le modèle a montré que les collisions à l'intérieur du noyau sont beaucoup plus fréquentes que dans l'espace vide. Le "milieu" (la présence des autres billes) agit comme un frein puissant qui réduit la probabilité de collision, un peu comme si la résistance de l'eau ralentissait vos mouvements.

4. Pourquoi est-ce important ?

En réussissant à faire correspondre leur modèle aux données réelles (en trouvant le bon réglage pour la "colle" et le bon réglage pour le "frein"), les scientifiques ont prouvé que :

1. Leur nouvelle méthode de calcul (le hasard intelligent) est bien meilleure que l'ancienne.
2. Ils peuvent maintenant utiliser les battements de cœur des noyaux pour mesurer des choses qu'on ne peut pas voir directement, comme la façon dont la matière nucléaire se comporte sous des conditions extrêmes (comme dans les étoiles à neutrons).

En résumé :
Les chercheurs ont remplacé un vieux modèle de billard trop rigide par un nouveau modèle basé sur la probabilité. Grâce à cela, ils ont pu expliquer pourquoi le "cœur" des atomes bat et s'arrête si vite. Ils ont découvert que l'intérieur d'un atome est un endroit très encombré où les particules se freinent mutuellement, ce qui nous aide à mieux comprendre la recette secrète de la matière qui compose notre univers.

Résumé technique

1. Problématique

L'étude des réactions d'ions lourds (HIC) et des excitations collectives nucléaires repose souvent sur des approches de transport. Un élément crucial de ces modèles est le terme de collision, qui nécessite l'utilisation de la section efficace de diffusion nucléon-nucléon (NN).

• Limitation des modèles existants : Le modèle original EQMD (Extended Quantum Molecular Dynamics), développé par Maruyama et al. en 1996, utilise une approche géométrique pour traiter les collisions binaires. Cette méthode présente des défauts, tels que la possibilité de collisions artificielles (spurious collisions) et une dépendance arbitraire à des paramètres de seuil (distance minimale, section efficace maximale).
• Le défi physique : La section efficace NN dans le vide ($\sigma_{NN}^{free}$) ne suffit pas à décrire les phénomènes au sein de la matière nucléaire dense. Les effets du milieu (in-medium effects) réduisent cette section efficace, mais l'ampleur exacte de cette réduction reste une question ouverte.
• Objectif : L'article vise à explorer les effets du milieu en analysant la largeur de la Résonance Géante Dipolaire Isoscalaire (GDR) dans le noyau $^{208}\text{Pb}$. La GDR est sensible aux collisions binaires qui amortissent l'oscillation collective. Les auteurs cherchent à déterminer si l'approche géométrique originale de l'EQMD est capable de reproduire les données expérimentales de la largeur de la GDR, et à évaluer l'impact d'une nouvelle méthode stochastique couplée à une correction du milieu.

2. Méthodologie

Les auteurs ont mis à jour le modèle EQMD en intégrant deux changements majeurs :

• Modèle de base « Soft » EQMD : Ils utilisent une version récente du modèle EQMD qui a corrigé le terme de champ moyen (potentiel effectif) pour éviter une surévaluation de l'incompressibilité nucléaire. Le potentiel est basé sur une fonctionnelle de densité d'énergie de type Skyrme (paramètres SLy7), incluant un terme d'énergie de symétrie ($c_s$) et des interactions dépendantes de l'impulsion.
• Nouvelle approche de collision (Stochastique) : Au lieu de l'approche géométrique originale, les auteurs implémentent une méthode stochastique (équivalente à la méthode du libre parcours moyen).
  • La probabilité de collision est proportionnelle au recouvrement de densité des paquets d'ondes des nucléons.
  • La section efficace utilisée est la section élastique NN libre ($\sigma_{NN}^{free}$) paramétrée à partir de données expérimentales (réf. [71]), couvrant une large gamme d'énergies.
  • Une correction de blocage de Pauli est appliquée pour assurer la conservation des propriétés fermioniques.
• Modélisation de l'effet du milieu : Pour simuler la réduction de la section efficace dans le milieu nucléaire, une paramétrisation phénoménologique est appliquée :
  $$\sigma_{NN}^* = \left(1 - \alpha_1 \frac{\rho}{\rho_0}\right) \sigma_{NN}^{free}$$
  où $\alpha_1$ est le paramètre de force de la correction du milieu, $\rho$ la densité locale et $\rho_0$ la densité de saturation.
• Excitation de la GDR : Une perturbation dipolaire est appliquée au Hamiltonien pour exciter l'oscillation collective protons-neutrons. La fonction de réponse temporelle est ensuite transformée par Fourier pour obtenir la fonction de force et la largeur de la résonance (FWHM).

3. Contributions Clés

1. Transition Géométrique vers Stochastique : L'article démontre que l'approche stochastique est supérieure à l'approche géométrique originale de l'EQMD pour traiter les collisions binaires, en éliminant les artefacts numériques et en permettant une description plus réaliste de l'amortissement collisionnel.
2. Sensibilité de la GDR : Confirmation que la position du pic et la largeur de la GDR dans le $^{208}\text{Pb}$ sont des observables extrêmement sensibles à la fois au coefficient d'énergie de symétrie ($c_s$) et à la section efficace NN dans le milieu ($\sigma_{NN}^*$).
3. Extraction des paramètres du milieu : L'étude permet de contraindre quantitativement la réduction de la section efficace NN à basse énergie. Les résultats indiquent qu'une réduction significative est nécessaire pour correspondre aux données expérimentales.

4. Résultats

• Échec de l'approche géométrique originale : Les simulations utilisant l'approche géométrique de l'EQMD (avec ou sans coupure basse énergie) sous-estiment systématiquement la largeur de la GDR d'au moins 1 MeV par rapport aux données d'évaluation. Cela suggère que le terme de collision original est inadéquat pour reproduire l'amortissement dans la région de basse énergie.
• Performance de l'approche stochastique :
  • L'utilisation de la section efficace libre ($\sigma_{NN}^{free}$) sans correction du milieu conduit à un amortissement trop fort (largeur de GDR surestimée).
  • L'application de la correction du milieu avec un paramètre $\alpha_1 \approx 0.57$ permet de reproduire avec précision la largeur expérimentale de la GDR ($\Gamma \approx 4.08 \pm 0.09$ MeV).
  • La position du pic de la GDR ($E_\gamma \approx 13.38$ MeV) est correctement reproduite avec un coefficient d'énergie de symétrie $c_s \approx 33.2$ MeV.
• Stabilité du noyau : Les tests de stabilité montrent que la méthode stochastique ne perturbe pas la stabilité de l'état fondamental du noyau $^{208}\text{Pb}$ sur de longues périodes d'évolution (jusqu'à 1000 fm/c), validant son utilisation pour des réactions nucléaires à long terme.

5. Signification et Conclusion

Cette étude valide l'utilisation de l'approche stochastique dans le cadre du modèle EQMD pour l'étude des réactions nucléaires à basse et moyenne énergie. Elle démontre que la largeur de la Résonance Géante Dipolaire est un outil puissant pour sonder la physique nucléaire dense.

Les résultats confirment l'hypothèse d'une réduction significative des sections efficaces de diffusion nucléon-nucléon à l'intérieur du milieu nucléaire (environ 40-60% de réduction à saturation), en accord avec d'autres modèles comme le LBUU. Cela fournit une contrainte importante pour l'équation d'état (EoS) de la matière nucléaire et la compréhension des effets de milieu dans les collisions d'ions lourds. L'article ouvre la voie à une utilisation plus fiable du modèle EQMD pour explorer les propriétés de la matière nucléaire exotique et les phénomènes astrophysiques.