# Formule de Weizsäcker

La formule de Weizsäcker, appelée aussi formule de Bethe-Weizsäcker, est une formule semi-empirique donnant une valeur approximative de l'énergie de liaison nucléaire B caractérisant la liaison entre les nucléons qui constituent le noyau des atomes (voir un résumé dans Modèle de la goutte liquide).

## Histoire

L'éponyme de formule de Weizsäcker est le physicien allemand Carl Friedrich von Weizsäcker (1912-2007) qui l'a proposée en 1935$^,$ dans un article publié dans le Zeitschrift für Physik$^,$. Les physiciens Hans Bethe (1906-2005) et Robert Bacher (1905-2004) en ont simplifié l'expression en 1936$^,$$^,$$^,$. Le physicien Eugene Wigner (1902-1995) l'a généralisée en 1937$^,$$^,$$^,$. En 1939, Niels Bohr (1885-1962) et John Wheeler (1911-2008) y ont introduit le terme $\delta$$^,$.

## Expression

    $B\left(A,Z\right)=a_{v}A-a_{s}A^{\frac {2}{3}}-a_{c}{\frac {Z\left(Z-1\right)}{A^{\frac {1}{3}}}}-a_{a}{\frac {(A-2Z)^{2}}{A}}\pm a_{p}A^{-{\frac {1}{2}}}$
où :
* B est l'énergie de liaison, A est le nombre de masse (ou nombre de nucléons contenus dans le noyau A = Z+N), Z est le nombre de protons ;
* $a_{v}A$ est le terme de volume$^,$$^,$ ;
* $a_{s}A^{2/3}$ est le terme de surface$^,$$^,$ ;
* $a_{c}Z\left(Z-1\right)A^{-1/3}$ est le terme coulombien$^,$$^,$ ;
* $a_{a}\left(A-2Z\right)^{2}A^{-1}$ est le terme d'asymétrie$^,$$^,$$^,$ ;
* $a_{p}A^{-1/2}$ est le terme de parité ou d'appariement$^,$.
Remarques :
   Dans le terme coulombien, $Z\left(Z-1\right)$ est souvent remplacé par $Z^{2}$$^,$$^,$ car $Z\left(Z-1\right)\approx Z^{2}$.
   Dans certaines formules, le terme d'appariement varie en $A^{-3/4}$ plutôt qu'en $A^{-1/2}$.
Les valeurs des constantes utilisées sont (en MeV) :
* a$_v$ = 15,56
* a$_s$ = 17,23
* a$_c$ = 0,7
* a$_a$ = 23,6
* a$_p$ = 11,2
Il existe différents jeux de paramètres pour la formule de Weizsäcker. Le choix du jeu de paramètres se définit en fonction des noyaux étudiés. Ainsi certains jeux de paramètres donneront des énergies de liaison plus précises pour les noyaux stables tandis que d'autres jeux donneront des résultats plus satisfaisant pour les noyaux exotiques.
Cette formule permet d'expliquer avec une bonne approximation la courbe d'Aston (ci-contre), qui représente l'énergie de liaison par nucléon en fonction du nombre de nucléons dans le noyau. C'est une courbe expérimentale car chaque énergie de liaison y est calculée à partir de la masse du noyau atomique mesurée par spectrométrie de masse.

## Explication des différents termes

La formule de Bethe-Weizsäcker fait apparaitre cinq termes :
$\,B=E_{v}+E_{s}+E_{c}+E_{a}+E_{p}$
* Les deux premiers sont dus au modèle de la goutte liquide du noyau.
* Le troisième exprime la répulsion électrostatique entre les protons.
* Les deux derniers sont d'ordre quantique.
Pour expliquer ces différents termes, il faut supposer que le noyau est sphérique, de rayon $R_{0}$. Et comme il est compact (son volume est proportionnel au nombre de nucléons A), alors $R_{0}$ est proportionnel à $A^{1/3}$.

### Énergie de volume

Pour expliquer le premier terme, on peut utiliser une analogie avec un gaz parfait pour lequel l'énergie interne est proportionnelle au nombre de particules constituant le gaz. Ainsi, on pose que cette énergie de volume $E_{v}$ est proportionnelle à A. Elle permet d'expliquer les forces nucléaires de courtes portées, et la saturation des forces nucléaires : $\,E_{v}=a_{v}A$.

