Structure of the $^8$B and $^8$Li nuclei and the astrophysical $S_{17}(0)$-factor of the $^7$Be($p,\gamma$)$^8$B direct capture process within a three-body model

En utilisant un modèle de cluster à potentiel à trois corps avec la méthode de Lagrange-mesh hypersphérique, cette étude calcule les propriétés structurelles et les coefficients de normalisation asymptotique des noyaux $^8$B et $^8$Li afin de dériver un facteur astrophysique $S_{17}(0)$ à énergie nulle précis de $22.492 \pm 0.014$ eV b pour la réaction $^7$Be($p,\gamma$)$^8$B, révélant que le canal de spin 2 domine le processus de capture.

L'essentiel

Imaginez le noyau atomique non pas comme une bille solide, mais comme une minuscule piste de danse chaotique où les particules tournent constamment et se tiennent par la main. Cet article est une étude détaillée de deux danseurs spécifiques sur cette piste : les noyaux Bore-8 (8B) et Lithium-8 (8Li).

Les auteurs, travaillant depuis l'Ouzbékistan, voulaient comprendre exactement comment ces noyaux sont construits et comment ils se comportent lorsqu'ils interagissent avec d'autres particules. Voici le détail de leur travail en termes simples.

1. Le Décor : Une Danse à Trois

La plupart des gens imaginent un noyau comme une masse unique, mais les auteurs traitent ces noyaux spécifiques comme un système à trois corps.

• Les Danseurs : Ils imaginent le noyau comme un groupe de trois parties distinctes : une particule Alpha (un groupe serré de 2 protons et 2 neutrons), un noyau Hélium-3 ou Tritium (un groupe plus petit), et un seul proton ou neutron.
• Le Modèle : Ils ont utilisé une « piste de danse » mathématique appelée méthode du maillage de Lagrange hypersphérique. Imaginez cela comme une grille 3D ultra-précise qui leur permet de calculer exactement comment ces trois parties se déplacent et s'accrochent les unes aux autres sans entrer dans des zones interdites (un concept appelé « principe d'exclusion de Pauli », qui est comme une règle disant que deux danseurs ne peuvent pas occuper exactement le même endroit en même temps).

2. L'Objectif : Mesurer la « Prise » (ANC)

La principale chose que les chercheurs voulaient mesurer est quelque chose appelé le Coefficient de Normalisation Asymptotique (ANC).

• L'Analogie : Imaginez que le noyau est un aimant. L'ANC mesure la force de l'attraction magnétique au tout bord de l'aimant, juste au moment où un morceau de fer est sur le point de s'y accrocher.
• Pourquoi c'est important : Dans le monde des étoiles, les noyaux tentent constamment de s'accrocher les uns aux autres pour créer de l'énergie. Pour savoir avec quelle probabilité ils vont s'accrocher, vous devez connaître exactement la force de cette « prise au bord ». Si la prise est trop faible, ils rebondissent ; si elle est juste, ils fusionnent.

L'équipe a calculé cette « force de prise » pour deux scénarios différents :

1. Bore-8 : À quel point un proton s'accroche-t-il fermement à un noyau Béryllium-7 ?
2. Lithium-8 : À quel point un neutron s'accroche-t-il fermement à un noyau Lithium-7 ?

Ils ont constaté que la « prise » varie selon le spin des particules (comme si les danseurs tournaient dans le sens des aiguilles d'une montre ou dans le sens inverse). Ils ont calculé ces valeurs avec une grande précision, s'assurant que leurs mathématiques convergent (cessaient de changer) lorsqu'ils ajoutaient suffisamment de détails au modèle.

3. La Grande Question : Le Thermostat Solaire

La raison ultime de cette étude est de résoudre un mystère concernant le Soleil.

• La Réaction : Le Soleil brille grâce à une réaction en chaîne où le Béryllium-7 attrape un proton pour devenir du Bore-8. Cette étape est le « goulot d'étranglement » du processus.
• Le Problème : Nous ne pouvons pas facilement mesurer cette réaction en laboratoire car le cœur du Soleil est incroyablement chaud, mais la réaction se produit à des énergies très faibles où la répulsion électrique entre les particules agit comme un mur massif.
• La Solution : En calculant parfaitement la « force de prise » (ANC) dans leur modèle, ils ont pu prédire le facteur astrophysique S. Imaginez le facteur S comme un « score de probabilité » indiquant à quelle fréquence cette fusion se produit.

4. Les Résultats : Un Nouveau Chiffre pour le Soleil

L'équipe a calculé un nombre spécifique pour cette probabilité : 22,492 eV b.

Voici comment leur résultat se compare aux « livres de règles » utilisés par les scientifiques :

• Solar Fusion II (L'ancien livre de règles) : Suggérait une valeur autour de 20,8. Le résultat des auteurs est un peu plus élevé que cela.
• Solar Fusion III (Le livre de règles plus récent) : Suggérait une valeur de 20,5. Le résultat des auteurs est nettement plus élevé que cela.
• Le « Meilleur » Modèle Solaire (BAR2M) : Fait intéressant, le modèle solaire moderne le plus réussi utilise actuellement une valeur de 22,4.

La Conclusion : Le calcul des auteurs (22,49) correspond presque parfaitement à la valeur utilisée dans le modèle solaire actuel le plus réussi (22,4). Cela suggère que leur façon de modéliser la danse à trois corps est très précise et soutient l'idée que la température interne du Soleil et sa production d'énergie pourraient être légèrement différentes de ce que le « livre de règles Solar Fusion III » suggère.

Résumé

En bref, les auteurs ont construit une simulation mathématique hautement détaillée de la manière dont les noyaux Bore-8 et Lithium-8 sont construits. En mesurant exactement à quel point leurs particules externes sont retenues, ils ont calculé une probabilité spécifique pour une réaction nucléaire qui alimente le Soleil. Leur chiffre correspond aux modèles solaires modernes les plus réussis, suggérant que notre compréhension actuelle du « moteur » du Soleil pourrait devoir être légèrement ajustée pour correspondre à leurs découvertes.

Résumé technique

1. Énoncé du problème

La réaction de capture radiative $^{7}\text{Be}(p, \gamma)^{8}\text{B}$ constitue un goulot d'étranglement critique dans la chaîne proton-proton (pp) solaire, déterminant le flux de neutrinos solaires de haute énergie. Un défi majeur en astrophysique nucléaire est la détermination précise du facteur astrophysique à énergie nulle, $S_{17}(0)$.

• Limites expérimentales : Une mesure directe aux énergies stellaires est impossible en raison de la barrière de Coulomb et de sections efficaces extrêmement faibles.
• Discordances théoriques : Les modèles théoriques existants et les extractions expérimentales indirectes (via des réactions de transfert ou de fragmentation) donnent des valeurs contradictoires pour les Coefficients de Normalisation Asymptotique (ANC) et le $S_{17}(0)$ qui en découle.
  • Solar Fusion II recommande $S_{17}(0) \approx 20.8 \pm 0.7 \text{ (th)} \pm 1.4 \text{ (exp)}$ eV b.
  • Solar Fusion III recommande une valeur légèrement plus basse de $20.5 \pm 0.70$ eV b.
  • Cependant, le modèle solaire moderne le plus abouti (BAR2M) repose sur une valeur de 22.4 eV b.
• Objectif : Les auteurs visent à résoudre ces discordances en calculant la structure de $^{8}\text{B}$ et de son noyau miroir $^{8}\text{Li}$ à l'aide d'un modèle à trois corps auto-cohérent, afin de dériver des valeurs précises d'ANC et le $S_{17}(0)$ correspondant.

2. Méthodologie

Les auteurs emploient un modèle de clusters à potentiel à trois corps traitant les noyaux comme suit :

• $^{8}\text{B}$ : $\alpha + ^{3}\text{He} + p$
• $^{8}\text{Li}$ : $\alpha + ^{3}\text{H} + n$

Composantes techniques clés :

• Méthode de Lagrange-Mesh hypersphérique : L'équation de Schrödinger est résolue en coordonnées hypersphériques. La fonction d'onde est développée sur une base d'harmoniques hypersphériques.
  • La convergence a été atteinte avec un moment hyperangulaire maximal de $K_{\text{max}} = 22$ pour les états fondamentaux ($2^+$) et $K_{\text{max}} = 28$ pour les états excités ($1^+$).
• Potentiels d'interaction : Des potentiels à deux corps réalistes ont été adoptés de la littérature :
  • Noyau-$\alpha$ ($\alpha$-p, $\alpha$-n) : Potentiel de Voronchev et al. (incluant le couplage spin-orbite et les états interdits par le principe de Pauli).
  • $\alpha$-$^{3}\text{He}$ / $\alpha$-$^{3}\text{H}$ : Potentiels gaussiens avec états interdits par le principe de Pauli.
  • $p$-$^{3}\text{He}$ / $n$-$^{3}\text{H}$ : Potentiels attractifs gaussiens + répulsifs exponentiels.
  • Note : Des pseudo-potentiels orthogonalisants (OPP) ont été utilisés pour éliminer les états interdits par le principe de Pauli.
• Extraction des ANC : Les valeurs d'ANC pour les désintégrations virtuelles ($^{8}\text{B} \to ^{7}\text{Be} + p$ et $^{8}\text{Li} \to ^{7}\text{Li} + n$) ont été déterminées en faisant correspondre l'intégrale de recouvrement de la fonction d'onde à trois corps avec la fonction de Whittaker asymptotique à grandes distances inter-clusters.
• Calcul du facteur astrophysique : Le facteur $S_{17}(0)$ a été calculé en utilisant la théorie asymptotique développée par D. Baye, qui relie le facteur S à la longueur de diffusion en onde S et au carré de l'ANC ($C^2$).

3. Contributions et résultats clés

A. Propriétés de la structure nucléaire

Le modèle a reproduit avec succès les énergies de liaison expérimentales sans nécessiter de forces à trois corps :

• Énergies de liaison :
  • $^{8}\text{B}$ ($2^+$) : Calculée $-1.718$ MeV vs Exp. $-1.724$ MeV.
  • $^{8}\text{Li}$ ($2^+$) : Calculée $-4.340$ MeV vs Exp. $-4.499$ MeV.
• Rayons de matière :
  • $^{8}\text{B}$ : Calculé $2.76$ fm (Exp. $2.58 \pm 0.06$ fm). Le modèle surestime légèrement le rayon d'environ 5 %, ce qui est cohérent avec sa nature de halo.
  • $^{8}\text{Li}$ : Calculé $2.56$ fm (Exp. $2.50 \pm 0.06$ fm), bien dans la marge d'erreur expérimentale.

B. Coefficients de Normalisation Asymptotique (ANC)

L'étude a fourni des valeurs d'ANC auto-cohérentes pour les canaux de spin $S=1$ et $S=2$ :

• $^{8}\text{B} \to ^{7}\text{Be} + p$ :
  • $C_{S=1} = 0.211 \text{ fm}^{-1/2}$
  • $C_{S=2} = 0.739 \text{ fm}^{-1/2}$
  • ANC carrée totale : $C^2 = 0.591 \text{ fm}^{-1}$.
• $^{8}\text{Li} \to ^{7}\text{Li} + n$ :
  • $C_{S=1} = 0.220 \text{ fm}^{-1/2}$
  • $C_{S=2} = 0.774 \text{ fm}^{-1/2}$
  • ANC carrée totale : $C^2 = 0.648 \text{ fm}^{-1}$.

Comparaison : La $C^2$ calculée pour $^{8}\text{B}$ ($0.591$) est supérieure au résultat NCSM ab-initio ($0.509$) mais s'aligne bien avec les valeurs extraites de la réaction $^{7}\text{Be}(d,n)^{8}\text{B}$ ($0.613 \pm 0.060$).

C. Facteur astrophysique S ($S_{17}(0)$)

En utilisant les ANC dérivées et les longueurs de diffusion expérimentales, le facteur astrophysique à énergie nulle a été calculé :

• Canal de spin 2 (Dominant) : $S^{(2)}_{17}(0) = 20.838 \pm 0.014$ eV b.
• Canal de spin 1 : $S^{(1)}_{17}(0) = 1.654 \pm 0.003$ eV b.
• Résultat total : $S_{17}(0) = 22.492 \pm 0.014$ eV b.

4. Signification et conclusion

• Résolution de la discordance du modèle solaire : La valeur calculée de 22.492 eV b est en excellent accord avec la valeur (22.4 eV b) utilisée dans le modèle solaire BAR2M, actuellement considéré comme le modèle le plus abouti pour les prédictions de neutrinos solaires.
• Comparaison avec les revues :
  • Le résultat est cohérent avec l'estimation Solar Fusion II ($20.8 \pm \dots$).
  • Il est légèrement supérieur à la recommandation Solar Fusion III ($20.5 \pm 0.70$ eV b), suggérant que la valeur plus basse dans SFIII pourrait sous-estimer la contribution du canal de spin 2 ou l'amplitude spécifique de l'ANC dérivée de cette approche rigoureuse à trois corps.
• Validation méthodologique : L'étude démontre que la méthode de Lagrange-Mesh hypersphérique, combinée à des potentiels à deux corps réalistes, fournit un cadre robuste et auto-cohérent pour décrire les noyaux légers de halo et extraire des facteurs astrophysiques sans forces à trois corps ajustables.
• Implication pour la physique des neutrinos : Les résultats soutiennent la valeur plus élevée de $S_{17}(0)$, ce qui a des implications directes sur le flux prédit de neutrinos solaires $^{8}\text{B}$ et l'interprétation des données de neutrinos solaires.

En résumé, cet article fournit une confirmation théorique rigoureuse que le facteur astrophysique S pour $^{7}\text{Be}(p, \gamma)^{8}\text{B}$ est probablement plus proche de 22.5 eV b, comblant le fossé entre les modèles théoriques de structure nucléaire et les exigences de la modélisation solaire moderne.