### Énergie de surface

La notion de tension de surface (ou superficielle) d'une goutte liquide peut être utilisée pour interpréter le second terme. Intuitivement, les nucléons à la surface du noyau sont en contact avec moins de nucléons que ceux du centre, l'énergie de liaison en est donc diminuée, à l'instar (sans que ce soit le même phénomène physique en jeu) de ce qu'il se passe à une interface liquide/gaz telle qu'une goutte d'eau dans l'air.
En introduisant E$_s$=−a$_s$A$^{2/3}${\textstyle \, E_{s}=-a_{s}A^{2/3}} comme le second terme de la formule, on prouve, en première approximation, que la surface du noyau est proportionnelle à E$_s${\textstyle E_{s}}.
Pour cela, on approxime que le volume du noyau est proportionnel au nombre de nucléons. Il s'agit d'une approximation courante :
$\,r=r_{0}A^{1/3}$Avec $\,r_{0}$ le rayon moyen d'un nucléon.
La surface d'une sphère de rayon r{\textstyle \, r} étant S=4πr$^2${\textstyle \, S=4\pi r^{2}}, en remplaçant r{\textstyle \, r} par son approximation, on obtient : $\,S=4\pi (r_{0}A^{1/3})^{2}=4\pi r_{0}^{2}A^{2/3}$D'où, $\,A^{2/3}\propto S$

### Répulsion électrostatique

Les protons étant tous chargés positivement, ils se repoussent mutuellement. Cela participe à diminuer l'énergie de liaison par un terme de répulsion électrostatique $E_{c}$. Dans une approximation grossière, le noyau peut être considéré comme une sphère avec une densité de charge uniforme. L'énergie potentielle d'une telle distribution de charge est donnée par :
    $E={\frac {3}{5}}\left({\frac {1}{4\pi \epsilon _{0}}}\right){\frac {Q^{2}}{R}}$
où $Q$ est la charge totale, $R$ le rayon de la sphère. En identifiant $Q$ à $Ze$ et en prenant le rayon proportionnel à $A^{1/3}$, on obtient la forme du terme coulombien. Cependant, la répulsion coulombienne n'existe que lorsqu'il y a plus d'un proton ce qui induit que $Z^{2}$ devient $Z(Z-1)$. La valeur de $a_{c}$ peut être calculée approximativement en utilisant l'équation ci-dessus :
    $R=r_{0}A^{\frac {1}{3}}$
    $Q=Ze$
    $Z^{2}=Z(Z-1)$
$E_{c}={\frac {3}{5}}\left({\frac {1}{4\pi \epsilon _{0}}}\right){\frac {Q^{2}}{R}}={\frac {3}{5}}\left({\frac {1}{4\pi \epsilon _{0}}}\right){\frac {(Ze)^{2}}{(r_{0}A^{\frac {1}{3}})}}={\frac {3e^{2}Z^{2}}{20\pi \epsilon _{0}r_{0}A^{\frac {1}{3}}}}={\frac {3e^{2}Z(Z-1)}{20\pi \epsilon _{0}r_{0}A^{\frac {1}{3}}}}=a_{c}{\frac {Z(Z-1)}{A^{1/3}}}$
L'énergie potentielle de la distribution de charge est donc :
    $E_{c}={\frac {3e^{2}Z(Z-1)}{20\pi \epsilon _{0}r_{0}A^{\frac {1}{3}}}}$
La constante $a_{c}$ du terme de répulsion électrostatique est :
    $a_{c}={\frac {3e^{2}}{20\pi \epsilon _{0}r_{0}}}$
Une autre valeur de $a_{c}$ peut être obtenue en utilisant la constante de structure fine :
    $a_{c}={\frac {3}{5}}\left({\frac {\hbar c\alpha }{r_{0}}}\right)$
où $\alpha$ est la constante de structure fine, $r_{0}A^{1/3}$, le rayon du noyau avec $r_{0}$ qui vaut approximativement 1,25 femtomètres. Cela donne une valeur théorique de $a_{c}$ de 0,691 MeV ce qui est peu éloigné des valeurs mesurées.
    $a_{c}=0,691MeV$

### Énergie d'asymétrie

La répulsion électrostatique étant en compétition avec l'interaction forte pour stabiliser le noyau, les noyaux lourds ont besoin d'un surplus de neutrons afin que cette interaction forte contrebalance l'effet de la répulsion électrostatique. Il y a donc une asymétrie du nombre de neutrons par rapport au nombre de protons. Cela n'a, a priori, aucun autre effet sur l'énergie de liaison que ceux qui ont été étudiés plus haut. En réalité, un effet quantique va jouer un rôle : les nucléons se trouvent sur des niveaux d'énergie, ce qui fait qu'un surplus de neutrons va augmenter leur énergie. On obtient alors que l'effet sur l'énergie de liaison s'écrit : $E_{a}=-a_{a}{\frac {(N-Z)^{2}}{A}}=-a_{a}{\frac {(A-2Z)^{2}}{A}}$.

### Énergie d'appariement

Un deuxième effet quantique joue un rôle dans l'énergie de liaison : les nucléons ayant un spin demi-entier ont tendance à s'apparier deux à deux, pour se grouper préférentiellement en nombre pair. Ainsi, un nombre impair de neutrons ou de protons sera moins stable. Une formule empirique permet de rendre compte de cet effet en ajoutant à l'énergie de liaison une énergie d'appariement (ou de parité) $E_{p}$ ayant différentes valeurs selon qu'il y a un nombre pair ou impair de nucléons, neutrons ou protons : $E_{p}=a_{p}\left\{{\begin{matrix}+A^{-{\frac {1}{2}}}&{\mbox{cas pair-pair (A pair, Z pair, N pair) }}\;\;\;\;\;\;\;\;\\0&{\mbox{cas pair-impair (A impair, Z et N de parité différente)}}\;\;\;\;\\-A^{-{\frac {1}{2}}}&{\mbox{cas impair-impair (A pair, Z impair, N impair) }}\end{matrix}}\right.$
Les noyaux ayant un nombre pair de nucléons, neutrons et protons sont plus stables que ceux ayant un nombre impair de nucléons, eux-mêmes plus stables que ceux ayant un nombre pair de nucléons et impair de neutrons et protons, donc l'énergie de liaison varie en conséquence.

## Utilisation de la formule
### Équation de la vallée de la stabilité

Pour un nombre de masse donné A, on s'aperçoit que la formule de Bethe-Weizsäcker fournit une équation quadratique en fonction de la charge Z. On a ainsi :
    $B_{A}(Z)=(-a_{c}A^{-{\frac {1}{3}}}-4a_{a}A^{-1})Z^{2}+(a_{c}A^{-{\frac {1}{3}}}+4a_{a})Z+(a_{v}-a_{a})A-a_{s}A^{\frac {2}{3}}\pm a_{p}A^{-{\frac {1}{2}}}$
Par définition, les noyaux stables sont définis comme étant des noyaux qui maximisent l'énergie de liaison $B_{A}(Z)$. Ainsi en cherchant les valeurs de $Z$ qui annulent la dérivée $B_{A}(Z)$ par rapport à $Z$, on peut obtenir une équation donnant les noyaux de la vallée de stabilité.
    ${\frac {\partial B_{A}(Z)}{\partial Z}}=0\rightarrow Z={\frac {a_{c}A^{-{\frac {1}{3}}}+4a_{a}}{2(a_{c}A^{-{\frac {1}{3}}}+4a_{a}A^{-1})}}\approx {1 \over 2}{A \over 1+A^{2/3}{\frac {a_{c}}{4a_{a}}}}$

## Améliorations

De nombreux modèles ont par la suite tenter de raffiner le modèle de la goutte liquide proposé initialement par Carl Friedrich von Weizsäcker en introduisant de nouveaux paramètres afin de reproduire plus finement les masses mesurées expérimentalement.

### Modèle de la gouttelette liquide

Ce modèle a été proposé par William D. Myers et Wladyslaw J. Swiatecki dans les années 1970. Il ajoute deux paramètres supplémentaires au modèle de la goutte liquide, à savoir la compressibilité de la matière nucléaire et l’asymétrie locale proton-neutron.

### Modèle de Duflo–Zuker

Proposé initialement en 1995 par Jean Duflo et Andres Zuker puis raffiné en 1999.